[发明专利]纤维增强复合材料等效弹性参数的计算方法

权利要求书

1.一种纤维增强复合材料等效弹性参数的计算方法，其特征在于：包括如下步骤

步骤1：将试件样品进行XCT扫描，获得逐层的灰度图像序列；

步骤2：假设试件断层扫描图片序列一共有M幅图片，每幅图片的像素尺寸为宽W像素、高H像素，每个像素的灰度范围是0-255，其中采用(j,i,k)（k=0，1，2～M-1；i=0,1,2～H-1；j=0,1,2～W-1）来表示第k+1幅图片，第i+1行，第j+1列像素，将试件断层扫描图片序列读入内存，建立像素的灰度数组color[]，数组中的元素color[W*H*k+W*i+j]表示像素(j,i,k)的灰度值；

步骤3：确定基体、碳纤维和孔隙的参考灰度值colorM、colorF和colorV，理论上，孔隙的灰度值为0，基体和碳纤维的灰度值分别开展XCT实验，获得纯基体和纯碳纤维的灰度值，且实验参数与试件扫描时的参数相同，纯基体和纯碳纤维试件采用反应烧结SiC板和碳纤维布，碳纤维布用树脂固定；

步骤4：创建单元阵列，其通过一个像素代表一个三维8节点正方体单元，按照步骤2的设定，创建的单元阵列的长、宽和高分别包含W、H和M个单元；

步骤5：确定单元的弹性参数，已知碳纤维和基体的弹性张量分别为和其中f和m分别表示纤维和基体，i,j,k,l=1,2,3，孔洞的弹性张量设定为一个接近零而不为零的数值；

采用线性插值函数来计算单元(j,i,k)的弹性张量，假设co=color[W*H*k+W*i+j]，在纤维束区域，如果colorV<co<=colorF，则单元(j,i,k)的弹性张量

Eijkl=(co-colorV)(colorF-colorV)Eijklf+(co-colorF)(colorV-colorF)Eijklv---(1)]]>

如果colorF<co<=colorM，则单元的弹性张量

Eijkl=(co-colorM)(colorF-colorM)Eijklf+(co-colorF)(colorM-colorF)Eijklm---(2)]]>

在基体区域，则单元的弹性张量可用下式计算：

Eijkl=(co-colorM)(colorV-colorM)Eijklv+(co-colorV)(colorM-colorV)Eijklm---(3)]]>

步骤6：施加边界条件，对步骤4和步骤5建立的单元阵列依次施加如下位移边界条件：

第1组：x＝0,u_x＝0;y＝0,u_y＝0;z＝0,u_z＝0;x＝w,u_x＝w*0.001  (4a)

第2组：x＝0,u_x＝0;y＝0,u_y＝0;z＝0,u_z＝0;y＝l,u_y＝l*0.001  (4b)

第3组：x＝0,u_x＝0;y＝0,u_y＝0;z＝0,u_z＝0;z＝h,u_z＝h*0.001  (4c)

第4组：x＝0,u_x＝u_y＝u_z＝0;x＝w,u_z＝w*0.001  (4d)

第5组：x＝0,u_x＝u_y＝u_z＝0;x＝w,u_y＝w*0.001  (4e)

第6组：y＝0,u_x＝u_y＝u_z＝0;y＝l,u_y＝l*0.001  (4f)

步骤7：计算出复合材料等效弹性参数，依次计算出步骤6所述边界条件下的单元阵列的应力分布和应变分布，然后采用方程（5）计算出平均应力和平均应变并写成公式（6）所示的向量形式：

σ‾ij=1V∫Ωσijdv,ϵ‾ij=1V∫Ωϵijdv---(5)]]>

σ‾11σ‾22σ‾33σ‾23σ‾13σ‾12T,ϵ‾11ϵ‾22ϵ‾33ϵ‾23ϵ‾13ϵ‾12T---(6)]]>

将六组边界条件的计算结果填入矩阵（7），即可得到复合材料的等效刚度矩阵对求逆后得到复合材料的等效柔度矩阵最后由等式(8)可计算出复合材料的等效弹性参数

S‾=[sij],i,j=1~6]]>

E₁₁＝1/s₁₁,E₂₂＝1/s₂₂,E₃₃＝1/s₃₃

μ₁₂＝-E₁*s₁₂,μ₁₃＝-E₁*s₁₃,μ₂₃＝-E₂*s₂₃

G₂₃＝1/s₄₄,G₁₃＝1/s₅₅,G₁₂＝1/s₆₆      (8)。

2.如权利要求1所述的纤维增强复合材料等效弹性参数的计算方法，其特征在于：所述步骤3中如果没有开展独立的实验来确定colorM和colorF，则通过在color[]数组中寻找最大灰度值colorMax和最小灰度值colorMin，认为colorM=colorMax，colorV=colorMin，在材料内部没有其他杂质时，打开一个断层图片，锁定一个纤维束区域，然后寻找纤维束区域的最大值和最小值，然后取其平均值就是colorF。