中低速压缩加载下不同截面构型复合材料薄壁结构吸能特性及失效分析

张欣玥; 惠旭龙; 葛宇静; 舒挽; 白春玉; 刘小川

doi:10.11883/bzycj-2021-0347

中低速压缩加载下不同截面构型复合材料薄壁结构吸能特性及失效分析

doi: 10.11883/bzycj-2021-0347

中国飞机强度研究所结构冲击动力学航空科技重点实验室，陕西西安 710065

基金项目: 民机专项科研（MJ-2017-F-15）

详细信息

作者简介:
张欣玥（1994－　），女，硕士，工程师，nwpuzhangxinyue@163.com

通讯作者:
刘小川（1983－　），男，博士，研究员，asri02@163.com

中图分类号: O347
计量
- 文章访问数: 751
- HTML全文浏览量: 333
- PDF下载量: 60
- 被引次数: 1
出版历程
- 收稿日期: 2021-08-17
- 录用日期: 2022-03-17
- 修回日期: 2022-01-24
- 网络出版日期: 2022-04-06
- 刊出日期: 2022-06-24

Energy absorption characteristics and failure analysis of composite thin-walled structures with different cross-sectional configurations under medium- and low-speed compression loading

Aviation Key Laboratory of Science and Technology on Structures Impact Dynamics, Aircraft Strength Research Institute of China, Xi’an 710065, Shaanxi, China

摘要

摘要: 为研究开剖面复合材料薄壁吸能结构的吸能特性，基于高速液压伺服试验系统，开展了开剖面复合材料薄壁结构轴向压缩试验，分析了截面构型、截面长宽比、触发模式及加载速度对其吸能特性的影响，揭示了其在压溃过程中的失效及吸能机理。研究结果表明，复合材料薄壁结构压溃过程中主要通过材料弯曲、分层、剪切破坏以及压溃区之间的摩擦吸能。截面构型对其吸能特性影响显著，其中，帽形及Ω形试件的平均压溃载荷较C形试件分别高出14.1%和14.6%，比吸能较C形试件分别高出14.3%和14.8%；截面长宽比对复合材料薄壁结构吸能特性的影响不如截面构型明显；触发模式主要影响吸能结构的初始压溃阶段，在降低峰值载荷方面，C形试件采用45°倒角触发效果更好，帽形试件采用15°尖顶触发效果更好；当加载速度从0.01 m/s提高到1 m/s时，C形、帽形及Ω形试件的平均压溃载荷分别下降了6.1%、10.9%和6.1%，比吸能分别下降了6.2%、11.0%和6.2%。
- 复合材料薄壁结构 /
- 吸能 /
- 截面构型 /
- 截面长宽比 /
- 触发模式 /
- 加载速度
Abstract: In order to study the energy absorption characteristics of open-section thin-walled composite structures, axial compression tests were carried out by using a high-speed hydraulic servo test system. The loading speed was set to 0.01, 0.1 and 1 m/s. A high-speed camera was used to record the deformation and failure of the test specimens. The effects of cross-section shape, section aspect ratio, trigger mechanism, and loading speed on the energy absorption characteristics of the composite structures are analyzed. The failure and energy absorption mechanism of the structure in the crushing process is revealed. The results show that the energy absorption is mainly attributed to material bending, delamination, shear failure and friction between crushing zones during the crushing process. The cross-section shape has a significant influence on its energy absorption capacity. The average crushing loads of the hat shaped and Ω- shaped specimens are 14.1% and 14.6% higher than that of the C-channel specimens, and their specific energy absorption (SEA) are 14.3% and 14.8% higher than that of C-channel specimens, respectively. The stress concentration of C-channel specimens leads to insufficient material damage, responsible to their lower energy absorption capacity. On the other hand, the section aspect ratio has less effect on the energy absorption capacity of composite thin-walled structures. The trigger mechanism mainly affects the initial crushing stage of the structures. For the C-channel specimens, 45° chamfer trigger is more effective in reducing the initial peak load; while for the hat shaped test piece, the 15° steeple trigger is better. When the loading speed was increased from 0.01 m/s to 1 m/s, the average crushing load of the C channel, hat shaped and Ω-shaped specimens were reduced by 6.1%, 10.9% and 6.1%, respectively; while the SAE were reduced by 6.2%, 11.0% and 6.2% respectively. The increase of loading speed leads to more debris flying out, which reduces the loading area and material utilization of the structure, and it reduces the friction energy absorption of the collapse zone, too.
- composite thin-walled structure /
- energy absorption /
- cross-section shape /
- section aspect ratio /
- trigger mechanism /
- loading speed

HTML全文

爆震波是激波和燃烧波的强烈耦合，涉及爆震燃烧的流动存在多个时间尺度，容易出现严重的刚性源项，使数值求解流动方程存在困难。

目前处理刚性源项问题的方法主要有分裂求解和耦合求解两类。采用显式耦合方法求解时，时间步长需要足够小以满足稳定性要求，否则可能导致计算不稳定甚至出现非物理解^[1]。Bussing等^[2-3]针对反应流中出现的多尺度问题，提出了点隐算法，将反应与流动的积分时间步长统一到伪时间步长上，消除了源项的刚性。Zhong^[4]推广了Rosenbrock^[5]从半隐出发求解刚性常微分方程组的方法，采用不同的方法分别求解方程的刚性项和非刚性项，增强了计算稳定性，提高了时间方向的精度，并在化学反应流的模拟中得到了成功应用。隐式处理刚性源项时，虽然很好地解决了刚性问题，但由于涉及矩阵求逆运算，随着反应组分增多，计算量急剧增大。

Leveque等^[6]利用时间分裂法和MacCormack预测校正方法求解了含刚性源项的双曲型方程，认为分裂算法在空间分辨率足够的情况下，计算效果比耦合算法稍好。Oran等^[7]对刚性源项方程的积分方法做了全面的论述。一步法虽然在求解刚性方程组时稳定性要求严格，但由于只需已知变量值，求解相对方便，并针对特征值为正且相当大的物理刚性问题时，使用一步法是比较有利的。Young等^[8]采用逼近法与预测校正思想来处理化学反应源项，对刚性和普通常微分方程组都具有较好的计算性能，且方便与时间步长估算方法^[9]配合使用，求解效率高，但容易出现质量分数不守恒的问题，需要在流场推进时不断修正。Mott等^[10-11]在逼近法的基础上发展了αQSS方法，通过迭代校正组分值来保证质量守恒，比逼近法更加稳定和精确，并针对碳氢燃料的燃烧和爆震模拟取得了令人满意的效果。

