Ansys：施加隨坐標(biāo)變化的壓力載荷

2017-03-02 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

在實(shí)際工作中,經(jīng)常需要在結(jié)構(gòu)表面上施加壓力,并且壓力不是均勻分別得,而是隨結(jié)構(gòu)表面位置變化而分布。如果在實(shí)體結(jié)構(gòu)表面施加不均勻分布?jí)毫?需要在實(shí)體表面先覆蓋一層表面效應(yīng)單元SURF153或者SURF154,然后在表面效應(yīng)單元上施加這種壓力,此時(shí)與在殼體上施加載荷一樣進(jìn)行處理。

eg:在一個(gè)長(zhǎng)度和寬度為10m,厚度為0.1m殼體上施加按下來(lái)函數(shù)p(x,y)=c0+c1*x**3+c2*y**3

!創(chuàng)建隨函數(shù)分布得壓力通用宏:press_funct.mac
*create,press_funct,mac
cm,e_sel,elem
!******提取所有單元編號(hào)******
*get,ecount,elem,,count
*dim,enum,array,ecount
*get,enum(1),elem,,num,min
*do,i,2,ecount
enum(i)=elnext(enum(i-1))
*enddo
!******按單元定義每個(gè)單元壓力******
!pn(x,y)=c0+c1*x**3+c2*y**3
c0=0.1
c1=5
c2=5
*do,i,1,ecount
esel,s,,,enum(i)
nsel,r
*get,n_e,node,,count
xn=nx(nelem(enum(i),1))
yn=ny(nelem(enum(i),1))
zn=nz(nelem(enum(i),1))
pn1=c0+c1*xn**3+c2*yn**3
xn=nx(nelem(enum(i),2))
yn=ny(nelem(enum(i),2))
zn=nz(nelem(enum(i),2))
pn2=c0+c1*xn**3+c2*yn**3
xn=nx(nelem(enum(i),3))
yn=ny(nelem(enum(i),3))
zn=nz(nelem(enum(i),3))
pn3=c0+c1*xn**3+c2*yn**3
pn4=0
*if,n_e,eq,4,then
xn=nx(nelem(enum(i),4))
yn=ny(nelem(enum(i),4))
zn=nz(nelem(enum(i),4))
pn4=c0+c1*xn**3+c2*yn**3
*endif
sfe,enum(i),1,pres,,pn1,pn2,pn3,pn4
*enddo
cmsel,s,e_sel
!******刪除宏內(nèi)部應(yīng)用的變量******
cmdele,e_sel
xn=
yn=
zn=
pn1=
pn2=
pn3=
n_e=
ecount=
enum(1)=
anum=
*end
!******
!主程序
finish
/filname,byb
/prep7
et,1,shell63
r,1,0.1
mp,ex,1,2e11
mp,nuxy,1,0.3
mp,dens,1,7800
rectng,-5,5,-5,5
esize,1
mshape,0,2d
mshkey,1
amesh,all
/psf,pres,norm,2,0
press_funct
eplot
allsel,all
finish
save