Dynamic response for functionally graded materials with penny-shaped cracks

This paper provides a method for studying the penny-shaped cracks configuration in functionally graded material(FGM) structures subjected to dynamic or steady loading. It is assumed that the FGMs are transversely isotropic and all the material properties only depend on the axial coordinatez. In the analysis, the elastic region is treated as a number of layers. The material properties are taken to be constants for each layer. By utilizing the Laplace transform and Hankel transform technique, the general solutions for the layers are derived. The dual integral equations are then obtained by introducing the mechanical boundary and layer interface conditions via the flexibility/stiffness matrix approach. The stress intensity factors are computed by solving dual integral equations numerically in Laplace transform domain. The solution in time domain is obtained by utilizing numerical Laplace inverse. The main advantage of the present model is its ability for treating multiple crack configurations in FGMs with arbitrarily distributed and continuously varied material properties by dividing the FGMs into a number of layers with the properties of each layer slightly different from one another. This work was supported by Failure Mechanics Laboratory of State Education Commission and the Post-doctor Research Fund of China.     

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题．假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向（与裂纹面垂直的方向）为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换祛,在Laplace域内推导出了控制问题的对偶积分方程组．利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率．研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响．     

Generalized thermoelastic functionally graded solid with a periodically varying heat source

This paper deals with the problem of thermoelastic interactions in a functionally graded isotropic unbounded medium due to the presence of periodically varying heat sources in the context of the linear theory of generalized thermoelasticity without energy dissipation (TEWOED). The governing equations of generalized thermoelasticity without energy dissipation (GN model type II) for a functionally graded materials (FGM) (i.e. material with spatially varying material properties)are established. The governing equations are expressed in Laplace–Fourier double transform domain and solved in that domain. Now, the inversion of the Fourier transform is carried out by using residual calculus, where poles of the integrand is obtained numerically in complex domain by using Laguerre’s method and the inversion of Laplace transform is done numerically using a method based on Fourier series expansion technique. The numerical estimates of the displacement, temperature, stress and strain are obtained for a hypothetical material. The solution to the analogous problem for homogeneous isotropic material is obtained by taking nonhomogeneity parameter suitably. Finally the results obtained are presented graphically to show the effect of nonhomogeneity on displacement, temperature, stress and strain.     

Modeling method for a crack problem of functionally graded materials with arbitrary properties—piecewise-exponential model

A new model, piecewise-exponential model (PE model), is developed to investigate the crack problem of the functionally graded materials (FGMs) with arbitrary properties. In the PE model, the functionally graded material is divided into some nonhomogeneous layers along the gradient direction of the properties, with each layer’s properties varying exponentially. By this way, the real material properties can be approached by a series of exponential functions. Since the real material properties are used on both surfaces of each nonhomogeneous layer, the nature of continuously varying properties of FGMs can be approached accurately. The influences of the local nonhomogeneity on the crack-tip fields can be fully considered. By using the new model, the fracture problem of a functionally graded strip with arbitrary properties and a crack vertical to the free surfaces is studied. The integral transform method, the theory of residues and the theory of singular integral equation are applied. Some representative samples with different kinds of nonhomogeneous properties are analyzed and the corresponding stress intensity factors (SIFs) are presented. It is shown that the PE mode is effective for investigating the crack problems of the FGMs with arbitrary properties.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

轴对称任意梯度圆筒的弹性动力学解

基于轴对称平面应变问题的运动方程及弹性梯度材料的应力和位移关系,通过将圆筒分层使材料性质离散为分段常数函数,同时在时域内应用有限差分格式,求得了材料性质沿径向梯度变化的圆筒弹性动力学解。本文解不仅适合任意梯度的弹性圆筒,而且容易满足多种形式的初始条件和边界条件。通过对材料性质沿径向为连续函数分布和分段函数分布的梯度圆筒数值分析,并与已有文献结果比较,得出本文解与已有文献的解吻合较好,验证了本文解的正确性和有效性。对材料性质为分段函数的三层组合圆筒分析发现,中间功能梯度层的指数分布因子对圆筒的径向位移和应力随时间变化都会产生显著影响。     

Hybrid graded element model for transient heat conduction in functionally graded materials

This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest’s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.     

Non-Fourier heat conduction in an exponentially graded slab

The present article investigates one-dimensional non-Fourier heat conduction in a functionally graded material by using the differential transformation method. The studied geometry is a finite functionally graded slab, which is initially at a uniform temperature and suddenly experiences a temperature rise at one side, while the other side is kept insulated. A general non-Fourier heat transfer equation related to the functionally graded slab is derived. The problem is solved in the Laplace domain analytically, and the final results in the time domain are obtained by using numerical inversion of the Laplace transform. The obtained results are compared with the exact solution to verify the accuracy of the proposed method, which shows excellent agreement.     

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.     

含功能梯度材料加强环的任意几何形状孔附近应力集中分析

