Research on the dynamic response of shallow-buried circular non-complete bonded tunnels under anti-plane line source loading
-
摘要: 为加深理解波源距离和非完全粘结对地震波散射的影响规律,结合位移不连续模型、波函数展开法、Graf公式和镜像方法推导了反平面线源荷载下浅埋圆形非完全粘结隧道动力响应的级数解,并通过衬砌内外边界条件残余量与级数解截断项数的关系校验了该解的精度。通过对该级数解进行参数分析,系统地探讨了衬砌与围岩的接触刚度、衬砌模量、衬砌厚度、隧道埋深和波源距离等因素对衬砌内表面位移和周向剪应力的影响。结果表明:衬砌与围岩的接触刚度对隧道的动力响应具有显著的影响,尤其在某些较小接触刚度情况下隧道动力响应幅值可能非常大;增大衬砌模量会减小位移,但同时会导致周向剪应力增加;增大衬砌厚度能同时减小位移和周向剪应力;增大隧道埋深会使最大位移和周向剪应力向隧道拱顶附近移动;增大线源与隧道的水平距离会使隧道背波侧相对幅值增大。Abstract: The scattering of seismic waves by shallow-buried underground structures has significant theoretical value in the engineering field. However, previous studies have mainly focused on the case of plane waves or the case of complete bonding between lining and surrounding rock, with little consideration of the effects of source distance and non-complete bonding between lining and surrounding rock. In order to deepen the understanding of the influence of source distance and non-complete bonding on seismic wave scattering, the series solution of the dynamic response of shallow-buried circular non-complete bonded tunnels under the loading of anti-plane line source was derived based on the displacement discontinuity model, wave function expansion method, Graf formula and mirror method. The accuracy of the obtained solution was verified by the relationship between the residuals of the inner and outer boundary conditions of the lining and the number of truncated terms in the series solution. By systematically analyzing the parameters of this series solution, the influence of factors such as the contact stiffness between lining and surrounding rock, lining modulus, lining thickness, tunnel depth and source distance on the displacement and circumferential shear stress on the inner surface of the lining was discussed. The results show that the contact stiffness between lining and surrounding rock has a significant influence on the dynamic response of the tunnel, especially in cases with relatively low contact stiffness, where the amplitude of the dynamic response of the tunnel can be very large. Increasing the lining modulus reduces the displacement but increases the circumferential shear stress. Increasing the lining thickness can simultaneously reduce the displacement and circumferential shear stress. As the tunnel depth increases, the maximum displacement and circumferential shear stress on the inner surface of the lining shifts towards the apex of the tunnel. Increasing the horizontal distance between the line source and the tunnel increases the relative amplitude of the tunnel's back wave side.
-
反恐侦察需要在最短时间内了解敌方的现场情况, 抛投式侦察机器人可以直接投放到敌方控制区域, 在最短时间内为指战员提供作战现场的实时图像信息。目前抛投式侦察机器人是由士兵抛投到房屋内, 投掷距离在几米至十几米之间, 投放人员易遭受隐蔽在室内的敌人攻击。采用枪榴弹或榴弹发射器等弹射方式, 其发射过载很大(几千至数万个重力加速度), 机器人内部机械结构和传感器精密部件极易损坏, 导致整个侦察机器人失效, 所以需要对抛投式机器人做缓冲减振处理。
目前使用最多的缓冲材料是泡沫塑料, 其密度小, 弹性形变大, 缓冲性能好并且容易加工。但是对泡沫塑料的研究主要集中在其自身动态力学性能和应力应变特性[1-4], 忽略了系统的整体性。被缓冲物的物理性质、接触面积与缓冲材料的相关性质进行耦合形成了独立的系统, 该系统的动力学特性与单独的缓冲材料的动力学特性可能不同[5-7]。
本文中, 针对抛投式侦察机器人远距离部署过程中存在的高发射过载情况下结构易损坏等问题, 对机器人所使用的缓冲材料在高过载情况下的共振吸能特性进行讨论。设计抛投式机器人保护机构, 分析其结构和缓冲材料的安装方式, 建立单自由度支座激励系统数学模型。利用弹射器内弹道加速度测量系统测量抛投式机器人的缓冲动态响应, 并结合所建立的数学模型对实验数据进行理论分析。
1. 理论分析
抛投式侦察机器人、缓冲材料(EVA, 乙烯-醋酸乙烯酯共聚物)和机器人保护壳组成抛投式侦察机器人缓冲机构, 如图 1所示。抛投式机器人放置在保护壳内, 其轴向垂直填充EVA缓冲材料。该机构采用活塞式弹射器将侦察机器人发射出膛, 在保护壳底部有推杆连接, 通过推杆传递弹射器的弹射力。
根据所设计的抛投式侦察机器人的缓冲机构建立了单自由度支座激励系统模型, 如图 2所示, M为保护壳(支座)质量, m为机器人质量, k为缓冲材料的倔强系数, C为缓冲材料的阻力系数。该系统将机器人与缓冲材料组成的弹簧系统固定在机器人保护壳上, 机器人保护壳受轴向力作用进行简谐运动, 该弹射力为:
y(t)={yMsin(ft)0⩽t⩽τ0t>τ (1) 式中:yM为支座振幅, f为激励频率, τ为冲击持续时间, τ=π/f。机器人和缓冲材料在外壳所传递的力的作用下进行受迫振动[8-9]。
根据牛顿第二定律, 机器人阻尼运动微分方程为:
m¨x+C˙x+kx=0 (2) 又有:
{f2n=k/m2n=C/m (3) 式中:fn为单自由度支座系统固有频率, fn=12π√km , n为阻尼系数。解方程得缓冲机构振动周期为:
T1=2π√f2n−n2 (4) 当n≥fn时, 在该大阻尼条件下机器人缓冲机构受弹射力冲击后缓慢回复到平衡位置, 不产生振动。当n < fn时, 根据ζ=n/fn, 则有:
T1=T√1−ζ2 (5) 式中:T为缓冲机构无阻尼振动周期, ζ为阻尼比, ζ=C2√km。小阻尼条件下, 衰减振动的周期增大。
针对机器人的受迫振动, 对式(2)求对应的通解和特解, 具体求解过程见文献[6], 最终求得该缓冲系统的放大系数方程为:
β=xmyM=√1+4ζ2λ2(1−λ2)2+4ζ2λ2 (6) 式中:xm为机器人振幅, λ为频率比, λ=f/fn。
则抛投式侦察机器人在缓冲保护机构中所受最大过载为:
¨xm=¨yMβ (7) 2. 实验
2.1 实验系统
根据实验需求和弹射器内弹道特点, 研发出一套弹射器内弹道加速度测量系统, 系统中将传感器分别安装在抛投机器人和机器人保护壳上, 可以直接测量机器人和外壳所承受的发射过载。利用该系统在某靶场进行了实弹发射, 具体实验系统如图 3所示, 图 4为机器人保护壳和缓冲材料的安装方式。
2.2 实验材料
根据抛投式侦察机器人结构的强度要求, 设定机器人脆值G=100, 质量m=615g, 底面直径φ2=60mm; 保护壳内空间的高h=170mm, 直径φ=80mm; 采用全面缓冲方式包覆。通过对以上参数的解算, 得到在发射距离为50m时机器人缓冲材料所受的最大应力, 并与多种材料的缓冲系数进行匹配[6]。选取乙烯-醋酸乙烯酯共聚物为缓冲材料, 密度ρ=0.089 5g/cm3, 在相同的弹射力作用下, 增大被缓冲物与缓冲材料之间的接触面, 使材料所受的应力降低, 有助于减小缓冲材料的厚度。因此, 在保护壳结构强度和外形尺寸允许的条件下, 使被缓冲物与缓冲材料之间的接触面积最大, 缓冲材料的直径均为φ1=80mm, 对3种缓冲厚度h进行实验, 分别为50、70和90mm。每次实验所使用的机器人及保护壳质量、弹射器所使用的发射药药量和高低压室结构均相同。
2.3 实验结果及分析
为了便于对缓冲材料的缓冲效果进行量化对比, 对未填充缓冲材料时抛投式机器人所受的过载进行了测试, 如图 5所示, 整个内弹道作用时间约18ms, 加速度峰值约380g, 加速度峰值出现在约2ms时。在机器人出膛口时(约18ms时), 由于机械振动和空气阻力等原因, 加速度传感器检测的信号呈现衰减振荡, 最终变为减速运动。表 1所示, 未加缓冲材料时的固有频率项为弹射器的激励频率, 在其他厚度情况下, 所组成缓冲系统的固有频率是根据缓冲系统的倔强系数, 并利用单自由度支座系统的固有频率公式得到的。
表 1 抛投机器人缓冲系统实验数据Table 1. Experimental data of robot cushioning systemh/mm fn/Hz a/g β 0 55.50 380 50 58.12 1 200 3.150 70 49.07 480 1.260 90 43.31 410 1.070 110 30.33 314 0.826 图 6中给出了缓冲材料厚度不同的情况下机器人的加速度响应。由图 6(a)可以看出, 缓冲材料厚度为50mm时, 机器人所受加速度峰值出现在约5.2ms时, 与图 5中实验所测量的加速度峰值相比, 延迟了约3ms, 这是由缓冲材料的阻尼效应导致的。侦察机器人和缓冲材料所组成的单自由度支座激励系统的固有频率与弹射器的激励频率接近, 产生了共振, 使加速度峰值达到了1 200g; 在这种情况下, 抛投式机器人极易损坏。而在该加速度峰值后又出现一个小的加速度峰值, 是由于缓冲材料产生了振荡。
由图 6(b)可看出, 缓冲材料厚度为70mm时, 机器人所受的加速度峰值出现在约6.5ms时。由于缓冲材料的厚度增加, 其倔强系数发生变化, 缓冲系统作用时间延长, 所以机器人所受的加速度峰值出现得比图 6(a)实验中的晚。抛投式机器人与缓冲材料组成系统的结构发生变化, 缓冲材料变厚, 该系统的固有频率降低, 根据式(6)得到其放大系数为1.26, 抛投式机器人所受的加速度峰值为480g。
由图 6(c)可知, 缓冲材料厚度为90mm时, 机器人所受加速度峰值出现在约7.6ms时, 由于缓冲系统固有频率变低, 其放大系数为1.07, 机器人所受加速度峰值为410g。
根据实验数据和数学模型预测, 如在继续增加缓冲材料厚度、降低缓冲材料刚度的情况下, 单自由度支座激励系统的固有频率为30.33Hz, 可使抛投式机器人得到更好的缓冲效果, 如表 1所示, 缓冲材料厚度增加到110mm时, 机器人所受的最大过载为314g, 放大系数为0.826。
从以上分析可以看出, 系统固有频率接近弹射器激励频率时, 机器人所受加速度明显增加。根据式(6)进一步计算可知, 侦察机器人和缓冲材料组成系统的固有频率与弹射器的激励频率相差较大时, 缓冲材料起缓冲作用; 缓冲系统的固有频率与弹射器的激励频率接近时, 系统产生共振, 如图 6(a)所示。分析式(6)可知, 通过调整缓冲系统的倔强系数和阻尼系数, 并与被缓冲机构的质量对应, 使整个系统的阻尼比和频率比发生变化, 最终使放大系数小于1, 从而缓冲材料对机构进行缓冲, 避免发生共振。
3. 结论
反恐作战中常需要将抛投式机器人弹射到敌方控制区域, 但由于受机器人机械结构强度的限制, 机器人不能承受大的发射过载, 因此需要对抛投式机器人进行缓冲减振。本文中讨论了缓冲材料在高过载情况下的共振吸能特性, 根据所设计的缓冲投放系统, 建立了单自由度支座激励系统模型, 并利用内弹道加速度测量系统测量了弹射器的激励曲线和机器人的缓冲动态响应曲线。实验数据和理论分析结果表明:缓冲材料和机器人组成系统的固有频率与弹射器的激励频率接近时, 发生共振, 机器人所受过载增加, 最高过载达1 200g, 放大系数为3.42;当两者频率相差较大时, 缓冲材料才能起缓冲作用。可通过调整缓冲系统的物理性质和结构参数, 改变系统固有频率, 使放大系数小于1, 避免产生共振。
-
-
[1] LEE V W. On deformations near circular underground cavity subjected to incident plane SH waves [C]// Proceedings of the Application of Computer Methods in Engineering Conference. Los Angeles, 1977: 951–962. [2] LEE V W, TRIFUNAC M D. Response of tunnels to incident SH-waves [J]. Journal of the Engineering Mechanics Division, 1979, 105(4): 643–659. DOI: 10.1061/JMCEA3.0002511. [3] 李海波, 马行东, 李俊如, 等. 地震荷载作用下地下岩体洞室位移特征的影响因素分析 [J]. 岩土工程学报, 2006, 28(3): 358–362. DOI: 10.3321/j.issn:1000-4548.2006.03.015.LI H B, MA X D, LI J R, et al. Study on influence factors of rock cavern displacement under earthquake [J]. Chinese Journal of Geotechnical Engineering, 2006, 28(3): 358–362. DOI: 10.3321/j.issn:1000-4548.2006.03.015. [4] XIA X, LI H B, LIU Y Q, et al. A case study on the cavity effect of a water tunnel on the ground vibrations induced by excavating blasts [J]. Tunnelling and Underground Space Technology, 2018, 71: 292–297. DOI: 10.1016/j.tust.2017.08.026. [5] 袁晓铭. 地表下圆形夹塞区出平面散射对地面运动的影响 [J]. 地球物理学报, 1996, 39(3): 373–381. DOI: 10.3321/j.issn:0001-5733.1996.03.011.YUAN X M. Effect of a circular underground inclusion on surface motion under incident plane SH waves [J]. Acta Geophysica Sinica, 1996, 39(3): 373–381. DOI: 10.3321/j.issn:0001-5733.1996.03.011. [6] 刘殿魁, 林宏. 浅埋的圆柱形孔洞对SH波的散射与地震动 [J]. 爆炸与冲击, 2003, 23(1): 6–12. DOI: 10.3321/j.issn:1001-1455.2003.01.002.LIU D K, LIN H. Scattering of SH-waves by a shallow buried cylindrical cavity and the ground motion [J]. Explosion and Shock Waves, 2003, 23(1): 6–12. DOI: 10.3321/j.issn:1001-1455.2003.01.002. [7] 王国庆, 刘殿魁. SH波对浅埋相邻多个圆孔作用的动力分析 [J]. 哈尔滨工程大学学报, 2003, 24(1): 108–113. DOI: 10.3969/j.issn.1006-7043.2003.01.026.WANG G X, LIU D K. Dynamic analysis for effect of SH-wave on shallow fill multiple circular cavities [J]. Journal of Harbin Engineering University, 2003, 24(1): 108–113. DOI: 10.3969/j.issn.1006-7043.2003.01.026. [8] 陈志刚, 刘殿魁. SH波冲击下浅埋任意形孔洞的动力分析 [J]. 地震工程与工程振动, 2004, 24(4): 32–36. DOI: 10.3969/j.issn.1000-1301.2004.04.006.CHEN Z G, LIU D K. Dynamic response on a shallowly buried cavity of arbitrary shape impacted by vertical SH-wave [J]. Earthquake Engineering and Engineering Vibration, 2004, 24(4): 32–36. DOI: 10.3969/j.issn.1000-1301.2004.04.006. [9] 陈志刚. 各向异性半空间中浅埋孔洞对地表反平面运动的影响 [J]. 地震学报, 2015, 37(4): 617–628. DOI: 10.11939/jass.2015.04.008.CHEN Z G. Effect of shallow buried cavity on anti-plane motion of ground surface in anisotropic half-space [J]. Acta Seismologica Sinica, 2015, 37(4): 617–628. DOI: 10.11939/jass.2015.04.008. [10] 李敏, 刘殿魁, 周瑞芬. 含孔半圆形凸起地形及多个孔洞对SH波的散射 [J]. 哈尔滨工程大学学报, 2008, 29(1): 78–84. DOI: 10.3969/j.issn.1006-7043.2008.01.016.LI M, LIU D K, ZHOU R F. Scattering of SH-waves by a semi-cylindrical hill with a hole and multiple cavities around it in half-space [J]. Journal of Harbin Engineering University, 2008, 29(1): 78–84. DOI: 10.3969/j.issn.1006-7043.2008.01.016. [11] 刘刚, 刘殿魁. SH波入射时浅埋圆孔附近等腰三角形凸起地形的地震动 [J]. 固体力学学报, 2007, 28(1): 60–66. DOI: 10.3969/j.issn.0254-7805.2007.01.011.LIU G, LIU D K. The ground motion of an isosceles triangular hill above a subsurface cavity with incident SH-waves [J]. Acta Mechanica Solida Sinica, 2007, 28(1): 60–66. DOI: 10.3969/j.issn.0254-7805.2007.01.011. [12] 齐辉, 赵春香, 黄敏. 出平面线源荷载作用下半空间内浅埋圆孔对半圆形凸起的圆柱形弹性夹杂的动力影响 [J]. 振动与冲击, 2013, 32(17): 109–112,122. DOI: 10.3969/j.issn.1000-3835.2013.17.021.QI H, ZHAO C X, HUANG M. Dynamic effect of a subsurface cavity in half space under out-of-plane line source load on a cylindrical elastic inclusion with a semi-cylindrical hill [J]. Journal of Vibration and Shock, 2013, 32(17): 109–112,122. DOI: 10.3969/j.issn.1000-3835.2013.17.021. [13] GAO Y F, DAI D H, ZHANG N, et al. Scattering of plane and cylindrical SH waves by a horseshoe shaped cavity [J]. Journal of Earthquake and Tsunami, 2017, 11(2): 1650011. DOI: 10.1142/s1793431116500111. [14] CHEN X, ZHANG N, GAO Y F, et al. Effects of a V-shaped canyon with a circular underground structure on surface ground motions under SH wave propagation [J]. Soil Dynamics and Earthquake Engineering, 2019, 127: 105830. DOI: 10.1016/j.soildyn.2019.105830. [15] ZHANG X P, JIANG Y J, CHEN L J, et al. Anti-plane seismic performance of a shallow-buried tunnel with imperfect interface in anisotropic half-space [J]. Tunnelling and Underground Space Technology, 2021, 112: 103906. DOI: 10.1016/j.tust.2021.103906. [16] LEE V W, KARL J. Diffraction of SV waves by underground, circular, cylindrical cavities [J]. Soil Dynamics and Earthquake Engineering, 1992, 11(8): 445–456. DOI: 10.1016/0267-7261(92)90008-2. [17] 梁建文, 张浩, LEE V W. 地下洞室群对地面运动影响问题的级数解答—P波入射 [J]. 地震学报, 2004, 26(3): 269–280. DOI: 10.3321/j.issn:0253-3782.2004.03.006.LIANG J W, ZHANG H, LEE V W. A series solution for surface motion amplification due to underground group cavities—incident P waves [J]. Acta Seismologica Sinica, 2004, 26(3): 269–280. DOI: 10.3321/j.issn:0253-3782.2004.03.006. [18] LIANG J W, ZHANG H, LEE V W. A series solution for surface motion amplification due to underground twin tunnels: incident SV waves [J]. Earthquake Engineering and Engineering Vibration, 2003, 2(2): 289–298. DOI: 10.1007/s11803-003-0012-x. [19] MEI W Q, XIA Y Y, HAN G S, et al. Theoretical responses of shallow-buried circular cavity subjected to transient P wave [J]. Computers and Geotechnics, 2021, 139: 104411. DOI: 10.1016/j.compgeo.2021.104411. [20] LIN C H, LEE V W, TODOROVSKA M I, et al. Zero-stress, cylindrical wave functions around a circular underground tunnel in a flat, elastic half-space: incident P-waves [J]. Soil Dynamics and Earthquake Engineering, 2010, 30(10): 879–894. DOI: 10.1016/j.soildyn.2010.01.010. [21] LIU Q J, ZHAO M J, WANG L H. Scattering of plane P, SV or Rayleigh waves by a shallow lined tunnel in an elastic half space [J]. Soil Dynamics and Earthquake Engineering, 2013, 49: 52–63. DOI: 10.1016/j.soildyn.2013.02.007. [22] LUCO J E, DE BARROS F C P. Dynamic displacements and stresses in the vicinity of a cylindrical cavity embedded in a half-space [J]. Earthquake Engineering & Structural Dynamics, 1994, 23(3): 321–340. DOI: 10.1002/eqe.4290230307. [23] LIU Q J, YUE C, ZHAO M J. Scattering of harmonic P1 and SV waves by a shallow lined circular tunnel in a poroelastic half-plane [J]. Soil Dynamics and Earthquake Engineering, 2022, 158: 107306. DOI: 10.1016/j.soildyn.2022.107306. [24] SON M, CORDING E J. Ground–liner interaction in rock tunneling [J]. Tunnelling and Underground Space Technology, 2007, 22(1): 1–9. DOI: 10.1016/j.tust.2006.03.002. [25] ACHENBACH J D. Wave propagation in elastic solids [M]. Amsterdam: Elsevier, 1973. DOI: 10.1016/c2009-0-08707-8. -