总结了结构在冲击荷载作用下的4种响应模式，并分析了各响应模式下的动力响应（冲击力、节点速度、节点位移、杆件应力）特点。通过改变冲击点位置以及下部柱高度等参数，分析这些参数对带下部钢管柱的单层球面网壳结构响应模式的影响，揭示了结构响应随各参数的分布规律。研究结果表明：水平侧向冲击作用下结构的冲击力主要为等腰三角形脉冲荷载；冲击柱顶点时，上部结构有杆件失效，结构的节点水平位移和杆件应力最大；冲击柱高度较低时，冲击力峰值较小，节点最大水平位移和杆件最大应力较大。
关键词：冲击荷载；网壳；钢管柱；响应模式；动力响应
中图分类号：TU312文献标志码：A
0引言
冲击荷载的形式很多，如高空坠物、飞射物撞击、爆炸、山洪泥石流等，虽然其属于偶然荷载，但在当今的特殊环境与形势下，由于战争、恐怖袭击以及施工中的失误等原因，发生的概率却在不断增加。9·11事件发生后，冲击荷载对结构的作用已经引起学者们的广泛关注。大跨空间结构多用于大型公共建筑、工业建筑、地标性建筑等重要的建筑工程领域中，此类结构如果遭到袭击出现严重破坏，必将带来生命及财产的重大损失。
Karagiozova等[15]主要对冲击荷载作用下的圆柱壳构件进行了系统的研究，分析了圆柱壳的冲击屈曲特性。Nagel等[6]对锥形截面的薄壁方钢管在斜向冲击荷载作用下的动力性能进行了研究。陆新征等[7]对飞机撞击世贸中心倒塌过程进行了仿真分析。刘锋等[8]对冲击荷载作用下框架结构的非线性动力响应进行了研究。目前关于冲击荷载作用下构件材料和框架的动力响应研究比较成熟，但对于大跨网壳结构方面的抗冲击研究较少，主要有：李海旺等[910]和史俊亮[11]他对40 m跨度K8型单层球面网壳结构在冲击荷载下的动力响应进行了研究；范峰等[1214]对60 m跨度K8型单层网壳在冲击荷载作用下的动力响应、失效模式和失效机理进行了研究；王秀丽等[1516]对基于不同加载点的冲击荷载作用下K6型单层球面网壳结构动力响应进行了研究。以上有关网壳的文献主要是关于对网壳上部结构的冲击，而对下部支撑结构受冲击的网壳研究则尚未开展。因此，本文研究了冲击荷载作用于下部结构时网壳动力响应及破坏类型。
1有限元分析模型
本文在有限元软件ANSYS/LSDYNA中建立跨度为70 m，矢跨比为0.2，环向分割频数为9的K8型带下部钢管柱的单层球面网壳结构数值分析模型[1718]。在网壳的8个主肋方向及2个主肋之间分别设置钢管柱与上部网壳（环梁）刚接，与地面完全固接，下部柱柱高为10 m。网壳模型见图1，其图1网壳模型
采用三节点梁单元Beam161，屋面荷载通过质量单元Mass166施加，冲击物采用八节点六面体实体单元Solid164，材料模型采用刚性体模型，下部柱及上部杆件所采用的材料模型为分段线性塑性模型，是ANSYS/LSDYNA中专门适用于钢材的的模型[1214]，考虑了应变率效应，所谓应变率效应就是指应变率提高后，材料内部发生复杂变化，导致其应力应变关系更为复杂，会使材料产生应力硬化现象。该模型可以自定义有效真应力与有效塑性应变的关系，以及应变率对屈服应力σy影响的比例因子，材料模型如下
2.1响应模式
带下部钢管柱的单层球面网壳结构当下部柱受冲击荷载作用时，其动力响应具有一定规律。通过变化冲击物的质量和速度对下部柱进行冲击，其中通过变化冲击物密度来改变冲击质量[1214]，因此本文不考虑冲击物尺寸变化对结果的影响。经大量数值模拟算例分析，以结构动力响应及变形为判断依据，可以发现其响应模式有以下4种：①模式1（结构轻度损伤），直接受冲击的下部柱进入塑性，其他构件基本处于弹性状态，结构无明显变形；②模式2（轻微局部坍塌），直接受冲击的下部柱进入塑性阶段且有较大塑性应变，上部结构塑性进一步发展，出现轻微局部坍塌；③模式3（局部坍塌），直接受冲击的下部柱破坏失效，上部部分杆件进入塑性，出现局部坍塌，比模式2严重；④模式4（结构剪切破坏），下部柱受高速冲击被切断，上部结构虽部分进入塑性阶段但塑性应变较小，结构无明显变形。4种响应模式如图2所示，以下通过4个典型算例来说明这4种响应模式。4个算例冲击位置均为下部柱柱高为10 m，柱高1/2处，模式1冲击物质量m=1 t，速度V=5 m·s-1；模式2冲击物质量m=20 t，速度V=30 m·s-1；模式3冲击物质量m=50 t，速度V=30 m·s-1；模式4冲击物质量m=100 t，速度V=300 m·s-1。冲击荷载下结构的响应模式分布见表3。
2.2动力响应特点
2.2.1冲击力
冲击力是指结构受到冲击荷载的大小，一定程度上反映了结构的损伤程度。4种响应模式的冲击力时程曲线见图3。
从图3中可知，冲击力主要为等腰三角形脉冲荷载。存在多次碰撞现象，但主接触过程持时很短，约为2～3 ms。从图3（a），（b）中可知，冲击力一般是首次冲击峰值很大，随后逐渐减弱。图3（c）中冲击力峰值在最后才达到最大，这是由于在低速冲击下模式3中下部柱失效导致的。一般情况下冲击力峰值是首次最大，随后逐渐减弱。
2.2.2节点速度
节点速度代表了结构的动能，速度变化在一定程度上反映了能量的变化规律。选取模式2来分析，模式2中冲击区节点0和与其相邻柱节点1的水平速度变化如图4所示。
图4中节点速度在很短时间内达到极值，然后衰减至平衡位置不断波动。从整体看，节点1速度比节点0的波动大，且其幅值都随时间呈衰减趋势。说明网壳被冲击物冲击后，冲击物动能转化为网壳结构的动能及内能，随着网壳结构阻尼的影响，节点速度逐渐趋于0。
2.2.3节点位移
节点位移表示的是结构的变形，网壳主肋向节点水平位移如表4所示。
从表4中可知，所有节点的水平位移从模式1到模式3增大，从模式3到模式4急剧减小。节点水平位移从模式1到模式3增大，这是由于在低速冲击下上部结构随着下部柱的变形而变形；从模式3到模式4急剧减小，这是由于较高速度冲击下下部柱直接剪切破坏，上部结构动力响应很小。由表4还可知，节点2的水平位移最大，这是由于节点2距冲击区最近的缘故。
2.2.4杆件应力
应力的变化反映了结构的内力变化，杆件应力的大小直接反映其受冲击影响的程度。与受冲击柱相连的网壳杆件应力如表5所示。
从表5中可知，与冲击柱相连的网壳杆件应力从模式1到模式3增大，从模式3到模式4减小。R1杆件应力值最大，这是因为R1杆件与受冲击柱直接相连；同是直接相连的R4杆件应力没有R1杆件的大，这是因为R4（环梁杆件）比R1（主肋杆件）的截面尺寸大。同时，在所有下部柱中，只有直接受冲击柱的应力较大，其他柱的应力都很小。
综上所述，4种响应模式的动力响应各有特点，
按照上文介绍的方法，对带下部钢管柱支撑的网壳结构进行参数分析，考察冲击点位置、受冲击柱高度等参数对结构冲击响应的影响。参数分析方案如表6所示。
从图5中可知：冲击柱中点（h/2）时，冲击力峰值比冲击其他位置时大；冲击柱顶点（h）时，冲击力峰值最小，这是由于冲击柱顶点时，该节点周围的杆件存在失效的情况而消耗较多能量，传递给其余杆件的能量较少，因此冲击力峰值较小。
从图6中可知：节点最大水平位移随着冲击点位置的升高而增大，从h/3到2h/3增幅较缓，从2h/3到h（顶点）增幅较陡；冲击柱顶点时，上部结构节点最大水平位移最大，这是由于柱顶点连接了柱以及2根环梁和1根主肋，相对刚度较柔，远离下部柱固接端。
从图7中可知：虽然冲击柱在2h/3时杆件最大应力较冲击柱中点时有所下降，但上部结构杆件最大应力总体趋势是随着冲击点位置的升高而升高；当冲击柱顶点时，杆件应力较大。
3.2冲击下部柱柱高变化对网壳结构动力响应的影响
选取质量10 t、冲击速度20 m·s-1的冲击物，对带柱高为6，8，10，12，14 m下部钢管柱的单层球面网壳结构在距柱底3 m处进行冲击。
带下部钢管柱的单层球面网壳结构的动力响应随冲击柱高不同的变化规律见图8～10。
从图8中可知，冲击柱高从6 m到12 m时，冲击力峰值呈现增长趋势，冲击柱高从12 m到14 m时，冲击力峰值略微下降，其下降的原因是到达一定高度后，随着柱高度的变化，结构的整体刚度相对减小，所以冲击力峰值会有所下降。
从图9中可知，节点最大水平位移随着冲击柱高的升高变化不大，处在0.7 m左右，柱高从6 m到10 m时节点最大水平位移有所下降，从10 m到14 m时略微提高。由于所选冲击位置较低，冲击荷载较小，所以结构整体位移反应较小。
从图10中可知：虽然冲击柱高14 m时，杆件最大应力较冲击柱高12 m时有所提高，但上部结构杆件最大应力总体趋势是随着冲击点位置的升高而降低，且柱高从8 m到10 m杆件最大应力下降较大；当冲击柱高度较低时，杆件应力较大。
4结语
（1）依据结构弹塑性变化及整体结构变形，可以将钢管柱单层球面网壳结构在冲击荷载下的响应模式分为4种：模式1，结构轻微损伤；模式2，整体结构轻微局部坍塌；模式3，整体结构有杆件失效的较严重局部坍塌；模式4，下部柱直接剪切破坏。
（2）水平侧向冲击作用下，结构的冲击力主要为等腰三角形脉冲荷载，冲击持时很短，持时为2～3 ms。
（3）当冲击柱顶点时，上部结构有杆件失效，此时这个结构的节点水平位移、杆件应力最大。除顶点特殊情况外，在冲击柱h/2（柱中点）处时出现整体结构最大冲击力峰值。此外还发现，冲击位置越靠上，柱更易于发生弯曲型破坏；反之，柱更易于发生剪切型破坏。
（4）当冲击柱高度较低时，冲击力峰值较小，节点最大水平位移、杆件最大应力较大。
（5）从大量的参数分析结果中可以发现，冲击破坏主要体现在直接承受冲击的下部柱及直接与其相连的网壳杆件，由于结构有足够的约束冗余度，即使是某个柱遭到严重的冲击破坏，也不会发生结构整体倒塌的严重后果。
参考文献：
[1]KARAGIOZOVA D，JONES N.Dynamic Elasticplastic Buckling of Circular Cylindrical Shells Under Axial Impact[J].International Journal of Solids and Structures，2000，37（14）：20052034.
[2]KARAGIOZOVA D，JONES N.Influence of Stress Waves on the Dynamic Progressive and Dynamic Plastic Buckling of Cylindrical Shells[J].International Journal of Solids and Structures，2001，38（38）：67236749.
[3]KARAGIOZOVA D，JONES N.On Dynamic Buckling Phenomena in Axially Loaded Elasticplastic Cylindrical Shells[J].International Journal of Nonlinear Mechanics，2002，37（7）：12231238.
[4]KARAGIOZOVA D，ALVES M.Transition from Progressive Buckling to Global Bending of Circular Shells Under Axial Impact — Part Ⅰ：Experimental and Numerical Observations[J].International Journal of Solids and Structures，2004，41（5/6）：15651580.
[5]KARAGIOZOVA D，ALVES M.Transition from Progressive Buckling to Global Bending of Circular Shells Under Axial Impact — Part Ⅱ：Theoretical Analysis[J].International Journal of Solids and Structures，2004，41（5/6）：15811604.
[6]NAGEL G M，THAMBIRATNAM D P.Dynamic Simulation and Energy Absorption of Tapered Thinwalled Tubes Under Oblique Impact Loading[J].International Journal of Impact Engineering，2006，32（10）：15951620.
[7]陆新征，江见鲸.世界贸易中心飞机冲击后倒塌过程的仿真分析[J].土木工程学报，2001，34（6）：810.
LU Xinzheng，JIANG Jianjing.Dynamic Finite Element Simulation for the Collapse of World Trade Center[J].China Civil Engineering Journal，2001，34（6）：810.
[8]刘锋，吕西林.冲击载荷作用下框架结构的非线性动力响应[J].振动工程学报，2008，21（2）：107114.
LIU Feng，LU Xilin.Nonlinear Dynamic Responses of Impulsive Loaded Frame Structure[J].Journal of Vibration Engineering，2008，21（2）：107114.
[9]李海旺，郭可，魏剑伟，等.撞击载荷作用下单层球面网壳动力响应模型实验研究[J].爆炸与冲击，2006，26（1）：3945.
LI Haiwang，GUO Ke，WEI Jianwei，et al.The Dynamic Response of a Singlelayer Reticulated Shell to Drop Hammer Impact[J].Explosion and Shock Waves，2006，26（1）：3945.
[10]郭可.单层球面网壳在冲击载荷作用下的动力响应分析[D].太原：太原理工大学，2004.
GUO Ke.Dynamic Analysis of Singlelayer Reticulated Shells Under Impact[D].Taiyuan：Taiyuan University of Technology，2004.
[11]史俊亮.K8型单层网壳在撞击荷载作用下的动力响应研究[D].太原：太原理工大学，2005.
SHI Junliang.Dynamic Response Research on K8 Singlelayer Reticulated Dome Under Impact[D].Taiyuan：Taiyuan University of Technology，2005.
[12]范峰，王多智，支旭东，等.K8型单层球面网壳抗冲击荷载性能研究[J].工程力学，2009，26（6）：7581.
FAN Feng，WANG Duozhi，ZHI Xudong，et al.Performance for Kiewitt8 Singlelayer Reticulated Domes Subjected to Impact Load[J].Engineering Mechanics，2009，26（6）：7581.
[13]王多智，支旭东，范峰，等.冲击荷载下K8单层球面网壳的破坏类型[J].工程力学，2008，25（增2）：144149.
WANG Duozhi，ZHI Xudong，FAN Feng，et al.Failure Patterns of Kiewitt8 Singlelayer Reticulated Domes Under Impact Loads[J].Engineering Mechanics，2008，25（S2）：144149.
[14]范峰，王多智，支旭东，等.冲击荷载下凯威特型球面网壳的失效模式及其判别方法[J].土木工程学报，2010，43（5）：5662.
FAN Feng，WANG Duozhi，ZHI Xudong，et al.Failure Modes and Discrimination Method for Kiewitt Reticulated Dome Under Impact Loads[J].China Civil Engineering Journal，2010，43（5）：5662.
[15]王秀丽，马肖彤，李俊杰.基于不同加载点的冲击荷载作用下单层球面网壳结构动力响应研究[J].空间结构，2014，20（2）：2328.
WANG Xiuli，MA Xiaotong，LI Junjie.Dynamic Response of Singlelayer Reticulated Shells Under Different Points Impact Loading[J].Spatial Structures，2014，20（2）：2328.
[16]王秀丽，马肖彤.冲击荷载作用下受损网壳结构全过程动力响应分析[J].建筑科学与工程学报，2013，30（3）：1419.
WANG Xiuli，MA Xiaotong.Whole Process of Dynamic Response of Damaged Reticulated Structure Shell Under Impact Load[J].Journal of Architecture and Civil Engineering，2013，30（3）：1419.
[17]王占飞，包龙生，隋伟宁，等.外置钢管补强圆形钢桥墩非线性有限元分析[J].筑路机械与施工机械化，2008，25（12）：7679.
WANG Zhanfei，BAO Longsheng，SUI Weining，et al.Nonlinear Finite Element Analysis for Tubular Steel Pier Reinforced with Outer Steel Tube[J].Road Machinery & Construction Mechanization，2008，25（12）：7679.
[18]沈世钊，陈昕.网壳结构稳定性[M].北京：科学出版社，1999.
SHEN Shizhao，CHEN Xin.Stability of Latticed Shells Structures[M].Beijing：Science Press，1999.
推荐访问： 荷载 球面 单层 钢管 响应
