案例三-接触面（均匀网格）

一、移来移去法：首先将需要建立接触面的两个网格分开建立，然后在一个网格的指定位置建立接触面，然后把另外一部分网格移到特定的位置；

缺点：在已完成的模型中显然不奏效；

zone create radial-cylinder point 0 (0,0,0) point 1 (8,0,0) ...
                            point 2 (0,0,-5) point 3 (0,8,0) ...
                            point 4 (8,0,-5) point 5 (0,8,-5) ...
                            point 6 (8,8,0) point 7 (8,8,-5) ...
                            point 8 (0.3,0,0) point 9 (0,0.3,0) ...
                            point 10 (0.3,0,-5) point 11 (0,0.3,-5) ...
                            size 3 10 6 15 ratio 1 1 1 1.15
zone create radial-cylinder point 0(0,0,-5) point 1 (8,0,-5) ...
                            point 2 (0,0,-8) point 3 (0 8 -5) ...
                            point 4 (8,0,-8) point 5 (0,8,-8) ...
                            point 6 (8,8,-5) point 7 (8,8,-8) ...
                            point 8 (0.3,0,-5) point 9 (0,0.3,-5) ...
                            point 10 (0.3,0,-8) point 11 (0,0.3,-8) ...
                            size 3 6 6 15 ratio 1 1 1 1.15 fill on
zone reflect dip-direction 270 dip 90
zone group 'clay'
;
zone interface '1' create by-face range cylinder end-1 0 0 0 end-2 0 0 -5.1 radius 0.31 ...
                                        cylinder end-1 0 0 0 end-2 0 0 -5.1 radius 0.29 not
zone interface '2' create by-face range cylinder end-1 0 0 -4.9 end-2 0 0 -5.1 radius 0.31
;
pause
zone create cylinder point 0 (0,0,6) point 1 (0.3,0,6) point 2 (0,0,1) ...
                     point 3 (0,.3,6) point 4 (0.3,0,1) point 5 (0,0.3,1) size 3 10 6
zone create cylinder point 0 (0,0,6.1) point 1 (0.3,0,6.1) point 2 (0,0,6) ...
                     point 3 (0,0.3,6.1) point 4 (0.3,0,6) point 5 (0,0.3,6) size 3 1 6
zone reflect dip-direction 270 dip 90 range position-z 1 6.1
zone group 'pile' range position-z 1 6.1
pause
zone gridpoint initialize position-z -6.0 add range group 'pile'
model save 'pile_geom.sav'

二、导来导去法：移来移去法在已经完成的模型中显然不能奏效，同时因为模型通常都较为复杂，建模本身花费的精力较大，如果修改模型的话，就会造成大量不必要的重复劳动；

导来导去法能解决上述问题

;;;;;;;;;;;;;;;;;建模&分组;;;;;;;;;;;;;;;;;;
model new
zone create brick size 3 3 3
zone group "2" range position-x 0 2 position-y 0 2 position-z 0 2
zone group "1" range group "2" not
model save '1.sav'

;;;;;;;;;;;;;;;;;删除接触面以外的实体&建立接触面;;;;;;;;;;;;;;;;;;
zone delete range group '1'
zone interface 'jiechumian' create by-face range group "2"
zone export 'test2-jiechumian.f3grid'
model save '2.sav'

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
model restore '1.sav'
zone delete range group '2'
zone import "test2-jiechumian.f3grid"
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

zone cmodel assign elastic
zone property bulk 3e7 shear 1e7 density 2000
zone face apply velocity-x 0 range position-z 0
zone face apply velocity-y 0 range position-z 0
zone face apply velocity-z 0 range position-z 0
model gravity 0 0 -10
;zone initialize-stresses ratio 0.3
zone interface 'jiechumian' node property stiffness-normal 20e6 stiffness-shear 20e6 cohesion 10e3 friction 15
zone apply force -200e3 range position-z 3; position-x 0 1 position-y 1 2
model solve

三、切割模型法：此方法是在已建好的模型上使用，先把group1与周围组分开，这个组的周围点虽说和周围相邻的组的点有相同的坐标，但是有不同的节点。然后在group1上建立接触面；

model new
zone create brick size 3 3 3 group '2'
zone group '1' range position-x 1 2 position-y 1 2 position-z 2 3
zone densify
zone interface '1' create by-face separate range group '1' group '2'

---

---

---

