›› 2006, Vol. 27 ›› Issue (4): 525-529.

• 基础理论与实验研究 • 上一篇    下一篇

岩土介质三维快速拉格朗日数值分析方法研究

刘建华,朱维申,李术才   

  1. 山东大学 岩土与结构工程研究中心,济南 250061
  • 收稿日期:2004-09-27 出版日期:2006-04-10 发布日期:2013-11-05
  • 作者简介:刘建华,男,1965年生,副教授,在职博士研究生,主要从事岩土与水利专业教学与科研工作

Study on fast lagrangian numerical analysis method for rock and soil medium in 3 dimensions

LIU Jian-hua, ZHU Wei-shen, LI Shu-cai   

  1. Geotechnical and Structural Engineering Research Center, Shandong University, Jinan 250061, China
  • Received:2004-09-27 Online:2006-04-10 Published:2013-11-05

摘要: 研究了FLAC3D(fast lagrangian analysis of continua in 3 dimensions)的特点,并与有限单元法作了比较。FLAC3D方法以结点运动方程为支配方程,追踪了介质从受荷到达到平衡状态的过程,而有限元法是根据介质力学平衡方程直接求解,这是二者主要区别。FLAC3D没有采用介质真实的阻尼特性和结点质量,给出的不是介质所经历的真实过程,不能正确反映过程的影响,因此给出的介质应力和变形计算结果的物理意义是不甚明确的。求解过程中的介质振动,是一种噪音,可引起弹塑性介质计算结果误差,而弹性介质的计算结果几乎不受影响。研究指出了FLAC3D方法的优缺点。还通过算例作了四个方面的研究:(1) Drucker-Prager屈服准则与Mohr-Coulomb屈服准则比较;(2) 膨胀角取值对计算结果的影响;(3) 大变形与小变形对计算结果的影响;(4) 精度设置对计算结果的影响。研究表明,Drucker-Prager准则与Mohr-Coulomb准则结果差异颇大;膨胀角取值对结果的影响是敏感和显著的;一般情况下,取小变形模式是合适的,计算精度取10-5是足够的。

关键词: 岩土介质, 拉格朗日数值分析方法, 有限单元法, 屈服准则, 膨胀角

Abstract: The characters of FLAC3D are studied and compared with finite element method. FLAC3D method uses node motion equation and traces the medium motion process from being loaded to reaching equilibrium state. While the finite element method solves from mechanical equilibrium equation directly. The above is the basic difference between the two methods. Because FLAC3D does not use the real medium damping property and node mass, the real process experiencing by the medium is not given and the influence of medium motion process on result cannot be reflected rightly. So the physical meaning of stress and deformation given by FLAC3D is not very clear. The medium vibration during solving process is a noise rather than a more real simulation to medium mechanical behavior and could cause calculation error for elastoplastic medium. The study shows that FLAC3D has a very high precision in determining the stress and deformation of elastic medium at equilibrium state and the vibration in solving process almost does not affect the final result. The advantages of FLAC3D method are that its mathematical operation is simple, and the solving process is convergent, and it is convenient in treating cases of medium large deformation, elastoplastic property, non-associated flow rule, excavation and supporting, etc. compared with finite element method. The defects of FLAC3D are that its calculation time is long, and the efficiency is low, and in some cases the number of time steps is high amazingly which may result in accumulating error. Studies in four aspects are also made through examples. (1) Comparisons between Drucker-Prager yield criterion and Mohr-Coulomb yield criterion. (2) Influence of dilation angle value on calculation result. (3) Influence of large deformation mode and small deformation mode on calculation result. (4) Influence of calculation precision setting on result.

Key words: rock and soil medium, Lagrangian numerical analysis method, finite element method, yield criterion, dilation angle

中图分类号: 

  • TU 452
