Torsion of rigid circular shaft of varying diameter embedded in an elastic half space

The axially symmetric torsion of rigid circular shaft of varying diameter embedded in an elastic half space is studied by line-loaded integral equation method (LLIEM), where the problem is formulated by distributions of ficitious fundamental loads PRCHS (point ring couple in half space) along the axis of symmetry in interval of the shaft and is reduced to a one-dimensional and non-singular Fredholm integral equation of the first kind and is easily solved numerically. Numerical examples of torsin of rigid conic, cylinder, conical-cylinder embedded in an elastic half space are given and compared with the known result obtained by the others. The exact solution of torsion of rigid half sphere embedded in an elastic half space is also presented. Project Supported by the National Science Foundation of China.     

Scattering of SH-wave from interface cylindrical elastic inclusion with a semicircular disconnected curve

Scattering of SH wave from an interface cylindrical elastic inclusion with a semicircular disconnected curve is investigated. The solution of dynamic stress concentration factor is given using the Green＇s function and the method of complex variable functions. First, the space is divided into upper and lower parts along the interface. In the lower half space, a suitable Green＇s function for the problem is constructed. It is an essential solution of the displacement field for an elastic half space with a semi-cylindrical hill of cylindrical elastic inclusion while bearing out-plane harmonic line source load at the horizontal surface. Thus, the semicircular disconnected curve can be constructed when the two parts are bonded and continuous on the interface loading the undetermined anti-plane forces on the horizontal surfaces. Also, the expressions of displacement and stress fields are obtained in this situation. Finally, examples and results of dynamic stress concentration factor are given. Influences of the cylindrical inclusion and the difference parameters of the two mediators are discussed.     

Scattering of SH-wave from interface cylindrical elastic inclusion with a semicircular disconnected curve

Scattering of SH wave from an interface cylindrical elastic inclusion with a semicircular disconnected curve is investigated.The solution of dynamic stress concentration factor is given using the Gteen's function and the method of complex variable functions.First,the space is divided into upper and lower parts along the interface.In the lower half space,a suitable Green's function for the problem is constructed.It is an essential solution of the displacement field for an elastic half space with a semi-cylindrical hill of cylindrical elastic inclusion while bearing out-plane harmonic line source load at the horizontal surface.Thus,the semicircular disconnected curve can be constructed when the two parts are bonded and continuous on the interface loading the undetermined anti-plane forces on the horizontal surfaces.Also,the expressions of displacement and stress fields are obtained in this situation.Finally,examples and results of dynamic stress concentration factor are given.Influences of the cylindrical inclusion and the difierence parameters of the two mediators are discussed.     

Dynamics of a prestressed stiff layer on an elastic half space: filtering and band gap characteristics of periodic structural models derived from long-wave asymptotics

Anti-plane and plane-strain, time-harmonic, small-amplitude vibrations of an elastic layer on an elastic half space are considered, superimposed upon a state of finite, uniform stress and strain. A (compressible) elastic material is considered, orthotropic with orthotropy axes aligned parallel and orthogonal both to the layer and the prestress principal directions. A non-uniform mass density is assumed in the layer. A formal long-wave asymptotic solution is derived under the assumptions of high contrast between the stiffnesses of the layer and the half space and between certain prestress components and the current elastic shear modulus.It is shown that (i) the layer asymptotically behaves as a beam subject to transversal and axial vibrations; (ii) the response of the half space can be found in a closed-form, under the assumption of plane wave motion (which becomes consistent when the density of the layer is uniform), otherwise it is represented by a hypersingular integral equation; (iii) if the nonlocality introduced by the hypersingular integral equation is restricted to an influence area of finite extent, the integral can be analytically approximated, so that a Winkler-type spring model representing the half space is rigorously derived. For uniform density of the layer, the constants defining the spring model are given as functions of the prestress and anisotropy parameters of the half space; and, finally, (iv) the asymptotic solution provides new analytical expressions for incremental displacement of the layer, which, compared to the exact numerical solution (also included), are shown to perform quite well, even for values of parameters much beyond the limits imposed by the asymptotic analysis.The asymptotic analysis allows us to explore, for the first time, dynamic properties of a periodic layer bonded to an elastic half space and subject to a uniform prestress state. We find that the system exhibits band gaps (ranges of forbidden frequencies) and that the prestress can be used as a parameter tuning the filtering properties of the structure, an effect which may have important consequences in the design of resonant devices.     

Stresses in a half space due to Newtonian gravitation

An efficient general solution is obtained for the problem of the elastic half space z > 0 with a traction-free surface experiencing gravitational attraction to an arbitrarily shaped body located in z < 0. Many components of the stress field can be written down immediately if the potential of the attracting body is known. Results are given for the case of attraction to a uniform sphere.     

Scattering of SH-wave by cracks originating at an elliptic hole and dynamic stress intensity factor

