正交异性统一粘塑性本构模型-中国力学学会.docVIP

正交异性统一粘塑性本构模型 及其在航空发动机热端部件上的应用 陈火红 ，何君毅 中国航空工业六二八所 摘 要：本文提出了一个能描述定向结晶、单晶高温合金材料特性的正交异性统一粘塑性本构模型，它克服了与时间无关的热弹塑加蠕变的传统分析方法中一些根本的缺点，更精确地反映了在变化中的高温及机械载荷作用下的结构材料的性态。通过大量的材料高温非常规试验拟合了某镍基合金该模型方程的材料参数。此粘塑性统一本构模型为一个非线性常微分组，在数学上被认为较为“刚硬”，难以直接用差分求解，为此本文针对所提出的本构模型方程并借鉴国外相关研究的进展发展了一种具有较高效率且稳定的积分方法。而后将所提出的本构模型及积分方法编成程序，通过大型通用非线型程序MARC的用户子程序HYPELA与MARC集成。最后利用MARC程序及由试验得到的材料参数，对一个实际航空发动机涡轮叶片进行了粘塑性分析，说明所提出的本构模型、积分方法及程序可以满足当前定向结晶、单晶燃气热端部件结构分析的要求。 关键词：正交异性；本构模型；粘塑性；涡轮叶片；燃气轮机热端部件 1.前 言 航空发动机热端部件，如涡轮叶片与涡轮盘、燃烧室、涡轮导气叶片等工作在严劣的环境下，承受严重的热、机械循环载荷，会产生较大的非弹性应变，工作寿命很短。目前，通常采用定向结晶或单晶高温合金材料来满足它们的安全性和耐久性的要性。但要对它们进行强度校核和寿命估算，必须先分析得到它们关键部位的准确的应力、应变历程。 传统分析高温与时间相关的结构采用热弹塑性加蠕变的方法，该法从理论上存在一些根本性缺点，如： 1.不能充分反映高温时材料性质与时间严重相关的事实，应变率对材料性质有较大影响。 2.不能反映周期载荷作用下的循环硬化（或软化）现象等。 3. 传统分析理论中蠕变只与初始应力水平有关，但实际上它还与蠕变前的应力历程相关。 统一粘塑性本构理论将非弹性应变统一成一个与时间相关非弹性应变，克服了人为将非弹性应变分为与时间无关的弹塑性应变及与时间相关的蠕变应变的牵强假定，以应变率的形式可考虑时间相关，因此与实际更接近，可克服传统理论的不足，故可较大地提高理论的精确度。 如果以参考资料[1]作为起点，对粘塑统一本构理论的研究已逾三十年，但适合涡轮叶片定向结晶或单晶正交异性材料的本构模型研究较少。根据对已有的研究成果，认为将各向同性Walker模型[2]推广到各向异性构成修正Walker模型较好，它对材料性质的力学描述比较清晰，对各种高温力学现象，如：循环硬（软）化、松驰、蠕变均能较好表达，而且Walker模型在各向同性情况下较成熟。 定向结晶材料与单晶材料从连续力学角度均可视为正交异性，可用本文所述方法分析。但近年来研究单晶材料多从微观滑移理论出发研究粘塑模型，这是由于单晶晶粒方向是已知的，易于由微观推至宏观。而对于定向结晶材料横向晶粒增长是随机的，只能采用连续力学方法宏观上讨论它的本构模型。2.正交异性Walker模型 本文约定：对定向结晶而言材料坐标轴3为结晶（纵向）方向，1、2轴为垂直于结晶方向的横向；对单晶而言三个方向均为结晶方向。对正交异性材料有如下应力与应变关系 {σ}=[C]（{ε}-{εI}-{α}Δθ） （1） 式中 {σ}T=[σ11 σ22 σ33 σ23 σ31 σ12]为真实应力。 由于材料正交异性质，三个轴线热胀系数为 {α}T=[α1 α2 α3 0 0 0] [C]为材料弹性系数矩阵， 由于定向结晶横向各向同性: {ε}为材料总应变， {εI}为非弹性应变， Δθ为当前温度与初始温度之差。 在粘塑性本构方程中一般可将方程分为两部分，即描述应力与非弹性应变率（为计及与时间相关一般采用率形式而不是应变本身）关系的流动方程与描述材料硬化（软化）与恢复的内部状态变量本身变化的演进方程 修正Walker模型的流动方程 式中 Mijkl为材料各向异性塑性张量。 aij为背应力（back stress）张量分量，它描述材料随动硬化及Bauschinger效应。 Sij为偏应力张量分量。 K为阻应力（drag stress）标量，它描述了材料初始及继屈服强化效应。 Walker模型之演进方程为 式中：n1、n2、n3、n6、n7、a0、k1、k1、n、m为方程常系数由试验确定。εI，a，K等变量用以描述材料硬化与恢复性质的变量，通常是无法直接测量得到，而采用微分方程表达它们的变化过程，这类变量通称内部状态变量。 本文将各向异性弹塑性Hill理论[5]推广至粘塑性。可以将二阶张量σij、εIij、aij写矢量，则方程（3）、（5）、（6）中四阶张量Mijkl写