一种基于分段权函数的涡轮盘裂纹扩展模拟方法技术

本发明专利技术涉及一种基于分段权函数的涡轮盘裂纹扩展模拟方法，步骤包括：(1)基于涡轮盘静强度分析，确定危险部位并得到相应的应力分布；(2)基于应力强度因子模拟结果，确定任意非线性(包括均布、一次和二次)应力分布权函数；(3)基于步骤(2)建立的权函数，开展涡轮盘裂纹扩展模拟，得到涡轮盘裂纹扩展规律。

【技术实现步骤摘要】
一种基于分段权函数的涡轮盘裂纹扩展模拟方法
本专利技术是一种针对航空发动机涡轮盘的裂纹扩展模拟方法，它是一种考虑涡轮盘裂纹扩展平面应力分布情况基于分段权函数求解裂纹强度应力从而模拟涡轮盘裂纹扩展的方法，属于航空航天发动机

技术介绍
根据断裂力学可知，应力强度因子与结构形状、尺寸、载荷及裂纹长度密切相关，为保证裂纹扩展分析的计算效率，工程上通常采用简化模型代替实际结构，以远端均布应力代替实际载荷，从而建立应力强度因子的经验公式，如现有文献"NewmanJC,RajuIS.Stressintensityfactorequationsforcracksinthree-dimensionalfinitebodiessubjectedtotensionandbendingloads[R].Hampton,Virginia,USA:NASALangleyResearchCenter,1984."针对三维矩形板裂纹模型，建立了多种类型裂纹在承受拉伸和弯曲载荷下的应力强度因子近似表达式。对于复杂结构而言，简化模型端面应力往往呈现非均匀应力分布，导致经验公式中的均布载荷假设不成立，从而使裂纹扩展预测结果存在一定偏差。因此，如何在非均布应力场下准确求解应力强度因子，是决定裂纹扩展模拟精度、效率的关键一环。现有文献"BuecknerHF.Anovelprincipleforthecomputationofstressintensityfactors[J].JournalofAppliedMathematicsandMechanics,1970,50(10):529-546."提出权函数方法，用以计算复杂应力分布下的应力强度因子。现有文献"KiciakA,GlinkaG,EmanM,etal.Weightfunctionsandstressintensityfactorsforcornerquarter-ellipticalcrackinfinitethicknessplatesubjectedtoin-planeloading[J].EngineeringFractureMechanics,1998,60(2):221-238."建立了复杂非均匀应力分布下三维矩形板裂纹体的应力强度因子通用权函数计算方法。然而该方法只适用于应力沿板厚度方向分布的情况，而对于航空发动机涡轮盘等部位，其裂纹面应力分布主要是沿着径向分布的周向应力，而将此裂纹进行简化成三维矩形板裂纹时，裂纹面上应力分布则沿矩形板宽度方向分布，因此该方法不太适用。
技术实现思路
本专利技术技术解决方案：克服现有涡轮盘裂纹扩展模拟方法的不足，提供一种任意非线性应力分布裂纹应力强度因子的权函数表达式方法，能够充分考虑涡轮盘裂纹扩展平面应力分布情况采用分段权函数模拟涡轮盘偏心孔处裂纹扩展，提高涡轮盘裂纹扩展模拟预测精度。本专利技术技术解决方案：一种基于分段权函数的涡轮盘裂纹扩展模拟方法，结合涡轮盘静强度分析结果，确定危险部位及相应的应力分布，根据应力强度因子模拟结果，确定任意非线性应力分布权函数，最后结合涡轮盘危险部位平面应力分布采用分段拟合权函数开展涡轮盘裂纹扩展模拟，实现步骤如下：(1)基于涡轮盘静强度分析，确定危险部位。为提高有限元计算效率、减少计算量，取涡轮盘一个扇区区域作为分析模型，施加周期对称边界条件；为便于网格划分将涡轮盘切分为榫槽、偏心孔、盘心和盘身多个子模块，采用六面体单元划分网格，提高计算精度；查阅高温合金手册获得材料参数，包括弹性模量、泊松比和密度，覆盖涡轮盘全部服役温度；载荷为涡轮盘设计工作状态的转速和叶片离心力(转换为榫槽接触面拉应力)，温度场由涡轮盘典型位置温度通过拟合获得；位移边界条件为约束安装边的轴向和周向位移；进行静强度分析，确定涡轮盘危险部位。(2)基于应力强度因子模拟结果，确定任意非线性应力分布权函数。对于孔边椭圆角裂纹，如图2(a)所示，记孔外表面裂纹前缘为A点(对应裂纹长度a)、孔内表面裂纹前缘为B点(对应裂纹长度b)。根据权函数法定义，A、B两点的应力强度因子可表示为：其中，σ(x)为A、B两点的应力，mA(x,a)、mB(x,a)分别为A、B两点的权函数，可表示为：其中，a为裂纹长度，M1A、M2A、M3A、M1B、M2B、M3B为代表结构特征的待定系数，通过参考载荷状态下的应力强度因子有限元解加以确定。分别选取沿简化模型x向呈均布、线性、二次分布的裂纹面正应力作为参考载荷状态，如图4所示，用以确定上述待定系数。其中，σ0为常数，a为裂纹长度，为几何修正因子，分别对应i＝0,1,2分别表示参考载荷为均布、线性分布和二次分布，Q为第2类完全椭圆积分近似表达式的平方，Q＝E(k)2，时，时，将图4中的应力分布及式(3)、式(4)带入式(1)、式(2)，积分得：其中，KiX(i＝0,1,2,X＝A,B)为A、B两点分别在参考载荷为均布、线性分布和二次分布时的应力强度因子，FiX(i＝0,1,2,X＝A,B)为几何修正因子。将M1A、M2A、M3A、M1B、M2B、M3B表示为F0A、F1A、F2A、F0B、F1B、F2B的形式：影响应力强度因子的几何参数主要有：裂纹长度a和b、矩形板厚度t，采用有限元奇异单元法计算孔边椭圆角裂纹分别在裂纹面承受均布、线性分布和二次分布应力时的应力强度因子以建立应力强度因子表达式中的几何修正因子。基于参考应力强度因子解，采用多项式函数分布拟合获得了形式较为统一的、高精度的参考应力强度因子表达式。取FiA、FiB(i＝0,1,2)为(b/t)的三次多项式，即：其中，αij、βij(i＝0,1,2、j＝0,L,3)为几何修正因子FiX(i＝0,1,2,X＝A,B)的拟合系数，可表示为(b/a)的二次多项式，即：其中，为αij、βij的拟合系数。Franc3D是目前唯一可以计算一般各向异性材料中三种断裂模式应力强度因子K的断裂力学软件，且其理论模型适用性强、软件参数设置灵活、计算精度较高，因此本专利技术选取Franc3D结合有限元软件Ansys进行断裂力学计算。通过有限元建立有限厚度(t＝7.8mm)平板角裂纹模型，在σ0＝100MPa下分别计算图4所示均布、线性、二次分布载荷下A、B两点的应力强度因子、几何因子及参数αij、βij。采用最小二乘法拟合确定pkij取值。同理，对于孔边穿透裂纹，在权函数法中，应力强度因子、几何因子、拟合参数之间的关系可表示为：其中，K为应力强度因子，σ(x)为应力，m(x,a)为权函数，可表示为：其中，M1、M2、M3为代表结构特征的待定系数，通过参考载荷状态下的应力强度因子有限元解加以确定。其中，Fi(i＝0,1,2)为几何修正因子。将M1、M2、M3表示为F0、F1、F2的形式：取Fi(i＝0,1,2)为(b/t)的二次多项式，即：其中，αik(i＝0,1,2)为Fi的拟合系数。通过有限元建立有限厚度(t＝7.8mm)平板穿透裂纹模型，在σ0＝100MPa下分别计算图4所示均布、线性、二次分布载荷下裂纹尖端应力强度因子、几何因子及拟合参数αik取值。以上，计算得到孔边椭圆角裂纹分别在裂纹面承受均布、线性分布和二次分布应力时的应力强度因子解及经验关系式，将其分别带入式(本文档来自技高网...

【技术保护点】
1.一种基于分段权函数的涡轮盘裂纹扩展模拟方法，其特征在于，包括以下步骤：(1)基于涡轮盘静强度分析，确定涡轮盘危险部位；取涡轮盘一个扇区作为分析模型，施加周期对称边界条件；将涡轮盘切分为榫槽、偏心孔、盘心及盘身多个子模块，采用六面体单元划分网格；根据高温合金手册获得材料参数，所述材料参数包括弹性模量、泊松比和密度，覆盖涡轮盘全部服役温度；施加机械载荷和温度场，所述机械载荷包括转速和叶片离心力，所述温度场由涡轮盘典型位置拟合得到；进行静强度分析，确定涡轮盘危险部位；(2)基于裂纹扩展应力强度因子模拟结果，确定任意非线性应力分布权函数；根据权函数的定义，基于孔边椭圆角裂纹和穿透裂纹简化模型，分别选取沿简化模型x向呈均布、线性、二次分布的裂纹面正应力作为参考载荷状态，采用有限元奇异单元法通过参考载荷状态下的裂纹扩展应力强度因子确定权函数中的待定系数，最终获得孔边椭圆角裂纹和孔边穿透裂纹的权函数；(3)基于步骤(2)建立的权函数，开展涡轮盘裂纹扩展模拟；分析由步骤(1)确定的涡轮盘危险部位应力特点，确定涡轮盘裂纹扩展方向，绘制相应的裂纹扩展平面应力分布并拟合应力分布的表达式，由拟合参数结合步骤(2)建立的权函数计算裂纹扩展应力强度因子；初始1/4圆形角裂纹半径选定1～3mm，每次计算间隔给定裂纹扩展增量选定0.05～0.2mm，得到涡轮盘危险部位表面裂纹长度与内表面裂纹长度随循环数N的变化关系，从而模拟获得涡轮盘裂纹扩展规律。...

【技术特征摘要】
1.一种基于分段权函数的涡轮盘裂纹扩展模拟方法，其特征在于，包括以下步骤：(1)基于涡轮盘静强度分析，确定涡轮盘危险部位；取涡轮盘一个扇区作为分析模型，施加周期对称边界条件；将涡轮盘切分为榫槽、偏心孔、盘心及盘身多个子模块，采用六面体单元划分网格；根据高温合金手册获得材料参数，所述材料参数包括弹性模量、泊松比和密度，覆盖涡轮盘全部服役温度；施加机械载荷和温度场，所述机械载荷包括转速和叶片离心力，所述温度场由涡轮盘典型位置拟合得到；进行静强度分析，确定涡轮盘危险部位；(2)基于裂纹扩展应力强度因子模拟结果，确定任意非线性应力分布权函数；根据权函数的定义，基于孔边椭圆角裂纹和穿透裂纹简化模型，分别选取沿简化模型x向呈均布、线性、二次分布的裂纹面正应力作为参考载荷状态，采用有限元奇异单元法通过参考载荷状态下的裂纹扩展应力强度因子确定权函数中的待定系数，最终获得孔边椭圆角裂纹和孔边穿透裂纹的权函数；(3)基于步骤(2)建立的权函数，开展涡轮盘裂纹扩展模拟；分析由步骤(1)确定的涡轮盘危险部位应力特点，确定涡轮盘裂纹扩展方向，绘制相应的裂纹扩展平面应力分布并拟合应力分布的表达式，由拟合参数结合步骤(2)建立的权函数计算裂纹扩展应力强度因子；初始1/4圆形角裂纹半径选定1～3mm，每次计算间隔给定裂纹扩展增量选定0.05～0.2mm，得到涡轮盘危险部位表面裂纹长度与内表面裂纹长度随循环数N的变化关系，从而模拟获得涡轮盘裂纹扩展规律。2.根据权利要求1所述的一...

【专利技术属性】
技术研发人员：王荣桥，胡殿印，胡如意，崔金铎，毛建兴，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11