Influence of water coupling coefficient on the blasting effect of red sandstone specimens

This study investigates the impact of different water coupling coefficients on the blasting effect of red sandstone.The analysis is based on the theories of detonation wave and elastic wave,focusing on the variation in wall pressure of the blasting holes.Using DDNP explosive as the explosive load,blasting tests were conducted on red sandstone specimens with four different water coupling coefficients:1.20,1.33,1.50,and 2.00.The study examines the morphologies of the rock specimens after blasting under these different water coupling coefficients.Additionally,the fractal dimensions of the surface cracks resulting from the blasting were calculated to provide a quantitative evaluation of the extent of rock damage.CT scanning and 3D reconstruction were performed on the post-blasting specimens to visually depict the extent of damage and fractures within the rock.Additionally,the volume fractal dimension and damage degree of the post-blasting specimens are calculated.The findings are then combined with numerical simulation to facilitate auxiliary analysis.The results demonstrate that an increase in the water coupling coefficient leads to a reduction in the peak pressure on the hole wall and the crushing zone,enabling more of the explosion energy to be utilized for crack propagation following the explosion.The specimens exhibited distinct failure patterns,resulting in corresponding changes in fractal dimensions.The simulated pore wall pressure–time curve validated the derived theoretical results,whereas the stress cloud map and explosion energy-time curve demonstrated the buffering effect of the water medium.As the water coupling coefficient increases,the buffering effect of the water medium becomes increasingly prominent.

Hydromechanical characterization of gas transport amidst uncertainty for underground nuclear explosion detection

Given the challenge of definitively discriminating between chemical and nuclear explosions using seismic methods alone,surface detection of signature noble gas radioisotopes is considered a positive identification of underground nuclear explosions(UNEs).However,the migration of signature radionuclide gases between the nuclear cavity and surface is not well understood because complex processes are involved,including the generation of complex fracture networks,reactivation of natural fractures and faults,and thermo-hydro-mechanical-chemical(THMC)coupling of radionuclide gas transport in the subsurface.In this study,we provide an experimental investigation of hydro-mechanical(HM)coupling among gas flow,stress states,rock deformation,and rock damage using a unique multi-physics triaxial direct shear rock testing system.The testing system also features redundant gas pressure and flow rate measurements,well suited for parameter uncertainty quantification.Using porous tuff and tight granite samples that are relevant to historic UNE tests,we measured the Biot effective stress coefficient,rock matrix gas permeability,and fracture gas permeability at a range of pore pressure and stress conditions.The Biot effective stress coefficient varies from 0.69 to 1 for the tuff,whose porosity averages 35.3%±0.7%,while this coefficient varies from 0.51 to 0.78 for the tight granite(porosity<1%,perhaps an underestimate).Matrix gas permeability is strongly correlated to effective stress for the granite,but not for the porous tuff.Our experiments reveal the following key engineering implications on transport of radionuclide gases post a UNE event:(1)The porous tuff shows apparent fracture dilation or compression upon stress changes,which does not necessarily change the gas permeability;(2)The granite fracture permeability shows strong stress sensitivity and is positively related to shear displacement;and(3)Hydromechanical coupling among stress states,rock damage,and gas flow appears to be stronger in tight granite than in porous tuff.

Three-Dimensional Sound Source Location Algorithm for Subsea Leakage Using Hydrophone

Leakages from subsea oil and gas equipment cause substantial economic losses and damage to marine ecosystem,so it is essential to locate the source of the leak.However,due to the complexity and variability of the marine environment,the signals collected by hydrophone contain a variety of noises,which makes it challenging to extract useful signals for localization.To solve this problem,a hydrophone denoising algorithm is proposed based on variational modal decomposition(VMD)with grey wolf optimization.First,the average envelope entropy is used as the fitness function of the grey wolf optimizer to find the optimal solution for the parameters K andα.Afterward,the VMD algorithm decomposes the original signal parameters to obtain the intrinsic mode functions(IMFs).Subsequently,the number of interrelationships between each IMF and the original signal was calculated,the threshold value was set,and the noise signal was removed to calculate the time difference using the valid signal obtained by reconstruction.Finally,the arrival time difference is used to locate the origin of the leak.The localization accuracy of the method in finding leaks is investigated experimentally by constructing a simulated leak test rig,and the effectiveness and feasibility of the method are verified.

