磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

高速移动荷载下黏弹性半空间体的动力响应

分别以移动荷载和黏弹性半空间体模拟运动列车荷载和地基,分析了地基在运动列车作用下的动力响应.首先采用Green函数法求解黏弹性半空间体在各种移动荷载模式作用下的动力响应的解析解,包括恒常和简谐移动点源、线源和面源荷载.然后采用IFFT算法和自适应数值积分算法计算解析解中的二维积分,得到了包括低音速、跨音速和超音速移动荷载作用下位移的数值结果.最后分析了速度对位移的分布和最大值的影响,发现当速度大于Rayleigh波速时,位移发生显著变化.     

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

State-space approach to two-temperature generalized thermoelasticity without energy dissipation of medium subjected to moving heat source

In this work,a model of two-temperature generalized thermoelasticity without energy dissipation for an elastic half-space with constant elastic parameters is constructed.The Laplace transform and state-space techniques are used to obtain the general solution for any set of boundary conditions.The general solutions are applied to a specific problem of a half-space subjected to a moving heat source with a constant velocity.The inverse Laplace transforms are computed numerically,and the comparisons are shown in figures to estimate the effects of the heat source velocity and the two-temperature parameter.     

Integration of Singularities in FE/BE Analyses of Soil-foundation Interaction with Non-homogeneous Elastic Soils

This paper addresses the derivation of the boundary integral equations for a non-homogeneous elastic half-space subjected to constant surface tractions on an arbitrarily shaped area on the basis of the respective Green's functions. The type of non-homogeneity considered is a power law variation of Young's modulus with depth below the surface of the half-space. Two different methods, a contour integral and an integration-free approach are presented, applicable for arbitrarily and rectangular shaped boundary elements, respectively. In the first one the divergence theorem is used in order to reduce the integration of a two-dimensional surface element to an integration over the element's confining boundary only. In the second approach no integration at all is needed since the solution is found simply by evaluating functions to be determined at the boundaries of the loaded rectangle.     

Two-dimensional generalized thermoelasticity problem for a half-space subjected to ramp-type heating

In this paper, we will consider a half-space filled with an elastic material, which has constant elastic parameters. The governing equations are taken in a unified system from which the field equations for coupled thermoelasticity as well as for generalized thermoelasticity can be easily obtained as particular cases. A linear temperature ramping function is used to more realistically model thermal loading of the half-space surface. The medium is assumed initially quiescent. Laplace and Fourier transform techniques are used to obtain the general solution for any set of boundary conditions. The general solution obtained is applied to a specific problem of a half-space subjected to ramp-type heating. The inverse Fourier transforms are obtained analytically while the inverse Laplace transforms are computed numerically using a method based on Fourier expansion techniques. Some comparisons have been shown in figures to estimate the effect of the ramping parameter of heating with different theories of thermoelasticity.     

Finite amplitude one-dimensional wave propagation in a thermoelastic half-space

The problem of finite wave propagation in a nonlinearly thermoelastic half-space is considered. The surface of the half-space is subjected to a time-dependent thermal and normal mechanical loading. The solution is obtained by a numerical procedure, which is shown to furnish accurate results, and linear dynamic thermoelastic problems are obtained as special cases. The accuracy of the results is checked by comparison with some known analytical solutions which can be obtained in some special cases of both the linear and the nonlinear problems. In those cases where the solution contains shocks, it is shown that the numerical results satisfy the necessary jumps conditions which need to hold across such discontinuities.     

Dynamic Response of a Multilayered Poroelastic Half-Space to Harmonic Surface Tractions

The propagator matrix method is developed to study the dynamic response of a multilayered poroelastic half-space to time-harmonic surface tractions. In a cylindrical coordinate system, a method of displacement potentials is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to those equations are obtained. After that, the propagator matrix method and the vector surface harmonics are employed to derive the solutions for a multilayered poroelastic half-space subjected to surface tractions. It is known that the original propagator algorithm has the loss-of-precision problem when the waves become evanescent. At present, an orthogonalization procedure is inserted into the matrix propagation loop to avoid the numerical difficulty of the original propagator algorithm. Finally, a high-order adaptive integration method with continued fraction expansions for accelerating the convergence of the truncated integral is adopted to numerically evaluate the integral solutions expressed in terms of semi-infinite Hankel-type integrals with respect to horizontal wavenumber. Furthermore, to validate the present approach, the response of a uniform poroelastic half-space is examined using the formulation proposed in this article. It is shown that the numerical results computed with this approach agree well with those computed with the analytical solution of a uniform half-space.     