刘瑜^[12]对比了VODE(variable-coefficient ordinary differential equation solver)、梯形公式^[13]和αQSS方法处理化学反应刚性源项的能力，认为在计算精度基本相同的情况下，梯形公式和αQSS方法效率更高。刘世杰等^[14]利用梯形公式和αQSS方法对爆震波胞格进行了模拟，认为使用αQSS方法时计算的化学反应更完全，放热量更大。目前对刚性源项处理方法在化学反应流动模拟中的适用性，以及各方法之间的联系方面还缺少一定的认识。本文中在对各源项处理方法的稳定性分析基础上，分析各种方法之间的关系以及对化学反应源项的处理能力，并通过具体算例比较不同处理方法的计算效率。

1. 计算方法介绍

1.1 控制方程

激波诱导燃烧为激波和燃烧波耦合的结果，可忽略黏性的影响，因此使用有限速率化学反应的欧拉方程描述为:

$\frac{\partial }{{\partial t}}\int_\mathit{\Omega } {\mathit{\boldsymbol{W}}{\rm{d}}\mathit{\Omega }} + \oint_{\partial \mathit{\Omega }} {\mathit{\boldsymbol{F}}{\rm{d}}S} + \int_\mathit{\Omega } {{\mathit{\boldsymbol{Q}}_{\rm{s}}}{\rm{d}}\mathit{\Omega }} = \int_\mathit{\Omega } {\mathit{\boldsymbol{Q}}{\rm{d}}\mathit{\Omega }}$

(1)

式中：W=(ρY₁, …, ρY_N, ρu, ρv, ρE)^T，F=(ρY₁V, …, ρY_NV, ρuV+n_xp, ρvV+n_yp, ρHV)^T， ${\mathit{\boldsymbol{Q}}_s} = \frac{\delta }{y}$ (ρY₁v, …, ρY_Nv, ρuv, ρv², ρHv)^T， $\mathit{\boldsymbol{Q}} = {({{\dot \omega }_1}, \ldots ,{{\dot \omega }_N},0,0,0,0)^{\rm{T}}}$ 。其中：W为守恒变量，F为对流通量，Q为化学反应源项，Q_s为轴对称源项，Ω为网格体积，∂Ω表示包含网格的封闭曲面，δ=0, 1, 2时分别代表了二维平面流动、轴对称和球对称的控制方程，Y_i为组分i的质量分数，N为组分个数，总密度ρ可通过质量守恒得到，p为压力，u，v分别为x与y方向的速度，V为逆变速度，n_x和n_y分别表示界面法向向量在x和y方向的分量。假定各组分为热完全气体，组分的定压比热c_pi、定容比热c_Vi都只是温度的函数，单位质量的内能e_i、焓h_i、总能E、总焓H、混合物的绝热指数γ以及组分生成速率 ${{\dot \omega }_i}$ 的表达式可参见文献[7]。

1.2 数值方法

对式(1)控制方程采用有限体积离散，界面两侧变量值使用2阶MUSCL(monotonic upwind scheme for conversation laws)重构，采用AUSMPW+(advection upstream splitting method by new pressure-based weighted function)分裂格式进行通量求解，时间推进采用具有TVD(total variation diminishing)性质的2阶Runge-Kutta法，化学反应采用Jachimowski 8组分19步反应机理^[15]。在涉及化学反应的流动中，组分i的反应特征时间τ_i^(r)为源项雅各比矩阵对角元素 $\frac{{\partial {{\dot \omega }_i}}}{{\partial {\rho _i}}}$ 的倒数，通常为10^-11~10^-9 s，流动特征时间τ^(f)为网格尺度与当地声速的比值，一般为10^-9~10^-5 s。采用时间尺度定义的刚性系数R_s为最大特征时间与最小特征时间的比值，范围在1~10⁶。当网格尺度特别小时，τ^(f)与τ^(r)相当，然而在绝大多数计算中R_s很大，刚性问题严重，因此分别采用Strang分裂和耦合求解方法来处理化学反应出现的刚性源项问题，其中一步法、逼近法、αQSS方法为Strang分裂中化学源项的求解算子，而点隐则与流动耦合求解。

具有二阶时间精度的Strang分裂格式表示为：

${\mathit{\boldsymbol{W}}^{\left( {n + 1} \right)}} = {L_{\rm{r}}}\left( {\frac{{\Delta t}}{2}} \right){L_{\rm{f}}}\left( {\Delta t} \right){L_{\rm{r}}}\left( {\frac{{\Delta t}}{2}} \right){\mathit{\boldsymbol{W}}^{\left( n \right)}}$

(2)

式中：算符L_f为流动求解算子，L_r为化学反应的求解算子，Δt为流场推进的时间步长，上标n为求解时间步数。

一步法为常微分方程初值问题的经典求解方法，求解变量n+1时刻的值时，只需n时刻的信息。解耦后的化学源项方程如下：

$\frac{{{\rm{d}}\mathit{\boldsymbol{W}}}}{{{\rm{d}}t}} = F\left( {t,\mathit{\boldsymbol{W}}} \right)$

(3)

由于Strang分裂对算子求解精度的要求，选取2阶Runge-Kutta法作为化学反应求解算子。

在利用逼近法求解化学反应源项时，将组分i的净生成率写成如下形式：

$\frac{{{\rm{d}}{\rho _i}}}{{{\rm{d}}t}} = {q_i} - \frac{{{\rho _i}}}{{{\tau _i}}}\;\;\;\;\;\;1 \le i \le N$

(4)

式中：q_i为组分i在反应中的生成项，τ_i为组分i的消耗特征时间。则刚性方程组中组分i在n+1时刻的密度为:

${{\bar \rho }_i} = \frac{{\rho _i^{\left( n \right)}\left( {2\tau _i^{\left( n \right)} - \Delta t} \right) + 2\Delta tq_i^{\left( n \right)}\tau _i^{\left( n \right)}}}{{2\tau _i^{\left( n \right)} + \Delta t}}$

(5)

$\rho _i^{\left( {n + 1} \right)} = \frac{{\rho _i^{\left( n \right)}\left( {{{\bar \tau }_i} + \tau _i^{\left( n \right)} - \Delta t} \right) + 0.5\Delta t\left( {{{\bar q}_i} + q_i^{\left( n \right)}} \right)\left( {{{\bar \tau }_i} + \tau _i^{\left( n \right)}} \right)}}{{{{\bar \tau }_i} + \tau _i^{\left( n \right)} + \Delta t}}$

(6)

其中式(5)为预测步，式(6)为校正步；τ_i、q_i表示预测步的求解结果。

αQSS方法在逼近法的基础上进一步改进其迭代形式，在求解时将组分净生成率写为：

$\frac{{{\rm{d}}{\rho _i}}}{{{\rm{d}}t}} = {q_i} - {p_i}{\rho _i}\;\;\;\;\;\;\;\;1 \le i \le N$

