Compensation excavation method control for large deformation disaster of mountain soft rock tunnel

In recent years,the mine tunneling method and the new Austrian tunneling method have been considered the main theories of tunneling approaches in China.It is difficult for the traditional technique to overcome the large deformation problems imposed by complex geological conditions of mountain soft rock tunneling.Hence,the compensation excavation method has been proposed to solve this issue under the consideration that all damage in tunneling originates from the excavation.It uses supportive strategies to counteract the excavation effects successfully.This paper provides an overview of the fundamental ideas of the compensation excavation method,methodologies,and field applications.The scientific validity and feasibility of the compensation excavation method were investigated through the practical engineering study of the Muzhailing and Changning tunnels.

Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Large deformation analysis of a cantilever beam made of axially functionally graded material by homotopy analysis method

Large deformation of a cantilever axially functionally graded (AFG) beam subject to a tip load is analytically studied using the homotopy analysis method (HAM). It is assumed that its Young’s modulus varies along the longitudinal direction according to a power law. Taking the solution of the corresponding homogeneous beam as the initial guess and obtaining a convergence region by adjusting an auxiliary parameter, the analytical expressions for large deformation of the AFG beam are provided. Results obtained by the HAM are compared with those obtained by the finite element method and those in the previous works to verify its validity. Good agreement is observed. A detailed parametric study is carried out. The results show that the axial material variation can greatly change the deformed configuration, which provides an approach to control and manage the deformation of beams. By tailoring the axial material distribution, a desired deformed configuration can be obtained for a specific load. The analytical solution presented herein can be a helpful tool for this procedure.

Frictional contact algorithm study on the numerical simulation of large deformations in deep soft rock tunnels

There exist three types of nonlinear problems in large deformation processes of deep softrock engineering, i.e., nonlin- earity caused by material, geometrical and contact boundary. In this paper, the numerical method to tackle the nonlinear eontact and large deformation problem in A Software on Large Deformation Analysis for Soft Rock Engineering at Great Depth was presented. In the software, based on Lagrange multiplier method and Coulomb friction law, kinematic constraints on contact boundaries were introduced in functional function, and the finite element equations was established for two incremental large deformation analyses models, polar decomposition model and additive decomposition model. For every incremental loading step, by searching for the contact points in the potential contact interfaces （the excavation boundaries）, the Lagrange multipliers, i.e., contact forces are cal- culated iteratively by Gauss-Seidel method, and justified through satisfy the inequalities of static constraint on contact boundaries. With the software, large deformation and frictional contact of a transport roadway were analyzed numerically by the two models. The numerical examples demonstrated the efficiency of the method used in the software.

Large deformation simulations of geomaterials using moving particle semi-implicit method

Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.

Dynamic analysis of a rotating rigid-flexible coupled smart structure with large deformations

Based on Hamilton＇s principle, a new kind of fully coupled nonlinear dynamic model for a rotating rigid-flexible smart structure with a tip mass is proposed. The geometrically nonlinear effects of the axial, transverse displacement and rotation angle are considered by means of the first-order approximation coupling （FOAC） model theory, in which large deformations and the centrifugal stiffening effects are considered. Three kinds of systems are established respectively, which are a structure without piezoelectric layer, with piezoelectric layer in open circuit and closed circuit. Several simulations based on simplified models are presented to show the differences in characteristics between structures with and without the tip mass, between smart beams in closed and open circuit, and between the centrifugal effects in high speed rotating state or not. The last simulation calculates the dynamic response of the structure subjected to external electrical loading.

Hybrid natural element method for large deformation elastoplasticity problems

We present the hybrid natural element method（HNEM） for two-dimensional elastoplastic large deformation problems. Sibson interpolation is adopted to construct the shape functions of nodal incremental displacements and incremental stresses. The incremental form of Hellinger–Reissner variational principle for elastoplastic large deformation problems is deduced to obtain the equation system. The total Lagrangian formulation is used to describe the discrete equation system.Compared with the natural element method（NEM）, the HNEM has higher computational precision and efficiency in solving elastoplastic large deformation problems. Some numerical examples are selected to demonstrate the advantage of the HNEM for large deformation elastoplasticity problems.

A dynamic large-deformation particle finite element method for geotechnical applications based on Abaqus

In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.

Elastoplastic Large Deformation Using Meshless Integral Method

