钢筋混凝土结构是土木工程中应用最普遍的结构形式之一，随着社会的不断发展，一些建构筑物失去价值后需利用控制爆破技术将其拆除，在城市建构筑物和人口密集、倒塌空间受限、周边环境复杂等情况下，常常采用数值模拟技术对爆破效果进行预测以优化爆破设计方案。目前，LS-DYNA软件中没有一种现成的材料模型能模拟钢筋混凝土，并同时考虑钢筋混凝土在受拉、受压以及率荷载下的力学性能，随着钢筋混凝土结构高度、体积、配筋率越来越大，采用分离式模型和组合式模型需耗费大量计算资源和时间，且钢筋与混凝土相互作用比较复杂，数值模拟结果往往难以达到预期爆破效果，为钢筋混凝土结构建构筑物的安全精确爆破埋下一定隐患，因此，钢筋混凝土材料模型一直是研究的重点与难点。

谣盐：2011年3月11日日本东北部地区突发9.0级大地震后，位于本州岛福岛的核电站发生爆炸并出现核泄漏；由于外界盛传服用碘盐可以抵抗核辐射，从而引发中国大陆民众大量抢购、囤积碘盐。类似“服用碘盐可以抗辐射”及“此后一段时间内生产出来的盐将受到核污染”的说法并无科学依据，只是一部分人制造的谣言。又因为此事与盐有关，并且“言”与“盐”谐音，故以“谣盐”代替“谣言”。

沈新晋等[2]以ABAQUS软件为平台，采用应变协调假设和强度等效假设，给出了一个描述钢筋混凝土材料与结构损伤塑性分析的本构模型，认为钢筋屈服后，立刻进入下降阶段，不考虑钢筋与混凝土相互作用引起的受拉刚化效应及受压箍筋约束效应。骈超等[2]借助ABAQUS软件中混凝土损伤本构模型，利用刚度等效原则将钢筋混凝土视为连续均质材料，对钢筋混凝土材料的本构关系进行分段线性简化，对钢筋与混凝土相互作用引起的受拉刚化效应进行了修正，但未考虑钢筋混凝土受压时的箍筋约束效应以及配筋率对受拉刚化效应的影响。宋伟等[3]在已有连续损伤理论的基础上，提出了综合考虑混凝土损伤、钢筋混凝土粘结滑移、以及受拉刚化效应等因素的钢筋混凝土拉伸本构模型。殴碧峰等[4]考虑钢筋和混凝土的相互作用，在Baant等提出的混凝土微平面模型的基础上通过引入钢筋的影响提出了一个钢筋混凝土动态本构模型。商霖等[5]提出了理想情况下钢筋混凝土本构关系可由混凝土材料黏弹性本构关系与一个依赖于增强材料特性的常量G的乘积确定的思想，得到了单向加筋、正交双向加筋和正交三向加筋混凝土的损伤型动态本构关系。孔丹丹等[6]在Zhao和Weng建立的纤维增强等效材料的等效模量理论的基础上，针对钢筋混凝土材料的力学特征，提出了钢筋单方向、三方向分布情况下钢筋混凝土材料的平均等效模量的简化公式。

大肠杆菌的基本特征：肠埃希氏菌通常称为大肠杆菌，是 Escherich在1886年发现的，在相当长的一段时间内，一直被当做正常肠道菌群的组成部分，认为是非致病菌。直到20世纪中叶，才认识到一些特殊血清型的大肠杆菌对人和动物有病原性，尤其对婴儿和幼畜（禽），常引起严重腹泻和败血症。

本文中，在已有研究的基础上，考虑钢筋混凝土受拉刚化效应、配筋率对受拉刚化效应的影响和钢筋混凝土受压箍筋约束效应，采用组合模量方法对钢筋混凝土复合材料进行模量等效，并对钢筋混凝土本构关系进行分段线性简化，借助LS-DYNA软件中混凝土损伤塑性模型，建立钢筋混凝土损伤塑性模型，在实验验证的基础上，应用此模型对一双切口钢筋混凝土烟囱同向延时爆破拆除效果进行数值模拟。

1 钢筋混凝土本构模型

1.1 钢筋混凝土拉伸力学行为

