本发明属于复合材料健康监测技术领域，具体涉及变截面三维五向编织复合材料构件的单胞区域划分方法。
背景技术：
三维五向编织复合材料具有较高的比强度、比刚度、冲击韧性和结构可设计性等优势，因此被广泛应用于多个领域，但是其应用的前提是进行充分的相关力学性能的预测，一般有试验方法、理论预测方法、宏观均质的有限元方法等。由于近年来计算机水平的提高以及商业有限元软件的发展，有限元方法普遍借助商业有限元软件进行三维五向编织复合材料的力学性能仿真预测，其普遍采用多尺度模型有限元分析方法，将细观单胞的力学性能赋给宏观构件，从而进行均质化的宏观构件力学响应分析。但是，仅选用内胞的力学性能代表整体构件的力学性能，忽略面胞和角胞的影响，存在宏观构件力学响应预测结果误差较大的问题。此外，对于具有变截面形状的三维编织复合材料构件的动态损伤分析，若是采用均质化的弹性性能作为损伤分析的基础，则会存在损伤分析结果准确性不足的问题，以此为基础的数值模拟结果对后续工程应用的帮助较为有限。
技术实现要素：
本发明的目的在于提供变截面三维五向编织复合材料构件的单胞区域划分方法，可针对具有变截面形状的三维五向编织复合材料构件进行精度较高的力学响应预测，得到构件在宏观尺度上的应力场、应变场的分布，指导后续构件的优化设计；此外，该方法在有限元软件abaqus中进行了二次开发，增加了通用性，利于工程问题的批量化处理。
本发明所采用的技术方案是：变截面三维五向编织复合材料构件的单胞区域划分方法，包括以下步骤
步骤1：将具有变截面形状的三维五向编织复合材料构件几何模型导入有限元软件abaqus，导入时几何模型坐标系需满足垂直横截面的方向与z轴一致；
步骤2：使用abaqus软件实现构件几何模型的四面体网格划分，网格单元尺寸应小于花节高度h；
步骤3：给定三单胞区域划分方式，根据所选纤维束线密度、体积密度、填充因子、内部编织角等参数确定单胞模型相关参数，以花节高度值h作为构件z轴方向相邻两截面划分差值的尺寸值，进行自上而下顺序的截面划分，将得到的截面依次标号为si，i表示截面编号，为正整数，取值范围为1～截面划分数量；
步骤4：建立划分截面上轮廓边缘的矩形包络框，借助软件内全局坐标系在该矩形包络框所在平面建立局部坐标系，以该矩形中心为坐标原点，取平行于全局坐标系x轴的方向为x轴，平行于全局坐标系y轴的方向为y轴，建立局部坐标系oli；并且读取截面上轮廓边缘的长度值l，以1/100l作为一个单元尺度，等距划分轮廓边缘，选取构件截面轮廓边缘与第一、第四象限间x轴的交点作为划分起始节点，从而得到该截面内所有划分节点，并且基于已建立的局部坐标系oli，从起始节点开始逆时针旋转依次对截面内所有节点进行编号排序，记作
步骤5：对于非带孔截面从外轮廓边缘等距离偏移wc得到内胞截面区域与面胞角胞截面区域的分界线，具体方法为对步骤4中节点坐标数据数组
步骤6：将相邻两截面上节点编号数字相同的节点连接起来，外分界线上节点相互连接，内分界线上节点相互连接，这样得到多组连线，外分界线对应节点连线记作
本发明的特点还在于，
所述的变截面三维五向编织复合材料构件的单胞区域划分方法，确定花节高度h与所选纤维束线密度λ、体积密度ρ、填充因子ε及结构设计的编织角等参数间的数学关系；步骤3的具体步骤为：在对三维五向编织复合材料样件横截面扫描的基础上，给出单胞划分方式，并且建立三维五向编织复合材料单胞模型相关参数间的数学关系式：
其中，sa、sb分别为轴纱、编织纱截面面积；λa、ρa、εa分别为轴纱线密度、体积密度、填充因子，λb、ρb、εb分别为编织纱线密度、体积密度、填充因子；δ、k、l分别为编织纱截面顶角角度、宽度、高度；e为轴纱截面边长；γ为内部编织角；h为花节高度值h。
本发明的有益效果是：本发明变截面三维五向编织复合材料构件的单胞区域划分方法克服了以往借助有限元软件分析三维五向编织复合材料构件力学响应时，仅以内胞力学性能代表构件整体力学性能，忽略了面胞和角胞的影响，得到的宏观构件力学响应准确性不足的问题；此外，使用该方法进行构件强度预测、损伤分析等可得到更准确有效的结果。
附图说明
图1是本发明变截面三维五向编织复合材料构件的单胞区域划分方法所选用的三单胞区域划分方式示意图；
图2是本发明以花节高度值h作为划分构件截面尺度的构件截面划分示意图；
图3是本发明轮廓边缘节点(
图4是本发明轮廓边缘节点偏移方式的示意图(图中内分界线、外分界线本应为首尾依次相连的多段线段，此处为简化仅表示为一条线)；
图5是本发明相邻截面分界线上节点编号相同的节点连线的示意图，此处仅示意了2条连线，具体为节点
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明是变截面三维五向编织复合材料构件的单胞区域划分方法具体如下：
步骤1：将用于有限元仿真分析的某型高温环境下工作的结构构件几何模型导入商业有限元软件abaqus，导入后使得该模型的编织方向与abaqus软件中坐标系z轴一致，便于后续该构件模型横截面的提取；
步骤2：在abaqus设置网格类型为c3d4，且保证网格单元尺寸应小于花节高度h，对构件几何模型划分网格单元；构件网格单元模型如图1所示；
步骤3：确定单胞划分方式见图2；根据对三维五向编织复合材料样件横截面电镜扫描图像建立的单胞模型几何参数与编织角、纤维束相关参数的数学关系，通过纤维束相关参数及编织角求得单胞几何模型相关参数单胞模型横截面尺寸wi、ws、wc、ti、ts、tc以及编织花节高度h等；本例某型高温结构构件编织纤维束和轴向纤维束均选用t300-6k，线密度为1.76g/m，体积密度为0.8g/cm3，填充因子为0.75，且内部编织角为35.26°；根据已构建的数学关系求得wi＝ws＝2.88mm，wc＝0.72mm，ti＝2.88mm，ts＝tc＝0.72mm，h＝4.08mm；依次进行横截面划分位置的选取，并且建立所有截面的局部坐标系，并将得到的截面依次标号为si(i取1～划分截面的数量)，图2示意了划分后的部分截面(s1、s2)；
步骤4：调用步骤3中的轮廓边缘节点坐标，借助软件内全局坐标系，对全局坐标系下x坐标最大值节点、x坐标最小值节点作过该点垂直于x轴的垂线，对全局坐标系下y坐标最大值节点、y坐标最小值节点作过该点垂直于y轴的垂线，四条垂线组成截面上轮廓边缘的矩形包络框；以该矩形中心为坐标原点，与全局坐标系x轴平行方向作为x轴，与全局坐标系y轴平行方向作为y轴建立局部坐标系oli；并且读取截面上轮廓边缘的长度值l，以1/100l作为一个单元尺度，等距划分轮廓边缘，选取构件截面轮廓边缘与第一、第四象限间x轴的交点作为划分起始节点，从而得到该截面内所有划分节点，并且基于已建立的局部坐标系oli，从起始节点开始逆时针旋转依次对截面内所有节点进行编号排序，记作
步骤5：对于非带孔截面，仅存在外轮廓边缘，对步骤4中外轮廓边缘划分节点的坐标数据数组
步骤6：将相邻两截面偏移后的对应节点按照编号顺序连接起来(例如，内分界线上节点
1.变截面三维五向编织复合材料构件的单胞区域划分方法，其特征在于，
步骤1：将具有变截面形状的三维五向编织复合材料构件几何模型导入有限元软件abaqus，导入时几何模型坐标系需满足垂直横截面的方向与z轴一致；
步骤2：使用abaqus软件实现构件几何模型的四面体网格划分，网格单元尺寸应小于花节高度h；
步骤3：给定三单胞区域划分方式，根据所选纤维束线密度、体积密度、填充因子、内部编织角参数确定单胞模型相关参数，以花节高度值h作为构件z轴方向相邻两截面差值的尺寸值，并且进行自上而下顺序的截面划分，将得到的截面依次标号为si，i表示截面编号，为正整数，取值范围为1～截面划分数量；
步骤4：建立划分截面上轮廓边缘的矩形包络框，借助软件内全局坐标系在该矩形包络框所在平面建立局部坐标系，以该矩形中心为坐标原点，取平行于全局坐标系x轴的方向为x轴，平行于全局坐标系y轴的方向为y轴，建立局部坐标系oli；并且读取截面上轮廓边缘的长度值l，以1/100l作为一个单元尺度，等距划分轮廓边缘，选取构件截面轮廓边缘与第一、第四象限间x轴的交点作为划分起始节点，从而得到该截面内所有划分节点，并且基于已建立的局部坐标系oli，从起始节点开始逆时针旋转依次对截面内所有节点进行编号排序，记作
步骤5：对于非带孔截面从外轮廓边缘等距离偏移wc，即面胞角胞厚度值，得到内胞截面区域与面胞角胞截面区域的分界线；具体方法为：
对步骤4中坐标值数组
步骤6：将相邻两截面上节点编号数字相同的节点连接起来，外分界线上节点相互连接，内分界线上节点相互连接，外分界线对应节点连线记作
2.根据权利要求1所述的变截面三维五向编织复合材料构件的单胞区域划分方法，其特征在于，确定花节高度h与所选纤维束线密度λ、体积密度ρ、填充因子ε及结构设计的编织角参数间的数学关系；步骤3的具体步骤为：在对三维五向编织复合材料样件横截面扫描的基础上，给出单胞划分方式，并且建立三维五向编织复合材料单胞模型相关参数间的数学关系式：
l＝k 2，
其中，sa、sb分别为轴纱、编织纱截面面积；λa、ρa、εa分别为轴纱线密度、体积密度、填充因子，λb、ρb、εb分别为编织纱线密度、体积密度、填充因子；δ、k、l分别为编织纱截面顶角角度、宽度、高度；e为轴纱截面边长；γ为内部编织角；h为花节高度值h。
技术总结