设压力分布函数为F=2RSIN(AF)COS(XT)
命令流
*do,i,1,N (共N个单元)
csys,2 (激活极坐标系)
esel,,,i (选择第i个单元)
xe=centrx(i) (xe为单元中心X坐标)
ye=centry(i) (ye为单元中心y坐标)
ze=centrz(i) (ze为单元中心z坐标)
*GET,AA,ELEM,I,AREA
f=2*sin(xe)*cos(ye)*cos(ze)
sfe,i,1,pres,,f/AA,,,
*enddo
=====================
先选要加载荷的面,get其最大单元号和最小单元号。再做上面提到的循环,循环从最小单
元号到最大单元号。循环过程中,判断所循环的单元号在所选的面内是否存在。
命令流:
asel,... (要加载荷的面)
esla,,,, (和面相联的单元)
*GET,eMAX,ELEM,0,NUM,MAX !获取最大单元号
*GET,eMIN,ELEM,0,NUM,MIN !获取最小单元号
*do,i,emin,emax
*if,esel(i),eq,1,then (该句语法请查帮助,这里凭记忆写出)
!!!*do,i,1,N (共N个单元)
csys,2 (激活极坐标系)
esel,,,i (选择第i个单元)
xe=centrx(i) (xe为单元中心X坐标)
ye=centry(i) (ye为单元中心y坐标)
ze=centrz(i) (ze为单元中心z坐标)
*GET,AA,ELEM,I,AREA
f=2*sin(xe)*cos(ye)*cos(ze)
sfe,i,1,pres,,f/AA,,,
*endif (判断结束)
*enddo