有限元计算注意事项

1、1、计算原理任何有限元模拟的第一步都是利用一个有限单元的集合离散结 构的实际几何形状，每个单元（element）代表这个实际结构的一个离 散部分.这些单元通过公用的节点（node）来连接.节点和单元的集合 成为网格（mesh.在一个特定网格中的单元数目称为网格密度（meshdensity）.在ansys计算过程中，程序以每个节点的每个自由度建立 平衡方程，以节点的位移作为未知量，利用矩阵求解节点的位移一旦 节点位移求出，整个结构的应力和应变都很容易计算出来.这种计算 的过程和方法，数学上称之为隐式方法.从上叙述来看，整个计算过程 中就是求解n个n元一次方程组（n表示节点数量），当计算模型复 杂

2、而且庞大时，隐式求解方法的计算量还是相当大的.与之相对应的，显式求解方法.显式求解方法是通过动态方法从 一个增量步前推到下一个增量步得到的.具体显式求解方法和隐式求 解方法例子如下:(1) 隐式求解节点a图1-4桁架问题的离散化榄型图1-5每个羽点的分离图图1 5显示了模型中毎个节点的分离图。在通常情况下，模型中的每个节点将承受 外部载荷P和内部载荷4治者是由于节点相连单元的应力引起的“对于一个处于静平 衡状态的模型，每个节点上的合力必须为零，即每 个节点上的外部载荷和内部载荷必须相互平衡。对于节点“其平衡方程可建立如下。假设圆杆的长度变化很小，则单元1的应变 可表示为其中 /利屮分别是节点尬

3、和&的位移丄是单元的初始长度。假设材料是弹性的，圆杆的应力可以通过应变乘以杨氏模量E给出：a = Em作用在端部节点的轴向力等于圆杆的应力乘以其横截面积人。这样，内力、材料性能和位移的关系为E =九 A Eg A 竽（t? ua）在节点ci处的平衡方程因此可以写成节点b的平衡必须同时考虑与该节点相连的两个单元的内力。对于节点乩单元1的 内力件用在反方向上，因此变为负值&平衡方程为Pf, 一 （U* 扯 * ） +r f 卅）R 0节点*的平衡方程为罠一竿w I?) = 0对于隐式方法这荃平衡方程需要同吋进行求解以获取每个节点的位移。以上求解 最好采用矩阵算法。因此，将内力和外力的贡献写成矩阵形

