Semi-analytical solution for mechanical analysis of tunnels crossing strike-slip fault zone considering nonuniform fault displacement and uncertain fault plane position

The tunnel subjected to strike-slip fault dislocation exhibits severe and catastrophic damage.The existing analysis models frequently assume uniform fault displacement and fixed fault plane position.In contrast,post-earthquake observations indicate that the displacement near the fault zone is typically nonuniform,and the fault plane position is uncertain.In this study,we first established a series of improved governing equations to analyze the mechanical response of tunnels under strike-slip fault dislocation.The proposed methodology incorporated key factors such as nonuniform fault displacement and uncertain fault plane position into the governing equations,thereby significantly enhancing the applicability range and accuracy of the model.In contrast to previous analytical models,the maximum computational error has decreased from 57.1%to 1.1%.Subsequently,we conducted a rigorous validation of the proposed methodology by undertaking a comparative analysis with a 3D finite element numerical model,and the results from both approaches exhibited a high degree of qualitative and quantitative agreement with a maximum error of 9.9%.Finally,the proposed methodology was utilized to perform a parametric analysis to explore the effects of various parameters,such as fault displacement,fault zone width,fault zone strength,the ratio of maximum fault displacement of the hanging wall to the footwall,and fault plane position,on the response of tunnels subjected to strike-slip fault dislocation.The findings indicate a progressive increase in the peak internal forces of the tunnel with the rise in fault displacement and fault zone strength.Conversely,an augmentation in fault zone width is found to contribute to a decrease in the peak internal forces.For example,for a fault zone width of 10 m,the peak values of bending moment,shear force,and axial force are approximately 46.9%,102.4%,and 28.7% higher,respectively,compared to those observed for a fault zone width of 50 m.Furthermore,the position of the peak internal forces is influenced by variations in the ratio of maximum fault displacement of the hanging wall to footwall and the fault plane location,while the peak values of shear force and axial force always align with the fault plane.The maximum peak internal forces are observed when the footwall exclusively bears the entirety of the fault displacement,corresponding to a ratio of 0:1.The peak values of bending moment,shear force,and axial force for the ratio of 0:1 amount to approximately 123.8%,148.6%,and 111.1% of those for the ratio of 0.5:0.5,respectively.

A SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR SOLVING 2-D SOUND-STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.

Semi-analytical solutions to one-dimensional consolidation for unsaturated soils with semi-permeable drainage boundary

The semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations （PDFs）, which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pres- sures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are con- ducted on the pore-air and pore-water pressures at different ratios （the air permeability coefficient to the water permeability coefficient） and depths.

General semi-analytical solutions to one-dimensional consolidation for unsaturated soils

This paper presents general semi-analytical solutions to Fredlund and Hasan＇s one-dimensional （1D） consolidation equations for unsaturated soils subject to different initial conditions, homogeneous boundaries and time-dependent loadings. Two variables are introduced to transform the two-coupled governing equations of pore-water and poreair pressures into an equivalent set of partial differential equations （PDEs）, which are solved with the Laplace transform method. The pore-water and pore-air pressures and settlement are obtained in the Laplace transform domMn. The Crump＇s method is used to perform inverse Laplace transform to obtain the solutions in the time domain. The present solutions are more general in practical applications and show good agreement with the previous solutions in the literature.

Consolidation of partially saturated ground improved by impervious column inclusion:Governing equations and semi-analytical solutions

This study focuses on the consolidation behavior and mathematical interpretation of partially-saturated ground improved by impervious column inclusion.The constitutive relations for soil skeleton,pore air and pore water for partially saturated soils are proposed in the context of partially-saturated ground improved by impervious column inclusion.Settlement equation and dissipation equations of excess pore air/water pressures for a partially saturated improved ground are then derived.The semi-analytical solutions for ground settlement and pore pressure dissipation are then obtained through the Laplace transform and validated by the existing solutions for two special cases in the literature and the numerical results obtained from the finite difference method.A series of parametric studies is finally conducted to investigate the influence of some key factors on consolidation of partially saturated ground improved by impervious column inclusion.Based on the parametric study,it can be found that a higher value of the area replacement ratio or modulus of the pile results in a longer dissipation time of excess pore air pressure(PAP),a shorter dissipation time of excess pore water pressure(PWP),and a lower normalized settlement.

