Surface Vibration Cavity Effect of Underpass Blasting in Urban Metro Liaison Channel
-
摘要: 研究城市地下工程爆破施工地表振动衰减规律对邻近建筑保护具有重要意义。以武汉地铁8号线二期联络通道爆破开挖工程为例,运用现场监测与ANSYS/LS-DYNA三维有限元数值模拟计算相结合的方法,分析了联络通道爆破开挖作用下地表振动的空洞效应,并预测其衰减规律。研究表明:地表质点的峰值振速随着与掌子面距离的增大而不断减小,成洞区域上方地表振速明显大于未开挖区域;空洞效应放大系数随着与爆源纵向距离的增大先增大后缓慢减小,沿通道两侧随着距离的增大而不断减小,空洞效应的影响减弱;在距爆源8 m处(掌子面后方6 m),放大系数达到最大,应重点在距离爆源2~8 m即空洞效应较大的地表区域内开展振动监测。在开挖区域,与介质和爆破条件相关的系数、振动衰减系数和层介质吸收系数分别为58.52、1.43和0.019,而在未开挖区域,则分别为152.09、1.74和0.023。Abstract: The study of surface vibration attenuation law of urban underground engineering blasting construction is of great significance for the protection of adjacent buildings. This paper takes the Wuhan Metro Line 8 Phase Ⅱ liaison channel blasting excavation project as an example and uses a combination of field monitoring and ANSYS/LS-DYNA 3D finite element numerical simulation to analyze the characteristics of surface vibration hollow effect under the liaison channel blasting excavation and predict its attenuation law. The results indicate that where the surface vibration speed is significantly greater than the unexcavated area, there is a “cavity effect”; with the increase of the longitudinal distance between the mass point and the source of the explosion, the cavity effect amplification coefficient increases rapidly until the extreme value and then slowly decreases. Along both sides of the channel with the increase in distance, the amplification coefficient decreases, the effect of the cavity effect is weakened, at a distance of 8 m from the source (6 m behind the palm surface) to reach the maximum, 2−8 m from the source of the cavity effect should focus on vibration monitoring within the surface area. The excavation area with the blasting conditions is related to the coefficient of 58.52, and vibration attenuation coefficient of 1.43, while the unexcavated area of 152.09, and vibration attenuation coefficient of 1.74, the absorption coefficient of the layer media are 0.019, 0.023, respectively.
-
Key words:
- access tunnel /
- blasting excavation /
- cavity effect /
- vibration velocity /
- attenuation prediction
-
TC4钛合金是一种
α+β 型中强度钛合金,具有较高的强度和较优异的塑性,在航空、航天、船舶以及兵器领域应用广泛[1]。鸟撞问题是飞机结构在起降过程中面临的主要威胁之一。鸟撞发生时,鸟体在毫秒级时间内瞬间冲击机体结构。高速冲击产生的巨大能量将导致机身结构严重损伤,从而引发伤亡事故。因此,国际适航标准要求所有向前部件在使用前必须分析其抗鸟撞性能。大量的研究表明,鸟体在高速冲击作用下表现出明显的流体流动飞溅特性。近年来,诸多学者针对飞机结构的抗鸟撞性能开展了大量研究工作。普遍认为,高速鸟撞冲击问题是一个应变率相关的流固耦合问题。目前结构抗鸟撞性能的数值分析方法主要有3种:拉格朗日有限元法(Lagrangian finite element)、任意拉格朗日-欧拉法(Arbitrary Lagrange-Euler,ALE)和光滑粒子流体动力学方法(Smooth particle hydro-dynamic,SPH)。采用拉格朗日有限元法分析高速鸟撞冲击时,因鸟体结构变形大,致使单元发生畸变,故拉格朗日有限元法只适合模拟低速鸟撞。SPH方法是一种基于拉格朗日技术的自适应无网格粒子法,将其与有限元方法进行耦合,可在流固耦合问题求解中展现显著的优势[2-3]。由于SPH粒子在空间相互独立,因此SPH法比拉格朗日有限元法和ALE法更适于解决高速鸟撞冲击问题[4]。例如:刘军等[5]通过对比鸟撞平板叶片实验结果和SPH法及拉格朗日有限元法数值分析结果,发现SPH方法与实验结果更接近;刘富等[6]采用SPH方法进行了2024-T3铝合金平板抗高速鸟撞冲击性能研究,得到了与实验结果相近的模拟结果;Liu等[7]通过不同速度的平板鸟撞冲击实验和数值分析,研究了适用于不同鸟撞速度的鸟体模型;姚小虎等[8]通过鸟撞圆弧风挡实验和数值计算,分析了风挡玻璃在鸟撞冲击过程中的损伤破坏。
本研究采用三维图像相关法(3 dimensional digital correlate,3D-DIC),分析TC4钛合金平板高速鸟撞过程中的变形场,基于SPH方法和TC4钛合金的Johnson-Cook动态损伤模型,建立TC4钛合金平板鸟撞数值模型,并将模拟结果与鸟撞实验进行对比验证。
1. 钛合金平板鸟撞实验
鸟撞实验装置由鸟弹发射系统、TC4钛合金靶板、速度测试系统、照明系统和高速摄像系统组成。实验装置如图1所示。本实验使用的鸟弹为长L = 228 mm、直径D = 114 mm的明胶弹。鸟弹由空气炮发射,利用激光测速仪记录发射速度,激光测速仪的系统误差小于0.5%。为了解析TC4钛合金靶板背面的三维变形场,在靶板背面设置两台I-SPEED 716型高速摄影机,拍摄帧率设置为104帧每秒。位于靶板正面的两台SA-X型高速摄影机记录鸟弹飞行轨迹和撞击靶板时的响应,保证鸟弹垂直撞击TC4钛合金靶板。实验开始之前,进行调焦、视场校准和同步设置。将4台高速摄影机的触发开关通过BNC线引至操作间,其中用于动态3D-DIC测量的两台相机使用转接头连接,以实现同步触发。高速摄影机布局如图2所示。
试件材料为TC4钛合金平板,尺寸为600 mm × 600 mm × 1.6 mm。通过均匀分布的16颗M10螺栓及4.0 mm厚的夹具,将试件固定在试验工装上,夹具尺寸与螺栓分布如图3所示。
鸟撞实验共设3个发射速度,分别为149、167和180 m/s。每组进行4次重复实验。图4显示了3种速度工况下鸟撞实验结果。图4中第1行的3幅图为平板正面高速摄影图像,可以看出:鸟弹包裹在弹托中由炮管发射,在空气阻力和实验舱入射口的作用下,鸟弹和弹托在撞击TC4钛合金平板前完全分离。弹托保证了鸟弹在发射过程中的整体形状和结构不受炮管内高压气体的破坏,弹托与鸟弹的完全分离消除了弹托对TC4平板鸟撞响应的影响。图4中第2行和第3行图像分别显示了TC4钛合金平板的正面和背面鸟撞冲击结果。发射速度为149 m/s的4次实验中,TC4钛合金平板均未发生破坏;发射速度为167 m/s的4次实验中,2次发生破坏,2次未发生破坏;而发射速度为180 m/s的4次实验中,平板均发生破坏。
图5为TC4钛合金平板破坏照片。鸟体撞击平板后产生的拉伸波向外传播,在螺栓处产生剪切作用,平板发生了剪切破坏。
2. 钛合金平板鸟撞数值计算
2.1 计算模型
数值计算采用的鸟体几何模型与实验相同,为两端半球状、中间圆柱体的胶囊状柱体,长径比L/D = 2,如图6所示。鸟体模型的质量为1.8 kg。采用SPH单元模拟高速鸟撞冲击过程中的鸟体流体状飞溅,鸟体材料参数列于表1。
表 1 鸟体材料参数Table 1. Material parameters of bird bodyDensity/(kg·m−3) Elastic modulus/GPa Poisson’s ratio Yield stress/MPa Failure strain Tangent modulus/MPa 928 0.068 0.49 0.69 1.25 5 高速鸟撞实验过程中,靶板夹具和支撑架的刚度足够大,夹具和支撑架只发生线弹性变形,因此采用钢材的线弹性本构模型描述。TC4钛合金平板在高速鸟撞冲击载荷作用下发生了大变形和损伤破坏。实验发现,TC4钛合金平板的主要破坏形式是剪切破坏,因此在数值仿真计算中需要考虑剪应力的影响。大量实验表明,钛合金材料具有拉压不对称性,需要对von Mises屈服准则进行修正。本研究将Johnson-Cook动态本构模型和Johnson-Cook损伤失效模型引入邹学韬等[9]提出的von Mises修正本构框架中。该本构可以表征TC4钛合金在强冲击载荷作用下的塑性流动应力和损伤破坏行为。Johnson-Cook动态本构模型的表达式为
σs=(A+Bεn)(1+Cln˙ε∗)(1−T∗m) (1) 式中:
σs 为塑性流动应力;A为参考应变率下的屈服应力;B和n为应变强化系数;ε 为等效塑性应变;C为应变率敏感系数;˙ε∗=˙ε/˙ε0 为无量纲应变率,˙ε0 为参考应变率;T∗=(T−Tr)/(Tm−Tr) 为无量纲温度,T为温度,Tm 为材料熔化温度,Tr 为参考温度;m为温度软化系数。考虑到TC4钛合金材料的拉压不对称性,引入拉压不对称因子
G(σ) 修正von Mises屈服面,得到屈服函数为ϕ=f(σ)G(σ)=1 (2) f(σ)=3J/σ2s (3) G(σ)=exp[−c(ξ+1)]=1 (4) 式中:
f(σ)=1 为von Mises屈服面函数;ξ=cos(3θ)=272⋅J3(3J2)3/2 为Lode参数,其中θ 为Lode角,J2为偏应力第二不变量,J3为偏应力第三不变量;c 为不对称系数,可以通过两种简单应力状态求解。本研究中,不对称系数c 取c=−2ln(σs√3τs)=−2lnα,α=σs√3τs (5) 式中:
τs 为纯剪切加载的屈服应力。于是,根据相关联的流动法则,即屈服面函数等于塑性势函数g(σ) ,得到该本构的增量表达式{dσij=De(dεij−dεpij)dεpij=dλ∂g∂σijdλ=g(σtrailij)+∂g∂˙εd˙ε∂g∂σijYij−∂g∂σpeqHYij=De3σ2sexp[2lnα(ξ+1)](∂J2∂σij+2J2lnα∂ξ∂σij)H=√23∂g∂σij∂g∂σij (6) 式中:De为弹性矩阵,
dλ 为塑性流动因子,εpij 为塑性应变,σtrailij 为试探应力,σpeq 为等效塑性应力。Johnson-Cook损伤失效模型为
εf=[D1+D2exp(D3σ∗)](1+D4ln˙ε∗)(1+D5T∗) (7) 式中:
εf 为失效应变;D1~D5为材料参数;σ∗=p/σeq ,其中p为静水压力,σeq 为等效应力。数值计算所使用的本构模型参数列于表2[10-11],其中E为弹性模量,ρ为密度,μ为泊松比。
TC4钛合金平板高速鸟撞的数值计算有限元模型如图7所示。TC4靶板、夹具和M10螺栓均采用C3D8R六面体八节点减缩积分单元模拟。通过建立一般接触,计算鸟体撞击TC4钛合金靶板以及螺栓和靶板之间的接触。夹具通过16颗M10螺栓固定在支架上,在数值计算中对螺栓进行固支约束。鸟体速度分别设置为149、167和180 m/s。
在TC4钛合金平板上选取6个具有代表意义的观测点,如图8所示,其中观测点S1、S2和S3沿轴向分布,S4、S5和S6沿对角线方向分布。
2.2 计算结果
图9为鸟撞速度为149 m/s时TC4钛合金平板的等效应力云图。鸟体撞击平板后,鸟体前端受到冲击压缩后解体并呈流体状飞溅,鸟体后端仍保持固体状态。平板受鸟体冲击后产生拉伸波,并向平板四周传播。1.66 ms时鸟体完全解体,鸟体撞击的冲击能量完全耗散,此时TC4钛合金平板的应力、应变和位移达到最大值,随后开始一定程度回弹。
图10为鸟撞击后TC4钛合金平板变形的数值计算结果和3D-DIC实验结果对比。图10中149 m/s和167 m/s工况下的最大位移(Smax)图像选自未破坏实验,180 m/s工况下的最大位移图像选自平板破坏飞出前(2.00 ms前)。鸟撞过程中,TC4钛合金平板的变形大,对角线方向隆起,隆起处亮度明显增大,使得平板部分区域被遮挡,同时也遮挡了高速摄影机,因此出现部分区域未追踪到变形场的问题。由图10可知,计算得到的最大位移场与实验结果吻合较好。3种工况下数值仿真和实验得到的观测点最大位移如表3所示。
表 3 TC4钛合金平板鸟撞最大位移Table 3. Maximum displacement of titanium alloy plate impacted by a birdVelocity/(m·s−1) Method Maximum displacement/mm S1 S2 S3 S4 S5 S6 149 Sim. 65 51 36 64 55 23 Exp. 68 60 38 63 53 25 167 Sim. 75 60 48 76 65 27 Exp. 80 59 47 76 63 24 180 Sim. 119 101 63 118 103 54 Exp. 117 90 53 103 95 50 图11为计算得到的180 m/s工况下TC4钛合金平板破坏过程中的等效塑性应变云图。从图11可以看出:0.40 ms时,位于轴线上的4颗螺栓附近开始出现裂纹;0.68 ms时,平板与夹具接触处进入塑性阶段;1.04 ms时,平板对角线和夹具接触处开始起裂,并沿着夹具边缘和对角线方向扩展;2.00 ms时,最先起裂的4颗螺栓孔处裂纹贯穿。对比可见,计算得到的螺栓孔处的损伤和破坏形式与实验结果基本一致。
计算与实验得到的位移-时间曲线对比如图12所示。从图12中可以看出,计算得到的位移变化趋势及大小与实验结果基本吻合,表明本研究使用的Johnson-Cook动态本构和损伤失效模型对于模拟TC4钛合金高速鸟撞冲击问题是比较准确的。
图13对比了6个观测点的应变时程曲线。数值计算得到的6个观测点应变与实验数据的整体吻合度较高。从图13中可以看出,最靠近鸟撞点的观测点S1和S4的等效应变在0~0.2 ms内增大,0.2~1.0 ms内保持平稳,1.0~1.4 ms再次增大,1.4 ms后再次保持不变,呈现双台阶模式。其余观测点均未表现出此双台阶模式。观测点S1和S2的等效应变出现双台阶的原因在于这两个点位于鸟弹半径范围之内。鸟弹撞击TC4钛合金平板瞬间,应变瞬间增大;0.2~1.0 ms内应力波向边界传递并在边界处反向,此时S1和S2区域内材料包裹着鸟弹运动,因此应变出现平台段;1.0 ms时,边界反射的应力波再次到达S1和S2区域,使得应变再次增大。
3. 结 论
通过3D-DIC实验和数值计算方法,研究了1.8 kg鸟体高速撞击1.6 mm厚TC4钛合金平板的动态响应和损伤破坏,得到了较精确、有效的有限元模型,并得到如下结论。
(1)1.6 mm厚的TC4钛合金在1.8 kg鸟体高速撞击下的临界破坏速度为167 m/s。撞击过程中平板内部未破坏,而螺栓和夹具处发生剪切破坏。
(2)3D-DIC测试技术能够比较准确地测定鸟撞冲击过程中TC4钛合金平板的变形场。高速冲击过程中平板的变形较大,易出现光线遮挡和反光,需要设置补充高速摄影机。
(3)实验表明,鸟撞冲击后TC4钛合金平板破坏主要为螺栓等边界处的剪切破坏。将修正的von Mises屈服准则引入Johnson-Cook动态本构和损伤模型中,在本构中同时考虑拉伸和剪切两种应力状态。该本构能够准确地模拟鸟撞平板问题。
-
表 1 上台阶爆破参数
Table 1. Blasting parameters of upper bench
No. Blast hole Blast hole depth/m Quantity of blast holes Single hole charge/kg 1 Empty hole 1.0 1 0 2 Cut hole 1.2 4 0.25 3 Auxiliary hole 1.0 11 0.20 4 Peripheral hole 1.0 9 0.20 5 Bottom hole 1.0 7 0.25 表 2 数值模拟参数
Table 2. Numerical simulation parameters
Material Density/(g·m−3) E/GPa μ c/MPa φ/(°) Fill stratum 1.98 0.044 0.28 Clay 1.98 0.039 0.35 0.035 15 Weathered sandstone 2.70 50 0.25 5.5 43 Lining 2.60 32 0.20 表 3 数值模拟与现场监测数据对比
Table 3. Comparison of numerical simulation and field monitoring data
Method Monitoring vx,max/(cm·s−1) vy,max/(cm·s−1) vz, max/(cm·s−1) vr,max/(cm·s−1) Field experiment C-1 0.301 0.252 0.247 0.313 C-2 0.146 0.083 0.183 0.240 C-3 0.057 0.078 0.123 0.156 Numerical simulation C-1 0.322 0.243 0.354 0.394 C-2 0.157 0.056 0.180 0.254 C-3 0.054 0.047 0.130 0.147 -
[1] 张继春, 曹孝君, 郑爽英, 等. 浅埋隧道掘进爆破的地表震动效应试验研究 [J]. 岩石力学与工程学报, 2005, 24(22): 4158–4163. doi: 10.3321/j.issn:1000-6915.2005.22.024ZHANG J C, CAO X J, ZHENG S Y, et al. Experimental study on surface vibration effect of blasting in shallow buried tunnel boring [J]. Journal of Rock Mechanics and Engineering, 2005, 24(22): 4158–4163. doi: 10.3321/j.issn:1000-6915.2005.22.024 [2] 郭得福, 黄博, 高文乐, 等. 浅埋隧洞爆破施工引起的振动效应 [J]. 工程爆破, 2017, 23(1): 71–76. doi: 10.3969/j.issn.1006-7051.2017.01.015GUO D F, HUANG B, GAO W L, et al. Vibration effects caused by blasting construction in shallow buried tunnels [J]. Engineering Blasting, 2017, 23(1): 71–76. doi: 10.3969/j.issn.1006-7051.2017.01.015 [3] 郭辉, 张继春, 俞定, 等. 爆破地震波作用下浅埋隧道地表振动的空洞效应研究 [J]. 路基工程, 2017(4): 151–154.GUO H, ZHANG J C, YU D, et al. Study on the cavity effect of surface vibration in shallow buried tunnels under the action of blasting seismic waves [J]. Roadbed Engineering, 2017(4): 151–154. [4] 刘光汉, 周建敏, 余红兵. 浅埋隧道掘进爆破空洞效应研究 [J]. 采矿技术, 2017, 17(5): 112–113. doi: 10.3969/j.issn.1671-2900.2017.05.037LIU G H, ZHOU J M, YU H B. Research on blast cavity effect in shallow buried tunnel boring [J]. Mining Technology, 2017, 17(5): 112–113. doi: 10.3969/j.issn.1671-2900.2017.05.037 [5] 李志堂, 尹荣申, 孟亚锋, 等. 土岩交错地层隧道爆破施工的振动响应及空洞效应分析 [J]. 隧道建设, 2018, 38(4): 588–593.LI Z T, YIN R S, MENG Y F, et al. Analysis of vibration response and cavitation effect of tunnel blasting construction in soil-rock interlaced strata [J]. Tunnel Construction, 2018, 38(4): 588–593. [6] 石连松, 高文学, 王林台. 地铁浅埋隧道爆破振动效应试验与数值模拟研究 [J]. 北京理工大学学报, 2018, 38(12): 1237–1243.SHI L S, GAO W X, WANG L T. Experimental and numerical simulation study of blasting vibration effects in shallow buried underpass tunnels [J]. Journal of Beijing University of Technology, 2018, 38(12): 1237–1243. [7] 刘志波. 莲花山隧道爆破振动效应试验研究 [J]. 爆破, 2020, 37(3): 78–84. doi: 10.3963/j.issn.1001-487X.2020.03.013LIU Z B. Experimental study of blasting vibration effect in Lotus Hill tunnel [J]. Blasting, 2020, 37(3): 78–84. doi: 10.3963/j.issn.1001-487X.2020.03.013 [8] 冯小冬. 地铁钻爆法施工对邻近建筑物的振动响应预测 [J]. 地下空间与工程学报, 2021, 17(2): 580–589.FENG X D. Prediction of vibration response of subway construction by drill-and-blast method on adjacent buildings [J]. Journal of Underground Space and Engineering, 2021, 17(2): 580–589. [9] 中国工程爆破协会. 爆破安全规程: GB 6722—2014 [S]. 北京: 中华人民共和国国家质量监督检验检疫总局, 2014. [10] ZHU B, JIANG N, ZHOU C, et al. Dynamic failure behavior of buried cast iron gas pipeline with local external corrosion subjected to blasting vibration [J]. Journal of Natural Gas Science and Engineering, 2021. [11] 中国建筑科学研究院. 混凝土结构设计规范: GB 50010—2015 [S]. 北京: 中华人民共和国住房和城乡建设部, 2015. [12] 孔丹丹, 赵颖华, 王萍, 等. 钢筋混凝土材料有限元分析中的等效模量方法 [J]. 沈阳建筑大学学报(自然科学版), 2005, 21(3): 200–203.KONG D D, ZHAO Y H, WANG P, et al. Equivalent modulus method in finite element analysis of reinforced concrete materials [J]. Journal of Shenyang University of Architecture (Natural Science Edition), 2005, 21(3): 200–203. [13] 吕国鹏, 周传波. 隧道断层带注浆加固围岩体爆破动力损伤特征 [J]. 岩石力学与工程学报, 2021, 40(10): 2038−2047.LÜ G P, ZHOU C B. Blasting dynamic damage characteristics of grouting reinforced surrounding rock mass in tunnel fault zone [J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(10): 2038−2047. [14] 罗杰峰. 隧道掘进爆破合理微差时间研究 [D]. 长沙: 长沙理工大学, 2014.LUO J F. Research of reasonable millisecond delay time with tunneling blasting [D]. Changsha: Changsha University of Science and Technology, 2014. [15] 蔡军, 苏莹, 邱秀丽. 爆破荷载作用下空洞效应对围岩振动速度的影响[J]. 矿冶工程, 2021, 41(5): 10−13, 17.CAI J, SU Y, QIU X L. Influence of cavity effect on vibration velocity of surrounding rock under blasting load [J]. Mining and Metallurgy Engineering, 2021, 41(5): 10−13, 17. -