地震作用下纵筋屈曲的研究进展

李晰; 杨昊

doi:10.3969/j.issn.0258-2724.20220549

地震作用下纵筋屈曲的研究进展

doi: 10.3969/j.issn.0258-2724.20220549

李晰^{1, 2,},
杨昊¹

1.
西南交通大学土木工程学院，四川成都 610031
2.
青岛理工大学土木工程学院，山东青岛 266033

基金项目: 国家自然科学基金项目（51508473）；山东省自然科学基金项目（ZR2021ME094）

详细信息

作者简介:
李晰（1984—），男，副教授，研究方向为桥梁抗震，E-mail：xi.li@qut.edu.cn

中图分类号: TU375.3
计量
- 文章访问数: 258
- HTML全文浏览量: 185
- PDF下载量: 55
- 被引次数: 0
出版历程
- 收稿日期: 2022-08-07
- 修回日期: 2023-01-08
- 网络出版日期: 2024-07-12
- 刊出日期: 2023-03-29

Research Progress on Buckling of Longitudinal Reinforcement Under Earthquake

LI Xi^{1, 2
,},
YANG Hao¹

1.
School of Civil Engineering, Southwest Jiaotong University, Chengdu 610031, China
2.
School of Civil Engineering, Qingdao University of Technology, Qingdao 266033, China

摘要

摘要:
在地震作用下，钢筋混凝土结构中的墩柱在承受竖向力的同时，通常还会受到往复荷载作用下的水平力，若结构中横向约束不足，纵筋可能会发生明显的横向变形，使得结构承载力和延性显著降低，修复困难. 因此，研究纵筋在受压时的屈曲机理对结构设计以及施工都有着重要的指导意义. 鉴于此，对基于单根钢筋直压试验的纵筋屈曲研究进行总结，分析得出钢筋屈曲的主要影响因素是长细比和屈服强度；概括总结基于RC （reinforced concrete）短柱直压试验的纵筋屈曲问题的研究进展，阐述在复杂相互作用下纵筋屈曲可能受到的各种影响因素；结合土木发展现状，对使用新材料/新构造的墩柱中的纵筋屈曲研究进行讨论. 现阶段对纵筋屈曲的研究中没有给出判断钢筋屈曲方向和范围的理论方法，对纵筋的整体屈曲以及新型结构中的纵筋屈曲问题仍有待进一步研究，且应在后续研究中综合考虑锈蚀、加载历程、截面形式以及各材料相互作用等因素对纵筋屈曲的影响.
- 桥梁工程 /
- 纵筋屈曲 /
- 综述 /
- 钢筋混凝土柱 /
- 抗震设计
Abstract:
During the outbreak of an earthquake, the pier columns in reinforced concrete (RC) structures are usually subjected to horizontal forces and vertical forces under reciprocating loads. When the transverse constraints in the structure are insufficient, the longitudinal reinforcement may have obvious transverse deformation. The bearing capacity and ductility of the structure are thus significantly reduced, and it is difficult to repair. Therefore, it is of great significance for structural design and construction to study the buckling mechanism of longitudinal reinforcements under compression. Firstly, the research on longitudinal reinforcement buckling based on the direct compression test of a single reinforcement was summarized, and the analysis shows that the main influencing factors of reinforcement buckling are the slenderness ratio and yield strength. Secondly, the research progress of reinforcement buckling based on the direct compression test of RC short columns was summarized, and the possible factors influencing longitudinal reinforcement buckling under complex interaction were described. Moreover, combined with the development status of civil engineering, the research progress of longitudinal reinforcement buckling in pier columns with new materials and structures was discussed. Finally, it is concluded that the current research on longitudinal reinforcement buckling lacks a theoretical method for determining the direction and extent of reinforcement buckling. Further investigations are needed to address the overall buckling issues of longitudinal reinforcements and the longitudinal reinforcement buckling concerns in novel structures. Subsequent studies should comprehensively consider the effect of factors such as corrosion, loading history, cross-sectional geometry, and interactions among materials on longitudinal reinforcement buckling.
- bridge engineering /
- longitudinal reinforcement buckling /
- review /
- reinforced concrete columns /
- aseismic design

HTML全文

泡沫铝作为多孔材料的一种，由于其低密度、高孔隙率且能够在稳定的应力水平下发生较大的塑性变形来耗散能量等优点^[1-3]，现已成为一种理想的缓冲吸能材料被广泛应用. 泡沫铝材料在工程应用时既能满足轻量化的设计需求，又能显著提高结构的抗冲击性能，目前更多地应用于汽车、航空等领域，提高载运工具的运行安全.

由于泡沫铝材料受压时表现出良好的吸能特性，引起研究人员对其压缩性能的关注. 目前已经有很多国内外研究学者对影响泡沫铝力学性能的基本因素展开了研究，包括基体材料属性、相对密度和应变率等：郭亚周等^[4]分析了泡沫铝在中低应变率、不同相对密度和不同泡沫铝基体材料下的力学承载性能变化，研究结果表明中低应变率下泡沫铝不具备应变率效应，此外，相对密度较小的泡沫铝具有更高的吸能效率；李忠献等^[5]采用改进的INSTRON高速动力加载系统对闭孔泡沫铝的动态压缩性能展开了研究，结果表明泡沫铝应力随着应变率的提高而增大，材料的能量吸收特性也随之提升；Wang等^[6]对闭孔泡沫铝进行了不同速度下的直接冲击式霍普金森压杆（DHPB）试验，结果表明当冲击速度超过临界值时平台应力动态增强，并将这种增强归因于惯性效应. 此外还有部分学者^[7-9]对于温度对泡沫铝压缩性能的影响展开了研究，李雪艳等^[7]通过进行不同温度下的泡沫铝压缩试验研究，表明闭孔泡沫铝具有温度软化效应，其压溃强度随温度升高而降低，且与温度变化成线性关系. 在泡沫铝的工程应用方面：郭亚周等^[10]通过LS-DYNA开展了机翼前缘未填充和局部填充泡沫铝材料抗鸟撞的分析，结果表明局部填充泡沫铝的机翼前缘能够比空机翼前缘在相同撞击工况下吸收更多的能量；翟希梅等^[11]将泡沫铝填充至6082-T6铝合金圆管中制备了耗能复合构件，通过轴压试验证明填充泡沫铝能够有效改善构件在轴压荷载下的变形能力，避免其发生不规则变形破坏，并提高构件耗能能力.

