接觸面建立方法

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

2、立方法，通過238版本后的impgrid，expgrid等命令進(jìn)行模型導(dǎo)入導(dǎo)出，花費(fèi)很少的時(shí)間就可以解決復(fù)雜模型的接觸面建立問題。下面用一個(gè)實(shí)例來做一下說明。工況：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 notsave 1.sav（1）將模型文件存為1.sav；（2）將除了需要加接觸面的單元以外的所有單元?jiǎng)h除，僅保存group 2 命令為：del ran grop 2 not；（3）在group

3、 2 周圍建立接觸面interface 1 face range x 1 y 1 2 z 1 2interface 1 face range x 2 y 1 2 z 1 2interface 1 face range x 1 2 y 1 z 1 2interface 1 face range x 1 2 y 1 2 z 1interface 1 face range x 1 2 y 1 2 z 2或interface 1 face也可以 建成后可以看到接觸面的效果： plo inter yel ske（4）這時(shí)候需要將上面建好的group 2 和接觸面進(jìn)行保存 save 2.sav（5）重新r

4、estore先前的模型，并將需要建立接觸面的實(shí)體刪除，再將整個(gè)模型用expgrid命令進(jìn)行導(dǎo)出： rest 1.savexpgrid 1.fac3d（6）寫到這兒，大家應(yīng)該明白了吧。下面的工作是這樣的：將已建好的group 2和接觸面的文件rest，并把其他的模型用impgrid命令進(jìn)行導(dǎo)入:impgrid 1.flac3d得到的就是已經(jīng)建立好接觸面的模型下面進(jìn)行簡(jiǎn)單的計(jì)算，以驗(yàn)證接觸面是否工作。model elafix x y z ran z 0ini den 20005set gra 0 0 -10interface 1 prop kn 20e6 ks 20e6 c 10e3 f 15ap

5、p nstr -200e3 ran x 0 1 y 1 2 z 3solve可以看到接觸面的接應(yīng)力分布結(jié)果：得到的結(jié)果還是蠻理想的3、用gen separate group1 int 1 wrap group1 group2此方法是在已建好的模型上使用，先把group1與周圍組分開，這個(gè)組的周圍點(diǎn)雖說和周圍相鄰的周圍組的點(diǎn)有相同的位置，但有不同的節(jié)點(diǎn)然后在group1上建立接觸面，此法我稱之為“分離法”ngen zone brick size 3 3 3group 2 range x 1 2 y 1 2 z 1 2group 1 range gr 2 notsave 1.savdel ran

6、group 2 notinterface 1 face range x 1 y 1 2 z 1 2interface 1 face range x 2 y 1 2 z 1 2interface 1 face range x 1 2 y 1 z 1 2interface 1 face range x 1 2 y 1 2 z 1interface 1 face range x 1 2 y 1 2 z 2;interface 1 face也可以plo inter yel skesave 2.savrest 1.savexpgrid 1.fac3drest 2.savimpgrid 1.flac3dp

7、lo 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ù)試驗(yàn)確定,一般可取接觸面兩側(cè)彈性模量較小值的10倍（手冊(cè)上所說）這個(gè)參數(shù)取值確實(shí)沒有太多經(jīng)驗(yàn)，最好的方法還是用手冊(cè)建議的10倍剛度的方法。如果這種方法的計(jì)算結(jié)果合理，那么最好不過了，因?yàn)槭謨?cè)中建議的方法自有它的道理，但具體是何道理，我也不了解。當(dāng)然，如果有試驗(yàn)，那么最好通過試驗(yàn)反分析的方法來確定參數(shù)。Ks切向剛度通常采用接觸面的剪切試驗(yàn)得到1 u( x. B0 C$ u2 K( * q1 l5 l# b" K& B  e/ C. jKn法向剛度通常設(shè)置一個(gè)較大的值，但是這個(gè)值不能太大，否則又會(huì)存在收斂慢的問題，我所見到的Ks在1e7-1e9之間5、就是記錄界面上的應(yīng)力隨時(shí)間的變化該用哪個(gè)命令我看例子上用hist interface i sstress x y zL