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

Dynamic 3D axisymmetric problems in continuously non-homogeneous piezoelectric solids

The meshless local Petrov-Galerkin (MLPG) method is used to analyze transient dynamic problems in 3D axisymmetric piezoelectric solids with continuously inhomogeneous material properties. Both mechanical and thermal loads are considered here. A 3D axisymmetric body is created by rotation of a cross section around an axis of symmetry. Axial symmetry of geometry and boundary conditions reduces the original 3D boundary value problem into a 2D problem. The cross section is covered by small circular sub-domains surrounding nodes randomly spread over the analyzed domain. A unit step function is chosen as test function, in order to derive local integral equations on the boundaries of the chosen sub-domains, called local boundary integral equations (LBIE). These integral formulations are either based on the Laplace transform technique or the time-difference approach. The local integral equations are non-singular and take a very simple form, despite of inhomogeneous and anisotropic material behaviour across the analyzed structure. Spatial variation of all physical fields (or of their Laplace transforms) at discrete time instants are approximated on the local boundary and in the interior of the sub-domain by means of the moving least-squares (MLS) method. The Stehfest algorithm is applied for the numerical Laplace inversion, in order to retrieve the time-dependent solutions.     

竖向荷载作用下半透水边界固结问题解

本文基于Biot固结理论,采用积分变换的方法,获得了竖向点荷载、环形荷载作用下半透水边界固结问题基本解,并求得圆形分布荷载作用下的固结解,计算分析了圆形分布荷载作用下的竖向位移、总应力和孔压的时间效应.     

饱和土体中压力隧洞衬砌-土相互作用解析研究

基于Biot理论，采用渗流一力学耦合模型研究分析了内水压力作用下饱和土体中压力隧洞衬砌一土的相互作用问题。假定衬砌和土体均为饱和多孔介质，且衬砌和土体完全接触。运用积分变换理论，在hplace变换域中得到了衬砌和土体中的应力、位移和孔隙水压力解答，并利用bplace数值逆变换得到时域中的解。文末的算例分析结果表明：(1)采用渗流一力学耦合模型能较好地反映隧洞衬砌与土体相互作用中应力和变形的随时间变化过程；(2)衬砌和土体的相对刚度对隧洞的计算结果有很大影响。     

Generalized solutions of transient thermal shock problem with bounded boundaries

In this work, the generalized thermoelastic solutions with bounded boundaries for the transient shock problem are proposed by an asymptotic method. The governing equations are taken in the context of the generalized thermoelasticity with one relaxation time (L–S theory). The general solutions for any set of boundary conditions are obtained in the physical domain by the Laplace transform techniques. The corresponding asymptotic solutions for a thin plate with finite thickness, subjected to different sudden temperature rises in its two boundaries, are obtained by means of the limit theorem of Laplace transform. In the context of these asymptotic solutions, two specific problems with different boundary conditions have been conducted. The distributions of displacement, temperature and stresses, as well as the propagations, intersections and reflections of two elastic waves, named as thermoelastic wave and thermal wave separately, are obtained and plotted. These results are agreed with the results obtained in the existing literatures.     

Dynamic analysis of beams on viscoelastic foundation

This study is intended to analyze dynamic behavior of beams on Pasternak-type viscoelastic foundation subjected to time-dependent loads. The Timoshenko beam theory is adopted in the derivation of the governing equation. Ordinary differential equations in scalar form obtained in the Laplace domain are solved numerically using the complementary functions method to calculate exactly the dynamic stiffness matrix of the problem. The solutions obtained are transformed to the real space using the Durbin's numerical inverse Laplace transform method. The dynamic response of beams on viscoelastic foundation is analyzed through various examples.     

The extended Durbin method and its application for piezoelectric wave propagation problems

Based on the well-known Durbin method, an efficient numerical method was developed for the inversion of the two-sided Laplace transform. The accuracy of the method was verified using examples. As an application of the method, transient elastic waves propagating in a two-layered piezoelectric medium subjected to anti-plane concentrated loading and in-plane electric displacement loading were investigated. One-sided and two-sided Laplace transforms were applied to determine the shear stresses and electric displacements in the double Laplace transform domain. Subsequently, the Durbin method for one-sided Laplace transform inversion and the extended Durbin method for two-sided Laplace transform inversion were used to implement the numerical inversions. Additionally, the numerical results of the transient stresses and electric displacements were evaluated and discussed. It showed that the arrival time of transient waves satisfies physical phenomena, and the transient solution oscillates near the static solution and rapidly approximates the static solution.     