4、式。如果两个单元的性质和 维数相同，平衡方程可以化简成如下形式：通常悄况下，毎个单元的刚度（本例中的EA/L项）是不同的。因此，可以将模型中 两个单元的刚度分别写成嵐和K“我们感兴趣的是获得平衡方程的解答，使施加的外 力P与产生的内力I达到平衡。若考虑到问题的收敛性和非线性，可以将平衡方程写成0tTJ=0Ks .迟K】Ki(K+KJ对于完粮的两单元三节点结构经过这种符号变换后，可以将平衡方程更写为LPJ L 0布隐武方法中例如应用在ABAQUS/Standard中，由此方程组能够解出3个未知变撤的 值：“及（在该问题中，护已指定为0.0）。一旦位移求出后就能利用位移返何计 算出桁架单元的应力当

5、每一个求解增超步结束时，隐式的有限元方法要求解一组方 程组。（2）显式求解本节提供从概念上理解当应用显式动力学算法时应力是如何在模型中传播的U在这 牛演示例子中，考虑应力波沿着一个由3个单元构成的杆件模型传播的过程，如图1-6所 永随着时间增灵的变化，我们将研究杆件的各个状态.】 2 图1-fi自由端炸用有集中力P的杆件的初蛇状念在第1个时间增量段哀施加在节点1的集中力F使节点I产生了一个加速度亦这 个加速度使节点1产生速度d在单元内引起应变速率钵“在第1个吋间増量段内对 应变速率进行积分就获得了单元的应变增駅山日“总应变灵是初始应变號和应变增 虽的和。在这个问题中的初始应变为零。一 U计算出

6、了单元应变就可以通过材料的本 构摸型求出单元的应力兀“对于线弹性材料，应力是弹性模址与总应变的乘积。这个过 程如图1-7所示。在第1个时间增量段甲，节点2和节点3因为没有力作用在其上*所以 役有移动*左第2个时间增鼠段，由单元的应力产生的单元内力械施加到了与单元相连的所 有节点上，如图1-8所示。然后应用这些单元应力计算节点1和节点2的动力平衡 方程*jj= to罚 二击 = Eetn珂=軌w图1-7自由端作用有烷中力的杆件在第1时间皓屋段第東时的状态人吐=破1皿讣1冃竺丁巴= 山lt =心出 ii +Hrnii=*rdl =Ef4|j師-巴护令二半+L曲图1-fi 杆件在第2个时间增邀段开始

7、时的状态这牛过程继续下去，到第3个时间增虽段开始时单元和单元中都已存在了应 力，而节点1,2和3则存在了作用力，如图1码所示.这个过程将一苴继续下去哀直釧总 的分折时间结束1 234 JP / /An Ai2卩图1-9杆杵在第3牛时间蜡It段开始时的狀态隐式求解中，计算的精度完全控制于计算步数，在一般的计算软件中（flac、abaqus）,软件均是利用不平衡力来控制计算步数（当 不平衡力10-5时，停止计算）不平衡力=A+B.A表示施加在节点上的 集中力；B表示：在n步数下,根据第n步计算出来的应力，求出节点 的内力.= Lpby+丄Flac软件中B,以上公式是根据虚功原 理推倒而得到具体推倒

8、过程见flac 原理2、Ansys计算注意事项：计算单位、参数、荷载、标准值、设计值，计算过程中系数的加入.（1） beam单元对于beam单元.Ansys软件中我们常用的有两种梁单元：beam188和beam4.这两种单元均是三维的梁单元，每个节点都具有6个自由度（ux、uy、uz、mx my mz）,并且单元坐标系 x轴是i点指向j点.Beam188单元是基于铁木辛哥理论的梁,beam4单元是我们常用的 经典结构力学梁.（铁木辛哥理论考虑了梁的剪切变形，而我们常用的 经典结构力学梁只考虑了弯矩对结构的变形影响）所以说 ,beam188 可以更精确的计算梁单元，因此我们结构计算中，一般都采用

9、beam188 单元.当然还有beam189单元,189单元属于三维二次的梁单元（beam188属于三维一次梁单元），精度比beam188更加高.定义beam188单兀， !定义单元/prep7et,1,beam188般米用如卜形式：!进入前处理!定义单元1 ；88号标号为1!定义材料属性mp,ex,1,2.55e7!定义弹性模量（kn/m 2）mp, nu xy,1,0.167!定义泊松比mp,de ns,1,2.5!定义密度（KN/N*KG/M)nummrg,all!合并重合节点nu mcmp,all!压缩编号!定义梁截面SECTYPE,1,BEAM,RECT,A1,0! 1表示梁编号；R

10、ECT表示是矩形梁（还有其他t型等等,具体见ansys帮助）；A1表示梁的名称；0表示薄壁梁单元网格划分精细程度（0 5).SECDATA,1,3,4,12!1表示梁b ; 3 表示梁h ; 4和12定义对应宽长等分份数.SECOFFSET,CENT!cent 质心；shrc剪切中心；origin原始中心;user用户定义;!注意：当梁单元和壳单元一起使用时 ，可以设置梁单元的偏心，使梁的一面和壳的一面共面(secoffset,user,offsety,offset z),如下图:!划分网格LSEL,S,1LATT,1,1,1,1LESIZE,ALL,0.2,1LMESH,ALLnCEHTXU

11、SEROFFSETYj OFFSET-ZASHRCORIGIN!选中编号为1的线.!mp,r,et,方向点 ”SECTYPE截面号.!0.2是单元大小,1是确认细分规则. !用beam单元离散模型，形成网格.!对于划分网格，空间的beam单元，由于需要确定b、h的方向，ansys软件利用方向点来控制b、h的方向.方向点的编号最好定义的很大，如果定义太小，会影响后面的加载具体方向点如何控制见上面的latt命令和ansys帮助.自己试两下就知道怎么用了Alls FINISH !加载加约束/SOLUACEL,9.8!重力加速度.注意方向，数值和整体坐标相反，比如重力指向z轴 负向，则为正值.SFCU

12、M,ALL,ADD!设置单兀荷载是叠加还是替代 ，只对加在单兀和节点上的荷载有 效,对于加在面、线上的荷载，都只有替代作用（对冋一个面 ，第 二次加的荷载替代第一次加的荷载）!对于beam单元,只能根据sfbeam命令增加均布荷载 等大小的均布荷载Lsel,s,1ESLL,S,1sfbeam,ALL,1,PRES,-161.5 !1 表示作用在beam单元的面上(如下图,面表示beam 单元的轴向，面表示单元侧面，面表示beam单元顶 面),-161.5表示均布荷载大小，正负号可以控制作用力的 方向Z梯形均布荷载Sfbeam命令是对每个单元进行加载 如果一根梁承受10100的梯形均布荷载，而且

