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

FORMULATION OF A SEMI-ANALYTICAL APPROACH BASEDON GURTIN VARIATIONAL PRINCIPLE FOR DYNAMICRESPONSE OF GENERAL THIN PLATES

4 semi-analytical approach for the dynamic response of general thin plates whichemployes finite element discretization in space domain and a series of representation intime demain is developed on the basis of Curtin variational principles.The formulationof time series is also investigated so that the dynamic response of plates with arbitraryshape and boundary constraints can be achieved with adequate accuracy.

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.

Semi-analytical solution for drained expansion analysis of a hollow cylinder of critical state soils

The expansion of a thick-walled hollow cylinder in soil is of non-self-similar nature that the stress/deformation paths are not the same for different soil material points.As a result,this problem cannot be solved by the common self-similar-based similarity techniques.This paper proposes a novel,exact solution for rigorous drained expansion analysis of a hollow cylinder of critical state soils.Considering stress-dependent elastic moduli of soils,new analytical stress and displacement solutions for the nonself-similar problem are developed taking the small strain assumption in the elastic zone.In the plastic zone,the cavity expansion response is formulated into a set of first-order partial differential equations(PDEs)with the combination use of Eulerian and Lagrangian descriptions,and a novel solution algorithm is developed to efficiently solve this complex boundary value problem.The solution is presented in a general form and thus can be useful for a wide range of soils.With the new solution,the non-self-similar nature induced by the finite outer boundary is clearly demonstrated and highlighted,which is found to be greatly different to the behaviour of cavity expansion in infinite soil mass.The present solution may serve as a benchmark for verifying the performance of advanced numerical techniques with critical state soil models and be used to capture the finite boundary effect for pressuremeter tests in small-sized calibration chambers.

A semi-analytical pressure and rate transient analysis model for inner boundary and propped fractures exhibiting dynamic behavior under long-term production conditions

The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.

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.

Enhancing perianal disease management with integrated physical and psychological approaches

This article provides a comprehensive analysis of the study by Hou et al,focusing on the complex interplay between psychological and physical factors in the postoperative recovery(POR)of patients with perianal diseases.The study sheds light on how illness perception,anxiety,and depression significantly influence recovery outcomes.Hou et al developed a predictive model that demonstrated high accuracy in identifying patients at risk of poor recovery.The article explores the critical role of pre-operative psychological assessment,highlighting the need for mental health support and personalized recovery plans in enhancing POR quality.A multidisciplinary approach,integrating mental health professionals with surgeons,anesthesiologists,and other specialists,is emphasized to ensure comprehensive care for patients.The study’s findings serve as a call to integrate psychological care into surgical practice to optimize outcomes for patients with perianal diseases.

Semi-analytical modeling of tokamak density evolution

Tokamak plasma density evolution is generally modeled by a diffusion-convection equation in cylindrical geometry. By using a semi-analytical approach, we solve such an equation for a given diffusion coefficient and inward convection velocity as an arbitrary function of the radial position. Through variable separation, a Sturm-Liouville-type eigenvalue problem is solved, thereby constructing a complete set of orthogonal eigenfunctions. Based on the decomposition of the solution, the initial function, and the source function in these eigenfunctions, several problems of practical interest about the density evolution are analyzed. They include the density evolution, with boundary density not being zero; the density profile with internal transport barrier; the damping profile during particle source being shut-down. Results are found to be qualitatively consistent with the tokamak experiments.

MODIFIED SCHEME BASED ON SEMI-ANALYTIC APPROACH FOR COMPUTING NON-PROBABILISTIC RELIABILITY INDEX

A new computation scheme proposed to tackle commensurate problems is devel- oped by modifying the semi-analytic approach for minimizing computational complexity. Using the proposed scheme, the limit state equations, usually referred to as the failure surface, are obtained from transformation of an interval variable to a normalized one. In order to minimize the computational cost, two algorithms for optimizing the calculation steps have been proposed. The monotonicity of the objective function can be determined from narrowing the scope of interval variables in normalized infinite space by incorporating the algorithms into the computational scheme. Two examples are used to illustrate the operation and computational efficiency of the approach. The results of these examples show that the proposed algorithms can greatly reduce the computation complexity without sacrificing the computational accuracy. The advantage of the proposed scheme can be even more efficient for analyzing sophistic structures.

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.

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.

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.

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.

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.

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.

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