对循环塑性模型的广义径向返回算法.doc

对循环塑性模型的广义径向返回算法 卿海1，刘小明1，杨卫1,2 1清华大学，工程力学系，北京100084 2 浙江大学，杭州 310027 摘要 本文描述了由广义径向返回法对含应变非强化区的率无关循环塑性模型进行有限元执行的完整算法。文中采用向后Euler隐式积分算法进行离散积分将多变量问题最终化简为一个非线性方程，并且推导得到了切向刚度算子。文中通过一个算例验证了广义径向返回法的精确性和稳定性。 关键词： 径向返回法，应变非强化区，循环塑性模型，切向刚度算子 引言 精确描述循环载荷下结构的力学行为需要一个合适的材料模型，这方面的一些成果可参见文献[1-4]。Ohno等[1]和Chaboche[2]在假设循环强化取决于塑性应变幅值的基础上提出应变非强化区，在应变非强化区内不发生各向同性强化。他们的模型可以表征材料的循环强化和弱化，饱和迟滞响应及不可恢复的记忆特征等力学特征。 大量研究表明，应用径向返回算法的向后Euler隐式积分可以高效地对材料模型进行离散积分。该算法可以将率形式的多变量方程组最终化为一个标量形式的非线性方程[5-9]。本文叙述有限元执行Ohno模型[1]，重点放在发展广义径向返回法及推导相应的切向刚度算子。文中假设材料处于小变形状态，同时采用Einstein 求和约定和固定的笛卡尔直角坐标系统。 含非强化区的双面循环塑性本构 按照应变率的加法分解，总应变率可以表示成弹性与塑性部分两部分之和 (1) 应力与弹性应变满足广义Hooke定律如下 (2) 其中， (3) 这里,和是Lame常数。下面概述Ohno循环塑性模型。 多轴应力状态下的限制面与屈服面(见图1)，其表达式分别为： (4) 其中，和分别为屈服面和限制面上的偏应力状态，和指两个面的中心(背应力)，而和分别指两个面的大小。两个面的背应力及限制面的半径的演化规律如下 (5) 这里，是与应变非强化区相关的量，将在后面介绍；指累积塑性应变增量，是与相关的一个标量， , (6) 图1 多轴应力状态下的限制面与屈服面示意图 图2 多轴应力状态下的非强化区 多轴应力状态下的应变非强化区(见示意图2)可以表示为 (7) 这里是背应变，为非强化区的大小(初值为)，其演化规律如下 , (8) 其中， ，且 ，其余情况 (9) 因此，应用有限元来实现Ohno模型时，一共需要确定11个材料参数，分别为，，，，，，，，，和。 隐式积分过程 完全隐式积分(向后欧拉法)在大规模计算 (包括大量时间步的计算) 中表现出很好的收敛性和精确度。因此，这里采用向后欧拉法来对Ohno塑性本构方程进行离散积分。 为了简化下面的表达式，本文中用上标(n)表示第n步的量，而无上标(n)的量为第n+1步的更新量。例如，和分别表示第n和第n+1步的应力值。由式(2)可得， (10) 这里，为偏应力空间里的弹性试应力，该应力是假设增量步内不发生塑性变形的情况下得到的。 首先，计算。若，该增量步内没有塑性变形，则有及，而其它量保持不变。反之，，该增量步有塑性变形发生 (11) 从式(5)可得 (12) 因而， (13) 将式(12)代入(13)，经过化简得 (14) 其中， (15) 定义在偏量空间内的Von-Mises距离函数 (16) 对式(14)两边分别双点积，可以得到 (17) 同时， (18) 由式(17)中第一等式所定义的Eq值是数值迭代中屈服状态的表征。在方程(17)中，有和两个变量。在求解时，可再次采用径向返回法来确定当前应变是否处于应变非强化区内。首先令，求解方程(17)，得到，此时的塑性应变称为塑性试应变。还需要检验，如果，则处于应变非强化区，成立，物理量作相应的更新。否则，，塑性应变等量还必须满足应变一致性条件， , (19) 通过推导，可以得到 (20) 由式(18)知，由唯一确定。因此，式(20)中也由唯一确定。求解由式(20)和式(17)组成的非线性方程组时，可以先将(20)代入(17)求解，再将代入(20)来求解。然后对所有的物理量作相应的更新。整个积分过程的流程见图3。 图3 隐式积分流程图 这里采用“PEGASUS”[10]方法来求解方程(17)。该方法不需要一阶导数，能够无条件稳定地收敛到问题的解，且收敛速度远高于割线法和二分法。应用“PEGASUS”方法需要确定上下界。可选时的解作为一个下界，当处于屈服时显然有。上界则必须使得，可以通过以下步骤得到。假定全加载步均导致塑性变形，那么应力卸除将引起弹性应变的减小。减小后的应力存在上界 (21)