本文介紹了ANSYS中任意面任意方向任意變化的壓力施加方法。

　　在某些特殊的應(yīng)用場合，可能需要在結(jié)構(gòu)件的某個面上施加某個坐標(biāo)方向的隨坐標(biāo)位置變化的壓力載荷，當(dāng)然，這在一定程度上可以通過ANSYS表面效應(yīng)單元實現(xiàn)。如果利用ANSYS的參數(shù)化設(shè)計語言，也可以非常完美地實現(xiàn)此功能，下面通過一個小例子描述此方法。

　　!!!在執(zhí)行如下加載命令之前,請務(wù)必用選擇命令asel將需要加載的幾何面選擇出來

　　!!!

　　finish

　　/prep7

　　et,500,shell63

　　press=100e6

　　amesh,all

　　esla,s

　　nsla,s,1

　　! 如果載荷的反向是一個特殊坐標(biāo)系的方向,可在此建立局部坐標(biāo)系,并將

　　! 所有節(jié)點坐標(biāo)系旋轉(zhuǎn)到局部坐標(biāo)系下.

　　*get,enmax,elem,,num,max

　　dofsel,s,fx,fy,fz

　　fcum,add !!!將力的施加方式設(shè)置為"累加",而不是缺省的"替代"

　　*do,i,1,enmax

　　*if,esel,eq,1,then

　　*get,ae,elem,i,area !此命令用單元真實面積，如用投影面積，請用下幾條命令

　　! *get,ae,elem,i,aproj,x !此命令用單元X投影面積，如用真實面積，請用上一條命令

　　! *get,ae,elem,i,aproj,y !此命令用單元Y投影面積

　　! *get,ae,elem,i,aproj,z !此命令用單元Z投影面積

　　xe=centrx !單元中心X坐標(biāo)(用于求解壓力值)

　　ye=centry !單元中心Y坐標(biāo)(用于求解壓力值)

　　ze=centrz !單元中心Z坐標(biāo)(用于求解壓力值)

　　! 下面輸入壓力隨坐標(biāo)變化的公式,本例的壓力隨X和Y坐標(biāo)線性變化.

　　p_e=(xe-10)*press (ye-5)*press

　　f_tot=p_e*ae

　　esel,s,elem,,i

　　nsle,s,corner

　　*get,nn,node,,count

　　f_n=f_tot/nn

　　*do,j,1,nn

　　f,nelem(i,j),fx,f_n !壓力的作用方向為X方向

　　! f,nelem(i,j),fy,f_n !壓力的作用方向為Y方向

　　! f,nelem(i,j),fz,f_n !壓力的作用方向為Z方向

　　*enddo

　　*endif

　　esla,s

　　*enddo

　　aclear,all

　　fcum,repl !!!將力的施加方式還原為缺省的"替代"

　　dofsel,all

　　allsel