Three-dimensional dynamic Green’s functions for a multilayered transversely isotropic half-space

With the aid of a method of displacement potentials, an efficient and accurate analytical derivation of the three-dimensional dynamic Green’s functions for a transversely isotropic multilayered half-space is presented. Constituted by proper algebraic factorizations, a set of generalized transmission–reflection matrices and internal source fields that are free of any numerically unstable exponential terms are proposed for effective computations of the potential solution. Three-dimensional point-load Green’s functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The present formulations and solutions are analytically in exact agreement with the existing solutions given by Pak and Guzina (2002) for the isotropic case. For the numerical computation of the integrals, a robust and effective methodology which gives the necessary account of the presence of singularities including branch points and poles on the path of integration is laid out. A comparison with the existing numerical solutions for multilayered isotropic half-space is made to confirm the accuracy of the numerical solutions.     

竖向力作用在粘弹性土体内部时的解答与应用

选择三参量固体模型描述土体的粘弹性本构关系,利用半空间体内部受竖向集中力的Mindlin弹性解,根据粘弹性理论中的准静态弹性-弹粘性对应原理,推导了竖向力作用在粘弹性半无限土体内部的理论解。通过对应力和位移解答进行Laplace逆变换,给出了应力与位移的时域解。作为解答的应用,建立了粘弹性半无限体内部矩形面积上作用有三角形分布、均匀分布荷载时的粘弹性沉降计算公式。为了便于计算与工程应用,根据粘弹性理论解编制了计算程序。结果验证和实例分析表明,文中理论解是正确的,并为工程实际应用提供了理论依据。     

A nonlinearly thermoelastic half-space under time-dependent normal and shear loading

A nonlinearly thermoelastic half-space is subjected to combined time-dependent normal and shear loading. The solution is obtained by a numerical method which is shown to yield accurate results by comparison with some known analytical solutions which can be obtained in some special cases. When shocks are involved, it is shown that the numerical results satisfy all the Rankine-Hugoniot jump conditions as well as the entropy condition across the shock.     

Two-temperature generalized thermoelasticity under ramp-type heating by finite element method

In this work, a general finite element model is proposed to analyze transient phenomena in thermoelastic half-space filled with an elastic material, which has constant elastic parameters. The governing equations are taken in the context of the two-temperature generalized thermoelasticity theory (Youssef in IMA J. Appl. Math. 71(3):383–390, 2006). A linear temperature ramping function is used to more realistically model thermal loading of the half-space surface. The medium is assumed initially quiescent. A finite element scheme is presented for the high accuracy numerical purpose. The numerical solutions of the non-dimensional governing partial differential equations of the problem have been shown graphically and some comparisons have been shown in figures to estimate the effect of the ramping parameter of heating and the parameter of two-temperature.     

Elasticity solutions for a piezoelectric cone under concentrated loads at its apex

