Numerical Simulation of Surrounding Rock Deformation and Grouting Reinforcement of Cross-Fault Tunnel under Different Excavation Methods

Tunnel construction is susceptible to accidents such as loosening, deformation, collapse, and water inrush, especiallyunder complex geological conditions like dense fault areas. These accidents can cause instability and damageto the tunnel. As a result, it is essential to conduct research on tunnel construction and grouting reinforcementtechnology in fault fracture zones to address these issues and ensure the safety of tunnel excavation projects. Thisstudy utilized the Xianglushan cross-fault tunnel to conduct a comprehensive analysis on the construction, support,and reinforcement of a tunnel crossing a fault fracture zone using the three-dimensional finite element numericalmethod. The study yielded the following research conclusions: The excavation conditions of the cross-fault tunnelarray were analyzed to determine the optimal construction method for excavation while controlling deformationand stress in the surrounding rock. The middle partition method (CD method) was found to be the most suitable.Additionally, the effects of advanced reinforcement grouting on the cross-fault fracture zone tunnel were studied,and the optimal combination of grouting reinforcement range (140°) and grouting thickness (1m) was determined.The stress and deformation data obtained fromon-site monitoring of the surrounding rock was slightly lower thanthe numerical simulation results. However, the change trend of both sets of data was found to be consistent. Theseresearch findings provide technical analysis and data support for the construction and design of cross-fault tunnels.

The Ductile Deformation Characteristics of Caledonian Intracontinental Orogeny in the Northeastern Jiangshan-Shaoxing Tectonic Zone:Insights from Magnetic Fabric Study and Its Geodynamic Implication

The Jiangshan-Shaoxing tectonic zone was the northeastern boundary between the Yangtze Block and the Cathaysia Block during the Neoproterozoic and was an intracontinental orogenic belt during late of the early Paleozoic. In this tectonic zone, there develops a lot of mylonite underwent strong ductile deformation and schist, gneiss, and amphibolite with medium and high grade metamorphism which was formed during the late of early Paleozoic. The research of geometry and kinematic of ductile deformation in Jiangshan-Shaoxing tectonic zone is very important to reveal the tectonic process of intracontinental orogeny. This paper uses the anisotropy of magnetic susceptibility （AMS） to determine the ductile deformation geometry and kinematic of Jiangshan-Shaoxing tectonic zone combing with the field survey. In this study, 190 specimens of 19 locations and 221 specimens of 23 locations from Wangjiazhai section and Lipu-Sizhai section were analyzed. The magnetic foliation over magnetic lineation in both Wangjiazhai and Lipu-Sizhai sections together with the field observations indicated a compressional deformation pattern. 3 and 4 strong ductile deformation zones can be established in the Wangjiazhai section and the Lipu-Sizhai section, respectively. According to the magnetic fabric and petro-fabric studies, the Northeastern Jiangshan-Shaoxing tectonic zone suffered two kinds of deformation patterns during the late early Paleozoic, i.e., the thrusting deformation followed by sinistral shear deformation.

On the calibration and verification of Voronoi-based discontinuous deformation analysis for modeling rock fracture

Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.

Non-contact monitoring and analysis system for tunnel surrounding rock deformation of underground engineering

It is very important to monitor surrounding rock deformation in tunnel construction. The principle, function, development and application of the system composed of a total station and computer for monitoring and analyzing surrounding rock deformation were discussed. The new methods of two free station of 3D measurement and its mathematic adjustment mode were presented. The development of software for total station on-board and post for computer were also described. Without centering it and measuring its height, the total station controlled by the software on-board can fulfill the whole measurements to target points. Monitoring data can be processed by the post software and results of regression analysis, forecasting information of the tunnel surrounding rock deformation can be provided in time. The practical use shows that this system is practicable, highly accurate and efficient. It satisfies the needs of safety and information construction in tunnel construction of underground engineering.

Performance analysis of empirical models for predicting rock mass deformation modulus using regression and Bayesian methods