关于泡沫铝的有限元表征模型，国内外不少学者都开展了一定研究，因为泡沫铝内部孔隙杂乱无章分布且不均匀，建模难度较大且内部孔隙结构难以重建. 泡沫铝三维几何模型通常采用连续的实体模型^[12-13]或者带有随机孔隙分布的实际结构模型^[14-15]：前者通过将泡沫铝压缩力学本构模型赋予实体模型以表征其力学响应，该方法建模简单且计算求解速度快，但无法模拟内部孔隙的准确变形过程，适用于泡沫铝充当吸能结构与其他部件的耦合仿真；后者建模时带有随机分布的孔隙，赋予基体材料力学属性并以此来反映泡沫铝的真实力学属性，较为常用的主要有Kelvin模型和Voronoi模型，此种建模方法难度较大，且计算时间较长，但是可以看出准确的孔隙变形特征.

目前建立的泡沫铝等效有限元模型的研究中尚未较好地考虑应变率效应的问题，本文通过试验和仿真结合的方式建立了能够反映应变率效应的两种密度泡沫铝的等效模型，并将其应用于直升机驾驶舱低空抗坠毁分析当中.

1. 泡沫铝压缩试验

1.1 试验介绍

本文中所用试验试样为两种不同孔隙率的闭孔泡沫铝，制备基体材料均为纯铝. 其中高孔隙率泡沫铝的孔径尺寸较大且分布不均，大部分为2～5 mm，相对密度 ${\rho }^{*}$ 较小，仅为0.080；低孔隙率泡沫铝孔径尺寸相对较小且均匀，大部分为2～3 mm，相对密度 ${\rho }^{*}$ 较大，为0.126. 为区分上述两种泡沫铝材料，称前者为低密度泡沫铝，称后者为高密度泡沫铝，此处的高低仅表示相对关系.

试验试件设计为圆柱形，统一采用线切割进行加工. 考虑到低密度泡沫铝的孔隙直径较大，部分可达5 mm以上. 根据试件设计规则，所有试样在每个方向上都至少有7个胞元^[16]，以避免压缩过程中可能产生的边界效应和尺寸效应，最终试件直径确定为60 mm，高度为15 mm，其中高孔隙率泡沫铝试件如图1所示.

图 1 高孔隙率泡沫铝试样及试件设计

Figure 1. Design of aluminum foam specimen with high porosity

下载: 全尺寸图片幻灯片

本文对两种不同密度的泡沫铝分别进行了0.001 /s、500 /s、1000 /s不同应变率下的压缩试验. 其中准静态试验在Zwick/Roell Z020型万能试验机上进行，试验装置如图2（a）所示. 试验过程中通过力传感器获取力-时间曲线，通过在试验机压头喷射散斑使用非接触测量（DIC）的方式来获取压头位移，进而获取力-位移曲线. 并根据应力 $\sigma =F/A$ 、应变 $\varepsilon =d/H$ 换算得到应力-应变曲线，其中：F为压力；A为试件横截面积；d为压头位移；H为试件高度.

图 2 试验装置

Figure 2. Test apparatus

下载: 全尺寸图片幻灯片

通过分离式霍普金森压杆（SHPB）进行中高速试验，试验装置如图2（b）所示，包括高压充气系统、装置控制系统、高速摄像及灯光系统、应变采集仪、试件固定等部分. 试验之前将试件胶粘固定于输出杆端，根据试验应变率估计气压. 入射杆直径为100 mm，弹性模量为3513 MPa，应力波经入射杆传播后在入射杆端发生反射和透射，通过应变片记录入射、透射和反射的信号得到试件两个端面上应力和应变率的相关信息. 采用经典二波法进行试验数据处理. 整个试验持续约2～3 ms，通过高速摄像机采集试验过程，采集的频率约为20000～40000帧/s.

1.2 试验结果

为保证试验结果的可靠性，每个工况都进行4次重复试验取平均值. 两种密度泡沫铝不同工况下的应力-应变曲线及压缩过程如图3、4所示. 由图3可以看出：泡沫铝压缩应力-应变曲线均表现出明显的三段特征，即线弹性、塑性平台段和密实段.

图 3 不同密度的泡沫铝压缩应力-应变曲线

Figure 3. Stress-strain curves of aluminum foams with different densities

下载: 全尺寸图片幻灯片

图 4 不同密度的泡沫铝准静态压缩过程

Figure 4. Quasi-static compressive processes of aluminum foams with different densities

下载: 全尺寸图片幻灯片

由图4可知：1）低密度泡沫铝试样孔壁更薄，实际加工过程中必然存在孔壁缺失、孔壁连接部位较弱等缺陷，在压缩过程中首先被压溃，压溃区域周围的孔壁应力升高并重新分布，变形带逐渐延伸，表现为“V”形剪切破坏，该过程重复进行，变形带基本上被压实，该过程应力快速升高，压缩进入密实化阶段. 由于整个过程中应力传递较为均匀，压缩的应力应变曲线比较平滑. 2）高密度泡沫铝压缩过程前期基本可以保持完整，但后期压缩变形较大时边界发生碎裂，有较多的基体碎屑散落在试样周围，该现象在高速工况下更为显著. 这是因为高密度泡沫铝内部胞壁的缺陷很少，只能通过孔壁的断裂而非塑性变形来吸能，导致压缩过程中应力传递不均匀，波动较大. 此外随着应变率的增加，平台段应力值提升更快，且平台段的波动更加明显. 这是由于压溃过程中，周围孔壁应力来不及发生重分布，坍塌带迅速扩展导致的. 随着应变率的提升，应力波动现象愈发明显.

对于泡沫铝等耗能缓冲构件，常采用比吸能和平台应力作为评价其承载吸能能力的重要指标. 其中比吸能（SEA）代表结构单位质量或单位体积下的吸能能力，比吸能越大，泡沫铝的承载吸能能力越好. 质量比吸能及体积比吸能的计算方法如式（1）.

$\begin{split} {S_m} = \dfrac{{{E_{\rm{A}}}}}{m},\\ {S_V} = \dfrac{{{E_{\rm{A}}}}}{V}, \end{split}$

(1)

