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.

New semi-analytical approach for ice lens heaving during artificial freezing of fine-grained material

The calculation of frost heaving with ice lens formation is still not standard for construction projects using artificial ground freezing(AGF).In fine-grained material,ice lenses may initiate and lead to significant heaving at the ground surface,which should be considered in advance.However,the complex processes during ice lens formation are still not fully understood and difficult to capture in a simple approach.In the past,the semi-analytical approach of Konrad and Morgenstern used one soil constant,the“segregation potential(SP)”.It has been mainly and most successfully applied to the heave calculation of natural-induced soil freezing in cold regions.Its application to AGF has been so far unsuccessful.To solve this,a new semi-analytical approach is presented in this paper.It includes AGF conditions such as bottom-up freezing,temperature gradients to reach great freezing velocities,and a distinction between two freezing states.One is the freezing-up state until a certain frost body thickness is reached(thermal transient state),and the other is a holding phase where the frost body thickness is kept constant(thermal quasi-steady state).To test its ability,the results are applied to another freezing direction,the top-down freezing.The new approach is validated using two different frost-susceptible soils and,in total,50 frost heave tests.In the thermal transient region,where the SP is applicable,the two semi-analytical approaches are compared,showing improved performance of the current method by about 15%.

Semi-analytical solution for internal forces of tunnel lining with multiple longitudinal cracks

Longitudinal cracks on the tunnel lining significantly influence the performance of tunnels in operation.In this study,we propose a semi-analytical method that provides a simple and effective way to calculate the internal forces of tunnel linings with multiple cracks.The semi-analytical solution is obtained using structural analysis considering the flexural rigidity for the cracked longitudinal section of the tunnel lining.Then the proposed solution is verified numerically.Using the proposed method,the influences of the crack depth and the number of cracks on the bending moment and modified crack tip stress are investigated.With the increase in crack depth,the bending moment of lining scetion adjacent to the crack decreases,while the bending moment of lining scetion far away from the crack increases slightly.The more the number of cracks in a tunnel lining,the easier the new cracks initiated.

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.

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.

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.

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation efectiveness and precision, is presented for solving the acoustic radiation from a submerged infnite non-circular cylindrical shell stifened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work, besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The efects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efcient than the existing mixed FE-BE method.

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.

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.

Semi-analytical Monte Carlo simulation of laser-induced fluorescence propagation in an optically participating spray

A semi-analytical Monte Carlo (SMC) simulation was developed to simulate the propagation of laser-induced fluo- rescence (LIF) in an optically participating spray, which simultaneously exhibits spectrally dependent emission, anisotropic scattering, absorption, and re-emission. The SMC simulation is described and then applied to an experimental configuration of a cloud of polydisperse droplets composed of water and sulforhodamine B dye. In the SMC simulation, the collected LIF flux on the remote receptor element is calculated as the global contribution from the emissive source, single, twice, … and nth collision events in any sequence. The effects on the fluorescence photons propagation of spray parameters like the dye concentration, droplets concentration, and droplets size are examined. Three spectral bands representing different optical properties are chosen to analyze the interference of absorption, scattering and re-emission on the detected LIF flux. The obtained spectral LIF flux distribution on the receptor demonstrates a “red shift” phenomenon.

Using a semi-analytical model to retrieve Secchi depth in coastal and estuarine waters

Secchi depth(SD,m)is a direct and intuitive measure of water's transparency,which is also an indicator of water quality.In 2015,a semi-analytical model was developed to derive SD from remote sensing reflectance,thus able to provide maps of water's transparency in satellite images.Here an in-situ dataset(338 stations)is used to evaluate its potential ability to monitor water quality in the coastal and estuarine waters,with measurements covering the Zhujiang(Pearl)River Estuary,the Yellow Sea and the East China Sea where measured SD values span a range of 0.2–21.0 m.As a preliminary validation result,according to the whole dataset,the unbiased percent difference(UPD)between estimated and measured SD is 23.3%(N=338,R^2=0.89),with about 60%of stations in the dataset having relative difference(RD)≤20%,over 80%of stations having RD≤40%.Furthermore,by excluding the field data which with relatively larger uncertainties,the semi-analytical model yielded the UPD of 17.7%(N=132,R^2=0.92)with SD range of 0.2–11.0 m.In addition,the semi-analytical model was applied to Landsat-8 images in the Zhujiang River Estuary,and retrieved high-quality mapping and reliable spatial-temporal patterns of water clarity.Taking into account the uncertainties associated with both field measurements and satellite data processing,and that there were no tuning of the semi-analytical model for these regions,these findings indicate highly robust retrieval of SD from spectral techniques for such turbid coastal and estuarine waters.The results suggest it is now possible to routinely monitor coastal water transparency or visibility at high-spatial resolutions from measurements,like Landsat-8 and Sentinel-2 and newly launched Gaofen-5.

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.

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.

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.