基于复变函数理论,结合保角变换技术研究含功能梯度材料（FGM）加强环的任意几何形状孔附近应力集中。采用分层均匀化方法,给出了远场均布载荷作用下材料参数沿孔周法线方向任意变化的FGM加强环内的复势函数解。通过数值算例,详细讨论了加强环内杨氏模量不同变化规律对三角形、正方形、矩形等各种几何形状孔附近应力分布的影响。结果表明：通过在孔周衬入FGM加强环并合理选择加强环内材料参数的递变规律,可以有效缓解各种几何形状孔附近的应力集中。同时通过一些特例与已有文献比对验证了本文结果的正确性。     

Dynamic investigation of a functionally graded layered structure with a crack crossing the interface

The dynamic response of a functionally graded layered structure with a crack crossing the interface is analyzed. The in-plane impact loading condition is considered. By using the Laplace and Fourier integral transforms, singular integral equation method and residue theory, the present problem is reduced to a singular integral equation in the Laplace transform domain. The influences of Young’s modulus ratio, thickness ratio, and crack length and location on the dynamic stress intensity factors (DSIFs) are investigated. Particularly, the DSIFs corresponding to different crack locations are shown in the case when the crack center moves from one layer to another layer through the interface. The peak and static values and overshoot characteristics of the DSIFs are analyzed. It is found that these values typically exhibit kinking behavior when the crack tips arrive at the interface. This study is different from previous other investigations in the following respects: (1) the dynamic response of a crack crossing the interface of a functionally graded structure is studied analytically, which has hardly been done in the past and (2) the present model can be reduced to some important problems, such as a functionally graded coating-substrate structure with a crack in the graded coating or homogeneous substrate or one intersecting the interface.     

The Anti-plane Dynamic Fracture Analysis of Functionally Graded Materials with Arbitrary Spatial Variations of Material Properties

Anti-plane dynamic fracture analysis is presented for functionally graded materials (FGM) with arbitrary spatial variations of material properties. The FGM with the material properties varying continuously in an arbitrary manner is modeled as a multi-layered medium with the elastic modulus and mass density varying linearly in each sub-layer and continuous at the interfaces between two adjacent sub-layers. With this linearly inhomogeneous multi-layered model, the problem of a crack in a graded interfacial zone bonded to two homogeneous half-spaces or in a coating bonded to a homogeneous half-space subjected to the anti-plane shear impact load is investigated. Laplace and Fourier transforms and transfer matrix are applied to reduce the associated mixed boundary value problem to a Cauchy singular integral equation which is solved numerically in the Laplace transformed domain. The dynamic stress intensity factors (DSIF) are obtained by using the numerical technique of Laplace inversion.     

压电涂层-功能梯度压电界面层-基底结构的二维接触问题研究

在多层压电元件中,由于界面处材料成分和性质的突变,常常导致界面处应力集中,使得界面处出现开裂或蠕变现象,从而大大缩短了压电元件的使用寿命。功能梯度压电材料作为界面层,可有效的缓解界面材料不匹配导致的破坏。本文主要研究利用功能梯度压电材料界面层连接压电涂层和基底,分析三层结构在圆柱型压头作用下的力电响应。利用傅里叶积分变换技术,本文将压电涂层-功能梯度压电层-基底结构在刚性圆柱压头作用下的二维平面应变接触问题转化为带有柯西核的奇异积分方程。运用高斯-切比雪夫积分公式,将奇异积分方程转化为线性方程组并对其进行数值求解,得到压电涂层-功能梯度压电层-基底结构在圆柱形压头作用下的应力分布和电位移分布。数值结果表明,梯度压电材料参数的变化对结构中的力电响应具有重要的影响。本文研究结果对于利用功能梯度压电界面层消除界面处的应力不连续导致的界面破坏具有重要的理论指导意义,研究结果可为功能梯度压电材料界面层的设计提供帮助。     

Dynamic analysis of functionally graded plates using the hybrid Fourier-Laplace transform under thermomechanical loading

In this study, the analytical solution is presented for dynamic response of a simply supported functionally graded rectangular plate subjected to a lateral thermomechanical loading. The first-order and third-order shear deformation theories and the hybrid Fourier-Laplace transform method are used. The material properties of the plate, except Poisson’s ratio, are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is subjected to a heat flux on the bottom surface and convection on the upper surface. A third-order polynomial temperature profile is considered across the plate thickness with unknown constants. The constants are obtained by substituting the profile into the energy equation and applying the Galerkin method. The obtained temperature profile is considered along with the equations of motion. The governing partial differential equations are solved using the finite Fourier transformation method. Using the Laplace transform, the unknown variables are obtained in the Laplace domain. Applying the analytical Laplace inverse method, the solution in the time domain is derived. The computed results for static, free vibration, and dynamic problems are presented for different power law indices for a plate with simply supported boundary conditions. The results are validated with the known data reported in the literature. Furthermore, the results calculated by the analytical Laplace inversion method are compared with those obtained by the numerical Newmark method.     

Coupled thermal/mechanical analysis for the fracture of functionally graded materials under transient thermal loading

