首页 岩土在线 评论员 岩土论坛 注册岩土 工程市场 岩土问答 岩土人才网 岩土领域 岩土企业
论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好！ 欢迎 登录 或 注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
热搜： 注册岩土规范勘察基坑
发帖 回复
返回列表
  • 3181阅读
  • 3回复

【求助】修正剑桥模型无法计算 总是报错 求指点 [复制链接]

 上一主题 下一主题
离线johnnyliu924
助理工程师
 

发帖
100
土币
484
威望
26
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2014-12-28
在算一个隧道开挖的模型，周边有桩与污水管
在地应力平衡时总是出错，调了快一周了都没知道如何解决，求坛子里的大神帮助！谢谢

在对管线和桩先行用model ela时 总是提示：
property ‘cam_cp’ not found in zone XXX（均是污水管或桩在模型边缘的单元）
于是先整个模型都用model cam-clay
但总是报错“Cam-clay:  Bulk modulus bound must be increased for stability”
我查了一下模型的源代码，如下
  /* --- bulk modulus --- */
         dBulk  = dMV * ps->working[Pav] / dKappa;
         if (dBulk > dBulkB) return("Cam-clay:  Bulk modulus bound must be increased for stability");
但我无论把bulk_bound怎么调，都还是报错
下面是我的地应力平衡代码：
复制代码
  1. rest 建模（含接触面）.sav
  4. ;—————————————以下为初平运算部分————————————
  6. ;设置模型参数
  8. model cam-clay
  9. prop shear 150000 bulk_bound 20e6
  10. prop mm 0.73 lambda 0.0938 kappa 0.0234
  11. prop mpc 0.395e6 mp1 1.0 mv_l 3.32
  12. ini dens 1870
  14. ;接触面参数赋值
  15. int 1 prop ks 4.52e7 kn 4.52e7 coh 4.96e3 fric 24
  16. ;int 2 prop ks 4.52e5 kn 4.52e5 coh 4.96e3 fric 24
  17. int 3 prop ks 4.52e6 kn 4.52e6 coh 4.96e2 fric 24
  19. ini szz 0 grad 0 0 18.7e3 ran z 0 -30
  21. ini sxx 0 grad 0 0 9.35e3 ran z 0 -30
  23. ini syy 0 grad 0 0 9.35e3 ran z 0 -30
  26. ;加位移边界条件
  27. fix x ran x -5.1 -4.9
  28. fix x ran x 54.9 55.1
  29. fix y ran y -0.1 0.1
  30. fix y ran y 39.9 40.1
  31. fix x y z ran z -29.9 -30.1
  34. ;———————-设置修正剑桥模型参数——————
  36. def camclay_ini_p
  37. pnt=zone_head
  38. loop while pnt # null
  39. OCR=1.2
  40. s1=-z_sxx(pnt)
  41. s2=-z_syy(pnt)
  42. s3=-z_szz(pnt)
  44. p0=(s1+s2+s3)/3.0
  46. z_prop(pnt,'cam_cp')=p0
  48. q0=sqrt(((s1-s2)*(s1-s2)+(s2-s3)*(s2-s3)+(s3-s1)*(s3-s1))*0.5)
  50. temp1=q0/(z_prop(pnt,'mm')*p0)
  52. pc=p0*(1.0+temp1*temp1)*OCR
  54. z_prop(pnt,'mpc')=pc
  56. pnt=z_next(pnt)
  58. endloop
  59. end
  61. camclay_ini_p
  63. ;————————————————————————————
  65. model ela ran group 桩
  66. prop bulk 2.4e9 shear 5e9 ran group 桩
  67. ini density 2410 ran group 桩
  69. model ela ran group 污水管
  70. prop bulk 2.4e8 shear 5e8 ran group 污水管
  71. ini density 2410 ran group 污水管
  73. set gravity 0 0 -10
  75. hist unbal
  76. set hist_rep 1
  78. set mech ratio 1e-7
  80. plot set rot 30 0 40
  81. plot con szz inter 1.5e5 ou on
  82. plot add hist 1
  84. set movie avi step 5 file 初平.avi
  85. movie start
  87. solve
  89. save 初平.sav
  91. call 隧道开挖（新）.txt
回复
举报
分享到 淘江湖 新浪 QQ微博 QQ空间 开心 人人 豆瓣 网易微博 百度 鲜果 白社会 飞信
离线johnnyliu924
助理工程师

发帖
100
土币
484
威望
26
原创币
0
只看该作者 1楼 发表于: 2014-12-28
参数我均是参考陈育民老师书中最后一张的隧道开挖的修正剑桥模型参数
不知道到底哪里出错了 求大神指点
回复
举报
离线johnnyliu924
助理工程师

发帖
100
土币
484
威望
26
原创币
0
只看该作者 2楼 发表于: 2015-01-04
求大神帮助啊 毕业的东西 搞不出来就over了
回复
举报
离线科费卡
实习生

发帖
18
土币
81
威望
4
原创币
0
只看该作者 3楼 发表于: 2016-08-01
请问大神解决了没？
回复
举报
发帖 回复
返回列表
快速回复
限100 字节
温馨提示：欢迎交流讨论，请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站