式中：m和V分别为泡沫铝的质量和体积；E_A为压缩开始至密实化应变区间内泡沫铝的总吸能量，本文取应力-应变曲线中对应1.3倍平台应力^[16]处的应变值作为密实化应变.

平台应力（ ${\sigma }_{{\rm{pl}}}$ ）作为评估泡沫铝压溃过程中吸收能量的关键参数，如式（2）.

${\sigma }_{{\rm{pl}}}=\frac{1}{{\varepsilon }_{2}-{\varepsilon }_{1}}{\int }_{{\varepsilon }_{1}}^{{\varepsilon }_{2}}\sigma \left(\varepsilon \right){\rm{d}}\varepsilon \text{，}$

(2)

式中： ${\varepsilon }_{1}$ 和 ${\varepsilon }_{2}$ 分别为泡沫铝轴向承载时平台段的起始应变和密实化应变.

不同压缩应变率工况下泡沫铝的压缩结果如表1所示，各指标均为4次压缩试验的均值.

表 1 不同密度泡沫铝压缩试验结果

Table 1. Test results of aluminum foams with different densities

试件	应变率/（s⁻¹）	平台应力 /MPa	密实化应变	质量比吸能/ （kJ·kg⁻¹）
低密度试件	0.001	0.95	0.57	1.88
	500	1.31	0.53	2.75
	1000	1.37	0.50	2.89
高密度试件	0.001	3.80	0.55	4.24
	500	4.73	0.49	5.36
	1000	4.85	0.48	5.42

下载: 导出CSV

| 显示表格

1.3 相对密度及应变率对承载能力的影响规律

泡沫铝材料的初始相对密度以及应变率对其力学性能承载指标影响显著，其影响规律如图5所示. 由图5可知：1）无论哪种相对密度的泡沫铝，平台应力以及质量比吸能都随应变率的增加而增加，表现出明显的应变率敏感性. 随着相对密度的增大，两个指标的提升现象变得不明显. 2）无论哪种相对密度的泡沫铝，其压缩过程中的密实化应变均随应变率的增加而呈现出下降趋势. 这是因为在泡沫铝的压溃过程中，孔隙壁的屈曲和坍塌是主要的吸能方式，在低应变率下，孔隙壁有足够的时间来进行旋转变形使得遗留孔隙最小，然而在高应变率下，在下一阶段孔隙壁变形之前，前一阶段孔隙壁没有足够的时间变形至最终模式，所以导致致密模式发生改变，引起了致密化应变的减小，这与文献[17]中的结论一致.

图 5 相对密度及应变率对不同密度泡沫铝承载吸能的影响规律

Figure 5. Influence of relative density and strain rate on the load-carrying and energy absorption of aluminum foams with different densities

下载: 全尺寸图片幻灯片

2. 考虑应变率效应的泡沫铝压缩等效模型建立

2.1 等效模型建立

为了评估泡沫铝在直升机驾驶舱低空坠落工况下的承载吸能效果，需要建立可用于驾驶舱垂直跌落仿真的泡沫铝模型. 考虑到泡沫铝内部由无规律分布的孔隙组成，如果建立其原样准确模型用于驾驶舱跌落的仿真，泡沫铝较小的网格尺寸与驾驶舱模型中的大网格尺寸的时间步长不匹配，导致耦合仿真计算时间难以接受. 因此，本文基于HyperMesh和LS-DYNA平台建立了可表征应变率效应的不同密度的泡沫铝有限元等效模型. 该等效模型在进行耦合仿真时可以大幅减少计算时间，并且可以真实反映出泡沫铝压缩试验中所表现出的力学特性以及应变率效应以保证仿真结果的准确性，能在实际工程中进行应用. 泡沫铝等效有限元模型如图6所示，该模型分为3部分：泡沫铝试件、顶部加载刚性板以及底部约束刚性板.

图 6 泡沫铝等效模型示意

Figure 6. Equivalent FE model of aluminum foam

下载: 全尺寸图片幻灯片

为了验证泡沫铝等效模型的有效性，利用上述建模方式建立了与试验尺寸信息一致的有限元模型. 泡沫铝模型为 $\phi 60\times 15$ mm的圆柱体. 泡沫铝模型置于上下两块刚性平板中间. 下刚性板施加全固定约束，上刚性板只有压缩方向1个移动自由度. 对上刚性板施加强制位移使其冲击泡沫铝模型以模拟压缩工况.

由于试验过程中泡沫铝试件发生横向膨胀变形后与加载刚性板之间的接触界面节点发生了改变，因此选用自动面-面接触来模拟泡沫铝与上下刚性板的接触，摩擦系数设为0.3. 对泡沫铝施加CONTACT_INTERIOR接触防止出现负体积. 上下刚性板以及泡沫铝模型均采用八节点六面体单元进行网格划分，积分方式为0号，沙漏类型为2号^[18].

2.2 材料模型

采用刚性材料（MAT_20）模拟两端刚性板. 泡沫铝采用MAT_MODIFIED_HONEYCOMB本构进行模拟，对于所有的法向应力和切向应力，可以分别定义其非线性弹塑性行为. 通过在3个方向上指定相同的材料属性可以实现各向同性.

该材料模型遵循式（3）所示的应力强化规律^[19].

$\sigma ={\sigma }_{{\rm{p}}} + \gamma \frac{\widehat{\varepsilon }}{{\varepsilon }_{{\rm{D}}}} + \alpha {\rm{ln}}\left[\frac{1}{1-{\left(\widehat{\varepsilon }/{\varepsilon }_{{\rm{D}}}\right)}^{\beta }}\right] \text{，}$

(3)

式中： $\sigma$ 和 $\widehat{\varepsilon }$ 分别为结构当前的应力和应变； ${\sigma }_{{\rm{p}}}$ 为初始平台应力；等式右边第二项为线性应变硬化项，其中线性应变硬化系数 $\gamma$ 为应力-应变曲线的斜率， ${\varepsilon }_{{\rm{D}}}$ 为密实化应变， ${\varepsilon }_{{\rm{D}}}=-\mathrm{l}\mathrm{n}\left({\rho }_{{\rm{f}}}/{\rho }_{{\rm{f0}}}\right)$ ， ${\rho }_{{\rm{f}}}$ 为泡沫铝的密度， ${\rho }_{{\rm{f0}}}$ 为泡沫铝基体材料的密度；等式右边最后一项表示非线性应变硬化，由系数 $\alpha$ （比例系数）和 $\beta$ （形状系数）表示.

