Numerical Simulation of Responses and Failure Modes of Reinforced Concrete Beams under Drop-Weight Impact Loadings
-
摘要: 钢筋混凝土在动态冲击下表现出与静态加载不同的结构响应,且破坏模式更为复杂。在钢筋三折线本构模型中引入应变率效应,利用ABAQUS显式动力分析模块,对钢筋混凝土梁在不同高度冲击下的结构响应进行了数值模拟。得到的冲击力和跨中挠度时程曲线与实验结果吻合较好,验证了模型的有效性。基于该模型,研究了配筋率分别为2.56%、2.66%和2.76%时,钢筋混凝土梁在不同冲击速度下的结构响应。结果表明:增大配筋率能够提高梁的承载能力;随冲击速度的增大,配筋率对梁抗变形能力的增强效果逐渐减弱;当冲击速度为4.85 m/s时,配筋率对梁破坏模式的影响微弱;当冲击速度大于4.85 m/s时,随配筋率的减小,破坏模式由剪切破坏转变为弯曲破坏。Abstract: Reinforced concrete structural members subjected to impact loads behave quite differently as compared to those subjected to quasi-static loading, with their failure mode becoming more complex.In this work, by introducing the strain rate effect of reinforcement in the trilinear model of the reinforcement, we simulated the structural responses of reinforced concrete beams under different impact loadings based on the dynamic analysis module of ABAQUS.The curves of impact-time and mid-point deflection-time were observed to agree well with those from the experiments.Based on this model, we simulated the responses of beams with the reinforcement ratios of 2.56%, 2.66%, and 2.76%, respectively.The comparison shows that the bearing capacity and deformation resistance of the beams increased with the increase of the reinforcement ratios; the enhancement effect of the reinforcement ratio weakens gradually as the impact velocity increases; when the impact velocity is 4.85 m/s, the reinforcement ratios have slight effect on the failure mode of beam at low impact velocities; in addition, when the impact velocity is higher than 4.85 m/s, the failure mode changes from shear failure to bending failure with the decrease of the reinforcement ratio.
-
Key words:
- drop weight impact /
- reinforcement ratio /
- strain rate effect /
- reinforced concrete
-
钢筋混凝土结构在建筑、桥梁工程中有广泛的应用。混凝土作为非弹性材料,在动态和静态加载条件下会表现出不同的力学特性,且存在明显的加载率效应。国内外许多研究表明,钢筋混凝土结构在高应变率加载条件下具有较高的变形抗力,结构的失效模式也会发生相应的变化。例如,Ožbolt等[1]将率相关微平面模型引入3D有限元模型中,分析了加载速率对锚杆拔出时混凝土失效模式的影响。Travaš等[2]模拟了素混凝土梁在准静态和动态加载下的结构响应,发现准静态和低速冲击(v=2 m/s)下梁发生弯曲破坏;当冲击速度大于2 m/s时,梁的破坏模式由弯曲破坏转为剪切破坏。付应乾等[3]对钢筋混凝土梁进行三点弯曲试验,探讨了不同冲击高度下钢筋混凝土结构的破坏机理。Saatci等[4]对钢筋混凝土梁进行不同质量锤头的落锤冲击试验,发现梁可承载的冲击力随抗剪强度的升高而增大。展婷变等[5]对钢筋混凝土梁进行落锤试验,分析了冲击速率对最大承载力和梁变形的影响。Zineddin等[6]对钢筋混凝土板进行不同高度的冲击试验,发现高速冲击下配筋率对于改善混凝土结构的层裂破坏有较大作用,但对最大承载力影响较小。范向前等[7]开展了非标准钢筋混凝土三点弯曲梁试件断裂试验,认为梁能够承载的最大载荷由配筋情况决定。沈培峰[8]研究了配筋率对钢筋混凝土梁断裂性能的影响以及钢筋在裂缝扩展过程中的限裂作用。落锤试验[9-12]是目前研究钢筋混凝土梁在低速冲击下破坏的有效手段。由于试件制备过程中浇筑技术和温度等外在条件影响,混凝土的力学性能具有一定的离散性。另一方面,实验过程中数据采集和加载条件的控制难免存在误差,因此最终实验结果的可靠性降低。
随着计算机软件的开发和有限元技术的不断成熟,数值仿真计算为实验结果的分析提供了有效的参考依据。李敏等[13]模拟了地震作用范围内加载速率对钢筋混凝土梁影响,姜华等[14]使用弹塑性断裂带模型(EPDC)模拟了钢筋混凝土梁在动态冲击下的结构响应。有限元建模存在本构模型不同、网格划分精细程度不一等差别,与实验结果进行对比后可验证其有效性。本研究基于ABAQUS有限元软件,模拟并分析了钢筋混凝土梁在不同高度冲击下的动态响应。通过模拟与实验结果对比,验证了模型的有效性,并讨论了不同冲击速度下,配筋率对钢筋混凝土梁承载力、抗变形能力和破坏模式的影响。
1. 有限元模型
付应乾等[3]利用落锤对钢筋混凝土梁进行三点弯曲试验,探讨了不同冲击高度下混凝土结构的破坏机理。本研究对该实验进行有限元建模,所用模型尺寸与参考文献[3]一致,如图 1所示。混凝土设计强度为C30,实验用梁长1.2 m,截面尺寸为120 mm×120 mm,纵向钢筋采用HRB335,箍筋采用HPB235,保护层厚度为20 mm。
模型中通过减小锤头的几何尺寸和增大密度来保持冲击能量不变。支座密度为7 800 kg/m3,弹性模量为200 GPa,泊松比为0.3。为减少计算时间,根据v=√2gh (1g=9.8 m/s2),落锤下落高度1.2、2.4和4.8 m,在预定义场中将初始速度分别设为4.85、6.86和9.60 m/s。对锤头约束在只有y方向存在位移,锤头与混凝土梁上表面接触为通用接触。支座设置为在x方向上(即混凝土梁纵向)存在位移。锤头与混凝土梁上表面初始距离为1 mm。有限元模型如图 2所示。模型中混凝土采用CDP(Crack Damage Plastic)模型,参数[15]如表 1所示。混凝土采用三维8节点缩减积分单元(C3D8R),网格尺寸为1.2 mm×1.2 mm×1.2 mm。
钢筋采用三折线模型,有限元模型中网格单元为三维8节点缩减积分单元(C3D8R),网格尺寸为1.2 mm×1.2 mm×1.2 mm,以Embed形式嵌入混凝土基体中。模型中不考虑钢筋与混凝土之间的握裹力和拉拔效应。试验表明[16],高应变率作用下建筑钢筋的弹性模量E保持不变,钢筋HPB235的弹性模量为210 GPa,HRB335和HRB400的弹性模量为200 GPa。钢筋的屈服强度和极限强度随应变率的增大而提高,其变化规律可表示为
fy,dfy,st=1+D1fy,stln˙ε˙ε0 (1) fu,dfu,st=1+D2fu,stln˙ε˙ε0 (2) 式中: ˙ε为钢筋的应变率,˙ε0=0.000 3 s-1,fy, d和fy, st分别为钢筋动态和静态加载下屈服强度,fu, d和fu, st分别为钢筋动态和静态加载下极限强度。
林峰等[16]通过实验,得到HPB235和HRB335的D1和D2分别为10.05、8.73和6.38、6.54。以三折线建筑钢筋模型为例,其应力-应变关系可表示为
σ={Eε0<ε≤εb,dfy,dεb,d<ε≤εe,dfy,d+ε−εe,dεu,d−εe,dfu,dεe,d<ε≤εu,d (3) 式中: εb, d和εe, d分别为屈服平台开始和结束时相应的应变,εb, d=fy, d/E,εe, d=εb, d+(εe, st-εb, st);εu, d为极限强度对应的应变,εu, d=εu, st。
钢筋HRB335选用˙ε=2.9 s-1时的屈服强度和极限强度,HPB235选用˙ε=6.5 s-1时的屈服强度和极限强度,见表 2。
2. 有限元模型验证
将模拟得到的冲击力、跨中挠度的最大值与实验结果[3]进行对比,如表 3所示。不同冲击速度下得到的冲击力和跨中挠度最大值与实验结果吻合较好,且略大于实验值。该现象的原因是实验中梁两端会架设钢板来保持其在冲击过程中不会产生太大位移,而有限元模型中出于简化模型没有建立这些约束条件。另一方面,下落过程中锤头与下滑轨道存在摩擦,导致锤头与混凝土梁上表面接触时的瞬时速度相比计算值较小。这两个方面最终导致实验结果比模拟结果偏小。
表 3 模拟与实验结果对比Table 3. Comparison of results between FEM and experimentMethod Drop height (m)/Mass (kg) Impact velocity/(m·s-1) Maximum impact force/kN Maximum deflection of mid-point/mm Experiment 1.2/124 104.0 28.3 FEM -/124 4.85 149.1 30.7 Experiment 2.4/124 94.7 51.2 FEM -/124 6.86 134.7 56.5 Experiment 4.8/124 173.9 90.4 FEM -/124 9.60 181.2 102.9 随着冲击速度的提高,钢筋在实验中的屈服强度和极限强度与有限元模型中的值接近,冲击力峰值误差逐渐减小。实验中通过激光位移传感器测量锤头位移,默认锤头位移与梁跨中挠度值相同,模拟中提取梁上表面中心点的位移时程曲线。实际工况中锤头位移达到最大值后锤头开始反弹,梁具有少量动能继续向下运动,跨中挠度继续增大,实验测得的锤头位移最大值比梁跨中挠度最大值偏小。随着冲击高度的增加,锤头反弹时梁具有的动能相对增大,跨中挠度增大值相对增大,导致模拟与实验的误差增大。
冲击高度为4.8 m(即v=9.60 m/s)时的冲击力模拟与实验结果对比如图 3所示。由图 3可知,模拟得到的冲击力时程曲线与实验结果变化趋势基本一致,峰值分别为181.2和173.9 kN,模拟值比实验值大4.2%,该误差在允许范围内。
图 4为不同冲击速度下跨中挠度时程曲线的模拟与实验结果对比。由图 4可知,随着冲击速度的提高,跨中挠度不断增大。当第一次冲击结束时,锤头与混凝土梁上表面分离,梁出现轻微回弹效应,跨中挠度达到极值后开始减小。由表 3可知,模拟得到的最大值比实验值偏大,分别增大8.5%、10.3%和13.8%,该结果是模拟中冲击速度比实际值偏大导致的。
由图 3和4可知,本研究模型可以很好地模拟出不同冲击速度下钢筋混凝土梁的结构响应。选取v=4.85 m/s和v=9.60 m/s时的冲击力-跨中挠度曲线与实验结果进行对比,如图 5所示。模拟与实验得到的曲线变化趋势相同,冲击力在锤头和垫块接触后迅速增大,随后冲击能量减小,一部分能量用于梁的破坏,另一部分被梁的惯性消耗。当冲击速度为4.85 m/s时,冲击力随跨中挠度的增大呈振荡发展,梁的破坏模式是弯曲失效;当冲击速度为9.60 m/s时,在跨中点形成塑性铰后,冲击力由振荡变化逐渐趋向稳定,梁的破坏模式由弯曲失效转变为局部剪切失效。
图 5中h=4.8 m时冲击力峰值的模拟值和实验值分别为173.9和181.2 kN,误差为4.20%。h=1.2 m时冲击力峰值的模拟值与实验值分别为149.1和104.0 kN,误差为43.36%。钢筋应变率随冲击高度的增大而增大,其屈服强度和极限强度均有所增大,且差别较小。有限元模型参考林峰等[16]的实验,钢筋HRB335选用˙ε=2.9 s-1时的屈服强度和极限强度,HPB235选用˙ε=6.5 s-1时的屈服强度和极限强度,见表 2。h=1.2 m时实验中钢筋的屈服强度和极限强度小于模型中的值,导致模型中梁的承载力较高,冲击力峰值误差增大。
实验中通过高速摄影记录了钢筋混凝土梁在不同冲击速度下的变形、损伤及破坏状况。图 6表示不同冲击高度下跨中挠度达到最大值时梁表面的损伤破坏情况,简图为对高速摄像记录的照片处理后的裂纹分布图,其下方分别是相对应冲击速度下梁的塑性应变云图。由图 6可知,随着冲击速度的提高,梁的挠度不断增大。冲击速度较低时,裂纹从梁下表面中心处产生且数量较多,此时梁发生弯曲破坏;冲击速度增大时,在跨中发生局部剪切破坏,梁的破坏模式发生改变。
3. 配筋率对梁承载力和破坏模式的影响
配筋率是钢筋混凝土构件中纵向受力钢筋面积与构件的有效面积之比,在桥梁工程中,一般指面积配筋率,即受拉钢筋面积与主梁面积之比。计算公式为P=As/A,其中As为纵向受力钢筋截面面积,A为主梁面积。上述实验中配筋率为2.66%,根据建筑工程钢筋配筋率规范[17]可知,混凝土强度为30 MPa,纵向受力钢筋为HRB335时,最大配筋率为2.62%。
基于本研究模型,通过改变单根钢筋直径得到配筋率为2.56%、2.66%和2.76%的梁在不同冲击速度下的冲击力时程曲线,如图 7所示,表 4为冲击力峰值结果。由图 7可知,不同配筋率下冲击力变化趋势相同,最大冲击力随着配筋率的增大而增大,且冲击力达到峰值的时间随配筋率的增大而减小。配筋率的增大,显著地提高了梁的承载能力。表 4中百分比代表配筋率为2.66%和2.76%的梁的冲击力相比配筋率为2.56%的增大百分数。当冲击速度为4.85 m/s时,以配筋率2.56%作为基准,其增大3.9%(P=2.66%)时,梁的承载力最大值增大30%;增大7.8%(P=2.76%)时,梁的承载力最大值增大49%。冲击速度为6.86 m/s时,配筋率增大3.9%和7.8%时梁的承载力最大值分别增大7.2%和60.8%。冲击速度为9.60 m/s时,配筋率增大3.9%和7.8%时梁的承载力最大值分别增大7.1%和31.3%。低速冲击下,梁的失效模式主要为弯曲破坏,承载能力随配筋率的增大显著提高。随着冲击速度的增大,梁的失效模式由弯曲破坏转向剪切破坏,钢筋对于提高混凝土梁承载力的效果逐渐减小。
表 4 冲击力峰值结果对比Table 4. Comparison of the results of maximum impact forceImpact velocity/(m·s-1) Reinforcement ratio:2.56% Reinforcement ratio:2.66% Reinforcement ratio:2.76% Max. impact force/kN Time/ms Max. impact force/kN Time/ms Max. impact force/kN Time/ms 4.85 114 0.66 149(+30.0%) 0.6 170(+49.0%) 0.5 6.86 125 0.50 134(+7.2%) 0.5 201(+60.8%) 0.4 9.60 169 0.33 181(+7.1%) 0.3 222(+31.3%) 0.3 不同配筋率梁在不同冲击速度下的跨中挠度时程曲线和具体数值如图 8和表 5所示。同一冲击速度下,梁的抗变形能力随配筋率的增大而提高,变形相对较小,跨中挠度达到最大值的时间随配筋率的增大而减小。表 5中百分数表示配筋率为2.66%和2.76%的梁跨中挠度相对配筋率为2.56%的减小百分比。冲击速度大于6.86 m/s时,配筋率对梁抗变形能力的增强效果随冲击速度的增大而逐渐减弱。
表 5 跨中位移峰值结果对比Table 5. Comparison of results of maximum mid-point deflectionImpact velocity/(m·s-1) Reinforcement ratio:2.56% Reinforcement ratio:2.66% Reinforcement ratio:2.76% Max. deflection/mm Time/ms Max. deflection/mm Time/ms Max. deflection/mm Time/ms 4.85 38.26 16.0 31.72(-17.10%) 13.32 27.70(-27.60%) 12.3 6.86 68.15 25.5 56.56(-17.00%) 17.30 49.73(-27.03%) 15.6 9.60 120.33 26.0 102.94(-14.43%) 25.00 91.37(-24.07%) 21.0 图 9为不同冲击速度下冲击力-跨中挠度曲线。由图 9可知,改变配筋率后,冲击力变化趋势一致。冲击速度为4.85 m/s时,不同配筋率梁的冲击力随跨中挠度增大呈振荡变化,梁的破坏模式为弯曲失效。冲击速度为6.86 m/s时,配筋率较小的梁冲击力在跨中挠度达到40 mm后逐渐稳定,梁的破坏模式由弯曲失效转为剪切失效,其他配筋率梁的冲击力随跨中挠度增大保持振荡变化,其发生弯曲破坏。冲击速度为9.60 m/s时,配筋率为2.56%和2.66%的梁发生剪切失效,配筋率为2.76%的梁发生弯曲失效。因此,随着配筋率的减小,梁在冲击过程中发生剪切破坏的可能性增大。
4. 结论
在钢筋的本构模型中引入应变率效应,对钢筋混凝土梁进行有限元建模。利用该模型模拟了在不同冲击高度下梁的结构响应,通过对比实验和模拟的冲击力、跨中挠度时程曲线和峰值,验证了模型的有效性。基于该模型模拟了不同配筋率钢筋混凝土梁在不同冲击速度下的结构响应,得到以下结论:
(1) 不同冲击速度下,配筋率的增大提高了梁的承载力,当冲击速度为4.85 m/s时,增大配筋率对梁承载力的增强效果更为明显;
(2) 冲击速度一定时,梁的抗变形能力随配筋率的增大而提高,当冲击速度较大时,随着冲击速度的增大,配筋率对梁抗变形能力的增强效果逐渐减小;
(3) 低速冲击下,配筋率对梁的破坏模式影响微弱,当冲击速度较大时,随着配筋率的增大,梁的破坏模式由剪切破坏转变为弯曲破坏。
-
Dilation angle/(°) Eccentricity fb0/fc0 K Viscosityparameter Density/(kg·m-3) E/GPa ν 30 0.1 1.16 0.666 7 0.000 5 2 400 26.48 0.167 Reinforced grades Diameter/mm ρ/(kg·m-3) E/GPa ν Yield strength/MPa Ultimate strength/MPa HPB235 8 7 862 210 0.3 400 540 HRB335 10 7 860 200 0.3 438 687 HRB335 12 7 850 200 0.3 438 687 表 3 模拟与实验结果对比
Table 3. Comparison of results between FEM and experiment
Method Drop height (m)/Mass (kg) Impact velocity/(m·s-1) Maximum impact force/kN Maximum deflection of mid-point/mm Experiment 1.2/124 104.0 28.3 FEM -/124 4.85 149.1 30.7 Experiment 2.4/124 94.7 51.2 FEM -/124 6.86 134.7 56.5 Experiment 4.8/124 173.9 90.4 FEM -/124 9.60 181.2 102.9 表 4 冲击力峰值结果对比
Table 4. Comparison of the results of maximum impact force
Impact velocity/(m·s-1) Reinforcement ratio:2.56% Reinforcement ratio:2.66% Reinforcement ratio:2.76% Max. impact force/kN Time/ms Max. impact force/kN Time/ms Max. impact force/kN Time/ms 4.85 114 0.66 149(+30.0%) 0.6 170(+49.0%) 0.5 6.86 125 0.50 134(+7.2%) 0.5 201(+60.8%) 0.4 9.60 169 0.33 181(+7.1%) 0.3 222(+31.3%) 0.3 表 5 跨中位移峰值结果对比
Table 5. Comparison of results of maximum mid-point deflection
Impact velocity/(m·s-1) Reinforcement ratio:2.56% Reinforcement ratio:2.66% Reinforcement ratio:2.76% Max. deflection/mm Time/ms Max. deflection/mm Time/ms Max. deflection/mm Time/ms 4.85 38.26 16.0 31.72(-17.10%) 13.32 27.70(-27.60%) 12.3 6.86 68.15 25.5 56.56(-17.00%) 17.30 49.73(-27.03%) 15.6 9.60 120.33 26.0 102.94(-14.43%) 25.00 91.37(-24.07%) 21.0 -
[1] OŽBOLT J, RAH K K, MEŠTROVIĈD.Influence of loading rate on concrete cone failure[J].International Journal of Fracture, 2006, 139(2):239-252. doi: 10.1007/s10704-006-0041-3 [2] TRAVAŠ V, OŽBOLT J, KOŽAR I.Failure of plain concrete beam at impact load:3D finite element analysis[J].International Journal of Fracture, 2009, 160(1):31-41. doi: 10.1007/s10704-009-9400-1 [3] 付应乾, 董新龙.落锤冲击下钢筋混凝土梁响应及破坏的实验研究[J].中国科学:技术科学, 2016, 46(4):400-406. doi: 10.1360/N092015-00337FU Y Q, DONG X L.An experimental study on impact response and failure behavior of reinforced concrete beam[J].Scientia Sinica:Technological, 2016, 46(4):400-406. doi: 10.1360/N092015-00337 [4] SAATCI S, VECCHIO F J.Effects of shear mechanisms on impact behavior of reinforced concrete beams[J].ACI Structural Journal, 2009, 106(1):78-86. [5] ZHAN T, WANG Z, NING J.Failure behaviors of reinforced concrete beams subjected to high impact loading[J].Engineering Failure Analysis, 2015, 56:233-243. doi: 10.1016/j.engfailanal.2015.02.006 [6] ZINEDDIN M, KRAUTHAMMER T.Dynamic response and behavior of reinforced concrete slabs under impact loading[J].International Journal of Impact Engineering, 2007, 34(9):1517-1534. doi: 10.1016/j.ijimpeng.2006.10.012 [7] 范向前, 胡少伟.不同配筋率对钢筋混凝土三点弯曲梁断裂韧度的影响[J].水电能源科学, 2013, 31(12):117-121. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=sdny201312031&dbname=CJFD&dbcode=CJFQFAN X Q, HU S W.Effects of various reinforcement ratio on fracture toughness of reinforced concrete for three-points bending beams[J].Water Resources and Power, 2013, 31(12):117-121. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=sdny201312031&dbname=CJFD&dbcode=CJFQ [8] 沈培峰.配筋率对混凝土断裂参数的影响[J].防灾减灾工程学报, 2013, 33(2):235-240. https://www.researchgate.net/profile/Shaowei_Hu/publication/272609888_Influence_of_Reinforcement_Ratios_on_Concrete_Fracture_Parameters/links/57e1083508ae52b3078c57fe.pdf?origin=publication_listSHEN P F.Influence of reinforcement ratio on concrete fracture parameters[J].Journal of Disaster Prevention and Mitigation Engineering, 2013, 33(2):235-240. https://www.researchgate.net/profile/Shaowei_Hu/publication/272609888_Influence_of_Reinforcement_Ratios_on_Concrete_Fracture_Parameters/links/57e1083508ae52b3078c57fe.pdf?origin=publication_list [9] BANTHIA N, MINDESS S, BENTUR A, et al.Impact testing of concrete using a drop-weight impact machine[J].Experimental March, 1989, 29(1):63-69. doi: 10.1007/BF02327783 [10] RAO M C, BHATTACHARYYA S K, BARAI S V.Behaviour of recycled aggregate concrete under drop weight impact load[J].Construction and Building Materials, 2011, 25(1):69-80. doi: 10.1016/j.conbuildmat.2010.06.055 [11] 任晓虎, 霍静思, 陈柏生.火灾下钢管混凝土梁落锤冲击试验研究[J].振动与冲击, 2012, 30(20):110-115. https://www.researchgate.net/profile/Jingsi_Huo2/publication/289149807_Anti-impact_behavior_of_concrete-filled_steel_tubular_beams_in_fire/links/5763f27f08ae1658e2ea2074.pdf?origin=publication_detailREN X H, HUO J S, CHEN B S.Anti-impact behavior of concrete-filled steel tubular beams in fire[J].Journal of Vibration and Shock, 2012, 30(20):110-115. https://www.researchgate.net/profile/Jingsi_Huo2/publication/289149807_Anti-impact_behavior_of_concrete-filled_steel_tubular_beams_in_fire/links/5763f27f08ae1658e2ea2074.pdf?origin=publication_detail [12] 展婷变, 宁建国, 王志华.冲击载荷下钢筋混凝土力学行为的研究[J].高压物理学报, 2016, 30(2):109-115. doi: 10.11858/gywlxb.2016.02.004ZHAN T B, NING J G, WANG Z H.Mechanical behavior of reinforced concrete under dynamic loading[J].Chinese Journal of High Pressure Physics, 2016, 30(2):109-115. doi: 10.11858/gywlxb.2016.02.004 [13] 李敏, 李宏男.钢筋混凝土梁动态试验与数值模拟[J].振动与冲击, 2015, 34(6):110-115. http://www.cnki.com.cn/Article/CJFDTotal-ZDCJ201506021.htmLI M, LI H N.Dynamic tests and numerical simulation of reinforced concrete beams[J].Journal of Vibration and Shock, 2015, 34(6):110-115. http://www.cnki.com.cn/Article/CJFDTotal-ZDCJ201506021.htm [14] JIANG H, WANG X, HE S.Numerical simulation of impact tests on reinforced concrete beams[J].Materials and Design, 2012, 39(15):111-120. https://www.sciencedirect.com/science/article/pii/S0261306912000751 [15] 刘巍, 徐明, 陈忠范.ABAQUS混凝土损伤塑性模型参数标定及验证[J].工业建筑, 2014(增刊1):167-171. http://www.docin.com/p-1535116684.htmlLIU W, XU M, CHEN Z F.Parameters calibration and verification of concrete damage plasticity model of ABAQUS[J].Industrial Construction, 2014(Suppl 1):167-171. http://www.docin.com/p-1535116684.html [16] 林峰, 顾祥林, 匡昕昕, 等.高应变率下建筑钢筋的本构模型[J].建筑材料学报, 2008, 11(1):14-20. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_jzclxb200801003LIN F, GU X L, KUANG X X, et al.Constitutive models for reinforcing steel bars under high strain rates[J].Journal of Building Materials, 2008, 11(1):14-20. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_jzclxb200801003 [17] 中国建筑科学研究院. 混凝土设计规范: GB 50010-2010[S]. 北京: 中国建筑工业出版社, 2010. -