13、这根梁 被分成了 10个beam单元，这样施加荷载就非常困难因此我将这种加载过程写成命令流，让软件自动进行加载ESLL,S,1命令流如下：beam单元!选中要加载的那根梁(线)!选中属于这根梁(线)的*GET,Nelem,ELEM,COUNT, , , ,!获得当前所选单元个数，赋予参数Nelem*GET,Ne,ELEM,NUM,MIN, , , ,!获得当前所选单元最小编号 ，赋予参数Ne*DO,l,1,Nelem!循环加载，循环次数=单元个数ESEL,S,NeNSLE,S,1*GET,Nnode,NODE,COUNT, , , ,!获得当前所选节点个数，赋予参数Nnode*GET,Nn,N

14、ODE,NUM,MIN, , , ,!获得当前所选节点最小编号，赋予参数NnNN1X=NX(NN)!将nn节点的x坐标赋予 NN1X(NX表示x坐标,NY表示y坐标)NN=NDNEXT(NN)!NN=当前所选节点的下一个编号NN2X=NX(NN)!将nn节点的x坐标赋予 NN2Xsfbeam,ALL,1,PRES,-1630.76/3.23*(3.23-NN1X),-1630.76/3.23*(3.23-NN2X)!以上荷载公式应根据实际情况进行调整ESLL,S,1NE=EL next(NE)!NE=当前所选单元的下一个编号*ENDDO!对于此命令流，根据不同的实际情况，ABC部分需要修改，其

15、他不需要修改!后处理（XY平面）（大拇指指向y，就是my）etable,lmY,smisc,2!显示弯距etable,JmY,smisc,15PLLS,IMY,JMYETABLE,IFX,SMISC,1!显示轴力ETABLE,JFX,SMISC,14PLLS,IFX,JFXETABLE,IFY,SMISC,5 !显示剪力ETABLE,JFY,SMISC,18PLLS,IFY,JFY!注意：beam单元的结果输出都是以单元坐标系输出的，且拉为正、压为负前面我们已经知道，单元坐标系x轴就是i点指向j点,其他坐标可以根据整体坐标系推出.详细内容见ansys帮助（2） Shell 单元对于shell单

16、元应用的范围，ansys软件并没有强制规定，只是从 字面上区分了薄壳和厚壳我以前看过一本电子教案仿真在线， 里面说一般规定壳体的主尺寸是厚度的10倍左右,都是可以用壳体来模拟的 一般高度与跨度之比（非与单元尺寸比较）1/15 & 1/10,可以当作厚壳来处理shell63是薄壳单元，他包含弯曲和薄膜效应，但是忽略横向剪切变 形；shell43,shell143,shell181,shell91,shell93和 shell99,都属于厚壳单元， 不仅有弯曲、薄膜效应，他也包含了横向剪切效应横向剪切被表示为 整个厚度上的常剪切应变这种一阶近似只适用于中等厚度壳体线形分析时，如果不包含横向剪切应变

17、，使用63,163单元;如果横 向剪切变形重要，则遵守以下原则：均匀材料使用43,93,143单元复合 材料使用91,99,181.我们土木工程中，一般利用shell43计算.!定义单元!进入前处理!定义单元43号标号为1!定义弹性模量（kn/m 2）!定义泊松比!定义密度（KN/N* *KG/M）/prep7 et,1,shell43!定义材料属性mp,ex,1,2.55e7mp, nu xy,1,0.167mp,de ns,1,2.5!定义墙体厚度!等厚度板R,1,2R,2,3Asel,s,1Aatt,1,1,1ESIZE,0.2MSHAPE,0,2D!1表示编号,2表示厚度（m!选中1号

18、面!mp,real,type!疋义单兀大小为 0.2左右!规定划分单元形状,0表示四边形（1表示三角形），2d表示划分面（3d 表示划分体MSHKEY,2!指定是自由划分还是映射划分,2表示：尽量用映射划分，不符合要求就自动使用自由划分，具体参见ansys帮助的eshkey命令.AMESH,ALL!划分面单元.注意：在网格剖面方面，最好全部用四边形，而且形状尽量规则、均匀！因为将来后处理内力提取的时候，提取出来的力和单元的大小有直接的关系。后处理部分将详细讲解。!渐变厚度的板对于渐变的板，首先选中要划分的面，然后进行网格划分但是在网格划分的 时候，赋予的厚度编号开始选用100 （前面没有定义的