The method of complex function and the method of Green‘s function are used to investigate the problem of SH-wave scattering by radial cracks of any limited length along the radius originating at the boundary of an elliptical hole, and the solution of dynamicstress intensity factor at the crack tip was given. A Green‘s function was constructed for the problem, which is a basic solution of displacement field for an elastic half space containing a half elliptical gap impacted by anti-plane harmonic linear source force at any point of its horizontal boundary. With division of a crack technique, a series of integral equations can be established on the conditions of continuity and the solution of dynamic stress intensity factor can be obtained. The influence of an elliptical hole on the dynamic stress intensity factor at the crack tip was discussed.     

Analysis of sloping elastic pile under arbitrary loads by line-loaded integral equation method

ＩｎｔｒｏｄｕｃｔｉｏｎＳｌｏｐｉｎｇｐｉｌｅｉｓｏｆｔｅｎｕｓｅｄｉｎｆｏｕｎｄａｔｉｏｎｓｏｆｐｉｅｒ,ｂｒｉｄｇｅ,ｏｆｓｈｏｒｅｓｔｒｕｃｔｕｒｅ,ｓｔｒｕｃｔｕｒｅａｎｄｂｕｉｌｄｉｎｇ．Ｔｈｅｓｔｕｄｙｈａｓｂｅｎｃｏｎｄｕｃｔｅｄｂｙｔ...     

Coupling of elastic and thermal deformations I

Flow, Turbulence and Combustion - The general solution for the stress and temperature distribution in an isotropic thermoelastic solid occupying the half space has been given when the classical...     

Scattering of SH-waves by an interacting interface linear crack and a circular cavity near bimaterial interface

An analytical method is developed for scattering of SH-waves and dynamic stress concentration by an interacting interface crack and a circular cavity near bimaterial interface. A suitable Green‘s function is contructed, which is the fundamental solution of the displacement field for an elastic half space with a circular cavity impacted by an out-plane harmonic line source loading at the horizontal surface. First, the bimaterial media is divided into two parts along the horizontal interface, one is an elastic half space with a circular cavity and the other is a complete half space. Then the problem is solved according to the procedure of combination and by the Green‘s function method. The horizontal surfaces of the two half spaces are loaded with undetermined anti-plane forces in order to satisfy continuity conditions at the linking section, or with some forces to recover cracks by means of crack-division technique. A series of Fredholm integral equations of first kind for determining the unknown forces can be set up through continuity conditions as expressed in terms of the Green‘s function. Moreover, some expressions are given in this paper, such as dynamic stress intensity factor (DSIF) at the tip of the interface crack and dynamic stress concentration factor (DSCF) around the circular cavity edge. Numerical examples are provided to show the influences of the wave numbers, the geometrical location of the interface crack and the circular cavity, and parameter combinations of different media upon DSIF and DSCF.     

爆炸荷载作用下饱和半空间中衬砌隧道弹性动力响应IBIEM模拟

基于Biot两相介质理论,采用一种高精度间接边界积分方程法（IBIEM）研究了饱和半空间中浅埋衬砌隧道在内部爆炸荷载作用下的瞬态弹性动力反应。通过典型算例,给出了爆炸荷载作用下隧道附近地表位移、衬砌动应力、围岩径向位移和衬砌表面孔隙水压的时程响应,并对比分析了饱和半空间和全空间中隧道动力响应的区别。研究表明：覆土层厚度对浅埋隧道-围岩整体动力响应特征具有明显影响;衬砌表面透水状态对爆炸荷载的时程响应的影响不显著;随半空间饱和介质孔隙率增加,围岩受隧道内部爆炸影响程度降低,衬砌承担的爆炸作用增大;当和直达波、衬砌内部反射波的峰值叠加作用时,半空间表面反射波对衬砌隧道拱顶附近响应影响显著,使得衬砌动应力幅值、径向位移相比深埋情况大幅度增加。     

Dynamic response of a pile embedded in elastic half space subjected to harmonic vertical loading

An analytical method is developed to investigate the dynamic response of a pile subjected to harmonic vertical loading.The pile is modeled as a one-dimensional(1D)elastic rod.The elastic soil is divided into a homogeneous half space underlying the base of pile and a series of infinitesimally thin layers along the vertical shaft of pile.The analytical solution for the soil-pile dynamic interaction problem is obtained by the method of Hankel transformation.The proposed solution is compared with the classical plane strain solution.Arithmetical examples are presented to demonstrate the sensitivity of the vertical impedance of the pile to relevant parameters.     

The method of solving axisymmetric problems in elastic space by complex function and some examples

This paper proves Love’s stress function of space axisymmetric problem can be represented by choosing two generalized analytic functions of complex variates reasonably^[1], and deduces the expressions of the components of stress displacements and boundary conditions in complex function. To present the feasibility of the method here and examining the truth of the formulae founded in this paper, the problem of circular shaft with globular cavity pressed on the side and pulled at the ends is solved by using power series and the result is the same as that solved by other methods. In the end, the problem of a cone sheared by uniform shear stress on the sideface is solved, and the solution of a cone acted on by gravity is given by converting constant body forces into surface forces.     

Analytical investigation of steady convection of a near-critical Van der Waals gas in a porous thin annular cylinder embedded in a heat-conducting medium

