Three-dimensional elasticity solution for bending of transversely isotropic functionally graded plates

This paper presents a three-dimensional elasticity solution for a simply supported, transversely isotropic functionally graded plate subjected to transverse loading, with Young’s moduli and the shear modulus varying exponentially through the thickness and Poisson’s ratios being constant. The approach makes use of the recently developed displacement functions for inhomogeneous transversely isotropic media. Dependence of stress and displacement fields in the plate on the inhomogeneity ratio, geometry and degree of anisotropy is examined and discussed. The developed three-dimensional solution for transversely isotropic functionally graded plate is validated through comparison with the available three-dimensional solutions for isotropic functionally graded plates, as well as the classical and higher-order plate theories.     

Three-dimensional analytical solution for a rotating disc of functionally graded materials with transverse isotropy

Based on the basic equations for axisymmetric problems of transversely isotropic elastic materials, the displacement components are expressed in terms of polynomials of the radial coordinate with the five involved coefficients, named as displacement functions in this paper, being undetermined functions of the axial (thickness) coordinate. Five equations governing the displacement functions are then derived. It is shown that the displacement functions can be found through progressive integration by incorporating the boundary conditions. Thus a three-dimensional analytical solution is obtained for a transversely isotropic functionally graded disc rotating at a constant angular velocity.The solution can be degenerated into that for an isotropic functionally graded rotating disc. A prominent feature of this solution is that the material properties can be arbitrary functions of the axial coordinate. Thus, the solution for a homogeneous transversely isotropic rotating disc is just a special case that can be easily derived. An example is finally considered for a special functionally graded material, and numerical results shows that the material inhomogeneity has a remarkable effect on the elastic field.     

均布荷载作用下功能梯度简支梁弯曲的解析解

利用应力函数半逆解法,研究了均布载荷作用下、材料属性在厚度上任意变化的功能梯度简支梁弯曲的解析解,给出了各向应力应变与位移的解析显式表达式.首先根据平面应力状态的基本方程,得出了功能梯度梁的应力函数应满足的偏微分方程,并根据应力边界条件得出了各应力分布的表达式;进而根据功能梯度材料的本构方程和位移边界条件,得出了应变和位移的分布.最后,通过将本文的解退化到均质各向同性梁并与经典弹性解比较,证明了本文理论的正确性,并求解了材料组分呈幂律分布的功能梯度梁的应力和位移分布,分析了上下表层材料的弹性模量比λ与组分材料体积分数指数n对应力和位移分布的影响.     

The radially nonhomogeneous elastic axisymmentric problem

The plane axisymmetric problem with axisymmetric geometry and loading is analyzed for a radially nonhomogeneous circular cylinder, in linear elasticity. Considering the radial dependence of the stress, the displacements fields and of the stiffness matrix, after a series of admissible functional manipulations, the general differential system solving the problem is developed. The isotopic radially inhomogeneous elastic axisymmetric problem is also analyzed. The exact elasticity solution is developed for a radially nonhomogeneous hollow circular cylinder of exponential Young’s modulus and constant Poisson’s ratio and of power law Young’s modulus and constant Poisson’s ratio. For the isotropic elastic axisymmetric problem, a general expression of the stress function is derived. After the satisfaction of the biharmonic equation and making compatible the stress field’s expressions, the stress function and the stress and displacements fields of the axisymmetric problem are also deduced. Applications have been made for a radially nonhomogeneous hollow cylinder where the stress and displacements fields are determined.     

POINT FORCE SOLUTION FOR A TRANSVERSELY ISOTROPIC ELASTIC LAYER

ＰＯＩＮＴＦＯＲＣＥＳＯＬＵＴＩＯＮＦＯＲＡＴＲＡＮＳＶＥＲＳＥＬＹＩＳＯＴＲＯＰＩＣＥＬＡＳＴＩＣＬＡＹＥＲＰＯＩＮＴＦＯＲＣＥＳＯＬＵＴＩＯＮＦＯＲＡＴＲＡＮＳＶＥＲＳＥＬＹＩＳＯＴＲＯＰＩＣＥＬＡＳＴＩＣＬＡＹＥＲ￥ＤｉｎｇＨａｏｊｉａｎｇ（丁...     