19、就可以），然后利用以下命 令进行厚度赋值asel,s,1Asel,s,1Esla,s,1Local,11,0,1,0,0Emodif,all,esys,11 注意：修改单元坐标系的时候 坐标系修改单元坐标系!划分好网格后，有一个很重要的事情：重新调整单元坐标系，因为好多单元的单 元坐标系是乱七八糟的，但是将来结果文件输出值的方向是单元坐标系的因此 必须把所有单元的单元坐标系调整的和整体坐标系一样（只是方向一样，不一定 x、y、z对应,因为壳体的z轴永远是壳体的法向）.命令如下：!选中1号面!选中属于1号面的所有单元!11局部坐标系编号；0表示笛卡尔坐标系；1表示x值!修改单元坐标系为局部坐标系

20、11,必须先定义一个局部坐标系因为emodif命令只能利用局部!将壳体法向反方向：有时候需要调整壳体法向的时候，可以利用如下命令：Asel,s,1AREVERSE,ALL !将面的方向反过来，ansys会提示你：是否也要将属于面的单元法向 也反过来？AllsFINISH!加载加约束/SOLUACEL,9.8!施加重力加速度，详细解释见上.等大小的均布荷载对壳体施加等大小的均布荷载比较简单，以下命令均可以实现，只要你选中相应的面、点或者单元.Sfa,ALL,1,PRES,100Sf,ALL,PRES,100!个人不推荐此命令，因为此命令对施加在实体面上的节点有限制.Sfe,ALL,1,PRES,

21、0,100 !主力推荐此命令.1表示作用壳体的哪个面（如下图所示）;0表示确实存在的压力（SURF153 or SURF154单元选 用2）;100表示压力（kn/m ）.梯形荷载SFGRAD,PRES,O,Z,5,-10*1.2!为下面的加载形成梯度。pres表示施加压力；0表示坐标系代号（笛卡尔坐标系）;z表示梯度沿z方向；5 表示：在0坐标系下z= 5的位置（下面命令的300就加 在这个位置）;-10*1.2表示施加的梯度，负号表示随 着z的增大，压力在减小.SFe,ALL,2,PRES,0,300!见上 sfe 命令解释.!后处理（如下图）!壳体单元输出的是每个单元在单元坐标系下，单位

22、长度的值。（这句话是ansys帮助里翻译而来的）有点含糊哈，大家凑合理解。我举个列子给你们听听。 如下例题：计算图选取单元图求基础和板之间的剪力？用屁股想都想的出来，精确解是100k n。如此简单的问题，我们选用shell63做计算（当然我们也可以选用平面单元， 这里只是为了试验shell63求解结果的含义）对这个模型，我们分两种情况计算，两种情况唯一不同之处就是网格大小， 第一种情况网格大小1*1,第二种情况网格大小0.5*0.5。取出靠近地面的最近的一排单元（如上图），利用单元table输出每个单元 的剪力值（单元坐标系下），然后相加。第一种情况，相加后 A= 99.8 kn第二种情况，相

23、加后 A= 198.9 kn由上可得：剪力值=A*单元大小这里可以看出，划分单元尽量规则、均匀的重要性，特别是将来你要提取内力的位置， 在划分网格的时候就要做好准备！而且计算前单元坐标系要统一也非常重要！本人编写了一个提取内力的通用程序，各位可以参考参考。只要选中你要提取内力的那排单元（如同上图的那排），然后输入以下命令即可：JLzo nghe=0 ! 剪力ZLZONGHEx=0 !X 方向轴力 ZLZONGHEy=0 !Y 方向轴力 mzo nghex=0 !X 方向弯矩 mzo nghey=O !Y 方向弯矩 etab,JL,smisc,3 ETAB,ZL2,SMISC,2 ETAB,ZL

24、1,SMISC,1 etab,m1,smisc,4etab,m2,smisc,5*GET,MAXN,ELEM,O,NUM,MAX !所选单元的最大编号给max n*GET,MINN,ELEM,O,NUM,MIN*DO,I,MINN,MAXN*IF,ESEL(I),EQ,1,THEN !判断单元号I是否属于所选范围*get,jia nli,elem,i,etab,JL*get,ZHOULIx,elem,i,etab,ZL1 *get,ZHOULIy,elem,i,etab,ZL2*get,wanjux,elem,i,etab,m1*get,wa nju y,elem,i,etab,m2JLz o

