基于主应力空间修正的Mohr-Coulomb准则塑性回映精确描述

1

作者 陈正峰 巨广宏 +1 位作者 刘高 符文熹 《西安建筑科技大学学报(自然科学版)》 北大核心 2025年第4期530-537,共8页

针对Mohr-Coulomb(M-C)准则未考虑岩土材料的弹脆性破坏和拉破坏,以及在三维应力空间中常常不收敛的问题.建立了带拉截断修正Mohr-Coulomb准则(修正的M-C准则),在主应力空间下,通过主应力塑性回映算法并结合屈服面的空间几何特征给出了... 针对Mohr-Coulomb(M-C)准则未考虑岩土材料的弹脆性破坏和拉破坏,以及在三维应力空间中常常不收敛的问题.建立了带拉截断修正Mohr-Coulomb准则(修正的M-C准则),在主应力空间下,通过主应力塑性回映算法并结合屈服面的空间几何特征给出了应力返回空间位置以及返回规则.依托ANSYS软件平台开发了修正Mohr-Coulomb准则用户子程序,并以深埋圆形洞室作为算例验证了提出的方法,验证了该方法的准确性. 展开更多

关键词 修正Mohr-Coulomb准则 主应力空间塑性回映算法 ANSYS 数值模拟

在线阅读 下载PDF 职称材料

基于Hoek-Brown准则的岩体高温-荷载耦合损伤模型及其验证

2

作者 姜谙男 黄炫茜 +1 位作者 许梦飞 蒋浩鹏 《沈阳工业大学学报》 CAS 北大核心 2024年第4期442-448,共7页

针对隧道岩体在高温和荷载共同作用下产生的耦合损伤问题,基于连续损伤力学理论,引入高温-荷载耦合损伤因子,建立了考虑高温作用的广义Hoek-Brown(H-B)准则的弹塑性损伤本构模型,解决了模型在数值积分过程中的奇异点问题,并利用Fortran... 针对隧道岩体在高温和荷载共同作用下产生的耦合损伤问题,基于连续损伤力学理论,引入高温-荷载耦合损伤因子,建立了考虑高温作用的广义Hoek-Brown(H-B)准则的弹塑性损伤本构模型,解决了模型在数值积分过程中的奇异点问题,并利用Fortran语言编写了模型的有限元程序。结果表明:模型计算的应力-应变曲线与试验所得到的曲线趋势一致且计算结果误差在5%以内,温度作用下隧洞围岩损伤逐渐加剧,验证了模型的准确性和程序的稳健性,可为火灾作用下隧道围岩安全稳定性评价提供一种有效计算方法。 展开更多

关键词 广义Hoek-Brown准则 耦合损伤模型 完全隐式应力返回映射算法 弹塑性模型 高温-荷载耦合 单轴压缩试验 应力-应变曲线 隧洞

在线阅读 下载PDF 职称材料

材料强化模型对板料回弹量的影响 被引量：11

3

作者 张冬娟 崔振山 +1 位作者 李玉强 阮雪榆 《上海交通大学学报》 EI CAS CSCD 北大核心 2006年第10期1671-1674,共4页

基于Lemaitre and Chaboche非线性随动强化理论、等向强化和Mises屈服准则,建立了复杂加载模式下非线性混合强化材料模型的弹塑性应力应变本构关系,并采用Backward Euler切向预测径向返回算法计算应力应变增量.基于ABAQUS开发式程序接口... 基于Lemaitre and Chaboche非线性随动强化理论、等向强化和Mises屈服准则,建立了复杂加载模式下非线性混合强化材料模型的弹塑性应力应变本构关系,并采用Backward Euler切向预测径向返回算法计算应力应变增量.基于ABAQUS开发式程序接口,编写了非线性混合强化材料模型用户子程序.以Numisheet’93板料U型弯曲考题为例,分析了不同材料强化模型对板料回弹量的影响.结果表明,线性随动强化因模拟板料成形后的应力最小而低估了回弹量,各向同性强化因模拟成形后的应力最大而使预测的回弹量偏大.与Numisheet’93实验值的比较可知,对于复杂加载问题,采用非线性混合强化材料模型预测板料回弹量的精度最高. 展开更多

