ansys提出实体单元建立弯矩

对于3D实体单元的内力计算，可采用面操作（应该是V8以上吧，没有细考！）。例如用悬臂梁采用3D实体单元，其某个截面的的内力计算如下，并具有详细解释。这种方法较原来的更方便，大家不放一用。

!******************************************************

!EX4.31 面操作及悬臂梁的内力计算

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

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方向应力

SUPL,SUZ2

SUPL,SUZ2,MYSZ !显示面本身、面项MYSZ

SUPL,SUZ2,MYSYZ

SUPL,SUZ2,MYS !显示面项MYSYZ，矢量显示应力

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轴的面积矩并赋给变量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 !对面项SZGCX在面上积分得到MY1

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

MX1=MX1-ZFORCE*MYYA !理论结果为80000，误差为0.08%

MY1=MY1-ZFORCE*MYXA !理论结果-32000，误差为0.2%