Hamiltonian parametric element and semi-analytical solution for smart laminated plates

Based on the Hellinger-Reissner （H-R） mixed variational principle for piezoelectric material, a unified 4-node Hamiltonian isoparametric element of anisotropy piezoelectric material is established. A new semi-analytical solution for the natural vibration of smart laminated plates and the transient response of the laminated cantilever with piezoelectric patch is presented. The major steps of mathematical model are as follows： the piezoelectric layer and host layer of laminated plate are considered as unattached three-dimensional bodies and discretized by the Hamiltonian isoparametric elements. The control equation of whole structure is derived by considering the compatibility of generalized displacements and generalized stresses on the interface between layers. There is no restriction for the side-face geometrical boundaries, the thickness and the number of layers of plate by the use of the present isoparametric element. Present method has wide application area.

Semi-analytical solution to one-dimensional consolidation in unsaturated soils

In this paper, a series of semi-analytical solutions to one-dimensional consolidation in unsaturated soils are obtained. The air governing equation by Fredlund for unsaturated soils consolidation is simplified. By applying the Laplace transform and the Cayley-Hamilton theorem to the simplified governing equations of water and air, Darcy＇s law, and Fick＇s law, the transfer function between the state vectors at top and at any depth is then constructed. Finally, by the boundary conditions, the excess pore-water pressure, the excess pore-air pressure, and the soil settlement are obtained under several kinds of boundary conditions with the large-area uniform instantaneous loading. By the Crump method, the inverse Laplace transform is performed, and the semi-analytical solutions to the excess pore-water pressure, the excess pore-air pressure, and the soils settlement are obtained in the time domain. In the case of one surface which is permeable to air and water, comparisons between the semi-analytical solutions and the analytical solutions indicate that the semi-analytical solutions are correct. In the case of one surface which is permeable to air but impermeable to water, comparisons between the semi-analytical solutions and the results of the finite difference method are made, indicating that the semi-analytical solution is also correct.

A Semi-Analytical Potential Solution for Wave Resonance in Gap Between Floating Box and Vertical Wall

Based on potential flow theory, a dissipative semi-analytical solution is developed for the wave resonance in the narrow gap between a fixed floating box and a vertical wall by using velocity potential decompositions and matched eigenfunction expansions. The energy dissipation near the box is modelled in the potential flow solution by introducing a quadratic pressure loss condition on the gap entrance. Such a treatment is inspired by the classical local head loss formula for the sudden change of cross section in channel flow, where the energy dissipation is assumed to be proportional to the square of local velocity for high Reynolds number flows. The dimensionless energy loss coefficient is calibrated based on experimental data. And it is found to be insensitive to the incident wave height and wave frequency. With the calibrated energy loss coefficient, the resonant wave height in gap and the reflection coefficient are calculated by the present dissipative semi-analytical solution. The predictions are in good agreement with experimental data. Case studies suggest that the maximum relative energy dissipation occurs near the resonant frequency, which leads to the minimum reflection coefficient. The horizontal wave forces on the box and the vertical wall attain also maximum values near the resonant frequency, while the vertical wave force on the box decreases abruptly there to a small value.

A semi-analytical solution for electric double layers near an elliptical cylinder

A theoretical analysis on the electric double layer formed near the surface of an infinite cylinder with an elliptical cross section and a prescribed electric potential in an ionic conductor was performed using the linearized Gouy–Chapman theory. A semi-analytical solution in terms of the Mathieu functions was obtained. The distributions of the electric potential, cations, anions, and electric field were calculated. The effects of various physical and geometric parameters were examined. The fields vary rapidly near the elliptical boundary and are nearly uniform at far field. Electric field concentrations were found at the ends of the semi-major and semi-minor axes of the ellipse. These concentrations are sensitive to the physical and geometric parameters.

A semi-analytical solution for frost heave prediction of clay soil

Frost heave is one of the main freezing problems for construction in permafrost regions.The Konrad-Morgenstern segregation potential（SP） model is being used in practice for frost heave using numerical techniques.However,the heat release from in-situ and migrated water in the freezing zone could result in some numerical instability,so the simulation of frost fringe is not ideal.In this study,a semi-analytical solution is developed for frost heave prediction of clay soil.The prediction results to the two tests with different freezing mode with clay soil agree well with the tested behavior,which indicates the feasibility of the solution.