关键词 非线性随动强化 屈服准则 本构关系 径向返回算法

在线阅读 下载PDF 职称材料

三维弹塑性模型在路堤软基固结分析中应用 被引量：16

4

作者 孙德安 甄文战 黄文雄 《岩土力学》 EI CAS CSCD 北大核心 2009年第3期669-674,共6页

针对基于SMP准则改进的剑桥模型(MCC-SMP),采用弹性预测和塑性修正两阶段的回映算法,给出了应力更新算法,根据其算法编制了非线性有限元子程序,进行了三轴压缩、三轴伸长及平面应变等单元试验问题的数值模拟,通过数值分析和试验结果的对... 针对基于SMP准则改进的剑桥模型(MCC-SMP),采用弹性预测和塑性修正两阶段的回映算法,给出了应力更新算法,根据其算法编制了非线性有限元子程序,进行了三轴压缩、三轴伸长及平面应变等单元试验问题的数值模拟,通过数值分析和试验结果的对比,验证了子程序的精度;采用该本构模型,对澳大利亚某试验段路堤软土地基固结进行了三维有限元分析,重点分析和对比了孔压、沉降及水平位移的变化规律。 展开更多

关键词 回映算法 三维 弹塑性模型 固结 曼德尔效应

在线阅读 下载PDF 职称材料

混凝土弹塑性损伤本构模型研究 Ⅱ:数值计算和试验验证 被引量：35

5

作者 吴建营 李杰 《土木工程学报》 EI CSCD 北大核心 2005年第9期21-27,共7页

针对建议的基于能量的混凝土弹塑性损伤本构模型[1],将损伤演化和塑性变形解耦,本文建立了模型的弹性预测-塑性修正-损伤修正数值分析框架,给出了无条件稳定的应力更新算法和相应的算法一致性切线模量。在弹性预测-塑性修正过程中,利用... 针对建议的基于能量的混凝土弹塑性损伤本构模型[1],将损伤演化和塑性变形解耦,本文建立了模型的弹性预测-塑性修正-损伤修正数值分析框架,给出了无条件稳定的应力更新算法和相应的算法一致性切线模量。在弹性预测-塑性修正过程中,利用谱分解回映算法建立了塑性流动因子和硬化参数的统一迭代格式,减小了有效应力更新的计算量。根据上述数值算法,编制非线性有限元分析程序,对单调和低周反复荷载作用下的混凝土材料和结构试验进行数值模拟,分析结果表明建议的弹塑性损伤本构模型可以较好地描述混凝土材料的典型非线性行为,其数值方法是有效的,为进一步的结构非线性分析应用奠定了基础。 展开更多

关键词 混凝土 非线性分析 本构模型 回映算法 后退欧拉法

在线阅读 下载PDF 职称材料

Hoek-Brown准则的主应力回映算法及其二次开发 被引量：7

6

作者 陈培帅 陈卫忠 +1 位作者 贾善坡 杨建平 《岩土力学》 EI CAS CSCD 北大核心 2011年第7期2211-2218,共8页

在弹塑性模型数值实现的过程中,需要进行应力更新的回映算法。针对三维应力空间回映算法在奇异点收敛性方面的不足,提出主应力空间的回映算法,讨论了算法实现过程的应力空间转化问题,分析了应力更新过程中确定回映区域的方法,建立了相... 在弹塑性模型数值实现的过程中,需要进行应力更新的回映算法。针对三维应力空间回映算法在奇异点收敛性方面的不足,提出主应力空间的回映算法,讨论了算法实现过程的应力空间转化问题,分析了应力更新过程中确定回映区域的方法,建立了相应的一致性刚度矩阵。基于大型有限元软件ABAQUS提供的用户材料子程序UMAT接口,利用Fortran编程语言,实现了Hoek-Brown准则主应力空间的隐式积分算法,利用开发的模型,进行了岩石常规三轴压缩试验的数值模拟,通过与ABAQUS内部的Mohr-Coulomb准则计算结果的对比,验证了模型和程序的可行性和准确性。 展开更多

