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

2、接觸面建立方法，通過238版本后的impgrid, expgrid等命令進行模型導人導出，花費很少的時間就可以解決復雜模型的接觸面建立問題。下面用一個實例來做一下說明。工況：3*3*3的模型，兩個 group ,模型正中間的一個單元為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)將除了需要加接觸面的單元以外的所有單元刪除，僅保存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)這時候需要將上面建好的group 2和接觸面進行保存save 2.sav(5)重新restore先前的模型，并將需要建立接觸面的實體刪除，再將整個模型用 命令進行導出：rest 1.savexpgrid 1.fac3d(6)寫到這兒，大家應該明白了吧。下面的工作是這樣的：將已建好的 group的文件rest,并把其他的模型用impgrid命令進彳T導入：impgrid 1.flac3dexpgrid2和接觸面得到的就是已經建立好接觸面的模型卜面進行簡單的計算，以驗證接觸面是否工作。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可以看到接觸面的接應力分布結果：得到的結果還是蠻理想的Interface Shear Stress(0.0040e+0001a Di。031 號0DOQ.QN0«*QQ01qi 1 OQCHKs+OW l.D040e*C03-lD 2m口加50g J 2.0M0e+0031O 3.0000*003 二 m.口的口無十口。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與周圍組分開，這個組的周圍點雖說和周圍相鄰的周圍組的點有相同的位置，但有不同的節點然后在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、接觸面參數的選擇kn ks應該根據13t驗確定，一般可取接觸面兩側彈性模量較小值的10倍（手冊上所說）這個參數取值確實沒有太多經驗，最好的方法還是用手冊建議的 10 倍剛度的方法。如果這種方法的計算結果合理，那么最好不過了，因為手冊中建議的方法自有它的道理，但具體是何道理，我也不了解。當然，如果有試驗，那么最好通過試驗反分析的方法來確定參數。Ks 切向剛度通常采用接觸面的剪切試驗得到* q1 l5l# b" K& B e/ C.jKn 法向剛度通常設置一個較大的值，但是這個值不能太大，否則又會存在收斂慢的問題，我所見到的 Ks 在 1e7-1e9 之間5、就是記錄界面上的應力隨時間的變化該用哪個命令我看例子上用 hist interface i ss