[转帖]ansys载荷为位置函数的一个例子

VibInfo

　　设压力分布函数为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,... (最大单元号)
*get,emin,... (最小单元号)

*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

因为手头没有ansys软件和电脑，只能在网吧上网。上面命令语法可能有错。但思路不会有
问题，我以前做过。

因为要加载荷的面可能划分了很多个单元，且其单元号基本上不会是从1到N.所以“雨人”
的程序适用性不是很强。改进的程序只要选择要加载荷的面即可，应该很适用。

usefully

来的好,学习一下,谢谢楼主！

mzj5000

谢谢！学习了！