Solution approximations for a mathematical model of relativistic electrons with beta derivative

The main aim of this paper is to obtain the exact and semi-analytical solutions of the nonlinear Klein-Fock-Gordon(KFG)equation which is a model of relativistic electrons arising in the laser thermonuclear fusion with beta derivative.For this purpose,both the modified extended tanh-function(mETF)method and the homotopy analysis method(HAM)are used.While applying the mETF the chain rule for beta derivative and complex wave transform are used for obtaining the exact solution.The advantage of this procedure is that discretization or normalization is not required.By applying the mETF,the exact solutions are obtained.Also,by applying the HAM semi-analytical results for the considered equation are acquired.In HAM?curve gives us a chance to find the suitable value of the for the convergence of the solution series.Also,comparative graphical representations are given to show the effectiveness,reliability of the methods.The results show that the m ETF and HAM are reliable and applicable tools for obtaining the solutions of non-linear fractional partial differential equations that involve beta derivative.This study can bring a new perspective for studies on fractional differential equations.On the other hand,it can be said that scientists can apply the considered methods for different mathematical models arising in physics,chemistry,engineering,social sciences and etc.which involves fractional differentiation.Briefly the results may cause a new insight who studies on relativistic electron modelling.

Semi-analytic Solution of Steady Temperature Fields During Thin Plate Welding

Usually, it is very difficult to find out an analytical solution to thermal conduction problems during high temperature welding. Therefore, as an important numerical approach, the method of lines （MOLs） is introduced to solve the temperature field characterized by high gradients. The basic idea of the method is to semi-discretize the governing equation of the problem into a system of ordinary differential equations （ODEs） defined on discrete lines by means of the finite difference method, by which the thermal boundary condition with high gradients are directly embodied in formulation. Thus the temperature field can be obtained by solving the ODEs. As a numerical example, the variation of an axisymmetrical temperature field along the plate thickness can be obtained.

Closed-form solution of mid-potential between two parallel charged plates with more extensive application

Efficient calculation of the electrostatic interactions including repulsive force between charged molecules in a biomolecule system or charged particles in a colloidal system is necessary for the molecular scale or particle scale mechanical analyses of these systems. The electrostatic repulsive force depends on the mid-plane potential between two charged particles. Previous analytical solutions of the mid-plane potential, including those based on simplified assumptions and modern mathematic methods, are reviewed. It is shown that none of these solutions applies to wide ranges of interparticle distance from 0 to 10 and surface potential from 1 to 10. Three previous analytical solutions are chosen to develop a semi-analytical solution which is proven to have more extensive applications. Furthermore, an empirical closed-form expression of mid-plane potential is proposed based on plenty of numerical solutions. This empirical solution has extensive applications, as well as high computational efficiency.

Solution to 1-D consolidation of non-homogeneous soft clay

In this work, semi-analytical methods were used to solve the problem of 1-D consolidation of non-homogeneous soft clay with spatially varying coefficients of permeability and compressibility. The semi-analytical solution was programmed and then verified by comparison with the obtained analytical solution of a special case. Based on the results of some computations and comparisons with the 1-D homogeneous consolidation (by Terzaghi) and the 1-D non-linear consolidation theory (by Davis et al.) of soft clay, some diagrams were prepared and the relevant consolidation behavior of non-homogeneous soils is discussed. It was shown that the result obtained differs greatly from Terzaghi’s theory and that of the non-linear consolidation theory when the coefficients of permeability and compressibility vary greatly.

Stress analysis of anisotropic thick laminates in cylindrical bending using a semi-analytical approach

Semi-analytical elasticity solutions for bending of angle-ply laminates in cylindrical bending are presented using the state-space-based differential quadrature method (SSDQM). Partial differential state equation is derived from the basic equations of elasticity based on the state space concept. Then, the differential quadrature (DQ) technique is introduced to discretize the longitu- dinal domain of the plate so that a series of ordinary differential state equations are obtained at the discrete points. Meanwhile, the edge constrained conditions are handled directly using the stress and displacement components without the Saint-Venant principle. The thickness domain is solved analytically based on the state space formalism along with the continuity conditions at interfaces. The present method is validated by comparing the results to the exact solutions of Pagano’s problem. Numerical results for fully clamped thick laminates are presented, and the influences of ply angle on stress distributions are discussed.