A solution for an isotropic sector under anti-plane shear loadings

In this study, the problem of an isotropic sector subjected to anti-plane shear loadings is investigated. The loadings were applied to the arc of the sector, and the radial edges of the sector were under traction-free or fixed conditions. Depending on these conditions, three problems, namely, free–free, fixed–free, and fixed–fixed edges were studied. A procedure using the finite Mellin transform combined with the Laplace transform was proposed for solving these problems. Explicit closed form solutions for the displacement and stress fields throughout the sector were obtained. The stress intensity factor (SIF) for each problem was analyzed using the obtained stress fields. It was determined that the SIF disappeared under the special condition of a fixed–fixed edge. Other special cases having anti-symmetric conditions were deduced from the derived solutions, and the results of these verified those cited in the literature as well as those obtained using finite element analysis (FEA).     

Laplace Transform and New Mathematical Theory of Viscoelasticity

In this paper generalized variational principles in thesense of the Laplace transform for viscoelastic problems arederived. Then mathematical theory of viscoelasticity ingeneralized Hardy spaces and in weighted anisotropic Sobolevspaces and spectral theory of corresponding non-selfadjointoperators is elaborated. Finally the Laplace transform--FEMfor numerical analysis of time-dependent problems of themathematical physics is proposed and analysed.     

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(1 D)consolidation for unsaturated soils with a semi-permeable drainage boundary are presented. 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 pressures 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 conducted on the pore-air and pore-water pressures at different ratios(the air permeability coefficient to the water permeability coefficient) and depths.     

Unsteady Couette flow of a micropolar fluid with slip

The unsteady Couette flow of an isothermal incompressible micropolar fluid between two infinite parallel plates is investigated. The motion of the fluid is produced by a time-dependent impulsive motion of the lower plate while the upper plate is set at rest. A linear slip, of Basset type, boundary condition on both plates is used. Two particular cases are discussed; in the first case we have assumed that the plate moves with constant speed and in the second case we have supposed that the plate oscillates tangentially. The solution of the problem is obtained in the Laplace transform domain. The inversion of the Laplace transform is carried out numerically using a numerical method based on Fourier series expansion. Numerical results are represented graphically for the velocity, microrotation, and volume flux for various values of the time, slip and micropolar parameters.     

Transient analysis of viscoelastic helical rods subject to time-dependent loads

The transient analysis of viscoelastic helical rods subject to time-dependent loads are examined in the Laplace domain. The governing equations for naturally twisted and curved spatial rods obtained using the Timoshenko beam theory are rewritten for cylindrical helical rods. The curvature of the rod axis, effect of rotary inertia and, shear and axial deformations are considered in the formulation. The material of the rod is assumed to be homogeneous, isotropic and linear viscoelastic. The viscoelastic constitutive equations are written in the Boltzmann–Volterra form. Ordinary differential equations in canonical form obtained in the Laplace domain are solved numerically using the complementary functions method to calculate the dynamic stiffness matrix of the problem. The solutions obtained are transformed to the real space using an appropriate numerical inverse Laplace transform method. Numerical results for quasi-static and dynamic response of viscoelastic models are presented in the form of graphics.     

Linear viscoelastic analysis of a semi-infinite porous medium

This paper considers the problem of a semi-infinite, isotropic, linear viscoelastic half-plane containing multiple, non-overlapping circular holes. The sizes and the locations of the holes are arbitrary. Constant or time dependent far-field stress acts parallel to the boundary of the half-plane and the boundaries of the holes are subjected to uniform pressure. Three types of loading conditions are assumed at the boundary of the half-plane: a point force, a force uniformly distributed over a segment, a force uniformly distributed over the whole boundary of the half-plane. The solution of the problem is based on the use of the correspondence principle. The direct boundary integral method is applied to obtain the governing equation in the Laplace domain. The unknown transformed displacements on the boundaries of the holes are approximated by a truncated complex Fourier series. A system of linear equations is obtained by using a Taylor series expansion. The viscoelastic stresses and displacements at any point of the half-plane are found by using the viscoelastic analogs of Kolosov–Muskhelishvili’s potentials. The solution in the time domain is obtained by the application of the inverse Laplace transform. All the operations of space integration, the Laplace transform and its inversion are performed analytically. The method described in the paper allows one to adopt a variety of viscoelastic models. For the sake of illustration only one model in which the material responds as the standard solid in shear and elastically in bulk is considered. The accuracy and efficiency of the method are demonstrated by the comparison of selected results with the solutions obtained by using finite element software ANSYS.     

考虑损伤时带微孔粘弹性体的应力分布