(7)

式中：p_i为特征消耗时间的倒数，为了形式简便略去下标，如果q_i与p_i在时间段(t, t+Δt)之内为常数时，则可积分式(7)得到精确解：

${\rho ^{\left( {n + 1} \right)}} = {\rho ^{\left( n \right)}}{{\rm{e}}^{ - p\Delta t}} + \frac{q}{p}\left( {1 - {{\rm{e}}^{ - p\Delta t}}} \right)$

(8)

整理式(8)得到：

${\rho ^{\left( {n + 1} \right)}} = {\rho ^{\left( n \right)}} + \frac{{\Delta t\left( {q - p{\rho ^{\left( n \right)}}} \right)}}{{1 + \alpha p\Delta t}}$

(9)

$\alpha \left( {p\Delta t} \right) = \frac{{1 - \left( {1 - {{\rm{e}}^{ - p\Delta t}}} \right)/\left( {p\Delta t} \right)}}{{1 - {{\rm{e}}^{ - p\Delta t}}}}$

(10)

式中：pΔt→∞时，α→1，表示无穷快的反应；而pΔt→0时，α→1/2，表示无穷慢的反应。αQSS方法中利用预测校正的方法获得一定的精度并保证质量分数守恒，而在计算中需要不断对校正步进行迭代，直到满足迭代标准，预测和校正步的详细实施过程见文献[10-11]。

使用点隐方法时，将源项写成n+1时间步的值，通过线性化消除源项刚性的问题，实现方式为：

${\mathit{\boldsymbol{Q}}^{\left( {n + 1} \right)}} \approx {\mathit{\boldsymbol{Q}}^{\left( n \right)}} + \frac{{\partial \mathit{\boldsymbol{Q}}}}{{\partial \mathit{\boldsymbol{W}}}}\Delta \mathit{\boldsymbol{W}}$

(11)

在与流动耦合，使用具有TVD性质的二阶Runge-Kutta法推进时，其形式可写为:

$\begin{array}{l} {\mathit{\boldsymbol{W}}^{\left( {n + \frac{1}{2}} \right)}} = {\mathit{\boldsymbol{W}}^{\left( n \right)}} - {\left( {\frac{1}{{\Delta t}}\mathit{\boldsymbol{I}} - \frac{{\partial \mathit{\boldsymbol{Q}}}}{{\partial \mathit{\boldsymbol{W}}}}} \right)^{ - 1}}\left( {\frac{{{\mathit{\boldsymbol{R}}^{\left( n \right)}}}}{\mathit{\Omega }} - {\mathit{\boldsymbol{Q}}^{\left( n \right)}}} \right)\\ {\mathit{\boldsymbol{W}}^{\left( {n + 1} \right)}} = \frac{1}{2}{\mathit{\boldsymbol{W}}^{\left( n \right)}} + \frac{1}{2}{\mathit{\boldsymbol{W}}^{\left( n \right)}} + \frac{1}{2}{\mathit{\boldsymbol{W}}^{\left( {n + \frac{1}{2}} \right)}} - \\{\left( {\frac{2}{{\Delta t}}\mathit{\boldsymbol{I}} - \frac{{\partial \mathit{\boldsymbol{Q}}}}{{\partial \mathit{\boldsymbol{W}}}}} \right)^{ - 1}}\left( {\frac{{{\mathit{\boldsymbol{R}}^{\left( {n + \frac{1}{2}} \right)}}}}{\mathit{\Omega }} - {\mathit{\boldsymbol{Q}}^{\left( n \right)}}} \right) \end{array}$

(12)

式中：I为单位质量矩阵；R为空间方程离散后的余项。源项雅各比矩阵推导见文献[2, 16]。

2. 源项处理方法的理论分析

2.1 源项处理方法的稳定性分析

考虑模型方程：

$\frac{{\partial u}}{{\partial t}} + a\frac{{\partial u}}{{\partial x}} = s\left( u \right)\;\;\;\;\;\;\;\;a > 0$

(13)

式中：a为常数，s(u)为反应源项。

为了简便，采用一阶迎风差分格式离散式(13)的对流项，得到如下形式:

$\frac{{u_j^{\left( {n + 1} \right)} - u_j^{\left( n \right)}}}{{\Delta t}} + a\frac{{u_j^n - u_{j - 1}^{\left( n \right)}}}{{\Delta x}} = s\left( {u_j^{\left( n \right)}} \right)$

(14)

式(14)中的源项可根据处理方法写为不同的形式。

考虑一步法，将式(14)改写为:

$\frac{{u_j^{\left( {n + 1} \right)} - u_j^{\left( n \right)}}}{{\Delta t}} + a\frac{{u_j^{\left( n \right)} - u_{j - 1}^{\left( n \right)}}}{{\Delta x}} = - \frac{{u_j^{\left( n \right)}}}{{{\tau ^{\left( {\rm{c}} \right)}}}}$

(15)

引入单波扰动u_j⁽ⁿ⁾=G⁽ⁿ⁾e^ijθ，代入式(15)中，整理得到误差放大因子G的表达式为:

$G = 1 - \lambda \Delta t$

(16)

其中，λ的具体形式为:

$\lambda = \frac{{a\left( {1 + \cos \theta } \right)}}{{\Delta x}} = \frac{1}{{{\tau ^{\left( {\rm{c}} \right)}}}} - i\frac{{a\sin \theta }}{{\Delta x}}$

(17)

为满足求解稳定性，要求误差放大因子G的模小于或等于1，解出满足稳定要求的Δt，其形式为：

$\Delta t \le \frac{{2{\mathop{\rm Re}\nolimits} \left( \lambda \right)}}{{{{\left( {{\mathop{\rm Re}\nolimits} \left( \lambda \right)} \right)}^2} + {{\left( {{\mathop{\rm Im}\nolimits} \left( \lambda \right)} \right)}^2}}}$

(18)

从式(17)可以看出，当τ^(c)非常小时，其倒数远大于其他项的值，式(17)变为λ≈1/τ^(c)。因此式(18)可简化为：

$\Delta t \le 2{\tau ^{\left( {\rm{c}} \right)}}$

(19)

可见，在使用一步法求解刚性问题严重的化学反应流动时，化学源项积分的时间步长需要小于或等于两倍最小特征反应时间。

逼近法的分析类似，假设在Δt内q和τ为常数，则可将式(14)改写为：

$\frac{{u_j^{\left( {n + 1} \right)} - u_j^{\left( n \right)}}}{{\Delta t}} + a\frac{{u_j^{\left( n \right)} - u_{j - 1}^{\left( n \right)}}}{{\Delta x}} = {q_j} - \frac{{u_j^{\left( {n + 1} \right)} - u_j^{\left( n \right)}}}{{2\tau _j^{\left( n \right)}}}$