关键词 霍克-布朗本构模型 ABAQUS~用户材料子程序 主应力空间回映算法

在线阅读 下载PDF 职称材料

考虑温度影响的混凝土弹塑性本构模型及其应用 被引量：7

7

作者 李亮 姜文华 +1 位作者 吴俊 杜修力 《应用力学学报》 CAS CSCD 北大核心 2019年第2期280-288,501-502,共10页

基于弹塑性理论的框架,引入热力学中的自由能概念,建立了考虑温度影响的混凝土弹塑性本构模型。假设混凝土为弹塑性解耦材料,通过采用包含温度的自由能函数建立了弹性应力-应变关系的表达式;基于Drucker-Prager强度准则构造屈服函数,以... 基于弹塑性理论的框架,引入热力学中的自由能概念,建立了考虑温度影响的混凝土弹塑性本构模型。假设混凝土为弹塑性解耦材料,通过采用包含温度的自由能函数建立了弹性应力-应变关系的表达式;基于Drucker-Prager强度准则构造屈服函数,以塑性剪应变作为硬化参量,通过在硬化规律中引入温度变量来反映温度对混凝土塑性性质的影响,得到了塑性应力-应变关系的表达式。基于ABAQUS有限元软件平台,应用半隐式回映算法编制了上述弹塑性本构模型的UMAT子程序,完成了该模型的程序实现。通过对C20及C40混凝土单元在不同温度条件下单轴压缩与双轴压缩的数值模拟,并将模拟结果与已有文献试验结果进行了对比,二者符合较好,验证了该模型的正确性。应用该模型模拟混凝土材料的力学特性,开展了不同温度下集中荷载作用的混凝土简支梁及混凝土轴心受压柱力学响应的数值模拟。计算结果表明,随着温度的升高,相同荷载作用下混凝土构件的力学响应逐渐增大,反映了混凝土材料的力学性能随着温度的升高逐渐退化,这符合不同温度条件下混凝土构件力学响应的基本规律。 展开更多

关键词 混凝土 弹塑性 本构模型 温度 半隐式回映算法

在线阅读 下载PDF 职称材料

不同应力路径下剪切带的数值模拟 被引量：12

8

作者 孙德安 甄文战 《岩土力学》 EI CAS CSCD 北大核心 2010年第7期2253-2258,2294,共7页

采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形... 采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明:超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。 展开更多

关键词 超固结黏土 剪胀 软化 回映算法 变形局部化 剪切带

在线阅读 下载PDF 职称材料

岩石弹塑性损伤本构模型建立及在隧道工程中的应用 被引量：25

9

作者 王军祥 姜谙男 《岩土力学》 EI CAS CSCD 北大核心 2015年第4期1147-1158,共12页

在实际隧道施工过程中,隧道开挖引起地下岩体应力重分布使得围岩的微裂纹扩展损伤,并伴随有塑性流动变形。在地下水环境中对于孔隙和微裂隙围岩介质受到应力作用时,在内部将产生高孔隙水压力影响岩石的力学性质,也改变了围岩的破坏模式... 在实际隧道施工过程中,隧道开挖引起地下岩体应力重分布使得围岩的微裂纹扩展损伤,并伴随有塑性流动变形。在地下水环境中对于孔隙和微裂隙围岩介质受到应力作用时,在内部将产生高孔隙水压力影响岩石的力学性质,也改变了围岩的破坏模式。为了研究损伤引起的刚度退化和塑性导致的流动两种破坏机制的耦合作用,从弹塑性力学和损伤理论的角度出发,同时引入修正有效应力原理来考虑孔隙水压力的作用,建立基于Drucker-Prager屈服准则的弹塑性损伤本构模型;针对该本构模型推导了孔隙水压力作用下弹塑性损伤本构模型的数值积分算法-隐式返回映射算法,分别对预测应力返回到屈服面的光滑圆锥面或尖点奇异处两种可能的情况给出了详细的描述,隐式返回映射算法具有稳定性和准确性的特点;大多数弹塑性损伤模型中涉及参数多且不易确定的问题,采用反分析方法获得损伤参数,解决了损伤参数不易确定的难题;采用面向对象的编程方法,使用C++语言编制了弹塑性损伤本构求解程序,并对所建立的弹塑性损伤模型和所编程序进行了试验和数值两个方面的验证;最后将其在吉林抚松隧道工程中进行应用,模拟了塑性区和损伤区的发展变化。研究结果表明:所建立的弹塑性损伤本构模型能够较好地描述岩石的力学性能、塑性和损伤变化趋势,所编程序能够进行实际工程问题的模拟,对现场施工给予一定的指导。 展开更多

