model newzone 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 onzone reflect normal 1 0 0 origin 0 0 0zone reflect normal 0 0 1 origin 0 0 0zone cmodel assign mohr-coulombzone property bulk 2.2e9 shear 1.3e9 friction 30 cohesion 1.3e6 tension 1.5e5zone initialize density 2000;zone gridpoint fix velocity-x range position-x -25.1 -24.9zone gridpoint fix velocity-x range position-x 24.9 25.1zone gridpoint fix velocity-y range position-y 49.9 50.1zone gridpoint fix velocity-z range position-z -25.1 -24.9zone 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 1structure cable create by-line 0 29 0 0 35 0 segments 6 id 1structure 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 1structure 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;删除端部的linkstructure 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 10model 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
;删除端部的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 newzone 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 onzone reflect normal 1 0 0 origin 0 0 0zone reflect normal 0 0 1 origin 0 0 0zone cmodel assign mohr-coulombzone property bulk 2.2e9 shear 1.3e9 friction 30 cohesion 1.3e6 tension 1.5e5zone initialize density 2000;zone gridpoint fix velocity-x range position-x -25.1 -24.9zone gridpoint fix velocity-x range position-x 24.9 25.1zone gridpoint fix velocity-y range position-y 49.9 50.1zone gridpoint fix velocity-z range position-z -25.1 -24.9zone 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 1structure 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 10structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e5 grout-stiffness 2e7 range component-id 11 17structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e8 grout-stiffness 2e10 range component-id 1 1struct cable apply tension value 60e3 range component-id 1 10model cycle 2000
;此处注意不同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
model newzone 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 onzone reflect normal 1 0 0 origin 0 0 0zone reflect normal 0 0 1 origin 0 0 0zone cmodel assign mohr-coulombzone property bulk 2.2e9 shear 1.3e9 friction 30 cohesion 1.3e6 tension 1.5e5zone initialize density 2000;zone gridpoint fix velocity-x range position-x -25.1 -24.9zone gridpoint fix velocity-x range position-x 24.9 25.1zone gridpoint fix velocity-y range position-y 49.9 50.1zone gridpoint fix velocity-z range position-z -25.1 -24.9zone 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 1structure cable create by-line 0 29 0 0 35 0 segments 6 id 1structure 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 1structure cable property young 2e10 yield-tension 310e3 cross-sectional-area 0.0004906 grout-cohesion 10e5 grout-stiffness 2e7 range component-id 11 17 id 1structure liner create by-face range position-y -.1 .1 position-x -1 1 position-z -1 1structure liner property isotropic 25e9 0.15 thickness 0.1structure 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.0structure link delete range component-id 1structure 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 10model cycle 2000
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 create on-nodeid 1 target node 18 group 'end'
图1 运行结果图
以上三种方法的模拟效果基本相同,陈老师觉得第一种方法较为简单、方便,我觉得还是第二种方法更简单方便,第三种最麻烦。