单晶叶片应力分析方法及应用.pdf

单晶叶片应力分析方法及应用 杨治国成晓鸣 尹泽勇 株〔洲航空动力机械研究所) 摘〔要〕采用晶体滑移理论，给出了单晶弹塑性姗变滑移本构模型。建立了全新的有限元列 式。并用于某气冷单晶祸轮叶片弹塑料性姗变应力分析。 关键词CA休滑移理论弹塑性姗变滑移本构模型有限元列勾 1.前言 由于单晶合金无晶界，多晶合金中常沮时的晶界强化元素在高温时成为拓变、疲劳损 伤滋的这一缺点，在单晶合金中得以遥免。因此，国外得以广泛应用，目前国内除小发外 也已用于某大发气冷涡轮叶片.然而，单晶合金与以往常用的多晶合金有很人不同。它是 一种各向异性材料，其物理尤其是机械性能其有方向性、不对称性等与各向同性材料迥然 不同的特点。另外，在高沮环境下，不随时间变化的塑性与随时间变化的场变同时出现. 更增加了问硬的复杂性。可以说。以往的分析方法难于适用单晶结构分析，有必要改进或 建立新的弹塑性蟠变本构模型和相应分析方法。 国内外自开发各类单晶合金以来.就十分注重开展对单晶合金机械性能及单晶结构 分析方法的研究。但有关单晶涡轮叶片结构分析未见报导。为保证采用这种先进材料的单 晶叶片研制工作的顺利进行，必须开展单晶叶片的分析研究。 2 基于滑移理论的弹塑性滑移本构模型 扮变形描述 图1所示，F为总变形梯度，F`为晶格畴变和刚体转动所产十的V,形梯度.F0140it 滑移产生的变形梯度，且F=FF。分别设m(a)、nla’为a滑移系变形前的滑移方向和滑 移面法向方向的单位矢绪，则晶格畸变后，新的滑移方向和滑移面的法线方向m n.’”‘:ma(l=Fma()·na()=[F{厂‘」na() 当前构形的速度梯度张量V V=V`+VP V`=F`(F`)一， VP=FFP(FP)一’(Fe)-=全m’沈“，。n’。，‘。T宁。，‘ 式中F为变形梯度率，,),(a)为。滑移系的剪应变率。 记变形率张量D、旋转率张量译为 D=D+DP W=W亡+WP 根据公式(2).有: D。一含{户F(·Y+[F( DP一Ep-(Y-) pa(，一IFm.()n(a()JT+na()际a()丫) ZL 、 ’ 、 ’」 一 88一 W0=2fF0Fy一!F‘(0r1]一’F(.Y}Wr二全W(a)Ya() w(。)=工 2 m[a()n((),7-n.(a)m((.)1TI 2)率无关弹塑性本构关系 设晶体的弹性性质不受滑移变形的影晌，弹性本构关系如下: i.二L:D` (4) 式中L为弹性模tI张tI,扩为以中间构形为羞础状态的Kirchhoff应力的Jaumann导数.且 奄.二t一W.T+，W. (5) 式中t为Kirchhoff应力的物质导数。相对初始构形状态的Kirchhoff应力Jaumann导数为 i二L:D一J(LP(a)+p‘“，)宁.‘， (6) p(以)=w()“卜rwa() (7) 在有限变形情况下，滑移系a上的分切应力，.‘)及其物质导数为 兮(。)二，:P(a) 8 ) ( t(a)mX(a卜D0