Blast Resistance of Containment Dome Reinforced Concrete Slab in NPP under Close-in Explosion
-
摘要: 安全壳是核电厂的最后一道防线,其穹顶采用60°配筋混凝土进行设计和建造,配筋方式特殊。借助ANSYS/LS-DYNA,采用CONWEP爆炸模型,建立60°和普通配筋的混凝土板有限元模型,研究了近场爆炸作用下60°配筋混凝土板的动态响应,参数化分析了板厚、药量、钢筋屈服强度和混凝土强度等因素对60°配筋钢筋混凝土板抗爆性能的影响规律;对比研究了普通配筋和60°配筋混凝土板的中心挠度、变形和应力云图,基于数值分析结果,拟合得到两种配筋方式混凝土板中心挠度最大值与药量之间的关系曲线,利用回归分析得到其计算公式。研究结果表明:在相同含钢量的条件下,60°配筋混凝土板中心挠度最大提高60.22%,抗爆性能更强,拟合公式可以较好地预测60°配筋混凝土板的挠度变化。
-
关键词:
- 核电厂安全壳 /
- 60°配筋 /
- 钢筋混凝土板 /
- 近场爆炸 /
- CONWEP爆炸模型
Abstract: The containment dome of the nuclear power plant (NPP) is the last barrier and designed with reinforced concrete in a 60° of configuration. This paper uses explicit dynamic finite element analysis software ANSYS/LS-DYNA to establish finite element model of the ordinary reinforced concrete slab and reinforced concrete slab of 60° of configuration (novel) with CONWEP explosion model. Dynamic responses of reinforced concrete slabs under close-in explosion are investigated. Base on the analysis of parameters, the effects of slab thickness, explosive charge, strength of concrete and reinforcing steel bar on the blast resistance of reinforced concrete slabs with 60° of configuration are analyzed. The central deflection, deformation and stress diagrams of ordinary and novel reinforced concrete slabs are compared, a relationship of the maximum central deflection of novel reinforced concrete slabs, the explosive charge and slab thickness, and the fitting formula are obtained and given. The results indicate that under the same explosive charge, the central deflection of novel reinforced concrete is increased by 60.22%, and the novel reinforced concrete slabs have better blast resistance, the fitting curve can preferably estimate the deflection of the novel reinforced concrete slab. -
冲击波加载多孔固体介质过程中形成的极端高温高压条件可激活常规温压条件难以发生的化学反应,为材料合成开辟了一条新途径,尤其是对合成条件苛刻的材料具有显著优势。然而,冲击波作用导致粉末从初始状态至孔隙湮灭的过程中涉及多种机械、物理和化学变化,因而冲击波压缩下的反应过程十分复杂[1]。冲击诱导化学反应取决于外部加载条件和粉末混合物的初始状态等多种因素[2-3],反应物粉末的初始孔隙率是控制冲击反应热力学和动力学的一个非常重要因素[4]。Cooper 等[5]应用理论模拟证实,冲击波和孔隙表面相互作用导致孔隙崩塌产生的大量能量集聚有利于激活化学反应,因此探讨粉末初始孔隙率对诱发超快冲击化学反应的影响极为重要。
铌硅二元金属间化合物因具有低密度和优良的高温强度等物理特性,可作为潜在的高温结构材料和金属含能材料[6-7]。根据铌-硅二元相图,铌-硅系存在Nb5Si3、Nb3Si和 NbSi2共3种金属间化合物[8]。其中Nb5Si3的熔点(2 515 ℃)最高,密度 (7.16 g/cm3)相对较低,因其具有优良的高温性能受到关注[9]。这类材料的熔点高、晶体结构复杂,采用常规方法难以合成,但该类材料的合成反应能够放热,理论上讲,一旦被激活便可自维持发生反应。由于该激活势垒高,传统方法难以达到激活条件,而冲击波加载产生的极端高温高压对于激活该类反应具有一定优势。目前,已经开展了铌硅化合物的冲击合成工作。Vecchio等[10]应用冲击波加载合成了NbSi2,发现其冲击诱导反应机制与传统固态反应机制不同;Meyers等[11]也通过铌硅粉末混合物的冲击回收,发现塑性变形可以影响铌硅粉末间的反应。此外,还有一些理论模拟研究对铌硅金属间化合物的冲击合成进行探讨[12]。但是,所有铌硅粉末混合物的冲击反应研究,无论是实验结果还是理论分析都仅局限于较低冲击强度条件(飞片速度在2.0 km/s以下),冲击产物也仅获得了较低熔点的NbSi2,而且这些结果对于冲击化学反应的程度和机理解释存在分歧。而且,高熔点金属间化合物(如Nb5Si3)这类传统方法较难合成的材料在高冲击强度范围(飞片速度在2.0 km/s以上)的冲击合成及反应机制未有报道。
课题组前期开展了铌硅粉末混合物的冲击回收实验,通过高冲击强度下的冲击回收获得了高熔点Nb5Si3化合物[13],发现相同组分配比的铌硅粉末混合物在不同的冲击强度下合成了不同组成的铌硅化合物,且反应特征不同。为更清楚地理解高冲击强度下铌硅粉末混合物的冲击合成及反应机制,在前期工作基础上,本研究进一步对不同孔隙率的铌硅粉末混合物样品进行高冲击强度下回收,通过对冲击回收产物进行表征分析,探讨不同初始孔隙率对铌硅混合物冲击反应的影响。
1. 实验材料与方法
1.1 冲击样品制备
实验使用高纯Nb(约325目,Alfa Aesar公司)和Si(约325目,Alfa Aesar公司)粉末原料。将铌粉和硅粉按原子比为5∶3称量,V型混料器中氩气气氛下将粉末混合2 h,通过不同压力将混合好的粉末冷压成直径约为12 mm、厚度约为3 mm的圆片状试样。冲击样品的准备过程(原料称量、混合前的真空封装、压片和铜回收盒的真空封装)均在水氧监测氩气气氛保护的手套箱中完成。将装有样品的铜制回收盒密封。样品孔隙率的计算方法参考文献[14]。
1.2 冲击回收实验
应用二级轻气炮加载装置对样片进行冲击回收,装置示意图和实物如图1所示。将铜飞片(
∅ 24 mm × 3.0 mm)粘在发射弹丸上,由磁测速系统测量飞片速度,在飞片速度约2.35 km/s的最高冲击速度下进行冲击回收实验。利用质量平均法计算混合物样品的雨贡纽参数[15-16],采用阻抗匹配法和疏松粉末状态方程计算样品的冲击压力和温度[15,17],表1列出了不同孔隙率粉末混合物的冲击回收实验参数和计算结果。表 1 不同孔隙度的粉末混合物的冲击回收实验参数Table 1. Shock loading conditions of Nb-Si powder mixtures with different porositySample Flyer velocity/
(km·s−1)Density/
(g·cm−3)Porosity/
%Shock
pressure/GPaSecond shock
pressure/GPaShock
temperature/KNb-Si-P1 5.505 10 45 60 1 173 Nb-Si-P2 2.35 ± 0.02 4.895 20 39 60 1 625 Nb-Si-P3 3.980 35 29 60 2 256 1.3 冲击回收产物表征
图2显示了完整回收的铜回收盒在冲击前、后的外形变化。用机械切割方式从铜制回收盒中将回收的冲击产物分离取出,切掉产物样品外围边缘的含铜部分,留下中心主体部分用于表征分析。将产物样块逐级打磨,并清洗、干燥,储存在真空干燥箱中待表征。
使用X 射线衍射仪(丹东浩元仪器有限公司,DX-2700型)对样品进行物相表征,采用CuKa 靶(λ = 0.154 06 nm) 射线,衍射角2θ为20°~80°,步长为0.03°,每步计数时间1 s。使用日本电子公司JSM-7001F型场发射扫描电子显微镜 (SEM)进行微观形貌观察,并用SEM配备的能谱仪分析各物相的成分。使用差示扫描量热计(DSC)测量冲击回收产物的热响应曲线(Ar气氛下,加热速率为10 ℃/min)。
2. 结果和讨论
2.1 冲击回收产物的XRD结果分析
图3为不同孔隙率铌硅粉末冲击产物的XRD图谱。从图3(a)中10%孔隙率冲击样品的XRD结果可见,产物主要为Nb和Si的特征衍射峰及较弱的NbSi2衍射峰,表明仍有大部分Nb和Si未发生反应,仅有极少部分NbSi2相生成。由图3(b)中20%孔隙率铌硅粉末冲击产物的XRD结果可见,该孔隙率产物中虽仍有部分未反应的Nb和Si,但生成NbSi2相的程度有所提高,初始孔隙率提高后,冲击化学反应更易进行。此外,图3(c)中35%孔隙率铌硅粉末混合物的冲击产物XRD图谱显示,其衍射峰仅由Nb5Si3(包括α-Nb5Si3相和β-Nb5Si3相)的特征峰组成。与低孔隙率样品相比,在相同冲击强度下高孔隙率粉末样品更容易获得高熔点Nb5Si3金属间化合物。因此,初始孔隙率对粉末混合物在冲击压缩下物相生成、固相化学反应的影响非常明显。参考表1估算的冲击平均温度也可看出,随着样品初始孔隙率增加,较高的冲击温度(2 256 K)更有利于铌硅金属间化合物发生反应。
2.2 冲击回收产物的SEM结果分析
图4为10%初始孔隙率铌硅粉末冲击产物的SEM结果,其中图4(b)为图4(a)的局部放大图,图4(c)为样品EDS能谱分析结果,w为质量分数。由图4(a)可以看出,样品在该冲击压力(约为45 GPa)下铌硅颗粒被明显压实,没有明显物相生成或发生化学反应。图4(b)中局部区域放大图显示,与初态未冲击样品相比,冲击后的粉末颗粒呈细粒化、等轴状,说明冲击加载作用导致颗粒发生崩裂和相互嵌合。EDS能谱显示图4(b)中亮灰色细粒子的为Nb,周围相对较暗细粒为Si,可能由于NbSi2量极少,在图中未能明确分辨。在较小孔隙率即相对密实的混合反应物样品中,冲击波更多作用于粉末颗粒的机械和物理变形。由于孔隙率低,孔隙崩塌所致的能量贮集弱,极少能够激活并促发冲击化学反应。
图5为20%初始孔隙率铌硅粉末冲击回收产物的SEM结果,其中图5(b)为图5(a)的局部放大图,图5(c)为样品EDS能谱分析结果。图5(a)中,经冲击加载后20%孔隙率回收产物与10%孔隙率回收产物明显不同,随着冲击作用下温度不断升高,由于Si的熔点较低,而Nb的屈服强度和熔点较高,Si会发生熔融,并分布于变形的Nb周围,与基体特征类似,并且其中弥散着新相。如图5(b)所示,对富Si区域高倍放大观察可见,Si基体中分布着球形、结节状粒子,表现出发生反应的微观特征。对图5(b)方框内粒子进行EDS元素分析发现,其化学组成与NbSi2化学式的原子个数比相符,与XRD结果中出现NbSi2衍射峰一致。由此可见,与10%孔隙率样品相比,在相同飞片速度冲击下,随着冲击温度升高,高初始孔隙率反应产物中的铌硅颗粒发生严重变形并伴随硅发生熔融,并且NbSi2产物增多,表明发生了明显的冲击化学反应。
图6为高孔隙率(35%)铌硅粉末混合物冲击回收样品的SEM结果。由图6(a)可知,样品已转化为相对均匀的单一微观组织结构,表明该孔隙率样品在相同冲击速度下已发生完全反应形成新的相结构。图6(b)为图6(a)中方框内区域局部放大图,图中的新相为枝条状结构,其间有微纳孔(孔隙率约为4.7%,孔隙尺寸介于0.2~4 μm之间)。EDS能谱分析结果显示,铌和硅的原子比约为63.02:36.98,与Nb5Si3化学式的原子个数比接近,表明冲击反应生成物为Nb5Si3化合物。35%孔隙率样品中,冲击加载创造了极高的温度条件(估算冲击样品中的平均温度达2 256 K),相同冲击速度作用于铌硅粉末在物相生成和反应程度上与前述低孔隙率样品相比显著提高,并且其物相形态有明显不同。
值得注意的是,高孔隙率样品在飞片速度约为2.35 km/s时就发生了完全反应,而前期对不同冲击速度下反应行为的研究发现,较低孔隙率(20%)样品[13]需要更高的冲击速度(飞片速度约为2.63 km/s)才能发生完全反应。因此,高孔隙率样品崩塌所致的高能量贮集产生的高温是发生完全反应的重要条件。
2.3 冲击回收产物的DSC分析
图7为不同孔隙率铌硅粉末冲击产物的DSC曲线。由图7可知,孔隙率为10%和20%冲击产物的DSC曲线分别在720 ℃和810 ℃附近出现了两个放热峰。结合冲击产物的XRD结果分析可知,720 ℃放热峰可能是冲击产物中未反应的Nb和Si之间进一步发生反应放热所致。可能因为未发生完全化学反应的冲击产物中,未反应的铌和硅粒子在比自蔓延燃烧反应开始温度低约200 ℃时就开始发生化学反应[18]。该异常是由于冲击作用引起粒子粉末发生塑性流动和混合,表面清洁以及样品自身缺陷导致在冲击压缩加载期间粒子粉末被活化。而810 ℃放热峰可能是由于冲击塑变的Nb颗粒发生再结晶。35%孔隙度冲击回收产物的DSC结果呈现光滑的曲线,在测试加热范围内未出现放热峰,表明冲击回收产物为稳定的化合物。这与前述的物相表征和形貌分析结果一致,也证实了35%孔隙度样品完全反应生成了Nb5Si3金属间化合物。
粉末颗粒混合物的初始状态涵盖颗粒自身组成结构及混合颗粒间的相互接触等条件[19-22]。冲击加载时,在冲击波的高压作用下,样品颗粒的变化不是孤立发生的,而受相互间直接作用和颗粒自身物理性质的限制,表现为样品内颗粒的整体协调变形及交互作用。这些作用对于颗粒的变形方式、变形程度和融合都会产生直接影响。更重要的是,不同变化方式和程度会导致不同的能量集聚,从而产生温度升高和局域热效应[23],这些变化又反过来影响物质的扩散、相互渗透甚至激活化学反应。研究表明[24-25],合适的孔隙率对于粉末颗粒反应物的冲击化学合成非常重要。
由表征分析结果可知,孔隙率对于铌硅粉末混合物冲击反应的影响,是导致其反应发生及反应机制变化的重要因素。20%孔隙率冲击反应伴随着铌硅颗粒间强烈的相互作用、较低的相互约束以及大面积的孔隙塌陷或湮灭,冲击加热效果显著提升,不仅导致硅粒子发生熔化,变形的铌粒子嵌入到熔融的硅中,而且以铌粒子所在位置为形核中心,铌与周围的硅粒子发生半“固-液”化学反应生成相应的铌硅化合物,不过该激活程度仍处于仅能激发NbSi2反应物生成的条件。在更高孔隙率(35%)下,铌硅颗粒交互作用、变形及孔隙塌陷引起的能量集聚很高,使得样品内部急剧升温,引发高熔点的铌熔化,熔融的铌硅原子发生自由迁移,在极高温度激活下,铌硅间发生“液相”化学反应获得相应的铌硅化合物,即全面激活了铌硅间的自维持反应,并获得高熔点的Nb5Si3金属间化合物且反应完全。值得注意的是,过大孔隙率会在样品中留下微纳尺度孔隙,可能会对样品的块体性能产生影响。此外,实验研究发现,若孔隙率过高,孔隙崩塌产生的过高温度可能使冲击回收样品盒局部熔化,出现破漏或样品内气体相较多导致回收盒爆裂,从而导致冲击回收样品失败。
3. 结 论
通过3种不同初始孔隙率(10%、20%和35%)铌硅粉末混合物的冲击化学反应实验,对冲击回收产物物相、微观形貌表征分析及冲击产物热响应行为进行研究分析,发现不同孔隙率铌硅粉末混合物冲击反应的发生情况、产物特征及反应程度有所不同。低孔隙率(10%)样品几乎不发生反应;当孔隙率提高至20%时,铌硅粉末混合物发生不完全化学反应,生成NbSi2化合物;当孔隙率增大到35%时,相同冲击速度下铌硅粉末混合物发生完全反应,并获得与初始配比一致的Nb5Si3化合物。实验结果证实,适当增大铌硅粉末的初始孔隙率有利于发生冲击化学反应。
-
表 1 试件几何尺寸和参数
Table 1. Geometry and parameters of specimens
No. Dimension/(cm×cm×cm) fcu.k/MPa Arrangement of reinforcement Reinforcement bar W/kg R/m Z/(m·kg–1/3) O1-1 120×120×4 30 ∅8@100 HPB300 1.36 0.4 0.361 O1-2 120×120×5 30 ∅8@100 HPB300 1.36 0.4 0.361 O1-3 120×120×6 30 ∅8@100 HPB300 1.36 0.4 0.361 O2-1 120×120×4 30 ∅8@100 HPB300 0.41 0.4 0.538 O2-2 120×120×4 30 ∅8@100 HPB300 6.52 0.4 0.214 O3-1 120×120×4 50 ∅8@100 HPB300 1.36 0.4 0.361 O3-2 120×120×4 60 ∅8@100 HPB300 1.36 0.4 0.361 O4-1 120×120×4 30 ∅8@100 HRB335 1.36 0.4 0.361 O4-2 120×120×4 30 ∅8@100 HRB400 1.36 0.4 0.361 N1-1 120×120×4 30 ∅8-60° HPB300 1.36 0.4 0.361 N1-2 120×120×5 30 ∅8-60° HPB300 1.36 0.4 0.361 N1-3 120×120×6 30 ∅8-60° HPB300 1.36 0.4 0.361 N2-1 120×120×4 30 ∅8-60° HPB300 0.41 0.4 0.538 N2-2 120×120×4 30 ∅8-60° HPB300 6.52 0.4 0.214 N3-1 120×120×4 50 ∅8-60° HPB300 1.36 0.4 0.361 N3-2 120×120×4 60 ∅8-60° HPB300 1.36 0.4 0.361 N4-1 120×120×4 30 ∅8-60° HRB335 1.36 0.4 0.361 N4-2 120×120×4 30 ∅8-60° HRB400 1.36 0.4 0.361 表 2 钢筋参数
Table 2. Parameters of reinforcement bar
ρ/(g·cm−3) E/GPa ν σ0/MPa Etan/MPa β C/s−1 7.89 200 0.3 335 1000 0 40 表 3 数值计算与试验结果对比
Table 3. Comparison of calculated and test results
R/m W/kg Z/(m·kg−1/3) Maximum deflection/mm Error/% Finite element simulation Test[16] 0.4 0.46 0.518 37.7 35.2 7.10 -
[1] YAO S, ZHANG D, CHEN X, et al. Experimental and numerical study on the dynamic response of RC slabs under blast loading [J]. Engineering Failure Analysis, 2016, 66: 120–129. doi: 10.1016/j.engfailanal.2016.04.027 [2] 汪维, 张舵, 卢芳云, 等. 钢筋混凝土楼板在爆炸荷载作用下破坏模式和抗爆性能分析 [J]. 兵工学报, 2010(Suppl 1): 102–106WANG W, ZHANG D, LU F Y, et al. Analysis for blast resistance and damage mode of reinforced concrete slab subjected to explosbe load [J]. Acta Armamentarii, 2010(Suppl 1): 102–106 [3] 王靖. 近距离爆炸下钢筋混凝土板局部破坏试验研究与数值模拟 [D]. 天津: 天津大学, 2017.WANG J. Experimental and numerical studies on thelocal damage of reinforced concrete slabsunder close-in explosions [D]. Tianjin: Tianjin University, 2017. [4] 王淞, 刘筱玲. 双向支撑钢筋混凝土板在爆炸荷载作用下破坏模式与影响因素分析 [J]. 四川建筑, 2017, 37(2): 160–162 doi: 10.3969/j.issn.1007-8983.2017.02.054WANG S, LIU X L. Analysis of influence factors and on the dynamic response of reinforced concrete slabs of double supported under blast loading [J]. Sichuan Architecture, 2017, 37(2): 160–162 doi: 10.3969/j.issn.1007-8983.2017.02.054 [5] PANDEY A K, KUMAR R, PAUL D K, et al. Non-linear response of reinforced concrete containment structure under blast loading [J]. Nuclear Engineering & Design, 2006, 236(9): 993–1002. [6] ZHAO C F, CHEN J Y. Damage mechanism and mode of square reinforced concrete slab subjected to blast loading [J]. Theoretical & Applied Fracture Mechanics, 2013(63/64): 54–62. [7] 赵春风, 陈健云. 内爆荷载作用下钢筋混凝土安全壳的非线性响应分析 [J]. 爆炸与冲击, 2013, 33(6): 667–672 doi: 10.3969/j.issn.1001-1455.2013.06.017ZHAO C F, CHEN J Y. Dynamic responses of reinforced concrete containment subjected to internal blast loading [J]. Explosion and Shock Waves, 2013, 33(6): 667–672 doi: 10.3969/j.issn.1001-1455.2013.06.017 [8] 赵春风. 强震及爆炸荷载作用下核岛厂房动力响应及减震抗爆措施研究 [D]. 大连: 大连理工大学, 2014: 1–20.ZHAO C F. Study on dynamic response and measures of seismic reduiction and anti-exposion of nuclear island subjected to strong earthquakes and blast loads [D]. Dalian: Dalian University of Technology, 2014: 1–20. [9] ZHAO C F, CHEN J Y, WANG Y, et al. Damage mechanism and response of reinforced concrete containment structure under internal blast loading [J]. Theoretical & Applied Fracture Mechanics, 2012, 61(1): 12–20. [10] 刘云飞, 王天运, 贺锋, 等. 核反应堆预应力钢筋混凝土安全壳内爆炸数值分析 [J]. 工程力学, 2007, 24(8): 168–172 doi: 10.3969/j.issn.1000-4750.2007.08.030LIU Y F, WANG T Y, HE F, et al. Numerical simulation for pre-stressed concrete containment under internal explosive loading [J]. Engineering Mechanics, 2007, 24(8): 168–172 doi: 10.3969/j.issn.1000-4750.2007.08.030 [11] ZHANG T, WU H, FANG Q, et al. Numerical simulations of nuclear power plant containment subjected to aircraft impact [J]. Nuclear Engineering & Design, 2017, 320: 207–221. [12] MALVAR L J, ROSS C A. A Review of strain rate effects for concrete in tension [J]. Construction and Building Materials, 1998, 95(6): 735–739. [13] CLARK L A. CEB-FIP model code 1990 [J]. Programs Usenix Unix Supplementary Documents, 2008, 40(95): 233–235. [14] HALLQUIST J O. LS-DYNA keyword user’s manual [M]. Livermore: Livermore Software Technology Corporation, 2007: 1–2. [15] MALVAR L J. Review of static and dynamic properties of steel reinforcing bars [J]. Materials Journal, 1998, 95(5): 609–616. [16] 汪维. 钢筋混凝土构件在爆炸载荷作用下的毁伤效应及评估方法研究 [D]. 长沙: 国防科学技术大学, 2012.WANG W. Study on damage effects and assessments method of reinforced concrete structural members under blast loading [D]. Changsha: National University of Defense Technology, 2012. -