FLAC3D接触面建立方法

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