Construction of Dynamic Alloy Interfaces for Uniform Li Deposition in Li-Metal Batteries

It is well accepted that a lithiophilic interface can effectively regulate Li deposition behaviors,but the influence of the lithiophilic interface is gradually diminished upon continuous Li deposition that completely isolates Li from the lithiophilic metals.Herein,we perform in-depth studies on the creation of dynamic alloy interfaces upon Li deposition,arising from the exceptionally high diffusion coefficient of Hg in the amalgam solid solution.As a comparison,other metals such as Au,Ag,and Zn have typical diffusion coefficients of 10-20 orders of magnitude lower than that of Hg in the similar solid solution phases.This difference induces compact Li deposition pattern with an amalgam substrate even with a high areal capacity of 55 mAh cm^(-2).This finding provides new insight into the rational design of Li anode substrate for the stable cycling of Li metal batteries.

Nonlinear Seebeck and Peltier effects in a Majorana nanowire coupled to leads

We theoretically study nonlinear thermoelectric transport through a topological superconductor nanowire hosting Majorana bound states(MBSs) at its two ends, a system named as Majorana nanowire(MNW). We consider that the MNW is coupled to the left and right normal metallic leads subjected to either bias voltage or temperature gradient. We focus our attention on the sign change of nonlinear Seebeck and Peltier coefficients induced by mechanisms related to the MBSs, by which the possible existence of MBSs might be proved. Our results show that for a fixed temperature difference between the two leads, the sign of the nonlinear Seebeck coefficient(thermopower) can be reversed by changing the overlap amplitude between the MBSs or the system equilibrium temperature, which are similar to the cases in linear response regime. By optimizing the MBS–MBS interaction amplitude and system equilibrium temperature, we find that the temperature difference may also induce sign change of the nonlinear thermopower. For zero temperature difference and finite bias voltage, both the sign and magnitude of nonlinear Peltier coefficient can be adjusted by changing the bias voltage or overlap amplitude between the MBSs. In the presence of both bias voltage and temperature difference, we show that the electrical current at zero Fermi level and the states induced by overlap between the MBSs keep unchanged, regardless of the amplitude of temperature difference. We also find that the direction of the heat current driven by bias voltage may be changed by weak temperature difference.

Mechanical behavior of rock under uniaxial tension:Insights from energy storage and dissipation

Many rock engineering projects show that the growth of tensile cracks is often an important cause of engineering disasters,and the mechanical behavior of rocks is essentially the transmission,storage,dissipation and release of energy.To investigate the tensile behavior of rock from the perspective of energy,uniaxial tension tests(UTTs)and uniaxial compression tests(UCTs)were carried out on three typical rocks(granite,sandstone and marble).Different unloading points were set before the peak stress to separate elastic energy and dissipated energy.The input energy density ut,elastic energy density ue,and dissipated energy density ud at each unloading point were calculated by integrating stress-strain curves.The results show that there is a strong linear relationship between the three energy parameters and the square of the unloading stress in UCT,but this linear relationship is weaker in UTT.The ue and ud increase linearly with the increase in ut in UCT and UTT.Based on the phenomenon that ue and ud increase linearly with ut,the applicability of W_(et)^(p) index in UTT was proved and the relative energy storage capacity and absolute energy distribution characteristics of three rocks in UCT and UTT were evaluated.The tensile behavior of marble and sandstone in UTT can be divided into two stages vaguely according to the energy distribution,but granite is not the case.In addition,based on dissipated energy,the damage evolution of three types of rocks in UCT and UTT was discussed.This study provides some new insights for understanding the tensile behavior of rock.

Scattering of Water Waves by Dual Symmetric Inclined Floating Porous Barriers Using the DBEM