In this paper, the meshless integral method based on the regularized boundary integral equation [1] has been extended to analyze the large deformation of elastoplastic materials. The updated Lagrangian governing integral equation is obtained from the weak form of elastoplasticity based on Green-Naghdi’s theory over a local sub-domain, and the moving least-squares approximation is used for meshless function approximation. Green-Naghdi’s theory starts with the additive decomposition of the Green-Lagrange strain into elastic and plastic parts and considers aJ2elastoplastic constitutive law that relates the Green-Lagrange strain to the second Piola-Kirchhoff stress. A simple, generalized collocation method is proposed to enforce essential boundary conditions straightforwardly and accurately, while natural boundary conditions are incorporated in the system governing equations and require no special handling. The solution algorithm for large deformation analysis is discussed in detail. Numerical examples show that meshless integral method with large deformation is accurate and robust.

Application of MLPG in large deformation analysis

Two-dimensional large deformation analysis of hyperelastic and elasto-plastic solids based on the Meshless Local Petrov-Galerkin method （MLPG） is presented. A material configuration based the nonlinear MLPG formulation is introduced for the large deformation analysis of both path-dependent and path-independent materials. The supports of the MLS approximation functions cover the same sets of nodes during material deformation, thus the shape function needs to be computed only in the initial stage. The multiplicative hyperelasto-plastic constitutive model is adopted to avoid objective time integration for stress update in large rota- tion. With this constitutive model, the computational formulations for path-dependent and path-independent materials become identical. Computational efficiency of the nonlinear MLPG method is discussed and optimized in several aspects to make the MLPG an O（N） algorithm. The numerical examples indicate that the MLPG method can solve large deformation problems accurately. Moreover, the MLPG computations enjoy better convergence rate than the FEM under very large particle distortion.

Numerical modeling of large deformation and nonlinear frictional contact of excavation boundary of deep soft rock tunnel

Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During the excavation process of a deep soft rock tunnel, the rock wall may be compacted due to large deformation. In this paper, the technique to address this problem by a two-dimensional (2D) finite element software, large deformation engineering analyses software (LDEAS 1.0), is provided. By using the Lagrange multiplier method, the kinematic constraint of non-penetrating condition and static constraint of Coulomb friction are introduced to the governing equations in the form of incremental displacement. The numerical example demonstrates the efficiency of this technology. Deformations of a transportation tunnel in inclined soft rock strata at the depth of 1 000 m in Qishan coal mine and a tunnel excavated to three different depths are analyzed by two models, i.e. the additive decomposition model and polar decomposition model. It can be found that the deformation of the transportation tunnel is asymmetrical due to the inclination of rock strata. For extremely soft rock, large deformation can converge only for the additive decomposition model. The deformation of surrounding rocks increases with the increase in the tunnel depth for both models. At the same depth, the deformation calculated by the additive decomposition model is smaller than that by the polar decomposition model.

Effects of liquefaction-induced large lateral ground deformation on pile foundations

The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,numerical methods involving finite difference approach of pile in liquefaction-induced lateral spreading ground were derived and implemented into a finite difference program.Based on the monotonic loading tests on saturated sand after liquefaction,the liquefaction lateral deformation of the site where group piles are located was predicted.The effects of lateral ground deformation after liquefaction on a group of pile foundations were studied using the fmite difference program mentioned above,and the failure mechanism of group piles in liquefaction-induced lateral spreading ground was obtained.The applicability of the program was preliminarily verified.The results show that the bending moments at the interfaces between liquefied and non-liquefied soil layers are larger than those at the pile＇s top when the pile＇s top is embedded.The value of the additional static bending moment is larger than the peak dynamic bending moment during the earthquake,so in the pile foundation design,more than the superstructure＇s dynamics should be considered and the effect of lateral ground deformation on pile foundations cannot be neglected.

Experimental analysis of crack tip fields in rubber materials under large deformation

A three-nested-deformation model is proposed to describe crack-tip fields in rubber-like materials with large deformation. The model is inspired by the distribution of the measured in-plane and out-of-plane deformation. The in- plane displacement of crack-tip fields under both Mode 1 and mixed-mode （Mode I-II） fracture conditions is measured by using the digital Moir6 method. The deformation character- istics and experimental sector division mode are investigated by comparing the measured displacement fields under differ- ent fracture modes. The out-of-plane displacement field near the crack tip is measured using the three-dimensional digital speckle correlation method.

Static-deformation based fault diagnosis for damping spring of large vibrating screen

