A method of geometry optimization for dynamic tensile specimen based on artificial neural network and genetic algorithm
-
摘要: 材料动态拉伸力学性能测试中,动态拉伸试样的几何尺寸对测试结果的准确性与有效性有着较大影响。为对动态拉伸试样的结构进行优化设计,以使其在动态拉伸过程中更好地满足一维应力与变形均匀等基本假设。首先,建立了量化的试样测量准确度指标,即应力平衡达到时间、变形均匀度、非轴向应力相对水平、过渡段相对变形。然后,对试样结构参数进行正交试验设计,通过数值模拟的方法得到了关于试样尺寸与测量准确度指标的正交试验数据库,并对正交试验数据库进行多目标正交试验矩阵分析,得到了试样结构参数对各测量准确度指标影响的主次顺序和规律。最后,以正交试验数据库为训练集,采用人工神经网络(artificial neural network, ANN)协同遗传算法(genetic algorithm, GA)的全局寻优方法对试样的结构尺寸进行优化设计,得到了试样的最优结构尺寸,并对最优尺寸的有效性进行了验证。结果表明,优化后的试样结构在材料动态拉伸力学性能测试精度上的表现明显得以提升。因此,采用ANN-GA协同优化的方法对动态拉伸试样的结构进行优化具有可行性和有效性。Abstract: The split Hopkinson tensile bar is one of the most commonly used apparatuses to test the dynamic tensile mechanical properties of materials at the high strain rates from 102 s−1 to 103 s−1, in which the specimens with a dog-bone shape are usually used. The dimensions of the specimen used are critical to ensure the basic assumptions during dynamic tensile process, such as one-dimensional stress state and uniform deformation of the specimen etc. And whether these assumptions can be satisfied would affect the measurement accuracy of the dynamic tensile properties directly. So, it is urgent to study the influence of the specimen structural parameters on the stress and deformation states of the specimen during the dynamic tensile tests. At the same time, developing and establishing an effective method which can realize the global optimization of specimen structural parameters in the entire parameter space is crucial. In order to actualize the above research objectives, indicators which can quantify the measurement accuracy of the dynamic tensile tests were firstly proposed, namely the time required to reach the stress equilibrium, the deformation uniformity, the relative level of the non-axial stress, and the relative deformation of the transition zones. Orthogonal tests with 6 factors and 5 levels were then designed for the important structural parameters of the dog-bone shaped sheet tensile specimens. According to the rules of the orthogonal test design, 25 dynamic tensile specimens with different structural dimensions were obtained. The commercial finite element software ABAQUS/Explicit was used to establish a finite element model of the split Hopkinson tensile bar, and dynamic tensile test simulations were performed on the dynamic tensile specimens obtained from the orthogonal test design. An orthogonal test dataset with the specimen structural parameters as the input and the measurement accuracy indicators as the output was then constructed. Multi-objective orthogonal test matrix analysis was carried out on the orthogonal test dataset to obtain the influence order as well as the influence law of the structural parameters of the tensile specimens on the measurement accuracy indicators of tests. Taking the orthogonal test dataset as the training dataset, an artificial neural network (ANN) model was used to fit the nonlinear relationship which can predict the measurement accuracy indicators of the test by using the structural parameters of the specimen, and then the fitness function in the genetic algorithm (GA) was established by using this model. Finally, the structural parameters of the dynamic tensile specimen were optimized using the ANN-GA collaborative optimization method, and the optimal structural dimensions of the dynamic tensile specimen were obtained as the result of the optimization. Finite element simulation results show that the optimal structural dimensions obtained by the ANN-GA optimization method are valid. The results of this study demonstrate the practicability and effectiveness of the ANN-GA method in the structural optimization of dynamic tensile specimens. On the other hand, it can provide guidance for the specimen design in the dynamic tensile mechanical properties tests of materials, and can also provide a reference for the validity analysis of the experimentally measured mechanical properties.
-
Key words:
- split Hopkinson bar /
- geometry effect /
- specimen design /
- machine learning /
- genetic algorithm
-
近年来,世界范围内意外爆炸事故频发,给人民的生命财产安全造成了严重危害。例如,2015年天津市滨海新区爆炸事故,2020年黎巴嫩贝鲁特港口爆炸事故,2021年湖北省十堰市燃气爆炸事故等,均对当地建筑物造成了严重破坏。为了保护生命和防止工程结构发生严重损坏,使用合理的防护手段对于缓解爆炸荷载对结构的危害至关重要[1-3]。
通过在结构表面设置牺牲包层来缓解爆炸荷载对结构的冲击,在近些年引起了学者们广泛的研究兴趣。Wu等[4]研究了泡沫铝包层对钢筋混凝土板的爆炸毁伤缓解性能;Rebelo等[5]探究了3D打印聚乳酸蜂窝结构作为牺牲包层可压碎芯时的吸能特性;Bohara等[6]研究了蜂窝夹芯板在近场和远场爆炸下对钢筋混凝土板的保护作用;范东宇等[7]对强动载荷下多孔泡沫牺牲层的动态压溃行为及缓冲吸能机理进行了研究;Zhao等[8]将泡沫水泥基材料用作隧道衬砌的牺牲包层并对其防护效果进行了评估。这些研究成果表明,在牺牲包层的选择上,大多更倾向于多孔固体结构。
非水反应发泡聚氨酯(polyurethane, PU)作为一种内部为泡孔结构、质量轻且力学性能优良的高聚物材料,也具有充当牺牲包层、缓解和抵抗爆炸荷载的潜力。目前,在静态力学领域,因其具有反应迅速可调节、抗压强度高、早强、耐腐蚀等[9]优点,已被广泛应用于地下混凝土管道脱空[10]、堤坝渗漏[11]、高速铁路无砟轨道沉降[12]等工程病害的非开挖修复中,均呈现出良好的修护效果。另一方面,在爆炸与冲击领域,高聚物材料的优良力学性能也引起了学者们的广泛关注。Wang等[13-14]通过落锤冲击试验和准静态压缩试验探究了填充聚氨酯泡沫吸能连接器的失效机理及能量吸收特性;Jamil等[15]研究了热塑性聚氨酯夹芯板作为牺牲包层的吸能效率;张勇[16]对聚氨酯泡沫铝复合结构进行了接触爆炸试验,探讨了聚氨酯泡沫铝的吸能性能。然而,现有研究成果主要集中于冲击、爆炸作用下聚氨酯及其复合结构的吸能特性与失效机理,而关于将非水反应发泡聚氨酯材料设置于结构表面充当牺牲包层时的防护性能、特别是针对钢筋混凝土结构的爆炸毁伤缓解效应的研究还比较少。
本文中,对带高聚物牺牲包层钢筋混凝土板(polyurethane-reinforced concrete slab, PU-RCS)开展接触爆炸试验,对比分析高聚物牺牲包层对钢筋混凝土板的爆炸毁伤缓解效应。在此基础上,建立现场试验的SPH-FEM耦合模型,参数化分析炸药量、高聚物牺牲包层密度、厚度对PU-RCS毁伤模式和吸能特性的影响。研究结果以期为新型非水反应发泡聚氨酯在工程结构抗爆防护领域的进一步研究与应用提供相关参考。
1. 试验概况
图1为PU试件的制备流程图。图2所示的试验所用PU牺牲包层是由异氰酸酯和多元醇按1∶1的质量比混合后反应生成的,制备的PU试件的尺寸为500 mm×500 mm×60 mm,密度为0.2 g/cm3,该密度的PU材料被认为具有良好的静态力学性能,且内部泡孔发育充分,其杨氏模量为35.81 MPa,弹性极限为1.54 MPa,屈服强度为2.04 MPa。
牺牲包层的整体能量吸收能力主要是由其应力应变-曲线中的屈服平台决定的。通过向特制的模具中注浆,共制备了密度为0.2 g/cm3、边长为70.7 mm的3个立方体试件,利用电液伺服万能试验机(WAW-300型)分别进行单轴压缩试验,加载速率为0.2 mm/min。通过3组重复性的单轴压缩试验,排除了试件自身所带缺陷和其他偶然性因素对测试结果的影响,从而得到了图3所示的PU试件的代表性应力-应变曲线。从图3中可以看到,非水反应发泡聚氨酯的应力-应变曲线主要分为3个阶段:坍塌前阶段、平台阶段和致密化阶段。在图3中,δy和εy分别为屈服应力和屈服应变,εd为致密化应变。不难发现,在屈服应变εy和致密化应变εd之间存在着一个相当长的屈服平台,这表明非水反应发泡聚氨酯具有较强的能量吸收能力,具备充当牺牲包层的潜力。
钢筋混凝土板的几何形状和钢筋分布如图4所示,其设计尺寸为500 mm×500 mm×60 mm。所选混凝土材料的强度等级为C40,所用钢筋的强度等级为HRB400。钢筋的直径为8 mm,采用单层双向布筋,钢筋的保护层厚度为20 mm。
接触爆炸试验装置如图5所示。在试验开始前,首先将钢筋混凝土板安装在特制的钢框架上,板的两边嵌入钢框架的卡槽内,然后将高聚物板覆盖在钢筋混凝土板的顶面充当牺牲包层,通过螺栓使PU-RCS处于两端固定的边界条件下。试验采用密度为1.05 g/cm3,爆速为4.2~5.0 km/s的岩石乳化炸药。称取50 g炸药置于高聚物牺牲包层迎爆面的中心处,药包形状采用球形装药。将普通非电导爆管雷管插入球形药包中心作为起爆装置,用于引爆炸药。在对照组中,炸药直接放置于钢筋混凝土板的中央,除此之外,其他试验条件均相同。
2. 数值模型
由于现场爆炸试验具有极高的危险性和试验成本,因此进行爆炸试验的机会是宝贵的。随着计算机技术的快速发展,数值模拟已成为研究爆炸荷载作用下结构动态响应的重要手段之一。本文中基于AUTODYN软件建立了接触爆炸作用下PU-RCS的三维精细化耦合模型,并通过数值模型深入探究高聚物牺牲包层对钢筋混凝土板的爆炸缓解效应。
2.1 SPH-FEM耦合模型
建立的SPH-FEM耦合模型如图6所示,耦合模型和现场爆炸试验中对应结构的物理尺寸相同。在SPH-FEM耦合模型中,混凝土由Lagrange网格建模,钢筋由Beam单元建模,炸药和高聚物牺牲包层由SPH粒子建模。Lagrange和Beam单元的网格尺寸均为5 mm,混凝土由120 000个Lagrange实体单元组成,钢筋由752个Beam单元组成。构成炸药和高聚物牺牲包层的SPH粒子的直径均为5 mm,其中50 g炸药所包含的SPH粒子数为344,而高聚物牺牲包层由120 000个SPH粒子组成。混凝土与钢筋之间通过共享节点连接,并且假设它们完全粘合,没有任何滑移,PU-RCS的两端均施加固定边界。
2.2 材料模型
2.2.1 混凝土
采用RHT(Riedel-Hiermaier-Thoma)动态损伤模型对混凝土材料进行建模[17]。该模型中采用了弹性极限面、失效面及残余强度面作为3个控制破坏面以描述混凝土材料的初始屈服强度、失效强度和残余强度。
RHT模型的失效面方程为:
Yfail=Y(p∗)R(θ)F(˙ε) (1) Y(p∗)=fc{A[p∗−p∗spallF(˙ε)N]} (2) 式中:Yfail为混凝土的失效强度, Y(p*)为压缩子午线等效应力强度,fc为单轴抗压强度,A为失效面常数,N为失效面指数,p*为由fc归一化的压应力,
p∗spall=ft/fc ,ft为单轴抗拉强度,F(˙ε) 为应变率的函数,R(θ)为任意应力角对应的子午线半径和压缩子午线半径之比,θ为罗德角。RHT模型的弹性极限面方程为:
Yelastic=YfailFcap(p)Felastic (3) 式中:Yelastic为混凝土的初始屈服强度,Felastic为弹性强度与失效面强度之比,Fcap(p)为盖帽函数,用于限制静水压力下的弹性偏应力。
RHT模型的残余失效面方程为:
Y∗res=Bp∗M (4) 式中:Y*res为混凝土的残余强度,B为残余强度面常数,M为残余强度面指数。具体材料参数见表1。
表 1 混凝土材料模型主要参数Table 1. Parameters of the concrete material model剪切模量G/GPa 体积模量A1/GPa 抗压强度fc/MPa ft/fc fs/fc 失效面常数A 16.7 35.27 40 0.06 0.18 1.6 失效面指数N 残余失效面常数B 残余失效面指数M 损伤常数D1 损伤常数D2 侵蚀应变 0.61 1.6 0.61 0.04 1.0 2.0 2.2.2 钢筋
采用Johnson-Cook模型[1]来描述钢筋的动态行为。该模型将屈服应力Y定义为:
Y=(A+Bεnp)(1+Clnε∗p)(1−TmH) (5) TH=(T−Troom)/(Tmelt−Troom) (6) 式中:
εp 为等效塑性应变,ε∗p 为归一化有效塑性应变率,A、B、C、m、n为材料常数,Tmelt为钢筋熔化温度,Troom为初始温度。具体材料参数见表2。表 2 钢筋材料模型主要参数Table 2. Parameters of the reinforcement steel密度ρ/(g·cm−3) 体积模量K/GPa 剪切模量G/GPa A/MPa B/MPa C m n Troom/K Tmelt/K 7.83 159 81.8 404 232.4 0.014 1.03 0.26 300 1793 2.2.3 炸药
采用状态方程[18-19]模拟炸药爆炸过程中压力和内能及爆轰产物相对体积之间的关系:
p=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωEV (7) 式中:p为爆轰压力,V为爆轰产物的相对体积,A、B、R1、R2和ω为材料常数,E为炸药体积内能。具体参数见表3。
表 3 炸药材料模型主要参数Table 3. Parameters of the explosiveρ/(g·cm−3) E/(GJ·m−3) A/GPa B/GPa R1 R2 ω 1.05 6.0 209.7 3.5 5.76 1.29 0.39 2.2.4 高聚物
采用Crushable foam模型[16]对高聚物材料进行建模。该模型需要输入密度、弹性模量、剪切模量、 最大拉伸应力以及应力-应变曲线等材料参数。为了提高计算效率,采用了瞬时几何应变侵蚀准则,当单元的动态响应达到给定的准则时,将自动从模型中删除畸变单元。本文中对不同的几何应变进行了试算,发现使用0.3的几何应变作为侵蚀标准,可以在保证计算不中断的同时获得可靠的结果。具体材料参数见表4。
表 4 高聚物材料模型主要参数Table 4. Parameters of the polymer material model密度/(g·cm−3) 屈服强度/MPa 最大拉伸应力/MPa 剪切模量/MPa 杨氏模量/MPa 侵蚀应变 0.2 2.04 1.77 18.7 35.8 0.3 3. 试验与数值结果
3.1 试验结果
图7为炸药量为50 g的接触爆炸下,PU-RCS的迎爆面和背爆面的毁伤结果。对于PU-RCS试件,将高聚物牺牲包层移走并清理碎屑后,发现其表现出了较轻的损伤程度,迎爆面中心区域没有明显的爆坑,只有深度极浅的混凝土压碎区域,其直径约为10 cm;背爆面中部出现了直径约13 cm的剥落区,但是混凝土剥落深度极浅。另一个明显的特征是背爆面分布着大约13条向四周发散的细长裂纹。
图8所示为对照组的试验结果。与PU-RCS的毁伤结果不同,RCS试件迎爆面中部有一个直径约为15 cm的圆形爆坑,最大深度约为1 cm;底面剥落区直径约为25 cm,最大深度约为4 cm。剥落区的外围分布着沿径向的发散状裂纹。从毁伤结果来看,接触爆炸下RCS试件发生了比较严重的局部毁伤。
表5详细地总结了上述试验后相应试件的毁伤情况。通过两者的毁伤结果对比表明,在相同的炸药量下,相较于RCS试件,PU-RCS的局部毁伤程度大大降低,这是由于高聚物牺牲包层有效地发挥了其缓解爆炸荷载的能力。
表 5 试验后试件的毁伤结果Table 5. Damage resuls of the specimens after the tests试件 炸药量/g 牺牲包层 板面 损伤区域直径/cm 损伤区最大深度/cm 裂纹数量/条 PU-RCS 50 有 迎爆面 10 <0.5 0 背爆面 13 <0.8 13 RCS 50 无 迎爆面 15 1 0 背爆面 25 4 10 为更加直观形象地描述接触爆炸下高聚物牺牲包层对钢筋混凝土板毁伤模式的影响,图9给出了炸药量为50 g时PU-RCS的毁伤模式示意图。炸药起爆后,首先是高聚物牺牲包层经历了压碎、破裂等大范围毁伤变形,这消耗了相当大的一部分爆炸能量,当剩余的能量作用于钢筋混凝土板时,其破坏作用已相当小,因此钢筋混凝土板的迎爆面只产生了毁伤深度极浅的压碎区,背爆面混凝土的剥落范围和深度也大大减小,随之产生的是充分发展的发散状裂纹。不难发现,高聚物牺牲包层发挥了其缓解爆炸荷载的能力,在一定程度上改变了钢筋混凝土板的毁伤模式,使得钢筋混凝土板的毁伤程度明显降低。
3.2 模型验证
采用适当的网格尺寸对于精确的模拟至关重要。为此,进行网格收敛性分析,以确定数值模型中的网格尺寸。建立了3个耦合模型,对PU-RCS分别采用10、5和4 mm的网格尺寸进行建模计算,装药量为50 g,其他条件均相同。图10所示为接触爆炸荷载下网格尺寸对钢筋混凝土板总能量的影响。从图10中可以看出不同网格尺寸下总能量的变化趋势非常相似,后两种网格尺寸对应的总能量峰值很相近,而前一种网格尺寸对应的总能量峰值较大,后两种网格的模拟结果相当一致,表明了网格的收敛性。进一步减小网格尺寸可能会得到更精确的数值结果,但会大大增加计算时间。因此,采用5 mm网格尺寸的平板模型是可行的。
为了验证所建立模型的准确性,将数值模拟结果与试验结果进行对比。如图11和表6所示。
表 6 模拟结果和试验结果参数对比Table 6. Parameters comparison of the simulation and experimental results方法 迎爆面破坏直径/cm 背爆面破坏直径/cm 迎爆面最大毁伤深度/cm 背爆面最大毁伤深度/cm 背爆面裂纹数量 试验 10 13 <0.5 <0.8 13 模拟 10.7 12.5 <0.6 <0.5 16 吻合度 93.46% 96.15% − − 81.25% 从图11和表6中可以看出,数值模拟结果不仅再现了带高聚物牺牲包层钢筋混凝土板迎爆面中心区域的压碎区和底面的层裂损伤,且底面发散性裂纹的扩散规律也基本保持一致。
图12为试验和数值模拟中高聚物牺牲包层的损伤结果对比。在现场试验中,高聚物牺牲包层首当其冲地遭受了接触爆炸荷载的冲击,这导致牺牲包层发生了比较严重的大范围毁伤,破碎情况较为严重。同时,数值模拟结果也非常准确地再现了高聚物牺牲包层的大范围压碎和贯穿性损伤。
对比以上数值模拟和试验的结果,二者的毁伤特征吻合度很高,本文中所建立的SPH-FEM耦合模型能够的准确地模拟接触爆炸下带高聚物牺牲包层钢筋混凝土板的毁伤,其结果合理并可靠。
3.3 PU-RCS动态响应过程
本文中模拟的爆炸粒子与PU-RCS的相互作用过程如图13所示。可以看到,在接触爆炸发生的瞬间(0.15 ms内),爆炸粒子快速向四周发散,球形炸药的体积急剧膨胀。显而易见的是,位于钢筋混凝土板上方的高聚物牺牲包层,首当其冲地受到了爆轰产物和高压冲击波的剧烈冲击。在t=0.01 ms时,高聚物牺牲包层已经发生了冲击破坏;在t=0.04 ms时,高聚物牺牲包层的中部已经发生贯穿破坏。正如在t=0.06, 0.15 ms时,在炸药粒子飞散的同时,高聚物牺牲包层的毁伤程度和范围不断增加,直到最终高聚物牺牲包层发生了非常严重的整体破碎,这也预示着高聚物牺牲包层已经基本完成了其所谓的“牺牲”任务,发挥了其抵抗并缓解爆炸荷载冲击的作用。
图14~15分别给出了应力波在RCS和PU-RCS的跨中截面处的传播过程。对于RCS,炸药爆炸后,产生的高压冲击波立即从装药位置向外传播,爆轰产物和高压冲击波直接作用在钢筋混凝土板的迎爆面上。在两者的猛烈作用下,装药周围的混凝土所受压应力急剧升高,使迎爆面中部混凝土被压碎成坑。压力传播到板内部后将以压缩波的形式向底面传播,在板背爆面压缩波反射为了拉伸波,并与压缩波相互作用,合成的拉伸应力高于混凝土的动态抗拉强度,从而导致背爆面混凝土发生震塌、剥落破坏。
而对于PU-RCS,因为高聚物牺牲包层的存在,产生的爆轰产物和高压冲击波将首先作用于牺牲包层,冲击波在高聚物内部的泡孔结构中发生反射、折射,同时高聚物牺牲包层的屈服变形能力缓解了爆轰产物和高压冲击波的剧烈冲击作用。二者与牺牲包层的相互作用过程消耗了相当大一部分能量和较长的时间,最终导致作用在钢筋混凝土板上的爆轰产物和冲击波压力明显减小,钢筋混凝土板的动态响应时刻明显延后,同时爆炸荷载被分散到除板中心外的其他板周区域。
与RCS相比,在应力波的传播过程中,PU-RCS顶面附近压应力值不仅明显降低,压应力作用区域也明显扩大,这说明高聚物牺牲包层的存在缓解了钢筋混凝土板中部的压应力集中效应,减小了爆炸瞬间产生的爆轰产物、高压冲击波对钢筋混凝土板中部的局部冲量作用。
4. 带高聚物牺牲包层钢筋混凝土板 (PU-RCS) 参数分析
基于所建立SPH-FEM耦合模型,探究了炸药量和高聚物牺牲包层密度、厚度对PU-RCS上下表面毁伤特性以及跨中位移、能量吸收变化趋势的影响。
4.1 炸药量
炸药量分别设置为0.10、0.13、0.16、0.19和0.22 kg,其余参数均保持不变,分析不同药量下PU-RCS的毁伤特性。图16为不同药量下PU-RCS上下表面的毁伤结果。随着药量的增加,板顶面中部的压碎区不断变大,产生的环形裂纹不断向顶面四周扩展,裂纹的密集程度不断加深。底面的发散性径向裂纹不断向周围扩展,且由于混凝土的抗拉强度较低,而作用于底面的拉应力却随着药量的增加而不断变大,因此板中心区域的混凝土剥落区也不断变大。当装药量达到0.16 kg时,无论是顶面还是底面,几乎都布满了裂纹。此后,随着药量从0.16 kg增加到0.22 kg,牺牲包层防护能力被明显削弱,板上下表面的毁伤程度大大加深。
图17所示为不同药量下PU-RCS的跨中截面损伤。从截面的角度看,随着药量的增加,PU-RCS的跨中位移和横向的毁伤范围不断变大。结合图18,当药量为0.10 kg时,PU-RCS跨中位移仅为15.6 mm,而当药量增加到0.16 kg时,PU-RCS已经出现了比较明显的弯曲破坏,此时跨中位移增加到32.2 mm,当药量继续增加至0.22 kg时,PU-RCS已经发生了比较严重的整体弯曲破坏,跨中位移达到42.5 mm。
图19为不同药量下高聚物牺牲包层(PU)和被保护钢筋混凝土板(RCS)的吸能曲线。随着起爆药量的增加,PU的吸能量不断增加,当药量为0.10 kg时,PU的吸能量约为7.4 kJ,而当药量最后增加到0.22 kg时,PU的吸能量已经达到15.6 kJ,增加了一倍以上;相比于PU,RCS的吸能量虽然也在增加,但其吸能量在4 kJ以下,始终维持在较低水平。
图20为不同药量下PU和RCS的吸能百分比。随着药量的增加,对于PU-RCS的总吸能量来说,PU的吸能量占比保持在80%以上,这意味着当PU-RCS遭受爆炸荷载冲击时,高聚物牺牲包层可以吸收大部分爆轰能量,使被保护钢筋混凝土板所受冲击显著削弱。也正是因为高聚物牺牲包层的吸能特性以及对接触爆炸荷载的分散作用,使得被保护钢筋混凝土板发生了明显的整体弯曲破坏。
4.2 高聚物牺牲包层密度
非水反应高聚物材料的密度范围一般在0.15~0.50 g/cm3。基于4.1节中的5种不同药量,将高聚物牺牲包层密度分别设置为0.20、0.25、0.30、0.35和0.40 g/cm3,分析不同高聚物牺牲包层密度下PU-RCS的毁伤特性。由于篇幅限制,这里只给出几组具有代表性的工况。图21是当药量为0.16 kg,高聚物牺牲包层密度分别为0.20、0.25、0.30、0.35和0.40 g/cm3时PU-RCS上下表面的毁伤结果图。可以看到,随着高聚物牺牲包层密度的不断增大,PU-RCS的顶面中部的压碎区面积不断减小,同时上下表面裂纹的密集程度有所降低,但底面中部的混凝土剥落区的面积变化并不明显。
图22为不同牺牲包层密度下PU-RCS的跨中截面处的毁伤结果。可以发现,随着高聚物牺牲包层密度的增大,PU-RCS跨中截面处的毁伤程度发生了比较明显的变化。结合图23可知,当牺牲包层密度仅为0.20 g/cm3时,PU-RCS发生了严重的整体弯曲破坏,跨中位移达到32.2 mm。当密度增大到0.30 g/cm3时,PU-RCS的整体弯曲破坏程度大大减轻,跨中位移减小到22.2 mm。最终,当密度增大到0.40 g/cm3时,PU-RCS几乎没有发生明显的弯曲破坏,跨中位移仅为15.1 mm。
以上结果表明,密度直接影响着高聚物牺牲包层对钢筋混凝土板的防护性能。相同药量下,随着高聚物牺牲包层密度的增大,PU-RCS的跨中位移明显减小,但减小的趋势逐渐变缓。
图24为上述5组不同牺牲包层密度下PU和RCS的吸能曲线。在相同药量下,随着密度的增大,PU的吸能量不断增多,当包层密度为0.20 g/cm3时,PU的吸能量约为12.1 kJ,当密度增大到0.40 g/cm3时,PU的吸能量达到18.0 kJ;另一方面,随着牺牲包层密度的增大,RCS的吸能量逐渐减少,这说明增大牺牲包层密度对缓解爆炸能量对RCS的传递和输入能起到正向积极作用。
图25为不同牺牲包层密度下PU和RCS的吸能百分比。相同药量下,随着牺牲包层密度的增大,PU吸能量在PU-RCS的总吸能量中所占比例不断增加,当密度为0.40 g/cm3时,PU吸能占比达到90.2%。从吸能曲线和两部分所占吸能百分比可以发现,牺牲包层密度是影响其爆炸缓解能力的重要因素,在一定范围内增大密度,有利于提高牺牲包层的能量吸收能力,使输入到被保护结构的能量明显降低,进而提升防护效果。
4.3 高聚物牺牲包层厚度
将牺牲包层密度保持为0.20 g/cm3不变,通过改变牺牲包层厚度,探究不同厚度对PU-RCS毁伤特性的影响。图26给出了当药量为0.16 kg,牺牲包层厚度分别为4、5、6、7和8 cm时,PU-RCS上下表面的毁伤结果。随着包层厚度的增大,钢筋混凝土板顶面的裂纹密集程度逐渐降低,板上下表面中部的局部毁伤面积明显减小,增大牺牲包层厚度能有效降低被保护结构的毁伤程度。
图27为不同牺牲包层厚度下,PU-RCS的跨中损伤结果。可以发现,当包层厚度为4 cm时,顶面中部混凝土出现明显的爆坑,底面中部混凝土因抗拉强度较低,几乎处于剥落状态,此时整个板块与接触爆炸下普通钢筋混凝土板的典型毁伤模式类似,以局部严重毁伤为主;而当厚度增大到8 cm时,板块不仅没有发生严重的局部毁伤,而且整体弯曲破坏也不明显,板的跨中位移很小。由此可见,不同于牺牲包层密度的变化,包层厚度的变化会改变钢筋混凝土板的毁伤模式,当包层厚度较小时,钢筋混凝土板以冲切成坑、剥落等典型局部毁伤为主;而当包层厚度较大时,钢筋混凝土板的损伤模式主要是大范围裂纹,整体弯曲破坏程度不明显。值得一提的是,因为包层厚度为5和7 cm时钢筋混凝土板的跨中截面毁伤模式分别与包层厚度为4和8 cm时的毁伤模式很接近,所以为了避免重复,这里并未给出。
图28为不同牺牲包层厚度下PU和RCS的吸能曲线。对于PU来说,随着其厚度的增大,吸能量逐渐增多,但是增长的幅度并不大,当厚度为4 cm时,吸能量约为10.9 kJ,厚度增加到8 cm时,吸能量只增加到了13.5 kJ,增长率仅为23.8%;而对于RCS来说,随着包层厚度的增加,输入到RCS的能量下降比较明显,当包层厚度从4 cm增加到8 cm时,输入到RCS的爆炸能量从5.6 kJ减小到1.4 kJ,减少了75%。另外,不同牺牲包层厚度下PU和RCS的吸能百分比如图29所示,可以看出增加高聚物牺牲包层的厚度能明显提升PU的吸能百分比。
5. 结 论
通过现场试验和数值模拟相结合的方法,探究了接触爆炸作用下带高聚物牺牲包层钢筋混凝土板 (PU-RCS) 的毁伤特征,对比分析了高聚物牺牲包层对钢筋混凝土板的爆炸毁伤缓解效应。参数化分析了炸药量和高聚物牺牲包层密度、厚度对PU-RCS毁伤特征及吸能特性的影响。可以得到以下结论。
(1) 高聚物牺牲包层可以有效地改变接触爆炸下钢筋混凝土板的毁伤模式。高聚物牺牲包层对爆炸荷载的缓解和分散作用减轻了板中部的应力集中效应,使钢筋混凝土板的局部毁伤程度明显减小。板中部的爆坑与混凝土剥落破坏转变为充分发展的发散性裂纹,在一定程度上减少了爆炸碎片的飞散。这表明高聚物牺牲包层在工程结构的抗爆防护中有着潜在的应用前景。
(2) 通过数值模拟得到的毁伤结果与现场爆炸试验结果基本一致,表明本文中所建立的SPH-FEM耦合模型能够有效地模拟接触爆炸下PU-RCS的毁伤特性。参数分析表明,随着药量的增大,跨中位移不断增大,板出现了明显的整体弯曲破坏现象;增大高聚物牺牲包层的密度和厚度有利于提高其防护性能,降低被保护钢筋混凝土板的弯曲破坏程度,但包层厚度的变化会引起被保护钢筋混凝土板毁伤模式的改变。
(3) 高聚物牺牲包层具有优良的吸能特性,药量在一定范围内增大时,高聚物牺牲包层依然能维持较高的吸能水平,吸能占比达到80%以上;牺牲包层密度、厚度的增大可以明显地增加包层的吸能量,提高包层的吸能占比,进而降低爆炸能量对被保护钢筋混凝土板的作用。
-
表 1 片状动态拉伸试样结构参数及其惯用值
Table 1. Structural parameters of the sheet specimen used for dynamic tensile tests and the reference values
lg/mm wg/mm lc/mm wc/mm R/mm δ/mm 12.0 4.0 18.0 16.0 3.0 2.0 表 2 试样结构参数正交表与模拟结果
Table 2. Orthogonal table of the specimen structural parameters and the simulation results
序号 结构参数/mm t/μs ˉσm/% s2/10−3 ˉDd/% lg wg lc wc R δ 01 6.0 2.0 16.0 14.0 0.5 1.0 18.05 4.55 2.241 0.57 02 6.0 3.0 17.0 15.0 1.0 1.5 18.50 6.11 2.176 3.54 03 6.0 4.0 18.0 16.0 2.0 2.0 19.00 5.90 1.776 7.23 04 6.0 5.0 19.0 17.0 3.0 2.5 18.00 5.76 1.874 11.44 05 6.0 6.0 20.0 18.0 4.0 3.0 19.50 6.17 1.412 15.58 06 8.0 2.0 18.0 17.0 1.0 3.0 18.50 2.94 1.930 2.01 07 8.0 3.0 19.0 18.0 2.0 1.0 19.50 2.59 0.865 4.83 08 8.0 4.0 20.0 14.0 3.0 1.5 20.00 3.33 0.765 8.47 09 8.0 5.0 16.0 15.0 4.0 2.0 20.00 4.29 0.753 12.17 10 8.0 6.0 17.0 16.0 0.5 2.5 21.01 8.42 4.844 1.90 11 10.0 2.0 20.0 15.0 2.0 2.5 22.00 1.68 0.807 3.73 12 10.0 3.0 16.0 16.0 3.0 3.0 22.00 2.63 0.651 6.54 13 10.0 4.0 17.0 17.0 4.0 1.0 23.12 1.91 0.414 9.65 14 10.0 5.0 18.0 18.0 0.5 1.5 21.99 6.32 3.203 1.19 15 10.0 6.0 19.0 14.0 1.0 2.0 21.00 6.88 2.925 2.34 16 12.0 2.0 17.0 18.0 3.0 2.0 23.50 1.39 0.357 5.17 17 12.0 3.0 18.0 14.0 4.0 2.5 24.00 1.99 0.316 8.14 18 12.0 4.0 19.0 15.0 0.5 3.0 21.32 5.21 1.972 1.06 19 12.0 5.0 20.0 16.0 1.0 1.0 21.99 4.58 1.852 1.70 20 12.0 6.0 16.0 17.0 2.0 2.0 23.00 4.46 1.651 4.18 21 14.0 2.0 19.0 16.0 4.0 1.5 26.00 0.94 0.146 6.99 22 14.0 3.0 20.0 17.0 0.5 2.0 24.51 3.53 1.295 0.74 23 14.0 4.0 16.0 18.0 1.0 2.5 25.11 3.93 1.205 1.77 24 14.0 5.0 17.0 14.0 2.0 3.0 25.03 3.95 0.805 4.00 25 14.0 6.0 18.0 15.0 3.0 1.0 25.50 3.39 0.750 5.18 表 3 模拟用材料模型参数
Table 3. Parameters of the material model used for simulations
材料 密度/(kg·m−3) 弹性模量/GPa 泊松比 Johnson-Cook 模型参数 A/MPa B/MPa N C 45钢 7 850 210 0.30 AA7075-T6 2 800 71 0.33 473 210 0.3813 0.033 表 4 试样结构参数对测量精度指标影响的主次顺序及规律
Table 4. The influence order and law of the specimen structural parameters on the measurement accuracy indicators
指标 主次顺序 关键因素 规律 t lg, R, lc, δ, wg, wc lg、R lg、R与t呈正相关 s2 R, lg, wg, δ, wc, lc R、lg R、lg与s2呈负相关 ˉσm wg, R, lg, δ, wc, lc wg、R wg与ˉσm呈正相关,R与ˉσm呈负相关 ˉDd R, lg, wg, δ, lc, wc R、lg R与ˉDd呈正相关,lg与ˉDd呈负相关 表 5 验证集
Table 5. Validation data set
编号 结构参数/mm t/μs ˉσm/% s2/10–3 ˉDd/% lg wg lc wc R δ 01 6.0 3.0 18.0 15.0 2.0 1.5 16.0 4.45 1.018 7.02 02 8.0 2.0 16.0 14.0 3.0 1.0 18.5 1.24 0.275 7.59 03 10.0 4.0 19.0 17.0 1.0 2.5 21.0 5.43 1.946 2.34 表 6 验证集真实值与预测值的对比
Table 6. Comparison of the true and predicted values of the validation dataset
编号 t ˉσm s2 ˉDd 真实值/μs 预测值/μs 误差/% 真实值/% 预测值/% 误差/% 真实值/10–3 预测值/10–3 误差/% 真实值/% 预测值/% 误差/% 01 16.00 15.62 2.38 4.45 4.87 9.44 1.018 1.072 5.30 7.02 6.73 4.13 02 18.50 17.36 6.16 1.24 1.37 10.48 0.275 0.261 5.09 7.59 8.18 7.77 03 21.00 20.78 1.05 5.43 4.99 8.10 1.946 1.778 8.63 2.34 2.43 3.85 表 7 优化前后试样测量准确度指标对比
Table 7. Comparison of the measurement accuracy indicators of specimens before and after optimization
t ˉσm s2 ˉDd 优化前/μs 优化后/μs 减小/% 优化前/% 优化后/% 减小/% 优化前/10−3 优化后/10−3 减小/% 优化前/% 优化后/% 减小/% 27.0 21.0 22.2 2.4 1.57 34.6 0.871 0.803 7.8 5.57 5.16 7.9 -
[1] 蔡明, 陈伟, 陈利强, 等. TC8钛合金的动态力学性能及本构关系 [J]. 机械工程材料, 2020, 44(12): 80–84. DOI: 10.11973/jxgccl202012015.CAI M, CHEN W, CHEN L Q, et al. Dynamic mechanical properties and constitutive relation of TC8 titanium alloy [J]. Materials for Mechanical Engineering, 2020, 44(12): 80–84. DOI: 10.11973/jxgccl202012015. [2] 杨建兴, 张江涛, 乔炎亮, 等. 硅橡胶动态压缩性能SHPB测试方法及其本构模型研究 [J]. 高分子通报, 2021(4): 27–34. DOI: 10.14028/j.cnki.1003-3726.2021.04.004.YANG J X, ZHANG J T, QIAO Y L, et al. SHPB test method for the dynamic compressive properties and dynamic constitutive model of silicon rubber [J]. Polymer Bulletin, 2021(4): 27–34. DOI: 10.14028/j.cnki.1003-3726.2021.04.004. [3] FENG F, HUANG S Y, MENG Z H, et al. A constitutive and fracture model for AZ31B magnesium alloy in the tensile state [J]. Materials Science and Engineering: A, 2014, 594: 334–343. DOI: 10.1016/j.msea.2013.11.008. [4] LIU Z, DONG Y, MAO P L, et al. Dynamic tensile and compressive properties of vacuum and ordinary die-casting AT72 magnesium alloy at high strain rates [J]. Journal of Magnesium and Alloys, 2013, 1(2): 150–162. DOI: 10.1016/j.jma.2013.07.004. [5] OGAWA K. Impact-tension compression test by using a split-Hopkinson bar [J]. Experimental Mechanics, 1984, 24(2): 81–86. DOI: 10.1007/BF02324987. [6] HARDING J, WOOD E O, CAMPBELL J D. Tensile testing of materials at impact rates of strain [J]. Journal of Mechanical Engineering Science, 1960, 2(2): 88–96. DOI: 10.1243/JMES_JOUR_1960_002_016_02. [7] HUH H, KANG W J, HAN S S. A tension split Hopkinson bar for investigating the dynamic behavior of sheet metals [J]. Experimental Mechanics, 2002, 42(1): 8–17. DOI: 10.1007/BF02411046. [8] VERLEYSEN P, DEGRIECK J, VERSTRAETE T, et al. Influence of specimen geometry on split Hopkinson tensile bar tests on sheet materials [J]. Experimental Mechanics, 2008, 48(5): 587–589. DOI: 10.1007/s11340-008-9149-x. [9] VERLEYSEN P, BENEDICT V, VERSTRAETE T, et al. Numerical study of the influence of the specimen geometry on split Hopkinson bar tensile test results [J]. Latin American Journal of Solids and Structures, 2009, 6(3): 285–298. DOI: 10.1007/s12205-009-0371-6. [10] XU W F, HUANG X C, HAO Z M, et al. Effect of the geometric shapes of specimens on impact tensile tests [J]. Journal of Zhejiang University: Science A, 2010, 11(10): 817–821. DOI: 10.1631/jzus.A1000139. [11] SUN X, SOULAMI A, CHOI K S, et al. Effects of sample geometry and loading rate on tensile ductility of TRIP800 steel [J]. Materials Science and Engineering: A, 2012, 541: 1–7. DOI: 10.1016/j.msea.2011.12.115. [12] ROTBAUM Y, RITTEL D. Is there an optimal gauge length for dynamic tensile specimens? [J]. Experimental Mechanics, 2014, 54(7): 1205–1214. DOI: 10.1007/s11340-014-9889-8. [13] PRABOWO D A, KARIEM M A, GUNAWAN L. The effect of specimen dimension on the results of the split-Hopkinson tension bar testing [J]. Procedia Engineering, 2017, 173: 608–614. DOI: 10.1016/j.proeng.2016.12.114. [14] NGUYEN K H, LEE C W, SHIN H, et al. A study on the effects of specimen geometry on measurement accuracy of dynamic constitutive properties of metals using SHTB [J]. International Journal of Precision Engineering and Manufacturing, 2020, 21(9): 1687–1695. DOI: 10.1007/s12541-020-00368-y. [15] 郭伟国, 李玉龙, 索涛. 应力波基础简明教程 [M]. 西安: 西北工业大学出版社, 2007: 131−132. [16] BAYRAMIN B, ŞIMŞIR C, EFE M. Dynamic strain aging in DP steels at forming relevant strain rates and temperatures [J]. Materials Science and Engineering: A, 2017, 704: 164–172. DOI: 10.1016/j.msea.2017.08.006. [17] JIANG B H, CAO L B, ZHU F. Dynamic tensile behavior of polypropylene with temperature effect [J]. Composites Part B: Engineering, 2018, 152: 300–304. DOI: 10.1016/j.compositesb.2018.08.133. [18] BARDELCIK A, SALISBURY C P, WINKLER S, et al. Effect of cooling rate on the high strain rate properties of boron steel [J]. International Journal of Impact Engineering, 2010, 37(6): 694–702. DOI: 10.1016/j.ijimpeng.2009.05.009. [19] SMERD R, WINKLER S, SALISBURY C, et al. High strain rate tensile testing of automotive aluminum alloy sheet [J]. International Journal of Impact Engineering, 2005, 32(1/2/3/4): 541–560. DOI: 10.1016/j.ijimpeng.2005.04.013. [20] 徐仲安, 王天保, 李常英, 等. 正交试验设计法简介 [J]. 科技情报开发与经济, 2002, 12(5): 148–150. DOI: 10.3969/j.issn.1005-6033.2002.05.084.XU Z A, WANG T B, LI C Y, et al. Brief introduction to the orthogonal test design [J]. Sci/Tech Information Development & Economy, 2002, 12(5): 148–150. DOI: 10.3969/j.issn.1005-6033.2002.05.084. [21] 陈建华, 周晨佳, 王雪, 等. 基于正交试验的高速井泵优化设计 [J]. 排灌机械工程学报, 2021, 39(5): 457–463. DOI: 10.3969/j.issn.1674-8530.19.0227.CHEN J H, ZHOU C J, WANG X, et al. Optimization design of high-speed well pump based on orthogonal experiment [J]. Journal of Drainage and Irrigation Machinery, 2021, 39(5): 457–463. DOI: 10.3969/j.issn.1674-8530.19.0227. [22] DAVIES E D H, HUNTER S C. The dynamic compression testing of solids by the method of the split Hopkinson pressure bar [J]. Journal of the Mechanics and Physics of Solids, 1963, 11(3): 155–179. DOI: 10.1016/0022-5096(63)90050-4. [23] VERLEYSEN P, DEGRIECK J. Non-homogeneous and multi-axial stress distribution in concrete specimens during split Hopkinson tensile tests [J]. Computers & Structures, 2000, 77(6): 669–676. DOI: 10.1016/S0045-7949(00)00022-5. [24] SONG B, CHEN W. Dynamic stress equilibration in split Hopkinson pressure bar tests on soft materials [J]. Experimental Mechanics, 2004, 44(3): 300–312. DOI: 10.1007/BF02427897. [25] ZHANG D N, SHANGGUAN Q Q, XIE C J, et al. A modified Johnson-Cook model of dynamic tensile behaviors for 7075-T6 aluminum alloy [J]. Journal of Alloys and Compounds, 2015, 619: 186–194. DOI: 10.1016/j.jallcom.2014.09.002. [26] 周玉珠. 正交试验设计的矩阵分析方法 [J]. 数学的实践与认识, 2009, 39(2): 202–207.ZHOU Y Z. A matrix analysis of orthogonal design [J]. Mathematics in Practice and Theory, 2009, 39(2): 202–207. [27] 何龙, 张冉阳, 赵刚要, 等. 基于BP神经网络的GH5188高温合金本构模型 [J]. 特种铸造及有色合金, 2021, 41(2): 223–226.HE L, ZHANG R Y, ZHAO G Y, et al. Constitutive model of GH5188 superalloy based on BP neural network [J]. Special Casting & Nonferrous Alloys, 2021, 41(2): 223–226. [28] NARAYANA P L, LEE S W, PARK C H, et al. Modeling high-temperature mechanical properties of austenitic stainless steels by neural networks [J]. Computational Materials Science, 2020, 179: 109617. DOI: 10.1016/j.commatsci.2020.109617. [29] 王小平, 曹立明. 遗传算法: 理论、应用与软件实现 [M]. 西安: 西安交通大学出版社, 2002: 1−15. 期刊类型引用(2)
1. 孔祥清,李若男,常雅慧,付莹. 泡沫填充负泊松比蜂窝夹层结构的抗爆性能数值模拟. 兵工学报. 2024(09): 3091-3104 . 百度学术
2. 李子玉,丁建国,周广盼. 基于ALE-FEM-SPH算法的混凝土板在爆炸荷载下响应分析. 南京理工大学学报. 2024(06): 773-780 . 百度学术
其他类型引用(0)
-