Deformation modulus of rock mass is one of the input parameters to most rock engineering designs and constructions.The field tests for determination of deformation modulus are cumbersome,expensive and time-consuming.This has prompted the development of various regression equations to estimate deformation modulus from results of rock mass classifications,with rock mass rating(RMR)being one of the frequently used classifications.The regression equations are of different types ranging from linear to nonlinear functions like power and exponential.Bayesian method has recently been developed to incorporate regression equations into a Bayesian framework to provide better estimates of geotechnical properties.The question of whether Bayesian method improves the estimation of geotechnical properties in all circumstances remains open.Therefore,a comparative study was conducted to assess the performances of regression and Bayesian methods when they are used to characterize deformation modulus from the same set of RMR data obtained from two project sites.The study also investigated the performance of different types of regression equations in estimation of the deformation modulus.Statistics,probability distributions and prediction indicators were used to assess the performances of regression and Bayesian methods and different types of regression equations.It was found that power and exponential types of regression equations provide a better estimate than linear regression equations.In addition,it was discovered that the ability of the Bayesian method to provide better estimates of deformation modulus than regression method depends on the quality and quantity of input data as well as the type of the regression equation.

Geological analysis of gravitational rock slope deformation：a case from Nujiang River, China

This paper presents a study on the gravityinduced rock slope deformation observed along the Nujiang River in China. We performed a comprehensive field investigation and analysis to identify the deformation pattern of the slope and its triggering factors. Moreover, a geologicalevolutionary model was developed, and it considers the effects of river incision and rock mass degradation caused by weathering and simulates the mechanisms underlying the initiation and progression of the slope deformation. The results support the proposed failure mechanism in which fractures within the slope are induced by rock mass degradation caused by weathering. Importantly, the modeling reveals that compressional deformation at the toe of the slope results in a tensile failure in the upper portion of the slope, demonstrating that the rock mass in the slope toe is the key factor inducing slope deformation. This analysis of slope deformation and its spatial and temporal correlations with rock weathering and river incision reveal the main triggering factors that control the evolution of the studied slope and provide insights into the deformation process.

FORMULATIONS OF THE THREE-DIMENSIONAL DISCONTINUOUS DEFORMATION ANALYSIS METHOD

This paper extends the original 2D discontinuous deformation analysis(DDA)method proposed by Shi to 3D cases,and presents the formulations of the 3D DDA.The formulations maintain the characteristics of the original 2D DDA approach.Contacts between the blocks are detected by using Common-Plane (C-P) approach and the non-smooth contact,such as of vertex-to-vertex,vertex- to-edge and edge-to-edge types,can be handled easily based on the C-P method.The matrices of equilibrium equations have been given in detail for programming purposes.The C program codes for the 3D DDA are developed.The ability and accuracy of the formulations and the program are verified by the analytical solutions of several dynamic examples.The robustness and versatility of the algorithms presented in this paper are demonstrated with the aid of an example of scattering of densely packed cubes.Finally,implications and future extensions are discussed.

Anisotropic NdFeB Magnet Fabricated by Single Stroke Hot Deformation

The general hot deformation process consists of two steps, hot pressing and die-upsetting in order to obtain the anisotropic NdFeB magnet. This is the first report that the high anisotropy NdFeB magnets can be fabricated by single stroke hot deforming the isotropic magnet. The magnetic properties of those materials are: coercivity iHc ～11 kOe, remanence Br ～12 kG, and the maximum energy product (BH)max ～28 MG.Oe.

Numerical simulation on deformation character of surrounding rock masses of Changjiashan tunnel through the gob of coalmine

Based on the construction project of the Changjiashan tunnel of the freeway,the variety rule of surrounding rock masses of the tunnel through the gob of coalmine wasstudied by using of finite element methed(FEM).The status of the stress and strain,thevariety of the plastic area were simulated in the whole rock mass before and after thetunnel was excavated.The characters of stress and deformation of surrounding rockmasses were analyzed when the tunnel was built.It concluded from the numerical simula-tion that the influence on the tunneling is great when the tunnel passing through the gob ofcoalmine is excavated,and the relative measures should be taken.

Characteristics of Metamorphic Rock Magnetic Fabrics in the Nyalam Area of the Southern Tibet and Its Geological Significance, China