Based on the statics theory, a novel and feasible twice-suspended-mass method(TSMM) was proposed to deal with the seldom-studied issue of fault diagnosis for damping springs of large vibrating screen(LVS). With the static balance characteristic of the screen body/surface as well as the deformation compatibility relation of springs considered, static model of the screen surface under a certain load was established to calculate compression deformation of each spring. Accuracy of the model was validated by both an experiment based on the suspended mass method and the properties of the 3D deformation space in a numerical simulation. Furthermore, by adopting the Taylor formula and the control variate method, quantitative relationship between the change of damping spring deformation and the change of spring stiffness, defined as the deformation sensitive coefficient(DSC), was derived mathematically, from which principle of the TSMM for spring fault diagnosis is clarified. In the end, an experiment was carried out and results show that the TSMM is applicable for diagnosing the fault of single spring in a LVS.

UNSYMMETRICAL LARGE DEFORMATION PROBLEM OF ORTHOTROPIC PLATES

Based upon the theory of ardsotropic plates, the unsymmetrical large deformation equations of orthotropic circular plates were derived. By using Fourier series, the partial differential equations of this problem can be transformed into sets of nonlinear differential equations . And the procedure to solve the problem using the iterative method is given .

Differentiai-algebraic approach to large deformation analysis of frame structures subjected to dynamic loads

A nonlinear mathematical model for the analysis of large deformation of frame structures with discontinuity conditions and initial displacements, subject to dynamic loads is formulated with arc-coordinates. The differential quadrature element method （DQEM） is then applied to discretize the nonlinear mathematical model in the spatial domain, An effective method is presented to deal with discontinuity conditions of multivariables in the application of DQEM. A set of DQEM discretization equations are obtained, which are a set of nonlinear differential-algebraic equations with singularity in the time domain. This paper also presents a method to solve nonlinear differential-algebra equations. As application, static and dynamical analyses of large deformation of frames and combined frame structures, subjected to concentrated and distributed forces, are presented. The obtained results are compared with those in the literatures. Numerical results show that the proposed method is general, and effective in dealing with disconti- nuity conditions of multi-variables and solving differential-algebraic equations. It requires only a small number of nodes and has low computation complexity with high precision and a good convergence property.

Analysis of dynamical response to large deformation of piles with initial displacements

In this paper, a nonlinear mathematical model for analyzing dynamical response to the large deformation of piles with initial displacements is firstly established with the arc-coordinate, and it is a set of nonlinear integral-differential equa- tions, in which, the Winkeler model is used to simulate the resistance of the soil to the pile. Secondly, a set of new auxiliary functions are introduced. The differential-integral equations are transformed into a set of nonlinear differential equations, and the differential quadrature method （DQM） and the finite difference method （FDM） are applied to discretize the set of nonlinear equations in the spatial and time domains, respectively. Then, the Newton-Raphson method is used to solve the set of discretization algebraic equations at each time step. Finally, numerical examples are presented, and the dynamical re- sponses to the deformation of piles, including configuration, bending moment and shear force, are graphically illuminated. In calculation, two types of initial displacements and dynamical loads are applied, and the effects of parameters on the dynamical responses of piles are analyzed in detail.

THE PROBLEM OF LARGE DEFORMATION IN SOFTROCK ENGINEERING AND PRACTICAL ANALYSIS OF FLOOR-HEAVING

At present, the mechanics theories studying softrock engineering generally depend onlinear small-deformation hypothesis of chosital mechanics. Although these theorics can be considered with the physical non-linear features of softrock, it is still an approximate theory of geometricsmall-deformation. Because of the specific characteristics of medium environment, the problem ofsoftrock engineering should be thought as large deformation. This article wili prove the advantagesof large deformation theory in solving softrock problem with an example of the No. 2 pit of NalongCoal Mine. This will provide a beneficai method for the studying of large deformation mechanics ofsoftrock engineering.

NON-SYMMETRICAL LARGE DEFORMATION OF A SHALLOW THIN SPHERICAL SHELL

By the modified iteration method,in this paper,non=symmetrical large deflectionof a shallow spherical shell is discussed. We solve the second-order approximate analytical solution of the deflection of a shallow spherical shell subjected to linearliquid loads. and portray the characteristic curves of load-deflection on a perturbing point. With this paper's method, the similar questions of other kind of shell can bediscussed. Through the examples.we discuss the large defection of a plane and shallow spherical shells with different initial deflections.

MICRO-EXPERIMENTAL STUDY OF LARGE PLASTIC DEFORMATION FOR HIGH-PURITY POLYCRYSTALLINE ALUMINUM

This paper aims at revealing various micro- deformation characteristics, such as crystalline slip and grain boundary slide, which are recorded under scanning electronic microscope for high-purity aluminum tensile specimen at room temperature. These experimental data provide us direct evidence for shear localization near the grain boundary network via multi- directional grain boundary slide. The nonuniform deformation induced in the grain interiors would have decisive effect on the plastic flow and failure of polycrystalline materials.