共查询到17条相似文献，搜索用时 411 毫秒
1.
2.
统一弹塑性本构模型在ABAQUS中的开发与应用 总被引：1，自引：1，他引：0
基于统一弹塑性本构模型的有限元理论格式,根据ABAQUS的UMAT格式要求,编制相应的接口程序,将统一弹塑性本构模型引入ABAQUS中。采用退化的统一强度模型（ 时,为Mohr-Coulomb模型）与ABAQUS自带的Mohr-Coulomb模型,对单轴试验和圆形硐室进行弹塑性分析,验证所开发材料子程序的正确性及高效性。考虑到统一强度模型的一般情况（ ）和屈服面硬化条件,对圆形硐室进行弹塑性计算,得到应力场的变化规律。所提出的研究思路具有普遍性,为采用ABAQUS平台进行本构模型的二次开发提供了借鉴和参考。 相似文献
3.
将SFG模型中的加载湿陷屈服面（LC屈服面）与BBM模型中提出的吸力增加屈服面（SI屈服面）统一成为单独光滑的屈服面,基于弹塑性理论框架,推导出相应的弹塑性本构模型。通过非饱和固结直剪试验对江西红土的模型计算参数的应力-应变曲线与试验结果进行对比。结果表明：模型对多种应力路径下非饱和三轴试验的计算结果与相应试验的测定结果较为吻合,验证了模型能较好地模拟非饱和土在多种应力路径下的体应变行为。统一屈服面简化了加载湿陷屈服面与吸力增加屈服面之间的耦合分析,解决了因双屈服面的不连续而造成计算不便问题,拓宽了SFG模型的适用性。 相似文献
4.
5.
Zienkiewicz-Pande屈服准则是Mohr-Coulomb准则的改进,在p-q子午面和π平面上都是光滑曲线,不存在尖点,不仅在数值迭代计算过程中易于处理,而且在一定程度上考虑了屈服曲线与静水压力的关系以及中主应力σ2对屈服的影响。该准则已应用于多个水电站地下洞室开挖支护计算中,但并未被大型数值计算软件所采用。根据弹塑性力学原理,详细推导了基于Zienkiewicz-Pande屈服准则的弹塑性本构模型在FLAC3D中的增量迭代计算格式;在FLAC3D提供的二次开发程序接口基础上,利用VC++程序语言实现了计算过程,并编译成动态链接库文件进行加载和调用。将自编的本构模型对一理想圆形隧洞进行模拟,计算结果与Mohr-Coulomb准则一致,验证了计算格式的正确性。将此本构模型应用于一大型地下洞室群的开挖围岩稳定分析中,得到的塑性区、位移场、应力场分布与Mohr-Coulomb准则计算结果分布规律相符,进一步验证了模型的适用性和可靠性。 相似文献
6.
高地应力下硬岩的本构模型研究 总被引：1，自引：0，他引：1
传统的本构模型在模拟高地应力下硬岩破坏的范围和深度方面并不理想。针对高应力下的拉西瓦花岗岩,通过真三轴压缩试验模拟了开挖时应力路径的演化。在摩尔-库仑强度准则、格里菲斯强度准则、德鲁克-普拉格准则、双剪统一强度理论和霍克-布朗强度准则的基础上,考虑十二面单元体主剪面上主剪应力、正应力和静水压力的共同作用,提出了一个有关高应力下硬岩的三剪强度准则。根据试验结果,利用遗传算法全局寻优的功能,搜索出了三剪强度准则的参数。预测样本和与摩尔-库仑强度准则、德鲁克-普拉格准则、双剪统一强度理论的比较表明,该三剪强度准则与试验结果吻合,适用于高地应力下的硬岩。根据试验得到的应力-应变关系曲线,建立了基于三剪强度准则和应变软化的弹-脆-塑性本构模型,并用FLAC3D提供的基于C++的用户自定义模型工具UDM,创建了用户自定义模型的动态链接库,嵌入了FLAC3D软件。采用遗传算法-FLAC方法,搜索得到了基于岩石试样和岩体的本构模型参数。计算结果与实测情况吻合,表明建立的本构模型适用于高地应力下的硬岩,为高地应力下硬岩地下工程的安全性和稳定性分析奠定了基础。 相似文献
7.
热黏弹塑性本构模型是描述土在温度（热）和时间（黏）耦合作用下的应力-应变关系的本构模型。在一些新型岩土工程诸如高放核废料地质处置、地热资源开发与贮存的建设中,需要同时考虑温度和时间对土的影响,所以建立一个热黏弹塑性本构模型具有理论和实际意义。将温度变化对黏土体积和强度参数的影响引入笔者之前提出的超固结土等向应力-应变-时间关系,建立了一个等向应力条件下的应力-应变-时间-温度关系。随后,基于该关系推导了屈服面硬化定律,并将其与超固结土统一硬化模型的屈服方程和流动法则结合,建立了超固结土的热黏弹塑性本构模型。最后,使用新模型预测室内试验,证明新模型能够反映时间和温度对土体积、一维压缩曲线和前期固结压力的耦合影响。 相似文献
8.
《岩土力学》2021,(2)
本文将SFG模型中的加载湿陷屈服面(LC屈服面)与BBM模型中提出的吸力增加屈服面(SI屈服面)统一成为单独光滑的屈服面,基于弹塑性理论框架,推导出相应的弹塑性本构模型。通过非饱和固结直剪试验对江西红土的模型计算参数进行测定。将计算参数输入模型程序并通过程序对已发表文献中非饱和三轴试验的应力路径进行计算,并将模型计算得到的应力-应变曲线与试验结果对比后得出结论:模型对多种应力路径下非饱和三轴试验的计算结果与相应试验的测定结果较为吻合,验证了模型能较好的模拟非饱和土在多种应力路径下的体应变行为。统一屈服面简化了加载湿陷屈服面与吸力增加屈服面之间的耦合分析,解决了因双屈服面的不连续而造成计算不便问题,拓宽了SFG模型的适用性。 相似文献
9.
考虑变形模量劣化的应变软化模型在FLAC3D中的开发与验证 总被引：2，自引：0，他引：2
FLAC3D软件中的应变软化本构模型是通过抗剪强度参数劣化来实现的,没有考虑变形模量劣化造成的影响。针对这一问题,采用了变形模量劣化的损伤本构模型,详细推导了劣化损伤本构模型的差分格式,研究该模型在FLAC3D中二次开发的流程,在VC++环境下实现了劣化损伤本构模型在FLAC3D中的二次开发;通过与深埋圆形洞室应变软化数值解算例对比,进一步验证了该模型弹塑性和损伤的力学性质,分析了变形模量劣化对围岩变形特性的影响。由于该模型能够反映围岩的劣化损伤特性和大变形特点,因此在工程应用中有一定的实用价值。 相似文献
10.
在季冻区高等级公路路基场地取土,制成重塑土试件进行冻融循环试验。而后对经历不同冻融循环次数的土体进行三轴固结排水（CD）剪切试验,测得应力、应变数据。在大量试验数据的基础上,遵循岩土本构变化规律,以岩土弹塑性力学为理论基础,在p-q空间,以椭圆方程拟合体变屈服面,以抛物线方程拟合剪切屈服面,采用塑性关联流动,得到考虑冻融循环条件下融土的弹塑性本构方程。通过编程计算后,对比计算与试验结果,证明该双屈服面本构模型能够较为准确地预估冻融循环土体的应力-应变关系,该本构模型对季冻区路基土体长期稳定性分析和工程预测具有重要参考意义。 相似文献
11.
复杂应力状态下土质高边坡稳定性分析 总被引：1，自引：1，他引：0
结合现有的实际滑坡场地,开展了关中盆地黄土塬边的高边坡稳定性研究。基于弹塑性理论和非关联流动法则,根据双剪统一强度理论,建立了双剪统一弹塑性本构模型的显式有限差分格式。针对陕西省泾阳县南塬高边坡滑坡的现场勘测和室内试验,分析了滑坡的规模和滑坡前后的地形特征。在此基础上,还原了滑坡前的边坡地形概状,建立三维计算模型并运用三维拉格朗日有限差分方法分析了中间主应力对高边坡稳定性的影响;然后考虑中间主应力的影响对不同地下水埋深和渗流时边坡的稳定性进行分析。结果表明,考虑中间主应力情况下,随着地下水位升高,坡顶水平向位移增大;当不考虑渗流时,坡顶水平位移增大到一定值后边坡达到稳定状态,而考虑渗流时坡顶水平位移持续增大,并且容易发生突然的破坏失稳。期望对泾阳南塬滑坡群的破坏机制研究提供一定参考。 相似文献
12.
针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager（D-P）强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法（NR-AL法）联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明：所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。 相似文献
13.
将解析反演和数值反演相结合,对试验毛洞的断面形状和所处应力场做出适当的简化,利用对应性原理求解出洞室位移的Kelvin粘弹性解,并以此去拟合试验洞不同位置处的实测位移而得到多组流变参数;分别将其作为FLAC3D数值模型蠕变模式的输入参数,得出各自对应的实际断面形状和实际应力场条件下试验洞的蠕变位移数值解。然后,利用BP神经网络,对各组流变参数和其对应的试验洞同一位置处的蠕变位移数值解进行训练,建立起两者之间的非线性关系;利用训练好的网络,依据实测蠕变位移值得出了围岩流变参数;最后,利用隧道实测位移数据对反演的流变参数进行了数值正分析验证。 相似文献
14.
堆石料三维边界面模型在FLAC3D中的开发与验证 总被引：1，自引：0，他引：1
堆石料三维边界面模型结合了统一非线性强度准则,可以反映三维应力空间堆石料的应变硬化、软化以及体积收缩和剪胀等传统本构模型难以反映的力学特性。基于三维边界面模型的基本理论,通过FLAC3D提供的二次开发平台,在VC++环境下实现了三维边界面本构模型的二次开发,并给出基于FLAC3D程序的边界面本构模型开发的关键步骤、编程要点和调试方法。利用开发的本构模型,开展了三轴压缩试验的数值模拟计算,并与理论值进行了对比分析。结果表明,嵌入在FLAC3D中的三维边界面本构模型能够较好地反映设定试验条件下的材料性能,而且三维边界面本构模型模型参数简单,都可从常规三轴试验获得,从而验证了三维边界面模型二次开发程序的优越性与合理性。 相似文献
15.
深海能源土是指含天然气水合物（俗称"可燃冰"）的深海沉积物,其本构特性的模拟对可燃冰的安全开采至关重要。首先分析了水合物对能源土强度、剪胀和软化等力学特性的影响机理,水合物饱和度越大,对能源土力学特性影响越显著。然后在修正剑桥模型的基础上,通过引入水合物的饱和度和统一硬化参数来修正屈服函数,以反映水合物对能源土强度、剪胀、软化等特性的影响,建立了能考虑天然气水合物胶结作用形成及退化影响的深海能源土弹塑性本构模型,推导了相应的弹塑性矩阵。最后,通过模拟结果与已有能源土三轴试验数据对比分析,表明模型能很好地预测能源土强度、剪胀和软化等特性,验证了模型的合理性和有效性。 相似文献
16.
滑坡形成是一个典型的岩土变形破坏时效过程,低速滑坡时效过程尤为显著。基于流变理论,建立反映滑坡变形破坏发展过程的本构模型、预测滑坡活动趋势,一直是国内外滑坡研究的基本途径和难点问题之一。然而,目前国内外已有的多数流变模型仅能反映滑坡岩土蠕变的第一、第二阶段,不能刻画滑坡岩土蠕变的第三阶段（加速蠕变阶段）。本文基于滑带在滑坡中的作用以及滑带岩土剪切蠕变发展的累进性和非线性特征,借助损伤理论,在FLAC3D内置的CVISC流变本构模型中引入非线性损伤黏塑性元件,构建了可描述滑坡加速蠕变过程的非线性损伤流变本构模型,依据类比原理建立了改进的CVISC三维差分模型,通过FLAC3D开放的用户接口实现了本构模型的二次开发,并将改进的CVISC模型用于长期缓慢滑移、伴随间歇性剧烈活动的甘肃舟曲泄流坡数值模拟中。模拟结果显示,该模型不仅呈现了滑坡的加速蠕变特征,而且揭示的滑坡活动特征与其曾经出现的活动特征基本一致,从而证实基于非线性损伤理论的改进CVISC模型具有较好的有效性。 相似文献
17.
底部滑脱层流变性和厚度变化对褶皱冲断带影响的数值模拟 总被引：1，自引：0，他引：1
滑脱层是影响褶皱冲断带的第一因素。为了探讨底部滑脱层流变性和厚度变化对褶皱冲断带的影响,本文通过二维有限差分软件FLAC开展了数值模拟研究。常见的泥岩滑脱层具摩擦性质,此类冲断带采用弹塑性模型;而盐岩滑脱层具蠕变性质,此类冲断带则采用黏弹塑性模型。结果表明,底部滑脱层流变性以及厚度对冲断带的变形影响很大。当滑脱层为泥岩时,发育典型的叠瓦状构造样式,变形序列为背驮式向前扩展的有序序列。而当滑脱层为盐岩时,发育背、向斜相间排列的侏罗山式构造样式,并以无序序列演化。底部滑脱层厚度的变化,对泥岩滑脱层的冲断带而言,主要改变其表面坡角,对构造样式和演化序列影响不大;但对盐岩滑脱层的冲断带影响很大,不管滑脱层最厚的地方分布在后陆还是前陆,变形都会最强烈,地表凸起最高。从后陆向前陆加厚的盐岩滑脱层结合由后陆向前陆略为减薄的上覆地层可产生罕见的后冲、后展式演化序列。这些变形的差异主要是由于盐岩比泥岩具有小得多的剪切强度所致。 相似文献