Natural convection of a near-critical Van der Waals gas in a horizontal porous thin annular cylinder embedded in a heat-conducting space with a temperature gradient given at infinity is considered. The two-dimensional problem in a plane orthogonal to the cylinder axis is investigated using analytical methods, accurate to quantities of the order of the annulus thickness. A criterion of the onset of an annulus-section-average flow is obtained. The critical Rayleigh-Darcy number is determined for the general case in which the physical properties of the gas can considerably vary throughout the medium. Several limiting cases are considered and the ranges of their applicability are discussed. It is shown that, as the thermodynamic critical point is approached, the asymptotics of the critical Rayleigh-Darcy number depend on the relation between non-Boussinesq parameters, such as hydrostatic compressibility criteria, the temperature difference, and the nearness to the critical point. In the case of steady convection, an analytical solution is also derived in the case in which the above-mentioned stability threshold is exceeded and the physical properties of the gas vary throughout space only slightly. A comparison with the case of a perfect gas is made.     

Determination of stress intensity factor with direct stress approach using finite element analysis

In this article,a direct stress approach based on finite element analysis to determine the stress intensity fac-tor is improved.Firstly,by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack,we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity fac-tor,which is nearly equal to the stress intensity factor given by the asymptotic solution.Secondly,the crack problem is solved numerically by the finite element method.Depending on the modeling capability of the software,we designed an adaptive mesh model to simulate the stress singularity.Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corre-sponding crack problem.Therefore,the stress intensity factor may be calculated from the stress distribution in the appro-priate interval,with a high accuracy.     

SH波对双相介质界面附近圆形孔洞的散射

建立了求解平面SH波对双相介质界面附近圆形孔洞散射与动应力集中的一种分析方法．利用复变函数与多极坐标的方法构造了一个Green函数，它是在含有圆形孔洞的弹性半空间的水平面上任一点上作用时间谐和的出平面线源荷载的位移解．利用“契合”模型，并根据界面上位移连续性条件，建立了求解SH波对双相介质界面附近圆形孔洞散射的具有弱奇异性的第一类Fredholm型积分方程．给出了圆孔周边上动应力集中系数的表达式．作为算例，分析了在界面一侧或界面两侧附近具有圆形孔洞时SH波的散射，并讨论了入射波波数、不同的材料组合以及孔心至界面的距离对动应力集中的影响．     

Dynamic stress concentration in an elastic half space with a semi-circular cavity excited by SH waves

Subject of the investigation is the stress distribution and the dynamic stress concentration factor at the surface of a semi-circular cavity in a half space excited by plane harmonic SH waves. Using wave function expansion for the incident wave and the reflected waves, a closed form solution is obtained. Numerical results are represented graphically.     

Analysis of three-dimensional edge cracks under tensile loading

Three-dimensional edge cracks are analyzed using the Self-Similar Crack Expansion (SSCE) method with a boundary integral equation technique. The boundary integral equations for surface cracks in a half space are presented based on a half space Green's function (Mindlin, 1936). By using the SSCE method, the stress intensity factors are determined by the crack-opening displacement over the crack surface. In discrete boundary integral equations, the regular and singular integrals on the crack surface elements are evaluated by an analytical method, and the closed form expressions of the integrals are given for subsurface cracks and edge crakcs. This globally numerical and locally analytical method improves the solution accuracy and computational effort. Numerical results for edge cracks under tensile loading with various geometries, such as rectangular cracks, elliptical cracks, and semi-circular cracks, are presented using the SSCE method. Results for stress intensity factors of those surface breaking cracks are in good agreement with other numerical and analytical solutions.     

Bounds for shakedown of cohesive-frictional materials under moving surface loads

In this paper, shakedown of a cohesive-frictional half space subjected to moving surface loads is investigated using Melan’s static shakedown theorem. The material in the half space is modelled as a Mohr–Coulomb medium. The sliding and rolling contact between a roller and the half space is assumed to be plane strain and can be approximated by a trapezoidal as well as a Hertzian load distribution. A closed form solution to the elastic stress field for the trapezoidal contact is derived, and is then used for the shakedown analysis. It is demonstrated that, by relaxing either the equilibrium or the yield constraints (or both) on the residual stress field, the shakedown analysis leads to various bounds for the elastic shakedown limit. The differences among the various shakedown load factors are quantitatively compared, and the influence of both Hertzian and trapezoidal contacts for the half space under moving surface loads is studied. The various bounds and shakedown limits obtained in the paper serve as useful benchmarks for future numerical shakedown analysis, and also provide a valuable reference for the safe design of pavements.     

SH波对内含裂纹衬砌结构的散射及动应力集中

当衬砌结构内含裂纹时 ,采用Green函数的方法 ,研究了SH波对裂纹的散射及其动应力集中 ,构造了在含有半圆形衬砌的弹性半空间上 ,在水平面上任一点承受时间谐和出平面线源载荷作用时的位移函数作为Green函数 ;推导了SH波对衬砌内有裂纹的散射定解积分方程组 ,进而求得裂纹尖端的动应力因子 ,重点讨论了衬砌及周围介质对裂纹尖端动应力因子的影响 ,给出了介质参数变化对裂纹尖端动应力因子的影响曲线 ,为工程设计提供了依据。     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.