文彦博，黄瑞源，李 平，马 剑，肖凯涛

（1. 南京理工大学瞬态物理国家重点实验室，江苏 南京 210094；

2. 安徽工业大学机械学院，安徽 马鞍山 243032；3. 江苏科技大学船舶与土木工程学院，江苏 镇江 215600；4. 防化研究院，北京 102205）

随着国民经济的不断发展，高层建筑、地铁以及其他一些混凝土材料建造的工程结构越来越多。除了设计范围内的正常载荷外，某些工程结构在使用期限内可能会因偶然事故受到冲击或发生爆炸，譬如：燃气爆炸、恐怖袭击中人为引起的爆炸[1-2]，导致工程结构中的混凝土材料同时处于高温、高应变率的环境中。因此，有必要对混凝土材料在高温、高应变率下的力学性能和损伤演化关系进行研究[3-5]。混凝土是由水泥、石子、细沙、水和减水剂等组成的一种多相复合材料，其成分复杂，且内部存在多种尺寸不等、形状各异的微裂纹和微孔洞[6]，在外部载荷作用下，微裂纹和微孔洞将演化发展，又会产生新的微裂纹和微孔洞，从而影响整个混凝土材料的力学性能[7]。目前，混凝土类材料的损伤研究一直是国内外学者关注的热点问题。Dougill[8]和Loland[9]在混凝土本构关系中引入了损伤的概念；Xie[10]研究了岩石材料中的裂纹扩展；Holmquist 等[11]将混凝土中的损伤定义为等效塑性应变和等效塑性体积应变的表达式；王春来等[12]根据Weibull 统计分布理论和等效应变假定原理给出了一种损伤本构模型；Bai 等[13]提出了一种微孔洞成核模型；王道荣等[7]采用新技术对混凝土材料实施损伤“冻结”实验，研究损伤演化与应变和应变率的关系。

虽然混凝土材料的损伤已有较多研究成果，但多集中于静态拉压下损伤演化的研究，或者根据某些假定给出其在特定工况下的损伤演化规律。而混凝土材料中损伤的演化发展与不同应力状态相关，作为一种常用的军事防护材料[14-15]，其在冲击、碰撞及爆炸过程中经历了大变形、高应变率、高压及高温作用[11,16-18]。因此，研究混凝土材料在高温、高应变率下的损伤演化方程对于防护工事的设计具有重要的指导意义。

本研究将对不同温度(20、200、400 ℃)下的C45 混凝土材料进行分离式霍普金森压杆（SHPB）实验，结合实验数据给出损伤变量关于塑性应变的关系式，通过不同温度和不同应变率下的实验数据确定损伤演化方程的相关参数，从而给出一种适用于高温、高应变率混凝土材料的损伤演化方程。

1 试件制备

制备试验的原材料：水泥、细骨料（砂）、粗骨料、粉煤灰、水和外加剂。水泥为张家港海螺水泥有限公司生产的P·O42.5 普通硅酸盐水泥；选用石英砂作为细骨料，颗粒级配采用筛分试验测定，细度模数为2.39，属于中砂；粗骨料选用张家港市的瓜子片，最大粒径8 mm；外加剂使用安徽省法尔胜科技有限公司生产的聚羧酸高性能减水剂。混凝土原材料混合后达到如表1 所示的配比要求，且水灰比为0.346，28 天抗压强度达到46.2 MPa。根据李胜林等[19]的研究结果，对于 ∅74 mm 大尺寸SHPB 试验装置，为满足“一维假定”和“应力均匀假定”，混凝土试件的长度应为30～74 mm，模具浇注形成试验试件，尺寸为 ∅70 mm× 35 mm。试件养护过程符合工程要求，并利用磨床对两端面进行研磨，试件尺寸误差小 于 ± 0.02 mm，不平行度小于0.02 mm。

表1 混凝土材料的配比Table 1 Mixture ratio of the concrete material kg/m3

2 实验装置