使用Origin自带的用户自定义拟合函数的非线性拟合工具对两种不同密度泡沫铝准静态压缩工况下的应力-应变曲线进行本构方程式（3）的拟合，具体的拟合方法是将初始平台应力作为定值，其余4个参数作为变量，采用Levenberg-Marquardt优化迭代算法进行迭代求得最优的拟合参数. 两种密度泡沫铝的本构拟合参数如表2所示，其中： ${E}_{{\rm{AAU}}}、 {E}_{{\rm{BBU}}}、{E}_{{\rm{CCU}}}$ 分别为3个正交方向下的弹性模量； ${G}_{{\rm{AAU}}}、 {G}_{{\rm{BBU}}}、{G}_{{\rm{CCU}}}$ 为3个方向下的剪切模量. 结果如图7所示.

表 2 不同密度的泡沫铝材料参数

Table 2. Material parameters of aluminum foam with different density

密度	${\sigma }_{{\rm{p}}}$ /MPa	$\gamma$	${\varepsilon }_{{\rm{D}}}$	$\alpha$	$\beta$	密度 $\rho /$ （kg/m³）	弹性模量/MPa ${E}_{{\rm{AAU}}}={E}_{{\rm{BBU}}}={E}_{{\rm{CCU}}}$	剪切模量/MPa ${G}_{{\rm{AAU}}}={G}_{{\rm{BBU}}}={G}_{{\rm{CCU}}}$	应力应变曲线
低	0.97	−1.4	0.77	1.28	1.17	216	300	300	图7（a）
高	4.34	−10.1	0.82	8.46	1.23	340	800	800	图7（b）

下载: 导出CSV

| 显示表格

此外该本构的杨氏模量和剪切模量从未开始压缩状态到压实状态随体积应变变化而线性变化，初始状态值需要用户给定. 两种不同密度泡沫铝的初始弹性模量由加载卸载试验来确定，试验结果表明低密度泡沫初始弹性模量为300 MPa，高密度泡沫初始弹性模量为800 MPa. 将两种密度泡沫铝准静态压缩下的应力-应变曲线作为基准应变率下的应力-应变曲线，两种密度泡沫铝的材料参数如表2所示.

图 7 泡沫铝本构拟合结果

Figure 7. Constitutive fitting results of aluminum foam

下载: 全尺寸图片幻灯片

该本构对于应变率效应的表征是通过首先确定不同应变率相较于基准应变率的缩放系数，然后对基准应变率下的应力-应变曲线进行缩放得到不同应变率下的力学响应曲线. 通过试验结果，将两种密度泡沫铝基准应变率（0.001 /s）下得到的压缩应力-应变曲线分别进行缩放以拟合相同密度泡沫铝在500 /s以及1000 /s应变率下的应力-应变曲线. 因为平台段作为泡沫铝压缩承载的关键阶段，因此，缩放标准要保证缩放曲线与原曲线的平台段基本一致，缩放拟合结果如图8所示，得到了两种不同密度泡沫铝在0.001 /s、500 /s、1000 /s应变率下的缩放系数.

图 8 对泡沫铝在0.001 /s下的应力-应变曲线缩放结果的拟合

Figure 8. Fitting on the scaling results of the stress-strain curve at 0.001 /s

下载: 全尺寸图片幻灯片

为了表征泡沫铝的应变率效应，需要在本构方程中采用准确的应变率项. 在确定了两种密度泡沫铝基准应变率（0.001 /s）、500 /s、1000 /s应变率下的数据点之后，使用Johnson-Cook（J-C）以及Cowper-Symonds（C-S）两种常用应变率模型进行标定. 其中，J-C模型应变率项采用对数函数表征应变率效应，而C-S模型采用幂函数来表征应变率效应，两者的标定结果如图9所示.

图 9 缩放系数随应变率变化标定结果

Figure 9. Calibrated results of scaling factors varying with strain rate

下载: 全尺寸图片幻灯片

由图9可以看出：1） J-C模型对高密度泡沫铝的应变率效应描述较好，但对于低密度泡沫铝应变率的描述则与试验结果有较大的偏差. 这是由于泡沫铝不同应变率下的缩放系数在对数坐标下并非线性，而J-C模型应变率项为线性函数导致的. 2）反观C-S模型的标定结果与试验结果吻合较好，相关系数达到0.99，故采用C-S模型作为最终泡沫铝应变率效应的表征函数，如式（4），标定系数如表3所示. 将该标定结果曲线作为泡沫铝等效模型应变率项的输入以准确表征泡沫铝的应变率效应.

表 3 应变率模型标定系数

Table 3. Calibrated coefficients of the strain-rate model

泡沫铝	b	e	相关系数（R²）
低密度	5.969 × 10⁶	0.485	0.999
高密度	2.112 × 10⁸	0.248	0.992

下载: 导出CSV

| 显示表格

${f}_{{\rm{rate}}}=1 + {\left(\dfrac{\dot \bar \varepsilon _{\rm{p}}^*}{b}\right)}^{e},$

(6)

式中： ${f}_{{\rm{rate}}}$ 为不同应变率下的缩放系数； ${\dot \bar \varepsilon _{\rm{p}}^*}$ 为等效塑性应变率； e和b $\mathrm{为}$ 需要拟合的材料参数.

2.3 等效模型结果验证

由于0.001 /s应变率下的仿真计算时间过长，因此本节只进行500 /s和1000 /s中高应变率下的仿真结果验证. 两种密度所对应的500 /s和1000 /s应变率下的仿真和试验应力-应变曲线的结果对比如图10所示. 由图10可以看到：除了仿真结果因为初始接触速度大而有所振荡之外，吸能的关键区域平台段曲线结果吻合较好，可以说明该模型具有较高的精度，可以代替真实泡沫铝模型进行仿真从而节约大量的计算时间.

图 10 仿真与实验结果对比

Figure 10. Simulation and experimental results for aluminum foam with low relative density

下载: 全尺寸图片幻灯片

3. 泡沫铝在直升机驾驶舱抗坠毁设计中的应用