均匀接触：zone interface '1' create by-face separate range group '1' group '2'

非均匀接触：zone interface '1' create by-face range position-z 1.1 4.5  position-x 0.1 2.9 position-y 0.1 2.9

注：用空间范围约定接触面建立的区域，接触面只能附在两个分组构建的表面，分组内部无法建立接触面

zone interface '1' node  property stiffness-normal ** stiffness-shear **  cohesion ** friction angle ** dilation angle ** tensile strength **

model new
model title "Real interface - slip and separation only"
fish automatic-create off
; Create Material Zones
zone create brick size 5 5 5 ...
         point 0 (0,0,0) point 1 (3,0,0) point 2 (0,3,0) point 3 (0,0,5) ...
         point 4 (3,3,0) point 5 (0,5,5) point 6 (5,0,5) point 7 (5,5,5)
zone create brick size 5 5 5 point 0 (0,0,5) edge 5.0
zone group 'Material'
; Create Bin Zones
zone create brick size 1 5 5 ...
         point 0 (3,0,0) point 1 add (3,0,0) point 2 add (0,3,0) ...
         point 3 add (2,0,5) point 4 add (3,6,0) point 5 add (2,5,5) ...
         point 6 add (3,0,5) point 7 add (3,6,5) group 'Bin1'
zone create brick size 1 5 5 ...
         point 0 (5,0,5) point 1 add (1,0,0) point 2 add (0,5,0) ...
         point 3 add (0,0,5) point 4 add (1,6,0) point 5 add (0,5,5) ...
         point 6 add (1,0,5) point 7 add (1,6,5) group 'Bin1'
zone create brick size 5 1 5 ...
         point 0 (0,3,0) point 1 add (3,0,0) point 2 add (0,3,0) ...
         point 3 add (0,2,5) point 4 add (6,3,0) point 5 add (0,3,5) ...
         point 6 add (5,2,5) point 7 add (6,3,5) group 'Bin2'
zone create brick size 5 1 5 ...
         point 0 (0,5,5) point 1 add (5,0,0) point 2 add (0,1,0) ...
         point 3 add (0,0,5) point 4 add (6,1,0) point 5 add (0,1,5) ...
         point 6 add (5,0,5) point 7 add (6,1,5) group 'Bin2'
; Assign models to groups
zone cmodel assign mohr-coulomb range group 'Material'
zone cmodel assign elastic range group 'Bin1' or 'Bin2'
; Separate zones, and create two interfaces
;zone separate by-face new-side group 'Interface' range group 'Material' group 'Bin1' or 'Bin2'
;zone interface 1 create by-face range group 'Interface' group 'Bin1'
;zone interface 2 create by-face range group 'Interface' group 'Bin2'
zone interface 1 create by-face separate range group 'Material' group 'Bin1'
zone interface 2 create by-face separate range group 'Material' group 'Bin2'
; Subdivide interface elements

; Assign properties
zone property shear 1e8 bulk 2e8 fric 30 range group 'Material'
zone property shear 1e8 bulk 2e8         range group 'Bin1' or 'Bin2'
zone property density 2000
zone interface 1 node property stiffness-shear 2e9 bonded-slip yes ...
                               stiffness-normal 2e9 friction 15
zone interface 2 node property stiffness-shear 2e9 bonded-slip yes ...
                               stiffness-normal 2e9 friction 15
; Assign Boundary Conditions
zone gridpoint fix velocity-x range union position-x 0 position-x 6
zone gridpoint fix velocity-y range union position-y 0 position-y 6
zone gridpoint fix velocity-z range position-z 0 group 'Bin1' or 'Bin2'
; Monitor histories
zone history displacement-z position (6,6,10)
zone history displacement-z position (0,0,10)
zone history displacement-z position (0,0,0)
; Settings
model largestrain on
model gravity 10
; Cycling
model step 4000
model save '20-5'

Question：

;zone separate by-face new-side group 'Interface' range group 'Material' group 'Bin1' or 'Bin2'
;zone interface 1 create by-face range group 'Interface' group 'Bin1'
;zone interface 2 create by-face range group 'Interface' group 'Bin2'
zone interface 1 create by-face separate range group 'Material' group 'Bin1'
zone interface 2 create by-face separate range group 'Material' group 'Bin2'

bonded-slip yes

zone interface 1 node property stiffness-shear 2e9 bonded-slip yes ...
                               stiffness-normal 2e9 friction 15

1