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

　　在某些特殊的應用場合，可能需要在結構件的某個面上施加某個坐標方向的隨坐標位置變化的壓力載荷，當然，這在一定程度上可以通過ANSYS表面效應單元實現。如果利用ANSYS的參數化設計語言，也可以非常完美地實現此功能，下面通過一個小例子描述此方法。

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

　　!!!

　　finish

　　/prep7

　　et,500,shell63

　　press=100e6

　　amesh,all

　　esla,s

　　nsla,s,1

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

　　! 所有節點坐標系旋轉到局部坐標系下.

　　*get,enmax,elem,,num,max

　　dofsel,s,fx,fy,fz

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

　　*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坐標(用于求解壓力值)

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

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

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

　　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