(20)

该形式与梯形公式相同，其稳定性分析见文献[1]，结果表明，逼近法只要时间步长满足对流项稳定条件即可。

点隐方法的分析与一步法类似，只需将式(15)中源项的变量换为n+1时刻的值，解出误差放大因子的模为:

$\left| G \right| = \left| {\frac{{1 - \left( {1 - {{\rm{e}}^{{\rm{i}}\theta }}} \right)a\Delta t/\Delta x}}{{1 + \Delta t/{\tau ^{\left( {\rm{c}} \right)}}}}} \right|$

(21)

式中：aΔt/Δx为CFL(Courant Friedrichs Lewy)数，而Δt/τ^(c)为刚性系数，可知只要满足对流项稳定条件，即CFL数小于1，误差放大因子的模是恒小于1的。

αQSS方法的稳定性分析可见文献[10]，在满足对流项稳定的条件下是绝对稳定的。

稳定性分析的结果表明:隐式方法在求解数学意义上的刚性问题时具有卓越的性能；一步法为了满足求解的稳定性，需要很小的时间步长；逼近法和αQSS方法中的特征消耗时间τ_i满足在求解时间步长内为常数的情况下，也具有很好的求解性能。

2.2 对化学反应流动模拟的适用性分析

化学反应流动存在多个时间尺度，但可能不完全符合刚性问题的数学定义^[17]，因为化学反应源项的雅各比矩阵特征值实部可能不全为负值^[7]。此时，特征值为正且相当大的情况表现为物理上的刚性问题，求解时必须保证时间步长足够小以捕捉到物理特征的变化。如果源项雅各比矩阵特征值实部全为负值，则对应的物理特征表现为组分快速变化到反应平衡状态。

基于以上分析，以点隐为首的隐式方法在求解数学意义上的刚性问题时具有优秀的性能，但针对物理意义上的刚性问题时，由于点隐对化学反应源项积分没有时间步长限制，可能会由于时间步长大而导致分辨不出物理特征变化的细节；另外，流场求解时网格数较多，使用隐式方法时，在每个网格单元都要进行矩阵的求逆运算，运算量随变量个数二次方增加^[1]，导致流场推进速度很慢。一步法受稳定性要求的严格限制，在刚性很严重的时候，时间步长取值造成的计算量需求将难以接受，但在求解物理和数学上都为刚性的问题时还是具有一定的优势，至少能够保证捕捉到变化细节所需要的时间步长。逼近法和αQSS方法在求解时简化了物理问题，将组分反应源项考虑为解耦形式，从而降低了求解的难度，但造成了在计算时容易出现质量不守恒的问题，使用逼近法时需要不断修正，而αQSS方法则利用组分迭代校正到一定精度来满足质量守恒。尽管如此，αQSS在积分化学反应源项时还是具有良好的计算性能，且对于化学反应的变化适应性很好，具体分析如下。

将组分i的化学反应净生成项写为式(7)形式，并将式(5)写成式(9)形式，得到：

${\rho ^{\left( {n + 1} \right)}} = {\rho ^{\left( n \right)}} + \frac{{\Delta t\left( {q - p{\rho ^{\left( n \right)}}} \right)}}{{1 + \frac{1}{2} \cdot p\Delta t}}$

(22)

逼近法对应αQSS方法中pΔt→0的情况，而梯形公式与逼近法本质上相同，只是缺少校正步。

对式(7)分别采用隐式和显式进行数值求解，得到如下两种形式:

$\frac{{{\rho ^{\left( {n + 1} \right)}} - {\rho ^{\left( n \right)}}}}{{\Delta t}} = q - p{\rho ^{\left( {n + 1} \right)}}$

(23)

$\frac{{{\rho ^{\left( {n + 1} \right)}} - {\rho ^{\left( n \right)}}}}{{\Delta t}} = q - p{\rho ^{\left( n \right)}}$

(24)

分别将式(23)、(24)改写成式(9)形式，得到：

${\rho ^{\left( {n + 1} \right)}} = {\rho ^{\left( n \right)}} + \frac{{\Delta t\left( {q - p{\rho ^{\left( n \right)}}} \right)}}{{1 + 1 \cdot p\Delta t}}$

(25)

${\rho ^{\left( {n + 1} \right)}} = {\rho ^{\left( n \right)}} + \frac{{\Delta t\left( {q - p{\rho ^{\left( n \right)}}} \right)}}{{1 + 0 \cdot p\Delta t}}$

(26)

由αQSS中参数α的定义，式(25)为解耦的隐式方法，对应pΔt→∞的情况，适合快速衰减的反应，与点隐方法类似。而式(26)则为一步法，对应pΔt→-∞的情况，表明组分消耗特征时间为负值，而αQSS方法不存在这种情况。因此，由式(22)~(26)可以看出，如果组分反应可表示为解耦的形式，如式(7)，则αQSS方法几乎涵盖了其他几种方法。但是，αQSS方法中的参数α变化只能表征组分衰减到平衡状态快慢的速度，即对应特征值小于零的情况，在生成项q远大于消耗项pρ时，只依赖参数α变化并不能捕捉到快速增长的化学反应特征，此时αQSS方法通过对比生成项和消耗项的大小，自动调整时间步长以适应化学反应的变化^[11]，其性质与一步法相同。从以上分析可以得到，对应物理模态快速衰减的刚性问题时，αQSS方法通过参数α的变化来适应不同的衰减速度；而对于物理模态快速增长的刚性问题时，αQSS方法通过改变积分时间步长来适应增长速度，因此对化学反应变化的适应性强。

3. 数值计算

Lehr^[18]进行了一系列弹道靶实验，将球形弹丸以接近爆震的速度射入一定条件下的氢气和空气预混气体中，观察到弹丸头部不稳定的燃烧现象，并拍摄了不同入射速度下弹丸头部振荡燃烧现象的阴影图，图 1(a)、(b)分别对应来流马赫数为4.48和4.79的工况。McVey等^[19]提出了解释激波诱导振荡燃烧的机理，其后Choi^[20]对该现象做了详尽的数值模拟，认为激波诱导周期振荡燃烧现象是验证非平衡化学流计算程序的典型算例。本文中使用该算例，从计算结果、计算效率等方面比较了源项处理方法的性能。

图 1 激波诱导的周期振荡燃烧实验阴影图^[18]

Figure 1. Experimental shadowgraph of periodically oscillating shock-induced combustion^[18]

下载: 全尺寸图片幻灯片

3.1 计算结果可靠性分析