基于建立的可反映应变率效应的不同密度泡沫铝的等效模型，以某型号直升机结构为建模依据，建立了简化后的直升机驾驶舱低空垂向跌落的有限元模型. 通过进行驾驶舱有限元模型和泡沫铝等效模型的耦合跌落仿真，研究泡沫铝密度对于驾驶舱垂向跌落的影响. 该模型由蒙皮、长桁、隔框、底板、横梁、纵梁、支撑板件、泡沫铝等效模型以及底部刚性墙组成，有限元模型如图11所示. 其中泡沫铝的尺寸是根据横纵梁交错位置的空间大小确定的，最终尺寸为180 mm × 176 mm × 300 mm，其安装位置是驾驶舱底部的横纵梁交接空隙处.

图 11 直升机驾驶舱有限元模型

Figure 11. Finite element model of helicopter cockpit

下载: 全尺寸图片幻灯片

该有限元模型各部分的连接关系为：长桁、隔框、底板、横纵梁等内部结构采用共节点接触；长桁、隔框、横纵梁等均与外部蒙皮采用MAT-100 Spotweld焊点连接来模拟铆接；内部泡沫铝等效模型与蒙皮和底板之间采用自动面面接触，驾驶舱模型整体采用自接触，并且与刚性墙之间采用自动面面接触. 对底部刚性墙进行完全固定约束，对驾驶舱模型施加大小为1、10、20 m/s的初始跌落速度以模拟驾驶舱模型低空跌落的工况. 3种跌落速度反映在泡沫铝模型上，对应的单元应变率分别为1 /s、10 /s、30 /s左右. 该模型中，驾驶舱结构均为铝合金材料，采用MAT 24本构进行模拟，统一采用四节点壳单元进行网格划分，泡沫铝等效模型建模与第2节所述一致.

无泡沫铝、置入低、高密度泡沫铝模型的低空跌落仿真变形不同. 以20 m/s的初始跌落速度为例，相同时刻下3种模型的仿真过程如图12所示. 由图12可知：无泡沫铝驾驶舱下部结构的蒙皮、横纵梁首先受到冲击，出现屈曲塑性变形，随着撞击过程继续，撞击能量继续向上传递，底部基本压溃，上部的隔框、长桁以及蒙皮开始发生塑性折弯变形来吸能；填充泡沫铝的驾驶舱模型，横纵梁发生屈曲塑性变形的同时，填充泡沫铝受压发生塑性压缩变形吸收了很大一部分碰撞动能；相较无泡沫铝的驾驶舱模型，相同时刻下底板下部结构变形减少，如红色椭圆部分所示.

图 12 3种模型跌落仿真过程对比

Figure 12. Comparison of simulation results for the three models under dropping conditions

下载: 全尺寸图片幻灯片

3种不同跌落速度下的最大变形量具体数据如图13所示，泡沫铝起到了明显的缓冲效果，底舱最大变形量得到有效减少. 此外随着跌落速度的增加，碰撞工况更加剧烈，底舱的变形量逐渐增大. 由于泡沫铝的填充，驾驶舱底部的承载能力提高，上部隔框等承载部件更早地进入塑性变形，变形加大，如图12黄色框线部分所示. 由于高密度泡沫铝承载能力高于低密度泡沫铝，因此相同时刻下其底部结构变形更小，但隔框等结构变形更严重.

图 13 驾驶舱底部最大变形量结果

Figure 13. Maximum deformation of the bottom of the cockpit

下载: 全尺寸图片幻灯片

驾驶舱结构跌落中主要依靠机身结构各部件发生不同程度的塑性变形来吸收碰撞动能，主要吸能形式是底部横纵梁组件的弯曲、内部泡沫铝的塑性压缩变形以及蒙皮、隔框的塑性折弯变形. 3种跌落速度工况下，内部部件的吸能规律基本保持一致，典型工况20 m/s跌落过程中各个主要部件的质量比吸能如图14所示. 可以看出，通过置入泡沫铝结构在压溃过程中吸收大量的碰撞动能，大大缓解了驾驶舱原有结构的承载压力，原有结构吸收能量减少. 这就意味着原有主要承载结构如横纵梁等的设计要求可以降低，承载能力留有一定的冗余，安全性能得到进一步保证. 随着跌落速度的增加，模型输入的初始碰撞动能增加. 跌落过程中，各个承载部件的变形更加严重，吸收能量也有所提高. 以泡沫铝结构为例，不同跌落速度下的泡沫铝质量比吸能如图14（b）所示，可以看出，随着初始跌落速度的增加，不论何种密度的泡沫铝的质量比吸能均随着跌落速度的增加而增加. 此外相对密度更大的泡沫铝在实际承载中体现出更高的能量吸收能力.

图 14 驾驶舱垂向跌落仿真结果

Figure 14. Simulation results of the cockpit under dropping conditions

下载: 全尺寸图片幻灯片

4. 结　论

本文首先进行了不同应变率、不同相对密度泡沫铝的压缩试验，研究了不同相对密度的闭孔泡沫铝的力学性能；然后，建立了可放映应变率效应的泡沫铝等效有限元模型并应用于直升机耐坠性设计中，得出了以下结论：

1）低密度泡沫铝压缩过程更为稳定，高密度泡沫铝压缩过程波动较大，且有碎屑脱落.

2）泡沫铝压缩时的平台应力和质量比吸能随相对密度及应变率的增加而增大，但密实化应变相反，随相对密度或应变率增加而减小.

3）基于MAT_MODIFIED_HONEYCOMB本构建立了可反映应变率效应的不同密度泡沫铝的等效有限元模型，经过试验结果验证具有较高的精度.

4）置入泡沫铝后驾驶舱原有各部件的承载压力和底舱的最大变形量均得到有效的减小，内部乘客的生存空间可以得到有效的保证，高密度泡沫铝相较低密度泡沫铝吸能效果更佳.

图 1 长细比对平均响应的影响

Figure 1. Effect of slenderness ratio on average response

下载: 全尺寸图片幻灯片

图 2 屈服强度对平均响应的影响

Figure 2. Effect of yield strength on average response

下载: 全尺寸图片幻灯片

图 3 Bae等^[9]的试验分析结果

Figure 3. Experimental and analysis results of Bae et al^[9].

下载: 全尺寸图片幻灯片

图 4 未锈蚀与锈蚀钢筋的屈曲机理^[15]

Figure 4. Buckling mechanism of uncorroded and corroded reinforcements ^[15]

下载: 全尺寸图片幻灯片

