›› 2015, Vol. 36 ›› Issue (2): 393-402.doi: 10.16285/j.rsm.2015.02.014
王军祥1, 2,姜谙男1, 2
WANG Jun-xiang1, 2, JIANG An-nan1, 2
摘要: 针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager(D-P)强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法(NR-AL法)联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明:所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。
中图分类号:
[1] | 程昊, 唐辉明, 吴琼, 雷国平. 一种考虑水力滞回效应的非饱和土弹塑性扩展 剑桥本构模型显式算法有限元实现[J]. 岩土力学, 2020, 41(2): 676-686. |
[2] | 金青, 王艺霖, 崔新壮, 王成军 , 张珂, 刘正银, . 拉拔作用下土工合成材料在风化料-废弃轮胎 橡胶颗粒轻质土中的变形行为研究[J]. 岩土力学, 2020, 41(2): 408-418. |
[3] | 邓子千, 陈嘉帅, 王建伟, 刘小文, . 基于SFG模型的统一屈服面本构模型与试验研究[J]. 岩土力学, 2020, 41(2): 527-534. |
[4] | 李潇旋, 李涛, 彭丽云, . 控制吸力循环荷载下非饱和黏性土 的弹塑性双面模型[J]. 岩土力学, 2020, 41(2): 552-560. |
[5] | 何鹏飞, 马巍, 穆彦虎, 黄永庭, 董建华, . 黄土−砂浆块界面剪切特性试验及本构模型研究[J]. 岩土力学, 2019, 40(S1): 82-90. |
[6] | 王伟, 陈国庆, 郑水全, 张广泽, 王栋, . 考虑张拉-剪切渐进破坏的边坡矢量和法研究[J]. 岩土力学, 2019, 40(S1): 468-476. |
[7] | 刘斯宏, 沈超敏, 毛航宇, 孙 屹. 堆石料状态相关弹塑性本构模型[J]. 岩土力学, 2019, 40(8): 2891-2898. |
[8] | 张超, 杨期君, 曹文贵. 考虑峰值后区应力跌落速率的 脆岩损伤本构模型研究[J]. 岩土力学, 2019, 40(8): 3099-3106. |
[9] | 张凌凯, 王睿, 张建民, 唐新军, . 考虑颗粒破碎效应的堆石料静动力本构模型[J]. 岩土力学, 2019, 40(7): 2547-2554. |
[10] | 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246. |
[11] | 汪俊敏, 熊勇林, 杨骐莱, 桑琴扬, 黄强. 不饱和土动弹塑性本构模型研究[J]. 岩土力学, 2019, 40(6): 2323-2331. |
[12] | 汪 杰, 宋卫东, 谭玉叶, 付建新, 曹 帅, . 水平分层胶结充填体损伤本构模型及强度准则[J]. 岩土力学, 2019, 40(5): 1731-1739. |
[13] | 孙逸飞, 陈 成, . 无状态变量的状态依赖剪胀方程及其本构模型[J]. 岩土力学, 2019, 40(5): 1813-1822. |
[14] | 徐 鹏, 蒋关鲁, 雷 涛, 刘 琪, 王智猛, 刘 勇, . 考虑填土强度的加筋土挡墙动位移计算[J]. 岩土力学, 2019, 40(5): 1841-1846. |
[15] | 杨骐莱, 熊勇林, 张 升, 刘干斌, 郑荣跃, 张 锋, . 考虑温度影响的软岩弹塑性本构模型[J]. 岩土力学, 2019, 40(5): 1898-1906. |
|