落锤冲击气囊施加均布动荷载的试验方法

潘腾　陈力　方秦

摘要：在结构构件上施加均布冲击（爆炸）荷载一直是动力试验领域的难题。提出了一种落锤冲击气囊施加均布动荷载的试验装置，和一个双自由度弹簧阻尼冲击加载简化分析模型；给出了该模型的基本力学方程、边界条件和求解方法，并进行了有限元验证；建立了所提出冲击加载试验装置的精细化有限元模型，分析了构件上各受荷区域的荷载均布情况，验证了所提出装置的可行性。基于所建立的简化分析模型，讨论了刚度、荷载作用时间和阻尼比等关键参数，对构件上动荷载的影响，并给出了基本规律。计算结果表明，在选定合适的参数后，落锤冲击气囊试验装置能够较好地模拟构件上的均布动荷载，为实验室中进行爆炸冲击试验提供了一种可能。

关键词：冲击；气囊；均布动荷载；爆炸荷载

中图分类号：TU317

文献标志码：A文章编号：16744764（2016）01012207

Abstract：

It is very difficult to apply a uniform impact （explosion） loads on the structural member in the laboratory. We presented an innovative loading apparatus using a drop hammer impacting on the airbag to simulate the uniform loads，and proposed a twoDOF （degreeoffreedom） analytical model. The fundamental mechanical equation， boundary conditions and solution method were given and verified. A fine finite element model was established and the loading effects on the divided area were compared to verify the feasibility of presented apparatus. Some key factors of the dynamic loads were also discussed， including stiffness， load duration and damping ratio.The results show that the presented apparatus can be applied to simulate the uniform dynamic loads with choosing parameters，which provides a possibility for experiments simulating blast loading in the laboratory.

Keywords：impact；airbag；uniform dynamic loads； explosion load

燃气爆炸、恐怖袭击等偶然性爆炸灾害频繁发生，常常造成工程结构损伤破坏，乃至坍塌。因此，结构构件在爆炸荷载作用下的动力响应一直是学术界的研究热点[14]。在常见比例爆距范围，通常认为构件上的爆炸荷载是均布动荷载[1]，而野外爆炸试验又有难度大、费用高的不足，所以人们一直致力于寻找一种在实验室范围内模拟结构上均布动荷载（爆炸）的试验方法。

目前，实验室模拟冲击（爆炸）动荷载的装置主要有：模爆器、爆坑、轻气炮，SHPB和落锤等。模爆器和爆坑试验费用高、危险性大、试验条件苛刻，所需炸药不易获得，且荷载持续时间长于实际情况。落锤虽然在荷载作用时间和量级上接近爆炸，但是无法实现均布加载。SHPB和轻气炮的试验只能针对材料，无法对结构构件施加荷载。由于气囊刚度低，把气囊作为分配力媒介能够很好地解决构件变形后的荷载均布问题。Wang等[5]将作动器作用于钢板约束气囊，实现了水箱上壁的均布静力加载；陈帅等[6]通过理论分析，证明了气囊能较好地对机翼施加均布荷载；王琦等[78]则通过液压油缸及柔性橡胶研制出一种柔性均布压力加载装置，真实模拟了深部岩土应力边界条件。课题组[9]在前期工作中亦研制了一种适用于梁、柱构件的气囊均布加载拟静力试验装置。Chen等[10]采用摆锤撞击约束气囊，首次实现了对双层多拱板施加均布动荷载，但是其对加载效果没有进行评估。

综上所述，气囊加载方法目前主要应用于结构构件的拟静力加载，对于冲击动载下的研究较少，其适用条件和关键影响因素也尚待进一步研究。本文结合落锤试验和气囊加载技术，提出了一种落锤冲击气囊模拟结构上均布动荷载的试验装置，给出了相应的简化理论分析模型，系统验证了冲击气囊均布加载装置的可行性。

1试验装置和方法

本文建立的落锤冲击气囊试验装置主要由落锤冲击加载系统、钢质分配梁、气囊、支座和PC板气囊挡板组成。如图1所示，该加载系统通过落锤冲击产生压力作用于大刚度钢质分配梁上，继而钢质分配梁传递力给气囊，气囊再将力均布分配与梁试件上，从而实现均布加载的目的。具体传力路径为：落锤钢质分配梁气囊试件。系统采用透明高强PC板作为气囊均布加载系统的围护装置，能够实现构件破坏形态的实时观测并保证安全。