图1给出了一个含有钢筋和混凝土的二维方形钢筋混凝土样本材料单元[3]。图1中正方形边长为单位长度1，钢筋截面与混凝土截面的面积比为μ(配筋率)，0<μ<1，应力荷载如图所示。钢筋所受应力为σt,steel，混凝土所受应力为σt,concr。图1中的右图为简化的钢筋混凝土材料单元模型。由于平衡关系，CD的对边AB上所受的载荷相同(为了简单起见，没标出AB边上的载荷)。

图1 钢筋混凝土等效材料单元体示意图 Fig.1 Schematic diagram of reinforced concrete equivalent material unit

模型采用应变协调假设，即假设钢筋和混凝土间粘合完好，二者的应变量相等:

获得双边匹配方案：云服务需求企业X1和云服务供给企业Y1匹配，云服务需求企业X2和云服务供给企业Y4匹配，云服务需求企业X3和云服务供给企业Y5匹配，云服务需求企业X4和云服务供给企业Y2匹配，云服务需求企业X5和云服务供给企业Y6匹配，云服务供给企业Y3和Y7未匹配。

根据武汉市2016-12-29早高峰时段的公交乘客IC卡刷卡数据和GPS数据，对模型进行试算. 其中上车站点识别模型试算发现采用模型2(基于IC卡刷卡时间识别)识别，站点序号基本对应不上，识别失败. 而模型1(基于GPS时间识别)识别情况是：总的有效刷卡记录为23 206条，占所有刷卡记录85.67%(27 088条)，识别成功的记录数为22 814条，识别率为98.31%. 所以最终选取模型1为识别模型.

εconcr=εsteel=ε

(1)

由静力关系可得钢筋混凝土等效材料的组合应力σt为:

σt=[σt,steelμ+σt,concr(1-μ)]

(2)

钢筋混凝土等效均质材料单轴受拉应力在达到混凝土开裂应力前，认为钢筋与混凝土均处于弹性工作状态，当混凝土开始出现裂缝时，近似认为混凝土达到峰值拉应力，此时开裂应力为：

(3)

式中：Ee为钢筋混凝土等效材料单元组合模量，Ee=[Esteelμ+Econcr(1-μ)]；Esteel为钢筋弹性模量，Econcr为混凝土弹性模量。

在等效材料单元应力大于开裂应力但未达到钢筋屈服强度时，钢筋仍处于弹性阶段，假定钢筋与混凝土之间无相对滑移，通过调整钢筋刚度来模拟受拉刚化效应，荷载全部由钢筋来承担。考虑钢筋受拉刚化效应，混凝土峰值应力为：

针对喀斯特山区大白菜肥料需求不同及合理施肥，通过栽培2+X施肥进行试验，探索裸地大白菜氮肥的施肥方法及数量，氮、磷、钾肥施用量对大白菜产量的影响不同，缺磷、钾肥或钾肥施用过多会影响大白菜的产量，氮、磷肥过多不会提高产量，反而增加成本。本研究可为喀斯特山区修正和优化施肥技术提供参考，促进大白菜高产优质栽培种植提供科学施肥参考依据。

σt2=εt,limψ μEsteel

(4)

式中：εt,lim为混凝土极限拉应变，由钢筋混凝土等效材料配筋率决定，配筋率小于1%时，可取1×10-3，配筋率大于1%时，可取σt,steel/Esteel[3,7]，本文中为考虑受拉刚化效应的不均匀系数，取值为5[8]。

当等效材料单元应力超过初始屈服应力后，混凝土开始出现塑性，而钢筋依然处于弹性阶段。当混凝土达到约束混凝土峰值应力，认为纵向钢筋达到屈服应变，钢筋应力维持在受压屈服强度，此时钢筋混凝土等效材料达到最大屈服应力σc2：

表2中，拉格朗德Ⅱ级水电站的第1期只安排其可靠出力的2倍为它的装机容量值，第1期和第2期合计安排可靠出力的倍比为4。

σt3=μσt,steel

(5)

假设钢筋达到屈服强度后，混凝土已经全部拉裂，钢筋进入颈缩阶段，单元发生较大变形，钢筋混凝土等效材料受拉应力进入下降阶段。钢筋混凝土等效均质材料单轴受拉应力应变曲线见图2(a)。

