首页/文章/ 详情

【复刻】FLAC3D 6.0创建非全长锚固预紧力锚杆(索)

11月前浏览7598
摘要:本文采用FLAC3D 6.0复刻了《FLAC、FLAC3D基础与工程实例》(陈育民,第二版)中关于非全长锚固预紧力锚杆(索)的建立方法,背景介绍见书。

1. 通过删除-建立link连接来模拟托盘



























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
2. 通过设置极大锚固剂参数模拟托盘






















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
3. 借助别的结构单元(如liner单元)来模拟托盘


























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

图1 运行结果图

4. 总结

以上三种方法的模拟效果基本相同,陈老师觉得第一种方法较为简单、方便,我觉得还是第二种方法更简单方便,第三种最麻烦。

来源:FLAC3D小技巧
FLAC3DGIDOrigin
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2023-12-26
最近编辑:11月前
FLAC3D小技巧
硕士 专注FLAC3D中的小技巧分享...
获赞 36粉丝 209文章 40课程 0
点赞
收藏
未登录
1条评论
戒
签名征集中
6月前
你好,我想咨询一下,我删除了component-id 17 的link 当我重新创建18节点到19节点的link时,他显示link已存在,也就是前面structure link delete range component-id 17命令没有成功,这时为什么?
回复 1条回复
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