2双自由度弹簧阻尼冲击加载分析模型

2.1理论分析模型

本试验装置中，由于钢质分配梁抗弯刚度较大，冲击过程中可近似认为整体无变形，各点位移随时间同时变化，可用x（t）来表示。根据爆炸荷载作用下梁的等效单自由度SDOF（single degreeoffreedom）模型理论[11]，由于气囊底部能够在试验过程中贴合试验梁上表面发生协调变形，因此，也可将气囊同样简化为一个等效单自由度来进行动力响应分析，参考点选为梁上表面跨中点，用位移y（t）表示。因此，落锤冲击气囊施加均布动荷载试验装置可简化为一个双自由度弹簧阻尼冲击加载分析模型，如图2所示，其中m、k、c分别表示质量、刚度和阻尼，x（t）和y（t）分别表示质体1和质体2的位移，F（t）为落锤荷载，下标1、2分别表示钢质分配梁和气囊。

从图4可以看出，基本力学方程求出的解与ABAQUS结果完全重合，证明了基本力学方程和求解方法的正确性。

3关键影响参数讨论

本文提出落锤冲击气囊试验装置的目的是尽可能模拟冲击、爆炸等均布动荷载。对于双自由度弹簧阻尼冲击加载分析模型来说，对构件施加的荷载就是简化分析模型的支反力RF（t），由于爆炸冲击荷载作用时间短，升压时间快，荷载峰值大，因此，亟需对支反力RF（t）的荷载特点进行评估。冲击荷载持续时间很短，一般情况下阻尼力还来不及吸收太多能量结构就达到了最大变形[15]，所以本文主要考虑冲击荷载下体系的无阻尼反应，而对于模拟冲击动荷载，仅仅需峰值响应时，例如最大荷载峰值时，可以忽略高阶振型的影响，只考虑低阶阵型。本文选取最低频率w1对应的第一阵型的周期作为双自由度弹簧阻尼冲击加载分析模型的周期T。

双自由度弹簧阻尼冲击加载分析模型可以通过质量和刚度正交条件解耦为两个独立的单自由度方程，并可分别求出模型内两质体的位移x（t）、y（t），而支反力RF（t）=k2·y（t），即支反力只与气囊位移y（t）以及其刚度k2有关，所以，支反力问题是个单自由度相关问题。

落锤荷载可等效为三角形脉冲荷载，如图5所示[16]。可表示为

F（t）=P0（1-tt1）0≤t≤t10t≥t1 （3）

其中：t为时间，t1为荷载作用时间，P0为三角形脉冲荷载的初始大小。

模拟爆炸冲击均布动荷载主要关注两个方面，一是动荷载升压时间及动荷载全周期时间，二是动荷载的峰值，以下将对这两个方面进行探讨，对于双自由度弹簧阻尼冲击系统的有限元模型，动荷载即为支反力RF（t），故可以通过研究支反力来得出影响冲击动荷载的影响因素。

1）动荷载升压时间和荷载全周期的影响因素：

由式（6）可知，结构的固有频率w为影响支反力周期的主要因素。本文取第一阵型对应的最小频率Wmin为固有频率，由W2=km，可以通过增大结构刚度或减少结构质量达到使w增大的效果使支反力的升压时间和荷载全周期时间减小。

2）支反力峰值的影响因素：

由式（7b）可知，支反力的影响因素主要为结构的固有频率W、荷载作用时间t1以及脉冲荷载的初始大小P0。可以通过增大结构刚度或减少结构质量达到使W增大的效果、通过增大荷载作用时间t1以及增大脉冲荷载的初始大小P0使支反力的峰值增大。

通过双自由度弹簧阻尼冲击加载分析模型，基于ABAQUS有限元软件Spring/Dashpot功能来验证影响支反力的关键参数。

从图6可以看出，荷载作用时间越长，升压时间变化不大，且支反力峰值越高；从图7可以看出，气囊刚度在1.2×107～1.2×109N/m范围内，刚度越大，升压时间越短，峰值越大；而当气囊刚度处于1.2×105～1.2×107N/m范围时，升压时间大体相同，峰值变化大，但是不明显，因为与钢质分配梁的刚度相比，气囊刚度太小，对支反力的升压时间和峰值影响很小。