This paper reports Precambrian rock magnetic fabrics in the Nyalam area of southern Tibet. The analytical results of magnetic fabrics show that the values of H are high （〉 10% in general）, so the ductile deformations of the Precambrian rock are strong. The orientation of the maximum principal stress inferred from the minimum magnetic susceptibility is nearly S-N, NE-SW and NW-SE. The Flinn diagram of the magnetic fabrics show that the strain pattern is oblate and constrictional type. Magnetic foliation of great majority of rock samples is well developed and the magnetic lineation is poor and the magnetic susceptibility ellipsoid is flattened. The magnetic lineation of the minority rock samples is well developed and the magnetic foliation is poor and the magnetic susceptibility ellipsoid is prolate. According to the geological field and the magnetic fabrics, there are 3 times tectonic stress field in SN directed extruding, NW-SE directed extruding, NW-SE directed extension. It shows that the Nyalam area has undergone process the orientation of SN, NW-SE nappe structure and NW-SE directed extension structure. The change of tectonic stress is reflected by the field characteristics of the Precambrian rock magnetic fabrics that is the direct responding result of the arc-continental, continent-continental collision between the India and Asian continents in the late part of the Late Cretaceous to Late Eocene and subsequently shifted to intra-continental convergent, the plateau uplifting and extension structure stage since the Late Eocene.

Effect of Seismic Permanent Deformation on Safety and Stability of Earth-Rock Dam Slope

In order to study the effect of seismic permanent deformation on the safety and stability of earth-rock dam, the permanent deformation is considered as the non-design permanent load, and the stress-strain hysteresis curve is also considered when the earth is under cyclic load. The research work can make the calculation results of plastic collapse more accurate by including the effect of the post-earthquake degree of plastic deformation on the stability of the earth-rock dam, and the dam safety factor decreases from 2.50 to 1.90 after the magnitude-8 earthquake. Moreover, the research work will also improve the design of the earth-rock dam under abnormal operating conditions.

Late Mesozoic Tectonic Evolution of Southwestern Fujian Province,South China:Constraints from Magnetic Fabric,Zircon U-Pb Geochronology and Structural Deformation

A combined study of magnetic fabrics, zircon U-Pb geochronology and structural deformation was carried out for Late Paleozoic sedimentary and Mesozoic magmatic rocks in the southwestern Fujian rift basin, South China, aiming at deciphering the tectonic evolution during Late Mesozoic. Field observations showed that the Late Mesozoic structure deformations in southwestern Fujian were categorized into four phases： NW-SE compression, ENE-WSW extension, NNE-SSW compression and NNW-SSE extension, se- quentially. Zircons picked out from Juzhou granite and WNW-trending diabase dykes showed complete crys- tal shapes and clear oscillatory zonings on their edges, and the U-Pb dating yielded ages of 132 and 141 Ma, respectively. The susceptibility ellipsoid magnitude parameters of the Juzhou granite are characterized by flaser type strain ellipsoid, with pole density center of K3 falling into the first and the third quadrants, these fea- tures revealed that the Juzhou granite formed in ENE-WSW compressional stress field, indicating the early stage of Early Cretaceous extrusion in southwestern Fujian. The late stage of Early Cretaceous NNE-SSW ex- tension was limited by the widespread WNW-trending diabase dykes, which were usually regarded as impor- tant indications for a regional extensional setting. On the basic of the previous researches, structural deforma- tion studies, and the deductions above, it can be concluded that southwestern Fujian experienced five main tectonic stages during Late Mesozoic： Early Jurassic extension, Middie-Late Jurassic thrusting, early stage of Early Cretaceous extension, late stage of Early Cretaceous compression and Late Cretaceous extension.

Static and dynamic stability analysis using 3D-DDA with incision body scheme

Discontinuous deformation analysis （DDA） provides a powerful numerical tool for the analysis of discontinuous media. This method has been widely applied to the 2D analysis of discontinuous deformation. However, it is hindered from analyzing 3D rock engineering problems mainly due to the lack of reliable 3D contact detection algorithms for polyhedra. Contact detection is a key in 3-D DDA analysis. The limitations and advantages of existing contact detection schemes are discussed in this paper, and a new approach, called the incision body （IB）, is proposed, taking into account the advantages of the existing methods. A computer code 3DIB, which uses the IB scheme as a 3D contact detection algorithm, was programmed with Visual C^＋＋. Static and dynamic stability analysis for three realistic engineering problems has been carried out. Furthermore, the focus is on studying the stability of a gravity dam on jointed rock foundation and dynamic stability of a fractured gravity dam subject to earthquake shaking. The simulation results show that the program 3DIB and incision body scheme are capable of detecting 3D block contacts correctly and hence simulating the open-close and slide process of jointed block masses. In addition, the code 3DIB could provide an effective tool for evaluating the safety of 3D dam structures, which is quite important for engineering problems.