(6)

对应的应变εt3为：

图2 等效材料应力-应变关系曲线示意图 Fig.2 Illustration of stress-strain relationship of the equivalent material

1.2 钢筋混凝土压缩力学行为

在2018年的调查中，61%的美国医生感受到了所在组织在采取针对倦怠的干预行动，52%的医生参与过相应训练，30%使用过心理门诊。尽管如此，美国医生在应对倦怠的11个行为中(人数占比90%)，有超过一半的行为都是不良行为或可能导致不良后果的行为，包括吸烟、暴食，甚至吸大麻。这说明对医生职业倦怠的干预应尽早采取行动，制定科学有效系统的干预计划，研究适合医生工作生活方式的干预手段，使医生倦怠问题控制在可控范围内，避免造成严重的医疗事故。

σc1=[Esteelμ+Econcr(1-μ)](0.3Fc,concr)/Econcr

(7)

式中：Fc,concr为素混凝土的抗压强度。

对应的应变εc1为：

εc1=σc1/E=0.3Fc,concr/Econcr

(8)

随着应变增加，钢筋应力逐渐增大，直到钢筋达到屈服应力。此时钢筋混凝土等效材料单元应力为：

(9)

式中：为约束混凝土峰值应力。

式中：为塑性应变率；为塑性乘数率；gp为塑性势。

(10)

式中：λ为箍筋综合约束指标，λ=μσt,steel/Fc,concr。

当钢筋混凝土等效材料轴压应变超过εc3时，由于箍筋的约束作用使得混凝土峰后应力应变曲线延性增加。根据《混凝土结构设计规范2010》可知，当混凝土强度等级在C30～C60时，混凝土的极限压应变εc3与峰值压应变εc2之比可近似认为是2，极限压应变对应的混凝土强度取0.15倍的峰值压应力，在此过程中认为钢筋一直保持屈服应力不变。钢筋混凝土等效材料在极限压应变时的应力σc3为：

前两个测点A、B相对于车辆悬浮架的坐标分别为(x1,z1)、(x2,z2)。通过A、B两点做直线，得到：

(11)

当钢筋混凝土等效材料应力小于极限压应力时，近似认为钢筋混凝土等效材料受压应力进入下降阶段。钢筋混凝土等效材料单轴受压应力应变曲线见图2(b)。

截至2017年12月，武隆区纳入低保保障对象共11504户21133人，其中城镇低保户3029户4716人，农村低保户8475户16417人；累计发放城乡低保金8697.6万元，其中城镇低保2509.6万元，农村低保6188万元。按照100元/人的标准，为47935个贫困人口购买了精准脱贫保，已赔付2633人160万元；健康扶贫医疗救助基金已赔付3324人150万元；解决大病医疗救助268.8万元，惠及贫困人口202人次。

2 钢筋混凝土损伤塑性模型

2.1 混凝土损伤塑性模型

混凝土损伤塑性模型是基于以下应力应变关系[10]：

(12)

式中：和分别为有效应力张量的正、负部分；ωt和ωc为两个标量损伤变量，取值范围0～1,0代表无损伤，1代表完全损伤破坏。

日日顺物流发展至今经历了四个阶段。自从2000年开始的企业物流，针对的是海尔的服务升级，从送货到送服务，到物流企业是卖服务和卖方案；第二阶段时开始进入快消品行业，同时也做了很多家电产业、原材料供应，积累了很多经验；第三阶段成为平台企业，从配送平台到交互平台，最核心的是跟用户的触点，最后一公里的服务。现阶段我们正在打造生态企业，完成从物流平台到生态平台的升级，进一步创造商业模式的价值。

有效应力可定义为：

(13)

式中：De为基于弹性模量E和泊松比ν的弹性刚度张量；ε为应变张量；εp为塑性应变张量。

式(13)中有效应力的正、负部分可通过有效主应力进行确定：

假定混凝土的弹性极限应力为0.3倍混凝土抗压强度，钢筋混凝土等效材料应力达到初始屈服应力σc1：

(14)

式中：<>+和<>-分别为正和负部分算子，<x>+=max(0,x),<x>-=min(0,x)，例如，对于一个结合拉伸和压缩应力状态与有效主应力分量正、负主应力分别为和

