三维四向编织复合材料刚度预报.docVIP

三维四向编织复合材料刚度预报 方国东，梁军(，王宝来 （哈尔滨工业大学复合材料与结构研究所 ，哈尔滨， 150001） 摘要：本文利用刚度均匀化模型，对三维四向编织复合材料的刚度进行模拟。分为三个层次对单胞进行刚度均匀化：从纤维束到子单胞对角微段，到子单胞，再到整体单胞。本文得到三维四向编织复合材料的刚度预报结果与文献的实验结果相符。同时本文讨论不同几何参数以及纤维束在不同弯曲状态下，对材料刚度的影响。 关键字：编织复合材料，均匀化方法，刚度 三维编织复合材料较传统层合复合材料具有优良的抗分层能力、高的损伤容限和抗冲击能力，广泛应用到航空航天结构中的承力部位。由于三维编织复合材料具有复杂的纤维构造，编织工艺参数、结构参数等诸多因素影响结构的力学性能，所以目前的研究工作者主要对三维编织复合材料的刚度进行预测，而对强度的预测还很少。 Ma, Yang and Chou提出了弹性应变能法[1] 和纤维倾角模型[2] 。Wu[3] 提出了表征编织结构的三胞模型。Wang和Wang[4] 和陈利等[5] 分别通过跟踪内部纤维和最小二乘法拟合得到编织件的拓扑结构模型。文献[6-10]利用选择平均化方法和有限元对三维四向编织复合材料进行刚度分析。庞宝君和曾涛等[11][12] 利用混合有限元对三维四向编织复合材料进行损伤和强度分析。刚度的预测是对材料很好的强度预测的前提。本文利用刚度均匀化方法，对三维四向编织复合材料进行刚度分析得到工程材料常数，与文献试验结果吻合。同时考虑几何参数以及纤维束在单胞中弯曲对材料刚度的影响。 1．单胞的模型和假设 （1）单胞的几何模型 三维四向11编织体的纤维束按四个方向对称分布，如图1 (a) 所示。取含有四向的八个子单胞为一个代表单胞，如图1（b、c）所示。由于编织件中纤维束在一定的牵引作用下编织，不同的牵引力得到的编织件中的纤维束的轨迹不同。在足够大的牵引力作用下，纤维束取直；在小牵引力的作用下，纤维束还具有一定的弯曲[13] 。在编织过程中纤维束的牵引和复合材料成型的过程中要对编织件填充基体并施加一定的压力，纤维束的横截面积在挤压力的作用下不再是理想的圆形[14] 。本文假设纤维的横截面积为椭圆形，如图1(d)所示。本文采用文献[14]中的体胞的几何关系。，，其中分别是纤维束椭圆形横截面的长轴和短轴，为一个单胞编织节的高度，为内部编织角。所以1/8个子单胞的边长为。 （2）单胞的分析模型 如图2所示，取1/4单胞。因为考虑纤维束的弯曲影响，并且能够清晰的利用数学描述单胞分析模型，建立三个层次的坐标系：子胞坐标系；子胞对角坐标系；子单胞对角微段坐标系。坐标系之间利用纤维束偏离角和内部编织角旋转得到。模型首先对子单胞对角微段进行刚度平均化，然后对子单胞纤维束进行刚度平均化，最后对子单胞进行刚度平均化。本文假设在坐标系中弯曲的纤维束利用描述，其中为子单胞对角线的长度，为纤维束弯曲的幅值。纤维束曲线的幅值大小在[0,]的范围内变化。 图表 1 三维四向编织体几何结构 图表 2 三维四向编织结构单胞分析模型 2．刚度均匀化方法和有效弹性模量的预测 （1）刚度均匀化方法 基于等应变和等应力假设，把不同组分材料的应力项和应变项分别分成两组，三个等应变项和三个等应力项，即 （1） 其中n代表等应变项，s代表等应力项，k代表不同的组分材料（纤维束和基体）。从而组分材料的刚度矩阵相应分为四个矩阵项，即 其中， （2） 所以组分材料的本构方程可以写为等应变项和等应力项的组合形式 （3） （4） 刚度平均化后整体应力和应变为 ，，， （5） 其中，和分别是刚度平均后整体应力的等应力部分和整体应变的等应变部分。和分别为刚度平均化后整体应力的非等应力部分和整体应变的非等应变部分。为组分材料的体积分数，有由于。 利用（4）式得 （6） 代入（3）式得 （7） 把式（6）和式（7）代入式（5）得 ， （8） 其中 由式（5）和式（8）得到均匀化后整体的刚度矩阵 （9） （2）有效弹性模量的预测 本文由纤维束与基体的刚度得到代表单胞的刚度需要通过三次刚度均匀化，每个层次的等应力项和等应变项的选取不同。 从纤维束和基体到子单胞对角微段的刚度均匀化。