Ansys：施加隨坐標變化的壓力載荷

2017-03-02  by:CAE仿真在線  來源:互聯網

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

eg:在一個長度和寬度為10m,厚度為0.1m殼體上施加按下來函數p(x,y)=c0+c1*x**3+c2*y**3

!創建隨函數分布得壓力通用宏:press_funct.mac
*create,press_funct,mac
cm,e_sel,elem
!******提取所有單元編號******
*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
!******按單元定義每個單元壓力******
!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
!******刪除宏內部應用的變量******
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

相關標簽搜索：Ansys：施加隨坐標變化的壓力載荷 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