关键词 孔隙水压力 弹塑性损伤本构模型 修正有效应力原理 Drucker-Prager屈服准则 隐式返回映射算法 反分析方法

在线阅读 下载PDF 职称材料

岩石摩尔-库仑弹塑性损伤本构模型及其主应力隐式返回映射算法研究 被引量：7

10

作者 刘扬 杨刚 +1 位作者 王军祥 姜谙男 《岩土力学》 EI CAS CSCD 北大核心 2017年第S1期418-428,共11页

在弹塑性损伤理论框架内考虑岩石的塑性变形机制和刚度退化,建立基于Mohr-Coulomb(M-C)屈服准则的弹塑性损伤模型,采用内变量即等效塑性应变表征岩石损伤变量的演化。由于M-C屈服准则在应力空间为一个六棱锥,在数值实施过程中六棱锥角... 在弹塑性损伤理论框架内考虑岩石的塑性变形机制和刚度退化,建立基于Mohr-Coulomb(M-C)屈服准则的弹塑性损伤模型,采用内变量即等效塑性应变表征岩石损伤变量的演化。由于M-C屈服准则在应力空间为一个六棱锥,在数值实施过程中六棱锥角点和棱线上的应力更新存在"奇异性"问题,角点光滑化方法可以处理该问题,但其不可避免的导致近似的计算结果。在M-C本构数值积分算法的基础上,推导弹塑性损伤本构方程的主应力空间隐式返回映射算法,包括弹性预测、塑性修正和损伤修正3个主要计算步骤。在塑性修正过程中,针对流动向量返回到主平面、左右棱线和尖点3种情况分别进行讨论,从主应力空间的角度出发解决"奇异性"问题。采用面向对象的编程方法,使用C++语言开发弹塑性损伤本构求解程序(RDM-C),并采用单轴压缩试验、地基和洞室算例对程序计算的结果进行分析和验证。研究结果表明,所建立的弹塑性损伤本构模型能够较好地描述岩石材料主要的力学和变形特性、塑性区和损伤区变化趋势。基于主应力空间的隐式积分算法所开发的程序可以进行岩土工程问题的数值分析,对现场施工提供指导和理论依据。 展开更多

关键词 岩石材料 隐式返回映射算法 主应力空间 MOHR-COULOMB屈服准则 弹塑性损伤本构模型

在线阅读 下载PDF 职称材料

岩石应变软化本构模型建立及NR-AL法求解研究 被引量：9

11

作者 王军祥 姜谙男 《岩土力学》 EI CAS CSCD 北大核心 2015年第2期393-402,共10页

针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager(D-P)强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳... 针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager(D-P)强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法(NR-AL法)联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明:所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。 展开更多

关键词 应变软化 本构模型 隐式返回映射算法 Newton-Raphson法 arc-length法 程序编制

在线阅读 下载PDF 职称材料

黏弹-双曲线Drucker-Prager塑性模型应力更新隐式算法 被引量：2

12

作者 闫富有 常键 刘忠玉 《岩土力学》 EI CAS CSCD 北大核心 2017年第6期1797-1804,共8页

