Ansys apdl：參數(shù)化建模：連續(xù)變厚板殼模型

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

實(shí)際中,某些板殼結(jié)構(gòu)中存在厚度不均勻的殼體,必須利用實(shí)常數(shù)定義殼的厚度,如果厚度按照某種規(guī)律進(jìn)行變化,這就意味著每個(gè)殼單元必須采用一個(gè)獨(dú)立的實(shí)常數(shù),并且每個(gè)殼單元的各節(jié)點(diǎn)具有不同的厚度。用人工方法定義這些實(shí)常數(shù)非常困難,也難將它們一一分配倒適當(dāng)?shù)奈恢?所以必須利用APDL和殼單元厚度函數(shù)功能批處理這些工作。

eg:一個(gè)10×10的矩形面劃分20×20的shell單元網(wǎng)格模型,其厚度為thick(x,y)=0.5+0.2*x+0.02*y**2

!初始化ansys環(huán)境
finsih
/clear
/filname,shell-thick
!創(chuàng)建殼體網(wǎng)格模型
/prep7
et,1,63
rect,,10,,10
eshape,2
esize,,20 +
amesh,1
eplot
!變厚度殼的實(shí)現(xiàn)
*SET,mxnode,ndinqr(0,14)
!提取節(jié)點(diǎn)總數(shù)mxnode
*dim,thick,,mxnode
!定義厚度數(shù)組thick(mxnode)
*do,node,1,mxnode
*if,ndinqr(node,1),eq,1,then
*SET,thick(node),0.5+0.2*nx(node)+0.02*ny(node)**2
*else
*SET,thick(node),0
*endif
*enddo
!根據(jù)厚度分布函數(shù)計(jì)算每個(gè)節(jié)點(diǎn)的厚度,并將每個(gè)節(jié)點(diǎn)厚度記錄在厚度數(shù)組thick(mxnode)
*SET,node,
*SET,mxnode,
!刪除變量node,mxnode
rthick,thick(1),1,2,3,4
!分配單元厚度
/eshape,1
eplot
save