在任意面施加任意方向任意變化的壓力
如果利用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