The scattering of normally incident water waves by two surface-piercing inclined perforated barriers in water with a uniform finite depth is investigated within the framework of linear water wave theory.Considering that thin barriers are zero-thickness,a novel numerical method involving the the coupling of the dual boundary element method(DBEM)with damping layers is applied.In order to effectively damp out the reflected waves,two damping layers,instead of pseudoboundaries are implemented near the two side boundaries of the computational domain.Thus,the modified linearized free surface boundary conditions are formulated and used for solving both the ordinary boundary integral equation as well as the hypersingular boundary integral equation for degenerate boundaries.The newly developed numerical method is validated against analytical methods using the matched eigenfunction expansion method for the special case of two vertical barriers or the inclined angle to the vertical being zero.The influence of the length of the two damping layers has been discussed.Moreover,these findings are also validated against previous results for several cases.After validation,the numerical results for the reflection coefficient,transmission coefficient and dissipation coefficient are obtained by varying the inclination angle and porosity-effect parameter.The effects of both the inclination angle and the porosity on the amplitudes of wave forces acting on both the front and rear barriers are also investigated.It is found that the effect of the inclination angle mainly shifts the location of the extremal values of the reflection and the transmission coefficients.Additionally,a moderate value of the porosity-parameter is quite effective at dissipating wave energy and mitigating the wave loads on dual barriers.

Prediction of Coefficient of Restitution for Impact Elastoplastic Spheres Considering Finite Plate Thickness

Collisions between objects are a relatively common phenomenon in nature.Analyses of collision processes can greatly contribute to solving problems such as impact-rub faults and particle impacts.The coefficient of restitution is a critical parameter in the analysis of collision processes.Many experiments have shown that the coefficient of restitution is closely related to the plate thickness,and the smaller the plate thickness,the more inaccurate the coefficient of restitution predicted by the existing model,which seriously affects the process of collision analysis.To remedy this shortcoming,this paper proposes a plate thickness influence factor with the ratio of sphere diameter to plate thickness as the variable.The plate thickness influence factor can optimize the coefficient of restitution model to effectively predict the coefficient of restitution of impacting elastoplastic spheres with finite plate thickness.Finally,the validity of the new model is verified using a large amount of experimental data.

A photogrammetric approach for quantifying the evolution of rock joint void geometry under varying contact states

Accurate measurement of the evolution of rock joint void geometry is essential for comprehending the distribution characteristics of asperities responsible for shear and seepage behaviors.However,existing techniques often require specialized equipment and skilled operators,posing practical challenges.In this study,a cost-effective photogrammetric approach is proposed.Particularly,local coordinate systems are established to facilitate the alignment and precise quantification of the relative position between two halves of a rock joint.Push/pull tests are conducted on rock joints with varying roughness levels to induce different contact states.A high-precision laser scanner serves as a benchmark for evaluating the photogrammetry method.Despite certain deviations exist,the measured evolution of void geometry is generally consistent with the qualitative findings of previous studies.The photogrammetric measurements yield comparable accuracy to laser scanning,with maximum errors of 13.2%for aperture and 14.4%for void volume.Most joint matching coefficient(JMC)measurement errors are below 20%.Larger measurement errors occur primarily in highly mismatched rock joints with JMC values below 0.2,but even in cases where measurement errors exceed 80%,the maximum JMC error is only 0.0434.Thus,the proposed photogrammetric approach holds promise for widespread application in void geometry measurements in rock joints.

Influence of High-Pressure Induced Lattice Dislocations and Distortions on Thermoelectric Performance of Pristine SnTe

As a sister compound of PbTe, SnTe possesses the environmentally friendly elements. However, the pristine SnTe compounds suffer from the high carrier concentration, the large valence band offset between the L and Σpositions and high thermal conductivity. Using high-pressure and high-temperature technology, we synthesized the pristine SnTe samples at different pressures and systemically investigated their thermoelectric properties.High pressure induces rich microstructures, including the high-density dislocations and lattice distortions, which serve as the strong phonon scattering centers, thereby reducing the lattice thermal conductivity. For the electrical properties, pressure reduces the harmful high carrier concentration, due to the depression of Sn vacancies.Moreover, pressure induces the valence band convergence, reducing the energy separation between the L and Σpositions. The band convergence and suppressed carrier concentration increase the Seebeck coefficient. Thus, the power factors of pressure-sintered compounds do not deteriorate significantly under the condition of decreasing electrical conductivity. Ultimately, for a pristine SnTe compound synthesized at 5 GPa, a higher ZT value of 0.51 is achieved at 750 K, representing a 140% improvement compared to the value of 0.21 obtained using SPS. Therefore, the high-pressure and high-temperature technology is demonstrated as an effectively approach to optimize thermoelectric performance.