IntroductionTheproblemofaconesubjectedtoconcentratedloadsatitsapexisaclassicalprobleminthetheoryofelasticity.AnumberofscholarshavestUdiedtheproblem.LovereportedthesolutionstotheproblemofanisotropicconeunderconcentfatCdforcesatitsapex['].Lur'estudiedthisclassofproblemssystematicallybymeansofPapkovich-Neubergeneralsolution[2].LekniskiiandHu,byusingtheirrespectivegeneralsolutions,studiedcompressionandbendingproblemsofatransverselyisotropicconesubjectedtoaxialconcentfatedforcesandtfansverseconc…     

Steady-State Response of a Thermoelastic Half-Space to the Rapid Motion of Surface Thermal/Mechanical Loads

Shear and normal tractions and a heat flux are applied to a largely arbitrary area that moves with constant subsonic speed over a half-space surface. The half-space is a coupled thermoelastic solid of the Jeffreys type, so that the governing steady-state equations involve three thermoelastic characteristic lengths and a dimensionless coupling constant. This constant and one of the lengths remain in the limit as the solid reduces to the standard coupled thermoelastic material. The problem is solved exactly in an integral transform space, and asymptotic expressions for the normal displacement and the temperature change induced on the half-space surface are extracted. These are in principle valid for large distances from the loading zone as measured along its line of travel but, because the scaling dimension is of O(10^-14)μ m, they are robust. Exact inversions are performed, and the results show marked dependence on both loading zone speed and thermoelastic parameters. Indeed, the role of the latter is enhanced as the speed is increased. Singular behavior is found, in particular, when the loading zone moves with the effective thermoelastic Rayleigh speed, an exact formula for which is also given. This revised version was published online in August 2006 with corrections to the Cover Date.     

Three-dimensional analytical elasticity solution for loaded functionally graded coated half-space

In the present paper, a three-dimensional problem of elasticity of normal and tangential loading of surface of the functionally graded coated half-space is considered. In case when Poisson's ratio is constant and the Young's modulus is a power or exponential function of the distance from the surface of the half-space, analytical solution using Fourier transform is obtained. Stress field due to Hertz contact pressure in an elliptical region are studied as a function of the parameter b/a (where a and b are axes of the contact ellipse) and coating thickness.     

Stokes’ first problem for a second grade fluid in a porous half-space with heated boundary

Based on a modified Darcy's law, Stokes’ first problem was investigated for a second grade fluid in a porous half-space with a heated flat plate. Exact solutions of the velocity and temperature fields were obtained using Fourier sine transforms. In contrast to the classical Stokes’ first problem, there is a steady-state solution for the second grade fluid in the porous half-space, which is a damping exponential function with respect to the distance from the flat plate. The well-known solutions for Newtonian fluids in non-porous or porous half-space appear in limiting cases of our solutions.     

Displacement discontinuity method for modeling axisymmetric cracks in an elastic half-space

This paper describes a displacement discontinuity method for modeling axisymmetric cracks in an elastic half-space or full space. The formulation is based on hypersingular integral equations that relate displacement jumps and tractions along the crack. The integral kernels, which represent stress influence functions for ring dislocation dipoles, are derived from available axisymmetric dislocation solutions. The crack is discretized into constant-strength displacement discontinuity elements, where each element represents a slice of a cone. The influence integrals are evaluated using a combination of numerical integration and a recursive procedure that allows for explicit integration of hyper- and Cauchy singularities. The accuracy of the solution at the crack tip is ensured by adding corrective stresses across the tip element. The method is validated by a comparison with analytical and numerical reference solutions.     

Anti-plane shear Green’s functions for an isotropic elastic half-space with a material surface

This paper presents analytical Green’s function solutions for an isotropic elastic half-space subject to anti-plane shear deformation. The boundary of the half-space is modeled as a material surface, for which the Gurtin–Murdoch theory for surface elasticity is employed. By using Fourier cosine transform, analytical solutions for a point force applied both in the interior or on the boundary of the half-space are derived in terms of two particular integrals. Through simple numerical examples, it is shown that the surface elasticity has an important influence on the elastic field in the half-space. The present Green’s functions can be used in boundary element method analysis of more complicated problems.     

Anti-plane shear Lamb''s problem treated by gradient elasticity with surface energy

The consideration of higher-order gradient effects in a classical elastodynamic problem is explored in this paper. The problem is the anti-plane shear analogue of the well-known Lamb's problem. It involves the time-harmonic loading of a half-space by a single concentrated anti-plane shear line force applied on the half-space surface. The classical solution of this problem based on standard linear elasticity was first given by J.D. Achenbach and predicts a logarithmically unbounded displacement at the point of application of the load. The latter formulation involves a Helmholtz equation for the out-of-plane displacement subjected to a traction boundary condition. Here, the generalized continuum theory of gradient elasticity with surface energy leads to a fourth-order PDE under traction and double-traction boundary conditions. This theory assumes a form of the strain-energy density containing, in addition to the standard linear-elasticity terms, strain-gradient and surface-energy terms. The present solution, in some contrast with the classical one, predicts bounded displacements everywhere. This may have important implications for more general contact problems and the Boundary-Integral-Equation Method.     

层状地基任意形状刚性基础动力响应求解

提出了基于积分变换、对偶方程与精细积分算法求解多层地基任意形状刚性基础的动力刚度问题. 首先在频率波数域内圆柱坐标体系中利用圆形微元的对称与反对称特性建立多层地基中格林影响函数的波动方程,然后将应力和位移关系表示成对偶形式进行精细积分求解以提高计算精度和稳定性. 再将任意形状刚性基础与地基的交界面离散化为一系列圆形微元,利用格林影响函数建立其平动与转动动力刚度的矩阵方程. 该求解方法高效、准确并且计算稳定,适于任意复杂多层地基任意形状基础动力刚度的计算.