黏弹-塑性模型是由黏弹性模型与塑性元件串联而成,可视为黏弹-黏塑性模型中黏塑性黏度参数趋于0的一种极限情况,为黏弹性材料的结构破坏分析提供了一个途径。把黏弹-塑性模型的应变增量分解为黏弹和塑性增量两部分,考虑黏弹性应变历史,... 黏弹-塑性模型是由黏弹性模型与塑性元件串联而成,可视为黏弹-黏塑性模型中黏塑性黏度参数趋于0的一种极限情况,为黏弹性材料的结构破坏分析提供了一个途径。把黏弹-塑性模型的应变增量分解为黏弹和塑性增量两部分,考虑黏弹性应变历史,把黏弹性积分型本构关系在一个时步内线性化,定义与时间增量相关的剪切模量和体积模量,导出应力递推公式,把黏弹-塑性本构积分转化为与弹塑性相似的形式。针对由黏弹性和双曲线Drucker-Prager塑性、各向同性硬化的黏弹-塑性模型,通过黏弹性预估和塑性校正"二步"算法实现对应力的更新,给出完全隐式算法和最终的计算公式。算例比较分析表明,由于迭代过程中仅需要简单的函数计算,该算法具有很好的收敛性。一般经过2次迭代运算后,屈服函数值已达到10-10的量级,应力点便返回到屈服面上。 展开更多

关键词 黏弹性 塑性 应力更新算法 一致切线算子 有限元法

在线阅读 下载PDF 职称材料

土体各向异性边界面模型隐式算法在ABAQUS软件中的实现 被引量：5

13

作者 钦亚洲 孙钧 《岩土力学》 EI CAS CSCD 北大核心 2012年第1期307-313,共7页

基于Wheeler土体各向异性旋转硬化法则,结合边界面理论,构造一个能够反映土体初始各向异性及加载后应力诱发各向异性的边界面本构模型,并借助ABAQUS软件提供的UMAT子程序接口,采用隐式积分算法——图形返回算法实现。通过对正常固结状态... 基于Wheeler土体各向异性旋转硬化法则,结合边界面理论,构造一个能够反映土体初始各向异性及加载后应力诱发各向异性的边界面本构模型,并借助ABAQUS软件提供的UMAT子程序接口,采用隐式积分算法——图形返回算法实现。通过对正常固结状态下(OCR=1)高岭土试样三轴不排水剪切试验进行模拟,并将模拟结果与ABAQUS自带的修正剑桥模型模拟结果进行了比较分析,表明本模型的模拟结果能够反映土体在偏压加载过程中产生的各向异性现象。在此基础上,采用本模型对中等超固结(OCR=4)高岭土试样三轴不排水剪切试验进行模拟,并再次与ABAQUS自带的修正剑桥模型模拟结果进行比较,表明本模型能够较好地反映中等超固结土在小应变情况下的非线性特性。相比于经典弹塑性模型,如修正剑桥模型,本模型的模拟结果更符合中等超固结土的变形特性。 展开更多

关键词 土体各向异性边界面模型 图形返回算法 ABAQUS软件 UMAT接口

在线阅读 下载PDF 职称材料

各向异性弹黏塑性模型在ABAQUS中的研发 被引量：1

14

作者 钦亚洲 李宁 许建聪 《岩土力学》 EI CAS CSCD 北大核心 2012年第4期1240-1246,共7页

通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载... 通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明:(1)对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2)对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。 展开更多

关键词 各向异性 弹黏塑性模型 回映算法 ABAQUS软件 UMAT子程序

在线阅读 下载PDF 职称材料

关于返回映射算法中应力的四阶张量值函数 被引量：6

15

作者 陈明祥 《力学学报》 EI CSCD 北大核心 2010年第2期228-237,共10页