计算模型和网格选取与文献[21]相同，取球头前方2.5 mm，上方10 mm，网格点数为251×201。计算工况与Lehr^[18]实验条件一致，来流温度293 K，压力42 665 Pa，预混气体来流速度分别为1 804 m/s(Ma=4.48)和1 931 m/s (Ma=4.79)，该实验条件下对应的爆震马赫数约为4.845，CFL数取为0.85。通过记录球头轴线驻点处压力随时间变化的数据，对其进行快速傅里叶变换得到振荡频率，对比计算振荡频率和实验所得振荡频率，验证化学反应模块的正确性。为了避免由于取值区间小而造成的振荡频率计算的波动，用于计算主频的压力时间曲线选取至少100个振荡周期。

为验证源项处理方法计算结果的正确性，对来流马赫数分别为4.48和4.79的工况进行了计算，其轴线驻点压力振荡频率见表 1。从表 1可以看出，来流马赫数分别为4.48和4.79时，不同源项处理方法得到的振荡频率与实验值相对误差都较小。在使用251×201的网格尺度且来流马赫数4.79时的计算条件下，计算得到的组分最小反应特征时间约为6×10^-11 s，流动特征时间为3.5×10^-9 s，刚性系数R_s约为58，从稳定性分析可知，采用的化学反应时间步长应小于等于2倍特征反应时间，但文献[2]表明，在实际计算中化学反应推进步长可取为10倍最小反应特征时间。表 1所列的一步法中，化学反应时间步长为6倍最小时间尺度，计算得到的振荡频率与实验值吻合，但从稳定性方面考虑，计算过程中可能存在一定的振荡。在Jachimowski化学反应模型中，组分H₂O₂的化学反应特征时间最小，因此出现振荡的可能性最大，图 2~3所示为对应表 1中Ma=4.79时，不同源项处理方法计算得到的球头轴线上组分变化曲线，其中图 2为不同积分时间步长的一步法计算结果，图 3(a)、(b)、(c)分别为αQSS、点隐和逼近法的计算结果。由于各组分之间的质量分数相差较大，为了对比的直观性，对图中显示的各组分质量分数进行归一化处理，归一化之后，H₂的归一化质量分数为0.03，H₂O的归一化质量分数采用0.25，H₂O₂的归一化质量分数为0.000 5，而N₂的归一化质量分数为0.76。

表 1 实验数据与理论模型结果对比

Table 1. Comparison of experimental and theoretical data

方法	f/Hz
方法	Ma=4.48	Ma=4.79
实验	425.0	712.0
αQSS	424.0	721.5
点隐法	426.1	725.7
一步法	424.1	721.5
逼近法	425.2	721.5

下载: 导出CSV

| 显示表格

图 2 一步法处理组分质量分数随轴线分布曲线

Figure 2. Mass fraction distribution of species along the axis by one step method

下载: 全尺寸图片幻灯片

3(a) αQSS方法处理组分质量分数随轴线分布曲线

3(a). Mass fraction distribution of species along the axis by αQSS method

下载: 全尺寸图片幻灯片

3(b) 点隐法处理组分质量分数随轴线分布曲线

3(b). Mass fraction distribution of species along the axis by point implicit method

下载: 全尺寸图片幻灯片

3(c) 逼近法处理组分质量分数随轴线分布曲线

3(c). Mass fraction distribution of species along the axis by asymptotic method

下载: 全尺寸图片幻灯片

图 2(a)表示使用一步法时积分步长为6倍最小化学反应特征时间的计算结果，与图 2(b)使用两倍化学反应特征时间的计算结果对比可发现，在燃烧阵面之后，H₂O₂变化曲线存在非物理振荡现象，与稳定性分析的预测结果一致，即积分时间步长超过稳定性要求的最大时间步长后，计算结果将出现不稳定或者发散现象，说明在使用一步法积分化学反应刚性源项时，要想获得准确的组分质量分数变化曲线，时间步长必须满足稳定性要求。积分时间步长取值满足稳定要求时，几种源项处理方法的结果相同，如图 3(a)~(c)所示，因此，在计算化学反应流动时，逼近法、αQSS和点隐在稳定性方面要优于一步法。

观察图 2(a)，发现振荡并没有随时间和空间变化而加剧，这是因为化学反应具有强非线性与随时间变化的特点，源项的雅各比矩阵特征值将在反应过程中不断变化，使得引入的小扰动在增长模态和衰减模态之间反复转变，因此扰动不会进一步增强。另外，从表 1可知，在使用一步法计算时，积分时间步长不满足稳定性要求时，也能获得正确的燃烧振荡频率。这是由于振荡频率主要由化学反应放热量决定的，H₂O₂组分质量分数很小，对于整个化学反应放热量的影响很微弱，因此其质量分数的小波动对于振荡频率几乎没有影响。此外，图 3中组分变化曲线表明这几种源项处理方法对本算例中化学反应特征的捕捉效果基本相同，其原因在于化学反应是一个动态的非线性系统，某组分在这一时刻呈快速的增长模态，在下一时刻变为快速衰减的模态或者趋于平衡；并且本算例的刚性系数不大，化学反应特征在单个流动时间步内的变化也很小，因此各个源项处理方法对于该算例的捕捉效果基本相同。

3.2 计算效率分析

在计算硬件条件相同且CFL数为0.85的情况下，以来流马赫数4.79且完全发展的振荡燃烧流场作为初始计算流场，比较几种源项处理方法的计算效率。表 2所示为流场推进物理时间473 μs时，不同源项处理方法所用的CPU时间，其中一步法的化学推进时间步长为6倍最小反应特征时间。

表 2 不同源项处理方法消耗CPU时间

Table 2. The CPU time in different stiff source term methods

方法	CPU time/s
αQSS	77 590
点隐法	130 007
一步法	154 774
逼近法	62 372

下载: 导出CSV

| 显示表格

从表 2中可以看到，点隐与一步法消耗的CPU时间最多，采用源项完全求逆的点隐所需时间约为αQSS的2倍，点隐方法需要计算源项的雅各比矩阵，并涉及复杂的矩阵求逆运算，求解时间随组分的增加呈指数增长，在使用基元反应模型求解爆震燃烧问题时，计算量很大；而使用一步法时，时间步长需取足够小以满足稳定性要求，极大降低了计算效率；αQSS计算时间与组分求解迭代时指定的误差ε有关，ε越小，所需计算时间越长；逼近法相对αQSS来说，缺少迭代过程，因此所用时间最少。由此可见，在计算结果精度基本相同的情况下，逼近法和αQSS的计算效率较其他两种处理方法高。

4. 结论