图 5 Montil-Nuti模型

Figure 5. Montil-Nuti model

下载: 全尺寸图片幻灯片

图 6 Dhakal-Maekawa模型

Figure 6. Dhakal-Maekawa model

下载: 全尺寸图片幻灯片

图 7 RC柱中的纵筋屈曲类型

Figure 7. Types of longitudinal reinforcement buckling in RC columns

下载: 全尺寸图片幻灯片

图 8 G-A本构模型

Figure 8. G-A constitutive model

下载: 全尺寸图片幻灯片

图 9 修正的G-A模型

Figure 9. Modified G-A model

下载: 全尺寸图片幻灯片

图 10 屈曲约束套管

Figure 10. Buckling-restrained casing

下载: 全尺寸图片幻灯片

图 11 屈曲约束套管^[62]

Figure 11. Buckling-restrained casing proposed^[62]

下载: 全尺寸图片幻灯片

参数 $\; \beta$ 的选取

Selection of parameter $\;\beta$

条件 $L / D \leqslant 5$ $5<L / D \leqslant 10$ $L / D>10$

取值 0.0050 0.0065 0.0125

下载: 导出CSV

表 2 各国规范对比

Table 2. Comparison of norms among different countries

规定箍筋
最小数量箍筋间距构造措施

欧洲规范
EN 1998 √ √ √

新西兰规范
NZS 3101 √ √ √

美国规范
ACI 318—05 √

我国规范 √ √

下载: 导出CSV

参考文献(62)