针对各向同性材料,基于一组相互正交的基张量,建立了一套有效的相关运算方法.基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴.根据张量函数表示定理,... 针对各向同性材料,基于一组相互正交的基张量,建立了一套有效的相关运算方法.基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴.根据张量函数表示定理,本构方程和返回映射算法中所涉及到的应力的二阶、四阶张量值函数及其逆都由这组基所表示.推演结果表明:这些张量之间的运算,表现为对应系数矩阵之间的简单关系.其中,四阶张量求逆归结为对应的3×3系数矩阵求逆,它对二阶张量的变换则表现为该矩阵对3×1列阵的变换.最后,对这些变换关系应用于返回映射算法的迭代格式进行了相关讨论. 展开更多

关键词 本构方程积分 应力更新 返回映射算法 各向同性 张量函数表示定理

在线阅读 下载PDF 职称材料

金属成形数值模拟算法及软件开发 被引量：1

16

作者 周洋靖 冯志强 +1 位作者 宁坡 汪昌勇 《西南交通大学学报》 EI CSCD 北大核心 2019年第4期855-862,共8页

以实现拥有部分算法特色的金属成形定制化CAE(computer-aided engineering)软件为目的,在自主开发的OmtDesk平台基础上,嵌入金属成形数值计算程序.首先,以双势理论为基础,运用Uzawa算法求解金属成形中的接触摩擦问题;其次,采用建立在更... 以实现拥有部分算法特色的金属成形定制化CAE(computer-aided engineering)软件为目的,在自主开发的OmtDesk平台基础上,嵌入金属成形数值计算程序.首先,以双势理论为基础,运用Uzawa算法求解金属成形中的接触摩擦问题;其次,采用建立在更新拉格朗日框架上的返回映射算法,求解金属成形中材料非线性问题;最后,通过压缩回弹和挤压成形两个数值算例,验证了自主开发的金属成形定制化CAE软件FEM/Form的准确性与通用性,并将求解结果与ANSYS对比.在压缩回弹算例中,采用ANSYS得到10–3 mm量级的穿透量,相同算例采用FEM/Form则能将穿透量控制在10–7 mm量级,说明采用双势接触模块能更加精确地控制金属成形数值模拟中的穿透量,使金属成形过程更加准确. 展开更多

关键词 金属成形 接触摩擦 双势算法 弹塑性 返回映射算法

在线阅读 下载PDF 职称材料

黏弹-Perzyna黏塑性有限元法应力更新隐式算法

17

作者 闫富有 崔昊 +1 位作者 张晓婉 刘忠玉 《计算力学学报》 EI CAS CSCD 北大核心 2018年第5期611-618,共8页

黏弹-黏塑性耦合模型的黏弹性部分由弹簧、黏壶和Kelvin链串联而成,黏塑性部分为双曲线型DruckerPrager屈服函数、各向同性硬化和Perzyna黏塑性流动模型。基于黏弹性蠕变柔度,通过定义与弹性问题相对应的与时间增量相关的黏弹性剪切模... 黏弹-黏塑性耦合模型的黏弹性部分由弹簧、黏壶和Kelvin链串联而成,黏塑性部分为双曲线型DruckerPrager屈服函数、各向同性硬化和Perzyna黏塑性流动模型。基于黏弹性蠕变柔度,通过定义与弹性问题相对应的与时间增量相关的黏弹性剪切模量和体积模量,导出增量递推形式的本构方程。为保证算法的收敛和稳定性,把Perzyna黏塑性流动方程转化为与弹塑性相似的一致性条件,建立黏塑性增量因子单侧逼近其收敛值的N-R迭代算法。最后,给出应力更新完全隐式算法和最终计算公式。分别采用黏弹性、黏弹-塑性和黏弹-黏塑性本构关系对一地基蠕变模型进行三维有限元分析和比较,结果表明,本文算法具有较高的计算效率和稳定性。 展开更多

关键词 黏弹性 Perzyna黏塑性 应力更新算法 一致切线算子 有限元法

在线阅读 下载PDF 职称材料

基于修正剑桥模型的光滑粒子流体动力学的土柱崩塌模拟

18

作者 周学君 陈丁 唐轶 《科学技术与工程》 北大核心 2017年第31期1-7,共7页