本文中从稳定性分析和数值计算方面对比一步法、逼近法、αQSS和点隐方法计算化学反应刚性源项问题的性能，得到结论如下：(1)一步法在积分刚性源项时，积分步长需小于或等于2倍最小时间尺度，而其他3种源项计算方法对时间步长取值没有影响；(2)一步法、逼近法实质上是αQSS的特例，αQSS方法可根据化学反应特征自动调整系数α和时间步长，适用范围广；(3)一步法在积分化学反应刚性源项时，可适当放宽稳定性要求的时间步长限制；(4)在计算结果相近以及和精度基本相同的前提下，逼近法和αQSS在求解刚性源项时具有较高的效率，并且αQSS能够适应化学反应的变化，是采用分裂法求解刚性源项时综合性能最好的方法。

图 1 试件横截面尺寸（单位：mm）

Figure 1. Cross-sectional dimensions of specimens (unit: mm)

下载: 全尺寸图片幻灯片

图 2 试件照片

Figure 2. Photos of the specimens

下载: 全尺寸图片幻灯片

图 3 高速液压伺服试验机

Figure 3. High speed hydraulic servo testing machine

下载: 全尺寸图片幻灯片

图 4 试件夹持方式

Figure 4. Clamping method of the specimens

下载: 全尺寸图片幻灯片

图 5 典型复合材料吸能试件渐进压溃载荷-位移曲线

Figure 5. Typical progressive crushing load-displacement curve of composite energy-absorbing specimens

下载: 全尺寸图片幻灯片

图 6 不同截面构型试件压溃过程中的典型载荷-位移曲线

Figure 6. Typical force-displacement curves of the specimenswith different cross-section shapes

下载: 全尺寸图片幻灯片

图 7 不同截面构型试件吸能特性对比

Figure 7. Comparison of energy-absorption characteristics of the specimens with different cross-section shapes

下载: 全尺寸图片幻灯片

图 8 不同截面构型试件加载过程中的破坏情况（加载速度1 m/s）

Figure 8. Failure modes of the specimens with different cross-section shapes during loading (loading speed: 1 m/s)

下载: 全尺寸图片幻灯片

图 9 加载速度为1 m/s时不同截面构型试件的破坏形貌

Figure 9. Failure morphology of the specimens with different cross-section configurations under loading speed of 1 m/s

下载: 全尺寸图片幻灯片

图 10 不同截面长宽比试件压溃过程中的典型载荷-位移曲线

Figure 10. Typical force-displacement curves of the specimens with different section aspect ratios

下载: 全尺寸图片幻灯片

图 11 不同长宽比试件吸能特性对比

Figure 11. Comparison of energy-absorption characteristics of the specimens with different section aspect ratios

下载: 全尺寸图片幻灯片

图 12 不同截面构型试件破坏形貌

Figure 12. Failure morphology of the specimens with different cross-section shapes

下载: 全尺寸图片幻灯片

图 13 不同触发模式试件典型载荷-位移曲线

Figure 13. Typical force-displacement curves of the specimens with different trigger methods

下载: 全尺寸图片幻灯片

图 14 不同触发模式试件吸能特性对比

Figure 14. Comparison of energy-absorption characteristics of the specimens with different trigger methods

下载: 全尺寸图片幻灯片

图 15 尖顶触发试件加载过程中破坏情况

Figure 15. Failure modes of the specimens with steeple trigger method

下载: 全尺寸图片幻灯片

图 16 不同加载速度下的典型载荷-位移曲线

Figure 16. Typical force-displacement curves of the specimens with different loading speeds

下载: 全尺寸图片幻灯片

图 17 不同加载速度的试件吸能特性对比

Figure 17. Comparison of energy-absorption characteristics of the specimens with different loading speeds

下载: 全尺寸图片幻灯片

图 18 不同构型试件比吸能随加载速度变化情况

Figure 18. Variation of specific energy absorption of the specimensof different section shapes with loading speed

下载: 全尺寸图片幻灯片

图 19 不同截面构型试件加载过程中破坏情况（加载速度0.01 m/s）

Figure 19. Failure modes of the specimens of different cross-section shapes during loading (loading speed: 0.01 m/s)

下载: 全尺寸图片幻灯片

图 20 加载速度为0.01 m/s时不同截面构型试件破坏形貌

Figure 20. Morphology of the specimens of different section shapes under loading speed of 0.01 m/s

下载: 全尺寸图片幻灯片

表 1 复合材料薄壁吸能结构压溃试验

Table 1. Composite thin-walled structures compression test matrix

试件构型	触发方式	加载速度/（m∙s⁻¹）	压缩行程内质量/g	试验组数
C1型	倒角触发	1	14.680	3
帽形	倒角触发	1	14.678	3
Ω形	倒角触发	1	14.677	3
C2型	倒角触发	1	14.680	3
C3型	倒角触发	1	14.680	3
C1型	尖顶触发	1	13.507	3
帽形	尖顶触发	1	14.029	3
C1型	倒角触发	0.01	14.680	3
C1型	倒角触发	0.1	14.680	3
帽形	倒角触发	0.01	14.678	3
帽形	倒角触发	0.1	14.678	3
Ω形	倒角触发	0.01	14.677	3
Ω形	倒角触发	0.1	14.677	3

下载: 导出CSV

参考文献(36)

