Numerical Simulation of Explosion in Twin Tunnel System

Rapid growth of urban population in Indian cities have led traffic congestion leading to demand for scientific utilization of underground space. Immediate underground level and deep level underground below the major arterial roads are the sustainable spaces available for meeting the demand of the future traffic/transport. Due to recent increased transit activities it has become one of the soft targets by terrorists or prone to catastrophic accidents in recent years which have increased the importance of rock structures study under explosive loading. In this paper, the response of a underground metro tunnel subjected dynamic loads have been investigated including explosive capacity (30 kg TNT), ground characteristics, liner thickness and blast pressure characteristics. Blast pressure representing CONWEP air blast loading model with positive over pressure phase was applied to lining of tunnel. A three dimensional explicit finite element method was used to analyze dynamic response and damage in twin tunnels of underground metro. It is found that liner of thickness 28 cm will start deforming at the explosive loading of more than 65 kg TNT.     

A coupling method for soil structure interaction problems

This paper describes a soil‐structure coupling method to simulate blast loading in soil and structure response. For the last decade, simulation of soil behavior under blast loading and its interaction with semi buried structure in soil becomes the focus of computational engineering in civil and mechanical engineering communities. In current design practice, soil‐structure interaction analysis often assumes linear elastic properties of the soil and uses small displacement theory. However, there are numerous problems, which require a more advanced approach that account for soil‐structure interaction and appropriate constitutive models for soil. In simplified approaches, the effect of soil on structure is considered using spring‐dashpot‐mass system, and the blast loading is modeled using linearly decaying pressure–time history based on equivalent trinitrotoluene and standoff distance, using ConWep, a computer program based on semi‐empirical equations. This strategy is very efficient from a CPU time computing point of view but may not provide accurate results for the dynamic response of the structure, because of its significant limitations, mainly when soil behavior is strongly nonlinear and when the buried charge is close to the structure. In this paper, both soil and explosive are modeled using solid elements with a constitutive material law for soil, and a Jones–Wilkins–Lee equation of state for explosive. One of the problems we have encountered when solving fluid structure interaction problems is the high mesh distortion at the contact interface because of high fluid nodal displacements and velocities. Similar problems have been encountered in soil structure interaction problems. To prevent high mesh distortion for soil, a new coupling algorithm is performed at the soil structure interface for structure loading. The coupling method is commonly used for fluid structure interaction problems in automotive and aerospace industry for fuel sloshing tank, and bird impact problems, but rarely used for soil structure interaction problems, where Lagrangian contact type algorithms are still dominant. Copyright © 2012 John Wiley & Sons, Ltd.     

Blast simulation of explosives buried in saturated sand

Shallow buried explosives pose a significant threat to lightweight vehicles and their onboard personnel. To date, designers of lightweight vehicles are limited in their knowledge of what occurs during the blast. The high intensity, short term loading imparted by the explosion is enormously complex and can be significantly affected by a number of parameters including the size, shape, type, detonation point and depth of burial (DOB) of the explosive and the type, density and water content of the soil. Recent advancements in numerical simulations have enabled the complex blast event to be accurately modelled by coupling Eulerian and Lagrangian analyses: the former is well suited to modelling the blast and while the latter, the structural response. Further validation of the modelling technique is considered in the current paper, which details simulations performed utilising the coupled Eulerian-Lagrangian analysis to study the blast output of explosives buried in saturated sand. These experiments varied explosive charge size, its depth of burial, the target stand-off (SO) distance and the dimensions of the target plate. The investigation concludes with a discussion of the accuracy of the numerical simulations when compared with the experimental observations.     

Numerical simulation of round robin numerical test on tunnels using a simplified kinematic hardening model