土柱崩塌是一种被广泛研究的岩土大变形模型。考虑非关联流动法则的修正剑桥(MCC)模型描述土体行为;并将其导入光滑粒子流体动力学(SPH)算法来模拟土柱崩塌。在算法设计中,利用有限元方法中常用的中心回映算法来计算应力积分,以提高算... 土柱崩塌是一种被广泛研究的岩土大变形模型。考虑非关联流动法则的修正剑桥(MCC)模型描述土体行为;并将其导入光滑粒子流体动力学(SPH)算法来模拟土柱崩塌。在算法设计中,利用有限元方法中常用的中心回映算法来计算应力积分,以提高算法的精度。应用所设计的SPH算法模拟初始高宽比分别为1.5、3和5.5的土柱崩塌模型;并与实验结果比较,表明SPH算法在研究土柱崩塌问题的可靠性。 展开更多

关键词 土柱崩塌 光滑粒子流体动力学(SPH) 修正剑桥(MCC)模型 非关联流动法则 中心回映算法

在线阅读 下载PDF 职称材料

一种适用于饱和黏土循环动力分析边界面塑性模型的隐式积分算法 被引量：1

19

作者 刘艳秋 胡存 刘海笑 《岩土力学》 EI CAS CSCD 北大核心 2013年第12期3617-3624,共8页

修正了传统隐式回映算法,建立了适用于饱和黏土循环动力分析的边界面塑性模型的完全隐式积分格式。该模型基于无弹性域概念和临界状态理论,采用各向同性、运动硬化准则、旋转的边界面,并引入表征土体结构损伤和重塑程度的损伤变量以反... 修正了传统隐式回映算法,建立了适用于饱和黏土循环动力分析的边界面塑性模型的完全隐式积分格式。该模型基于无弹性域概念和临界状态理论,采用各向同性、运动硬化准则、旋转的边界面,并引入表征土体结构损伤和重塑程度的损伤变量以反映循环载荷作用下饱和黏土的各向异性、刚度、强度软化及塑性变形累积等特征。针对等压固结0K?1.0和偏压固结0K?0.67的饱和高岭黏土的不排水三轴试验进行模拟,采用不同的应变增量步长进行计算,并与试验数据对比,结果表明,修正隐式回映算法应用于该类边界面模型的合理性、积分格式的精确性和稳定性;另外,结合有限元软件自动时间步长的增量迭代解法,对饱和黏土应力控制的不排水动三轴试验进行预测,结果表明,修正的适用于该边界面的塑性模型隐式回映算法可以得到比较合理的数值分析结果,能够反映饱和黏土的循环刚度的退化和强度的弱化等动力特性。 展开更多

关键词 边界面塑性模型 隐式回映算法 循环动力特性 饱和黏土

在线阅读 下载PDF 职称材料

有限元实现含应变非强化区的双面循环塑性模型

20

作者 卿海 刘小明 +1 位作者 杨卫 吕坚 《应用力学学报》 EI CAS CSCD 北大核心 2008年第4期693-697,741,共5页

在有限单元方法的框架下,由广义径向返回法来实现含有应变非强化区的循双面环塑性本构模型。在小变形的假设下,本文详细描述了有限元方法实现率无关的双面循环塑性本构模型的完整算法。采用向后Euler隐式算法进行离散积分,其最终可化简... 在有限单元方法的框架下,由广义径向返回法来实现含有应变非强化区的循双面环塑性本构模型。在小变形的假设下,本文详细描述了有限元方法实现率无关的双面循环塑性本构模型的完整算法。采用向后Euler隐式算法进行离散积分,其最终可化简为一个非线性方程,并且由Peg-asus方法来求解。对径向返回法进行线性化,可以得到切向刚度算子。最后,文中通过两个算例表征了离散积分的精确性和材料模型的正确性。 展开更多

关键词 数值积分 径向返回法 应变非强化区 循环塑性模型 切向刚度算子

在线阅读 下载PDF 职称材料