Simulation of failure process of jointed rock

A modified discontinuous deformation analysis (DDA) algorithm was proposed to simulate the failure behavior of jointed rock. In the proposed algorithm, by using the Monte-Carlo technique, random joint network was generated in the domain of interest. Based on the joint network, the triangular DDA block system was automatically generated by adopting the advanced front method. In the process of generating blocks, numerous artificial joints came into being, and once the stress states at some artificial joints satisfy the failure criterion given beforehand, artificial joints will turn into real joints. In this way, the whole fragmentation process of rock mass can be replicated. The algorithm logic was described in detail, and several numerical examples were carried out to obtain some insight into the failure behavior of rock mass containing random joints. From the numerical results, it can be found that the crack initiates from the crack tip, the growth direction of the crack depends upon the loading and constraint conditions, and the proposed method can reproduce some complicated phenomena in the whole process of rock failure.

Simulation of blast induced crater in jointed rock mass by discontinuous deformation analysis method

Rock blasting is a dynamic process accom panied with the propagations of shock waves and the dispersion of the explosion gas.This paper adopts the discontinuous deformation analysis(DDA)method to simulate the rock blasting process.A dynamic parameter adjustment and the non-reflecting boundary condition are implemented in the DDA method.The sub-block DDA method to simulate fracture problems is used.The blasting process in jointed rock mass is simulated by application of the explosion gas pressure on the expanding borehole walls and induced connected fracture surfaces around the boreholes.The blast craters with different overburdens are derived.The whole process including the explosion gas dispersion,borehole expansion,rock mass failure and cast,and the formation of the final blasting piles in rock blasting are well reproduced numerically.Parametric study for different overburdens is carried out,and the results are analyzed and discussed.

A rapid micro-magnetic resonance imaging scanning for three-dimensional reconstruction of peripheral nerve fascicles

The most common methods for three-dimensional reconstruction of peripheral nerve fascicles include histological and radiology techniques. Histological techniques have many drawbacks including an enormous manual workload and poor image registration. Micro-magnetic resonance imaging（Micro-MRI）, an emerging radiology technique, has been used to report results in the brain, liver and tumor tissues. However, micro-MRI usage for obtaining intraneural structures has not been reported. The aim of this study was to present a new imaging method for three-dimensional reconstruction of peripheral nerve fascicles by ~1T micro-MRI. Freshly harvested sciatic nerve samples from an amputated limb were divided into four groups. Two different scanning conditions（Mannerist Solution/GD-DTPA contrast agent, distilled water） were selected, and both T1 and T2 phases programmed for each scanning condition. Three clinical surgeons evaluated the quality of the images via a standardized scale. Moreover, to analyze deformation of the two-dimensional image, the nerve diameter and total area of the micro-MRI images were compared after hematoxylin-eosin staining. The results show that rapid micro-MRI imaging method can be used for three-dimensional reconstruction of the fascicle structure. Nerve sample immersed in contrast agent（Mannerist Solution/GD-DTPA） and scanned in the T1 phase was the best. Moreover, the nerve sample was scanned freshly and can be recycled for other procedures. MRI images show better stability and smaller deformation compared with histological images. In conclusion, micro-MRI provides a feasible and rapid method for three-dimensional reconstruction of peripheral nerve fascicles, which can clearly show the internal structure of the peripheral nerve.

基于Bayes-LSTM的公路隧道围岩变形预测方法研究

在公路隧道施工过程中,围岩的稳定性对隧道施工的影响较大。因此公路隧道围岩变形的监控量测与准确预测是保障隧道施工安全的关键。针对当前隧道围岩变形的预测精度较低以及泛化能力较差等问题,该文提出一种基于贝叶斯(Bayes)优化长短期记忆网络(LSTM)的方法,该方法首先对拱顶沉降和周边收敛的原始监测数据进行预处理,而后构建公路隧道拱顶沉降与周边收敛的初始LSTM模型,并利用Bayes优化模型中的超参数,最终得出预测结果。利用该模型对某公路隧道拱顶沉降和周边收敛进行预测,将预测结果以均方根误差为评价指标与神经网络(CNN)和支持向量回归(SVR)进行对比。预测拱顶沉降时,Bayes-LSTM模型的平均预测精度相较于CNN与SVR模型分别提高了1.0与1.26;预测周边收敛时,Bayes-LSTM模型平均精度相较于CNN与SVR分别提高了0.3与0.32。表明Bayes-LSTM模型的预测精度较高,同时其能在训练模型过程中对历史信息进行判断和取舍,极大地提高了时序数据处理的效率,为公路隧道围岩变形预测提供了新的思路和探索。

