A remote sensing imagery-based model for assessment of building damage induced by large-equivalent explosions
-
摘要: 为了研究大当量爆炸建筑物毁伤评估问题,基于遥感影像解译和大数据分析构建了大当量爆炸建筑物毁伤评估模型。首先,基于大当量爆炸的具体历史案例构建了毁伤数据集,具体指基于遥感影像提取建筑物毁伤信息,辅助大数据信息补充毁伤细节,利用地理信息系统空间分析数字化毁伤信息,构成毁伤数据集。然后,基于毁伤数据集中的训练样本修正经验模型参数,构建了适用于大当量爆炸的针对不同类型建筑物的毁伤评估模型,并基于毁伤数据集中的验证样本测试了模型性能。实验证明:所构建模型拟合优度高于96%,检验样本准确度高于84%,整体误差在可接受范围内。所构建模型在一定精度要求下可为大当量爆炸事故评估提供参考。Abstract: To address challenges in the field of large-scale explosive building damage assessment, where the explosion process is too complex for high-precision numerical simulation, and relying solely on change detection from remote sensing imagery cannot capture detailed internal information and lacks the capability of predicting in advance, this paper establishes a building damage assessment model for large-scale explosive events by coupling empirical mechanics models with remote sensing image interpretation and big data analysis. The study initially constructs a damage dataset based on specific historical cases of large-scale explosions. This involves extracting building damage information (including building types and damage levels) from remote sensing imagery and supplementing damage details with additional big data sources such as collected online images, videos, and news reports to enhance the precision of the sampled data. Geographic information systems spatial analysis is employed to digitize the damage information, obtaining data on building types, damage levels, and the distance from the target building to the explosion center, forming the damage dataset. Subsequently, the empirical model parameters are refined based on the training samples from the damage dataset, creating damage assessment models applicable to different building types for large-scale explosive events. The performance of the model is then tested using validation samples from the damage dataset. Experimental results demonstrate a model fitting goodness of over 96%, accuracy on validation samples exceeding 84%, and an overall error within an acceptable range. The model, under certain accuracy requirements, can provide guidance for site selection of storage locations for chemicals and hazardous materials, emergency evacuation of people in the event of a risk of large-scale explosions, critical equipment evacuation during an emergency, resource dispatching for rescue and relief after an accident, and building damage assessment.
-
Key words:
- large equivalent explosion /
- remote sensing imagery /
- damage assessment /
- big data analysis /
- building
-
薄壁结构质量轻并且具有较好的能量吸收特性,常被用来生产各种能量吸收或缓冲装置,已被广泛应用于航天航空、车辆工程及国防装备等领域。薄壁结构在轴向冲击载荷作用下通过胞壁的塑性变形吸收大量的冲击能量,所以研究薄壁结构的轴向压缩特性具有重要的工程意义。
McFarland[1]最先提出正六边形多胞结构轴向准静态平均压缩应力的计算方法,Wierzbicki[2]修正了McFarland提出的基本折叠模式,并根据超折叠单元理论[3]对正六边形多胞结构的轴向平均压缩应力和折叠波长进行了理论求解;Tran等[4]通过理论和数值手段获得了多种薄壁管在冲击载荷作用下的轴向压缩应力理论模型;尹汉锋等[5]基于超折叠单元理论对几种胞元构型蜂窝的平均压缩应力进行求解,并开展耐撞性优化设计。
薄壁结构的几何构型对其轴向冲击载荷作用下的吸能特性有着较大的影响。为提高普通构型(三角形、正六边形、圆形、正方形等)薄壁管的吸能能力,许多学者对一些新型层级薄壁结构产生了兴趣。Sun等[6]将规则六边形多胞结构的每个顶点替换为一个较小的六边形结构,并分析了其轴向冲击载荷作用下的耐撞性能。Mousanezhad等[7]提出了具有高刚度和高韧性的蜘蛛网层级多胞结构,并研究了结构参数对其面内及面外力学性能的影响。Sun等[8]用同性子结构替换了规则六边形多胞结构的实心胞壁,研究发现该层级特性的引入能大幅度提高其面内刚度。张越等[9]运用数值模拟方法研究了二阶自相似四边形蜂窝结构参数对其面外动态压缩性能的影响。于国际等[10]研究了二阶层级六边形蜂窝的面内动态压缩性能。赖燕辉等[11]基于多级蜂窝构型法分析了层级蜂窝的弹性模量等力学参数。
受Sierpinski三角形的启发,本文中将Sierpinski分形特性引入薄壁吸能管的层级设计,从而提出一种具有Sierpinski层级特性的新型薄壁多胞管(Sierpinski hierarchical tube, SHT),并对SHTs在轴向冲击作用下的变形模式和吸能特性进行数值模拟分析。进一步基于能量守恒原理对SHTs在轴向压缩下的平均压缩应力理论模型进行求解,以期为新型薄壁构型轴向缓冲吸能装置的设计提供指导。
1. SHTs的层级结构设计
分形结构可以通过多种方法产生,Sierpinski[12]于1916年提出了一个典型的自相似集,并命名为Sierpinski三角形。如图1所示,Sierpinski垫是一个由若干个三角形构成的自相似几何结构。
基于Sierpinski[12]的研究工作,本文研究的具有Sierpinski层级特性的多胞管是通过胞壁连接规则三角形薄壁管胞壁中点形成的。在每个SHT中有3i个边长为
li 的单位三角形,通过重复这一迭代过程能够得到更高阶的SHTs,i表示层级数,同时普通三角形薄壁管可以定义为零阶SHTs。第零~三阶SHTs的结构示意图如图1所示。由图1中(i−1)阶SHTs发展为i阶SHTs的典型过程可知,层次结构参数可以定义为:
γ=lili−1 (1) 式中:
li 和li−1 表示每个单位三角形的胞元壁长。值得注意的是,SHTs的结构参数γ将保持在0.5不变。l0 为普通三角形薄壁管的胞壁长。因此,i阶SHTs的所有胞壁长度之和为:Li=3i+12il0 (2) 本文中假设这些结构的胞壁厚度是均匀的。作为薄壁结构的重要结构参数,根据面积比理论,i阶SHTs的表观密度可以表示为:
ρi=4√3⋅(32)itil0ρs (3) 式中:
ti 为i阶SHTs中单位三角形的胞壁厚度,ρs 为基体材料密度。因此,SHTs的相对密度由下式给出:
¯ρi=ρiρs=4√3⋅(32)itil0 (4) 构成SHTs的单位三角形薄壁长度随着层级参数i的增大而变短,这就会给数值模拟带来很大的麻烦,因此本文中只考虑到三阶SHTs。本文中
l0 为90 mm,多胞管轴向长度h为250 mm。考虑3组壁厚,即t0 =1.0,1.3,1.6 mm,薄壁管的相对密度分别为0.077、0.100和0.123。由于Sierpinski层级结构特征的引入,本文中所研究的SHTs结构变得更复杂。然而,SHT作为一种典型的薄壁结构,其仍可以被划分为若干个典型单元。如图2所示,所有的SHTs均可看作由2种基本单元组成,即V形角单元和K形角单元。
2. SHTs的轴向压缩模拟分析
2.1 SHTs轴向压缩有限元模型
利用通用显式非线性有限元分析软件LS-DYNA模拟轴向冲击载荷作用下SHTs的动态压缩特性,计算模型如图3所示。构建上下端面刚性板,当顶端刚性板沿轴向以10 m/s的恒定速度冲击薄壁试件时,底端刚性板固定。为了准确地模拟薄壁管的大变形,胞壁采用Beltschko-Tsay四边形壳单元,单元厚度方向采用五点积分,面内采用单点积分。基体材料为铝合金AA6060T4,材料的杨氏模量
EY=68.2GPa ,屈服应力σy=80MPa ,极限应力σu=173MPa ,密度ρs=2700kg/m3 ,泊松比μ=0.3 ,幂指强化因数n=0.23 [13]。由于铝合金对应变率不敏感,本文中未考虑应变率的影响。薄壁管胞壁采用自动单面接触算法来考虑自身变形产生的接触;薄壁管与刚性板之间采用自动点-面接触算法。数值模型中的静摩擦因数及动摩擦因数均取0.2。2.2 有限元模型验证
为验证该模型的可靠性,首先对薄壁方形管[14]的压缩行为进行模拟。通过有限元计算后处理可以提取冲击端应力应变曲线,名义应力
σ 可表示为:σ=F(d)/S (5) 式中:
F(d) 为冲击端压力,d为冲击端压缩距离,S为薄壁结构所占面积。应变
ε 可表示为:ε=d/h (6) 动态平均压缩应力
σdm 则表示为:σdm=1εd∫εd0σdε (7) 式中:
εd 为密实化应变。计算出应力σdm 后,动态平均压缩力则可表示为:Pdm=σdm⋅S (8) 模拟结果与实验数据[14]的比较见表1,误差均小于10%,表明模拟结果与实验结果吻合较好。
2.3 变形模式和塑性机理
图4给出了
ˉρ=0.077 的不同SHTs结构在轴向冲击载荷作用下不同时刻的典型变形图。由于多胞管的底部被固定,胞壁均从冲击端开始发生塑性变形,随着冲击端的不断下压,胞元逐层发生折叠变形。研究结果还表明,具有Sierpinski层级特性的多胞管在轴向压缩过程中均呈现出轴对称渐进屈曲模式。轴对称渐进屈曲是理想的屈曲模式,已有研究表明薄壁结构在受到轴向冲击载荷作用时也多发生轴对称的向外屈曲变形[2-5],并且塑性变形主要集中在管壁的连接处。与此同时,Sierpinski层级特征的引入导致胞元节点段数的增加。为保持结构的相对密度不变,胞壁厚度会变小,从而在某种程度上削弱了薄壁结构的弯曲变形能。与规则三角形管相比,引入Sierpinski分级特性的多胞管变形吸能最显著的区别就是折叠单元的数量。SHTs完全折叠单元的数量会随着层级数的增加而增加。Sierpinski分级特性的引入大大缩短了胞壁变形的半折叠波长,这也就意味着在压缩过程中更多的塑性变形能被耗散,增强了薄壁结构的抗压缩能力。
Abramowicz等[15]指出,薄壁管在塑性变形过程中有2类基本的变形单元,即非延展性基本单元和延展性基本单元,如图5所示。图5中2H为折叠单元高度,图5(a)中延展吸能区域阴影面积为S*,图5(b)中延展吸能区域阴影面积为S#。
前文介绍SHTs结构时提及到作为一种典型的薄壁结构,其仍可以被划分为若干个典型单元,即如图2所示,所有的SHTs均可看作由2种基本单元组成,即V形角单元和K形角单元。
图6给出了2种角单元的典型变形轮廓图,通过对变形模式的仔细观察能得到SHT结构在整个变形过程中的2种典型变形模式。这些变形模式主要取决于相邻胞壁的相对运动。对于V形角单元,在每一个折叠单元变形过程中,2个胞壁都向同一方向移动,其变形轮廓如图6(a)所示,这属于传统的非延展性变形模式(Ⅰ型)。
对于K形单元,4个胞壁在折叠变形过程中均向外移动,如图6(b)所示,此时相邻胞壁间的变形机理属于延展性变形模式(Ⅱ型)。
3. SHTs的轴向压缩应力计算
图7为薄壁结构的胞壁受轴向冲击载荷作用时的变形示意图,其中H被称为胞元折叠单元半波长塑性铰长度,
δ 为轴向压缩长度,q 为冲击载荷。轴对称渐进屈曲是理想的屈曲模式,已有研究[2-5]表明薄壁结构在受到轴向冲击载荷作用时也多发生轴对称渐进屈曲变形。学者们基于简化超折叠单元(simplified super folded element, SSFE)理论求解了大量该变形模式下薄壁管的塑性坍塌问题[3-5]。
SHTs是一种典型的薄壁结构,通过数值模拟和分析可知其在压缩载荷作用下发生轴向渐进屈曲模式变形,这就启发我们运用SSFE理论来推导SHTs的轴向压缩力。为确定SHTs的轴向压缩应力,还需作如下假设:
(1)薄壁管材料塑性好,可视为理想刚塑性材料;
(2)如图7所示,薄壁管在轴向压缩载荷作用下发生向外的轴对称渐进屈曲,变形过程中各个折叠的塑性铰长度相等,为2H;
(3)在变形过程中,同一胞元上下表面保持平行,即胞元的各边变形量相同;
(4)压缩过程中胞壁间粘接强度足够大,不发生破裂,可忽略粘接对薄壁结构力学性能的影响。
考虑到冲击能量主要由薄壁结构胞壁的塑性变形耗散,根据能量守恒原理,外力做功等于薄壁结构塑性变形所耗散的能量,即弯曲耗能Ebending和拉伸耗散能
Emembrane 。相应的能量平衡方程可表示为:Pmη(2H)=Ebending+Emembrane (9) 式中:
Pm 为准静态平均压缩力。事实上,折叠单元在塑性变形过程中无法完全被压实。因此,这里给出了考虑有效压缩的修正因数。η 为有效压缩距离δe 与半折叠波长H 的比值。Wierzbicki等[15]提出η 可能在0.70~0.75之间变化。为简单起见,本文中采用η=0.70 。3.1 SHTs的弯曲耗能计算
Chen等[16]指出通过将折叠单元的水平铰线的弯曲耗散能量相加,可以估算出压缩过程中的全部弯曲耗能,即:
Ebending=j∑i=1M0φiLi (10) 式中:
M0 为胞元单位长度的塑性极限弯矩,M0=σ0t2i/4 ;φi 表示第i 条铰线的旋转角度,j 表示铰线的数量;Li 是所有胞壁的总长度;σ0 为蜂窝基体材料的流动应力[13]。σ0 的表达式为:σ0=√σyσu/(1+n) (11) 式中:
σy 和σu 分别为SHTs基体材料的屈服应力和极限应力,n 为硬化指数。如图8所示,假设薄壁管胞壁被完全压实,水平铰线的弯曲角度即为
π/2 、π 和π/2 。因此,弯曲耗散能可表示为:Ebending=2πM0Li (12) 3.2 角单元的拉伸耗散能
3.2.1 V形角单元
图6给出了简化超折叠单元理论中的基本折叠单元延展吸能示意图,本文中考虑了SSFE理论中的2种拉伸变形模式。图5(a)所示为变形模式,延展吸能区由3个三角形阴影部分组成,通过对阴影面积积分可求得该变形模式下单个胞壁的拉伸耗散能:
Easymmembrane_f=∫S∗σ0tidS∗=12σ0tiH2=2M0H2ti (13) Chen等[16]指出2个互相连接的胞壁对角单元的拉伸形变能具有相似的贡献。因此,在非延展变形的情况下,直角单元的拉伸耗散吸能可表示为单个角折叠单元拉伸耗散吸能的2倍:
Easymmembrane_r=2Easymmembrane_f=4M0H2ti (14) Tran等[4]在角单元的拉伸耗散吸能方面做了大量的研究,指出与直角单元相比,V形角单元的拉伸耗散吸能较小。如图9所示,V形角单元在完全塑性坍塌过程中的拉伸耗散吸能可表示为:
EV-shapedmembrane=Easymmembrane_rcosα=4M0H2ticosα (15) 式中:
α 为V形角单元两胞壁的夹角。图9中β 为V形角单元胞壁与直角单元胞壁之间的夹角。3.2.2 K形角单元
通过对变形轮廓图的分析,不难发现K形角单元的变形机理比直角角单元的变形机理复杂。将K形角单元的薄膜耗散吸能简化为由一个V形角单元和两个附加平面的薄膜耗散吸能组成。K形角单元中的每个胞壁都具有相似的变形模式。K形角单元中包含的V形角单元与单个V形角单元具有相同的抗压强度,需注意的是根据仿真结果(图7)可知此时的V形角单元处于延展性变形模式,如图6(b)所示。阴影部分为延展吸能区,通过积分阴影区域来评估完全塌陷期间每个胞壁的薄膜耗散吸能:
Esymmembrane_f=∫S#σ0tidS#=σ0tiH2=4M0H2ti (16) 图10给出了在K形角单元的拐角处形成的延展吸能三角形单元。K形角单元中包含的V形角单元的薄膜耗散能可表示为:
EV-shaped′membrane=2Esymmembrane_fcosβ=8M0H2ticosβ (17) 另外两个附加胞壁的薄膜耗散吸能为:
E2-additionalmembrane=2Esymmembrane_f=8M0H2ti (18) 因此,在一个折叠波长压缩期间,K形角单元的薄膜耗散吸能由一个V形角单元和两个附加胞壁的薄膜耗散吸能组成:
EK-shapedmembrane=EV-shaped′membrane+E2-additionalmembrane=8M0H2ti(1+1cosβ) (19) 3.3 SHTs的动态平均压缩应力计算
普通三角形薄壁管,即零阶SHTs是由3个V形角单元组合而成。将式(12)和(15)代入式(9),平均压缩力
Pm,0th 的理论计算方程可写为:Pm,0th(2H)η=Ebending+3EV-shapedmembrane=2πM0L0+2M0H2t0(6cosα) (20) 对式(20)进行转换可得:
Pm,0thηM0=πL0H+Ht0(6cosα)=πL0H+Ht0F(α) (21) 根据准静态压缩过程中的能量最低原理可得
∂Pm,0th/∂H=0 ,对式(21)求导得:0=−πL0H2+F(α)t0⇒H=√πL0t0F(α)=√0.518πl0t0 (22) 将式(22)代入式(20),可求得普通三角形薄壁管的准静态平均压缩载荷为:
Pm,0th=πM0L0ηH+M0Hηt0F(α)=π0.5σ0t1.50L0.50√F(α)2η=√17.387π2ησ0t1.50l0.50 (23) 式中:
F(α)=6cosα 。而薄壁管所占面积为:
S=√3l20/4 (24) 所以,普通三角形薄壁管的轴向准静态平均压缩应力为:
σm,0th=Pm,0th/S=2√5.796πησ0(t0l0)3/2 (25) 同理,一阶SHTs可看作由3个V形角单元和3个K形角单元组合而成;二阶SHTs可看作由3个V形角单元和12个K形角单元组合而成;三阶SHTs可看作由3个V形角单元和39个K形角单元组合而成。对3种SHTs的轴向准静态平均压缩应力进行求解,可得轴向准静态平均压缩应力分别为:
σm,1st=2√62.693πησ0(t1l0)3/2 (26) σm,2nd=2√337.04πησ0(t2l0)3/2 (27) σm,3rd=2√1599.06πησ0(t3l0)3/2 (28) 研究表明动态压缩应力一般比相应的准静态压缩应力高[5],而上文的理论推导过程中忽略了动态压缩(包括惯性和应变率)的影响。考虑到铝的应变速率效应可以忽略不计,采用增强因子
λ 来考虑惯性效应。Hanssen等[17]指出λ可在1.3~1.6范围内变化。为简化计算,取10 m/s冲击速度下的增强因子为1.3。因此,可以得到零~三阶SHTs的动态压缩应力理论计算值分别为:σdm,0th=2λ√5.796πησ0(t0l0)3/2 (29) σdm,1st=2λ√62.693πησ0(t1l0)3/2 (30) σdm,2nd=2λ√337.04πησ0(t2l0)3/2 (31) σdm,3rd=2λ√1599.06πησ0(t3l0)3/2 (32) 3.4 理论模型模拟验证
本文中
t0 分别取1.0、1.3和1.6 mm,胞壁长l0 取90 mm,分别对这些新型SHT多胞结构进行轴向压缩的数值模拟。用式(29)~(32)计算出动态压缩应力的理论值,并将其与模拟结果进行比较,如表2所示,当
t0=1.0 mm(ˉρ=0.077 )时,一阶、二阶及三阶SHTs的动态压缩应力较普通三角形薄壁管分别升高了85.8%、138.2%和183.8%。这表明:将Sierpinski层级特性应用到薄壁管的设计中,可以有效提高薄壁管的吸能特性。4种薄壁管的动态压缩应力模拟值和理论值的差异分别为−2.64%~3.42%、−5.23%~1.54%、−2.77%~3.14%和−1.19%~3.67%,结果基本吻合,进一步验证了本文动态压缩应力理论推导方法的可行性,说明理论计算结果具有工程应用价值,可用于指导新型薄壁结构轴向缓冲吸能装置的设计。表 2 各采样点有限元结果与理论预测值对比Table 2. Comparison of finite element results and theoretical predictions for all the sampling points层级级数 ˉρ ti/mm li/mm σdm/MPa 相对误差/% 模拟 理论 0th 0.077 1.00 90 1.925 1.968 −2.18 0.100 1.30 2.840 2.917 −2.64 0.123 1.60 4.118 3.982 3.42 1st 0.077 0.67 45 3.577 3.523 1.54 0.100 0.87 4.948 5.221 −5.23 0.123 1.07 7.046 7.129 −1.17 2nd 0.077 0.44 22.5 4.586 4.446 3.14 0.100 0.58 6.407 6.590 −2.77 0.123 0.71 9.186 8.998 2.09 3rd 0.077 0.30 11.25 5.464 5.271 3.67 0.100 0.39 7.720 7.813 −1.19 0.123 0.47 10.972 10.669 2.84 4. 结 论
将Sierpinski分形结构引入到薄壁吸能管的层级设计中,提出了一种具有Sierpinski层级特性的新型薄壁多胞管(SHT),并通过数值模拟和理论方法研究其受轴向冲击载荷作用下的变形模式和能量吸收特性,研究结果表明:
(1)具有Sierpinski层级特性的SHTs在轴向压缩过程中均呈现出理想的轴对称渐进屈曲模式。SHTs完全折叠单元的数量会随着层级数的增加而增加。Sierpinski分级特性的引入大大缩短了胞壁变形的半折叠波长,压缩过程中更多的塑性变形能量被耗散,增强了薄壁结构的抗压缩能力。
(2)SHTs均可看作由两种基本单元组成,即V形角单元和K形角单元。V形角单元在每一个折叠单元变形过程中,两个胞壁都向同一方向移动,属于传统的非延展性变形模式(Ⅰ型)。K形角单元4个胞壁在折叠变形过程中均向外移动,相邻胞壁间的变形机理属于延展性变形模式(Ⅱ型)。
(3)基于能量守恒理论和塑性铰理论建立了SHTs的塑性屈曲理论模型,获得了压溃过程中的塑性弯曲耗散能、薄膜耗散能、轴向压缩应力的理论表达式,模拟结果与理论计算结果基本吻合,验证了本文中动态压缩应力理论推导方法是可行的,理论计算结果具有工程应用价值。
(4)在相同的相对密度下,一阶、二阶及三阶SHTs的动态压缩应力较普通三角形薄壁管分别增高了85.8%、138.2%和183.8%。将Sierpinski层级特性引入到薄壁管的设计中,能够有效提高薄壁管的耐撞性能,这可为新型吸能元件的研究和设计提供参考。
-
表 1 建筑物目视解译毁伤等级标准
Table 1. Standard for visual interpretation of damage grade of buildings
建筑物类型
与毁伤等级毁伤前影像 毁伤后影像 标准细则 砖混结构
中度毁伤(1)屋顶和楼板出现明显破损或塌陷,
形成明显洞口或凹陷;
(2)墙体出现明显裂缝;
(3)结构柱明显变形,墙体倾斜砖混结构
重度毁伤(1)屋顶或楼板出现严重破损,楼层部分
或全部垮塌;
(2)墙体部分或全部崩塌;
(3)结构柱断裂或严重变形,墙体严重倾斜
或移位工业厂房
中度毁伤(1)屋顶破损但整体结构未受到影响;
(2)钢结构轻微变形或扭曲,但并未断裂
或明显变形;
(3)主体结构轻微倾斜或不平衡工业厂房
重度毁伤(1)屋顶严重破损;
(2)钢结构出现明显脱离、断裂
或腐蚀等情况;
(3)主体结构受到严重破坏,部分
或整体倒塌表 2 样本详情
Table 2. Sample details
样本 事故详情 发生时间 经纬度 TNT当量/t 爆炸前成像时间 爆炸后成像时间 训练样本1 河北省张家口市盛华化工有限公司爆炸 2018.11.28 40.75°N 115.00°E 0.5 2018.11 2018.12 训练样本2 山东省青岛市中石化东黄输油管道爆炸 2013.11.22 36.05°N 120.21°E 3 2013.03 2014.02 训练样本3 美国德克萨斯州韦斯特镇的韦斯特化工厂爆炸 2013.04.17 31.81°N 97.09°W 10 2012.10 2013.04 训练样本4 河南省三门峡市义马气化厂爆炸 2019.07.19 34.74°N 111.84°E 30 2019.03 2020.04 训练样本5 江苏省盐城市陈家港镇工业园区爆炸 2019.08.04 34.32°N 119.79°E 260 2018.03 2019.04 训练样本6 黎巴嫩贝鲁特港口爆炸 2020.08.04 33.90°N 35.52°E 620 2020.07 2020.08 检验样本1 江苏省连云港市堆沟港镇化工园爆炸 2017.12.09 34.42°N 119.78°E 10 2017.09 2018.03 检验样本2 天津市滨海新区天津港爆炸 2015.08.12 39.04°N 117.74°E 450 2015.05 2015.08 表 3 训练样本数据
Table 3. Training sample data
样本 TNT当量/t 砖混结构中度毁伤 砖混结构重度毁伤 工业厂房中度毁伤 工业厂房重度毁伤 y/m ŷ/m p/% y/m ŷ/m p/% y/m ŷ/m p/% y/m ŷ/m p/% 1 0.5 45 38 84.44 30 28 93.33 2 3 141 70 49.65 50 45 90.00 3 10 192 129 67.19 105 82 78.10 263 171 65.02 154 124 80.52 4 30 252 223 88.49 180 142 78.89 343 296 86.30 195 214 90.26 5 260 733 655 89.36 468 417 89.10 980 872 88.98 735 631 85.85 6 770 1065 1128 94.08 677 717 94.09 1417 1500 94.12 1025 1085 94.15 R2/% 96.97 97.80 97.74 97.86 注:y为原始值,ŷ为预测值;p为准确度,p=1−|相对误差|;
拟合优度由R2度量,R2=1−残差平方和/总平方和=1−Σ(原始值−预测值)2/Σ(原始值−原始值均值)2。表 4 检验样本数据
Table 4. Test sample data
样本 TNT当量/t 砖混结构中度毁伤 砖混结构重度毁伤 工业厂房中度毁伤 工业厂房重度毁伤 y/m ŷ/m p/% y/m ŷ/m p/% y/m ŷ/m p/% y/m ŷ/m p/% 1 10 130 129 99.23 96 82 85.42 157 171 91.08 107 124 86.29 2 450 895 862 96.31 473 548 84.14 988 1147 83.91 716 830 86.27 pmean/% 97.77 84.78 87.49 86.28 注:pmean为平均准确度。 -
[1] 张连玉, 汪令羽, 苗瑞生. 爆炸气体动力学基础 [M]. 北京: 北京工业学院出版社, 1987: 425–426. [2] 周子钦. 基于多时相图像的打击效果评估技术研究 [D]. 武汉: 华中科技大学, 2020: 40–45. DOI: 10.27157/d.cnki.ghzku.2020.006392.ZHOU Z Q. Study of damage effect assessment based on multi-temporal image [D]. Wuhan: Huazhong University of Science and Technology, 2020: 40–45. DOI: 10.27157/d.cnki.ghzku.2020.006392. [3] 张军. 多时相图像检测方法及其在毁伤评估系统中的应用 [D]. 上海: 上海交通大学, 2008: 17–18.ZHANG J. Processing method of multi-temporal remote sensing images and the use in damage assessment [D]. Shanghai: Shanghai Jiao Tong University, 2008: 17–18. [4] SIRMACEK B, UNSALAN C. Damaged building detection in aerial images using shadow information [C]//Proceedings of the 4th International Conference on Recent Advances in Space Technologies. Istanbul: IEEE, 2009: 249–252. DOI: 10.1109/RAST.2009.5158206. [5] 王威. 基于图像理解的打击效果评估系统研究与实现 [D]. 武汉: 华中科技大学, 2015: 9–13.WANG W. Research and implementation of battle damage assessment system based on image understanding [D]. Wuhan: Huazhong University of Science and Technology, 2015: 9–13. [6] 勾涛. 基于图像分析的毁伤评估系统关键技术研究 [D]. 长春: 吉林大学, 2019: 5–6.GOU T. Research on key technologies of damage assessment system based on image analysis [D]. Changchun: Jilin University, 2019: 5–6. [7] 孔祥锡, 秦闻远, 苏飘逸, 等. 基于深度学习及模糊层次分析的毁伤评估算法 [J]. 航空学报, 2023, 40: 1–18.KONG X X, QIN W Y, SU P Y, et al. Damage assessment algorithm based on deep learning and fuzzy analytic hierarchy process [J]. Acta Aeronautica et Astronautica Sinica, 2023, 40: 1–18. [8] 中国人民解放军总装备部军事训练教材编辑工作委员会. 核爆炸物理概论 [M]. 北京: 国防工业出版社, 2003: 51–55. [9] 李秦超, 姚成宝, 程帅, 等. 神经网络状态方程在强爆炸冲击波数值模拟中的应用 [J]. 爆炸与冲击, 2023, 43(4): 044202. DOI: 10.11883/bzycj-2022-0222.LI Q C, YAO C B, CHENG S, et al. Application of the neural network equation of state in numerical simulation of intense blast wave [J]. Explosion and Shock Waves, 2023, 43(4): 044202. DOI: 10.11883/bzycj-2022-0222. [10] GLASSTONE S, DOLAN P J. The effects of nuclear weapons [R]. USA: Defense Technical Information Center, 1977: 453–501. DOI: 10.21236/ada087568. [11] ZHANG H M, WANG M C, ZHANG Y X, et al. TDA-Net: a novel transfer deep attention network for rapid response to building damage discovery [J]. Remote Sensing, 2022, 14(15): 3687. DOI: 10.3390/rs14153687. [12] 孙家抦. 遥感原理与应用 [M]. 3版. 武汉: 武汉大学出版社, 2013: 174–176.SUN J B. Principles and applications of remote sensing [M]. 3rd ed. Wuhan: Wuhan University Press, 2013: 174–176. [13] 全国地震标准化技术委员会. 中华人民共和国地震行业标准: DB/T 75—2018 [S]. 北京: 中国标准出版社, 2019. [14] 卢芳云, 李翔宇, 田占东, 等. 武器毁伤与评估 [M]. 北京: 科学出版社, 2021: 93–152. 期刊类型引用(7)
1. 黄翠萍,邓小林,蔡珍珍. 自相似嵌套金属多胞管耐撞性研究. 机械强度. 2024(03): 568-578 . 百度学术
2. 杨永宝,奚杰,李志超. 阶梯多胞方管横向弯曲响应特性研究. 农业装备与车辆工程. 2024(12): 98-105 . 百度学术
3. 蔡玮雯,马其华,甘学辉. 低速冲击载荷下CFRP-Al多胞薄壁管的耐撞性. 塑性工程学报. 2023(04): 187-196 . 百度学术
4. 覃善甘,邓小林,刘鑫炎,刘夫云. 新型层级梯度六边形管的耐撞性多目标优化. 塑性工程学报. 2023(07): 194-204 . 百度学术
5. 孔志成,胡俊,郭智平. 星形混合多胞管的耐撞性数值与理论研究. 爆炸与冲击. 2023(08): 43-55 . 本站查看
6. 邓敏杰,刘志芳. 仿马尾草薄壁结构的设计与耐撞性研究. 高压物理学报. 2022(03): 111-120 . 百度学术
7. 贺璞,邓庆田,李新波. 层合多孔圆柱壳的轴向冲击吸能特性. 高压物理学报. 2022(04): 145-154 . 百度学术
其他类型引用(12)
-