1、第 卷 第 期 年 月福建工程学院学报 ：长柱侧向冲击响应及刚塑性分析潘爽，郭金龙，郭晓（福建农林大学 交通与土木工程学院，福建 福州；福建江夏学院 工程学院，福建 福州；中建海峡建设发展有限公司，福建 福州）摘要：采用有限元软件对圆中空夹层钢管混凝土（）长柱构件进行侧向冲击数值模拟，揭示了 长柱构件冲击响应全过程工作机理，探究了空心率、长细比、冲击能量等参数对构件动力响应的影响。结果表明：长柱构件的冲击力发展过程大致分为接触、震荡和下降 个阶段。动力冲击抗弯承载力提高系数随空心率、长细比的增大而降低，随冲击能量的增大而提高。基于刚塑性理论对构件冲击过程进行简化处理，提出了 长柱构件冲击抗弯承
2、载力提高系数和跨中最大挠度计算公式。关键词：长柱；冲击响应；刚塑性分析；冲击抗弯承载力；跨中最大挠度中图分类号：文献标志码：文章编号：（），（，；，；，）：（），：，：；收稿日期：基金项目：福建省自然科学基金项目（）；福建江夏学院国家基金培育项目（）第一作者简介：潘爽（），男，湖北麻城人，硕士研究生，研究方向：钢混凝土组合结构。通信作者：郭金龙（），男，福建福州人，副教授，博士，研究方向：钢混凝土组合结构。近年来，中空夹层钢管混凝土构件在输电塔、桥梁等工程中得到了越来越多的应用，不少学者开展了该类构件受弯、耐火、抗震、压扭等方面的实验和理论研究。但实际工程中冲击事故频发，因此有必要对工程中的撞
3、击作用及响应机理进行分析。王潇宇等对钢管混凝土试件进行冲击试验，发现边界条件不同会影响构件的破坏模式。王路明等提出了冲击后构件的损伤评估方法。王丙斌等对外包不锈钢中空夹层钢管混凝土构件进行冲击试验，发现空心率超过一定范围时，其抗冲击性能呈下降趋势。史艳莉等试验发现空心第 期潘爽，等：长柱侧向冲击响应及刚塑性分析率为的方套圆中空夹层钢管混凝土构件在冲击作用下主要发生局部凹陷变形。综合上述，现有研究中针对圆中空夹层钢管混凝土（以下简称）长柱构件的抗冲击性能研究较少，已有研究多采用试验或数值模拟分析法，未见研究者利用刚塑性理论探究侧向冲击下 长柱构件的抗冲击设计方法。本文利用 软件针对 的 长柱构件
4、建立有限元分析模型，模拟侧向冲击试验，并针对该类构件基于刚塑性理论建立简化分析模型。有限元模型建立 材料模型（）混凝土本构模型核心混凝土强度等级为，其中材料单轴受压、受拉应力应变关系采用钢管对混凝土约束的本构模型。混凝土材料属性设置采用塑性损伤模型，参数取值参考文献为：泊松比 取；密度 取 ；流动势偏心率 取；双轴极限抗压强度与单轴极限抗压强度之比取；拉伸子午面和压缩子午面上的第二应力不变量的比值 取 ；黏性系数 取；膨胀角 取，混凝土弹性模量按 计算。在冲击作用下，核心混凝土的受压动力强度放大值采用式（）进行计算：，（）式中，为应变率 时的混凝土动力抗压强度，；为混凝土静力抗压强度，；为静态
5、受压应变率值，取值 ；（），取值 。混凝土受拉强度放大值的计算公式如式（）所示：，（）式中，为应变率 时的混凝土动力抗拉强度，；为混凝土静力抗拉强度，；为静态受拉应变率值，通常取；（），取值 。（）钢材本构模型钢材单轴应力（）应变（）关系采用五阶段模型。参数如下：密度 ，弹性模量 ，泊 松 比 ，屈 服 强 度 为。采用 模型考虑钢材在冲击作用下的应变率效应，如式（）所示：（）（）式中，为钢材在应变率 时的屈服强度，；为拉伸实验测得的钢材屈服强度，；和 为材料参数，取 ，取。单元类型、网格划分和界面接触图 为有限元模型及截面示意图，采用三维实体单元 创建混凝土和钢管部件；将冲击体简化为刚性体，
6、运用刚体壳单元模拟。在接触属性设置中，钢管和混凝土、冲击体和钢管的相互作用关系均可处理为表面与表面接触，其中钢管与混凝土切向摩擦因数取值；冲击体和钢管切线方向摩擦因数为，两类接触的法线方向均采用硬接触。落锤置于构件跨中正上方 ，指定冲击体初始速度，整个模型施加重力场，采用显式动力学分析步进行计算。图 有限元模型及截面示意图 有限元模型验证利用文献中的 落锤试验进行对比验证，试件、探究了不同冲击能量下的冲击响应，数字、表示冲击速度分别为、。图 为试验与有限元模拟结果对比曲线，可见两组曲线在持续时段、整体发展趋势上均较为接近，而冲击力、挠度峰值的误差福建工程学院学报第 卷均在 以内，两组曲线吻合度
7、较高。图 为文献中试件试验与有限元的破坏模态对比，试验中随着冲击能量的增大，构件整体弯曲变形增大，跨中冲击部位凹陷增大。构件 整体弯曲和凹陷均不明显，有明显的弯曲变形，的整体变形和局部凹陷最为明显，原因是冲击能量的增大使得构件需要吸收更多的塑形应变能，从而发生较大变形。模拟中构件跨中同样出现局部凹陷，整体呈弯曲变形，随着冲击能量的增大，破坏现象越明显，有限元变形结果与试验破坏现象基本一致。综上，根据以上参数设置所建立的有限元模型能对 构件侧向冲击进行较为准确的模拟，可用于本文后续对 长柱构件的侧向冲击模拟。图 试验与有限元结果对比 图 有限元与文献的构件整体破坏模态对比 动力响应机理分析根据上
8、述 参数设置，建立 长柱构件有限元模型，初始模型基本参数为：外钢管直径 为 ，内钢管直径 为 ，外钢管壁厚 和内钢管壁厚 均为 ，构件长度 为 ，冲击速度 为 ，落锤质量 为 ，长细比 为，空心率为。通过改变内钢管直径、构件长度、冲击速度、冲击质量等条件，共设计了 个构件，分别探究空心率、长细比、冲击能量对侧向冲击作用下两端固定约束的 长柱构件动力响应的影响。以初始模型为典型构件，对构件冲击全过程、破坏模态及跨中截面应力应变发展进行分析。冲击全过程分析图 为典型构件的各时程曲线，包括冲击力、挠度、构件速度、落锤速度，根据冲击力时程曲线，将整个冲击过程分为 个阶段，即接触段（）、震荡段（）、下降
9、段（）。接触段（）：落锤速度减小，构件速度急速增大且构件的冲击力迅速增大，时刻达到峰值，此时构件和落锤的速度几乎一致，随后构件与落锤一同减速运动，冲击力又迅速下降。震荡段（）：构件和落锤继续做减速运动，此时冲击力较峰值大幅降低，点时，构件和落锤的速度降低为，此时构件跨中挠度达到最大，冲击力开始进入下降段。下降段（）：冲击力迅速下降，跨中挠度有少许回弹，构件和落锤开始反向加速。冲击全过程在 内完成。破坏模态分析图 为典型构件冲击过程应力云图，时刻，构件没有明显变形，冲击部位迅速形成应力集中。时刻，构件由于向下运动与支座处拉扯，支座处也形成应力集中，此时构件有轻微弯曲，随着构件进一步向下运动。时刻
10、，冲击部位和支座处发第 期潘爽，等：长柱侧向冲击响应及刚塑性分析图 冲击力、跨中挠度和速度的时程曲线 ，生塑形变形，冲击部位有明显凹陷，支座下部也有鼓屈现象，过程中，应力不断增大。时刻，应力伴随构件的向上回弹有些许下降，最终，应力主要集中在跨中和两端支座处，构件有明显弯曲。图 构件冲击过程应力云图 应力应变发展分析图 为跨中各点位的应力应变曲线，、点位表示外钢管各处的标记点，、点位表示内钢管各处的标记点，观察发现应力应变曲线总体有以下特点：各点位的应力应变曲线大致经历了弹性、弹塑性、强化和卸载 个阶段，弹性阶段应力随应变的增大基本呈线性增加，弹塑性阶段应力水平保持平稳，强化阶段应力值随应变的增
11、大继续增大，卸载阶段应力值出现陡然下降。靠近冲击部位的 点位和 点位受压应变值为负，冲击部位侧向 点位和 点位及对向 点位和 点位受拉应力值为正，各点位的应力值最终都超过了钢材的屈服应力，由于钢管与混凝土良好的协同作用，钢材充分发挥了其材料强度。图 跨中截面应力应变曲线 动力响应参数研究 空心率在初始模型的基础上，通过改变内钢管直径设计了、的 个构件，空心率由 提高至 过程中，构件冲击力峰值明显降低，平台值反而逐步提高，与 时相比，、时，冲击力峰值分别下降、，平台值分别增大、，该结果反映出空心率的增大削减了构件的接触刚度，但可以提高构件整体抗弯性能。构件跨中挠度随着空心率的增大而减小，主要是因
12、为大空心率构件发生局部凹陷破坏。长细比保持 的截面尺寸不变，改变构件长度设计了、的构件，长细比的变化对构件冲击力峰值无影响，但冲击力平台值随着长细比的增大而大幅降低，同时跨中最大挠度福建工程学院学报第 卷有所增大，与 的构件相比，冲击力平台值和跨中最大挠度的最大变化幅度分别高达、，这是因为长细比的增大使得支座对构件的约束能力降低。冲击能量设计了两组模拟试验，从冲击速度和冲击质量两个方面进行分析，冲击质量分别为、，冲击速度分别为、，冲击能量 采用式（）计算（）随着冲击速度的增大，构件的冲击力峰值也逐渐增大，这可能是因为材料应变率效应的影响，相比于静力作用，动态作用下材料强度有一定提高。冲击质量的
13、变化对冲击力无显著影响。两组模拟中冲击持续时间和跨中最大挠度均随冲击能量的增大而增大，且跨中最大挠度与冲击能量呈正相关。刚塑性理论简化分析 刚塑性简化模型两端固定约束边界条件下，构件在冲击荷载作用下通常会在跨中和两端产生 个“塑性铰”，构件主要表现为对称式的跨中受弯破坏，内钢管与外钢管破坏模态极为相似，核心混凝土由内外两层钢管包裹，没有发生严重的断裂或破碎。由破坏模态可知，冲击响应过程中，构件支座端和跨中部位出现应力集中现象，且主要发生塑性变形，刚塑性理论认为大部分冲击能量主要转化为构件“塑性铰”区域的塑性应变能。根据等效塑性铰方法将冲击破坏模型简化为塑性铰模型，如图 所示，为塑性铰的旋转角度
14、（）。本模型为两端固支约束，故认为冲击破坏沿构件长度方向呈对称分布，所以，根据几何关系进一步将 简化，为净跨长度，采用构件跨中最大挠度，相应的假定构件截面两端和跨中转角处的动力塑性弯矩 完全相等，各塑性铰处的耗能可以表示为每处的动态塑性弯矩 与转角 的乘积，构件整体塑性变形耗能 即为 个转角处耗能之和，如式（）所示，构件整体变形消耗的塑性变形耗能 可由软件提取，采用跨中挠度变形量为 时对应的消耗值。图 圆中空夹层钢管混凝土塑性铰简化模型 （）进而动态塑形弯矩 可表示为：（）参考文献中 柱受弯破坏时的静态极限弯矩 表达式为：，（）其中，为外钢管和夹层混凝土的组合抗弯承载力，；，为内钢管抗弯承载力
15、，；如式（）（）。，（），（）式中，为外钢管和夹层混凝土的组合毛截面模量，（）（）；为内钢管的毛截面模量，（）（）；、为截面抗弯承载力系数，对于圆形截面，（）（），为名义约束效应系数标准值，。动态弯矩 与静态极限弯矩 存在一定比例关系，为冲击抗弯强度提高系数：（）实用计算公式如图 所示为不同参数对冲击抗弯强度提高系数、跨中最大挠度 的影响，其中空心率、长细比与冲击抗弯强度提高系数 呈线性负相第 期潘爽，等：长柱侧向冲击响应及刚塑性分析关，冲击能量与冲击抗弯强度提高系数 呈正相关，在本文所研究参数范围内，取值 ，跨中最大挠度 随长细比、冲击能量的增大而增大，随空心率的增大而减小。图 不同参数对、
16、的影响 以上参数是影响冲击抗弯强度提高系数、跨中最大挠度 的关键性因素，考虑空心率（）、长细比（）、冲击能量（）对构件冲击性能的影响，在为、为、为 的范围内，通过各参数结果回归分析，得到 的实用计算公式：（）（）（）（）式中，（），（），（）。的实用计算公式：（）（）（）（）式中，（），（），（）。图 为有限元计算与公式计算对比，其相对误差均在 以内，公式精度较高。图 有限元计算与公式计算对比 结论）长柱构件的冲击力发展过程大致分为接触、震荡、下降 个阶段，曲线几乎没有明显的平台段。）动力冲击抗弯承载力提高系数 与空心率、长细比呈线性负相关，与冲击能量呈正相关，在本文所研究参数范围内，取值；跨中最大挠度随长细比、冲击能量的增大而增大，随空心率的增大而减小。）利用刚塑性简化模型对冲击响应进行分析，发现冲击能量主要转化为构件冲击部位和两端支座处“塑性铰”消耗的塑性应变能，提出了 长柱构件冲击抗弯承载力提高系数 和跨中最大挠度 的计算公式。福建工程学院学报第 卷参考文献：谌扬宇，宁寄慧，张永鑫，等 中空夹层圆钢管混凝土短柱的压弯承载力计算法 建筑钢结构进展，（）：熊明祥，胡琪东，刘博元，等