Effect of surfactant frequently used in soil flushing on oxygen mass transfer in micro-nano-bubble aeration system

In-site soil flushing and aeration are the typical synergetic remediation technology for contaminated sites.The surfactant present in flushing solutions is bound to affect the aeration efficiency.The purpose of this study is to evaluate the effect of surfactant frequently used in soil flushing on the oxygen mass transfer in micro-nano-bubble(MNB)aeration system.Firstly,bio-surfactants and chemical surfactants were used to investigate their effects on Sauter mean diameter of bubble(dBS),gas holdup(ε),volumetric mass-transfer coefficient(kLa)and liquid-side mass-transfer coefficient(kL)in the MNB aeration system.Then,based upon the experimental results,the Sardeing's and Frossling's models were modified to describe the effect of surfactant on kL in the MNB aeration.The results showed that,for the twenty aqueous surfactant solutions,with the increase in surfactant concentration,the value of dBS,kLa and kL decreased,while the value ofεand gas-liquid interfacial area(a)increased.These phenomena were mainly attributed to the synergistic effects of immobile bubble surface and the suppression of coalescence in the surfactant solutions.In addition,with the presence of electric charge,MNBs in anionic surfactant solutions were smaller and higher in number than in non-ionic surfactant solutions.Furthermore,the accumulation of surfactant on the gas-liquid interface was more conspicuous for small MNB,so the reduction of kL in anionic surfactant solutions was larger than that in non-ionic surfactant solutions.Besides,the modified Frossling's model predicted the effect of surfactant on kL in MNB aeration system with reasonable accuracy.

Mechanical response and dilatancy characteristics of deep marble under different stress paths:A sight from energy dissipation

Dilatancy is a fundamental volumetric growth behavior observed during loading and serves as a key index to comprehending the intricate nonlinear behavior and constitutive equation structure of rock.This study focuses on Jinping marble obtained from the Jinping Underground Laboratory in China at a depth of 2400 m.Various uniaxial and triaxial tests at different strain rates,along with constant confining pressure tests and reduced confining pressure tests under different confining pressures were conducted to analyze the mechanical response and dilatancy characteristics of the marble under four stress paths.Subsequently,a new empirical dilatancy coefficient is proposed based on the energy dissipation method.The results show that brittle failure characteristics of marble under uniaxial compression are more obvious with the strain rate increasing,and plastic failure characteristics of marble under triaxial compression are gradually strengthened.Furthermore,compared to the constant confining pressure,the volume expansion is relatively lower under unloading condition.The energy dissipation is closely linked to the process of dilatancy,with a rapid increase of dissipated energy coinciding with the beginning of dilatancy.A new empirical dilatancy coefficient is defined according to the change trend of energy dissipation rate curve,of which change trend is consistent with the actual dilatancy response in marble under different stress paths.The existing empirical and theoretical dilatancy models are analyzed,which shows that the empirical dilatancy coefficient based on the energy background is more universal.

Quantitative Analysis of Frosting Characteristics on the Drip Tray of A Cryogenic Valve for LNG Ships

Cryogenic valves play a crucial role in the production and transportation of liquefied natural gas(LNG),and are primarily responsible for efficiently controlling the inflow and outflow of LNG and regulating pressure.However,due to their operation in low-temperature and high-humidity environments,crucial components such as drip trays are susceptible to frosting,which may lead to LNG leakage,thereby causing severe safety incidents.In this study,the user-defined function(UDF)is employed to redevelop Fluent,which integrates the frost growth model with the Eulerian multiphase flow model,to conduct a quantitative analysis of frosting on drip trays of cryogenic valves.The effects of environmental parameters,such as wind speed,ambient temperature,air humidity,and cold surface temperature on the growth of the frost layer were analyzed.This study reveals a limiting wind speed between 1 m/s and 2 m/s.Upon reaching this limit speed,the growth of the frost layer reaches its maximum,and further increases in the wind speed have no significant effect on the growth of the frost layer.Furthermore,the influence of the change in the flow field on droplet impingement and freezing during the growth of the frost layer is considered through the coupling method of the kinematic characteristics of water droplets and the collection coefficient of water droplets.This study identifies the influence of different parameters on the droplet impact efficiency,leading to the modification of the frost layer on the drip tray.

