一种二维编织陶瓷基复合材料氧化后剩余刚度预测方法
【专利摘要】本发明涉及一种二维编织陶瓷基复合材料氧化后剩余刚度预测方法。刚度指材料在受力时抵抗弹性变形的能力。通过分析材料各组分的刚度，即可确定材料内部的应力应变分布。因此本发明提供了一种能准确预测二维编织陶瓷基复合材料氧化后剩余刚度的方法。提出了考虑纤维氧化的动力学模型，在此基础上建立了考虑纤维氧化的微观尺度模型和二维编织陶瓷基复合材料中尺度模型。采用有限元法，通过施加周期性边界条件，计算了材料氧化后的剩余刚度。本发明能够精确的预测出该材料在不同氧化时间、不同氧化温度区间的剩余刚度，不需要通过实验花费大量的人力、物力去测试，因此节约了大量的试验成本。
【专利说明】
一种二维编织陶瓷基复合材料氧化后剩余刚度预测方法
技术领域
[0001] 本发明涉及一种二维编织陶瓷基复合材料氧化后剩余刚度预测方法。
【背景技术】
[0002] 二维编织陶瓷基复合材料具有高比强、高比模、耐高温、耐腐蚀和低密度等优良性 能，在空天飞行器高温防护系统具有广泛的需求。材料在其使用过程中，由于受到高温环境 因素的影响，会逐渐产生氧化损伤，导致材料力学性能下降，进而严重影响工程构件的使用 寿命与安全。刚度指材料在受力时抵抗弹性变形的能力，研究二维平纹编织陶瓷基复合材 料氧化后剩余刚度对其应用有着重要意义。
[0003] 由于单向陶瓷基复合材料存在非纤维方向力学性能弱等缺点，其应用范围受到了 限制。二维编织结构陶瓷基复合材料的出现，克服了单向复合材料的缺点，同时在厚度方向 上纤维束整体化更高，增加了材料层间剪切强度，减少了分层现象，并提高了复合材料抗冲 击性能和弯曲疲劳性能，因此大大扩展了陶瓷基复合材料的应用范围。
[0004] 然而由于二维编织陶瓷基复合材料是一种新型结构材料，国内外还没有高效的方 法预测其氧化后的剩余刚度，也未见公开的发明专利。杨成鹏（杨成鹏，矫桂琼，王波，等 .2D-C/SiC复合材料的氧化损伤及刚度模型[J].复合材料学报，2009，26(3): 175-181.)采 用实验的方法测试了2D C/SiC复合材料在700°C环境下的剩余刚度，并基于细观结构的变 化建立了计算公式，计算值与实验值较为吻合。但通过实验的方式要消耗大量的实验资金， 其提出的计算模型也只能计算离散的特定温度下的剩余刚度。
[0005] 当前，如何准确的预测二维编织陶瓷基复合材料氧化后的剩余刚度是本技术领域 重要而难以解决的问题。
【发明内容】
[0006] 发明目的：针对上述现有技术，提出一种能有效预测二维编织陶瓷基复合材料氧 化后剩余刚度预测方法。
[0007] 技术方案:一种二维编织陶瓷基复合材料氧化后剩余刚度预测方法，包括如下步 骤：
[0008] (1)，基于质量损失率理论和纤维退化规律假设，建立氧化动力学模型；
[0009] (2)，基于氧化动力学模型，采用有限元软件，建立氧化后的微观尺度的单胞模型；
[0010] (3)，施加微观尺度的单胞模型的周期性边界条件；
[0011] (4)，计算该单胞模型6个方向的弹性参数；
[0012] (5)，采用有限元软件，建立二维编织陶瓷基复合材料单胞模型；
[0013] (6)，将计算得到的氧化后微观尺度的单胞模型6个方向的弹性参数作为纱线的基 本属性，带入二维编织陶瓷基复合材料单胞模型；
[0014] (7)，施加所述二维编织陶瓷基复合材料单胞模型的周期性边界条件；
[0015] (8)，计算得到二维编织陶瓷基复合材料轴向的剩余弹性模量。
[0016] 作为本发明的优选方案，所述步骤（1)中，所述质量损失率理论分为两个温度区 间：
[0017] 1)当温度在400°C~700°C区间时，公式如下：
[0019] 其中，Ar是复合材料的质量损失率，W是复合材料的质量，A W是材料的质量变化 量，Ko是与氧化速率相关的常数，是氧气的体积分数，P是大气气压，M。是碳纤维的摩尔质 量，R是气体常数，T是环境温度，E r是氧化反应活化能，t是氧化时间，Srff是碳的有效反应面 积;其中，Srff=yW，y是碳的反应有效系数；
[0020] 2)当温度在700°C~900°C区间时，公式如下：
[0022]其中，Nc是碳的摩尔密度，A是与初始状态有关的常数，Tc是基体开裂温度，L c是涂 层厚度；
[0023]所述纤维退化规律假设:假设纤维以圆形的规律在高温下进行退化，公式如下：
[0025] 其中，S是纤维的氧化长度，p#Pp。分别表示纤维和复合材料的密度，L是复合材料 的长度，H是复合材料的高度，n是碳的物质的量，N f是单位面积内纤维的数量；
[0026] 将所述复合材料的质量损失率Ar带入所述纤维的氧化长度S的计算公式（3)，得 到：
[0027] 1)当温度在400°C~700°C区间时：
[0029] 2)当温度在700°C~900°C区间时：
[0031 ]由此可得氧化后纤维的剩余半径为：
[0032] 1)当温度在400°C~700°C区间时，公式如下：
[0034] 2)当温度在700°C~900°C区间时，公式如下：
[0036]其中，rf是纤维氧化后的剩余半径，rfo是纤维未氧化时的初始半径。
[0037] 作为本发明的优选方案，所述步骤(2)中，假设氧化为均匀贯穿式氧化，氧化后纤 维各处半径相等。
[0038] 作为本发明的优选方案，所述步骤(3)中，所施加的边界条件满足位移的连续性和 在模型相反两个平面应力分布的一致性。
[0039] 作为本发明的优选方案，所述步骤(4)中，所述6个方向的弹性参数包括x、y、z三个 方向的弹性模量￡\上上，"、^、72三个方向的剪切模量6#^ 2為2，以及泊松比化、^、 Vyz 〇
[0040] 作为本发明的优选方案，所述步骤(7)中，施加周期性边界条件为：
[0041 ] uf+ = 3iXjZt +u] (8)
[0042] +t( (9)
[0043] 其中，Z+和Z-分别表示垂直于Z轴的两个相反的边界表面，》f+为在Z+边界表面上 的位移，为在Z-边界表面上的位移， Xlz+为Z+表面上节点的位移量，x产为Z-表面上节点 的位移量，<为在边界表面上位移的周期性部分，h为周期性结构的平均应变张量。
[0044]有益效果:本发明提供的二维编织陶瓷基复合材料氧化后剩余刚度预测方法，基 于质量损失率模型和纤维退化规律假设，提出了考虑纤维氧化的动力学模型。基于氧化动 力学模型，采用有限元法建立了考虑纤维氧化的微观尺度模型和二维编织陶瓷基复合材料 单胞尺度模型，预测了材料的剩余刚度。本发明提出的预测模型充分考虑了纤维随氧化时 间、温度的退化规律，因此能够精确的预测出二维编织陶瓷基复合材料氧化后的剩余刚度 且节约了大量的实验成本。
【附图说明】
[0045] 图1是裂纹尖端附近碳纤维氧化扫描电镜图；
[0046] 图2是陶瓷基复合材料二维平面模型；
[0047] 图3是单向陶瓷基复合材料氧化示意图；
[0048]图4是氧化后的微观模型；
[0049]图5是在氧化后微观模型边界条件示意图；
[0050]图6是二维平纹编织陶瓷基复合材料单胞模型；
[0051 ]图7是二维平纹编织陶瓷基复合材料单胞模型边界条件示意图；
[0052]图8是该预测模型的具体流程图；
[0053]图9是空气环境下700°C时二维平纹编织C/SiC复合材料轴向的剩余弹性模量预测 值与的实验值的对比曲线；
[0054]图10是空气环境下800°C时二维平纹编织C/SiC复合材料轴向的剩余弹性模量预 测值与的实验值的对比曲线；
[0055] 图11是空气环境下850°C时二维平纹编织C/SiC复合材料轴向的剩余弹性模量预 测值与的实验值的对比曲线；
[0056] 图12是空气环境下900°C时二维平纹编织C/SiC复合材料轴向的剩余弹性模量预 测值与的实验值的对比曲线。
【具体实施方式】
[0057]下面结合附图对本发明做更进一步的解释。
[0058] 本实施例以二维平纹编织C/SiC复合材料为例，对700 °C~900°C区间范围内材料 氧化后的剩余刚度进行预测，其中材料性能参数如表1所示。
[0059]表 1
[0061 ]如图8所示，本方法具体步骤如下：
[0062] (1 )，基于质量损失率理论和纤维退化规律假设，建立氧化动力学模型。其中，质量 损失率理论分为两个温度区间：
[0063] 1)当温度在400°C~700°C区间时，复合材料的质量损失率公式如下：
[0065]其中，Ar是复合材料的质量损失率;W是复合材料的质量；A W是复合材料的质量损 失率的质量变化量；Ko是与氧化速率相关的常数；是氧气的体积分数，在这里取 20.95% ;P是大气气压，在这里取101.325KPa;MC是碳纤维的摩尔质量，在这里取12 X 103kg/ mol ;R是气体常数，在这里取8.3145J/(mol ? K) ;T是环境温度;Er是氧化反应活化能；t是氧 化时间；Srff是碳的有效反应面积，Srff与试样的质量有关，表示为Srff = iiW，ii是碳的反应有 效系数，y取决于试样的微裂纹面积和孔隙横截面积以及试样密度，可通过实验测得。
[0066] 2)当温度在700°C~900°C区间时，复合材料的质量损失率公式如下：
[0068] 其中，N。是碳的摩尔密度;A是与初始状态有关的常数;T。是基体开裂温度，在这里 取1030°C;U是涂层厚度。
[0069] 纤维退化规律假设:根据图1所示的扫描电镜照片和图2所示的陶瓷基复合材料二 维平面模型，假设纤维以圆形的规律在高温下进行退化，如图3所示，纤维的氧化长度即为 图中BD段的长度;ro即0D的距离，指的是拟合的氧化区域的圆心0到纤维未氧化时表面的距 离;r'即为0C的距离，指的是拟合的氧化区域的半径;a即为ro和r'之间的夹角。纤维退化规 律公式如下：
[0071]其中，S是纤维的氧化长度;Pf和P。分别表示纤维和复合材料的密度;L是复合材料 的长度，H是复合材料的高度，如图2所示;n是碳的物质的量;Nf是单位面积内纤维的数量。
[0072] 将复合材料的质量损失率Ar带入纤维的氧化长度S的计算公式(3)，得到：
[0073] 1)当温度在400°C~700°C区间时：
[0075] 2)当温度在700°C~900°C区间时：
[0077]若纤维未氧化时的初始半径为rfQ，纤维氧化后的剩余半径为rf，根据关系rf = rf〇-S可得氧化后纤维的剩余半径为：
[0078] 1)当温度在400°C~700°C区间时，公式如下：
[0080] 2)当温度在700°C~900°C区间时，公式如下：
[0082] (2)，基于上述步骤（1)中的氧化动力学模型即式(6)和（7)，采用ANSYS软件，建立 氧化后的微观尺度的单胞模型，如图4所示。其中假设氧化为均匀贯穿式氧化，氧化后纤维 各处半径相等。
[0083] (3)，施加微观尺度的单胞模型的周期性边界条件:所施加的边界条件满足位移的 连续性和在微观尺度的单胞模型相反两个平面应力分布的一致性，如图5所示;所施加的周 期性边界条件如表2所示。
[0084] 表 2
[0086]表中：S(x'y，z)、S(x+，y，z)分别为x方向坐标最小和最大的两个平面，S(x，y'z)、 S(x，y+，z)分别为y方向坐标最小和最大的两个平面，S(x，y，z1、S(x，y，z+)分别为z方向坐 标最小和最大的两个平面。u x、uy、uz分别为x、y、z方向的位移约束，const表示平面内所有节 点位移耦合。
[0087] (4)，计算该单胞模型6个方向的弹性参数，6个方向的弹性参数包括x、y、z三个方 向的弹性模量Ex、Ey、Ez，xy、XZ、yz三个方向的剪切模量GXy、GXZ、Gy Z，以及泊松比VXy、VXZ、VyZ。
[0088] (5)，采用ANSYS软件，建立二维编织陶瓷基复合材料单胞模型，如图6所示;模型的 尺寸参数如表3所示。
[0089] 表3
[0091]其中，a为纱线宽度，b为同向纱线间间距，hf为纱线厚度，hb为基体层厚度。
[0092 ] (6 )，将计算得到的氧化后微观尺度的单胞模型6个方向的弹性参数作为纱线的基 本属性，带入二维编织陶瓷基复合材料单胞模型；
[0093] (7)，施加二维编织陶瓷基复合材料单胞模型的周期性边界条件，如图7所示，其表 达式为：
[0094] ^ = 3txtz+ +u* ：C8)
[0095] (9)
[0096]其中，Z+和Z-分别表示垂直于Z轴的两个相反的边界表面，w.f+为在Z+边界表面上 的位移为在Z-边界表面上的位移，Xlz+为Z+表面上节点的位移量，x产为Z-表面上节点 的位移量，〃)为在边界表面上位移的周期性部分，I为周期性结构的平均应变张量。
[0097] Z+和Z-两个边界面之间的相对位移量表示为：
[0098] wf+ - -.=我(xf+ -
[0099] 其中，Ajf表示Z+和Z-两个边界面上节点的位移改变量。
[0100] (8)，根据公式
，EL为材料的轴向弹性模量，^为单元轴向平均应力，^为 单元轴向平均应变，即可计算得到二维平纹编织C/SiC复合材料轴向的剩余弹性模量。
[0101] 图9给出了 700°C环境温度下二维平纹编织C/SiC复合材料轴向的剩余弹性模量预 测值与的实验值的对比曲线。图10给出了 800°C环境温度下二维平纹编织C/SiC复合材料轴 向的剩余弹性模量预测值与的实验值的对比曲线。图11给出了 850°C环境温度下二维平纹 编织C/SiC复合材料轴向的剩余弹性模量预测值与的实验值的对比曲线。图12给出了 900°C 环境温度下二维平纹编织C/SiC复合材料轴向的剩余弹性模量预测值与的实验值的对比曲 线。通过对比可见，本发明的方法能够有效预测二维编织陶瓷基复合材料氧化后的剩余刚 度。
[0102] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种二维编织陶瓷基复合材料氧化后剩余刚度预测方法，其特征在于，包括如下步 骤： (1 )，基于质量损失率理论和纤维退化规律假设，建立氧化动力学模型； (2)，基于氧化动力学模型，采用有限元软件，建立氧化后的微观尺度的单胞模型； (3 )，施加微观尺度的单胞模型的周期性边界条件； (4) ，计算该单胞模型6个方向的弹性参数； (5) ，采用有限元软件，建立二维编织陶瓷基复合材料单胞模型； (6) ，将计算得到的氧化后微观尺度的单胞模型6个方向的弹性参数作为纱线的基本属 性，带入二维编织陶瓷基复合材料单胞模型； (7) ，施加所述二维编织陶瓷基复合材料单胞模型的周期性边界条件； (8) ，计算得到二维编织陶瓷基复合材料轴向的剩余弹性模量。2. 根据权利要求1所述的二维编织陶瓷基复合材料氧化后剩余刚度预测方法，其特征 在于:所述步骤(1)中，所述质量损失率理论分为两个温度区间： 1) 当温度在400 °C~700 °C区间时，公式如下：其中，是复合材料的质量损失率，W是复合材料的质量，AW是材料的质量变化量，Ko是 与氧化速率相关的常数，是氧气的体积分数，Ρ是大气气压，Μ。是碳纤维的摩尔质量，R是 气体常数，Τ是环境温度，Er是氧化反应活化能，t是氧化时间，Srff是碳的有效反应面积;其 中，Srff=yW，y是碳的反应有效系数； 2) 当温度在700 °C~900 °C区间时，公式如下：其中，N。是碳的摩尔密度，λ是与初始状态有关的常数，T。是基体开裂温度，L。是涂层厚 度； 所述纤维退化规律假设:假设纤维以圆形的规律在高温下进行退化，公式如下：其中，δ是纤维的氧化长度，ρ#Ρρ。分别表示纤维和复合材料的密度，L是复合材料的长 度，Η是复合材料的高度，η是碳的物质的量，Nf是单位面积内纤维的数量； 将所述复合材料的质量损失率k带入所述纤维的氧化长度δ的计算公式(3)，得到： 1) 当温度在400 °C~700 °C区间时：2) 当温度在700 °C~900 °C区间时：由此可得氧化后纤维的剩余半径为： 1) 当温度在400 °C~700 °C区间时，公式如下：2) 当温度在700 °C~900 °C区间时，公式如下：其中，rf是纤维氧化后的剩余半径，rfo是纤维未氧化时的初始半径。3. 根据权利要求1所述的二维编织陶瓷基复合材料氧化后剩余刚度预测方法，其特征 在于:所述步骤(2)中，假设氧化为均匀贯穿式氧化，氧化后纤维各处半径相等。4. 根据权利要求1所述的二维编织陶瓷基复合材料氧化后剩余刚度预测方法，其特征 在于:所述步骤(3)中，所施加的边界条件满足位移的连续性和在模型相反两个平面应力分 布的一致性。5. 根据权利要求1所述的二维编织陶瓷基复合材料氧化后剩余刚度预测方法，其特征 在于:所述步骤(4)中，所述6个方向的弹性参数包括x、y、z三个方向的弹性模量E x、Ey、Ez， xy、xz、yz三个方向的剪切模量Gxy、Gxz、Gyz，以及泊松比v xy、vxz、Vyz。6. 根据权利要求1所述的二维编织陶瓷基复合材料氧化后剩余刚度预测方法，其特征 在于:所述步骤(7)中，施加周期性边界条件为：其中，Z+和Z-分别表示垂直于Z轴的两个相反的边界表面，为在Z+边界表面上的位 移，为在Z-边界表面上的位移，Xlz+为Z+表面上节点的位移量，X产为Z-表面上节点的位 移量，< 为在边界表面上位移的周期性部分，專为周期性结构的平均应变张量。
【文档编号】G06F17/50GK105930579SQ201610244709
【公开日】2016年9月7日
【申请日】2016年4月19日
【发明人】孙志刚, 陈西辉, 宋迎东
【申请人】南京航空航天大学