The paper summarizes the numerical simulation of the round robin numerical test on tunnels performed in Aristotle University of Thessaloniki. The main issues of the numerical simulation are presented along with representative comparisons of the computed response with the recorded data. For the simulation, the finite element method is implemented, using ABAQUS. The analyses are performed on prototype-scale models under plane strain conditions. While the tunnel behavior is assumed to be elastic, the soil nonlinear behavior during shaking is modeled using a simplified kinematic hardening model combined with a von Mises failure criterion and an associated plastic flow rule. The model parameters are adequately calibrated using available laboratory test results for the specific fraction of sand. The soil–tunnel interface is also accounted and simulated adequately. The effect of the interface friction on the tunnel response is investigated for one test case, as this parameter seems to affect significantly the tunnel lining axial forces. Finally, the internal forces of the tunnel lining are also evaluated with available closed-form solutions, usually used in the preliminary stages of design and compared with the experimental data and the numerical predictions. The numerical analyses can generally reproduce reasonably well the recorded response. Any differences between the experimental data and the numerical results are mainly attributed to the simplification of the used model and to differences between the assumed and the actual mechanical properties of the soil and the tunnel during the test.     

Dynamic Analysis of Subway Structures Under Blast Loading

Public transit systems have become one of the targets of terrorist attacks using explosives, examples of which are the 1995 attack on Paris subway and the 2004 attack on Moscow subway. Considering the intense threats of terrorist attacks on subway systems in metropolitan areas, explicit three-dimensional Finite Element method was used to investigate the dynamic response and damage of subway structures under internal blast loading. The study was motivated by the fact that explosion in subway structure may not only cause direct life loss, but also damage the subway structure and lead to further loss of lives and properties. The study based on the New York subway system, and investigated the influences of various factors on the possible damage of subway tunnel, including weight of explosive, ground media, burial depth and characteristics of blast pressure. A mitigation measure using grouting to improve ground stiffness and strength was also analyzed. Considering the amount of explosive terrorists may use, the present study focused on small-diameter single-track tunnels, which are more vulnerable to internal blast loading and are common in New York City. Blast pressure from explosion was applied to lining surface assuming triangle pressure–time diagram, and the elasto-plasticity of ground and lining as well as their nonlinear interaction was taken into account in the numerical model. It is found from the numerical study that maximum lining stress occurred right after explosion, before the blast air pressure reduced to the atmospheric one, and it was more dependent on the maximum magnitude of air pressure than on the specific impulse, which is the area below the pressure–time curve. Small tunnels embedded in soft soil, with small burial depth, might be permanently damaged even by modest internal explosion that may be perpetuated by terrorists.     

移动环形荷载作用下饱和土中圆形衬砌隧洞动力响应研究

基于Biot理论,研究了饱和土中带有衬砌的圆形隧洞在移动环形荷载作用下的动力响应。假定衬砌为弹性体,土体为饱和多孔介质,引入两类势函数来表示土体、孔隙水和衬砌的位移,使隧洞的控制方程解耦。结合边界条件及连续条件,通过傅立叶变换得到频率-波数域中衬砌和土体的应力、位移和孔隙水压力解答,最后用傅立叶积分逆变换得到时-空域中的数值解。计算并比较了3种隧洞模型（弹性土体隧洞、饱和土体隧洞和饱和土衬砌隧洞）的动力响应分析。数值分析结果说明：（1）移动荷载速度对3种隧洞动力响应均具有较大影响;（2）弹性土体隧洞和饱和土体隧洞的动力响应具有明显区别,所以在富水地区的隧洞动力响应中土体应该视为饱和土体;（3）衬砌对隧洞动力响应有较大影响,故隧洞的动力分析中不能忽略衬砌作用。     

冲击荷载作用下软土隧道结构热-流-固耦合动力响应分析

爆炸对结构物的作用可以看成是冲击力和热的共同作用结果。基于Biot热弹性波动理论,利用饱和多孔介质热- 流-固耦合动力响应模型,研究了饱和软黏土中隧道结构内受指数衰减热/力冲击荷载作用下的动力响应。利用Flügge薄壳理论,得到隧道结构的运动方程。考虑土-结构接触面上的协调条件并利用热-流-固耦合动力响应模型的通解,得到了热/力冲击荷载作用下应力、位移和孔隙水压力响应在Laplace变换域中的解。利用Laplace数值逆变换技术得到数值计算结果,探讨了冲击荷载和热荷载作用下隧道结构刚度和厚度对应力、位移和孔隙水压力响应的影响。结果表明：在爆炸冲击荷载作用下,隧洞衬砌对周围土体起到很好的屏蔽作用,且衬砌刚度越大,其保护效果越好;衬砌与土体接触面上应力随衬砌刚度增大而减小,并且刚度越大应力衰减越快;衬砌刚度和厚度对热冲击的位移响应有明显影响,而对其温度响应的影响较小。     

