采用ansys软件分析时，有时候想知道提取截面的内力—弯矩，利用ansys单元库的梁单元分析自然可以直接提取，但是采用三维实体单元分析时就不一定能提取了，此时可以采用ansys后处理中的面操作功能。下面看一个实例吧。

!悬臂梁的内力计算
finish
/clear
/prep7
et,1,solid95
mp,ex,1,2e11
mp,prxy,1,0.3
blc4,2,3,0.2,0.3,4
da,2,all
fk,1,fy,-2e4
fk,3,fy,-2e4
fk,3,fx,0.8e4
fk,4,fx,0.8e4
sfa,1,1,pres,1e6
esize,0.05
vmesh,all
finish
/solu
solve
finish
/post1
wpoff,,,2
sucr,suz2,cplane,3      !创建面suz2
sumap,mysx,s,x           !映射X方向应力
sumap,mysy,s,y           !映射Y方向应力
sumap,mysz,s,z           !映射Z方向应力
sumap,mysxy,s,xy        !映射xy方向应力
sumap,mysyz,s,yz        !映射yz方向应力
sumap,mysxz,s,xz        !映射xz方向应力
supl,suz2                       !显示面本身
supl,suz2,mysz             !显示面项mysz
supl,suz2,mysyz            !显示面项mysyz
supl,suz2,mys               !矢量显示应力

supr,all,mysz                 !列表显示面项mysz

sueval,xforce,mysxz,intg        !求截面上FX，理论结果为-16000,误差1%
sueval,yforce,mysyz,intg        !求截面上FY理论结果为40000，误差0.5%
sueval,zforce,mysz,intg          !求截面上FZ理论结果为-6000，没有误差

sueval,mya,da,sum                 !求截面面积并赋給变量mya
sueval,myya,gcy,intg               !求关于x轴的面积矩并赋给变量myya
myya=myya/mya                      !得到面积重心到x轴的距离=面积矩/面积
sueval,myxa,gcx,intg               !求关于y轴的面积矩并赋给变量myya
myxa=myxa/mya                      !得到面积重心到y轴的距离=面积矩/面积

sucalc,szgcy,mysz,mult,gcy     !计算mysz乘以gcy，并赋给面项szgcy
sueval,mx1,szgcy,intg              !对面项szgcy在面上积分得到mx1        
sucalc,szgcx,mysz,mult,gcx     !计算mysz乘以gcy，并赋给面项szgcx
sueval,my1,szgcx,intg              !对面项szgcx在面上积分得到my1

!!!!!!!!!上述弯矩基于总体直角坐标系原点而言的，应对面积重心取矩，将内力简化到面积重心上。

mx1=mx1-zforce*myya   !理论结果80000,0.08%
my1=my1-zforce*myxa   !理论结果-32000,0.2%