Exploring the hygroscopic behavior of highly energetic oxidizer ammonium dinitramide(ADN)at different temperatures and humidities using an innovative hygroscopic modeling

Ammonium dinitramide(ADN)is a new type of green energetic oxidizer with excellent energy density and low pollution combustion characteristics.However,the strong hygroscopicity has a significant impact on its practical application.To assist in the research on moisture-proof modification of ADN materials,an innovative hygroscopic modeling approach was proposed to evaluate the hygroscopicity of ADN at various temperatures and humidities.By investigating the diffusion coefficient of water molecules in molecular dynamics processes,a visual insight into the hygroscopic process of ADN was gained.Furthermore,analyzing the non-covalent interactions between ADN and water molecules,the hygroscopicity of ADN could be evaluated qualitatively and quantitatively.The energy analysis revealed that electrostatic forces play a dominant role in the process of water adsorption by ADN,whereas van der Waals forces impede it.As a whole,the simulation results show that ADN presents the following hygroscopic law:At temperatures ranging from 273 K to 373 K and relative humidity(RH)from 10%to 100%,the hygroscopicity of ADN generally shows an increasing trend with the rise in temperature and humidity based on the results of three simulations.According to the non-hygroscopic point(298 K,52%RH)of ADN obtained by experiment in the literature,a non-hygroscopic range of temperature and humidity for ADN can be depicted when the simulation results in relative hygroscopicity is less than or equal to 17%.This study can provide effective strategies for screening anti-hygroscopic modified materials of ADN.

Kinematic analysis of flexible bipedal robotic systems

In spite of its intrinsic complexities,the passive gait of bipedal robots on a sloping ramp is a subject of interest for numerous researchers.What distinguishes the present research from similar works is the consideration of flexibility in the constituent links of this type of robotic systems.This is not a far-fetched assumption because in the transient(impact)phase,due to the impulsive forces which are applied to the system,the likelihood of exciting the vibration modes increases considerably.Moreover,the human leg bones that are involved in walking are supported by viscoelastic muscles and ligaments.Therefore,for achieving more exact results,it is essential to model the robot links with viscoelastic properties.To this end,the Gibbs-Appell formulation and Newton's kinematic impact law are used to derive the most general form of the system's dynamic equations in the swing and transient phases of motion.The most important issue in the passive walking motion of bipedal robots is the determination of the initial robot configuration with which the system could accomplish a periodic and stable gait solely under the effect of gravitational force.The extremely unstable nature of the system studied in this paper and the vibrations caused by the impulsive forces induced by the impact of robot feet with the inclined surface are some of the very serious challenges encountered for achieving the above-mentioned goal.To overcome such challenges,an innovative method that uses a combination of the linearized equations of motion in the swing phase and the algebraic motion equations in the transition phase is presented in this paper to obtain an eigenvalue problem.By solving this problem,the suitable initial conditions that are necessary for the passive gait of this bipedal robot on a sloping surface are determined.The effects of the characteristic parameters of elastic links including the modulus of elasticity and the Kelvin-Voigt coefficient on the walking stability of this type of robotic systems are also studied.The findings of this parametric study reveal that the increase in the Kelvin-Voigt coefficient enhances the stability of the robotic system,while the increase in the modulus of elasticity has an opposite effect.

Effect of non-uniform swelling on coal multiphysics during gas injection: The triangle approach

In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.

TOEPLITZ DETERMINANTS IN ONE AND HIGHER DIMENSIONS

In this study,we derive the sharp bounds of certain Toeplitz determinants whose entries are the coefficients of holomorphic functions belonging to a class defined on the unit disk U.Furthermore,these results are extended to a class of holomorphic functions on the unit ball in a complex Banach space and on the unit polydisc inℂn.The obtained results provide the bounds of Toeplitz determinants in higher dimensions for various subclasses of normalized univalent functions.

Study of deep transportation and plugging performance of deformable gel particles in porous media