本文根据粘弹性理论和带空沿材料的线性理论，利用Laplace变换及其逆变换，给出了粘弹性损伤材料的一种本构模型。应用这种本构关系。讨论了带损伤的圆环形板或圆筒体在内外压力作用下的平面问题，得到了圆孔边缘附近的应力场和损伤增量场的分布情况，同时根据Laplace变换的终值定理，得到圆环或圆筒的终态应力和损伤增量的分布，并分析了工程中的一个典型的例子。得到了最终损伤的一种近似分布。可为工程设计提供一定的参考。     

空间杆系结构非线性动力稳定性分析的一种实用算法

本文对任意空间杆系结构进行非线性动力稳定性分析。     

The transient fracture behavior for a functionally graded layered structure subjected to an in-plane impact load

The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally graded interlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors (DSIFs) in time domain are studied and some meaningful conclusions are obtained.The project supported by the National Science Foundation for Excellent Young Investigators (10325208), the National Natural Science Foundation of China (10432030) and the China Postdoctoral Science Foundation (2004036018)The English text was polished by Ron Marshall.     

A semi-analytical approach to the study of an elastic circular cylinder confined in a cylindrical fluid domain subjected to small-amplitude transient motions

This paper deals with the transient motions experienced by an elastic circular cylinder in a cylindrical fluid domain initially at rest and subjected to small-amplitude imposed displacements. Three fluid models are considered, namely potential, viscous and acoustic, to cover different fluid–structure interaction regimes. They are derived here from the general compressible Navier–Stokes equations by a formal perturbation method so as to underline their links and ranges of validity a priori. The resulting fluid models are linear owing to the small-amplitude-displacement hypothesis. For simplicity, the elastic flexure beam model is chosen for the circular cylinder dynamics. The semi-analytical approach used here is based on the methods of Laplace transform in time, in vacuo eigenvector expansion with time-dependent coefficients for the transverse beam displacement and separation of variables for the fluid. Moreover, the viscous case is handled with a matched asymptotic expansion performed at first order. The projection of the fluid forces on the in vacuo eigenvectors leads to a fully coupled system involving the modal time-dependent displacement coefficients. These coefficients are then obtained by matrix inversion in the Laplace domain and fast numerical inversion of the Laplace transform. The three models, written in the form of convolution products, are described through the analysis of their kernels, involving both the wave propagation phenomena in the fluid domain and the beam elasticity. Last, the three models are illustrated for a specific imposed motion mimicking shock loading. It is shown that their combination permits coverage of a broad range of motions.     

Theoretical analysis,numerical calculation and experimental measurement of transient elastic waves in strips subjected to dynamic loadings

In this study, the transient response of an elastic strip subjected to dynamic in-plane loadings on the surface is investigated in detail. One of the objectives of this study is to develop an effective analytical method for determining transient solutions in a strip. By applying Laplace transform, the analytical solution in the transformed domain is derived and expressed in matrix form. The solution is then decomposed into infinite wave groups in which the multiple reflected waves with the same reflection are involved. Each multi-reflected wave can be identified by a coding method and be verified by the theory of generalized ray. The inverse transform is performed by using the well-known Cagniard method. The transient solutions in time domain for stresses and displacements are expressed in a closed form and are discussed in detail by an example. The experimental results show that the early time transient responses of displacements on the surface agree very well with the numerical calculations based on the theoretical solutions.     

Thermoelastic transient response of an infinitely long annular multilayered cylinder

Thermoelastic transient response of multilayered annular cylinders of infinite lengths subjected to known temperature at traction-free inner and outer surfaces are considered. A method based on the Laplace transformation and finite difference method has been developed to analyze the thermoelasticity problem. Using the Laplace transform with respect to time, the general solutions of the governing equation are obtained in transform domain. The solution is obtained by using the matrix similarity transformation and inverse Laplace transform. Solutions for the temperature and thermal stress distributions in a transient state were obtained. It was found that the temperature distribution, the displacement and the thermal stresses change slightly as time increases. There is no limit of number of annular layers of the cylinder in the presented computational procedures.     

含阻尼效应的黏弹性杆稳态响应的解析解

基于格林函数法和叠加原理,给出了简谐激励下开尔文黏弹性杆稳态响应的解析解。运用拉普拉斯变换法得到了任意边界条件下黏弹性杆强迫振动的格林函数解。设计了一个统一方案,将杆的各种边界表示为矩阵形式,便于确定格林函数中的待定常数。算例表明,材料黏性阻尼和外部阻尼导致杆位移在空间和时间上不可分离。