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 notsave 1.sav（1）將模型文件存為1.sav；（2）將除了需要加接觸面的單元以外的所有單元?jiǎng)h除，僅保存group 2命令為： del ran grop2 not ；（3

3、）在 group 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.sa

4、v（5）重新restore 先前的模型，并將需要建立接觸面的實(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

5、10e3 f 15app nstr -200e3 ran x 0 1 y 1 2 z 3solve可以看到接觸面的接應(yīng)力分布結(jié)果：得到的結(jié)果還是蠻理想的3、用gen separate group1int 1 wrap group1 group2此方法是在已建好的模型上使用，先把group1與周圍組分開(kāi)，這個(gè)組的周圍點(diǎn)雖說(shuō)和周圍相鄰的周圍組的點(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

6、.savdel ran 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.savimpg

7、rid 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ù)試驗(yàn)確定,一般可取接觸面兩側(cè)彈性模量較小值的10 倍（手冊(cè)上所說(shuō)）這個(gè)參數(shù)取值確實(shí)沒(méi)有太多經(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 l5 l# b K&B e/ C. jKn 法向剛度通常設(shè)置一個(gè)較大的值，但是這個(gè)值不能太大，否則又會(huì)存在收斂慢的問(wèn)題，我所見(jiàn)到的 Ks 在 1e7-1e9 之間5、就是記錄界面上的應(yīng)力隨時(shí)間的變化該用哪個(gè)命令我看例子上用hist interface i sstress x y zLakewa