A theoretical method for calculating spatial periodic distribution of deformation of a spherical shell under explosive loading
-
摘要: 早期研究提出了对振动叠加应变增长现象的解剖式分析方法,进而发现爆炸加载下带扰动源球壳上的弯曲波和壳体变形呈空间周期分布的规律。参考Timoshenko梁的弯曲理论,基于平截面假定和壳体发生较小的弯曲变形的假设,推导出球壳上弯曲波波速和波长的关系,计算得到最短弯曲波和与膜振动频率相近的弯曲波的波速,还结合早期研究提出的壳体变形分布周期与弯曲波波速的关系,计算得到了壳体变形空间分布的周期。结果表明:(1)理论计算结果与数值仿真结果基本吻合,其中弯曲波波速的计算结果与数值仿真结果相差在15%以内,壳体变形空间分布周期的计算结果与数值仿真结果相差在12%以内;(2)弯曲波波长越短,波速越快,当波长无限短时,波速趋于极限值,约为声速的0.574倍。本计算方法为解剖式分析方法提供了一定的理论依据。Abstract: The strain growth, caused by vibration superposition, has been anatomized by the membrane strain and the bending strain in existing studies, and the bending wave and deformation spatial periodic distribution of a spherical shell under explosive loading have been found. By referring to the theoretical method for Timoshenko beam bending, based on a plane-section assumption and a small-deformation limit, the relation between the velocity and the wavelength of bending wave was deduced, and the velocities of the shortest bending wave and the bending wave with a frequency similar to that of the membrane vibration were calculated. By combining the relation between the deformation spatial distribution period and the bending wave velocity presented in existing studies, the deformation spatial distribution period was calculated. The main conclusions are as follows: (1) The theoretical results are in good agreement with the numerical results, in which the difference between the numerical and theoretical results of bending wave velocity is within 15%, and the difference between the numerical and theoretical results of the deformation spatial distribution period is within 12%. (2) The shorter the wavelength, the higher the wave velocity, when the wavelength is infinite short, the bending wave velocity tends to the limit value, about 0.574 times the speed of sound. The theoretical method presented in this paper provides a certain theoretical support for anatomizing strain growth.
-
Key words:
- explosive vessel /
- strain growth /
- bending wave /
- spatial periodic distribution
-
1976年Buzukov[1]实验发现了爆炸容器弹性变形范围内的应变增长现象,在这种现象中容器壳体的最大变形没有出现在第一个振动周期内,而是出现在后续周期内。应变增长现象对容器安全来说是不利的,因为传统的容器设计方法一般以第一个周期内的变形峰值为依据。近几年,应变增长现象的研究有不少新发现。Dong等[2]、Li等[3]发现壳体不稳定性导致弯曲振动进而形成应变增长现象,并研究了相关机理;刘文祥等[4-5]实验观察到壳体塑性变形下的应变增长现象,获得了应变增长系数超过6的数据,甚至推测存在系数超过10的情况。
振动叠加是形成应变增长现象的重要原因[2, 6-7]。Duffey等[7]基于振动力学中“拍”形成机理,结合壳体振动频率分析结果以及出现应变增长现象的应变曲线的特征,推测频率相近的振动叠加形成应变增长现象。Liu等[8]提出了对应变增长现象解剖式分析的方法,针对带扰动源球壳上的应变增长现象,分离出弯曲应变(弯曲波引起的壳体弯曲变形导致的应变)和膜应变(壳体拉伸变形导致的应变),直接证明了频率相近的膜振动和弯曲振动叠加是导致应变增长现象的主要原因,还发现了球壳变形呈空间周期分布。
解剖式分析是研究振动叠加形成应变增长现象的新方法,早期研究虽提出了壳体弯曲波、变形空间周期性分布的相关机理,但缺乏理论的支撑。本文中参考Timoshenko梁的弯曲理论[9-10],理论分析球壳上弯曲波的波速以及壳体变形空间分布的周期,并与数值仿真结果进行对比,以验证理论计算的可靠性。
1. 球壳上的弯曲波和变形空间周期分布
球形爆炸容器往往存在人孔门、观察窗等功能部件,在爆炸作用下这些部件与球壳的运动存在差异,由扰动源处产生弯曲波并在球壳上传播,这些部件被称为扰动源。完全约束扰动源是指绝对静止状态的扰动源,其导致的弯曲波比其他约束条件扰动源的更强,导致的应变增长现象也更严重。Liu等[8]利用数值仿真分析了内径523 mm、厚度3 mm的完全约束扰动源下球壳的响应情况。壳体模型见图1,其中扰动源半径为L,应变输出位置与旋转对称轴的距离为r′,相应的夹角为α,应变输出位置与扰动源正对位置的弧线长度为x。
图2为L=10 mm时球壳上扰动源正对位置外壁和内壁的应变曲线,可见该位置出现了明显的应变增长现象。Liu等[8]的研究表明,由于扰动源处产生的弯曲波在扰动源正对位置处汇聚,导致此位置出现了球壳上最严重的应变增长现象。
Liu等[8]根据薄壳变形机理(见图3),提出了由壳体外壁应变和内壁应变计算得到弯曲应变和膜应变的方法,其中弯曲应变为弯曲波引起壳体弯曲振动导致的应变,膜应变是壳体膜振动引起的壳体切向拉压变形导致的应变。
图4为分离出来的L=10 mm时壳体上不同位置的弯曲应变曲线,可以观察到弯曲波由扰动源产生并在球壳上的传播。其中弯曲波A为频率与球壳膜振动频率相近的第一个弯曲波,弯曲波波群的头部为波长最短的弯曲波。由于弯曲波频率越高,波速越快,波长短的波将逐渐超过波长长的波。根据弯曲波到达的壳体位置以及相应时刻,可计算出弯曲波A以及最短弯曲波的传播速度,如表1所示,其中弯曲波A的传播速度为420~450 m/s,最短弯曲波的传播速度为2 900~3 200 m/s。
表 1 弯曲波的速度Table 1. Velocities of bending waves propagating along shell.α/(°) 与膜振动频率相近的弯曲波A 最短弯曲波 到达时间/μs 平均波速/(m·s−1) 到达时间/μs 平均波速/(m·s−1) 44.2 436.3 427.1 63.4 2 964.8 88.0 909.7 427.1 131.6 3 069.1 136.0 1 436.6 443.4 203.8 3 188.7 180.0 1 894.7 267.5 Liu等[8]研究发现,弯曲波运动到扰动源正对位置后反射,而反射返回的弯曲波与入射的弯曲波由于频率相近,两者相互叠加将形成“驻波”现象,具备“驻波”的弯曲振动与膜振动叠加导致球壳上应变呈周期性空间分布。图5为L=10,62.5 mm时球壳上应变峰值的空间分布情况,可见两种情况下离扰动源正对位置较近时空间分布周期约为41 mm,离扰动源正对位置较远时空间分布周期约为82 mm。
2. 理论分析
Timoshenko在研究梁弯曲时考虑了剪切效应,推导了梁弯曲波波长与波速之间的关系[9-10]。本文中将参考相关理论,推导球壳弯曲波波长和波速的关系。取图1中球壳的微元段dx为对象,其弯曲变形如图6所示,Q表示作用在任一截面的切力,M表示相应的弯矩,α为球壳的微元段转角。需要说明的是,球壳微元除了受到外力Q、M,还受到切向拉压力以及垂直于纸面的拉压力的作用,但这两个外力对微元弯曲变形没有贡献,因此未显示在图中。
进行了以下假设和限制,将问题简化为可以进行较易求解的数学模型:
(1)弯曲波理论分析是基于平截面假定的,即假设与图6中虚线垂直的平截面在变形后仍为平面,并保持和变形后的虚线垂直。该假设可把平面弯曲问题简化为一维问题。
(2)壳体发生较小的弯曲变形,即tan α远小于1。该限制可以将后续分析的某些方程简化为线性方程。
按照动量定理和动量矩定理,可得到两个动力学公式:
∂Q∂x=ρ0A0∂V∂t (1) ∂M∂x−Q=ρ0I∂ω∂t (2) 式中:ρ0为变形前壳体的密度,A0为变形前壳体截面面积,I为截面对中性轴的转动惯量,V为壳体的横向移动速度,ω为截面转动的角速度。
在平截面假定下,可得到:
V=∂w∂t (3) tanα=∂w∂x (4) ω=∂α∂t (5) K=−∂2w∂x2[1+(∂w∂x)2]3/2 (6) 式中:w为壳体的横向位移,K为曲率。
当壳体仅发生较小的弯曲变形时,即tan α远小于1,可得到:
α=∂w∂x (7) ω=∂2w∂t∂x (8) K=−∂2w∂x2 (9) 另外,弯矩M按其定义是截面上法应力对中性轴的合力矩,为:
M=EIK (10) 式中:E为材料弹性模量。
Timoshenko认为,当弯曲波波长较短时,必须考虑剪切效应。对弯曲变形考虑剪切应力时,原来的平截面将因为剪切应变而发生挠曲,即壳体的横向位移w实际上是由于弯矩M作用下的转角α所对应的wM和由于切力Q作用下的剪切应变γ所对应的wQ两部分组成:
w=wM+wQ (11) 切力Q与切应变γ存在关系:
Q=GηA0γ (12) 式中:G为剪切弹性模量,η为壳体截面的形状系数。矩形面η的计算公式为[11]:
η=10(1+υ)12+11υ (13) 泊松比υ=0.29,则η=0.85。
把式(3)和(7)~(11)代入式(1)和式(2),则得到:
GA0η(∂2w∂x2−∂2wM∂x2)=ρA0∂2w∂t2 (14) EI∂3wM∂x3+GA0η(∂w∂x−∂wM∂x)=ρI∂3w∂t2∂x (15) 两式消去含wM的项,得到四阶偏微分方程:
EIρ0A0∂4w∂x4−IA0[1+E0ρ0ρ0Gη]∂4w∂x2∂t2+IA0ρ0Gη∂4w∂t4+∂2w∂t2=0 (16) 设
C0=√Eρ0 ,R=√IA0 ,CQ=√ηGρ0 ,则式(16)可转化为:∂4w∂x4−[1C20+1C2Q]∂4w∂x2∂t2+1C201C2Q∂4w∂t4+1C20R2∂2w∂t2=0 (17) 式(17)的通解形式为:
w=Dcos(pt−ψx) (18) 式中:D为弯曲波振幅,
ψ=2π/λ 为波数,p=ψC 为圆频率,λ 为波长,C为波速。则得到:1C201C2QC4−[1C20+1C2Q+1C20R2ψ2]C2+1=0 (19) 式(19)存在两个解,即:
C=√[C2Q+C20+C2Qλ24R2π2]±√C4Q+C40+C4Qλ416R4π4−2C2QC20+C20C2Qλ22R2π2+C4Qλ22R2π22 (20) 当λ无限小时,式(20)中右边根号内取正号“+”时,得到C=C0;右边根号内取负号“−”时,得到C=CQ。显然,C=CQ才是真解,因此式(19)的真解取:
C=√[C2Q+C20+C2Qλ24R2π2]−√C4Q+C40+C4Qλ416R4π4−2C2QC20+C20C2Qλ22R2π2+C4Qλ22R2π22 (21) 对于矩形截面,旋转半径计算公式为[12]:
R=√IA0=√3h6 (22) 取波速的无量纲数为
ˉC=CC0 ,波长的无量纲参数为ˉλ=λh ,h为球体厚度,则式(21)可变换为:ˉC=√12+12(CQC0)2+32π2ˉλ2(CQC0)2−12√1+(CQC0)4+9π4ˉλ4(CQC0)4−2(CQC0)2+6π2(CQC0)2ˉλ2+6π2(CQC0)4ˉλ2 (23) 式(23)给出了球壳上弯曲波波速和波长的关系。
3. 理论计算结果
由公式(23)得到球壳上弯曲波波速与弯曲波波长的关系曲线,见图7。
可见弯曲波的波长越短,波速越快;当波长λ无限短时,弯曲波波速趋于极限值。取弹性模量E=200 GPa,密度ρ=78 00 kg/m3,泊松比υ=0.29,其计算为:
C=CQ=√ηGρ=√ηE2ρ(υ+1) (24) 可见,最短弯曲波的波速为声速的0.574倍,约为2 895 m/s,与表1中由数值仿真得到的结果2 965~3 188 m/s相近。
另外,波的波速和波长还存在关系:
C=fλ (25) 式中:f为波的频率。
壳体膜振动主频的表达式为:
f=12πa√2Eρ(1−υ) (26) 弯曲波A为膜振动频率相近的第一个弯曲波,其频率近似取为壳体膜振动主频,取球壳半径a=261.5 mm,即可计算得到弯曲波A的频率约为5 173 Hz。
由式(23)和式(25)均可绘制出弯曲波A的波速和波长之间的关系曲线,如图8所示,两曲线的交点可以得到弯曲波A的波速。
由图8可知,频率与呼吸振动频率相近的弯曲波A的波速约为376 m/s。而表1中数值仿真得到的波速为427~443 m/s,计算结果与数值仿真结果相差(11.9~15.1)%。
Liu等[8]给出了壳体变形空间周期分布的计算公式,离扰动源正对位置较近时其周期为:
l=C2f (27) 离扰动源正对位置较远时其周期为:
l=Cf (28) 式中:C为波速,取弯曲波A的计算结果376 m/s;f为相应的频率,取计算结果5 173 Hz。则式(27)~(28)分别计算得到在离扰动源正对位置较近时周期为36.3 mm,离扰动源正对位置较远时其周期为72.7 mm,与数值仿真结果41、82 mm相差约11.5%。
可见,理论分析结果与数值仿真结果存在一定的差异,该差异主要来源于两个因素:一个因素是对理论分析的简化处理,理论分析的两个主要假设或者限制条件在数值仿真中并不完全成立;另一个因素是对数值仿真数据的读取方法,计算弯曲波波速要读取弯曲波在球壳上运动的距离和时间,壳体变形空间分布周期要读取球壳上每个周期的长度,这些数据的读取方法本身必然带来误差。
4. 结 论
参考Timoshenko梁的弯曲理论,基于平截面假定和壳体发生较小的弯曲变形的假设,推导出球壳上弯曲波波速和波长的关系,计算得到了最短弯曲波和与膜振动频率相近的弯曲波的波速,再结合早期研究提出的壳体变形分布周期与弯曲波波速的关系,计算得到了壳体变形空间分布的周期。主要结论有:
(1)理论计算结果与数值仿真结果基本吻合,其中弯曲波波速的计算结果与数值仿真结果相差在15%以内,壳体变形空间分布周期的计算结果与数值仿真结果相差在12%以内。
(2)弯曲波的波长越短,波速越快。当波长无限短时,波速趋于一极限值,约为声速的0.574倍。
本文的理论研究为解剖式分析应变增长现象的新方法提供了一定的理论依据,进一步验证了该分析方法的合理性。
-
表 1 弯曲波的速度
Table 1. Velocities of bending waves propagating along shell.
α/(°) 与膜振动频率相近的弯曲波A 最短弯曲波 到达时间/μs 平均波速/(m·s−1) 到达时间/μs 平均波速/(m·s−1) 44.2 436.3 427.1 63.4 2 964.8 88.0 909.7 427.1 131.6 3 069.1 136.0 1 436.6 443.4 203.8 3 188.7 180.0 1 894.7 267.5 -
[1] BUZUKOV A A. Characteristics of the behavior of the walls of explosion chambers under the action of pulsed loading [J]. Combustion, Explosion, and Shock Waves, 1976, 12(4): 549–554. DOI: 10.1007/BF00741150. [2] DONG Q, LI Q M, ZHENG J Y. Further study on strain growth in spherical containment vessels subjected to internal blast loading [J]. International Journal of Impact Engineering, 2010, 37(2): 196–206. DOI: 10.1016/j.ijimpeng.2009.09.001. [3] Li Q M, DONG Q, ZHENG J Y. Strain growth of the in-plane response in an elastic cylindrical shell [J]. International Journal of Impact Engineering, 2008, 35(10): 1130–1153. DOI: 10.1016/j.ijimpeng.2008.01.007. [4] 刘文祥, 张庆明, 钟方平, 等. 球壳塑性变形下的应变增长现象 [J]. 爆炸与冲击, 2017, 37(5): 893–898. DOI: 10.11883/1001-1455(2017)05-0893-06.LIU W X, ZHANG Q M, ZHONG F P, et al. Strain growth of spherical shell subjected to internal blast loading during plastic response [J]. Explosion and Shock Waves, 2017, 37(5): 893–898. DOI: 10.11883/1001-1455(2017)05-0893-06. [5] 刘文祥, 张德志, 程帅, 等. 球形爆炸容器应变增长现象的极限情况 [J]. 爆炸与冲击, 2017, 37(6): 901–906. DOI: 10.11883/1001-1455(2017)06-0901-06.LIU W X, ZHANG D Z, CHENG S, et al. Limit of strain growth in a spherical explosion vessel [J]. Explosion and Shock Waves, 2017, 37(6): 901–906. DOI: 10.11883/1001-1455(2017)06-0901-06. [6] ZHU W H, XUE H L, ZHOU G Q, et al. Dynamic response of cylindrical explosive chambers to internal blast loading produced by a concentrated charge [J]. International Journal of Impact Engineering, 1997, 19(9-10): 831–845. DOI: 10.1016/S0734-743X(97)00022-5. [7] DUFFEY T A, ROMERO C. Strain growth in spherical explosive chambers subjected to internal blast loading [J]. International Journal of Impact Engineering, 2003, 28(9): 967–983. DOI: 10.1016/S0734-743X(02)00169-0. [8] LIU W X, ZHANG Q M, ZHONG F P, et al. Further research on mechanism of strain growth caused by superposition of different vibration modes [J]. International Journal of Impact Engineering, 2017, 104: 1–12. DOI: 10.1016/j.ijimpeng.2017.01.025. [9] 王礼立. 应力波基础[M]. 2版. 北京: 国防工业出版社, 2005. [10] KOLSKY H. Stress waves in solids [M]. Oxford: Clarendon Press, 1953. [11] COWPER G R. The shear coefficient in Timoshenko’s beam theory [J]. Journal of Applied Mechanics, 1966, 33(2): 335–340. DOI: 10.1115/1.3625046. [12] 刘鸿文. 材料力学[M]. 5版. 北京: 高等教育出版社, 2011. 期刊类型引用(1)
1. 余同希,朱凌,许骏. 结构冲击动力学进展(2010-2020). 爆炸与冲击. 2021(12): 4-64 . 本站查看
其他类型引用(0)
-