任意梯度分布功能梯度圆环的热弹性分析

对沿径向任意变化的材料参数的功能梯度圆环进行了热弹性分析.与以前关于该问题的分析不同,既不需要预先给定具体的梯度变化形式,也不需要对结构进行细分.给出一种新的有效解法将问题转换为求解Fredholm积分方程,从而通过Fredholm积分方程的解给出热应力和位移的分布情况.最后通过算例分析了内外表面受不同温度作用时,材料参数呈现梯度变化对圆环的应力和位移变化的影响,计算结果表明某些特定的材料梯度可有效缓解圆环内的热应力分布.该文得到的结果对功能梯度圆环在结构安全设计方面有重要的理论指导意义.     

功能梯度圆板的轴对称自由振动

基于三维弹性理论,用直接位移法求解了横观各向同性功能梯度圆板的轴对称自由振动,其材料特性沿板厚度方向按指数形式变化,在弹性简支和刚性滑动两种边界条件下得到了三维精确解,即它逐点满足基本方程和边界条件,最后给出了数值算例并与以前的工作进行了对比.该方法也可推广应用于材料特性沿板厚度任意变化情形.     

功能梯度厚壁中空圆柱体弹性动力学平面应变响应的解析解

研究了边界表面受均布动压力作用的功能梯度(FGM)厚壁中空圆柱体,给出了其平面应变响应下的弹性动力学解．假设材料性能(除Poisson比外)随厚度按幂律函数变化．为了得到一个精确解,将动力径向位移分为准静力部分和动力部分,导出了每个部分的一个解析解．先由Euler方程得到准静力学部分的解,再由分离变量法和正交展开法得到动力学部分的解．在不同动荷载作用下,对不同的FGM中空圆柱体,画出径向位移和应力图,并对本方法的优点进行了讨论．该解析解适用于中空圆柱体各种组合的FGM,厚度可以是任意的,初始条件也可以是任意的,壁面上均匀分布着任意形式的动压力．     

功能梯度板的柱面弯曲弹性力学解

在推广后的England-Spencer功能梯度板理论基础上,研究了功能梯度板在不同荷载作用下的柱面弯曲问题.采用该理论中的位移展开公式,并且材料参数沿板厚方向可以任意连续变化,并将材料由各向同性推广到正交各向异性.假设板在y方向无限长,最终建立了一个从弹性力学理论出发的正交各向异性功能梯度板在横向分布荷载作用下柱面弯曲问题的板理论.通过算例分析,讨论了边界条件、材料梯度及板厚跨比等因素对功能梯度板静力响应的影响.     

正交各向异性功能梯度层合板稳态热传导问题的精确解

对轴对称正交各向异性功能梯度层合圆板稳态热传导问题进行精确分析.假设材料热传导率沿板厚方向按指数函数形式梯度分布,从正交各向异性功能梯度圆板稳态热传导的基本方程出发,利用分离变量法,获得了在上、下表面作用任意热分布情况下的精确解.通过数值算例的分析,指出材料性质的梯度变化、板厚边界条件等分析了对温度场分布的影响.所获得的精确结果,可以作为评价其它近似方法的标准解答.     

功能梯度板柱面弯曲的弹性力学解

利用推广后的Main和Spencer功能梯度板理论,研究了功能梯度矩形板在均布荷载作用下的柱面弯曲问题.采用该理论中的位移展开公式,并且材料参数沿板厚方向可以任意连续变化,但将材料由各向同性推广到正交各向异性,以及由不考虑板的横向荷载作用发展到受横向均布荷载作用.假设板在y方向无限长,从而得到了一个从弹性力学理论出发的正交各向异性功能梯度板在横向均布荷载作用下柱面弯曲问题的板理论.通过算例分析,讨论了边界条件和梯度变化程度对功能梯度板静力响应的影响.     

变厚度非均匀粘弹性复合圆柱体的旋转应力

在平面应变的假设下,给出了两个复合弹性圆柱体旋转时的解析解．外柱是由厚度按公式变化的正交各向异性材料所组成,它包裹着一个等厚度纤维增强粘弹性均匀各向同性的实心圆柱体．外圆柱体的厚度和弹性性质按半径方向的幂函数变化．应用边界和连续条件,确定复合圆柱体旋转时的径向位移和应力,应用等效模量和Illyushin逼近法,得到问题的粘弹性解．讨论了各向异性、厚度变化、本构参数以及时间参数,对径向位移和应力的影响．     

具有功能梯度加强环的有限尺寸开孔板应力集中问题

