FLAC3D接觸面建立方法

1、帖子：復(fù)雜模型內(nèi)部接觸面的建立方法 首先為什么要建立接觸面: 如果考慮接觸界面對(duì)問(wèn)題的影響或接觸界面的變化特征時(shí)就的設(shè)置接觸單元。如樁與土之間時(shí)，由于樁土，地連墻與土，由于剛度相差太大，位移不連續(xù)，有限元中尤其 如此；設(shè)置接觸可以反映接觸面之間得剪切滑移或侵入行為，如節(jié)理面，斷層的剛度相差較大， 所以它們之間的接觸單元定義后可以減少由剛度差別較大產(chǎn)生的誤差影響(如位移) 當(dāng)然在定義接觸單元后降低了剛度的影響也會(huì)提高軟件的運(yùn)行速度。供參考。接觸面的建立方法：1、就是大家通常用的“移來(lái)移去”法，手冊(cè)上也有例子，一般模型推薦用此方法。2、lakewater所推薦的建立復(fù)雜模型的方法推薦一種更加實(shí)用的

2、接觸面建立方法，通過(guò)238版本后的impgrid, expgrid等命令進(jìn)行模型導(dǎo)人導(dǎo)出，花費(fèi)很少的時(shí)間就可以解決復(fù)雜模型的接觸面建立問(wèn)題。下面用一個(gè)實(shí)例來(lái)做一下說(shuō)明。工況：3*3*3的模型，兩個(gè) group ,模型正中間的一個(gè)單元為gropu 2 ,其他為group 1gen zone brick size 3 3 3group 2 range x 1 2 y 1 2 z 1 2group 1 range gr 2 not save 1.sav(1)將模型文件存為1.sav;(2)將除了需要加接觸面的單元以外的所有單元?jiǎng)h除，僅保存group 2命令為：del ran grop2 not;(3

3、)在group 2周圍建立接觸面interface 1 facerange x 1y1 2 z1 2interface 1 facerange x 2y1 2 z1 2interface 1 facerange x 12y 1 z1 2interface 1 facerange x 12y 1 2z 1interface 1 facerange x 12y 1 2z 2或 interface 1 face 也可以建成后可以看到接觸面的效果：plo inter yel ske fix x y z ran z 0 ini den 20005 set gra 0 0 -10Interface Loc

4、ationsSketchMl 白 gfac O.OOOe+ ODOLinestyle (4)這時(shí)候需要將上面建好的group 2和接觸面進(jìn)行保存save 2.sav(5)重新restore先前的模型，并將需要建立接觸面的實(shí)體刪除，再將整個(gè)模型用 命令進(jìn)行導(dǎo)出：rest 1.savexpgrid 1.fac3d(6)寫到這兒，大家應(yīng)該明白了吧。下面的工作是這樣的：將已建好的 group的文件rest,并把其他的模型用impgrid命令進(jìn)彳T導(dǎo)入：impgrid 1.flac3dexpgrid2和接觸面得到的就是已經(jīng)建立好接觸面的模型卜面進(jìn)行簡(jiǎn)單的計(jì)算，以驗(yàn)證接觸面是否工作。model elain

5、terface 1 prop kn 20e6 ks 20e6 c 10e3 f 15app nstr -200e3 ran x 0 1 y 1 2 z 3solve可以看到接觸面的接應(yīng)力分布結(jié)果：得到的結(jié)果還是蠻理想的Interface Shear Stress(0.0040e+0001a Di。031 號(hào)0DOQ.QN0«*QQ01qi 1 OQCHKs+OW l.D040e*C03-lD 2m口加50g J 2.0M0e+0031O 3.0000*003 二 m.口的口無(wú)十口。310 4 dabDe*o3 二 4.00Q0c*003iD 5.0000e*OQ3 二 5.00fl0

6、e+fl031o 6 0000*00-3(G.CllMCletdOaiD 7 0口口the*0ba J 7.0IM0e*(l031& B.OOOO-b*ODS S.D(JflCI®+0031O SI 0000+003 I芻0附加鋪口q打骷lDM inwrvalT 1 dei/iEMB3、用 gen separate group1 int 1 wrap group1 group2此方法是在已建好的模型上使用，先把 group1與周圍組分開，這個(gè)組的周圍點(diǎn)雖說(shuō)和周圍相鄰的周圍組的點(diǎn)有相同的位置，但有不同的節(jié)點(diǎn)然后在group1上建立接觸面，此法我稱之為“分離法”ngen zone

7、 brick size 3 3 3group 2 range x 1 2 y 1 2 z 1 2group 1 range gr 2 notsave 1.savdel ran group 2 notinterface 1 facerange x1y 12z 12interface 1 facerange x2y 12z 12interface 1 facerange x12 y1z 12interface 1 facerange x12 y12 z1interface 1 facerange x12 y12 z2;interface 1 face 也可以plo inter yel ske sa

8、ve 2.savrest 1.savexpgrid 1.fac3drest 2.savimpgrid 1.flac3dplo inter yel skepaumodel elafix x y z ran z 0ini den 20005set gra 0 0 -10interface 1 prop kn 20e6 ks 20e6 c 10e3 f 15app nstr -200e3 ran x 0 1 y 1 2 z 3solveplo int sst4、接觸面參數(shù)的選擇kn ks應(yīng)該根據(jù)13t驗(yàn)確定，一般可取接觸面兩側(cè)彈性模量較小值的10倍（手冊(cè)上所說(shuō)）這個(gè)參數(shù)取值確實(shí)沒有太多經(jīng)驗(yàn)，最好的方法還是用手冊(cè)建議的 10 倍剛度的方法。如果這種方法的計(jì)算結(jié)果合理，那么最好不過(guò)了，因?yàn)槭謨?cè)中建議的方法自有它的道理，但具體是何道理，我也不了解。當(dāng)然，如果有試驗(yàn)，那么最好通過(guò)試驗(yàn)反分析的方法來(lái)確定參數(shù)。Ks 切向剛度通常采用接觸面的剪切試驗(yàn)得到* q1 l5l# b" K& B e/ C.jKn 法向剛度通常設(shè)置一個(gè)較大的值，但是這個(gè)值不能太大，否則又會(huì)存在收斂慢的問(wèn)題，我所見到的 Ks 在 1e7-1e9 之間5、就是記錄界面上的應(yīng)力隨時(shí)間的變化該用哪個(gè)命令我看例子上用 hist interface i ss