简谐荷载作用下饱和土体中圆形衬砌隧道三维动力响应分析

研究了全空间饱和土体中圆形衬砌隧道在径向简谐点荷载作用下的三维动力响应,将衬砌用无限长圆柱壳来模拟,土体用Biot饱和多孔介质模型来模拟,引入两类势函数来表示土骨架的位移和孔隙水压力,并利用修正Bessel方程来求解各势函数,结合边界条件,得到频率-波数域内衬砌和土骨架位移、孔隙水压力的解答,最后进行Fourier逆变换得到时间-空间域内的响应。通过算例分析了荷载振动频率和土体渗透性对土体和衬砌位移响应及土体孔压的影响。结果表明,饱和土体和弹性土体的位移响应具有明显区别。随着荷载频率的增大,土体和隧道位移幅值减小,土体孔压幅值增大;随着土体渗透性增大,土体位移及孔压幅值减小。     

An operator‐split ALE model for large deformation analysis of geomaterials

Analysis of large deformation of geomaterials subjected to time‐varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator‐split arbitrary Lagrangian–Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid–fluid coupling and strong material non‐linearity. Each time step of the operator‐split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one‐dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical modeling of combined effects of upward and downward propagation of shear bands on stability of slopes with sensitive clay

Modeling of progressive development of zones of large inelastic shear deformation (shear band) that results from strain‐softening behavior of sensitive clays could explain the failure mechanisms of large landslides. Because of toe erosion, a shear band can be initiated and propagated upward (inward) from the river bank. On the other hand, upslope surcharge loading could generate shear bands that might propagate down towards the river bank. In the present study, upward and downward propagation of shear bands and failure of sensitive clay slopes are modeled using the Coupled Eulerian Lagrangian approach in Abaqus finite element (FE) software. It is shown that the formation and propagation of shear bands are significantly influenced by kinematic constraints that change with displacements of the soil masses, and therefore the propagation of an existing shear band might be stopped and new shear bands could be formed. The main advantages of the present FE modeling are: (i) extremely large strains in the shear bands can be successfully simulated without numerical issues, (ii) a priori definition of shearing zones is not required to tackle severe strains; instead, the FE program automatically identifies the critical locations for shear band formation and propagation. Toe erosion could significantly increase the slope failure potential because of upslope surcharge loading. FE analyses with a thick and thin sensitive clay layers show that the global failure could occur at lower surcharge loads in the former as compared to the latter cases. Copyright © 2016 John Wiley & Sons, Ltd.     

饱和半空间瞬态加载衬砌隧道的动力响应研究

工程中,地下衬砌隧道会遇到水压破裂压力、爆炸及突然开挖等瞬态荷载作用,将这些荷载理想化为作用在衬砌内边界上的均布瞬态荷载,研究圆柱形衬砌隧道在突加荷载、阶跃荷载和三角形脉冲荷载作用下的动力响应规律。根据Biot波动理论推导出半空间饱和介质的控制方程;视衬砌结构为弹性材料导出衬砌结构的控制方程。用极大半径凸圆弧近似半空间直边界,采用Graff加法公式进行坐标变换,将直角坐标表示的通解转化为极坐标表示的通解。根据边界条件导出衬砌和土体的位移、应力和孔隙压力的Laplace变换域的解答。利用反Laplace变换数值计算方法,将解答转换为时域解,得出3种瞬态荷载作用下衬砌隧道地面位移峰值、衬砌应力和孔隙压力的分布规律。     

地铁循环荷载作用下交叉隧道沉降分析

