workbench實(shí)例計(jì)算solid186單元與solid185單元性能對比的案例
solid186與solid185單元結(jié)果對比下載
EDP材料參數(shù)計(jì)算
EC1=6*sin(alph1)/(3-sin(alph1))
EC2=6*Coh1*cos(alph1)/(3-sin(alph1))
tb,Edp,3,1,2,LYFUN
tbdata,1,EC1,EC2
tb,Edp,3,1,1,LFpot
tbdata,1,EC1
!內(nèi)部挖空圍巖
mp,ex,4,3.6e9
mp,prxy,4,0.32
mp,dens,4,2200
tb,Edp,4,1,2,LYFUN
tbdata,1,EC1,EC2
tb,Edp,4,1,1,LFpot
tbdata,1,EC1
!===============
結(jié)果對比：
一、等效應(yīng)力場
Plane42單元等效應(yīng)力場
Plane182單元等效應(yīng)力場
二、Y方向位移場
Plane42單元Y方向位移場
Plane182單元Y方向位移場
三、襯砌彎矩
Plane42單元襯砌彎矩
Plane182單元襯砌彎矩
通過對比可發(fā)現(xiàn)，兩者計(jì)算結(jié)果誤差較小，說明上述參數(shù)等效方法能較好的實(shí)現(xiàn)EDP材料模型參數(shù)的輸入，故而同學(xué)們在做類似巖土類模擬時(shí)可選擇如下方法進(jìn)行：
一、采用低級單元Plane42、Solid45，材料模型采用經(jīng)典DP模型；
二、采用高級單元Plane182、Solid185，材料模型采用EDP模型，模型參數(shù)可按本文所述方法進(jìn)行計(jì)算
下載地址：solid186與solid185單元結(jié)果對比
展開 ANSYS中的節(jié)點(diǎn)解與單元解是怎么回事？附solid186與solid185單元結(jié)果對比文檔下載
而20節(jié)點(diǎn)單元縮減積分后，有7個(gè)積分點(diǎn)，應(yīng)該輸出7個(gè)單元解，經(jīng)過計(jì)算如圖8所示：
圖8
圖8正好是7個(gè)輸出解。
Abaqus的計(jì)算表明單元輸出解果然是輸出單元積分點(diǎn)的值，采用完全積分和縮減積分單元輸出解不一樣，求解精度不一樣。
那么為什么ANSYS則沒有這種規(guī)律呢？
其實(shí)后臺程序計(jì)算是肯定是按照理論上走的，也就是先得到節(jié)點(diǎn)的位移，再得到單元積分點(diǎn)的應(yīng)力應(yīng)變，再外推得到各個(gè)單元節(jié)點(diǎn)的應(yīng)力應(yīng)變，最后平均得到節(jié)點(diǎn)解。
ANSYS之所以顯示的單元解不是單元積分點(diǎn)的解，而是各個(gè)節(jié)點(diǎn)的解，是因?yàn)锳NSYS已經(jīng)在得到單元積分點(diǎn)的解之后經(jīng)過外推得到了單元各個(gè)角節(jié)點(diǎn)的解，但是還沒有做平均。
也就是，ANSYS的單元解，其實(shí)不能完全看作單元解，筆者稱之為單元角節(jié)點(diǎn)解。
下載地址：solid186與solid185單元結(jié)果對比文檔下載
展開 