A comprehensive treatment of fracture of functionally graded materials (FGMs) is provided. It is assumed that the material properties depend only on the coordinate perpendicular to the crack surfaces and vary continuously along the crack faces. By using a laminated composite plate model to simulate the material non-homogeneity, an algorithm for solving the system based on Laplace transform and Fourier transform techniques is presented. Unlike earlier studies that considered certain assumed property distributions and a single crack problem, the current investigation studies multiple crack problem in the FGMs with arbitrarily varying material properties. Transient thermal stresses are presented. Project supported by the National Natural Science Foundation of China (Nos 10102004 and 19902003).     

Dynamic fracture mechanics analysis for composite material with material non-homogeneity in thickness direction

The problem considered here is the response of a non-homogeneous composite material containing some cracks subjected to dynamic loading. It is assumed that the composite material is orthotropic and all the material properties depend only on the coordinatey (along the thickness direction). In the analysis, the elastic region is divided into a number of plies of infinite length. The material properties are taken to be constants for each ply. By utilizing the Laplace transform and Fourier transform technique, the general solutions for plies are derived. The singular integral equations of the entire elastic region are obtained and solved by the virtual displacement principle. Attention is focused on the time-dependent full field solutions of stress intensity factor(SIF) and strain energy release rate. As a numerical illustration, the dynamic stress intensity factor of a substrate/functionally graded film structure with two cracks under suddenly applied forces on cracks face are presented for various material non-homogeneity parameters.     

A moving interface crack between two dissimilar functionally graded piezoelectric layers under electromechanical loading

The dynamic propagation of an interface crack between two dissimilar functionally graded piezoelectric material (FGPM) layers under anti-plane shear is analyzed using the integral transform method. The properties of the FGPM layers vary continuously along the thickness. The properties of the FGPM layers vary differently and the two layers are connected weak-discontinuously. A constant velocity Yoffe-type moving crack is considered. Fourier transform is used to reduce the problem to two sets of dual integral equations, which are then expressed to the Fredholm integral equations of the second kind. Numerical values on the dynamic energy release rate (DERR) are presented for the FGPM to show the effects on the electric loading, gradient of material properties, crack moving velocity, and thickness of layers. Followings are helpful to increase of the resistance of the interface crack propagation of FGPM: (a) certain direction and magnitude of the electric loading; (b) increase of the gradient of material properties; (c) increase of the material properties from the interface to the upper and lower free surface; (d) increase of the thickness of FGPM layer. The DERR increases or decreases with increase of the crack moving velocity.     

热冲击载荷作用下的含球形空腔的广义热弹性功能梯度球形各向同性体

This paper is concerned with the determination of thermoelastic displacement, stress and temperature in a functionally graded spherically isotropic infinite elastic medium having a spherical cavity, in the context of the linear theory of generalized thermoelasticity with two relaxation time parameters （Green and Lindsay theory）. The surface of cavity is stress-free and is subjected to a time-dependent thermal shock. The basic equations have been written in the form of a vector-matrix differential equation in the Laplace transform domain, which is then solved by an eigenvalue approach. Numerical inversion of the transforms is carried out using the Bellman method. Displacement, stress and temperature are computed and presented graphically. It is found that variation in the thermo-physical properties of a material strongly influences the response to loading. A comparative study with a corresponding homogeneous material is also made.     

Active transient elasto-acoustic response damping of a thick-walled liquid-coupled piezolaminated cylindrical vessel

The linear 3D piezoelasticity theory in conjunction with the versatile transfer matrix approach and the wave equation for the internal acoustic domain are employed for active non-stationary vibroacoustic response control of an arbitrarily thick, tri-laminate, fluid-filled, simply supported, piezocomposite cylindrical tank, excited by arbitrary (non-axisymmetric) time-dependent on-surface mechanical loads. The smart structure is composed of a supporting core layer of functionally graded orthotropic material perfectly bonded to inner and outer spatially distributed radially polarized functionally graded piezoceramic sensor and uniform force actuator (FGPM) layers. Active vibration damping is implemented by transferring the accumulated voltage on the sensor layer to the piezoelectric actuator layer in context of proportional and derivative control laws. Durbin's numerical inverse Laplace transform scheme is utilized to calculate the time response histories of the relevant interface displacement/stress components, center-point acoustic pressure, and actuator voltage, for selected loading configurations (i.e., concentrated step, impulse, and moving external loads). Numerical simulations demonstrate the effectiveness of the adopted distributed sensing/actuation configuration together with the active damping control strategy in suppressing the vibroacoustic response of a three-layered (Ba₂NaNb₅O₁₅/Al/PZT4) water-filled piezoelastic cylindrical tank. Limiting cases are considered and the validity of results is established by comparison with the available data as well as with the aid of a commercial finite element package.     

界面特性对功能梯度智能梁静动态响应的影响研究

采用状态空间法分析了两边简支的含压电夹层的功能梯度梁的静力弯曲和自由振动问题．为了考虑中间压电层与上、下功能梯度层之间的粘结效果,采用线性弹簧模型以模拟界面性能．假设上下功能梯度层的材料参数沿厚度连续变化,而压电层则是均匀材料,并且它们都是正交各向异性的．由于功能梯度梁的不均匀性使得直接求解比较困难,文中用层合模型来进行近似．数值算例中,分别考虑了压电层用于传感器或作动器的情形,分析了粘结界面完美程度对组合梁静力弯曲和自由振动频率的影响．