[1]	张敏政. 近年地震震害的几点启示[J]. 工程抗震,2001,23(1): 11-15,29. ZHANG Minzheng. Some Revelations of the earthquake disaster in recent years[J]. Earthquake Resistant Engineering, 2001, 23(1): 11-15,29.
[2]	文明秀. 地震工程与抗震原理探究[J]. 现代交际,2015(2): 90. doi: 10.3969/j.issn.1009-5349.2015.02.070
[3]	MONTI G, NUTI C. Nonlinear cyclic behavior of reinforcing bars including buckling[J]. Journal of Structural Engineering, 1992, 118(12): 3268-3284. doi: 10.1061/(ASCE)0733-9445(1992)118:12(3268)
[4]	DHAKAL R P, MAEKAWA K. Modeling for postyield buckling of reinforcement[J]. Journal of Structural Engineering, 2002, 128(9): 1139-1147. doi: 10.1061/(ASCE)0733-9445(2002)128:9(1139)
[5]	BAE S, MIESES A M, BAYRAK O. Inelastic buckling of reinforcing bars[J]. Journal of Structural Engineering, 2005, 131(2): 314-321. doi: 10.1061/(ASCE)0733-9445(2005)131:2(314)
[6]	COSENZA E, PROTA A. Experimental behaviour and numerical modelling of smooth steel bars under compression[J]. Journal of Earthquake Engineering, 2006, 10(3): 313-329.
[7]	陈昉健,易伟建. 高强纵筋受压屈曲试验研究[J]. 建筑科学,2019,35(1): 75-81. CHEN Fangjian, Yi Weijian. Experimental study of high-strength rebar buckling under compression[J]. Building Science, 2019, 35(1): 75-81.
[8]	KORENTZ J. The effect of yield strength on inelastic buckling of reinforcing bars[J]. Mechanics and Mechanical Engineering, 2010, 14(2): 247-255.
[9]	KORENTZ J, MARCINOWSKI J. Effect of mechanical parameters of steel on inelastic buckling of reinforcing bars[J]. International Journal of Structural Stability and Dynamics, 2016, 16(7): 1550028.1-1550028.12.
[10]	MASSONE L M, MORODER D. Buckling modeling of reinforcing bars with imperfections[J]. Engineering Structures, 2009, 31(3): 758-767. doi: 10.1016/j.engstruct.2008.11.019
[11]	史庆轩,牛荻涛,颜桂云. 反复荷载作用下锈蚀纵筋混凝土压弯构件恢复力性能的试验研究[J]. 地震工程与工程振动,2000,20(4): 44-50. SHI Qingxuan, NIU Ditao, YAN Guiyun. Experimental research on hysteretic characteristics of corroded RC member with flexural and compressive axial loads under repeated horizontal loading[J]. Earthquake Engineering and Engineering Vibration, 2000, 20(4): 44-50.
[12]	KASHANI M M, BARMI A K, MALINOVA V S. Influence of inelastic buckling on low-cycle fatigue degradation of reinforcing bars[J]. Construction and Building Materials, 2015, 94: 644-655. doi: 10.1016/j.conbuildmat.2015.07.102
[13]	KASHANI M M, CREWE A J, ALEXANDER N A. Nonlinear cyclic response of corrosion-damaged reinforcing bars with the effect of buckling[J]. Construction and Building Materials,2013,41:388-400.
[14]	KASHANI M M, LOWES L N, CREWE A J, et al. Phenomenological hysteretic model for corroded reinforcing bars including inelastic buckling and low-cycle fatigue degradation[J]. Computers & Structures,2015,156:58-71.
[15]	KASHANI M M, CREWE A J, ALEXANDER N A. Nonlinear stress-strain behaviour of corrosion-damaged reinforcing bars including inelastic buckling[J]. Engineering Structures, 2013, 48: 417-429. doi: 10.1016/j.engstruct.2012.09.034
[16]	陈武雄. FRP加固锈蚀RC柱力学性能及其锈蚀纵筋屈曲特性的试验研究[D]. 深圳:深圳大学,2015.
[17]	IMPERATORE S, RINALDI Z. Experimental behavior and analytical modeling of corroded steel rebars under compression[J]. Construction and Building Materials, 2019, 226: 126-138. doi: 10.1016/j.conbuildmat.2019.07.109
[18]	DHAKAL R P, MAEKAWA K. Path-dependent cyclic stress-strain relationship of reinforcing bar including buckling[J]. Engineering Structures, 2002, 24(11): 1383-1396. doi: 10.1016/S0141-0296(02)00080-9
[19]	FILIPPOU F C, POPOV E P, BERTERO V V. Effects of bond deterioration on hysteretic behavior of reinforced concrete joints[R]. Berkeley: Earthquake Engineering Research Center, University of California, 1983.
[20]	LIMBERT J, AFSHAN S, KASHANI M M, et al. Compressive stress-strain behaviour of stainless steel reinforcing bars with the effect of inelastic buckling[J]. Engineering Structures, 2021, 237: 112098.1-112098.15.
[21]	陈进可. 考虑屈曲及混凝土横向膨胀影响的钢筋本构模型改进方法研究[D]. 重庆:重庆大学,2014.
[22]	PANTAZOPOULOU S J. Detailing for reinforcement stability in RC members[J]. Journal of Structural Engineering, 1998, 124(6): 623-632. doi: 10.1061/(ASCE)0733-9445(1998)124:6(623)
[23]	BAYRAK O , SHEIKH S A. Plastic hinge analysis[J]. Journal of Structural Engineering,2001,127(9):1092-1100.
[24]	DHAKAL R P, MAEKAWA K. Reinforcement stability and fracture of cover concrete in reinforced concrete members[J]. Journal of Structural Engineering,2002,128(10):1253-1262.
[25]	HILSON C W, SEGURA C L, WALLACE J W. Experimental study of longitudinal reinforcement buckling in reinforced concrete structural wall boundary elements[C]//Proceedings of the 10th National Conference in Earthquake. Anchorage: Earthquake Engineering Research Institute, 2014.
[26]	FRAGIADAKIS M, PINHO R, ANTONIOU S. Modelling inelastic buckling of reinforcing bars under earthquake loading[C]//ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering. Greece: ECCOMAS , 2007.
[27]	MANDER J B, PANTHAKI F D, KASALANATI A. Low-cycle fatigue behavior of reinforcing steel[J]. Journal of Materials in Civil Engineering, 1994, 6(4): 453-468.
[28]	ZONG Z, KUNNATH S. Simulation of reinforcing bar buckling in circular reinforced concrete columns[J]. Aci Structural Journal,2013,110(4):607-616.
[29]	ZONG Z Y, KUNNATH S, MONTI G. Material model incorporating buckling of reinforcing bars in RC columns[J]. Journal of Structural Engineering, 2014, 140(1): 04013032.1-04013032.10.
[30]	邓宗才,姚军锁. 高强箍筋约束超高性能混凝土柱轴压性能[J]. 复合材料学报,2020,37(10): 2590-2601. DENG Zongcai, YAO Junsuo. Axial compression behavior of ultra-high performance concrete columns confined by high-strength stirrups[J]. Acta Materiae Compositae Sinica, 2020, 37(10): 2590-2601.
[31]	MASSONE L M, LÓPEZ E E. Modeling of reinforcement global buckling in RC elements[J]. Engineering Structures, 2014, 59: 484-494. doi: 10.1016/j.engstruct.2013.11.015
[32]	蒋连接,常明丰,文兆全. 锈蚀钢筋力学性能的试验研究[J]. 低温建筑技术,2011(3): 62-64. JIANG Lianjie, CHANG Mingfeng, WEN Zhaoquan. Experimental study on mechanical properties of corroded steel bars[J]. Low Temperature Architecture Technology, 2011(3): 62-64.
[33]	李凤兰,侯维玲,侯朋兵. 锈蚀钢筋的力学性能试验研究[J]. 华北水利水电学院学报,2013,34(4): 61-64. LI Fenglan, HOU Weiling, HOU Pengbing. Experimental study on mechanical properties of corroded steel bars[J]. Journal of North China Institute of Water Conservancy and Hydroelectric Power, 2013, 34(4): 61-64.
[34]	MA Y, CHE Y, GONG J X. Behavior of corrosion damaged circular reinforced concrete columns under cyclic loading[J]. Construction and Building Materials,2012,29:548-556.
[35]	LAVORATO D, FIORENTINO G, PELLE A, et al. A corrosion model for the interpretation of cyclic behavior of reinforced concrete sections[J]. Structural Concrete, 2020, 21(5): 1732-1746. doi: 10.1002/suco.201900232
[36]	SUDA K. Buckling behavior of longitudinal reinforcing bars in concrete columns subjected to reversed lateral loading[C]//Proceedings of the 11th World Conference on Earthquake Engineering. Mexico: Elsevier Science Ltd.,1996.
[37]	RODRIGUEZ M E, BOTERO J C, VILLA J. Cyclic stress-strain behavior of reinforcing steel including effect of buckling[J]. Journal of Structural Engineering, 1999, 125(6): 605-612. doi: 10.1061/(ASCE)0733-9445(1999)125:6(605)
[38]	MOYER M J, KOWALSKY M J. Influence of tension strain on buckling of reinforcement in concrete columns[J]. ACI Structural Journal, 2003, 100(1): 75-85.
[39]	ZONG Z, KUNNATH S. Buckling of reinforcing bars in concrete structures under seismic loads[C]//The 14th World Conference on Earthquake Engineering. Beijing: [s.n.], 2008:1-8.
[40]	CEN. Eurocode 8: design of structures for earthquake resistance-Part 2: Bridges:EN 1998-2:2005[S]. London:BSI Standards Limited,2005.
[41]	Standards New Zealand. Concrete structures standard. part 1: the design of concrete structures: NZS 3101.1: 2006[S]. New Zealand: New Zealand Standards Executive, 2006.
[42]	中华人民共和国住房和城乡建设部. 建筑抗震设计规范:GB 50011—2010[S]. 北京:中国建筑工业出版社,2010.
[43]	美国混凝土学会(ACI). 美国房屋建筑混凝土结构规范: ACI 318R-05[S]. 张川,白绍良,钱觉时译. 重庆:重庆大学出版社,2007.
[44]	招商局重庆交通科研设计院有限公司. 公路桥梁抗震设计规范:JTG/T 2231-01—2020[S]. 北京:人民交通出版社.
[45]	仇建磊,贡金鑫. 地震作用下钢筋混凝土柱纵筋屈曲研究进展[J]. 土木工程学报,2016,49(5): 50-62. QIU Jianlei, GONG Jinxin. Research on buckling of longitudinal bars in RC columns under earthquake excitation: state of the art[J]. China Civil Engineering Journal, 2016, 49(5): 50-62.
[46]	GOMES A, APPLETON J. Nonlinear cyclic stress-strain relationship of reinforcing bars including buckling[J]. Engineering Structures, 1997, 19(10): 822-826. doi: 10.1016/S0141-0296(97)00166-1
[47]	BROWN J R, KUNNATH S. Low-cycle fatigue behavior of longitudinal reinforcement in reinforced concrete bridge columns[R]. New York: State University of New York, 2000.
[48]	BOURNAS D A, TRIANTAFILLOU T C. Bar buckling in RC columns confined with composite materials[J]. Journal of Composites for Construction, 2011, 15(3): 393-403. doi: 10.1061/(ASCE)CC.1943-5614.0000180
[49]	BAI Y L, DAI J G, OZBAKKALOGLU T. Cyclic stress-strain model incorporating buckling effect for steel reinforcing bars embedded in FRP-confined concrete[J]. Composite Structures, 2017, 182: 54-66. doi: 10.1016/j.compstruct.2017.09.007
[50]	SATO Y, KO H. Modeling of reinforcement buckling in RC columns confined with FRP[J]. Journal of Advanced Concrete Technology, 2008, 6(1): 195-204. doi: 10.3151/jact.6.195
[51]	BAI Y L, DAI J G, TENG J G. Buckling of steel reinforcing bars in FRP-confined RC columns: an experimental study[J]. Construction and Building Materials, 2017, 140: 403-415. doi: 10.1016/j.conbuildmat.2017.02.149
[52]	BAI Y L, DAI J G, TENG J G. Monotonic stress-strain behavior of steel rebars embedded in FRP-confined concrete including buckling[J]. Journal of Composites for Construction, 2017, 21(5): 04017043.1-04017043.11.
[53]	BOURNAS D A, LONTOU P V, PAPANICOLAOU C G, et al. Textile-reinforced mortar versus fiber-reinforced polymer confinement in reinforced concrete columns[J]. ACI Structural Journal, 2007, 104(6): 740-748.
[54]	BOURNAS D A, TRIANTAFILLOU T C, PAPANICOLAOU C G. Retrofit of seismically deficient RC columns with textile- reinforced mortar (TRM) jackets[C]//4^th Colloquium on Textile Reinforced Structures (CTRS4). Dresden:[s.n.], 2009:471-490.
[55]	BOURNAS D A, TRIANTAFILLOU T C, ZYGOURIS K, et al. Textile-reinforced mortar versus FRP jacketing in seismic retrofitting of RC columns with continuous or lap-spliced deformed bars[J]. Journal of Composites for Construction, 2009, 13(5): 360-371. doi: 10.1061/(ASCE)CC.1943-5614.0000028
[56]	HUNG C C, YEN C H. Compressive behavior and strength model of reinforced UHPC short columns[J]. Journal of Building Engineering, 2021, 35: 102103.1-102103.14.
[57]	胡锐,方志,许宝丹. 不同配箍率和钢纤维掺量UHPC柱抗震性能试验[J]. 中国公路学报,2021,34(8): 65-77. doi: 10.3969/j.issn.1001-7372.2021.08.007 HU Rui, FANG Zhi, XU Baodan. Experimental study on seismic behavior of UHPC columns with different stirrup ratios and steel fiber contents[J]. China Journal of Highway and Transport, 2021, 34(8): 65-77. doi: 10.3969/j.issn.1001-7372.2021.08.007
[58]	LI X, CHEN K D, HU P, et al. Effect of ECC jackets for enhancing the lateral cyclic behavior of RC bridge columns[J]. Engineering Structures, 2020, 219: 110714.1-110714.17.
[59]	RUANGRASSAMEE A, SAWAROJ A. Seismic enhancement of reinforced-concrete columns by rebar-restraining collars[J]. Journal of Earthquake and Tsunami, 2012, 6(3): 1250015.1-1250015.19.
[60]	LUKKUNAPRASIT P, TANGBUNCHOO T, RODSIN K. Enhancement of seismic performance of reinforced concrete columns with buckling-restrained reinforcement[J]. Engineering Structures, 2011, 33(12): 3311-3316. doi: 10.1016/j.engstruct.2011.08.019
[61]	伍云天,潘毅,陈博文,等. 采用屈曲约束纵筋的钢筋混凝土短柱轴压性能试验研究[J]. 土木工程学报,2013,46(增1): 190-194. WU Yuntian, PAN Yi, CHEN Bowen, et al. Behavior of axially loaded RC columns with buckling restrained longitudinal rebars[J]. China Civil Engineering Journal, 2013, 46(S1): 190-194.
[62]	张文涛,贾俊峰,曹艳辉. 基于纵筋屈曲约束构造的RC柱轴压性能试验研究[J]. 土木工程学报,2018,51(增2): 62-67. ZHANG Wentao, JIA Junfeng, CAO Yanhui. Experimental study on axial compression behavior of RC columns based on buckling restrained reinforcements[J]. China Civil Engineering Journal, 2018, 51(S2): 62-67.