Finite-part integral and boundary element method to solve three-dimensional crack problems in piezoelectric materials

Using the hypersingular integral equation method based on body force method, a planar crack in a three-dimensional transversely isotropic piezoelectric solid under mechanical and electrical loads is analyzed. This crack problem is reduced to solve a set of hypersingular integral equations. Compare with the crack problems in elastic isotropic materials, it is shown that for the impermeable crack, the intensity factors for piezoelectric materials can be obtained from those for elastic isotropic materials. Based on the exact analytical solution of the singular stresses and electrical displacements near the crack front, the numerical method of the hypersingular integral equation is proposed by the finite-part integral method and boundary element method, which the square root models of the displacement and electric potential discontinuities in elements near the crack front are applied. Finally, the numerical solutions of the stress and electric field intensity factors of some examples are given.     

Three-dimensional elastic solutions for functionally graded circular plates

Three-dimensional elastic solutions are obtained for a functionally graded thick circular plate subject to axisymmetric conditions. We consider a isotropic material where the Young modulus depends exponentially on the position along the thickness, while the Poisson ratio is constant. The solution method utilises a Plevako’s representation form which reduces the problem to the construction of a potential function satisfying a linear fourth-order partial differential equation. We write this potential function in terms of Bessel functions and we pointwise assign mixed boundary conditions. The analytic solution is obtained in a general form and explicitly presented by assuming transversal load on the upper face and zero displacements on the mantle; this is done by superposing the solutions of problems with suitably imposed radial displacement. We validate the solution by means of a finite element approach; in this way, we highlight the effects of the material inhomogeneity and the limits of the employed numerical method near the mantle, where the solution shows a large sensitivity to the boundary conditions.     

Exact Solutions for a Thick Elastic Plate with a Thin Elastic Surface Layer

A procedure has been developed in previous papers for constructing exact solutions of the equations of linear elasticity in a plate (not necessarily thin) of inhomogeneous isotropic linearly elastic material in which the elastic moduli depend in any specified manner on a coordinate normal to the plane of the plate. The essential idea is that any solution of the classical equations for a hypothetical thin plate or laminate (which are two-dimensional theories) generates, by straightforward substitutions, a solution of the three-dimensional elasticity equations for the inhomogeneous material. In this paper we consider a thick plate of isotropic elastic material with a thin surface layer of different isotropic elastic material. It is shown that the interface tractions and in-plane stress discontinuities are determined only by the initial two-dimensional solution, without recourse to the three-dimensional elasticity theory. Two illustrative examples are described.     

Exact Lévy-Type Solutions for Plate Bending Exist for Transversely Isotropic but Not for General Monoclinic Materials

An exact three-dimensional Lévy-type solution for the bending of an elastic slab is one in which there are edge loads only and the unknown displacement and stresses have very simple polynomial dependence on the thickness coordinate. Lévy obtained such a solution in 1877 for a linearly elastic, isotropic, plate-like body. The most general material that allows bending and stretching to be uncoupled is monoclinic (13 elastic constants). However, it is shown that only transversely isotropic materials (5 elastic constants) admit exact solutions having polynomial dependence in the thickness direction. Such solutions are listed explicitly.     

压电介质中受拉伸与弯曲联合作用的圆币形裂纹问题

以弹性位移分量和电势函数为基本未知量时，横观各向同性压电介质非轴对称三维问题的控制微分方程是四个二阶线性偏微分方程相联立的方程组。本文导出了用四个调和函数表示位移及电势的该方程组的势函数通解。作为通解的应用举例，文中求解了压电陶瓷材料中受拉伸与弯曲联合作用的圆币形裂纹问题，得到了裂纹尖端附近应力场及电位移场的解析表达式。结果表明裂尖场以及应力强度因子和电位移强度因子均表现出复杂的机－电耦合行为。     

