FLac3D接触面简单实例

有接触面模型

gen zone brick p0 0 0 0 p1 5 0 0 p2 0 5 0 p3 0 0 2 size 10 10 4 group 1

gen zone brick p0 0 0 2 p1 5 0 2 p2 0 5 2 p3 0 0 4 size 10 10 4 group 2

gen zone brick p0 0 0 4 p1 5 0 4 p2 0 5 4 p3 0 0 6 size 10 10 4 group 3

gen separate 2

int 1 wrap 1 2

gen separate 3

int 2 wrap 2 3

plot block group int red

fix z range z -0.1 0.1

fix x range x -0.1 0.1

fix x range x 4.9 5.1

fix y range y -0.1 0.1

fix y range y 4.9 5.1

model mohr

prop bulk 2.4e9 shear 1.5e9 coh 2.5e6 fri 32 tension 1.2e6 ran group 1

prop bulk 2.4e9 shear 1.5e9 coh 2.5e6 fri 32 tension 1.2e6 ran group 2

prop bulk 2.4e9 shear 1.5e9 coh 2.5e6 fri 32 tension 1.2e6 ran group 3

ini dens 2660 ran group 1

ini dens 2660 ran group 2

ini dens 2660 ran group 3

inter 1 prop kn 0.88e11 ks 0.88e11 fri 15 coh 2e6 ten 1e6

inter 2 prop kn 0.88e11 ks 0.88e11 fri 15 coh 2e6 ten 1e6

set grav 0 0 -10

ini szz -7e6 range z 0 6

ini sxx -8.4e6 range z 0 6

ini syy -8.4e6 range z 0 6

app nstress -7e6 range z 5.9 6.1 x 0 5 y 0 5

plot con zd

plot add int black

solve

无接触面对比模型

gen zone brick p0 0 0 0 p1 5 0 0 p2 0 5 0 p3 0 0 2 size 10 10 4 group 1

gen zone brick p0 0 0 2 p1 5 0 2 p2 0 5 2 p3 0 0 4 size 10 10 4 group 2

gen zone brick p0 0 0 4 p1 5 0 4 p2 0 5 4 p3 0 0 6 size 10 10 4 group 3

plot block group int red

fix z range z -0.1 0.1

fix x range x -0.1 0.1

fix x range x 4.9 5.1

fix y range y -0.1 0.1

fix y range y 4.9 5.1

model mohr

prop bulk 2.4e10 shear 1.5e10 coh 2.5e10 fri 32 tension 1.2e10 ran group 1

prop bulk 2.4e10 shear 1.5e10 coh 2.5e10 fri 32 tension 1.2e10 ran group 2

prop bulk 2.4e10 shear 1.5e10 coh 2.5e10 fri 32 tension 1.2e10 ran group 3

ini dens 2660 ran group 1

ini dens 2660 ran group 2

ini dens 2660 ran group 3

set grav 0 0 -10

ini szz -7e6 range z 0 6

ini sxx -8.4e6 range z 0 6

ini syy -8.4e6 range z 0 6

app nstress -7e6 range z 5.9 6.1 x 0 5 y 0 5

plot con zd

plot add int black

solve