施引文献

附加材料(0)

访问统计

点击查看大图

图(11) / 表(2)

计量

文章访问数: 258
HTML全文浏览量: 185
PDF下载量: 55
被引次数: 0

1. 泡沫铝压缩试验
1.1 试验介绍
1.2 试验结果
1.3 相对密度及应变率对承载能力的影响规律
2. 考虑应变率效应的泡沫铝压缩等效模型建立
2.1 等效模型建立
2.2 材料模型
2.3 等效模型结果验证
3. 泡沫铝在直升机驾驶舱抗坠毁设计中的应用
4. 结　论

1. 泡沫铝压缩试验
1.1 试验介绍
1.2 试验结果
1.3 相对密度及应变率对承载能力的影响规律
2. 考虑应变率效应的泡沫铝压缩等效模型建立
2.1 等效模型建立
2.2 材料模型
2.3 等效模型结果验证
3. 泡沫铝在直升机驾驶舱抗坠毁设计中的应用
4. 结　论

参考文献(62)

施引文献

附加材料(0)

访问统计

条件	$L / D \leqslant 5$	$5<L / D \leqslant 10$	$L / D>10$
取值	0.0050	0.0065	0.0125

规定	箍筋最小数量	箍筋间距	构造措施
欧洲规范 EN 1998	√	√	√
新西兰规范 NZS 3101	√	√	√
美国规范 ACI 318—05		√
我国规范	√	√

地震作用下纵筋屈曲的研究进展

doi: 10.3969/j.issn.0258-2724.20220549

作者简介: 李晰（1984—），男，副教授，研究方向为桥梁抗震，E-mail：xi.li@qut.edu.cn

计量

出版历程