同时，本文还考虑了气囊阻尼比对支反力的影响，从图8可以看出，阻尼比主要影响荷载的峰值，阻尼比越大，荷载峰值越小，但是阻尼比太小时支反力曲线会出现明显的震荡线性，这对模拟爆炸冲击均布动荷载是不利的。

综上所述，比较发现，升压时间t1取2 ms、气囊刚度取1.2×108 N/m、阻尼比取30%时能较好地兼顾升压时间和峰值两个影响因素，且能全周期震荡较少，与爆炸荷载模拟较好。

4荷载均布情况分析

4.1有限元模型的建立

采用气囊加载的目的是施加结构所需的均布动荷载，因此，必须对这种加载方式产生的荷载均布情况进行评估。基于ABAQUS软件平台建立了该试验装置的精细化有限元模型，如图9所示。

气囊内部设置为空腔，通用气体常数（Universal gas constant）为8.314。采用关键字*Fluid cavity描述空腔，特定环境压力（Specify ambient pressure）为101.36 kPa，理想气体分子量（Ideal gas molecular weight）为0.044。通过在气囊四周施加法向位移约束来模拟气囊四周的PC板围护。为准确模拟梁构件的简支约束条件，分别在两端设置小垫块。采用通用接触关键字*General Contact来模拟重锤与钢质分配梁、钢质分配梁与气囊、构件与支座垫块之间的接触；采用通用关键字*Tie来模拟气囊与构件的贴合。

为简化计算，本文近似将锤头放置于大刚度钢质分配梁顶部，并给定个初速度v0来模拟落锤下落冲击过程；通过ABAQUS/Explicit求解器求解。气囊壳体采用四节点缩减积分壳单元S4R离散，其余组成部分采用三维八节点缩减积分单元C3D8R离散，单元尺寸大小均为50 mm。

4.2计算结果验证

采用气囊加载的目的是施加结构所需的均布动荷载。因此，将受荷构件上表面划分为8个区域，通过比较各区域荷载时程曲线来验证荷载均布情况，如图10所示。根据对称性，图11为落锤冲击初速度15 m/s时，试验梁上表面区域1～4上的荷载时程曲线。图12为2 ms时刻，试件上的Mises应力分布云图。

从图11可以看出，区域2～4上的荷载时程曲线吻合较好，说明区域2～4荷载均布情况较好。区域1的荷载与区域2～4相差较大，其主要原因是，随着落锤撞击力的传递，区域2～4所在的试验梁会发生挠曲，而区域1受支座限制，不能变形，所以受荷较其余区域大得多。模拟结果表明，支座之间区域的所受动荷载是均匀的，本文提出的落锤冲击气囊试验装置能较好地对构件施加均布动荷载。

5结论

本文提出了一种落锤冲击气囊施加均布动荷载的试验装置，建立了相应的双自由度弹簧阻尼冲击加载理论分析模型，给出了基本动力方程和求解方法，并进行了有限元验证。讨论了刚度、荷载作用时间和阻尼比等关键参数对构件上动荷载时程曲线的影响。建立了落锤冲击气囊施加均布动荷载试验装置的精细化有限元模型，对构件表面受荷均布情况进行了讨论。

1）等效双自由度弹簧阻尼冲击加载分析模型的计算结果表明，刚度、荷载作用时间和阻尼比等是影响构件上动荷载时程曲线的关键参数。在一定范围内，支反力升压时间随刚度的增大而缩短，而其峰值变大。荷载作用时间越长，升压时间越长，且支反力峰值越高。阻尼比越大，支反力峰值越小。

2）落锤冲击气囊试验装置能较好地对支座内试件区域施加均布动荷载，但在支座附近的荷载均布效果不理想。

3）分析给出了落锤冲击气囊施加均布动荷载试验装置最佳材料参数，为实验室模拟构件爆炸冲击试验提供了一种可能。

参考文献：

[1] 方秦， 柳锦春， 张亚栋，等. 爆炸荷载作用下钢筋混凝土梁破坏形态有限元分析[J]. 工程力学，2001，18（2）：18.

Fang Q， Liu J C， Zhang Y D， et al. Finite element analysis of failure modes of blastloaded R/C beams [J]. Engineering Mechanics， 2001，18（2）： 18.（in Chinese）