Analysis method of planar interface cracks of arbitrary shape in three-dimensional transversely isotropic magnetoelectroelastic bimaterials

Using the fundamental solutions for three-dimensional transversely isotropic magnetoelectroelastic bimaterials, the extended displacements at any point for an internal crack parallel to the interface in a magnetoelectroelastic bimaterial are expressed in terms of the extended displacement discontinuities across the crack surfaces. The hyper-singular boundary integral–differential equations of the extended displacement discontinuities are obtained for planar interface cracks of arbitrary shape under impermeable and permeable boundary conditions in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. An analysis method is proposed based on the analogy between the obtained boundary integral–differential equations and those for interface cracks in purely elastic media. The singular indexes and the singular behaviors of near crack-tip fields are studied. Three new extended stress intensity factors at crack tip related to the extended stresses are defined for interface cracks in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. A penny-shaped interface crack in magnetoelectroelastic bimaterials is studied by using the proposed method.The results show that the extended stresses near the border of an impermeable interface crack possess the well-known oscillating singularity r^?1/2±iε or the non-oscillating singularity r^?1/2±κ. Three-dimensional transversely isotropic magnetoelectroelastic bimaterials are categorized into two groups, i.e., ε-group with non-zero value of ε and κ-group with non-zero value of κ. The two indexes ε and κ do not coexist for one bimaterial. However, the extended stresses near the border of a permeable interface crack have only oscillating singularity and depend only on the mechanical loadings.     

Green’s functions of an exponentially graded transversely isotropic half-space

By virtue of a complete set of displacement potential functions and Hankel transform, the analytical expressions of Green’s function of an exponentially graded elastic transversely isotropic half-space is presented. The given solution is analytically in exact agreement with the existing solution for a homogeneous transversely isotropic half-space. Employing a robust asymptotic decomposition technique, the Green’s function is decomposed to the closed-form Green’s function corresponding to the homogeneous transversely isotropic half-space and grading term with strong decaying integrands. This representation is very useful for numerical methods which are based on boundary-integral formulations such as boundary-element method since the numerically evaluated part is not responsible for the singularity. The high accuracy of the proposed numerical scheme is confirmed by some numerical examples.     

Stress concentration factor due to a functionally graded ring around a hole in an isotropic plate

The aim of this work is to present an analytical solution to reduce the stress concentration factor (SCF) around a circular hole in an isotropic homogeneous plate subjected to far-field uniaxial loading. In this paper the elastic response of an inhomogeneous annular ring made of functionally graded material (FGM), inserted around a hole of a homogeneous plate, is studied. By assuming that Young’s modulus varies in the radial direction with power law and that Poisson’s ratio is constant, the governing differential equations for plane stress conditions are obtained. Using stress function a general solution in explicit closed form is presented and the SCF investigated to highlight the inhomogeneity effects. Furthermore, the explicit solution for an inner homogeneous ring, with different properties with respect to those of the plate, is explicitly obtained and numerical results are compared between homogeneous ring and FGM ring.     

横观各向同性材料三维裂纹问题的数值分析

严格从三维横观各向同性材料弹性空间问题的Green函数出发,采用Hadamard有限部积分概念,导出了三维状态下单位位移间断(位错)集度的基本解.在此基础上,将三维任意形状的片状裂纹问题归结为求解-组以未知位移间断表示的超奇异积分方程;并给出了边界元离散形式.对方程中出现的超奇异积分,采用了Had-alnard定义的有限部积分来处理.论文最后给出了若干典型片状裂纹问题的数值算例,数值结果表明了本文方法是非常有效的.     

Finite and transversely isotropic elastic cylinders under compression with end constraint induced by friction