基于复变函数理论,结合最小二乘边界配点法,对具有功能梯度加强环的有限尺寸开孔板在任意均布载荷作用下的应力集中问题进行了研究.首先,采用分层均匀化方法,给出了材料参数沿径向任意变化的功能梯度加强环内的复势及孔边应力的半解析解;然后,通过几组数值算例,讨论了组分梯度、加强环厚度、板相对尺寸及偏心率的变化对孔边应力集中的影响.结果表明,通过合理选择功能梯度加强环内材料参数的递变规律及加强环的厚度,可以有效缓解有限尺寸开孔板内的应力集中.     

饱和地基上弹性基础的竖向振动特性研究

采用解析的方法研究了饱和地基上受一简谐竖向荷载作用下弹性基础的动力响应．在分析中,首先利用积分变换技术获得了饱和介质基本控制方程的变换解,然后基于基础-半空间完全放松接触、半空间表面完全透水或不透水的假设,建立了该动力混合边值问题的对偶积分方程,并把该对偶积分方程进一步化为易于数值求解的第二类Fredholm积分方程A·D2文末数值算例给出了动力柔度系数、位移和孔隙水压力随振动频域和土-基础体系物理力学参数特性的变化曲线．结果表明:饱和地基上弹性基础的动力响应完全不同于饱和地基上刚性圆板的动力响应．所用方法可用于研究波的传播、土-结构动力相互作用等许多问题．     

面内功能梯度三角形板等几何面内振动分析

基于平面应变理论,利用等几何有限元方法分析了弹性边界条件下面内功能梯度三角形板的面内振动特性.板的材料属性沿厚度方向呈均匀分布,而在面内方向呈任意指数梯度变化.采用非均匀有理B样条(NURBS)基函数对三角形结构进行等几何建模和位移描述,实现了三角形板几何设计和振动分析的无缝衔接.在三角形板边界上引入虚拟弹簧约束并通过调节虚拟弹簧刚度,实现任意边界条件的施加.通过不同的单元细化方案和对比算例,验证了等几何方法的灵活性、准确性和快速收敛性.系统研究了边界条件、材料属性和几何参数对三角形板振动特性的影响.同时给出了弹性边界条件下面内功能梯度三角形板的振动特性解,具有重要参考价值.     

Rotation Stability of Anisotropic Disks

The critical state caused by the loss of stability of coils in rotating anisotropic disks is investigated by a numerical method. The investigation is performed for different ratios between the circumferential and radial elastic moduli and different relative thicknesses of the disks, Poisson ratios, and the ratios of radial and shear moduli. It is shown that, if a disk is made of materials compliant in the radial direction, the coils can lose their stability and shift according to the first harmonic (with displacement of the center of mass of the disk) or the second and third harmonics of vibrations (without displacement of the center of mass) within the range of its rotational velocities. These vibrations can be real for some materials, e.g., for unidirectional composites with a compliant (polyurethane) binder, and can precede the failure of disks from the action of circumferential stresses.     

Concerning an approach to identifying the Lamé parameters of an elastic functionally graded cylinder

Based on the general linear elasticity relations, an axisymmetric problem on the steady-state oscillations of a functionally graded hollow cylinder is formulated. The Lamé parameters are considered variable in radial coordinate. Oscillations are caused by the distributed load applied to the outer part of the cylinder boundary. Using the variable separation method, the direct problem on determining the radial and longitudinal components of the displacement field is investigated. The influence of the laws of variation for the Lamé parameters on acoustic characteristics is analysed. The inverse coefficient problem on the identification of the variable Lamé parameters from the data on the amplitude-frequency characteristic is stated. Based on the weak formulation of the problem for an elastic inhomogeneous body, a general linearised relation for the desired and given characteristics is obtained. A system of the Fredholm integral equations of the first kind is formulated with respect to two unknown corrections to the restored laws of the Lamé parameters change. The solution is built by means of an iterative process. A reconstruction of various laws of changing the Lamé parameters is carried out. The accuracy of the presented algorithm is estimated, and recommendations for the most efficient implementation of the reconstruction procedure are proposed.     

Investigation on a rotating FGPM circular disk under a coupled hygrothermal field

In this paper, the distributions of the temperature, moisture, displacement and stress of a functionally graded piezoelectric material (FGPM) circular disk rotating around its axis at a constant angular velocity under a coupled hygrothermal field are presented by a numerical method. The material properties of the FGPM circular disk are assumed to vary along the radial coordinate exponentially. First, the coupled hygrothermal field along the radius of a rotating circular disk is achieved by solving the coupled hygrothermal equations, and then the dynamic equilibrium is solved by utilizing the finite difference method. Finally, numerical results show the effects of functionally graded index, inner radius, angular speed and hygrothermal index on the hygrothermal behaviors of the FGPM circular disk. The results can be useful for the optimal design of rotating FGPM circular disks under a coupled hygrothermal field.     

An axisymmetric problem of an embedded mixed-mode crack in a functionally graded magnetoelectroelastic infinite medium

