一种陶瓷基复合材料的疲劳寿命的多尺度预测方法与流程

本发明涉及一种陶瓷基复合材料的疲劳寿命的多尺度预测方法，用以解决航空发动机热端部件使用陶瓷基复合材料时的疲劳寿命预测，属于航空发动机领域。

背景技术：

与传统金属材料相比，陶瓷基复合材料耐温能力提升近20％，密度下降约2/3，在1300～1400℃下无需冷却可保持较高的强度，覆盖热障涂层后无需冷却可在1600℃下长时间使用。陶瓷基复合材料应用于航空发动机热端部件，可直接减轻重量60％左右，还可以降低冷却气体用量15～25％，大大降低耗油率和co2排放量。因此，陶瓷基复合材料成为先进军民用航空发动机热端部件的理想材料。

疲劳失效是航空发动机构件最主要的失效模式之一。陶瓷基复合材料应用于航空发动机，首先需要解决的是其在疲劳载荷下的强度及寿命预测问题。

目前传统的疲劳寿命预测方法及相关模型主要用于金属材料，对于复杂的编织陶瓷基复合材料，并没有发展出成熟的疲劳寿命预测方法。

技术实现要素：

本发明提出了一种陶瓷基复合材料的疲劳寿命的多尺度预测方法，用以解决目前针对航空发动机热端部件使用陶瓷基复合材料时的疲劳寿命预测问题。

本发明采用的技术方案为：一种陶瓷基复合材料的疲劳寿命的多尺度预测方法，主要分为四个步骤，分别如下：

步骤一：1.1输入计算需要的编织陶瓷基复合材料的纤维和基体材料参数。

步骤一：1.2建立纤维和基体的疲劳失效准则，其表达式如下：

基体疲劳损伤累积速率：

其中dm表示基体损伤累积量，n表示循环数，am、bm和cm为材料常数。δσeq,m表示基体的等效应力。σult,m表示基体的静态极限应力。

当用于拉拉疲劳时，δσeq,m＝δσeq,mt；当用于压压疲劳时，δσeq,m＝δσeq,mc；当用于拉压疲劳时，δσeq,m＝δσeq,mt+δσeq,mc。

当p<0时为拉伸状态，当p>0时为压缩状态。δσeq,mt表示基体等效拉伸应力，δσeq,mc表示基体等效压缩应力。其中和分别表示基体的拉伸强度和压缩强度，p和σvon分别表示静水压力和vonmises等效应力。

纤维疲劳损伤累积速率：

其中di表示纤维束损伤累积量，n表示循环数，ai、bi和ci为材料常数。δfi表示静态损伤量。

δf1取值定义如下：

当用于纤维束轴向拉拉疲劳时，δf1＝δf1t；当用于纤维束轴向压压疲劳时，δf1＝δf1c；当用于纤维束轴向拉压疲劳时，δf1＝δf1t+δf1c。

δf2＝δf3取值定义如下：

当用于纤维束横向拉拉疲劳时，当用于纤维束横向压压疲劳时，当用于纤维束横向拉压疲劳时，

步骤二：基于扫描电镜的陶瓷基复合材料细观几何特征参数获取，对于编织陶瓷基复合材料，通过扫描电镜获取编织角，体积分数，纤维截面形状等参数。

步骤三：根据获得的几何特征参数信息建立陶瓷基复合材料周期结构rve。

本步骤主要基于步骤二中获取的编织陶瓷基复合材料的细观几何特征参数建立包括二维和三维的编织陶瓷基复合材料的rve模型。

对于二维编织陶瓷基复合材料，获得纤维束中心线曲线方程，其中x和y分别为局部坐标。

横截面曲线方程其中λ用来表征纤维束截面的扁平程度，a、b分别为纤维束宽度和厚度的一半。

体积分数vf可表示为

对于三维编织陶瓷基复合材料，中心线为直线；横截面曲线方程，其中c，d分别表示椭圆长短轴。

步骤四：周期性跳跃算法完成对陶瓷基复合材料疲劳寿命跨越式计算。具体包括以下两步：

一方面：依据疲劳载荷情况，使用雨流法提取出疲劳载荷谱。

另一方面：利用获得的疲劳载荷谱，使用周期跳跃方法进行寿命计算。

