载荷为位置函数的一个例子

来源:百度文库 编辑:神马文学网 时间:2024/04/30 11:31:13
设压力分布函数为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.所以“雨人”的程序适用性不是很强。改进的程序只要选择要加载荷的面即可,应该很适用。