This paper derives an exact solution for the non-uniform stress and displacement fields within a finite, transversely isotropic, and linear elastic cylinder under compression with a kind of radial constraint induced by friction between the end surfaces of the cylinder and the loading platens. The main feature of the present work is the introduction of a general solution form for Lekhnitskii’s stress function such that the governing equation and all end and curved boundary conditions of the cylinder are satisfied exactly. Two different solutions were obtained corresponding to the real or complex characteristic roots of the governing equation, depending on the combination of the elastic material constants. The solution by Watanabe [Watanabe, S., 1996. Elastic analysis of axi-symmetric finite cylinder constrained radial displacement on the loading end. Structural Engineering/Earthquake Engineering JSCE 13, 175s–185s] for isotropic cylinders under compression test can be recovered as a special case. Our numerical results show that both the non-uniform stress distribution and the difference between the apparent and the true Young’s moduli of the cylinder are very sensitive to the anisotropy of Young’s moduli, Poisson’s ratios and shear moduli. A more distinct bulging shape of the cylinder is expected when anisotropy in shear modulus is strong, the cylinder is relatively short, and the end constraint is large. The bulging shape, however, does not depend strongly on anisotropy of either Poisson’s ratio or Young’s modulus.     

胶补加筋板应力强度因子及胶层剪应力计算

 在铆钉力法的基础上，提出计算胶接修补加筋板结构应力强度因子及胶层剪应力的解析法, 对分析胶接修补的加筋板强度具有现实意义. 利用影响系数，可以得到筋条、补片及蒙皮位 移的表达式，根据位移协调条件，建立联立方程，从而得到胶层剪应力及铆钉力. 通过解析 式，计算裂纹尖端的应力强度因子，分析补片对结构修补的效果，综合考虑了筋条及补片对 应力强度因子的影响.     

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.     

Surface effects in the deformation of an anisotropic elastic material with nano-sized elliptical hole

We examine the effect of surface energy on an anisotropic elastic material weakened by an elliptical hole. A closed-form, full-field solution is derived using the standard Stroh formalism. In particular, explicit expressions for the hoop stress, normal, in-plane tangential and out-of-plane displacement components along the edge of the hole are obtained. These expressions clearly demonstrate the effect of elastic anisotropy of the bulk material on the corresponding field variables. When the material becomes isotropic, the hoop stress agrees with the well-known result in the literature while both the in-plane tangential and out-of-plane displacements vanish and the normal displacement is constant along the entire boundary of the elliptical hole.     

Benchmark trigonometric and 3-D elasticity solutions for an exponentially graded thick rectangular plate

The bending problem of a transverse load acting on an isotropic inhomogeneous rectangular plate using both two-dimensional (2-D) trigonometric and three-dimensional (3-D) elasticity solutions is considered. In the present 2-D solution, trigonometric terms are used for the displacements in addition to the initial terms of a power series through the thickness. The effects due to transverse shear and normal deformations are both included. The form of the assumed 2-D displacements is simplified by enforcing traction-free boundary conditions at the faces of the plate. No transverse shear correction factors are needed because a correct representation of the transverse shearing strain is given. The plate material is exponentially graded, meaning that Lamé’s coefficients vary exponentially in a given fixed direction (the thickness direction). A wide variety of results for the displacements and stresses of an exponentially graded rectangular plate are presented. The validity of the present 2-D trigonometric solution is demonstrated by comparison with the 3-D elasticity solution. The influence of aspect ratio, side-to-thickness ratio and the exponentially graded parameter on the bending response are investigated.     

绕径旋转横观各向同性圆球壳的变形和应力分析

本文用三维弹性力学理论研究绕径旋转横观各向同性圆球壳的变形和应力分布，给出任意厚度圆球壳的位移和应务的封闭形式解以及数值结果。研究表明球壳上的位移和应力大小与旋转角速度的平方成正比，球面内的正应力远大于径向正应力和剪应力。