节理化炭质页岩地层隧道围岩大变形及控制技术研究

某高速铁路XHS隧道穿越节理化炭质页岩地层,在施工过程中围岩大变形、失稳坍塌现象显著,现场采用强支护和仰拱加深等措施后围岩变形控制效果不佳。针对XHS隧道节理化炭质页岩地层地质条件,结合现场监测手段、离散-连续耦合数值模拟分析围岩大变形及破坏特征,基于数值模拟提出以采取地层预加固为主的围岩变形控制措施,并通过现场试验探讨该控制措施的应用效果。研究结果表明:隧道开挖后,围岩变形具有变形量大、变形速率快的特点,围岩拱部沉降量大于水平收敛量且变形具有非对称的分布特征;受三台阶法多次开挖扰动影响,围岩卸荷范围动态发育并不断向全环扩展,松动区逐渐由浅部围岩向深部转移,并呈现出非对称的破坏特征,最终引发围岩大变形;采取地层预加固后,模型中围岩变形和松动区范围明显减小,围岩非对称变形破坏也得到了有效控制;在采用地层预加固、管棚超前支护、三台阶临时仰拱法开挖的控制措施后,围岩变形得到控制,施工效果良好,隧道恢复正常施工,保证了隧道的顺利贯通;以地层预加固为主的控制措施是此类节理化炭质页岩地层围岩变形控制的有效手段。

嵌岩桩水平承载力地基反力Matlock法分析

嵌岩桩在海上风电、跨海大桥等结构中的应用十分广泛,但在嵌岩桩水平承载特性方面的研究仍有较多不足。弹性地基反力法中的Matlock法在中国是用于计算水平承载桩承载力时应用最多的方法,但地基反力Matlock法针对嵌岩桩缺少合理的取值方法。基于现场试验数据建立并验证有限差分数值模型,通过数值模拟与参数分析,研究了岩体力学参数和桩基参数等因素对桩岩相互作用的影响,提出了Matlock法中地基反力系数比例系数的拟合计算公式。通过现场试验数据对提出的方法进行了验证,结果表明计算结果与实际情况较为相符,与其他方法比较,所提出的修正Matlock法相对非线性的p-y曲线法具有取值计算较为简便的优势。

深部岩体变形模量数值反演计算方法及敏感性分析

岩体变形模量是地下岩石工程设计和建造的关键参数之一。获取准确的岩体变形模量对深入分析深部岩体工程围岩变形行为至关重要。基于此,系统梳理归纳了常用的岩体变形模量计算方法,并利用有限元方法模拟原位试验进行岩体变形模量反演计算,与原位试验结果进行对比验证分析。针对深竖井实际工程,建立了深部地层不同岩性岩体变形模量反演分析数值计算模型,开展了不同强度准则、泊松比和岩体黏聚力对岩体变形模量的敏感性分析,将数值模拟反演结果与经验公式计算结果进行对比分析和讨论。研究结果表明,基于刚性平板载荷试验计算得到的岩体变形模量值与有限元数值反演得到的数值高度接近,具有较好的一致性;虽然泊松比对岩体变形模量有一定影响,但总体上岩体变形模量对泊松比的敏感性较低;当岩体黏聚力较小时,强度准则的选择对变形模量的反演结果影响显著,随着岩体黏聚力的增大,强度准则的选择对变形模量反演结果的影响逐渐减小;对于不同深度岩性的地层,岩体变形模量对岩体黏聚力变化敏感性不同;通过数值反演求得的岩体变形模量和一些经验公式的计算结果具有很好的相关性。因此,深部岩体工程规划设计和开发过程中可综合运用数值模拟和经验公式进行岩体变形模量反演计算,进而在更大程度上保障工程的安全性和经济性。