[1]	REN Y R, JIANG H Y, JI W Y, et al. Improvement of progressive damage model to predicting crashworthy composite corrugated plate [J]. Applied Composite Materials, 2017, 25(1): 45–66. DOI: 10.1007/s10443-017-9610-z.
[2]	刘小川, 周苏枫, 马君峰, 等. 民机客舱下部吸能结构分析与试验相关性研究 [J]. 航空学报, 2012, 33(12): 2202–2210. LIU X C, ZHOU S F, MA J F, et al. Correlation study of crash analysis and test of civil airplane sub-cabin energy absorption structure [J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(12): 2202–2210.
[3]	冯振宇, 解江, 李恒晖, 等. 大飞机货舱地板下部结构有限元建模与适坠性分析 [J]. 航空学报, 2019, 40(2): 522394. DOI: 10.7527/S1000-6893.2018.22394. FENG Z Y, XIE J, LI H H, et al. Finite element modeling and crashworthiness analysis of large aeroplane sub-cargo structure [J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(2): 522394. DOI: 10.7527/S1000-6893.2018.22394.
[4]	冯振宇, 程坤, 赵一帆, 等. 运输类飞机典型货舱地板下部结构冲击吸能特性 [J]. 航空学报, 2019, 40(9): 222907. DOI: 10.7527/S1000-6893.2019.22907. FENG Z Y, CHENG K, ZHAO Y F, et al. Energy-absorbing charateristics of a typical sub-cargo fuselage section of a transport category aircraft [J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(9): 222907. DOI: 10.7527/S1000-6893.2019.22907.
[5]	李松晏, 郑志军, 虞吉林. 高速列车吸能结构设计和耐撞性分析 [J]. 爆炸与冲击, 2015, 35: 164–170. DOI: 10.11883/1001-1455(2015)02-0164-07. LI S Y, ZHENG Z J, YU J L. Energy-absorbing structure design and crashworthiness analysis of high-speed trains [J]. Explosion and Shock Waves, 2015, 35: 164–170. DOI: 10.11883/1001-1455(2015)02-0164-07.
[6]	SUBBARAMAIAH R, PRUSTY B G, PEARCE G M K, et al. Crashworthy response of fibre metal laminate top hat structures [J]. Composite Structures, 2017, 160: 773–781. DOI: 10.1016/j.compstruct.2016.10.112.
[7]	JIANG H Y, REN Y R, GAO B H, et al. Design of novel plug-type triggers for composite square tubes: enhancement of energy-absorption capacity and inducing failure mechanisms [J]. International Journal of Mechanical Sciences, 2017, 131: 113–136. DOI: 10.1016/j.ijmecsci.2017.06.050.
[8]	REN Y R, ZHANG H Y, XIANG J W. A novel aircraft energy absorption strut system with corrugated composite plate to improve crashworthiness [J]. International Journal of Crashworthiness, 2018, 23(1): 1–10. DOI: 10.1080/13588265.2017.1301082.
[9]	冯振宇, 周坤, 裴惠, 等. 复合材料薄壁方管准静态轴向压缩失效机理及吸能特性 [J]. 高分子材料科学与工程, 2019, 35(8): 94–104. DOI: 10.16865/j.cnki.1000-7555.2019.0220. FENG Z Y, ZHOU K, PEI H, et al. Failure mechanism and energy-absorption characteristics of composite thin-walled square tube under quasi-static axial compression load [J]. Polymer Materials Science and Engineering, 2019, 35(8): 94–104. DOI: 10.16865/j.cnki.1000-7555.2019.0220.
[10]	解江, 马骢瑶, 霍雨佳, 等. 纤维铺层角度对复合材料薄壁圆管轴向压溃吸能特性影响研究 [J]. 振动与冲击, 2018, 37(20): 205–211. DOI: 10.13465/j.cnki.jvs.2018.20.031. XIE J, MA C Y, HUO Y J, et al. Effect of ply orientations on energy-absorbing characteristics of composite thin-walled circular tubes under axial compression [J]. Journal of Vibration and Shock, 2018, 37(20): 205–211. DOI: 10.13465/j.cnki.jvs.2018.20.031.
[11]	JIANG H Y, REN Y R, GAO B H. Research on the progressive damage model and trigger geometry of composite waved beam to improve crashworthiness [J]. Thin-Walled Structures, 2017, 119: 531–543. DOI: 10.1016/j.tws.2017.07.004.
[12]	KIM J S, YOON H J, SHIN K B. A study on crushing behaviors of composite circular tubes with different reinforcing fibers [J]. International Journal of Impact Engineering, 2011, 38(4): 198–207. DOI: 10.1016/j.ijimpeng.2010.11.007.
[13]	KAKOGIANNIS D, YUEN S C K, PALANIVELU S, et al. Response of pultruded composite tubes subjected to dynamic and impulsive axial loading [J]. Composites Part B: Engineering, 2013, 55: 537–547. DOI: 10.1016/j.compositesb.2013.07.022.
[14]	黄建城, 王鑫伟, 卞航. SMA薄弱环节对复合材料圆管耐撞性影响的试验研究 [J]. 工程力学, 2011, 28(10): 222–227. HUANG J C, WANG X W, BIAN H. Effect of SMA trigger on the crashworthiness of composite tubes [J]. Engineering Mechanics, 2011, 28(10): 222–227.
[15]	王振, 宋凯, 朱国华, 等. 单向碳纤维复合材料锥管轴向吸能特性研究 [J]. 振动与冲击, 2018, 37(7): 172–178. WANG Z, SONG K, ZHU G H, et al. Axial energy absorption characteristics of unidirectional carbon-fiber composite cone tubes [J]. Journal of Vibration and Shock, 2018, 37(7): 172–178.
[16]	邓亚斌, 任毅如, 蒋宏勇. 复合材料吸能圆管在半圆凹槽触发机制下的斜向压溃失效行为 [J]. 复合材料学报, 2022, 39(4): 1796–1804. DENG Y B, REN Y R, JIANG H Y. Oblique crushing failure behaviors of composite energy-absorbing circu-lar tube under the semi-circular cavity triggering mechanism [J]. Acta Materiae Compositae Sinica, 2022, 39(4): 1796–1804.
[17]	谭丽辉, 徐涛, 崔晓梅, 等. 带有圆弧形凹槽金属薄壁圆管抗撞性优化设计 [J]. 爆炸与冲击, 2014, 34(5): 547–553. DOI: 10.11883/1001-1455(2014)05-0547-07. TAN L H, XU T, CUI X M, et al. Design optimization for crashworthiness of metal thin-walled cylinders with circular arc indentations [J]. Explosion and Shock Waves, 2014, 34(5): 547–553. DOI: 10.11883/1001-1455(2014)05-0547-07.
[18]	殷之平, 李玉龙, 黄其青. 含诱导缺陷薄壁圆管耐撞性优化设计 [J]. 爆炸与冲击, 2011, 31(4): 418–422. DOI: 10.11883/1001-1455(2011)04-0418-05. YIN Z P, LI Y L, HUANG Q Q. Optimal crashworthiness design of thin-walled circular tubes with triggering holes [J]. Explosion and Shock Waves, 2011, 31(4): 418–422. DOI: 10.11883/1001-1455(2011)04-0418-05.
[19]	黄建城. 含薄弱环节复合材料圆管轴向吸能特性研究 [D]. 南京: 南京航空航天大学, 2011: 46–96. HUANG J C. On the axial energy absorption behaviour of composite tubes with crush triggers [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2011: 46–96.
[20]	DEEPAK S. Crashworthy design and analysis of aircraft structures [D]. Philadelphia: Drexel University, 2013: 87–227.
[21]	蒋宏勇. 复合材料薄壁结构的损伤耦合破坏模型及其吸能机理的研究 [D]. 长沙: 湖南大学, 2017: 33–36. JIANG H Y. Research on the damage coupling destruction model and energy-absorbing mechanism of composite thin-walled structure [D]. Changsha: Hunan University, 2017: 33–36.
[22]	ZHAO X, ZHU G H, ZHOU C Y, et al. Crashworthiness analysis and design of composite tapered tubes under multiple load cases [J]. Composite Structures, 2019, 222: 110920. DOI: 10.1016/j.compstruct.2019.110920.
[23]	郝文乾, 卢进帅, 黄睿, 等. 轴向冲击载荷下薄壁折纹管的屈曲模态与吸能 [J]. 爆炸与冲击, 2015, 35(3): 380–385. DOI: 10.11883/1001-1455-(2015)03-0380-06. HAO W Q, LU J S, HUANG R, et al. Buckling and energy absorption properties of thin-walled corrugated tubes under axial impacting [J]. Explosion and Shock Waves, 2015, 35(3): 380–385. DOI: 10.11883/1001-1455-(2015)03-0380-06.
[24]	解江, 张雪晗, 宋山山, 等. CFRP薄壁C型柱轴向压缩破坏机制及吸能特性 [J]. 复合材料学报, 2018, 35(12): 3261–3270. DOI: 10.13801/j.cnki.fhclxb.20180319.002. XIE J, ZHANG X H, SONG S S, et al. Failure mechanism and energy-absorbing characteristics of CFRP thin-walled C-channels subject to axial compression [J]. Acta Materiae Compositae Sinica, 2018, 35(12): 3261–3270. DOI: 10.13801/j.cnki.fhclxb.20180319.002.
[25]	JIANG H Y, REN Y R. Crashworthiness and failure analysis of steeple-triggered hat-shaped composite structure under the axial and oblique crushing load [J]. Composite Structures, 2019, 229: 111375. DOI: 10.1016/j.compstruct.2019.111375.
[26]	JOOSTEN M W, DUTTON S, KELLY D, et al. Experimental and numerical investigation of the crushing response of an open section composite energy absorbing element [J]. Composite Structures, 2011, 93(2): 682–689. DOI: 10.1016/j.compstruct.2010.08.011.
[27]	RICCIO A, RAIMONDO A, CAPRIO F D, et al. Experimental and numerical investigation on the crashworthiness of a composite fuselage sub-floor support system [J]. Composites Part B: Engineering, 2018, 150: 93–103. DOI: 10.1016/j.compositesb.2018.05.044.
[28]	RICCIO A, SAPUTO S, SELLITTO A, et al. On the crashworthiness behaviour of a composite fuselage sub-floor component [J]. Composite Structures, 2020, 234: 111662. DOI: 10.1016/j.compstruct.2019.111662.
[29]	JACKSON A, DUTTON S, GUNNION A J, et al. Investigation into laminate design of open carbon-fibre/epoxy sections by quasi-static and dynamic crushing [J]. Composite Structures, 2011, 93(10): 2646–2654. DOI: 10.1016/j.compstruct.2011.04.032.
[30]	解江, 宋山山, 宋东方, 等. 复合材料C型柱轴压失效分析的层合壳建模方法 [J]. 航空学报, 2019, 40(2): 522395. DOI: 10.7527/S1000-6893.2018.22395. XIE J, SONG S S, SONG D F, et al. Stacked shell modeling menthod for failure analysis of composite C-channel subject to axial cmpression [J]. Acta aeronautica et Astronautica Sinica, 2019, 40(2): 522395. DOI: 10.7527/S1000-6893.2018.22395.
[31]	汪洋, 吴志斌, 刘富. 复合材料货舱地板立柱压溃响应试验 [J]. 复合材料学报, 2020, 37(9): 2200–2206. DOI: 10.13801/j.cnki.fhclxb.20200111.001. WANG Y, WU Z B, LIU F. Crush experiment of composite cargo floor stanchions [J]. Acta Materiae Compositae Sinica, 2020, 37(9): 2200–2206. DOI: 10.13801/j.cnki.fhclxb.20200111.001.
[32]	FERABOLI P. Development of a corrugated test specimen for composite materials energy absorption [J]. Journal of Composite Materials, 2008, 42(3): 229–56. DOI: 10.1177/0021998307086202.
[33]	WAIMER M, KOHLGRÜBER D, HACHENBERG D, et al. Experimental study of CFRP components subjected to dynamic crash loads [J]. Composite Structures, 2013, 105: 288–299. DOI: 10.1016/j.compstruct.2013.05.030.
[34]	PATEL S, VUSA V R, SOARES C G. Crashworthiness analysis of polymer composites under axial and oblique impact loading [J]. International Journal of Mechanical Sciences, 2019, 156(1): 221–234. DOI: 10.1016/j.ijmecsci.2019.03.038.
[35]	WAIMER M, SIEMANN M H, FESER T. Simulation of CFRP components subjected to dynamic crash loads [J]. International Journal of Impact Engineering, 2017, 101(1): 115–131. DOI: 10.1016/j.ijimpeng.2016.11.011.
[36]	蒋宏勇, 任毅如, 袁秀良, 等. 基于非线性渐进损伤模型的复合材料波纹梁耐撞性能研究 [J]. 航空学报, 2017, 38(6): 220717. DOI: 10.7527/S1000-6893.2016.220717. JIANG H Y, REN Y R, YUAN X L, et al. Crashworthiness of composite corrugated beam based on nonlinear progressive damage model [J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(6): 220717. DOI: 10.7527/S1000-6893.2016.220717.

施引文献

期刊类型引用(0)
其他类型引用(1)

资源附件(0)

访问统计

图(20) / 表(1)

计量

文章访问数: 751
HTML全文浏览量: 333
PDF下载量: 60
被引次数: 1

1. 计算方法介绍
1.1 控制方程
1.2 数值方法
2. 源项处理方法的理论分析
2.1 源项处理方法的稳定性分析
2.2 对化学反应流动模拟的适用性分析
3. 数值计算
3.1 计算结果可靠性分析
3.2 计算效率分析
4. 结论

中低速压缩加载下不同截面构型复合材料薄壁结构吸能特性及失效分析

doi: 10.11883/bzycj-2021-0347

作者简介: 张欣玥（1994－ ），女，硕士，工程师，nwpuzhangxinyue@163.com

通讯作者: 刘小川（1983－ ），男，博士，研究员，asri02@163.com

计量

出版历程

Energy absorption characteristics and failure analysis of composite thin-walled structures with different cross-sectional configurations under medium- and low-speed compression loading

1. 计算方法介绍

1.1 控制方程

1.2 数值方法

2. 源项处理方法的理论分析

2.1 源项处理方法的稳定性分析

2.2 对化学反应流动模拟的适用性分析

3. 数值计算

3.1 计算结果可靠性分析

3.2 计算效率分析

4. 结论

期刊类型引用(0)

其他类型引用(1)

计量

出版历程

目录

1. 计算方法介绍

1.1 控制方程

1.2 数值方法

2. 源项处理方法的理论分析

2.1 源项处理方法的稳定性分析

2.2 对化学反应流动模拟的适用性分析

3. 数值计算

3.1 计算结果可靠性分析

3.2 计算效率分析

4. 结论

作者简介:
张欣玥（1994－　），女，硕士，工程师，nwpuzhangxinyue@163.com

通讯作者:
刘小川（1983－　），男，博士，研究员，asri02@163.com