[发明专利]一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法

权利要求书

1.一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法，其特征在于，包括以下步骤：

S1、建立连续纤维增强陶瓷基复合材料的宏观尺度构件有限元模型；

建立宏观尺度模型，修改其边界从而预设一条裂纹；

S2、建立包含任意复杂编织预制体的纤维束尺度有限元模型；

S21、建立预制体参数化模型；

结合增强纤维和纤维束截面形状及编织特点，开展了复杂预制体参数化精确建模，实现一维、二维、2.5维、三维增韧相参数化建模，同时完成不同形状纤维束及预制体设计；

S22、预制体水平集函数描述和自动离散；

结合水平集函数φ⁽ⁱ⁾(x,y,z)的值和预制体参数化模型，建立三维四面体节点处的水平集函数值和距离值；若单元同时存在函数f_R(x,y,z)0和函数f_R(x,y,z)0的节点，则该单元穿过增韧相边界Г，被定义为相交单元；则通过单元自动离散方法对该相交单元进行细分和编号及存储更新；

S3、通过导入CT扫描得到的孔隙信息，建立其形状、位置、尺寸、分布的统计模型，并在纤维束尺度模型引入自然缺陷特征；

S4、建立纤维丝尺度有限元模型，并通过统计模型引入自然缺陷特征；

S41、通过基于水平集函数的自动建模方法，建立纤维丝尺度模型；

S42、通过统计模型引入纤维丝尺度的孔隙特征，通过基于水平集函数的自动建模建立孔隙模型并将孔隙模型内的单元材料模量设为0；

S5、引入并行计算，计算宏观尺度下所有单元的高斯点应变，作为纤维束尺度下模型的边界条件；

S6、引入并行计算，计算纤维束尺度下所有单元的高斯点应变，以作为纤维丝尺度下模型的边界条件；

S7、采用CPU+GPU并行计算算法，在纤维丝尺度下计算平均应力和平均应变，返回纤维束尺度计算其纤维束尺度位移；

S8、在纤维束尺度下计算平均应力和平均应变，返回宏观尺度计算其宏观尺度位移；

S9、在宏观尺度下，引入上一步计算处的位移场，利用相场法求解相场并更新宏观尺度的边界条件，迭代整个过程，从而预测损伤演化。

2.根据权利要求1所述的一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法，其特征在于，所述步骤S3具体实现方法为：

S31、对连续纤维增韧陶瓷基复合材料进行断层CT扫描，智能识别连续纤维增韧陶瓷基复合材料内的所有自然孔隙，获得自然孔隙的图像及CT扫描孔隙数据；

S32、提取CT扫描孔隙数据，并产生有限元模型，引入如下的孔隙量化描述函数：

式中：x、y、z为孔隙笛卡尔坐标系中的坐标；R_x、R_y、R_z分别表示为x、y、z方向上的半径；

同时根据如下的孔隙随机分布参数D_v来描述其分布特征和E_v来量化孔隙逼近效果；

式中：V_p为孔隙的实际体积；V_f为拟合得到的实际体积；V_c为孔隙的实际区域和拟合得到的实际区域重合的体积；

S33、通过基于水平集函数的自动建模方法，建立孔隙结构；

S34、筛选出所有自然孔隙区域包含的所有四面体单元并记录；

S35、将所有自然孔隙区域的四面体单元材料中的模量设为0。

3.根据权利要求1所述的一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法，其特征在于，所述步骤S5具体实现方法为：

S51、通过以下公式计算宏观尺度下所有单元的高斯点的应变；

ε＝B·q

式中：q为节点位移矩阵；B为应变矩阵；ε为应变；

S52、将宏观尺度下计算得到的高斯点应变，作为纤维束尺度的边界条件，其对应公式为；

式中：表示在x轴正方向且与x轴垂直的有限元模型外表面上节点x⁺所对应的位移矢量；表示在x轴负方向且与x轴垂直的有限元模型外表面上节点x^-所对应的位移矢量；

S53、将宏观尺度模型的所有高斯点分成N份，采用并行命令对N份高斯点同时进行处理，设置下层纤维束尺度的边界条件。