塑性模型是基于有效应力，与损伤无关。塑性模型[11]可以通过屈服函数，流动法则，硬化变量的演化法则以及加卸载条件进行描述。

屈服函数的形式可以表示为：

(15)

式中：qh1(κp)和qh2(κp)为控制屈服面大小和形状演化的无量纲函数，κp为硬化变量。

流动法则可以表述为：

(16)

约束混凝土峰值应力和峰值应变为[9]:

硬化变量κp的速率与塑性应变率的演化规律有关，加卸载条件为：

惨了！这两个字钻进步凡的脑海，他高速地运转自己的脑袋，想找出一个逃离的方法。可他还没来得及想出什么好办法，一个“龙死党”已经摁住了他的胳膊，并捂住了他的嘴巴。

(17)

损伤模型[12]可以通过损伤加载函数，加卸载条件以及拉伸和压缩损伤变量的演化规律进行描述。对于拉伸损伤方程可以表示为：

(18)

对于压缩损伤方程可以表示为：

(19)

式中：fdt和fdc为加载函数；和为等效应变；κdt,κdt1,κdt2,κdc,κdc1,κdc2为损伤历史变量；αc为区分拉伸加载和压缩加载条件的一个变量。

混凝土损伤塑性模型中，应变率效应通过引入因子αr(αr>1)， 值越大，损伤延迟的发生和强度提高越大。对于混凝土应变率效应，采用动力增大系数曲线体现动荷载作用下材料屈服强度的提高。混凝土抗压强度的动态增大系数采用欧洲规范CEB建议的公式[13]计算，混凝土抗拉强度增大系数采用修正的CEB公式[14]计算。因此，αr可以定义为：

【美国北极星医用放射性同位素公司网站2018年10月1日报道】 2018年10月1日，美国北极星医用放射性同位素公司(NorthStar Medical Radioisotopes)宣布获得美国能源部(DOE)国家核军工管理局(NNSA)颁发的杰出成就奖。该奖项主要表彰北极星即将成为美国近30年来的第一家钼-99生产商。

αr=(1-αc)αrt+αcαrc

(20)

式中：αrt和αrc为拉伸应变率效应因子和压缩应变率效应因子；αc为一标量(0<αc<1)，αc=1表示仅压缩，αc=0表示仅拉伸。

2.2 模型参数修改

LS-DYNA中混凝土损伤塑性模型参数及部分模型参数默认值[10]见表1。表1中：MID代表材料编号，Type代表损伤类型，ρ为密度，E为弹性模量，ν为泊松比，Ecc为偏心距，QH0为初始硬化参数，Ft为单轴抗拉强度，Fc为单轴抗压强度，Hp为硬化参数，Ah、Bh、Ch、Dh为硬化延性参数，As、Bs为软化延性参数，Df为流动法则参数，fco为应变率相关参数，Wf为线性损伤公式下拉伸阀值位移，Wf1、Ft1为双线性损伤公式下第二部分对应的拉伸位移阀值和拉伸强度阀值，Sf为应变率标识，Ff为失效方式标识。混凝土损伤塑性模型参数分别控制着应力应变曲线的弹性、硬化和软化三个阶段。

表1 LS-DYNA混凝土损伤塑性模型默认参数 Table 1 Default parameters of concrete damage plasticity model

MIDρEνEccQH0FtFc273 NoneNone 0.2AUTO0.3 NoneNone HpAhBhChDhAsDffco0.50.0080.0032.010-615.00.85AutoTypeBsWfWf1Ft1SfFf21None0.15Wf0.3Ft00

本文中，依据钢筋混凝土在单轴受拉、受压作用下的等效本构关系，基于LY-DYNA混凝土损伤塑性模型，对弹性模量E、抗拉强度Ft、抗拉强度阀值Ft1、断裂位移Wf、硬化参数Bh、抗压强度Fc、损伤延性参数As等8个参数进行修改，得到钢筋混凝土损伤塑性模型。具体的[11-12]：

E=Ec;Fc=σc2;Ft=σt2;Ft3=σt3;Wf=εt3hx/0.15,h为单元网格尺寸，x为钢筋混凝土与混凝土的等效抗拉强度之比；Bh=0.007 82ln [(εc2+1.242)/2.107]；Esoft=38.12As-1.11，Esoft为受压软化模量，由钢筋混凝土等效材料受压应力应变曲线中对应的特征点2和3应力应变数据计算得出。

