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

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

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

eg:一個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
!變厚度殼的實現(xiàn)
*SET,mxnode,ndinqr(0,14)
!提取節(jié)點總數(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ù)計算每個節(jié)點的厚度,并將每個節(jié)點厚度記錄在厚度數(shù)組thick(mxnode)
*SET,node,
*SET,mxnode,
!刪除變量node,mxnode
rthick,thick(1),1,2,3,4
!分配單元厚度
/eshape,1
eplot
save