This paper investigates the problem of an axisymmetric penny shaped crack embedded in an infinite functionally graded magneto electro elastic medium. The loading consists of magnetoelectromechanical loads applied on the crack surfaces assumed to be magneto electrically impermeable. The material’s gradient is parallel to the axisymmetric direction and is perpendicular to the crack plane. An anisotropic constitutive law is adopted to model the material behavior. The governing equations are converted analytically using Hankel transform into coupled singular integral equations, which are solved numerically to yield the crack tip stress, electric displacement and magnetic induction intensity factors. A similar problem but with a different crack morphology, that is a plane crack embedded in an infinite functionally graded magneto electro elastic medium, was considered by the authors in a previous work (Rekik et al., 2012) [25]. While the overall solution schemes look similar, the axisymmetric problem resulted in more mathematical complexities and let to different conclusions with respect to the influence of coupling between elastic, electric and magnetic effects. The main focus of this paper is to study the effect of material non-homogeneity on the fields’ intensity factors to understand further the behavior of graded magnetoelectroelastic materials containing penny shaped cracks and to inspect the effect of varying the crack geometry.     

Two-dimensional adhesion mechanics of a graded coated substrate under a rigid cylindrical punch based on a PWEML model

The paper aims at the contact mechanics of functionally graded coated substrate by taking into the adhesion effect. The coating-substrate structure is indented by a cylindrical punch to form a contact region where the adhesion forces are described by using the Maugis adhesion model. A piece-wise exponential multi-layered (PWEML) model is used to simulate the functionally graded materials with arbitrary spatial variation of material properties. This model divided the functionally graded coating into several sub-layers in which the elastic parameter varies as exponential form. Using the Fourier transform technique and the Transfer matrix method, the boundary value problem for adhesive contact of graded coated substrate is reduced to the singular integral equation. Some numerical results are presented to analyze the influence of gradient index on the pull-out force, contact stresses and adhesion region. The results can be applied to improve the performance of the coating by adjusting the gradient index.     

Boussinesq indentation of a transversely isotropic half-space reinforced by a buried inextensible membrane

The normal indentation of a rigid circular disk into the surface of a transversely isotropic half-space reinforced by a buried inextensible thin film is addressed. By virtue of a displacement potential function and the Hankel transform, the governing equations of this axisymmetric mixed boundary value problem are represented as a dual integral equation, which is subsequently reduced to a Fredholm integral equation of the second kind. Two important results of the contact stress distribution beneath the disk region as well as the equivalent stiffness of the system are expressed in terms of the solution of the Fredholm integral equation. When the membrane is located on the surface or at the remote boundary, exact closed-form solutions are presented. For the limiting case of an isotropic half-space the results are verified with those available in the literature. As a special case, the elastic fields of a reinforced transversely isotropic half-space under the action of surface axisymmetric patch loads are also given. The effects of anisotropy, embedment depth of the membrane, and material incompressibility on both the contact stress and the normal stiffness factor are depicted in some plots.     

An interface crack moving between magnetoelectroelastic and functionally graded elastic layers

A constant crack moving along the interface of magnetoelectroelastic and functionally graded elastic layers under anti-plane shear and in-plane electric and magnetic loading is investigated by the integral transform method. Fourier transforms are applied to reduce the mixed boundary value problem of the crack to dual integral equations, which are expressed in terms of Fredholm integral equations of the second kind. The singular stress, electric displacement and magnetic induction near the crack tip are obtained asymptotically and the corresponding field intensity factors are defined. Numerical results show that the stress intensity factors are influenced by the crack moving velocity, the material properties, the functionally graded parameter and the geometric size ratios. The propagation of the moving crack may bring about crack kinking, depending on the crack moving velocity and the material properties across the interface.     

Moving contact problems involving a rigid punch and a functionally graded coating

Frictional contact mechanics analysis for a rigid moving punch of an arbitrary profile and a functionally graded coating/homogeneous substrate system is carried out. The rigid punch slides over the coating at a constant subsonic speed. Smooth variation of the shear modulus of the graded coating is defined by an exponential function and the variation of the Poisson's ratio is assumed negligible. Coulomb's friction law is adopted. Hence, tangential force is proportional to the normal applied force through the coefficient of friction. An analytical method is developed utilizing the singular integral equation approach. Governing partial differential equations are derived in accordance with the theory of elastodynamics. The mixed boundary value problem is reduced to a singular integral equation of the second kind, which is solved numerically by an expansion-collocation technique. Presented results illustrate the effects of punch speed, coefficient of friction, material inhomogeneity and coating thickness on contact stress distributions and stress intensity factors. Comparisons indicate that the difference between elastodynamic and elastostatic solutions tends to be quite larger especially at higher punch speeds. It is shown that use of the elastodynamic theory provides more realistic results in contact problems involving a moving punch.