3 钢筋混凝土损伤塑性模型验证

3.1 钢筋混凝土轴压短柱实验

图3 钢筋混凝土等效轴压柱应力应变曲线 Fig.3 Stress-strain curve of reinforced concrete equivalent axial compression column

图4 钢筋混凝土等效简支梁荷载位移曲线 Fig.4 Load displacement curve of reinforced concrete simple supported beam

钢筋混凝土短柱尺寸为150 mm×150 mm×450 mm，混凝土抗压强度fc=63.1MPa，配筋率为0.84%，纵筋屈服强度为548.1 MPa，箍筋采用HPB235，配箍率为0.82%，配箍特征值为0.03，具体见文献[15]。

根据以上实验参数，在数值模型中设置：钢筋混凝土等效材料的组合模量为38.52 GPa；εc2=2.6×10-3；Fc=σc2=68.7 MPa；εc3=2εc2=5.2×10-3；σc3=58.6 MPa；Ft=σt2=8.4 MPa；Ft1=σt3=4.49 MPa；Bh=4.7×10-3；As=7.7；Wf=25εt3σt2/(0.15σt,concr)=0.79 mm。

图3给出了钢筋混凝土轴压短柱采用钢筋混凝土损伤塑性模型和默认参数混凝土损伤塑性模型的数值模拟与实验结果，可以看出，钢筋混凝土损伤塑性模型数值模拟结果与实验结果较为吻合。

3.2 简支梁压弯实验

钢筋混凝土简支梁梁长1 400 mm，支座间距为1 200 mm，梁截面为矩形，宽150 mm，高200 mm。混凝土标号为C25，保护层厚度为25 mm。钢筋混凝土梁受拉和受压均配置2跟∅10钢筋(HRB335)，箍筋选用∅6@150(HPB235)，梁受拉与受压配筋率和体积配箍率分别为0.522 5%和0.628%。实验采用两点加载，加载速率为10 mm/s，单调加载，简支梁内混凝土的应变率将达到0.01 s-1以上，具体见文献[16]。

根据以上实验参数，在数值模型中设置：钢筋混凝土等效材料的组合模量为38.52 GPa，εc2=2.4×10-3；Fc=σc2=31.6 MPa；εc3=4.8×10-3；σc3=27.4 MPa；Ft=σt2=5.225 MPa；Bh=4.6×10-3；As=15.8；Ft1=σt3=1.75 MPa；Wf=25εt3σt2/(0.15σt,concr)=0.58 mm。

图4给出了钢筋混凝土简支梁采用钢筋混凝土损伤塑性模型和默认参数混凝土损伤塑性模型的数值模拟计算结果与实验结果，可以看出，钢筋混凝土损伤塑性模型数值模拟结果与实验结果基本一致。

4 应 用

图5 钢筋混凝土烟囱有限元模型 Fig.5 Finite element of reinforced concrete chimney

某钢筋混凝土烟囱高210 m，烟囱底部筒体外半径为9.45 m，壁厚600 mm，烟囱顶部外半径为3.35 m，壁厚250 mm，混凝土标号C30，竖向筋截面配筋率在0.42%～0.66%之间，环向筋配筋率平均为0.45%。烟囱上部爆破切口布置在90 m处，切口角度210°，切口高度为2.5 m，正梯形切口，定向窗角度为30°，底部爆破切口角度为216°，切口高度为3.5 m，正梯形切口，定向窗角度为30°，爆破方案设计上切口先于下切口3 s起爆。

