- 积分
- 18
- 注册时间
- 2009-11-29
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 bananaliuchao 于 2011-11-10 22:41 编辑
有日子没发技术性主题了,因此积分也少得可怜,作为一个aba版斑竹,我承认我很懒,基础很不扎实。
今天要和坛油探讨的一个问题是隧道施工模拟的一种传统的方法——应力释放法。
应力释放法模拟隧道开挖主要有两种实现途径,大家可以参考帖子:
http://forum.simwe.com/thread-765338-1-1.html
简单的说,主要就是参数弱化法和施加节点反力法。这两种方法各有优劣。
参数弱化法:
优点:1、便于实现;2、便于控制应力释放过程;
缺点:1、无法给出明确的应力释放率;2、无法用在剑桥模型中,因为找不到一个合适的参数进行弱化。
施加节点反力法:
优点:1、能够明确给出应力释放率;2、物理意义较为明确;
缺点:实现过程不是很方便
下面跟大家汇报一下我的做法,希望有经验的童鞋给出建议,各种板砖,放心拍吧~
一、参数弱化法
参数弱化方法其实没有什么难度,版内已经讨论了无数次了。各种牛们也是不厌其烦的跟新手们说了一遍又一遍,可惜有些童鞋天生喜欢求助,上来就发帖问,完全置搜索功能于不顾。所以我在这要再次呼吁,搜索处要更加醒目,加红字,加粗,加浮动,另外加旁注:亲,如果你想快速入门aba,搜索吧,这是闷骚技术人新一代的学习方法,包会哦!
参数弱化,主要就是弱化弹性模量,这个可以通过关键字*field来实现,具体肿么实现,搜索吧。下面给出计算后的云图,抛砖引玉。
地应力平衡,勉强可以吧:
待开挖区域弹性模量折减20%:
如何添加衬砌以及如何设置追踪单元我想再单开一贴和大家进行总结探讨。上面位移矢量图隐去了待开挖部分,是为了矢量显示清晰,实际上这一步待开挖(中心土体)部分实际还是存在的,只不过弹性模量折减过了。从图中可以看出,这种方法计算出的位移场并不是很符合实际情况。隧道开挖土体位移场可以参见帖子:http://forum.simwe.com/thread-985421-1-1.html
在这里也顺便给这个帖子做个个人的解答吧。关于该贴提出的强制位移法,我想在后面专开一贴和大家讨论总结。
有经验的童鞋可否告知你们采用折减法做出的位移场是否也和上面一样的?
二、施加节点反力法
这一方法的技术难点在于,如何获得地应力平衡后隧道周边土体节点反力。通过搜索和查看manual,基本可以锁定两种输出量:NFORC & RF。其中,RF只能在有边界约束处输出,NFORC根据manual我猜测应该是由节点所涉及单元的应力进行外插平均后得到的节点力,参见帖子:http://forum.simwe.com/thread-883192-1-1.html10#。不管怎样,总的思路是,通过对开挖后隧道相邻土体节点施加非均布节点力,使得此时的模型(土体)保持初始地应力和位移准零状态,然后施加原节点荷载*(1-应力释放率),以此来模拟不同的应力释放率下隧道及土体的力学响应。
以前的帖子对这个问题讨论了很多,有代表性的列举如下(话说论坛的收藏是个好功能):
http://forum.simwe.com/thread-872620-1-1.html
该贴lz给出了四种获取节点力的方法,都很有参考价值。最后采用了NFORC,但得到的效果不是很理想。我做了一下,发现的确如该lz所说,NFORC输出的节点力是不能直接施加到原节点上的,即使乘以2倍,也只能是接近实际的节点力。所得到的位移场在初始状态就有很大的位移,不符合前面提出的总体思路。该贴5#和7#都提出采用约束洞周位移获取RF。事实证明,这个方法是可行的,但不同的操作方式得到的结果千差万别。该贴7#提出的思路是对的,但给出的文件中将洞周位移约束贯穿了整个施工过程,这样造成的后果是,其计算起始点已经不再是green field的情况,得到的应力场自然不符合实际情况,那么在错误的应力场上做的计算怎么能对呢?
http://forum.simwe.com/thread-854591-1-1.html
该贴给出了施加洞周约束获取反力的正确思路,即施加洞周位移约束要在地应力平衡后进行。这样能保证施加位移约束时的地应力初始状态。下面就我做的例子详细说明。
模型:
材料参数:
修正剑桥模型
*Material, name=soil-3
*Clay Plasticity, intercept=1.45
0.11, 1.27, 0., 1., 1., , 1.
*Density
1770.,
*Porous Elastic
0.009, 0.32, 0.
剑桥模型参数的设置参见:http://forum.simwe.com/thread-968882-1-1.html
这里再啰嗦几句,最近很多人在问为什么用剑桥模型的时候会提示“THE SUM OF INITIAL PRESSUE AND ELASTIC TENSILE STRENGTH FOR THE POROUS ELASTIC MATERIAL IN ELEMENT...",如果你确定有地应力平衡,那么就请你检查是否是用的直接加重力平衡的吧。多孔介质弹性模型是需要在initial conditions中定义初始地应力的。
地应力平衡:
平衡效果一般,主要是做实例,只平衡了一次,到了10E-4就没再管了。
施加洞周边界约束,同时remove掉待开挖部分
可以看出,这一步进行完后,剩余土体的应力基本没有发生改变,也就是基本保持了初始地应力,在这种情况下导出节点力再倒入折减后支护力才是正确的。可以看出施加洞周位移约束后土体位移分布发生了变化,但数量级还是很小,基本没变,可以忽略。最后一张图是洞周节点反力云图,在这种情况下就可以输出节点反力了。下面介绍一下我的做法。
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
2
查看全部评分
-
|