[1] 刘 扬,杨 刚,王军祥,姜谙男,. 岩石摩尔-库仑弹塑性损伤本构模型及其主应力隐式返回映射算法研究[J]. , 2017, 38(S1): 418-428.
[2] 罗 滔,E. T. Ooi,A. H. C Chan,傅少君,. 一种模拟堆石料颗粒破碎的离散元-比例边界有限元结合法[J]. , 2017, 38(5): 1463-1471.
[3] 张 锟,徐 青,王一凡,阿湖宝. 自适应差分进化算法在边坡滑面搜索中的应用[J]. , 2017, 38(5): 1503-1509.
[4] 管华栋,周晓敏. 基于围岩相互作用的冻结壁弹塑性分析对比研究[J]. , 2017, 38(3): 649-655.
[5] 刘新荣,王子娟,傅 晏,张 梁,袁 文,缪露莉, . 考虑干湿循环作用泥质砂岩的强度与破坏准则研究[J]. , 2017, 38(12): 3395-3401.
[6] 蒋 欢 ,王水林 ,王万军,. 一种求解应变软化岩体中球腔问题的数值方法[J]. , 2016, 37(S2): 697-705.
[7] 周 辉 ,李 震 ,朱国金 ,李 云 ,卢景景 ,张冬冬,. 基于岩石统一能量屈服准则的硬岩损伤模型[J]. , 2016, 37(3): 609-615.
[8] 周 辉,卢景景,徐荣超,张传庆,陈 王君,孟凡震. 硬脆性大理岩拉剪破坏特征与屈服准则研究[J]. , 2016, 37(2): 305-314.
[9] 金晓飞 ,梁书亭 ,朱筱俊 ,张玉良,. 软土地区深基坑连续墙最大变形简易计算方法[J]. , 2015, 36(S2): 583-587.
[10] 苏 辉,杨石飞,顾国荣. 基于旁压试验土体弹塑性本构模型初探[J]. , 2015, 36(S1): 131-136.
[11] 刘远锋 ,郑 刚,. 地基排水控制法的沉降和孔压数值分析[J]. , 2015, 36(S1): 592-596.
[12] 蔡燕燕 ,张建智 ,俞 缙 ,陈士海,. 考虑围岩蠕变全过程与扩容的深埋隧洞非线性位移解[J]. , 2015, 36(7): 1831-1839.
[13] 吴顺川 ,耿晓杰 ,高永涛 ,赵国军 ,李 健 ,严 琼,. 基于广义Hoek-Brown准则的隧道纵向变形曲线研究[J]. , 2015, 36(4): 946-952.
[14] 王军祥 ,姜谙男 , . 岩石弹塑性损伤本构模型建立及在隧道工程中的应用[J]. , 2015, 36(4): 1147-1158.
[15] 赵明华,张 锐,刘 猛. 下限分析有限单元法的非线性规划求解[J]. , 2015, 36(12): 3589-3597.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[2] 闻世强,陈育民,丁选明,左威龙. 路堤下浆固碎石桩复合地基现场试验研究[J]. , 2010, 31(5): 1559 -1563 .
[3] 张常光,张庆贺,赵均海. 非饱和土抗剪强度及土压力统一解[J]. , 2010, 31(6): 1871 -1876 .
[4] 杨天鸿,陈仕阔,朱万成,刘洪磊,霍中刚,姜文忠. 煤层瓦斯卸压抽放动态过程的气-固耦合模型研究[J]. , 2010, 31(7): 2247 -2252 .
[5] 郝冬雪,陈 榕,栾茂田,武科. SBPT测定饱和黏土不排水强度的数值分析[J]. , 2010, 31(7): 2324 -2328 .
[6] 胡秀宏,伍法权. 岩体结构面间距的双参数负指数分布研究[J]. , 2009, 30(8): 2353 -2358 .
[7] 邴 慧 ,何 平. 不同冻结方式下盐渍土水盐重分布规律的试验研究[J]. , 2011, 32(8): 2307 -2312 .
[8] 李卫超,熊巨华,杨 敏. 分层土中水泥土围护结构抗倾覆验算方法的改进[J]. , 2011, 32(8): 2435 -2440 .
[9] 范书立 ,陈健云 ,张俊清. 波浪荷载作用下斜向抗拔桩的承载特性分析[J]. , 2012, 33(1): 301 -306 .
[10] 查甫生,许 龙,崔可锐. 水泥固化重金属污染土的强度特性试验研究[J]. , 2012, 33(3): 652 -658 .