在某些特殊的應(yīng)用場(chǎng)合,可能需要在結(jié)構(gòu)件的某個(gè)面上施加某個(gè)坐標(biāo)方向的隨坐標(biāo)位置變化的壓力載荷,當(dāng)然,這在一定程度上可以通過ANSYS表面效應(yīng)單元實(shí)現(xiàn)。如果利用ANSYS的參數(shù)化設(shè)計(jì)語言,也可以非常完美地實(shí)現(xiàn)此功能,下面通過一個(gè)小例子描述此方法。
!!!在執(zhí)行如下加載命令之前,請(qǐng)務(wù)必用選擇命令asel將需要加載的幾何面選擇出來
!!!
finish
/prep7
et,500,shell63
press=100e6
amesh,all
esla,s
nsla,s,1
! 如果載荷的反向是一個(gè)特殊坐標(biāo)系的方向,可在此建立局部坐標(biāo)系,并將
! 所有節(jié)點(diǎn)坐標(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 !此命令用單元真實(shí)面積,如用投影面積,請(qǐng)用下幾條命令
! *get,ae,elem,i,aproj,x !此命令用單元X投影面積,如用真實(shí)面積,請(qǐng)用上一條命令
! *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 !壓力的作用方向?yàn)閄方向
! f,nelem(i,j),fy,f_n !壓力的作用方向?yàn)閅方向
! f,nelem(i,j),fz,f_n !壓力的作用方向?yàn)閆方向
*enddo
*endif
esla,s
*enddo
aclear,all
fcum,repl !!!將力的施加方式還原為缺省的"替代"
dofsel,all
allsel