2021年
第1期
【摘要】综合考虑机织复合材料成型及力学性能特点，设计了一体式结构3D 机织碳纤维复合材料避震塔，并基于多尺
度优化方法，综合考虑材料和结构参数，实现避震塔强度性能优化设计。为获取准确的仿真分析用材料参数，基于显微镜观测得到的机织参数建立多尺度预测模型。通过避震塔强度性能仿真分析与试验结果的对比，验证了考虑工艺过程影响的材料参数预测结果及结构性能仿真分析方法的正确性，与原高强钢避震塔相比，优化得到的3D 机织碳纤维复合材料避震塔实现了50%的轻量化效果。
主题词：3D 机织碳纤维复合材料避震塔强度性能多尺度优化设计一体式结构
中图分类号：U465.6
文献标识码：A DOI:10.19620/jki.1000-3703.20191334
Multi-Scale Optimization Design on Strength Performance of Carbon
Fiber Reinforced Composite Shock Absorption Tower
Lu Jiahai,Wang Xia,Kou Hongbin
（Pan Asia Technical Automotive Center Co.,Ltd.,Shanghai 201201）
【Abstract 】An integrated structure 3D woven carbon fiber reinforced composite shock absorption tower is designed considering molding and mechanical property characteristics of composite material,and multi-scale optimization is implemented on strength performance considering material and structure parameters simultaneously.Multi-scale prediction model of 3D woven carbon fiber reinforced composite is established based on weaving parameters obtained through electron microscope observation for accurate material parameter input of simulation analysis.Correlation results between strength test and simulation results show correctness of the modeling method.The optimized 3D woven carbon fiber reinforced
composite shock tower achieves 50%weight reduction comparing to original high strength steel structure.
Key words:3D woven carbon fiber reinforced composite,Shock absorption tower,Strength
performance,Multi-scale optimization design,Integrated structure
卢家海
汪霞
寇宏滨
（泛亚汽车技术中心有限公司，上海201201）
碳纤维复合材料避震塔强度性能多尺度优化设计
汽车技术·Automobile Technology
1前言
碳纤维复合材料因比刚度、比强度和比吸能高等特点，广泛应用于汽车零部件，在保证产品性能的同时，可实现汽车的轻量化设计[1-2]，但是碳纤维复合材料复杂的力学特性提高了产品设计难度。同时，碳纤维复合材料具有材料-结构一体成型的特性，相对于传统金属材料主要通过结构优化满足产品性能的特点，碳纤维复合材料可以在材料参数和结构参数两方面开展设计，得到最优的材料和结构组合，将产品的轻量化效果最大化。
对于层铺碳纤维复合材料，其力学性能与铺层方式有密切联系，许多学者在铺层参数优化方面开展了大量的研究工作[3-4]。对于机织碳纤维复合材料，机织参数
是影响其宏观力学性能并进而影响产品宏观结构性能的重要因素。结合多尺度建模和结构优化设计方法开展复合材料零部件产品多尺度设计，是实现机织碳纤维复合材料产品材料-结构一体设计的重要途径[5-6]。当前针对机织复合材料结构的多尺度优化设计研究中，对于介观尺度机织结构，均基于理想单胞进行多尺度建模。而基于介观观测结果可知，在成型过程中，单胞形态会发生变化，导致复合材料中实际机织参数与理想机织参数发生偏离[7-8]。为了更好地预测机织复合材料的宏观力学性能和产品宏观结构性能，需要在多尺度设计过程中考虑工艺过程对机织参数的影响。
本文旨在综合考虑机织复合材料的成型工艺和力学性能特点，设计3D 机织碳纤维复合材料汽车避震塔
【引用格式】卢家海，汪霞，寇宏滨.碳纤维复合材料避震塔强度性能多尺度优化设计[J].汽车技术,2021(1):51-56.
LU J H,WANG X,KOU H B.Multi-Scale Optimization Design on Strength Performance of Carbon Fiber Reinforced Composite Shock Absorption Tower[J].Automobile Technology,2021(1):51-56.
-
-
51
结构，并基于多尺度优化方法，在满足避震塔强度性能的基础上实现其轻量化设计。为考虑成型工艺对复合材料的影响，通过显微镜观测得3D机织碳纤维复合材料介观机织参数，并基于计算细观力学理论，建立考虑工艺过程影响的宏观各向异性多尺度预测模型及其本构模型。考虑到避震塔的载荷特性及复合材料成型的特点，设计一体式结构3D机织碳纤维复合材料避震塔，并进行强度性能试验和仿真分析对标。在此基础上，基于多尺度建模和优化设计理论，综合考虑介观机织参数及避震塔宏观结构参数，对避震塔强度性能进行多尺度优化设计。
2力学性能试验及机织参数观测
2.1样件制作及力学性能试验
本文所采用的3D机织碳纤维复合材料首先由纱线编织为预制件，之后通过树脂传递模（Resin Transfer Molding，RTM）工艺与树脂复合加工成复合材料样板。编织预制件采用的纱线为东丽T700-6k，机织结构为三维正交，理想单胞如图1所示，具体机织参数如表1所示。复合材料成型所用树脂为环氧树脂。最终加工得到的样板厚度为7.7mm，碳纤维体积分数为39%。
图13D正交机织理想单胞示意
表13D机织碳纤维复合材料机织参数
本文分别参照标准ASTM D638、ASTM D6641和ASTM D5379，规划3D机织碳纤维复合材料经向和纬向的拉伸和压缩力学性能试验，以及面内剪切试验，结果如表2所示。由表2可知，3D机织碳纤维复合材料面内力学性能具有典型的各向异性和拉压非对称性。
2.2样件机织参数观测
图2所示为样件经向和纬向断面的显微镜观测结果，由图2可知，加工过程对样板中纱线的形状和位置产生了影响，经向纱和纬向纱均呈平行四边形，Z向纱沿厚度方向发生了偏转，角度为55°。
表23D机织碳纤维复合材料面内力学性能试验结果
（a）经向断面示意
（b）纬向断面示意
图2样件断面显微观察结果
最终测量得到的介观尺度参数，即经向纤维束、纬向纤维束和Z向纤维束的截面尺寸如表3所示。
表3不同纤维束介观尺寸参数测量结果mm
3宏观力学性能预测及本构建模
3.1复合材料多尺度建模
由显微镜观测结果可知，由于加工过程中模流的作用，实际样板中的单胞形式及参数与理想单胞均有偏差。为更加准确地预测实际产品中3D机织碳纤维复合材料的宏观力学性能，本文基于显微镜观察得到的机织参数，分别建立细观和介观尺度模型。
基于表3测量的纤维束截面尺寸，可以分别计算细观尺度纤维束的截面面积S T。东丽T700碳纤维直径为6.8μm，每根纤维束内含有6000根纤维，每根纤维束中
碳纤维的截面积为S f，则细观尺度参数，纤维束中碳纤
经纱
Z向纱
纬纱
纬向
经向
机织参数经纱间距/mm 纬纱间距/mm 层数/层取值
1.0
3.3
8
加载方式
经向拉伸
经向压缩
纬向拉伸
纬向压缩
面内剪切
强度/MPa
638.4
433.4
189.9
215.9
49.9
弹性模量/GPa
62.5
54.6
23.1
19.2
3.4
Z向纱
经纱
经纱
纬纱
纱线方向
经向
纬向
Z向
高度
0.52
0.30
0.28
宽度
0.72
1.71
1.08
间距
1.00
0.33
1.00
-
-52
汽车技术
2021年
第1期
维的体积分数V f 为：
V f =S f /S T （1）
最终得到复合材料中经向、纬向和Z 向纤维束中碳纤维体积分数分别为63.09%、47.90%和79.66%。
基于表3和计算得到的纤维束中碳纤维的体积分数，分别对3D 机织碳纤维复合材料细观尺度和介观尺度进行几何模型和有限元网格模型重建，最终得到的细
观及介观尺度代表性单胞有限元模型分别如图3和图4所示。在细观尺度中，假设碳纤维六边形排布于纤维束内部。基于重建得到的模型，计算得到复合材料中碳纤维体积分数为38.5%，与样板中碳纤维体积分数的39%基本一致。
图3
细观尺度代表性单胞有限元模型
图4介观尺度代表性单胞有限元模型
3.2预测结果及对标
在细观尺度预测中，碳纤维和环氧树脂力学性能参
数均由材料生产商提供，具体数值如表4所示，其中，E 1T 、E 1C 分别为纤维轴向方向的拉伸和压缩弹性模量，E 2为垂直于纤维轴向的弹性模量，G 12、G 23分别为面内和面外剪切模量，v 12为面内泊松比。在仿真分析过程中，碳纤维力学性能采用横观各向同性本构模型描述，并考虑拉压非对称性的影响，仿真分析中采用ABAQUS 的用户
自定义（User-Defined Material ，UMAT ）子程序进行数值计算，环氧树脂采用德鲁克-普拉格（Drucker-Prager ）准则描述其力学行为。
基于图3所建立的细观尺度有限元模型和表4中的组成相力学性能参数，预测得到不同纤维束的力学性能参数如表5所示。其中，v 23为面外泊松比，X T 、X C 分别为纤维轴向拉伸和压缩强度，Y T 为垂直于纤维轴向的拉伸强度，S L 、S T 方分别为面内和面外的剪切强度。
表4碳纤维和环氧树脂力学性能参数
表5
预测得到的不同方向纤维束力学性能参数结果
为预测3D 机织碳纤维复合材料宏观强度性能，介
观尺度仿真分析中需要建立合适的纤维束强度准则。参考目前针对单向碳纤维复合材料的研究结果[9]，基于Pinho 准则描述纤维束复杂应力状态下的强度性能，并采用ABAQUS 的UMAT 子程序将上述准则数值化。将表5中的细观预测结果带入到介观尺度模型，预测得到的3D 机织碳纤维复合材料宏观力学性能结果如表6所示。与试验结果相比，建立的考虑工艺过程影响的多尺度预测方法可以很好地预测3D 机织碳纤维复合材料轴向拉伸及压缩载荷下的宏观力学性能。
表6
3D 机织碳纤维复合材料宏观力学性能预测结果
3.33D 机织碳纤维复合材料本构建模
考虑到3D 机织碳纤维复合材料宏观力学性能的拉
压非对称性和各向异性，建立其宏观本构模型：
参数E 1T /GPa E 1C /GPa
E 2/GPa v 12当当改装网
纤维23019915
0.2基体3
0.35
参数G 12/GPa G 23/GPa 拉伸强度/MPa
压缩强度/MPa
纤维24
5.03
49002400
基体
72
103
参数
E 1T /GPa E 1C /GPa E 2/GPa v 12v 23
G 12/GPa 经向145.6126.687.660.240.43
4.57
纬向
111.2196.8846.350.270.403.26
Z 向
182.82159.2110.270.220.467.46参数
G 23/GPa
X T /MPa X C /MPa Y T /MPa S L /MPa S T /MPa
经向2.68
3051.541554.2556.6953.0261.94
纬向2.14
2350.121205.8752.2251.2356.70
Z 向3.59
3841.711934.7263.5256.5366.77
加载方式经向拉伸
纬向拉伸
经向压缩
纬向压缩
面内剪切
力学性能模量/GPa 强度/MPa 模量/GPa 强度/MPa 模量/GPa 强度/MPa 模量/GPa 强度/MPa 模量/GPa 强度/MPa
试验值62.5634.823.1189.954.6433.419.2215.93.4
49.9
预测值60.5719.522.2218.553.5455.020.0210.73.2
41.3
误差/%-3.412.7
-4.015.1-1.94.84.0-2.5-5.9
-17.2-
-53
汽
车
技
术
ìíîïïüýþïïε1ε2ε12=ìíî
ïïïïïïïï
üýþ
ïïïïïïïï1
E W
-
v 12E W 0-v 12E W 1E F 000
12G WF ìíîïïüý
þïïσ
1σ2σ12（2）
式中，ε1、ε2分别为经向和纬向的线应变；ε12为剪应变；E W 、E F 分别为表6中所示经向和纬向的弹性模量；G WF 为表6中所示面内剪切模量；σ1、σ2分别为经向和纬向的正应力；σ12为剪应力。
考虑到材料力学性能的拉压非对称性，当σi ≥0时，
E W 、E
F 分别取表6中所示的拉伸弹性模量，当σi <0时，E W 、E F 分别取表6中所示的压缩弹性模量。
对于其强度性能，采用最大应力准则进行建模。所
建立的本构模型基于ABAQUS 的UMAT 子程序数值化。
4避震塔结构设计及强度性能对标
4.1
避震塔结构设计
考虑到3D 机织碳纤维复合材料的加工及避震塔所承受的载荷特性，设计了如图5a 所示的一体式结构避震塔，避震塔的高度方向为机织过程中的经向，环向为机织过程中的纬向。沿避震塔高度方向，该结构分为上部、中部、下部3个区域，如图5b 所示。在初始设计中，3个区域采用相同的经密和纬密，但机织层数不同，如表7所示。基于上述机织参数，首先将碳纤维纱线编织为预制件，之后通过RTM 工艺与环氧树脂结合，形成一体式结构碳纤维复合材料避震塔本体。在避震塔的加工过程中，所采用的工艺参数与样件加工所用参数完全相同。
（a ）一体式结构避震塔示意
（b ）避震塔剖视图
图5
一体式结构避震塔及其剖视图
表7
一体式结构避震塔机织参数
4.2避震塔强度性能试验
避震塔强度性能试验台架如图6所示，其强度性能
要求为载荷加载至56.1kN 时，避震塔不发生破坏。
图6避震塔强度性能试验台架
试验过程中，当载荷加载至40.3kN 时，避震塔发生破坏，其失效形式如图7所示，避震塔首先沿顶部圆孔的纬向发生破坏，
之后沿经向进行扩展。
图7避震塔强度破坏结果
4.3避震塔强度性能仿真分析及对标
建立的避震塔强度性能有限元分析模型采用共节
点处理不同厚度的连接关系。考虑到3D 机织碳纤维复合材料的各向异性，在避震塔有限元模型中建立局部坐标系描述材料方向。在40.3kN 加载下，避震塔CAE 分析得到的应力云图如图8所示。由图8可知，避震塔顶部圆孔周围的纬向拉应力达到拉伸强度极限，继续加载，将会发生沿经向的裂纹。
图8避震塔强度仿真分析应力结果示意
由避震塔强度性能对标结果可知，所建立的材料参
数预测流程、本构模型及CAE 分析方法可以很好地描述3D 机织碳纤维复合材料的强度力学性能。
机织参数
经纱间距/mm 纬纱间距/mm 上部层数/层
参数1.03.38
机织参数中部层数/层下部层数/层
参数64
经向
纬向
上部中部
下部
裂纹
纬向应力/MPa
181.3126.371.316.3-38.7
-
-
54
2021年
第1期
5避震塔多尺度优化设计
考虑到3D 机织复合材料机织参数的可设计性，针对避震塔强度性能，分别选择避震塔上部、中部和下部3个区域的经纱间距、纬纱间距和层数等介观机织参数，以及顶部配合面圆孔直径等宏观结构参数作为设计变量，进行多尺度优化设计，具体变量及其范围如表8
所示。多尺度优化设计流程如图9所示。
表8多尺度优化设计变量及其范围
图9复合材料避震塔多尺度优化设计流程
针对介观尺度设计变量，基于试验设计（Design of Experiment ，DOE ）结果，采用前述的多尺度建模方法，分别建立相应的介观尺度预测模型，计算得到与参数相对应的刚度和强度性能宏观材料参数。通
过将上述参数带入到避震塔强度性能CAE 仿真分析模型中，并结合宏观尺度设计变量，最终得到避震塔宏观质量以及经向和纬向等应力性能结果。在仿真分析得到的避震塔宏观性能的基础上，基于Kriging 建模方法，建立设计变量与避震塔强度和质量性能之间的近似模型。并基于此，建立如下的优化问题：
min M
<
σ1
σ1T,Stren ≤0.7,σ1≥0，
-σ1
σ1C,Stren ≤0.7,σ1<0，σ2
σ2T,Stren ≤0.7,σ2≥0，-σ2
σ2C,Stren
≤0.7,σ2<0，|
|σ12σ12,Stren
≤0.7
（3）
式中，M 为避震塔质量，σ1T,Stren 、σ1C,Stren 分别为经向拉伸、压缩强度；σ2T,Stren 、σ2C,Stren 分别为纬向拉伸、压缩强
度；σ12,Stren 为面内剪切强度。
基于式（3），采用遗传算法进行求解，最终优化得到
的参数组合如表9所示。由表9可知，优化得到的机织参数经纱间距有所增加，纬纱间距减小，提高了复合材料纬向的承载能力。
表9
优化后的设计参数与原设计参数对比
基于表9中的参数，建立优化后的避震塔强度性能仿真分析模型，最终得到的强度因子与原结构强度因
子对比结果如表10所示。由表10可知，优化后的结构满足避震塔强度性能要求。最终得到的3D 机织碳纤维复合材料避震塔的质量为2.28kg ，相对于原高强钢避震塔，实现了50%的轻量化效果。
表10优化后避震塔强度性能仿真分析结果对比
6结束语
本文综合考虑机织复合材料成型及力学性能特点，设计了一体式结构3D 机织碳纤维复合材料避震塔，并基于多尺度优化方法，综合考虑材料和结构参数，实现了避震塔强度性能优化设计。优化得到的复合材料避震塔相对于原高强钢避震塔实现了50%的轻量化效果。
变量名称
介观机织参数
宏观结构参数
经纱间距x 1/mm
上部纬纱间距x 2/mm
中下部纬纱间距x 3/mm 上部经纱层数x 4/层中部经纱层数x 5/层下部经纱层数x 6/层
上部圆孔直径x 7/mm
变量范围[1.2,2.2][1.5,2.5][2.0,3.3][8,11][6,9][4,7]
[45,65]
介观尺度有限元建模介观尺度几何建模
介观机织参数
DOE
多
尺度建模
宏观力学性能参数
宏观结构参数
避震塔宏观几何建模
避震塔宏观有限元建模
避震塔强度性能仿真基于Kriging 的近似建模优化问题定义及求解
结果验证
结束
N
开始Y
设计变量x 1/mm x 2/mm
x 3/mm x 4/层
优化结果1.21.62.210
原设计参数1.03.33.38
设计变量x 5/层x 6/层
x 7/mm 优化结果
75
45原设计
参数
64
66载荷方向经向拉伸经向压缩纬向拉伸纬向压缩面内剪切
初步设计强度因子
0.830.611.330.330.23
优化结果强度因子
0.660.440.670.160.11
-
-55
发布评论