HAMILTONIAN SYMPLECTIC SEMI-ANALYTICAL METHOD AND ITS APPLICATION IN INHOMOGENEOUS ELECTROMAGNETIC WAVEGUIDES

Dual vectors are applied in Hamilton system of applied mechanics. Electric and magnetic field vectors are the dual vectors in electromagnetic field. The Hamilton system method is introduced into the analysis of electromagnetism waveguide with inhomogeneous materials. The transverse electric and magnetic fields are regarded as the dual. The basic equations are solved in Hamilton system and symplectic geometry. With the Hamilton variational principle, the symplectic semi_analytical equations are derived and preserve their symplectic structures. The given numerical example demonstrates the solution of LSE (Longitudinal Section Electric) mode in a dielectric waveguide. (

Analytical solution and semi-analytical solution for anisotropic functionally graded beam subject to arbitrary loading

Analytical and semi-analytical solutions are presented for anisotropic functionally graded beams subject to an arbitrary load,which can be expanded in terms of sinusoidal series.For plane stress problems,the stress function is assumed to consist of two parts,one being a product of a trigonometric function of the longitudinal coordinate(x) and an undetermined function of the thickness coordinate(y),and the other a linear polynomial of x with unknown coefficients depending on y.The governing equations satisfied by these y-dependent functions are derived.The expressions for stresses,resultant forces and displacements are then deduced,with integral constants determinable from the boundary conditions.While the analytical solution is derived for the beam with material coefficients varying exponentially or in a power law along the thickness,the semi-analytical solution is sought by making use of the sub-layer approximation for the beam with an arbitrary variation of material parameters along the thickness.The present analysis is applicable to beams with various boundary conditions at the two ends.Three numerical examples are presented for validation of the theory and illustration of the effects of certain parameters.

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.

Theoretical investigation on axial cyclic performance of monopile in sands using interface constitutive models

Cyclic loads generated by environmental factors,such as winds,waves,and trains,will likely lead to performance degradation in pile foundations,resulting in issues like permanent displacement accumulation and bearing capacity attenuation.This paper presents a semi-analytical solution for predicting the axial cyclic behavior of piles in sands.The solution relies on two enhanced nonlinear load-transfer models considering stress-strain hysteresis and cyclic degradation in the pile-soil interaction.Model parameters are calibrated through cyclic shear tests of the sand-steel interface and laboratory geotechnical testing of sands.A novel aspect involves the meticulous formulation of the shaft loadtransfer function using an interface constitutive model,which inherently inherits the interface model’s advantages,such as capturing hysteresis,hardening,degradation,and particle breakage.The semi-analytical solution is computed numerically using the matrix displacement method,and the calculated values are validated through model tests performed on non-displacement and displacement piles in sands.The results demonstrate that the predicted values show excellent agreement with the measured values for both the static and cyclic responses of piles in sands.The displacement pile response,including factors such as bearing capacity,mobilized shaft resistance,and convergence rate of permanent settlement,exhibit improvements compared to non-displacement piles attributed to the soil squeezing effect.This methodology presents an innovative analytical framework,allowing for integrating cyclic interface models into the theoretical investigation of pile responses.

Dynamic characteristics of multi-span spinning beams with elastic constraints under an axial compressive force

A theoretical model for the multi-span spinning beams with elastic constraints under an axial compressive force is proposed.The displacement and bending angle functions are represented through an improved Fourier series,which ensures the continuity of the derivative at the boundary and enhances the convergence.The exact characteristic equations of the multi-span spinning beams with elastic constraints under an axial compressive force are derived by the Lagrange equation.The efficiency and accuracy of the present method are validated in comparison with the finite element method(FEM)and other methods.The effects of the boundary spring stiffness,the number of spans,the spinning velocity,and the axial compressive force on the dynamic characteristics of the multi-span spinning beams are studied.The results show that the present method can freely simulate any boundary constraints without modifying the solution process.The elastic range of linear springs is larger than that of torsion springs,and it is not affected by the number of spans.With an increase in the axial compressive force,the attenuation rate of the natural frequency of a spinning beam with a large number of spans becomes larger,while the attenuation rate with an elastic boundary is lower than that under a classic simply supported boundary.