Deformable gel particles(DGPs) possess the capability of deep profile control and flooding. However, the deep migration behavior and plugging mechanism along their path remain unclear. Breakage, an inevitable phenomenon during particle migration, significantly impacts the deep plugging effect. Due to the complexity of the process, few studies have been conducted on this subject. In this paper, we conducted DGP flow experiments using a physical model of a multi-point sandpack under various injection rates and particle sizes. Particle size and concentration tests were performed at each measurement point to investigate the transportation behavior of particles in the deep part of the reservoir. The residual resistance coefficient and concentration changes along the porous media were combined to analyze the plugging performance of DGPs. Furthermore, the particle breakage along their path was revealed by analyzing the changes in particle size along the way. A mathematical model of breakage and concentration changes along the path was established. The results showed that the passage after breakage is a significant migration behavior of particles in porous media. The particles were reduced to less than half of their initial size at the front of the porous media. Breakage is an essential reason for the continuous decreases in particle concentration, size, and residual resistance coefficient. However, the particles can remain in porous media after breakage and play a significant role in deep plugging. Higher injection rates or larger particle sizes resulted in faster breakage along the injection direction, higher degrees of breakage, and faster decreases in residual resistance coefficient along the path. These conditions also led to a weaker deep plugging ability. Smaller particles were more evenly retained along the path, but more particles flowed out of the porous media, resulting in a poor deep plugging effect. The particle size is a function of particle size before injection, transport distance, and different injection parameters(injection rate or the diameter ratio of DGP to throat). Likewise, the particle concentration is a function of initial concentration, transport distance, and different injection parameters. These models can be utilized to optimize particle injection parameters, thereby achieving the goal of fine-tuning oil displacement.

Enhanced piezoelectric properties and depolarization temperature in textured(Bi_(0.5)Na_(0.5))TiO_(3)-based ceramics via homoepitaxial templated grain growth

Enhanced piezoelectric response was usually achieved in(Bi_(0.5) Na_(0.5))TiO_(3)(BNT)-based ceramics with sacrifice of depolarization temperature T_(d),seriously limiting their usage range in electromechanical applications.In this work,we propose to explore piezoelectric anisotropy and domain engineering in composition&microstructure-controlled textured ceramics to resolve this issue.[001]c-textured 0.94(Bi_(0.5) Na_(0.5))TiO_(3)–0.06BaTiO_(3)(0.94BNT-0.06BT)ceramics with Lotgering factor F_(001)-91%were fabricated through homoepitaxial templated grain growth(TGG)via using 0.94BNT-0.06BT microplatelet templates.The textured samples exhibited more ordered domains with facilitated domain switching behavior,being consistent with saturated high polarization achieved at lower electric fields.Increasing F_(001) to above 60%enables rapid enhancement of piezoelectric response.Notably,compared to non-textured counterpart,the maximally textured ceramics exhibited-236%enhanced piezoelectric coefficient(d_(33)-302 pC/N)and-280%enhanced piezoelectric voltage coefficient(g_(33)-49.8×10^(−3)Vm/N),together with slightly increased depolarization temperature(T_(d)-106℃).Moreover,those values are approaching or even higher than the single-crystal values.This work not only provides important guidelines for design and synthesis of novel textured ceramics with improved comprehensive electrical properties,but also can expand application fields of BNT-based ceramics.

Global stability coefficient of large underground caverns under static loading and earthquake wave condition

Underground energy and resource development,deep underground energy storage and other projects involve the global stability of multiple interconnected cavern groups under internal and external dynamic disturbances.An evaluation method of the global stability coefficient of underground caverns based on static overload and dynamic overload was proposed.Firstly,the global failure criterion for caverns was defined based on its band connection of plastic-strain between multi-caverns.Then,overloading calculation of the boundary geostress and seismic intensity on the caverns model was carried out,and the critical unstable state of multi-caverns can be identified,if the plastic-strain band appeared between caverns during these overloading processes.Thus,the global stability coefficient for the multi-caverns under static loading and earthquake was obtained based on the corresponding overloading coefficient.Practical analysis for the Yingliangbao(YLB)hydraulic caverns indicated that this method can not only effectively obtain the global stability coefficient of caverns under static and dynamic earthquake conditions,but also identify the caverns’high-risk zone of local instability through localized plastic strain of surrounding rock.This study can provide some reference for the layout design and seismic optimization of underground cavern group.