Nonlinear dynamic response and buckling of laminated cylindrical shells with axial shallow groove based on a semi-analytical method

The effect of axial shallow groove on the nonlinear dynamic response and buckling of laminated cylindrical shells subjected to radial compression loading was investigated. Based on the first-order shear deformation theory （FSDT）, the nonlinear dynamic equations involving the transverse shear deformation and initial geometric imperfections were derived with the Hamilton philosophy. The axial shallow groove of the laminated composite cylindrical shell was treated as the initial geometric imperfections in the dynamic equations. A semi-analytical method of expanding displacements and loads along the circumferential direction and employing the finite difference method along the axial direction and in the time domain is used to solve the governing equations and obtain the dynamic response of the laminated shell. The B-R criterion was employed to determine the critical loads of dynamic buckling of the shell. The effects of the parameters of the shallow groove on the dynamic response and buckling were discussed in this paper and the results show that the axial shallow grooves greatly affect the dynamic response and buckling.

Construction of Crack Perturbation Model and Forward Semi-analytical Model of Attached Eddy Current Sensor

Fatigue damage monitoring is critical metallic structure health monitoring of aircraft.The sensor should be high sensitive,easy to be integrated into structure and well adaptable for poor working conditions.Therefore,an attached eddy current sensor with flexible plane is put forward and its characteristics are analyzed.By extracting material′s conductivity as the crack features,forward semi-analytical model is established and parameter optimizations are carried out.Crack perturbation model of attached eddy current sensor is constructed,and perturbation voltages of sensing channels under three-dimension structural crack are obtained.To verify the sensor′s performance,monitoring experiment on crack extension is conducted under condition of 3 MHz frequency.The validation experimental results show that perturbation model of 2A12-T4 aluminum alloy agrees well with experiment results,and perturbation model errors of four sensing channels are within 25%.The attached eddy current sensor is capable of testing the crack nondestructively and measuring the crack extension quantitatively with the accuracy of 1mm.

SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR DYNAMIC ANALYSIS OF FOUNDATION

A semi-analytical and semi-numerical method is proposed for the dynamic analysis of foundations. The Lamb＇s solution and the approximate formulae were used to establish the relation of the contact force and deflection between the foundation and soil. Therefore, the foundation can be separated from soil and analyzed by FEM as for the static cases. The plate can be treated as that the known forces are acting on the upper surface, and the contact pressure from soil can be represented as the deflection. So that only the plate needs to be divided into elements in the analysis. By this method, a series of vibration problems, including various shapes and rigidities of foundations, different excitation frequencies, were analyzed. Furthermore, it can be used for the embedded foundation. The numerical examples show that this method has simplicity, highly accurate and versatile. It is an effective method for the dynamic analysis of foundations.

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.

Numerical Study on Motion Responses of Two Parallel Ships During Underway Replenishment by the Semi-Analytical HOTP Method

Motion responses of two ships advancing parallel in waves with hydrodynamic interactions are investigated in this paper. Within the framework of the frequency-domain potential flow theory, a semi-analytical higher-order translating-pulsating source(HOTP) method is presented to solve the problems of coupled radiation and diffraction potential. The method employs nine-node bi-quadratic curvilinear elements to discretize the boundary integral equations(BIEs) constructed over the mean wetted surface of the two ship hulls. In order to eliminate the numerical oscillation, analytical quadrature formulas are derived and adopted to evaluate the integrals related to the Froudedependent part of the Green’s function along the horizontal direction in the BIEs. Based on the method, a numerical program is originally coded. Through the calculations of hydrodynamic responses of single ships, the numerical implementation is proved successful. Then the validated program is applied in the investigations on the hydrodynamic interactions of two identical Wigley Ⅲ hulls and the underway replenishment of a frigate and a supply ship in waves with and without stagger, respectively. The comparison between the present computed results with experimental data and numerical solutions of other methods shows that the semi-analytical HOTP method is of higher accuracy than the pulsating source Green’s function method with speed correction and better stability than the traditional HOTP method based on Gauss quadrature. In addition, for two ships with obviously different dimensions,the influence of hydrodynamic interactions on the smaller ship is found to be more noticeable than that on the larger ship, which leads to the differences between the motions of frigate with and without the presence of supply ship.

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. (