ANSYS中杆梁壳单元内力的输出
ANSYS提取单元内力的方法
ANSYS提取单元内力的方法ANSYS提取单元或节点内力的方法方法1：节点荷载（List Results→Nodal Loads）方法2：节点合力计算（Nodal Cals→Sum @ Each Node）方法3：单元解中的节点解（List Results→Element Solution→Structural Forces & Moments）方法4：支座反力（List Results→Reaction Solu）方法5：单元表（List Results→Elem Table Data）上述各方法提取的结果关系如下：（1）方法1和方法2提取的结果完全相同，但结果为0的项在方法1的结果列表中不显示，而方法2的结果列表则会全部显示。
（2）方法3提取的结果是每个单元各节点在该单元中的内力，针对同一节点，将其在各个单元中的内力求和，其累加结果与方法1和2得到的结果一致。
（3）方法4提取的结果只显示有施加位移约束的节点反力，其数值大小与方法1和2得到的结果相差一个正负号，即节点内力和节点反力刚好是一对作用力与反作用力。
（4）方法5提取的结果是单元的内力，如果单元的形函数为线性（如BEAM188单元设置“KEYOPT(3)=0”），则ANSYS会取单元中点作为积分点并将其数值代替单元内的线性变化，因此其输出结果的绝对值等于方法3中对应单元的各节点相应内力绝对值的平均值；如果单元的形函数为非线性（如BEAM188单元设置“KEYOPT(3)=2”），则单元各节点的内力不同，其结果与方法3得到的结果一致。
（5）方法1~4提取的结果都是默认基于整体坐标系的，而方法5提取的结果是基于单元坐标系的，因此提取结果的方向和正负号需特别注意。
有限元中力的方向和结构力学中的方向是有区别的，不论是什么结果坐标系，力的正方向取为对应结果坐标的正方向，弯矩则是对应坐标轴的顺时针为正。
几种梁单元用法与结果——ansys命令流
!EX6.3 几种梁单元用法与结果!EX6.3A-BEAM4单元finish/clear/prep7et,1,beam4mp,ex,1,2.1e11mp,prxy,1,0.3r,1,1032e-5,158936e-9,1947744e-11,0.18,0.3 rmore,,110976e-11r,2,1032e-5,1947744e-11,158936e-9,0.3,0.18 rmore,,110976e-11k,1k,2,,4k,3,3,4k,4,3,4,-2l,1,2l,2,3l,3,4dk,1,alldk,4,ux,,,,uy,uzfk,3,fy,-15000fk,3,fz,8000lesize,all,,,10lsel,s,,,1latt,1,1,1lmesh,allesll,ssfbeam,all,2,pres,3000lsel,s,,,2latt,1,1,1lmesh,allesll,ssfbeam,all,2,pres,2000lsel,s,,,3latt,1,2,1lmesh,allesll,ssfbeam,all,1,pres,-1000 allsel/solusolve/post1pldisp,1etable,fxi,smisc,1 etable,fxj,smisc,7 etable,fyi,smisc,2 etable,fyj,smisc,8 etable,fzi,smisc,3 etable,fzj,smisc,9 etable,mxi,smisc,4 etable,mxj,smisc,10 etable,myi,smisc,5 etable,myj,smisc,11 etable,mzi,smisc,6 etable,mzj,smisc,12 etable,smini,nmisc,2 etable,sminj,nmisc,4 !单元坐标系中的结果plls,fxi,fxjplls,fyi,fyjplls,fzi,fzjplls,mxi,mxjplls,myi,myjplls,mzi,mzjplls,smini,sminj prnsol,uprrsol!------------------------------!EX6.3B-BEAM24单元finish/clear/prep7et,1,beam24,1mp,ex,1,2.1e11mp,prxy,1,0.3r,1,-0.09,-0.14,0,0.09,-0.14,0.02 rmore,0,-0.14,0,0,0.14,0.012 rmore,-0.09,0.14,0,0.09,0.14,0.02 k,1k,2,,4k,3,3,4k,4,3,4,-2k,51,-3,2k,52,2,6k,53,0,4,-1l,1,2l,2,3l,3,4dk,1,alldk,4,ux,,,,uy,uzfk,3,fy,-15000fk,3,fz,8000lesize,all,,,10lsel,s,,,1latt,1,1,1,,,51lmesh,allesll,ssfbeam,all,1,pres,3000lsel,s,,,2latt,1,1,1,,,52lmesh,allesll,ssfbeam,all,1,pres,2000lsel,s,,,3latt,1,1,1,,,53lmesh,allesll,ssfbeam,all,1,pres,1000allsel/solusolve/post1pldisp,1prnsol,uprrsoletable,smini,nmisc,2etable,sminj,nmisc,4plls,smini,sminj!------------------------------!EX6.3C-BEAM44单元finish/clear/prep7et,1,beam44mp,ex,1,2.1e11mp,prxy,1,0.3sectype,1,beam,isecoffset,centsecdata,0.18,0.18,0.3,0.02,0.02,0.012 k,1k,2,,4k,3,3,4k,4,3,4,-2k,51,-3,2k,52,2,6k,53,0,4,-1l,1,2l,2,3l,3,4dk,1,alldk,4,ux,,,,uy,uzfk,3,fy,-15000fk,3,fz,8000lesize,all,,,10lsel,s,,,1latt,1,,1,,,51,1 lmesh,allsfbeam,all,1,pres,3000 lsel,s,,,2latt,1,,1,,,52,1 lmesh,allesllsfbeam,all,1,pres,2000 lsel,s,,,3latt,1,,1,,,53,1 lmesh,allesllsfbeam,all,1,pres,1000 allsel/solusolve/post1pldisp,1etable,smini,nmisc,2 etable,sminj,nmisc,4 plls,smini,sminj prnsol,uprrsol!------------------------------!EX6.3D-BEAM188单元finish/clear/prep7et,1,beam188mp,ex,1,2.1e11mp,prxy,1,0.3sectype,1,beam,isecoffset,centsecdata,0.18,0.18,0.3,0.02,0.02,0.012 k,1k,2,,4k,3,3,4k,4,3,4,-2k,51,-3,2k,52,2,6k,53,0,4,-1l,1,2l,2,3l,3,4dk,1,alldk,4,ux,,,,uy,uzfk,3,fy,-15000fk,3,fz,8000lesize,all,,,10lsel,s,,,1latt,1,,1,,,51,1lmesh,allsfbeam,all,1,pres,3000lsel,s,,,2latt,1,,1,,,52,1lmesh,allesllsfbeam,all,1,pres,2000lsel,s,,,3latt,1,,1,,,53,1lmesh,allesllsfbeam,all,1,pres,1000allsel/solusolve/post1pldisp,1prnsol,uprrsol/eshape,1plnsol,s,x!------------------------------!EX6.3E-BEAM189单元finish/clear/prep7et,1,beam189mp,ex,1,2.1e11mp,prxy,1,0.3sectype,1,beam,isecoffset,centsecdata,0.18,0.18,0.3,0.02,0.02,0.012 k,1k,2,,4k,3,3,4k,4,3,4,-2k,51,-3,2k,52,2,6k,53,0,4,-1l,1,2l,2,3l,3,4dk,1,alldk,4,ux,,,,uy,uzfk,3,fy,-15000fk,3,fz,8000lesize,all,,,10lsel,s,,,1latt,1,,1,,,51,1 lmesh,allsfbeam,all,1,pres,3000 lsel,s,,,2latt,1,,1,,,52,1 lmesh,allesllsfbeam,all,1,pres,2000 lsel,s,,,3latt,1,,1,,,53,1 lmesh,allesllsfbeam,all,1,pres,1000 allsel/solusolve/post1pldisp,1prnsol,uprrsol/eshape,1plnsol,s,x。
ANSYS梁壳单元及其应用
ANSYS梁壳单元及其应用一梁单元特性及应用梁单元用来建立三维结构的数学上理想化的一维有限元模型，与实体单元和壳单元相比，使用梁单元可以提高求解效率。
ANSYS提供了多种梁单元库以适应不同的需要，其中Beam44为3-D 渐变非对称截面梁，Beam188和Beam189为3-D有限应变梁，ANSYS的梁单元在非线性分析方面具有先进性和鲁棒性的独特优势。
1．3D真实描述梁单元在空间上是一维的线单元，单元特性和截面属性是相互独立的，通过指定截面编号，一维的梁单元就可以描述真实的三维空间结构，并且ANSYS可以以三维的形式显示分析结果。
ANSYS提供了11种常用的截面形状，并允许用户定制截面形状，用户可以利用二维建模的方式创建新截面，并可以把定制的截面形状保存在截面形状库中。
2．变截面梁ANSYS允许定义任意截面形状，允许单元的每一端具有不同的不对称几何形状，并允许其端节点从梁的中心轴偏移。
3．梁单元的预应力ANSYS的梁单元可以考虑预应力产生的应力刚化效应。
所谓应力刚化效应对于梁单元来说就是轴向应力引起的垂直轴向的刚度变化。
4．复合材料截面ANSYS可以定义任意几何形状由多个各向同性材料组成的横截面，可以用来模拟层状复合材料梁，长纤维增强复合材料梁和传感器等。
（图2）5．考虑剪切变形和翘曲的影响ANSYS的梁单元基于铁木辛柯梁理论，在平面假设的基础上可以考虑剪切变形的影响。
ANSYS的梁单元还可以考虑非圆截面梁扭转时产生的翘曲影响，这时每个端节点有7个自由度，包括3个平动，3个转动和一个翘曲自由度。
Shell188和Shell189单元不仅能模拟直梁的弯曲剪切响应，而且能模拟横向—扭转屈曲行为（特征屈曲和非线性崩塌）（图3）。
6．支持非线性材料本构模型ANSYS的梁单元支持弹性（线弹性和非线性弹性），塑性和蠕变等材料本构模型。
二壳单元特性及应用壳单元用来创建三维结构的二维理想化模型，用壳单元模拟壳状结构比用实体单元具有更高的计算效率。
ANSYS中不同单元之间的连接问题总结
MP,PRXY,1,0.3 !定义泊松比
R,1 !定义实体单元实常数
R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数
BLC4,,,20,7,10 !创建矩形块为实体模型
WPOFFS,0,3.5 !将工作平面向Y方向移动3.5
ANSYS中不同单元之间的连接问题
总结一下不同单元之间的连接问题
论坛里常有人问不同单元之间的连接问题，我自己也一直被这个问题所困绕，最近从ANSYS工程分析进阶实例上知道了ANSYS中不同单元之间的连接原则。感觉收收获不小，现把它上传与大家共享。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约事方程。例如：
MPC即Multipoint Constraint,多点约束方程，其原理与前面所说的方程的技术几乎一致，将不连续、自由度不协调的单元网格连接起来，不需要连接边界上的节点完全一一对应。
MPC能够连接的模型一般有以下几种。
solid模型-solid模型
shell模型-shell模型
solid模型-shell模型
单元类型的选择问题
初学ANSYS的人，通常会被ANSYS所提供的众多纷繁复杂的单元类型弄花了眼，如何选择正确的单元类型，也是新手学习时很头疼的问题。
单元类型的选择，跟你要解决的问题本身密切相关。在选择单元类型前，首先你要对问题本身有非常明确的认识，然后，对于每一种单元类型，每个节点有多少个自由度，它包含哪些特性，能够在哪些条件下使用，在ANSYS的帮助文档中都有非常详细的描述，要结合自己的问题，对照帮助文档里面的单元描述来选择恰当的单元类型。
总结一下不同单元之间的连接问题
总结一下不同单元之间的连接问题论坛里常有人问不同单元之间的连接问题，我自己也一直被这个问题所困绕，最近从ANSYS 工程分析进阶实例上知道了ANSYS中不同单元之间的连接原则。
感觉收收获不小，现把它上传与大家共享。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约事方程。
例如：（1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。
（2）梁与壳有公共节点怒可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
（3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
（4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
例如：杆与梁、壳、体单元有公共节点即可，不需要约束方程。
梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
举例：有一长为100mm的矩形截面梁，截面为10X1mm，与一规格为20mmX7mmX10mm的实体连接，约束实体的端面，在梁端施加大小为3N的y方向的压力，梁与实体都为一材料，弹性模量为30Gpa，泊松比为0.3。
本例主要讲解梁与实体连接处如何利用耦合及约束方程进行处理。
命令流如下：FINI/CLE/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名/PREP7ET,1,SOLID95 !定义实体单元类型为SOLID95ET,2,BEAM4 !定义梁单元类型为BEAM4MP,EX,1,3E4 !定义材料的弹性模量MP,PRXY,1,0.3 !定义泊松比R,1 !定义实体单元实常数R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数BLC4,,,20,7,10 !创建矩形块为实体模型WPOFFS,0,3.5 !将工作平面向Y方向移动3.5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPOFFS,0,5 !将工作平面向Y方向移动5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPCSYS,-1 !将工作平面设为与总体笛卡儿坐标一致K,100,20,3.5,5 !创建关键点K,101,120,3.5,5 !创建关键点L,100,101 !连接关键点生成梁的线实体LSEL,S,LOC,X,21,130 !选择梁线LATT,1,2,2 !指定梁的单元属性LESIZE,ALL,,,10 !指定梁上的单元份数LMESH,ALL !划分梁单元VSEL,ALL !选择所有实体VATT,1,1,1 !设置实体的单元属性ESIZE,1 !指定实体单元尺寸MSHAPE,0,2D !设置实体单元为2DMSHKEY,1 !设置为映射网格划分方法VMESH,ALL !划分实体单元ALLS !全选FINI !退出前处理/SOLU !进入求解器ASEL,S,LOC,X,0 !选择实体的端面DA,ALL,ALL !约束实体端面ALLS !全选FK,101,FY,-3.0 !在两端施加Y向压力CP,1,UX,1,21 !耦合节点1和节点21X方向自由度CP,2,UY,1,21 !耦合节点1和节点21Y方向自由度CP,3,UZ,1,21 !耦合节点1和节点21Z方向自由度CE,1,0,626,UX,1,2328,UX,-1,1,ROTY,-ABS(NZ(626)-NZ(2328)) !设置约束方程CE,2,0,67,UX,1,4283,UX,-1,1,ROTZ,-ABS(NY(67)-NY(4283)) !设置约束方程CE,3,0,67,UZ,1,4283,UZ,-1,1,ROTX,-ABS(NY(67)-NY(4283)) !设置约束方程ALLS !全选SOLVE !保存FINI !退出求解器/POST1 !进入通用后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, S,EQV, 0,1.0 !显示等效应力ETABLE,ZL1,SMISC,1 !读取梁单元上I节点X方向的力ETABLE,ZL2,SMISC,7 !读取梁单元上J节点X方向的力ETABLE,MZ1,SMISC,6 !读取梁单元上I节点Z方向的力矩ETABLE,MZ2,SMISC,12 !读取梁单元上J节点Z方向的力矩PLETAB,ZL1 !显示梁单元X方向的力PLETAB,MZ1 !显示梁单元Z方向力矩上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的，有一定的局限性，只适用于小位移，下面介绍一种支持大位移算法的方法，MPC法。
ANSYS中杆梁壳单元内力的输出
杆梁壳单元内力的输出计算完成后，就可以对杆梁壳单元进行后处理，显示位移、应力（需打开单元开关显示开关）的操作和普通的单元的后处理是一样的。
但其内力的输出和显示需要进行特殊的处理。
首先要保证在定义单元类型时，打开了内力输出的选项；其次，在完成计算后需要通过定义单元表的方法来提取。
一、定义输出选项Beam4单元说明的input summary中可以看到关于内力输出的开关选项说明：二、定义单元表查看对应项和序号：在单元说的后的单元输出表中查到MFOR(X,Y,Z)和MMOM(X,Y,Z)对应的项和序列号如下表，注意I和J节点上对应的序列号不同。
GUI方式如下图：命令流过程如下：（I节点处结果）Etable,forx_i,smisc,1 ！x方向剪力（即轴力）Etable,fory_i,smisc,2 ！y方向剪力Etable,forz_i,smisc,3 ！z方向剪力Etable,momx_i,smisc,4 ！x方向弯矩Etable,momy_i,smisc,5 ！y方向弯矩Etable,momz-i,smisc,6 ！z方向弯矩（J节点处结果）Etable,forx_j,smisc,7 ！x方向剪力（即轴力）Etable,fory_j,smisc,8 ！y方向剪力Etable,forz_j,smisc,9 ！z方向剪力Etable,momx_j,smisc,10 ！x方向弯矩Etable,momy_j,smisc,11 ！y方向弯矩Etable,momz-j,smisc,12 ！z方向弯矩三、绘制内力图对线单元和面单元，其绘制单元表的结果的命令是不一样的。
线单元的内力云图的绘制为：PLLS. 面单元的内力云图的绘制为：PLETAB四、注意事项1.弯矩图的绘制与显示比例在绘制弯矩图时（其他内力图也一样），可以给定一个显示比例系数Fact（默认等于1），当为负值时，反向显示，数值为放大系数。
同时还须注意的是，弯矩图的显示还与单元坐标系的方向有关，因此，如果同一段梁中各个单元的单元坐标系不一致（可以通过改变线的方向来改变线单元的单元坐标系），还将导致弯矩图显示不连续。
ANSYS中SHELL181单元理解和参数详解
ANSYS中SHELL181单元参数详解SHELL181单元说明:SHELL181单元适合对薄的到具有一定厚度的壳体结构进行分析。
它是一个4结点单元，每个结点具有6个自由度：x,y,z方向的位移自由度和绕X,Y,Z轴的转动自由度。
(如果应用了薄膜选项的话，那该单元则只有移动自由度了)。
简并三角形选项只在该单元做为充填单元进行网格划分时才会用到。
Shell181单元非常适用于分析线性的，大转动变形和非线性的大形变。
壳体厚度的变化是为了适应非线性分析。
在该单元的应用范围内，完全积分和降阶积分都是适用的。
SHELL181单元阐明了以下（荷载刚度）分布压强的效果。
SHELL181单元可以应用在多层结构的材料，如复合层压壳体或者夹层结构的建模。
在复合壳体的建模过程中，其精确度取决于第一剪切形变理论（通常指明德林-雷斯那壳体理论）在解决许多有收敛困难的问题上，SHELL181单元可以用来替代SHELL43单元。
参见ANSYS理论参考中的SHELL181单元介绍以了解有关该单元的更多细节问题。
SHELL181 输入数据其几何特性，节点的位置及坐标系参见图181.1 SHELL 181单元的几何图。
该单元由四个节点定义而成：I，J，K，L。
该单元等式是基于对数应变和实际应力的测量的。
从运动学来讲，该单元支持有限的薄膜应变（伸展）。
但是，在一定时间内的曲线变化必须是小变化。
你可以用实常数或者横截面定义来定义该单元的厚度或者其他的一些参数。
用实常量来定义其参数只限于定义单层壳体。
如果同时用实常数和一个有效的壳体横截面来定义SHELL181单元的话，实常数定义将被忽略。
SHELL181单元同样适用用预整合的壳体横截面类型。
当该单元用GENS类型横截面来定义时，厚度或者材料的定义就都不需要了。
如果想了解更多的信息，请参见预整合通用壳体横截面的使用。
用实常数来定义厚度壳体单元的厚度可以在其各个节点定义，在该单元内，其厚度的变化被认为是光滑的。
总结一下ANSYS中不同单元之间的连接问题
论坛里常有人问不同单元之间的连接问题，我自己也一直被这个问题所困绕,最近从ANＳＹＳ工程分析进阶实例上知道了AＮSYS中不同单元之间的连接原则。
感觉收收获不小，现把它上传与大家共享。
一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可,不需要约束方程,否则需要耦合自由度与约事方程。
例如：(1)杆与梁、壳、体单元有公共节点即可，不需要约束方程。
ﻫ（2)梁与壳有公共节点怒可,也不需要约束写约束方程;壳梁自由度数目相同,自由度也相同，尽管壳的ｒoｔz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
(3)梁与体则要在相同位置建立不同的节点,然后在节点处耦合自由度与施加约束方程。
(４)壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的,有一定的局限性,只适用于小位移，下面介绍一种支持大位移算法的方法，ＭＰC法。
MPC即Mｕｌtiｐoint Conｓtraｉｎt,多点约束方程,其原理与前面所说的方程的技术几乎一致,将不连续、自由度不协调的单元网格连接起来,不需要连接边界上的节点完全一一对应。
MPC能够连接的模型一般有以下几种。
ﻫsolid 模型-ｓoｌid模型ﻫsｈelｌ模型-ｓheｌｌ模型ｓolｉｄ模型-shell 模型ｓolｉd 模型-beaｍ模型shelｌ模型-bｅam模型ﻫ在ＡNSＹS中,实现上述MPC技术有三种途径。
ﻫ(1）通过MPC18４单元定义模型的刚性或者二力杆连接关系。
定义MＰC1８4单元模型与定义杆的操作完全一致，而MPC单(2）利用约束元的作用可以是刚性杆（三个自由度的连接关系)或者刚性梁(六个自由度的连接关系)。
ﻫ方程菜单路径Ｍain Meｎｕ>pｒeproceｓsor>Coupliｎｇ/Cｅqn>shelｌ／ｓｏｌid Ｉｎｔerface创建壳与实体模型之间的装配关系。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
杆梁壳单元内力的输出
计算完成后，就可以对杆梁壳单元进行后处理，显示位移、应力（需打开单元开关显示开关）的操作和普通的单元的后处理是一样的。
但其内力的输出和显示需要进行特殊的处理。
首先要保证在定义单元类型时，打开了内力输出的选项；其次，在完成计算后需要通过定义单元表的方法来提取。
一、定义输出选项
Beam4单元说明的input summary中可以看到关于内力输出的开关选项说明：
二、定义单元表
查看对应项和序号：在单元说的后的单元输出表中查到MFOR(X,Y,Z)和MMOM(X,Y,Z)对应的项和序列号如下表，注意I和J节点上对应的序列号不同。
GUI方式如下图：
命令流过程如下：
（I节点处结果）
Etable,forx_i,smisc,1 ！x方向剪力（即轴力）
Etable,fory_i,smisc,2 ！y方向剪力
Etable,forz_i,smisc,3 ！z方向剪力
Etable,momx_i,smisc,4 ！x方向弯矩
Etable,momy_i,smisc,5 ！y方向弯矩
Etable,momz-i,smisc,6 ！z方向弯矩
（J节点处结果）
Etable,forx_j,smisc,7 ！x方向剪力（即轴力）
Etable,fory_j,smisc,8 ！y方向剪力
Etable,forz_j,smisc,9 ！z方向剪力
Etable,momx_j,smisc,10 ！x方向弯矩
Etable,momy_j,smisc,11 ！y方向弯矩
Etable,momz-j,smisc,12 ！z方向弯矩
三、绘制内力图
对线单元和面单元，其绘制单元表的结果的命令是不一样的。
线单元的内力云图的绘制为：PLLS. 面单元的内力云图的绘制为：PLETAB
四、注意事项
1.弯矩图的绘制与显示比例
在绘制弯矩图时（其他内力图也一样），可以给定一个显示比例系数Fact（默认等于1），当为负值时，反向显示，数值为放大系数。
同时还须注意的是，弯矩图的显示还与单元坐标系的方向有关，因此，如果同一段梁中各个单元的单元坐标系不一致（可以通过改变线的方向来改变线单元的单元坐标系），还将导致弯矩图显示不连续。
壳单元的后处理中更加要注意单元坐标系的方向，而且如果面上的壳单元划分不整齐，还将导致提取出来的内力图方向不统一，没有参考价值。
2.划分数的影响
要特别注意梁单元划分数对弯矩图的影响，除了因为有限元法是一种把连续结构体离散化的近似求解方法，还因为弯矩图是以梁单元的I和J节点结果进行连线绘制的。
可通过例子中用10份，2份，100份进行对比。