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

雨人载荷为位置函数的一个例子????设压力分布函数为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坐标） ??f=2*xe*sin(ye)*cos(ze) ??sfe,i,1,pres,,f,,, ??*enddo ?上贴有误，面载荷没有除以面积，修改如下： ??????*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 ? fat_head_fish 好方法。 ??增加一点建议。 ??先选要加载荷的面，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.所以“雨人”的程序适用性不是很强。改进的程序只要选择要加载荷的面即可，应该很适用。 ??希望能与斑竹多交流（多给我加点积分，呵呵）。还有，“雨人”的*GET,AA,ELEM,I,AREA 命令前也应有我的这个 asel,s,... 。不然f/aa就问题大了：）? 雨人 *GET,AA,ELEM,I,AREA 是获取单元的面积，与面无关的。chenguanci 斑竹大哥和fat_head_fish 大哥你们好！ 非常感谢你们的例子和建议。我想请教个问题： 我想在实体面上加载连续移动的面载，这个面载是在一个圆形的面上， 用NSEL命令可以选定这个圆形区域内的节点，但我害怕这样加载不准确，但如何才能选中一个圆形区域，用ASEL命令好像选不中，因为事先并没有建立这个圆形的面，只是选中一个已存在面的一个小区域。 希望不胜赐教！无限感激！