钢筋混凝土烟囱有限元模型由3部分组成：筒体、钢筋和爆破切口，建模过程中混凝土采用Solid164单元，钢筋采用Beam161单元，由于烟囱竖向筋和环向筋配筋率不同，只能对环向配筋进行简化等效，竖向筋则按剩余配筋率正常建立，竖向筋与筒体之间采用分离式共节点的处理方式，钢筋与混凝土共计472 440个单元，钢筋混凝土烟囱有限元模型见图5。钢筋材料采用弹塑性硬化模型，密度为7 850 kg/m3，屈服强度为335 MPa，弹性模量为200 GPa，泊松比为0.3。筒体采用钢筋混凝土损伤塑性模型，钢筋混凝土损伤塑性模型参数见表2。地面采用刚体模型，密度为2 000 kg/m3，弹性模量为30 GPa，泊松比为0.2。烟囱重力采用关键字*LOAD_BODY_Y(Y方向)来实现。钢筋与刚体地面接触采用关键字*CONTACT_NODES_TO_SURFACE接触定义，防止钢筋触地穿透地面，筒体与刚体地面采用关键字*CONTACT_ERODING_SINGLE_SURFACE接触定义。爆破切口的形成采用关键字*MAT_ADD_EROSION中的时间控制参数定义，爆破切口范围内单元一次性控制删除，不考虑炸药爆炸以及毫秒延时对切口形成的影响。钢筋和混凝土材料失效均由应变控制，钢筋材料失效应变设为0.05，混凝土材料失效应变设为0.004。

表2 新钢筋混凝土损伤塑性模型参数 Table 2 Parameters of new concrete damage plasticity model

ρ/(kg·m-3)E/GPaνQH0Ft/MPaFc/MPaHpAhBhCh245029.50.20.35.0310.50.0080.0032.0DhAsDfTypeBsWf/mmWf1/mmFt1/MPaSfFf10-613.00.852150.751.510

采用钢筋混凝土损伤塑性模型，图6给出高耸双切口钢筋混凝土烟囱同向延时爆破拆除数值模拟倒塌过程与实际高速摄影观测结果，可以看出，实际高速摄影观测结果不论在倒塌运动状态还是在倒塌时间上均与数值模拟结果相符，说明采用钢筋混凝土损伤塑性模型模拟钢筋混凝土烟囱爆破拆除不仅简化了建模工作量，提高了计算效率，且计算结果与实际较为吻合，同时也证实了钢筋混凝土损伤塑性模型参数计算选取的正确性及合理性。

从“套路贷”的产生背景及概念的界定来看，“套路贷”脱胎于传统的民间借贷，是以民间借贷为诱饵，通过制造种种“合规”流程，设置重重诈骗陷阱而实施的违法犯罪行为。“套路贷”可以说是高利贷的高级发展形态：一方面，“套路贷”不同于传统的高利贷，借贷行为已由违法型向犯罪型转变，另一方面，借贷中的套路手段贯穿始终，借贷方式已由索利向侵财转变，不法分子分工明确，层层设套，令借款人防不胜防。总的来说，“套路贷”犯罪是指嫌疑人利用借款人急用钱款且自我保护意识不强等弱势情形，诱使借款人签订虚高金额的借款合同，随后通过恶意转贷垒高债务、肆意认定违约、非法索债、虚假诉讼等手段，最终达到非法占有借款人财物的犯罪行为。

图6 烟囱爆破拆除实际倒塌与数值模拟过程对比 Fig.6 Comparison between numerical simulation and actual collapse process of chimney blasting demolition

5 结 论

(1)借助LS-DYNA有限元软件中混凝土损伤塑性模型，采用组合模量的方法给出了一个钢筋混凝土损伤塑性模型，该模型综合了混凝土受拉和受压损伤、应变率效应，考虑了钢筋混凝土受拉刚化效应、受压箍筋约束效应以及配筋率对受拉刚化效应的影响，在保证钢筋混凝土受力性能的前提下，极大简化了建模工作量。(2)钢筋混凝土结构数值模拟采用钢筋混凝土材料模型受力更合理，在不影响计算精度的前提下，提高了计算效率。(3)拆除爆破工程具有唯一性和一定的高危性，数值模拟作为一种重要的分析手段发挥的作用愈加明显，本文采用钢筋混凝土材料模型模拟结果与实际高速摄影观测结果较为吻合，给出的钢筋混凝土损伤塑性模型可为钢筋混凝土结构爆破拆除工程数值模拟提供参考和借鉴。

参考文献:

[1] 沈新晋,王琛元,周琳.一个钢筋混凝土损伤塑性本构模型及工程应用[J].工程力学,2007,24(9):122-128.

SHEN Xinjin, WANG Chenyuan, ZHOU Lin. A damage plastic constitutive model for reinforced concrete and its engineering application[J]. Engineering Mechanics, 2007,24(9):122-128.

