model new
zone create radial-tunnel point 0 0 0 0 point 1 25 0 0 point 2 0 50 0 point 3 0 0 25 size 4 25 4 10 dimension 4 4 4 4 ratio 1 1 1 1.1 fill on
zone reflect normal 1 0 0 origin 0 0 0
zone reflect normal 0 0 1 origin 0 0 0
zone cmodel assign mohr-coulomb
zone property bulk 2.2e9 shear 1.3e9 friction 30 cohesion 1.3e6 tension 1.5e5
zone initialize density 2000
;
zone gridpoint fix velocity-x range position-x -25.1 -24.9
zone gridpoint fix velocity-x range position-x 24.9 25.1
zone gridpoint fix velocity-y range position-y 49.9 50.1
zone gridpoint fix velocity-z range position-z -25.1 -24.9
zone gridpoint fix velocity-z range position-z 24.9 25.1
;
structure cable create by-line 0 0 0 0 29 0 segments 10 id 1
structure cable create by-line 0 29 0 0 35 0 segments 6 id 1
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 1 grout-stiffness 1 grout-perimeter 0.0785 range component-id 1 10 id 1
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e5 grout-stiffness 2e7 range component-id 11 17 id 1
;删除端部的link
structure link delete range component-id 1
;在原处建立link并分组为’end'
structure link create on-nodeid 1 target zone group 'end'
;设置link的自由度值
structure link attach x rigid y rigid z rigid rotation-x rigid rotation-y rigid rotation-z rigid range group 'end'
;施加预应力
struct cable apply tension value 60e3 range component-id 1 10
model cycle 2000
model new
zone create radial-tunnel point 0 0 0 0 point 1 25 0 0 point 2 0 50 0 point 3 0 0 25 size 4 25 4 10 dimension 4 4 4 4 ratio 1 1 1 1.1 fill on
zone reflect normal 1 0 0 origin 0 0 0
zone reflect normal 0 0 1 origin 0 0 0
zone cmodel assign mohr-coulomb
zone property bulk 2.2e9 shear 1.3e9 friction 30 cohesion 1.3e6 tension 1.5e5
zone initialize density 2000
;
zone gridpoint fix velocity-x range position-x -25.1 -24.9
zone gridpoint fix velocity-x range position-x 24.9 25.1
zone gridpoint fix velocity-y range position-y 49.9 50.1
zone gridpoint fix velocity-z range position-z -25.1 -24.9
zone gridpoint fix velocity-z range position-z 24.9 25.1
;
structure cable create by-line 0 0 0 0 29 0 segments 10 id 1
structure cable create by-line 0 29 0 0 35 0 segments 6 id 1
;此处注意不同component的锚固参数值
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 1 grout-stiffness 1 grout-perimeter 0.0785 range component-id 2 10
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e5 grout-stiffness 2e7 range component-id 11 17
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e8 grout-stiffness 2e10 range component-id 1 1
struct cable apply tension value 60e3 range component-id 1 10
model cycle 2000
model new
zone create radial-tunnel point 0 0 0 0 point 1 25 0 0 point 2 0 50 0 point 3 0 0 25 size 4 25 4 10 dimension 4 4 4 4 ratio 1 1 1 1.1 fill on
zone reflect normal 1 0 0 origin 0 0 0
zone reflect normal 0 0 1 origin 0 0 0
zone cmodel assign mohr-coulomb
zone property bulk 2.2e9 shear 1.3e9 friction 30 cohesion 1.3e6 tension 1.5e5
zone initialize density 2000
;
zone gridpoint fix velocity-x range position-x -25.1 -24.9
zone gridpoint fix velocity-x range position-x 24.9 25.1
zone gridpoint fix velocity-y range position-y 49.9 50.1
zone gridpoint fix velocity-z range position-z -25.1 -24.9
zone gridpoint fix velocity-z range position-z 24.9 25.1
;
structure cable create by-line 0 0 0 0 29 0 segments 10 id 1
structure cable create by-line 0 29 0 0 35 0 segments 6 id 1
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 1 grout-stiffness 1 grout-perimeter 0.0785 range component-id 1 10 id 1
structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e5 grout-stiffness 2e7 range component-id 11 17 id 1
structure liner create by-face range position-y -.1 .1 position-x -1 1 position-z -1 1
structure liner property isotropic 25e9 0.15 thickness 0.1
structure liner property coupling-stiffness-normal 8e8 coupling-stiffness-shear 8e8 coupling-yield-normal 0.0 coupling-cohesion-shear 0.0 coupling-cohesion-shear-residual 0.0 coupling-friction-shear 0.0
structure link delete range component-id 1
structure link create on-nodeid 1 target node 18 group 'end'
structure link attach x rigid y rigid z rigid rotation-x rigid rotation-y rigid rotation-z rigid range group 'end'
struct cable apply tension value 60e3 range component-id 1 10
model cycle 2000
图1 运行结果图
以上三种方法的模拟效果基本相同,陈老师觉得第一种方法较为简单、方便,我觉得还是第二种方法更简单方便,第三种最麻烦。