25、n ghe=JLz on ghe+jia nliZLzon ghex=ZLz on ghex+ZHOULIxZLzon ghey=ZLz on ghey+ZHOULIymzon ghex=mz on ghex+wa njuxmzon ghey=mz on ghey+wa njuy*ENDIF*ENDDO注意：(1)以上命令求解出来的 Jlz on ghe 等等还没有乘以单元大小。(2)以上求解出来的X方向、Y方向的轴力和弯矩都是单元坐标系下的值(具体方向 如下图，SX、SY就是单元坐标系)，单元坐标系的 X Y方向可以通过视图显示出来，命令如 下：/psymb,esys,1蓝线表示Z轴，白线表

26、示 X轴，青线表示 Y轴 SXTOP) SX MlDl SX(BOD直M = Element x-axis if ESYS is not supplied.x = Element x-axis if ESYS is supplie d.（3）Solid 单元实体单元其实是ansys最简单的单元，但是也是ansys最难操作的单元。 实体单元从建模到形成计算书，70%的时间用在划分网格上，20%的时间用 在后处理上，还有10%就是一些小杂事了。前期操作方面，有两种方法可以完成：（1）利用ansys建模同时划分网格。我有一个很牛的同事，做窝壳方面很有研究，他曾经做过一个工程的窝壳， 写了整整2万行的

27、命令流，那窝壳真是漂亮，不得不让人佩服。不过这种方法需 要非常非常熟悉ansys,不是一般人可以的哦（汗颜！我为什么不行，55555）。（2）利用三维软件建模（个人推荐solidwork,和ansys的兼容性非常ok）， 然后导入ansys。这种方法比较通用，但是一般solidwork导入的都是一个整体块，需要人工 不停的用wp （工作面）切割，切到可以映射划分为止。这个切割技术还是比较 难的，需要自己多去试试。一定要自己动手去切割、去划分，只要认真做两个例 题，你就能自学成才了。上面废话一堆，实体单元的前处理我就不说了，重点说一下我们土建专业比 较关注的后处理：内力（剪力、轴力、弯矩）提取。

28、由于solid单元最后输出的都是应力，无法在工程中利用（当然也可以根据规 范提到的应力配筋， 但个人感觉不太稳妥，我曾经比较过：同一截面，应力配筋和提取内力的结果，内力配筋结果都大于应力配筋），但是我们可以利用这些应力积分得到我们 所需要的内力。积分的原理就不说了，哪个兄弟要是忘了就去再看看结构力学书 吧。首先将工作平面的xy面移到你要计算内力的位置，然后你要判断:(1) 你准备计算内力的截面属于整体坐标系的 xy平面SUCR,SUZ2,CPLANE,3 ! 移动工作平面、创建面 SUZ2 SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射 X和 Y方向应力 SUMAP,MYS

29、Z,S,Z$SUMAP,MYSXY,S,XY 映射 Z 和 XY方向应力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ 映射 YZ和 XZ 方向应力 SUEVAL,XFORCE,MYSXZ,INTG 求截面上 FX SUEVAL,YFORCE,MYSYZ,INTG 求截面上 FY SUEVAL,ZFORCE,MYSZ,INTG !求截面上 FZSUEVAL,MYA,DA,SUM !求截面面积并赋给 MYA变量 SUEVAL,MYYA,GCY,INTG !求关于X轴的面积矩并赋给变量MYYAMYYA=MYYA/MYA !得到面积重心到 X轴的距离=面积矩/面积 SUEVAL

30、,MYXA,GCX,INTG !求关于Y轴的面积矩并赋给变量 MYXA MYXA=MYXA/MYA !得到面积重心到 Y轴的距离=面积矩/面积 SUCALC,SZGCY,MYSZ,MULT,GCY计算 MYSZ GCY 并赋给面项 SZGCY SUEVAL,MX1,SZGCY,INTG !对面项 SZGCY在面上积分得到 MX1 SUCALC,SZGCX,MYSZ,MULT,GCX计算 MYSZ GCZ 并赋给面项 SZGCX SUEVAL,MY1,SZGCX,INTG !对面项SZGC淮面上积分得到 MY1!上述弯矩基于总体直角坐标系原点而言的，应对面积重心取矩，将内力简化到面积重心上 MX

31、1=MX1-ZFORCE*MYYA !MY1=MY1-ZFORCE*MYXA !(2) 你准备计算内力的截面属于整体坐标系的yz平面SUCR,SUZ2,CPLANE,3 ! 移动工作平面、创建面 SUZ2 SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射 X和 Y方向应力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY 映射 Z 和 XY方向应力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ 映射 YZ和 XZ 方向应力 SUEVAL,XFORCE,MYSx,INTG 求截面上 FX SUEVAL,zFORCE,MYSXZ,INTG 求截

32、面上 Fz SUEVAL,yFORCE,MYSXY,INTG !求截面上 FySUEVAL,MYA,DA,SUM !求截面面积并赋给 MYA变量 SUEVAL,MYZA,GCZ,INTG!求关于X轴的面积矩并赋给变量MYYAMYZA=MYZA/MYA !得到面积重心到 Y轴的距离=面积矩/面积 SUEVAL,MYYA,GCY,INTG !求关于Z轴的面积矩并赋给变量 MYYA MYYA=MYYA/MYA !得到面积重心到 Z轴的距离=面积矩/面积 SUCALC,SZGCZ,MYSX,MULT,GCZ 计算 MYSX GCZ 并赋给面项 SZGCZ SUEVAL,MY1,SZGCZ,INTG !

33、对面项 SZGCZ面上积分得到 MY1 SUCALC,SZGCY,MYSX,MULT,GC丫计算 MYSX GCY 并赋给面项 SZGCY SUEVAL,MZ1,SZGCY,INTG !对面项 SZGCY在面上积分得到 MZ1!上述弯矩基于总体直角坐标系原点而言的，应对面积重心取矩，将内力简化到面积重心上 MY1=MY1-XFORCE*MYZA !MZ1=MZ1-XFORCE*MYYA !(3) 你准备计算内力的截面属于整体坐标系的xz平面SUCR,SUZ2,CPLANE,3 ! 移动工作平面、创建面 SUZ2 SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y !映射 X和 Y方向应

34、力 SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY 映射 Z 和 XY方向应力 SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ 映射 YZ和 XZ 方向应力 SUEVAL,XFORCE,MYSxy,INTG 求截面上 FX SUEVAL,zFORCE,MYSYZ,INTG 求截面上 Fz SUEVAL,yFORCE,MYSy,INTG !求截面上 Fy SUEVAL,MYA,DA,SUM !求截面面积并赋给 MYA变量 SUEVAL,MYZA,GCZ,INTG !求关于X轴的面积矩并赋给变量 MYYA MYZA=MYZA/MYA !得到面积重心到 X轴的距离

35、=面积矩/面积 SUEVAL,MYXA,GCX,INTG !求关于Z轴的面积矩并赋给变量 MYXA MYXA=MYXA/MYA !得到面积重心到 Z轴的距离=面积矩/面积 SUCALC,SZGCZ,MYSy,MULT,GCZ 计算 MYSY GCZ 并赋给面项 SZGCZ SUEVAL,MX1,SZGCZ,INTG !对面项 SZGCZ面上积分得到 MX1 SUCALC,SZGCX,MYSy,MULT,GCX 计算 MYSY GCX 并赋给面项 SZGCX SUEVAL,MZ1,SZGCX,INTG !对面项 SZGC淮面上积分得到 MY1!上述弯矩基于总体直角坐标系原点而言的，应对面积重心取

36、矩，将内力简化到面积重心上 MX1=MX1-yFORCE*MYZA !MZ1=MZ1-yFORCE*MYXA !(1) XFORCE YF0RGEZF0RCEMX1 MY1 MZ1为最后的计算结果，而且是在 整体坐标系下的结果。(2) 对于xy平面，ZFORC为该断面的轴力；XFORCEYFORC分别表示x、y 方向的剪力；MX1 MY1分别表示右手准则，大拇指分别指向整体坐标系 下的x、y轴。其他yz、xz平面类似。(3) |这里强调一下求得内力的方向问题。 我们求解内力的时候，其实就是对各应力进行面积积分，所得内力正负值的意义和各应力一样。在ansys中各 应力的正负的意义是和有限元定义正负号的方法是一致的， 和单纯的整体 坐标系的正负号是不一样的。举例如下：比如你求一个属于xy平面的截面的内力，求出剪力=100kn。那么对于 这个剪力方向是什么呢？根据有限元正负号定义原则！可参照有限元原理(4) 以前我一个同事发现了一个命令，很简单的一个命令输入，然后就可以输 出以上结果。但是对于此命令得到结果，本人不太明白他输出的意思，而 且你所选的范围不一样，结果也不一样，所以我就不推荐大家了。（4）接触单元三维接触问题个人理解，我们平常用的最多的就是面-面接触和mpc装配。