首页/文章/ 详情

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

4月前浏览3295
摘要:本文采用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-normal8e8 coupling-stiffness-shear 8e8 coupling-yield-normal0.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
        最近编辑:4月前
        FLAC3D小技巧
        硕士 专注FLAC3D中的小技巧分享...
        获赞 24粉丝 117文章 40课程 0
        点赞
        收藏
        未登录
        还没有评论

        课程
        培训
        服务
        行家

        VIP会员 学习 福利任务 兑换礼品
        下载APP
        联系我们
        帮助与反馈