設壓力分布函數為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