如图1 所示，SHPB 是测量材料动态力学性能的重要仪器，由发射系统、波形输入/输出系统、波形存储系统和数据处理系统组成， εi为入射应变信号， εr为反射应变信号， εt为透射应变信号。本研究是在 ∅74 mm 口径的SHPB 装置上完成的，其撞击杆、入射杆和透射杆的长度分别为500、5 461 和3 485 mm，装置实物如图2 所示。为了得到混凝土材料在高温下的应力-应变曲线，采用管式高温炉装置加热混凝土试件[20]，高温炉和温度控制器如图3 所示。将试件加热到预定温度值后，保持恒温30 min 以确保试件的内部温度均匀，将发射装置的空气压力调整至预定值，然后使用特定夹具将加热的试样快速安装到指定位置，进行冲击测试。这样做的目的是为了避免加热炉加热试验系统，影响试验结果。整个过程在多人协作下完成，从组装开始到试验完成仅需7～10 s，大幅减少了混凝土表面的热损失。文献[21]显示，10 s 内混凝土的热损失在5 ℃左右，为了抵消实验过程中的热量损失，设置预定温度值比实验温度高5 ℃。由于混凝土是一种热惰性材料，传热系数非常低，最终试件的温度偏差控制在3 ℃以内。通过SHPB 实验装置及高温加热装置对C45 混凝土开展了3 种温度（20、200、400 ℃）下的动态压缩实验 ，在同一温度下对试件进行3 次实验，并通过改变子弹初速度调整应变率，相关测试信息见表2。

图1 SHPB 装置示意图Fig. 1 Schematic diagram of the SHPB device

图2 SHPB 装置Fig. 2 Image of the SHPB device

图3 高温加热装置Fig. 3 High temperature heating device

表2 测试相关信息Table 2 Summary of the test information

3 混凝土材料的应力-应变曲线

采用黄铜材料作为整形器，利用三波法[22]得到混凝土材料的应力-应变曲线，通过文献[23]中的方法确定应变率。图4 给出了C45 混凝土材料在不同温度(20、200、400 ℃)、不同应变率下的应力-应变曲线。可以看出，随着应变率增加，混凝土材料在不同温度下的峰值应力和峰值应变都随之增大。

图4 不同温度、应变率下混凝土材料的应力-应变曲线Fig. 4 Stress-strain curves of concrete materials with different temperatures and strain rates

为了进一步研究温度和应变率对峰值应力的影响，图5 给出了不同温度(20、200 、400 ℃)下峰值应力与应变率的关系。可见，随着温度不断增加，峰值应力随之增大，混凝土材料发生温度硬化和应变率硬化现象。受混凝土制备工艺的影响，试件本身含有大量的初始微裂缝和微孔洞。当应变率较低时，损伤主要沿初始微裂缝扩展、增大；随着应变率提高，损伤来不及扩展，从而衍生出大量的微裂纹，而微裂纹的产生需要大量能量，因此只能以提高强度或增加应变来抵消外部冲量，此时就表现出明显的应变率效应。不同温度下混凝土的实际破坏结果见图6。由于混凝土材料具有惯性效应和损伤滞后效应，表现为混凝土的强度随着应变率的增大而增大，试件的最终破坏结果表现为破碎程度随应变率的增加逐渐增大。由于混凝土试件内部结构随着高温的作用变得越来越松散，加载下其表面和内部的微裂纹更易演化发展，因此温度效应表现为在高温条件下混凝土试件的粉碎破坏程度更严重。

图5 不同温度下峰值应力和应变率的关系Fig. 5 The relationship between peak stress and strain rate at different temperatures

图6 温度不同时混凝土试件的破碎情况Fig. 6 Smashed concrete specimens at different temperatures

4 损伤与塑性应变的关系

从微观角度看，混凝土中存在着大小不一、形状各异的微裂纹，随着混凝土材料进入塑性阶段，内部的微裂纹开始演化发展并产生新的裂纹，混凝土材料的内部损伤随着塑性变形的程度增加不断变大。因此，混凝土材料内部的损伤演化与塑性应变的发展是相关的[11]。

4.1 损伤演化方程

Lemaitre[24]提出受损材料的本构关系可以从无损材料的本构关系推导得到，李永池等[25]、王春来等[12]在此基础上提出了相关的损伤本构模型。在单轴压缩情况下，混凝土材料损伤后的本构关系[12]可表示为

式 中： σ为 应力， ε为应变，E 为弹性模量，D 为损伤变量。对式（1）变换得到损伤变量

应 变 ε可表示为

式 中： εe为弹性应变， εp为 塑性应变。弹性应变 εe可 由应力 σ得到

塑性应变