[2] 师燕超，李忠献.爆炸荷载作用下钢筋混凝土柱的动力响应与破坏模式[J]. 建筑结构学报，2008，29（4）：112117.

Shi Y C， Li Z X. Dynamic responses and failure modes of RC columns under blast load [J]. Journal of Building Structures， 2008，29（4）： 112117.（in Chinese）

[3] 阎石，张亮，王丹，等.钢筋混凝土板在爆炸荷载作用下的破坏模式分析[J]. 沈阳建筑大学学报，2005，21（3）：177180.

Yan S， Zhang L， Wang D，et al. Failure mode analysis for RC slab under explosive loads [J].Journal of Shenyang Jianzhu University， 2005，21（3）： 177180.（in Chinese）

[4] 王礼立. 爆炸与冲击荷载下结构和材料动态响应研究的新进展[J]. 爆炸与冲击，2001，21（2）：8188.

Wang L L. Progress in studies on dynamic response of structures and materials under explosive/impact load [J]. Explosive and Shock Waves， 2001，21（2）： 8188.（in Chinese）

[5] Wang Y H， Liew J Y R ， Lee S C. Performance of water tank under static and dynamic pressure loads[C]∥ 6th International Conference on Protection of Structures against Hazard， Tianjin， 2014： 588597.

[6] 陈帅，杨智春，李斌.结构试验中均布载荷的气囊加载系统设计方法[J]. 工程力学，2012，29（6）：146150.

Chen S，Yang Z C，Li B. Design method of uniform loading system using airbags in structure test [J]. Engineering Mechanics， 2012，29（6）： 146150.（in Chinese）

[7] 王琦，王汉鹏， 李术才，等. 柔性均布压力加载装置的研制及试验分析[J]. 岩石力学与工程学报，2012，31（1）：133139.

Wang Q， Wang H P， Li S C， et al. Development and test analysis of flexible uniform pressure loading device [J]. Chinese Journal of Rock Mechanics and Engineering，2012，31（1）：133139.（in Chinese）

[8] 李术才，王琦，李为腾， 等.柔性均布压力加载装置在模型试验中的应用研究[J]. 岩土力学，2014，35（1）：6166.

Li S C， Wang Q， Li W T， et al. Research on application of flexible uniform pressure loading device to model test [J]. Rock and Soil Mechanics， 2014，35（1）：6166.（in Chinese）

[9] 中国人民解放军理工大学. 一种拟静力试验装置[P].中国，204679269U. 20150930.

PLA University of science and technology. A kind of quasistatic loading apparatus [P]. China， 20204679269U. 20150930.（in Chinese）

[10] Chen W S. Experimental and numerical studies of novel protective panels against blast and impact loadings[D]. Perth： The University of Western Australia， 2013.

[11] 方秦， 柳锦春. 地下防护结构[M]. 北京：中国水利水电出版社，2010.

Fang Q， Liu J C. Underground defensive project[M]. Beijing： China Water & Power Press， 2010. （in Chinese）

[12] 拉夫 R，彭津 J. 结构动力学[M].王光远， 译 .北京：高等教育出版社， 2013.

Ray C J. Structural dynamics [M].Beijing： Higher Education Press，2013. （in Chinese）

[13] Abaqus Inc.Abaqus Example Problems Manual[M].USA：Abaqus Inc， 2013.

[14] 潘腾，陈力，方秦，等. 钢筋混凝土梁气囊拟静力加载数值试验分析[J]. 实验室研究与探索，2015.

Pan T， Chen L， Fang Q， et al. Numerical testing analysis of quasistatic loading using airbag on reinforced concrete beam [J].Research and Exploration in Laboratory，2015.（in Chinese）

[15] 相恒波， 方秦，王玮， 等.磁流变阻尼器抗爆隔震性能的数值模拟[J]. 振动与冲击，2007，26（11）：106111.

Xiang H B， Fang Q， Wang W， et al. Numerical simulation on shock isolation of blast resistant structure with magneto rheological dampers [J]. Journal of Vibration and Shock， 2007，26（11）：106111.（in Chinese）

[16] 喻忠操 .基于落锤试验机的砌体填充墙抗爆加固试验研究[D]. 长沙：湖南大学，2013.

Yu Z C. Study on blast resistance of reinforced masonry infill wall based on impact of drop hammer machine [D].Changsha： Hunan University，2013.（in Chinese）

（编辑胡玲）