摘要
为了探究复合材料波纹腹板梁在坠撞过程中的吸能性能,对复合材料波纹腹板梁试验件进行了动态冲击压溃试验研究,得到了试验件的损伤破坏形貌以及载荷和能量的历程曲线。基于ABAQUS/Explicit平台二次开发得到了模拟复合材料波纹腹板梁冲击压溃过程的仿真分析模型,模型采用了改进的Hashin损伤判定准则和Choi‑Chang准则综合判断单元失效,并结合Cohesive界面单元,可较为真实地反映所研究复合材料层合结构的各向异性和渐进损伤特性。通过数值模拟得到了能量评估参数比吸能(Specific energy absorption, SEA)和平均压溃载荷,与试验结果进行了对比分析。基于仿真分析模型,本文进一步研究了不同波数的复合材料波纹腹板梁的吸能能力。数值模拟结果表明:波纹腹板梁在冲击载荷作用下发生了渐进压溃失效;平均压溃载荷的相对误差较小,验证了模型的有效性;在腹板单波长度不变的情况下,波纹腹板梁的长度对抗冲击吸能的影响较大。当波纹腹板梁长度较小时,结构容易失稳,无法有效地吸收能量。
复合材料具有结构可设计性,与组分材料相比,有更高的比强度和比模量(刚度),更广泛的使用温度范围,优良的化学稳定性(陶瓷基复合材料以及树脂基复合材料具有良好的耐腐蚀性),良好的抗疲劳、冲击和蠕变性能以及较高的断裂韧性等。在航空航天领域,复合材料结构已成为直升机抗坠毁结构首选的吸能材料,被广泛运用于机身底部的地板结构中,如

图1 波纹腹板梁吸能结构
Fig.1 Corrugated web beam in energy absorption structure
复合材料波纹腹板梁作为飞行器吸能结构的重要组成元件,研究其在冲击压溃过程中的平均载荷和比吸能(Specitic energy absorption,SEA)等吸能参数,对提高飞行器的抗坠毁能力具有重要指导意义。目前,国内外研究人员对复合材料吸能结构的研究多集中于准静态加载的条件
国内外研究人员提出了多种数值模拟方法对复合材料吸能结构耐撞性能进行研究。苏
还有研究表明,复合材料抗坠毁结构中不同形式的薄弱环节对构件的吸能水平有较大影响,对薄弱环节进行优化设计可改进结构的吸能形式,提高耐撞性
本文所研究的波纹腹板梁结构形式如

