- 积分
- 0
- 注册时间
- 2014-3-13
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2014-5-2 10:35:45
|
显示全部楼层
来自 北京
给你个提取solid45单元的命令流例子吧。【注:这是以往论坛上的例子,并非本人创作,但很有用的哦】
/PREP7
ET,1,SOLID45
MP,EX,1,2E11
MP,PRXY,1,0.3 !定义单元类型、材料特性
BLC4,2,3,0.2,0.3,4 !通过两个角点或Z方向的深度产生一个矩形面或块体
DA,1,ALL
!施加荷载
FK,8,FY,2E4 !在关键点施加集中荷载
FK,7,FY,2E4
FK,7,FX,0.8E4
FK,6,FX,0.8E4
SFA,4,1,PRES,1E6 !在选择的面上加表面荷载PRES压力
SFA,2,1,PRES,1E6
ESIZE,0.05
VMESH,ALL
FINISH
/SOLU
SOLVE !生成有限元模型并求解
FINISH
/POST1 !进入后处理层
WPOFF,,,2 !偏移工作平面
SUCR,SUZ2,CPLANE,3 !移动工作平面、创建面SUZ2
!映射X、Y和Z方向应力
SUMAP,MYSX,S,X !映射结果数据到所选择的面上
SUMAP,MYSY,S,Y
SUMAP,MYSZ,S,Z
!映射XY、YZ和XZ方向应力
SUMAP,MYSXY,S,XY
SUMAP,MYSYZ,S,YZ
SUMAP,MYSXZ,S,XZ
!显示面本身和面上的应力
SUPL,SUZ2
SUPL,SUZ2,MYSZ
SUPL,SUZ2,MYSYZ
SUPL,SUZ2,MYS !矢量显示应力
!列表显示MYSZ面项
SUPR,ALL,MYSZ
!求截面上FX,理论结果为16000,计算结果为15999.9983 误差1%
SUEVAL,XFORCE,MYSXZ,INTG !Integral of MYSXZ over surface积分
!求截面上FY,理论结果为-360000,计算结果为 -360000.069
SUEVAL,YFORCE,MYSYZ,INTG
!求截面上FZ,理论结果为-60000,计算结果为-60001.2485
SUEVAL,ZFORCE,MYSZ,INTG
!求截面面积并赋给MYA变量
SUEVAL,MYA,DA,SUM
!求关于X轴的面积矩并赋给变量MYYA
SUEVAL,MYYA,GCY,INTG
!得到面积重心到X轴的距离=面积矩/面积
MYYA=MYYA/MYA
!求关于Y轴的面积矩并赋给变量MYXA
SUEVAL,MYXA,GCX,INTG
!得到面积重心到Y轴的距离=面积矩/面积
MYXA=MYXA/MYA
!计算MYSZ*GCY,并赋给面项SZGCY
SUCALC,SZGCY,MYSZ,MULT,GCY
!对面项SZGCY在面上积分得到MX1
SUEVAL,MX1,SZGCY,INTG
!计算MYSZ*GCZ,并赋给面项SZGC
SUCALC,SZGCX,MYSZ,MULT,GCX
!对面项SZGCX在面上积分得到MY1
SUEVAL,MY1,SZGCX,INTG
!上述弯矩基于总体直角坐标系原点而言的,应对面积重心取矩,将内力简化到面积重心上
!理论结果为320000,计算结果为320437.403
MX1=MX1-ZFORCE*MYYA
!理论结果-32000, 计算结果为-32062.5187
MY1=MY1-ZFORCE*MYXA |
|