采用FLAC3D软件建立了交叉隧道的三维有限差分模型,以人工数定激励力模拟列车荷载,计算列车动荷载作用下土体的变形和应力。采用修正动偏应力长期沉降计算模型,结合分层总和法计算软土的沉降,预测交叉隧道的长期沉降规律。分析3种不同行车工况下地基长期沉降,工况1为一辆列车单独通过1#隧道,工况2为一辆列车单独通过2#隧道,工况3为两辆列车同时通过1#和2#隧道。对比不同的列车行驶速度、隧道衬砌刚度及厚度对隧道地基的长期影响。结果表明,由列车荷载引起的地基沉降主要集中在距隧道中心轴20 m范围内;在工况3下土体沉降大于工况1、2沉降之和;车速越快,沉降越小;衬砌刚度和厚度越小,沉降越大,且衬砌厚度对沉降的影响较大,衬砌刚度影响较小。     

隧道施工数值模拟及衬砌强度安全系数分析

以某隧道为背景 ,采用二维弹塑性有限元分析方法 ,模拟隧道开挖和支护过程 ,计算衬砌和围岩内力及衬砌的强度安全系数 ,研究衬砌厚度变化对其安全系数和结构稳定性的影响 ,探讨最小衬砌安全厚度 ,为类似隧道设计施工提供科学依据     

围岩弹性抗力对隧道结构受力的影响分析

运用荷载结构模型及有限元分析手段,以新建200km时速城际铁路隧道衬砌为研究对象,分析了围岩弹性抗力对二衬结构受力的影响。随着围岩弹性抗力的增大,拱顶弯矩及拱肩弯矩绝对值、拱顶竖向位移变小,而拱顶轴力及安全系数变大。说明,围岩越坚硬,衬砌结构受力越有利。随着二衬厚度的增大,拱顶弯矩及拱肩弯矩绝对值变大,而拱顶轴力、拱顶竖向位移及安全系数变小。随着混凝土强度等级的提高,拱顶弯矩、拱肩弯矩及安全系数变大,拱顶轴力和拱顶竖向位移变小,但这些量值的变化幅度较小,说明混凝土弹性模量对二衬受力影响很小。     

隧道衬砌结构三轴压缩蠕变特性试验研究

使用MTS-815液压伺服系统对隧道衬砌的混凝土试件进行了三轴压缩蠕变试验,当压缩蠕变强度达到混凝土试件抗压强度的85%时进入加速蠕变阶段。分析了隧道衬砌混凝土三轴压缩蠕变变形的特点,确定采用摩尔-库仑准则和Burgers模型串联组合而成的黏弹塑性模型可以合理地描述衬砌混凝土的蠕变特性。组合模型的黏弹性状态与Burgers模型一致,塑性状态与Mohr-Coulomb模型一致,它不仅能够反映衬砌结构蠕变的全过程,且能模拟材料的黏弹塑性变形特性,解决了典型Burgers蠕变模型只能描述黏弹性蠕变的缺陷。同时也研究了三轴试验条件下组合蠕变模型参数的求解方法,得到隧道衬砌混凝土的蠕变参数。     

Predicting Overbreak from Blast Vibration Monitoring in a Lake Tap Tunnel - A Success Story

Tunnels are required to be constructed for meeting different human needs such as power generation, transportation, underground storage, sewage etc. The predominant method of excavation, world over, is drilling and blasting owing to its capability to meet changing geo-technical conditions. Irrespective of the purpose for which the tunnels are driven, all are plagued by overbreak problems. Tunnels driven for water conveyance in hydroelectric power projects, in particular, need to be excavated with minimum overbreak to minimise cost of permanent concrete lining. Thus, predicting overbreak assumes significant importance to design site-specific blasts for minimizing rock damage. This paper presents a brief review of existing PPV (Peak Particle Velocity) based blast-induced rock damage estimation criteria and attempts to outline the ground vibration threshold levels for overbreak/rock damage in a tunnel driven through compact basalt. Rock damage manifested as overbreak is measured and correlated with the possible threshold levels of PPV. Also, the PPV levels for crack initiation and widening are proposed. The case pertains to a lake tap horizontal tunnel of Koyna Hydro-electric Power Project, India which is a water feeder tunnel for a fully underground hydroelectric power project. The tunnel was driven under a shallow rock cover of average depth ranging from 20 to 25m beneath a fully charged water body. The parting rock is mainly compact basalt. Blasting was carried out in two rounds in a controlled manner, i.e., by limiting the maximum charge per delay based on ground vibration monitoring. Ground vibration generated with free face (in second round) has been modeled and a new ground vibration propagation equation is proposed for tunnel blasting including the effect of an extra free face. The threshold limits of PPV for different degrees of overbreak/rock damage are proposed from extrapolated vibration predictor equation. The actual overbreak in the tunnel, measured using a Planimeter, varied from 2.45 per cent to 17.75 per cent of the finished tunnel area. The predicted overbreak from extrapolated PPV measurements is compared against the measured overbreak to validate the proposed blast-induced rock damage (BIRD) assessment model. The PPV level for overbreak was found to exceed 2050 mm/s in compact basalt. A linear relationship between the overbreak and maximum charge per delay is also established to design a tunnel blast in similar formations.     