图2 复合材料波纹腹板梁示意图
Fig.2 Schematic diagram of composite beam with corrugated web
波纹腹板梁的铺层形式为[(±45)3/(0, 90)/ (±45)3],铺层所用的碳纤维单向材料属性由
密度ρ/(g⋅c | 拉伸模量/GPa | 剪切模量/GPa | ||||
---|---|---|---|---|---|---|
2.53 | 120 | 9.5 | 9.5 | 5 | 5 | 3.2 |
强度/MPa | 泊松比 | |||||
1 000 | 40 | 700 | 180 | 70 | 0.31 | 0.52 |
N/MPa | S/MPa | T/MPa |
---|---|---|
21 | 28 | 28 |
GⅠC/(J⋅ |
GⅡC/(J⋅ |
GⅢC/(J⋅ |
306 | 632 | 632 |
铺层形式的波纹腹板梁中,预置的薄弱环节有效地降低冲击峰值载荷,明确了破坏起始位置(薄弱环节处),并能参与引导波纹腹板梁的压溃过
在数值研究中所建立的结构等效模型中,对波纹腹板梁的薄弱环节也进行了等效处理。将厚度为的薄弱环节(圆弧区),简化为自腹板下端竖直向上高度为15 mm的等效薄弱区域,并对该区域的材料刚度进行折减等效。假设薄弱环节处的等效材料与旁边的复合材料层板具有相同的轴向压缩刚度,当施加一定的轴向压缩载荷,圆弧形薄弱环节与等效薄弱环节的压缩变形量相同,由此可求得长15 mm的等效薄弱环节的厚度值。
(1) |
式中:L代表波纹腹板梁宽度,为等效薄弱环节的厚度,为数值模型中等效薄弱环节的高度,r代表实际圆弧形薄弱环节的半径。已知所研究的圆弧形薄弱环节,则通过
采用数值方法对波纹腹板梁的抗冲击性能进行仿真分析时,对模型底端高度为15 mm的区域进行等效处理。在厚度不变的情况下,对材料刚度进行折减,折减因子为。
波纹腹板梁的坠撞过程可被视为结构的轴向冲击压溃,试验件采用碳纤维单向布铺设复合而成,试验件尺寸与铺层形式同前节,坠撞试验选用了三波波纹腹板梁试验件(

图3 三波波纹腹板梁试验件照片
Fig.3 Photo of specimen of beam with three‑wave corrugated web
复合材料波纹腹板梁坠撞试验中,试验件投放高度为2.12 m,投放质量为108 kg。试验开始后,记录试验件在坠撞过程中的载荷时间历程;用高速摄像机对试验过程进行摄像记录,位移(即试件压缩量)、碰撞后峰值速度及加速度时间历程均通过对高速摄像数据处理后得到。
基于坠撞试验数据,并通过观察波纹腹板梁的破坏现象与失效模式,分别得到峰值载荷、平均载荷、吸收能量等数据。

图4 轴向冲击载荷‑时间曲线
Fig.4 Axial impact load‑time curve

图5 吸收能量‑时间曲线
Fig.5 Energy absorption‑time curve
试验件编号 | 3A‑1 | 3A‑2 | 3A‑3 | 3A‑4 |
---|---|---|---|---|
峰值载荷/kN | 82.66 | 94.02 | 48.29 | 55.82 |
平均载荷/kN | 13.54 | 12.84 | 16.96 | 16.28 |
峰值加速度/g | 13.80 | 15.85 | 18.64 | 16.39 |
峰值速度/(m· | 6.42 | 6.41 | 6.44 | 6.45 |
压缩量L/mm | 180.65 | 189.19 | 126.85 | 130.22 |
总功量Mg/kJ | 2.45 | 2.46 | 2.39 | 2.40 |
吸收功量/kJ | 2.45 | 2.43 | 2.15 | 2.12 |
坠撞试验结果表明,复合材料波纹腹板梁的压溃破坏过程主要分为两个阶段:初始局部破坏阶段和稳定压溃破坏阶段。波纹腹板梁的抗失稳能力较高,破坏过程稳定,有较高的吸能性能。本试验有效地模拟了坠撞冲击压溃过程,对于该型复合材料层合结构具有代表意义。
层合板可视为是由多个单层板构成,把多层复合材料单层板转化为具有相同截面面积的单层板,其强度也是通过单层板的强度来预测。复合材料拥有多种失效准则,在有限元数值计算过程中,单层板达到了失效条件,材料即发生失效。通常在宏观上可以认为复合材料层合板是正交各向异性材料。
单层板在材料主方向的应力‑应变关系为
(2) |
本文通过Hashin准则来判断单层板的失效,采用Choi‑Chang准
基体拉伸损伤判定
(3) |
式中:和代表2、3方向的应力;S则表示基体剪切强度;代表基体拉伸强度;和代表剪切应力。
基体压缩损伤判定
(4) |
式中代表压缩强度。
纤维拉伸损伤判定
(5) |
式中:代表纤维纵向的拉伸强度;代表纤维纵向的应力;和代表剪切应力。
纤维压缩损伤判定
(6) |
式中:代表纤维压缩强度;代表纤维纵向的应力。
分层判定标准如下
(7) |
式中:、为上铺层的横向剪应力;代表下铺层的横向正应力。式中的参数K随铺层厚度方向的拉和压而不同,具体如下
(8) |
式中:代表界面性能参数,其值等于2;、分别表示上、下铺层沿厚度方向上的正应力。
本文除了采用Choi‑Chang准则模拟单层板内部的分层损伤,还引入了Cohesive单元,模拟单层板之间的分层损伤。Cohesive单元分层破坏判定公式为
(9) |
式中:、、分别表示模式Ⅰ、Ⅱ、Ⅲ的层间强度;、、表示层间应力。当
当复合材料内部有损伤萌生并扩展,为达到模拟材料损伤的目的,需要对模型的刚度进行折减。Camanho
失效模式 | 刚度折减系数 | |||||
---|---|---|---|---|---|---|
Exx | Eyy | Ezz | Gxy | Gxz | Gyz | |
基体拉裂 | — | 0.2 | — | 0.2 | — | 0.2 |
基体压裂 | 0.4 | — | — | 0.4 | 0.4 | — |
纤维拉裂 | 0.01 | — | — | 0.01 | 0.01 | — |
纤维压断 | 0.01 | — | — | 0.01 | 0.01 | — |
分层失效 | — | — | 0.01 | — | 0.01 | 0.01 |
本文所研究的波纹腹板梁是由16层单向复合材料铺成,单层厚度仅有0.1 mm,为保证数值计算仿真过程中的单元质量,同时提高计算效率、降低计算成本,对波纹腹板梁有限元模型进行了简化处理。基于层合板理论,将4层单向复合材料合并等效为一个实体层,厚度方向共4层实体单元,每层单元厚度为0.4 mm。在每两层实体单元之间布置一层Cohesive单元,共计7层,有限元模型如

图6 模型厚度方向单元类型示意图
Fig.6 Element types along thickness direction of FE model
由于在实际压溃过程中波纹腹板梁上、下缘条的吸能能力几乎可忽略不计,所以在仿真模型中去掉了上、下缘条。将试验平台简化为两个刚性平板,定义上平台质量为108 kg。
经过等效处理和单元离散后,波纹腹板梁数值仿真模型共130 928个单元,其中实体单元74 816个,Cohesive单元56 112个。模型如

图7 三波波纹腹板梁有限元模型
Fig.7 FE model of a composite beam with three‑wave corrugated web



图8 复合材料波纹腹板梁有限元模型计算结果及试验后试验件照片
Fig.8 Numerical results of FE model of composite corrugated web beam and a photo of a specimen after test
数值仿真计算结果显示(
仿真结果/kN | 试验结果/kN | ||
---|---|---|---|
峰值载荷 | 平均载荷 | 峰值载荷 | 平均载荷 |
86.0 | 20.9 | 84.0 | 18.6 |
(1)通过降低刚度和强度所得到的等效薄弱环节,无法完全反映圆弧形薄弱环节实际的破坏与吸能过程。
(2)模型中选用的刚度折减系数也是在基于先前折减方法的基础上确定的,不一定能完全反应真实情况。
(3)不论是波纹腹板梁与刚性板之间的接触,还是波纹腹板梁自身接触的定义,都进行了理想化假设。另外,模型只定义了摩擦因数一个参量来引入摩擦导致的能量耗散,与实际接触形式相比较为简单。
(4)模型中所采用的Cohesive单元本构需要进一步修正,层间断裂韧性会影响结构的初始损伤状态。
SE
5种波数波纹腹板梁的复合材料本构模型在受到冲击载荷时的应力云图由

图9 不同波数的波纹腹板梁有限元模型压溃应力云图
Fig.9 Stress nephogram of FE models of corrugated web beams with different wavenumbers when collapsed
5种不同波数的波纹腹板梁有限元仿真计算所得的位移载荷曲线如

图10 不同波数的波纹腹板梁模型的载荷‑位移曲线
Fig.10 Load‑displacement curves of corrugated web beams with different wavenumbers
(1)5个模型受冲击后载荷瞬间达到峰值,薄弱环节首先被压溃,随后载荷大幅下降,并逐渐趋于稳定,可视为收敛于固定值,该值与材料的断裂韧性有关。在当结构进入压溃过程,载荷‑位移曲线趋于稳定,每种结构的SEA略有差异,平均压溃载荷基本稳定在10~12 kN之间。
(2)当结构位移行进至15 mm左右时,曲线会出现较小峰值,可能是由于等效薄弱环节与腹板之间的刚度不连续导致的。
(3)与其他波数模型相比,单波波纹腹板梁长度过短,其结构稳定性较差,当压溃到一定位移后(18 mm),承载能力发生明显下降,且压溃过程中载荷波动较为剧烈。两波波纹腹板梁峰值载荷出现得略晚(21 mm),但整个压溃过程载荷维持在一定范围,没有明显降低。
由此可知,当波纹腹板梁长度较短时,会导致结构的局部屈曲,从而改变腹板上的载荷分布状态,使得结构在压溃吸能过程中出现横向位移,进而影响其抗冲击性能,结构的SEA低于平均水平。因此,一般避免选用单波或两波的波纹腹板梁,而多选用三波、四波甚至更多波数的波纹腹板梁应用于抗坠撞的结构中。
(1)对复合材料波纹腹板梁进行了坠撞试验研究,得到了试验件的损伤破坏形貌,以及载荷和能量的历程曲线。结果表明,复合材料波纹腹板梁的压溃破坏过程主要分为两个阶段:初始局部破坏阶段和稳定压溃破坏阶段,波纹腹板梁断口形成于初始局部破坏阶段,并且其形貌主要依赖于薄弱环节的设置。
(2)基于ABAQUS/Explicit平台二次开发得到了模拟复合材料波纹腹板梁受冲击压溃过程的仿真分析模型,采用改进的Hashin损伤判定准则和Choi‑Chang准则综合判断单元失效,并结合Cohesive界面单元,可较为真实地反映所研究复合材料层合结构的各向异性和渐进损伤特性。
(3)建立了不同波数的复合材料波纹腹板梁的有限元模型。由数值计算结果可以得出,当腹板单波长保持不变时,腹板的波数(波纹腹板梁的长度)对抗冲击吸能的影响较大。当波纹腹板梁长度较大时(L≫H),结构才能够稳定有效地吸收能量。因此,在抗坠撞结构设计时,应综合考虑波纹腹板梁的结构参数。
参考文献
龚俊杰, 王鑫伟. 复合材料波纹梁吸能能力的数值模拟[J]. 航空学报, 2005, 26(3): 298‑302. [百度学术]
GONG Junjie, WANG Xinwei. Numerical simulations of energy absorption capability of composite components[J]. Acta Aeronautica et Astronautica Sinica, 2005, 26(3): 298‑302. [百度学术]
苏璇. 复合材料波纹梁吸能特性分析及优化设计[D]. 天津: 中国民航大学, 2020. [百度学术]
SU Xuan. Energy absorption characteristics and optimization design of composite semi‑circular corrugated beam[D]. Tianjin: Civil Aviation University of China, 2020. [百度学术]
SOKOLINSKY V S, INDERMUEHLE K C, HURTADO J A. Numerical simulation of the crushing process of a corrugated composite plate[J]. Composites Part A, 2011, 42(9): 1119‑1126. [百度学术]
HAMZA S S, AL EMRAN I, YUHAZRI M Y, et al. Hybridization effect on crashworthiness parameters of natural composite[J]. Frontiers in Materials, 2021. DOI:10.3389/fmats.2021.619245. [百度学术]
黄建城, 王鑫伟. 薄弱环节对复合材料吸能元件轴向压溃性能的影响综述[J]. 船舶力学, 2011, 15(8): 930‑939. [百度学术]
HUANG Jiancheng, WANG Xinwei. Overview on effect of trigger geometry on the axial crushing behavior of composite components for energy absorption[J]. Journal of Ship Mechanics, 2011, 15(8): 930‑939. [百度学术]
MCGREGOR C J, VAZIRI R, POURSARTIP A, et al. Simulation of progressive damage development in braided composite tubes under axial compression[J]. Composites Part A: Applied Science & Manufacturing, 2007, 38(11): 2247‑2259. [百度学术]
MA J, YAN Y. Quasi‑static and dynamic experiment investigations on the crashworthiness response of composite tubes[J]. Polymer Composites, 2013, 34(7): 1099‑1109. [百度学术]
龚俊杰. 复合材料结构耐撞性及可控薄弱环节研究[D]. 南京: 南京航空航天大学, 2007. [百度学术]
GONG Junjie. On the crashworthiness of composite structures and controllable trigger mechanism[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2007. [百度学术]
陶然. 复合材料波纹梁耐撞性研究[D]. 南京: 南京航空航天大学, 2015. [百度学术]
TAO Ran. Crashworthiness analysis of composite waved beams[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2015. [百度学术]
CHANG F K, CHANG K Y. A progressive damage model for laminated composites containing stress concentration[J]. Journal of Composite Materials, 1987, 19(9): 834‑855. [百度学术]
CAMANHO P P, MATTHEWS F L. A progressive damage model for mechanically fastened joints in composite laminates[J]. Journal of Composite Materials, 1999,33: 2248‑2280. [百度学术]
MCGREGOR C, VAZIRI R, POURSARTIP A. Simulation of progressive damage development in braided composite tubes undergoing axial crushing[J]. Composites Part A: Applied Science and Manufacturing, 2007, 38(11): 2247‑2259. [百度学术]