如何在Ansys/lsdyna中給Shell 163 賦值變厚度(注：不是Ansys下的Shell63號單元)
瀏覽：80815
這個問題困擾了小弟好久，盼有人能予以解答~
小弟一直想做一個變厚度鋼管，在網上看了好多Ansys變厚度板的代碼，也就是這個：
“FINISH
/clear
/prep7
et,1,63
rect,,10,,10
eshape,2
esize,,20
amesh,1
eplot
mxnode=ndinqr(0,14)
*dim,thick,,mxnode
*do,node,1,mxnode
*if,ndinqr(node,1),eq,1,then
thick(node)=0.5+0.2*nx(node)+0.02*ny(node)**2
*else
thick(node)=0
*endif
*enddo
node=
mxnode=
rthick,thick(1),1,2,3,4“
但是我在用Ansys/lsdyna做的時候遇到了極大的問題！！！問題如下：
1. 代碼中的單元是Shell 63號單元，此單元在Ansys/lsdyna下是沒有的，Ansys/lsdyna是Explicit的，所以只能用Shell 163號單元來做。不過使用相同的代碼，我在用Ansys/Lsdyna下已經用163號單元劃好網格后，將代碼相應代入命令流中(從mxnode=ndinqr(0,14)開始代入)，完全實現不了Shell 163號單元的命令使用。這樣賦值之后，Ansys里面模型的每個單元仍舊是等厚度的。我只是想做一個長600的空心方管，同一個x軸坐標的節點是一個厚度，譬如x=1的節點都為0.5，x=2的節點都為0.7厚度，但是這樣用上面的代碼實現不了。這個讓小弟非常不知道如何解決。
2. 另外，這樣賦值之后，不僅每個單元是等厚度的，而且每個單元的獨立成一個Part，這樣在我做比較細的模型的時候（我的管一共13000多個單元），每一個單元獨立成為一個Part，這樣計算的時候，也不是很穩定的。因為我的模型等X值的單元其實屬性是完全相同的，所以能否設定每一個X值的單元是一個Part呢？
補充問題：有一次給小弟逼急了，我拿Hypermesh畫個管，一層一層的賦值t1 t2 t3 t4，但是這樣之后又有一個問題，就是同一個X值的單元雖然屬于一個Component，但是由于是按照Component賦 t1 t2 t3 t4的，并不能保證每個單元上的節點方向都是一樣的，也就是說，有的 t1=t2>t3=t4，但是有的t2=t3>t1=t4，這樣一層一層賦，也是連續變厚度板。
小弟只想知道能不能在Ansys/LSDYNA中合理的實現命令，非常想知道，關于補充問題能解決就好，不能解決也沒關系~希望有大神能 予以幫助，小弟已經快被搞死了
小弟一直想做一個變厚度鋼管，在網上看了好多Ansys變厚度板的代碼，也就是這個：
“FINISH
/clear
/prep7
et,1,63
rect,,10,,10
eshape,2
esize,,20
amesh,1
eplot
mxnode=ndinqr(0,14)
*dim,thick,,mxnode
*do,node,1,mxnode
*if,ndinqr(node,1),eq,1,then
thick(node)=0.5+0.2*nx(node)+0.02*ny(node)**2
*else
thick(node)=0
*endif
*enddo
node=
mxnode=
rthick,thick(1),1,2,3,4“
但是我在用Ansys/lsdyna做的時候遇到了極大的問題！！！問題如下：
1. 代碼中的單元是Shell 63號單元，此單元在Ansys/lsdyna下是沒有的，Ansys/lsdyna是Explicit的，所以只能用Shell 163號單元來做。不過使用相同的代碼，我在用Ansys/Lsdyna下已經用163號單元劃好網格后，將代碼相應代入命令流中(從mxnode=ndinqr(0,14)開始代入)，完全實現不了Shell 163號單元的命令使用。這樣賦值之后，Ansys里面模型的每個單元仍舊是等厚度的。我只是想做一個長600的空心方管，同一個x軸坐標的節點是一個厚度，譬如x=1的節點都為0.5，x=2的節點都為0.7厚度，但是這樣用上面的代碼實現不了。這個讓小弟非常不知道如何解決。
2. 另外，這樣賦值之后，不僅每個單元是等厚度的，而且每個單元的獨立成一個Part，這樣在我做比較細的模型的時候（我的管一共13000多個單元），每一個單元獨立成為一個Part，這樣計算的時候，也不是很穩定的。因為我的模型等X值的單元其實屬性是完全相同的，所以能否設定每一個X值的單元是一個Part呢？
補充問題：有一次給小弟逼急了，我拿Hypermesh畫個管，一層一層的賦值t1 t2 t3 t4，但是這樣之后又有一個問題，就是同一個X值的單元雖然屬于一個Component，但是由于是按照Component賦 t1 t2 t3 t4的，并不能保證每個單元上的節點方向都是一樣的，也就是說，有的 t1=t2>t3=t4，但是有的t2=t3>t1=t4，這樣一層一層賦，也是連續變厚度板。
小弟只想知道能不能在Ansys/LSDYNA中合理的實現命令，非常想知道，關于補充問題能解決就好，不能解決也沒關系~希望有大神能 予以幫助，小弟已經快被搞死了