本发明使用的方法获得的编织陶瓷基复合材料的疲劳寿命相较于使用传统金属材料疲劳方法预测编织陶瓷基复合材料寿命精度提高30％；除此之外，使用周期跳跃算法同时兼顾计算精度和计算量，通过较少的增量步模拟得到经历了较高循环数后的结构疲劳渐进损伤扩展情况，使得计算效率大为提高。

附图说明

图1为周期跳跃算法框图；

图2为本发明预测的编织陶瓷基复合材料寿命与传统方法的对比；

图3为本发明总体流程图。

具体实施方式

下面结合附图以及具体实施方式方式进一步说明本发明。

本发明涉及一种陶瓷基复合材料的疲劳寿命的多尺度预测方法，针对编织陶瓷基复合材料疲劳寿命的预测，具体实施方式如下：

主要分为四个步骤，分别如下：

步骤一：1.1输入计算需要的编织陶瓷基复合材料的纤维和基体材料参数。

步骤一：1.2建立纤维和基体的疲劳失效准则，其表达式如下：

基体疲劳损伤累积速率：

其中dm表示基体损伤累积量，n表示循环数，am、bm和cm为材料常数。δσeq,m表示基体的等效应力。σult,m表示基体的静态极限应力。

当用于拉拉疲劳时，δσeq,m＝δσeq,mt；当用于压压疲劳时，δσeq,m＝δσeq,mc；当用于拉压疲劳时，δσeq,m＝δσeq,mt+δσeq,mc。

当p<0时为拉伸状态，当p>0时为压缩状态。δσeq,mt表示基体等效拉伸应力，δσeq,mc表示基体等效压缩应力。其中和分别表示基体的拉伸强度和压缩强度，p和σvon分别表示静水压力和vonmises等效应力。

纤维疲劳损伤累积速率：

其中di表示纤维束损伤累积量，n表示循环数，ai、bi和ci为材料常数。

δf1取值定义如下：

当用于纤维束轴向拉拉疲劳时，δf1＝δf1t；当用于纤维束轴向压压疲劳时，δf1＝δf1c；当用于纤维束轴向拉压疲劳时，δf1＝δf1t+δf1c。

δf2＝δf3取值定义如下：

当用于纤维束横向拉拉疲劳时，当用于纤维束横向压压疲劳时，当用于纤维束横向拉压疲劳时，

步骤二：基于扫描电镜的陶瓷基复合材料细观几何特征参数获取，对于编织陶瓷基复合材料，通过扫描电镜获取编织角，体积分数，纤维截面形状等参数。

步骤三：根据获得的几何特征参数信息建立陶瓷基复合材料周期结构rve。

本步骤主要基于步骤二中获取的编织陶瓷基复合材料的细观几何特征参数建立包括二维和三维的编织陶瓷基复合材料的rve模型。

对于二维编织陶瓷基复合材料，获得纤维束中心线曲线方程，其中x和y分别为局部坐标。

横截面曲线方程其中λ用来表征纤维束截面的扁平程度，a、b分别为纤维束宽度和厚度的一半。

体积分数vf可表示为

对于三维编织陶瓷基复合材料，中心线为直线；横截面曲线方程，其中c，d分别表示椭圆长短轴。

步骤四：周期性跳跃算法完成对陶瓷基复合材料疲劳寿命跨越式计算。具体包括以下两步：

一方面：依据疲劳载荷情况，使用雨流法提取出疲劳载荷谱。

另一方面：利用获得的疲劳载荷谱，使用周期跳跃方法进行寿命计算。

周期跳跃方法可以分为以下四步：

步骤1、在abaqus中使用多尺度渐进损伤模型，执行当前循环的强耦合计算，进行损伤扩展与宏细观应力重分布；

步骤2、根据第1步的计算结果，采用步骤一定义的纤维和基体的疲劳损伤累积模型评估损伤速率；

步骤3、在疲劳计算中，进行njump(i)个周期跳跃，更新基体与纤维束损伤状态；

步骤4、判断循环数或结构状态是否已经到达设定限值，如果是，则停止计算；否则回到第1步。

通常来说，数值模拟时会遇到起始循环时损伤速率较大的情况。例如对于位移控制的循环加载。因此，起始的周期跳跃数njump(i)应较小以保证计算精度。同时，计算过程中njump(i)逐渐增大，以保证计算效率。作为参考，逐次跳跃的周期循环数可取为njump(i)＝1,1,2,4,8,…,32768,65536,10⁵,10⁵,…。