有限应变下蠕变2弹塑性2损伤耦合本构模型的数值模拟Ξ.pdfVIP

第 14 卷第 2 期 计 算 力 学 学 报 V o l. 14 N o. 2 1997 年 5 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS M ay 1997 有限应变下蠕变弹塑性损伤耦合本 构模型的数值模拟 李锡夔 (大连理工大学工程力学研究所, 大连, 116023) 摘　要　本文提出了一个模拟有限应变下各向同性蠕变弹塑性损伤耦合本构行为的数值模型。模 型基于不可逆过程的热动力学和内状态变量理论, 以及应力应变本构关系的超弹性形式描述和变形 梯度的弹性及非弹性部分的乘式分解。导出了计算应力和内状态变量的一致性积分算法和考虑有限 应变下蠕变弹塑性损伤完全耦合效应的一致性切线矩阵。 关键词　有限应变; 蠕变; 弹塑性; 损伤; 耦合数值模型; 一致性算法 分类号　 242. 1 O 1　引　言 80 年代后期以来, 有限应变下材料非线性本构关系数值模拟及其算法的研究已受到了广 泛的关注。其主要特点是应用应力应变本构关系的超弹性描述和变形梯度的弹塑性部分乘式 〔1〕 分解 。相对于过去的利用应力应变关系亚弹性描述以及应变张量和式分解的途径, 这个途径 〔2 〕 〔3 〕 〔4 〕 的优点已有广泛的讨论 。在这个途径下, 、 和 分别提出了有限应变下 Sim o W eber A nand 〔5 〕 弹塑性及粘弹塑性本构数值模型。 和 采用主对数应变作为应变度量建立了有限 Sim o T aylo r 〔6 〕 弹性的本构数值模型, 近来, Sim o 和M eschke 把它推广到弹塑性本构模型中。 基于这个新途径, 本文发展了一个有限应变下蠕变弹塑性损伤耦合本构数值模型, 其目 的在于试图模拟高填充高分子一类材料的本构行为。这个本构模型的主要特点可概括为: ( 1) 材料在拉伸和压缩方向呈现不同的损伤效应, 定义两个损伤内变量以计及杨氏弹性模量的降 低和泊松比随杨氏模量降低所引起的变化; (2) 材料在拉伸和压缩方向具有不同的塑性屈服强 ( ) 度, 虽然材料的初始屈服强度很低, 但在应变软化 呈现材料不稳定 前具有很高的应变硬化; (3) 应力应变关系的滞后现象。 2　各向同性, 蠕变弹塑性损伤模型 选择物体的参考构形 R 3。构形 中材料点定义为 ∈ 。 时刻构形中的点表示为 x t x = ( , ) 。对于一个固定时刻 , 由 ( , ) 所描述的变形梯度 、速度梯度 和变形率 可分 x t tn x tn F l d 别定义为: x - 1 1 T F = ; 　l = F F ; 　d = ( l + l ) ( 1)