drcker-prave广义中点法在各向同性硬化条件下的应用.docxVIP

drcker-prave广义中点法在各向同性硬化条件下的应用

ruker和rough在1952年提出了标准。在最简单的数学形式上，它类似于莫霍-科茨的标准，用来设计土力学的弹性极限。这是著名的ruker-roll的标准。从其提出至今近60年的时间,Drucker-Prager弹塑性本构关系在岩土力学与工程中取得了极为广泛的应用与深远的发展。本构关系数值积分是指微分型本构关系在有限时间增量上的积分,是非线性数值分析的核心内容。基于DP弹塑性模型优越的数学形式与广泛的应用,值得对其数值积分进行独立研究。

总应变包括弹性应变和塑性应变:

弹性应变与应力满足广义Hooke定律:

其中,C为各向同性弹性柔度张量。

存在塑性势函数g及流动法则:

各向同性硬化条件:

其中:f(σ,uf06b)为屈服函数;H为硬化模量。

Kuhn-Tucker条件:

的求解方法主要分为两类:

第一类方法是基于常应变率假设,推导本构方程式(1)的积分精确解。Krieg等推导了vonMises理想弹塑性材料本构积分精确解;Carosio等将Krieg的方法应用于理想J2粘塑性材料;我国学者王勗成和常亮明推导了运动硬化vonMises材料本构关系的精确积分,并将其推广应用到各向同性硬化与混合硬化;Loret和Prevost以及RezaieePajand和Nasirai采用不同方法推导了DP理想弹塑性材料本构积分精确解。与vonMises材料精确解不同,DP材料精确解并非完全解析解,而是需要数值求解复杂的一元超越方程。本构方程积分精确解为数值积分方法的验证和精度评价提供了依据,但是对于更复杂的本构关系,精确解几乎无法得到。

第二类是数值积分方法,其核心是利用积分中值定理,将微分型本构方程转化为增量型非线性方程,通过数值方法进行求解。Ortiz和Popov系统研究了弹塑性数值积分算法的精度和稳定性,将常见的方法归纳为两类:广义梯形法(GeneralizedTrapezoidalMethod,GTM)和广义中点法。通过控制参数的不同取值,这两类方法涵盖了常用的显式和隐式积分方法。Borja等给出三种三主不变量型弹塑性模型(Mohr-Coulomb,Lade-Duncan和Matsuoka-Nakai)的向后欧拉法(CPPM别称);Hjiaj等给出非关联流动DP材料的应力调整方法,其中考虑了屈服面尖点的奇异;Tamagnini等研究了一般各向同性硬化材料的向后欧拉法,并用于模拟带状岩体的风化。

杨强等[11―12]提出关联流动理想弹塑性DP材料基于线性预测-径向校正策略(ElasticPredictorRadialCorrectorScheme)的显式解析解,并证明该解在积分策略上完全等价于隐式的CPPM,因此在效率,精度和稳定性三方面具有优越性。张小寒进一步证明,这种解法可由CPPM推广到更一般的GMM,同时基于精确解,分析加载步长和加载偏离角度对计算精度的影响。

岩土工程中的极限分析对应于采用关联流动法则的理想弹塑性分析,可采用文献建议的解析解。然而考虑非关联流动以及各向同性硬化的DP弹塑性材料同样非常重要。例如Rudnichi和Rice提出一种压力相关弹塑性本构关系描述散粒材料和脆性岩石的变形局部化,被后来学者大量引用。该本构关系推广到一般三维情况,就是非关联流动各向同性硬化DP弹塑性材料。本文采用GMM对这种本构关系进行数值积分,给出调整后最终应力的解析解,并对其精度进行评价。

1本构关系积分的方法选择

图1所示为一典型的弹塑性加载过程。域内一点初始应力状态满足f(σ0,k0)≤0;对该点施以应变增量,对应于弹性试应力:

其中,D=C-1为各向同性弹性刚度张量。如果导致塑性加载。弹性加载和塑性加载的分界点一般按线性条件确定:

后文仅讨论塑性加载,并将σc视为已知。

弹塑性本构关系积分对应于一个满足流动法则和硬化条件的应力调整过程。由式(1a)和式(1b)可知,经过调整的最终应力:

满足一致性条件f(σ1,k1)=0。

利用积分中值定理,流动法则式(1c)可写为如下增量形式::

其中,。精确满足式(5)的塑性势梯度mm存在但无法确定,需要基于一定的假设选取,这也是不同积分方法的主要差异。

本文采用GMM,其塑性势梯度mm取:

当且仅当控制参数时,GMM退化为显式的向前欧拉法,其余情况均为隐式方法。隐式方法的优点是自动满足一致性条件,并且具有良好的计算精度与数值稳定性;缺点是对于一般本构关系需要迭代求解,计算效率较低,并且迭代过程可能不收敛。特别地,为文献中讨论的CPPM,其本身具有一阶精度。文献证明,只要,GMM就是无条件稳定的。

2多态性弹塑性cppm的应力调整策略

本文沿用Rudnichi和Rice的描述,DP屈服准则为:

如果假定硬化模量