图7 给出了不同温度、应变率下损伤变量D 与塑性应变 εp的关系，其中损伤变量D 是先根据应力-应变曲线的线弹性段求得弹性模量E，再将应力、应变代入式(2)计算得到的。从图7 中可以看出，随着应变率增加，在相同塑性应变下，其对应的损伤变量相对较小，混凝土材料内部微裂纹的扩展随着应变率的增加受到一定的抑制作用。因此，可以认为应变率对混凝土力学行为的影响主要体现在其对损伤变量D 的演化发展上，所以损伤变量D 可以表示为含应变率 ε˙ 和 塑性应变 εp的函数，由于损伤变量D 与塑性应变 εp之 间大致为指数函数关系，并且当塑性应变 εp为零时材料的损伤变量D 也为零，在此基础 上可以给出损伤演化方程

图7 实验和拟合得到的损伤变量-塑性应变曲线对比Fig. 7 Comparison of experimental and fitted results of damage variable-plastic strain curves

采用式(6)拟合不同温度、应变率下的 D- εp曲线，拟合曲线与实验曲线的对比如图7 所示。可以看出，拟合曲线与实验曲线的重合度非常高，得到的不同温度、应变率下材料参数A 和B 的值如表3 所示。假定混凝土材料进入塑性时对应的应力为峰值应力的70%[12]。从实验结果还可以得出，在温度和应变率不同时，C45 混凝土破坏时刻对应的损伤变量介于0.7～0.8 之间。

由于在不同温度下得到的材料参数A 和B 与应变率有关，为了确定不同温度下材料参数A 和B 与应 变率的关系，采用式(7)对其进行拟合，图8 给出了材料参数A 和B 与应变率关系的拟合曲线。

表3 不同温度、应变率条件下的材料参数A 和BTable 3 Material parameters A and B at different temperatures and strain rates

图8 材料参数A、B 与应变率的关系Fig. 8 The relationship between material parameters A，B and strain rate

式中：a、b、c、d、e、f 均为拟合参数。因此，不同温度、应变率下混凝土材料损伤演化方程的具体形式为

4.2 本构关系验证

图9 弹性模量与应变率的关系Fig. 9 Relationship between elastic modulus and strain rate

混凝土材料的损伤本构关系可表示为

弹性模量随着应变率的增加而增加[26]，图9给出了不同温度下弹性模量与应变率的关系，通过拟合发现弹性模量与应变率可描述为

将式(8)中不同温度、应变率下混凝土材料损伤演化方程的具体形式代入式(9)混凝土材料的损伤本构关系中，图10 比较了不同温度、应变率下实验和计算的应力-应变曲线，可以看出，采用本研究给出的损伤演化方程可以有效描述混凝土材料在高温、高应变率下的损伤演化过程，对工程分析和设计具有一定的参考价值。

图10 实验和计算的应力-应变曲线比较Fig. 10 Comparison of experimental and calculated results of stress-strain curves

5 结 论

通过采用 ∅74 mm 大口径分离式霍普金森压杆装置对C45 混凝土材料在不同温度(20、200、400 ℃)下进行动态力学性能实验，得到以下结论。

(1)在20～400 ℃之间，随着温度和应变率的增加，C45 混凝土材料的峰值应力和峰值应变呈增加趋势，并且在高温条件下混凝土试件的破坏程度比常温下更高。这是由于试件本身含有大量的初始微裂缝和微孔洞，当应变率较低时，损伤主要沿初始微裂缝扩展、增大，而随着应变率提高，损伤来不及扩展，从而衍生出大量的微裂纹，而微裂纹的产生需要大量能量，因此只能以提高强度或增加应变来抵消外部冲量。

(2)根据理论推导给出了损伤变量D 与塑性应变 εp的关系，结合实验得到的损伤变量D 与塑性应变 εp曲线确定了混凝土材料的损伤演化方程，并通过SHPB 实验数据确定了不同温度及应变率下损伤演化方程参数。实验结果还表明在不同温度及不同应变率下，假定混凝土材料进入塑性时对应的应力为峰值应力的70%，那么C45 混凝土破坏时刻对应的损伤变量介于0.7～0.8 之间。

(3)将研究中的损伤演化方程应用于混凝土材料的本构关系中，对本实验中不同温度、应变率下应力-应变曲线进行预测，预测结果与实验数据具有较好的一致性，进一步证明了该损伤演化方程形式可以有效描述混凝土材料在高温、高应变率下的损伤演化过程。