[2] 骈超,刘柱,陶云芳,等.修正的钢筋混凝土等效材料塑性损伤本构模型及其在抗震分析中的应用[J].建筑结构,2015,45(8):19-24.

PIAN Chao, LIU Zhu, TAO Yunfang, et al. Revised plasticity damage constitutive model of reinforced concrete equivalent material and its application in seismic analysis[J]. Building Structure, 2015,45(8):19-24.

[3] 宋伟,袁勇.单调轴向荷载下钢筋混凝土拉伸损伤本构研究[J].工程力学,2005,22(6):142-147.

SONG Wei, YUAN Yong. Damage constitutive law of reinforced concrete specimens under monotonic uniaxial tensile load[J]. Engineering Mechanics, 2015,45(8):142-147.

[4] 欧碧峰,王君杰.钢筋混凝土微平面动态本构模型[J].振动与冲击,2007,26(12):74-80.

OU Bifeng, WANG Junjie. Microplane dynamic model for reinforced concrete[J]. Journal of Vibration and Shock, 2007,26(12):74-80.

[5] 商霖,宁建国,孙远翔.强冲击载荷作用下钢筋混凝土本构关系的研究[J].固体力学学报,2005,26(2):175-181.

SHAN Lin, NING Jian guo, SUN Yuan xiang. The constitutive relationship of reinforced concrete subjected to shock loading[J]. Acta Mechanica Solida Sinica, 2005,26(2):175-181.

[6] 孔丹丹,赵颖华,王萍,等.钢筋混凝土材料有限元分析中的等效模量方法[J].沈阳建筑大学学报,2005,21(3):200-203.

KONG Dandan, ZHAO Yinghua, WANG Ping, et al. Equivalent moduli method in the finite element analysis of reinforced concrete structures[J]. Journal of Shenyang Architectural and Civil Engineering Institute, 2005,21(3):200-203.

[7] 龙渝川,周元德,张楚汉.钢筋混凝土相互作用效应的钢筋刚化模拟[J].清华大学学报,2007,47(6):793-796.

LONG Yuchuan, ZHOU Yuande, ZHANG Chuhan. Reinforcement stiffening model for reinforced concrete interaction effect[J]. Journal of Tsinghua University, 2007,47(6):793-796.

[8] 过镇海.钢筋混凝土原理[M].北京:清华大学出版社,1999,157-166.

[9] SAATCIOGLU M, RAZVI S R. Strength and ductility of confined concrete[J]. Journal of Structural Engineering, ASCE, 1992,118(6):1590-1607.

[10] Livermore Software Technology Corporation. LS-DYNA keyword user’s manual volume II: material models[M]. California: LSTC, 2016:1294-1302.

[11] GRASSL P, JIRASEK M. Damage-plastic model for concrete failure[J]. International Journal of Solids & Structures, 2006,43(22):7166-7196.

[12] GRASSL P, XENOS D, NYSTROM U, et al. CDPM2: A damage-plasticity approach to modelling the failure[J]. International Journal of Solids and Structures, 2013,50(24):3805-3816.

[13] Comite Euro-International du beton. CEB-FIP model code 1990[M].Trowbridge, Wiltshire, UK: Redwood Books, 1993:48-51.

[14] MALVAR L J, ROSS C A. Review of strain rate effects for concrete in tension[J]. ACI Material Journal, 1998,98(6):735-739.

[15] 刘立新,张艳丽,李千,等.500MPa钢筋混凝土短柱轴压下受力性能研究[J].郑州大学学报,2006,41(S2):1-5.

LIU Lixing, ZHANG Yanli, LI Qian, et al. The study performance of axle loaded reinforced concrete puncheons used 500MPa[J]. Journal of Zhengzhou University, 2006,41(S2):1-5.

[16] 肖诗云,曹闻博,潘浩浩.不同加载速率下钢筋混凝土梁力学性能试验研究[J].建筑结构学报,2012,33(12):142-146.

XIAO Shiyun, CAO Wenbo, PAN Haohao. Experimental study on mechanical behavior of reinforced concrete beams at different loading rates[J]. Journal of Building Structures, 2012,33(12):142-146.