real,1
esys,0 (或 aatt) 激活单元类型 2,材料号 2,实常数 1,单元坐标系
vsweep,all,,, 扫掠网格
csys,4 激活坐标系 4
/prep7
numstr,kp,100 !define the following keypoint number start with with the 100
l,1,2,4 !如果 CSYS=0 则生成直线,如果 CSYS=1 则生成弧线,这个命令与当前的坐标系统有
关而楼上师兄的 LSTR 则始终生成直线.
lsel , !取线
wprof,,12 !移坐标
alsv !拾取一选定实体上的所有面
nsla !同理,拾取一选定面上的所有节点
aatt,1,1,1 !等效于楼上的 MAT,1 TYPE,1 REAL, 1 对面定义属性
mshke,0
!网格格划分进行限定:采用 FREE 进行划分;网格形状为 四 边形或六面体
mshape,1,2d
vmesh ,2 !划分实体网格,后面的参数是实体编号如:2
/solu !进入求解过程
antype,static !选择求解类型为静力分析
asel,s,loc,x,
nsla
d,all,uy,,,,,roty,rotz !对选定的面上的所有节点施加 UY ROTY ROTZ 的对称约束.
allsel !恢复全部选择等效于:ASELL,ALL ESEL,ALL NSEL,ALL
asel,s,,,1
sfa,all,1,press,1000 !对选定的面 1 施加均布力 1000
allsel
/stat,slou !显示求解状况
solve
/post1 !进入后处理
set,list !列出求解的步数及相关信息
set,last !读取最后一步结果
plns,s,eqv,,1 !绘出节点的等效应力云图