东巨寺沟隧道湿喷纤维混凝土支护结构和围岩稳定性位移模糊判别

湿法喷射混凝土可有效地改善隧道内作业环境、反弹量小、喷射质量高,是隧道支护结构施工值得推广的工艺。首次在宝鸡－兰州铁路复线东巨寺沟铁路隧道选取试验段,实施湿喷工艺喷射聚丙烯纤维混凝土,作为支护结构兼作永久性支护,免除二次衬砌工作,并进行水平位移收敛监测。采用非线性最小二乘法对该隧道水平收敛曲线进行回归分析,同时采用隧道稳定性模糊概率分析理论对不同断面位置处的隧道围岩稳定性进行判别分析。分析结果表明,采用湿喷纤维混凝土支护的隧道围岩处于稳定状态,可靠度指标高。     

Post-installation pore-pressure changes around spudcan and long-term spudcan behaviour in soft clay

This paper presents a dual-stage Eulerian–Lagrangian analysis for modelling the entire process of spudcan installation in soft clay, followed by consolidation and working load operation. The analysis consists of three components, namely undrained effective stress Eulerian analysis of spudcan installation, mesh-to-mesh variable mapping and coupled-flow Lagrangian analysis for the post-installation spudcan working behaviour. The results show good agreement with centrifuge model data but also highlight the importance of replicating the hysteretic behaviour of the soil. The findings also show that while a wished-in-place approach was able to model the long-term bearing response of the spudcan, rotational stiffness was over-estimated. This is due to the fact that, while the wished-in-place analysis was able to model the hardening of the soil ahead of the spudcan, it was unable to model the softening of back-flowed soil behind spudcan. The latter influences the spudcan fixity significantly, but not bearing response. Although the analyses were conducted using ABAQUS, they can, in principle, be conducted using other codes.     

爆炸成腔后土体物理力学性能的变化规律

在6块不同的土质场地完成了21次爆炸成腔试验,试验用炸药埋深为1.20～12.5 cm,起爆药量为0.05～51 kg。对其中的5个腔体进行了土工试验。对爆后土体的天然密度、孔隙比、饱和度、压缩系数、压缩模量、内摩擦角、黏聚力的分布规律进行了实测。试验结果表明,爆炸作用对土体的影响半径是爆炸成腔半径的4～6倍。土体的天然密度、干密度、饱和度、压缩模量和黏聚力等指标随土体距爆炸中心距离的增大而降低。土体的孔隙比、压缩系数和渗透系数等指标随土体距爆炸中心距离的增大而提高。未在腔壁附近发现密度低于原状土的低密度区。爆炸冲击作用对土体内摩擦角的影响规律不明显。     

爆炸地震波作用下盾构隧道动力响应分析

盾构隧道衬砌由于各种类型接头的存在而与整体式衬砌的力学特性存在较大差异。将盾构隧道衬砌结构看作由弹塑性铰链连接的刚性管片组成,考虑围岩介质的黏弹性,提出了爆炸地震波作用下盾构隧道动力分析的简化计算方法。采用该方法对南京地铁盾构段典型横断面进行了动力分析,研究了爆炸地震波入射角度、围岩介质特性及管片厚度对结构受力与变形的影响规律。分析结果表明:波入射角度对盾构隧道有很大影响,斜入射时结构的动力响应要大于垂直入射时结构的动力响应;围岩介质等级越高,围岩对隧道结构的约束越强,隧道的抗爆性能越好;管片厚度的增大会增大结构的内力,合理设置管片厚度有利于提高盾构隧道抗爆性能。