Thermal stresses analysis of ceramic/metal functionally gradient material cylinder

ＩｎｔｒｏｄｕｃｔｉｏｎＣｅｒａｍｉｃ／ｍｅｔａｌｆｕｎｃｔｉｏｎａｌｙｇｒａｄｉｅｎｔｍａｔｅｒｉａｌ（ＦＧＭ）ｈａｓｇｒａｄｉｅｎｔｉｎｇｒｅｄｉｅｎｔｄｉｓｔｒｉｂｕｔｉｏｎ,ｓｏｉｔｓｐｒｏｐｅｒｔｉｅｓｏｆｔｈｅｈｅａｔｔｒａｎｓｆｅｒａｎ...     

Anti-plane shear crack in a functionally gradient piezoelectric material

The main objective of this paper is to study the singular nature of the crack-tip stress and electric displacement field in a functionally gradient piezoelectric medium having material coefficients with a discontinuous derivative. The problem is considered for the simplest possible loading and geometry, namely, the anti-plane shear stress and electric displacement in-plane of two bonded half spaces in which the crack is parallel to the interface. It is shown that the square-root singularity of the crack-tip stress field and electric displacement field is unaffected by the discontinuity in the derivative of the material coefficients. The problem is solved for the case of a finite crack and extensive results are given for the stress intensity factors, electric displacement intensity factors, and the energy release rate. Project supported by the National Natural Science Foundation of China (No. 10072041), the National Excellent Young Scholar Fund, of China (No. 10125209) and the Teaching and Research Award Program for Outstanding Young Teachers in Higher Education Institutions of MOE, P. R. C..     

ANTI-PLANE FRACTURE ANALYSIS OF FUNCTIONALLY GRADIENT MATERIAL INFINITE STRIP WITH FINITE WIDTH

The special case of a crack under mode III conditions was treated, lying parallel to the edges of an infinite strip with finite width and with the shear modulus varying exponentially perpendicular to the edges. By using Fourier transforms the problem was formulated in terms of a singular integral equation. It was numerically solved by representing the unknown dislocation density by a truncated series of Chebyshev polynomials leading to a linear system of equations. The stress intensity factor (SIF) results were discussed with respect to the influences of different geometric parameters and the strength of the non-homogeneity. It was indicated that the SIF increases with the increase of the crack length and decreases with the increase of the rigidity of the material in the vicinity of crack. The SIF of narrow strip is very sensitive to the change of the non-homogeneity parameter and its variation is complicated. With the increase of the non-homogeneity parameter, the stress intensity factor may increase, decrease or keep constant, which is mainly determined by the strip width and the relative crack location. If the crack is located at the midline of the strip or if the strip is wide, the stress intensity factor is not sensitive to the material non-homogeneity parameter.     

功能梯度材料板件三维分析的半解析梯度有限元法

将半解析有限元与梯度有限元相结合,形成一种半解析梯度有限元来求解功能梯度材料板件问题。该方法兼有有限元法的适应性强、程序统一,半解析有限元法的节省单元与计算工作量,梯度有限元法的适应构件内部材料性能任意梯度分布等特点,并实现用一维数值计算给出构件三维分析结果。算例分析表明了方法的精度、功能与上述特点,充分揭示了功能梯度材料板件力学响应的三维形态。半解析梯度有限元法可推广应用到其他功能梯度材料面结构的各类分析中。     

Thermal post-bunkling analyses of functionally graded material rod

The non-linear governing differential equations of immovably simply supported functionally graded material (FGM) rod subjected to thermal loads were derived. The thermal post-buckling behaviors of FGM rod made of ZrO2 and Ti-6Al-4Vwere analyzed by shooting method. Firstly, the thermal post-buckling equilibrium paths of the FGM rod with different gradient index in the uniform temperature field were plotted, and compared with the behaviors of the homogeneous rods made of ZrO2 and Ti-6A1-4V materials, respectively. For given value of end rotation angles, the influence of gradient index on the thermal post-buckling behaviors of FGM rod was discussed. Secondly, the thermal post-buckling characteristics of the FGM rod were analyzed when the temperature difference parameter is changed while the bottom temperature parameter remains constant, and when the bottom temperature parameter is changed while the temperature difference parameter remains constant, and compared with the characteristics of the two homogeneous material rods.     

热环境中功能梯度圆形薄板的混沌运动

针对陶瓷-金属功能梯度圆板,同时考虑几何非线性、材料物性参数随温度变化且材料组分沿厚度方向按幂律分布的情况,应用虚功原理给出了热载荷与横向简谐载荷共同作用下的非线性振动偏微分方程。在固支无滑动的边界条件下,通过引入位移函数,利用伽辽金方法得到了达芬型非线性动力学方程。利用Melnikov方法,给出了热环境中功能梯度圆板可能发生混沌运动的临界条件。通过数值算例,给出了不同体积分数指数和温度的同宿分岔曲线,平面相图和庞加莱映射图,讨论其对临界条件的影响,证实了系统混沌运动的存在。通过分岔图和与其相对应的最大李雅普诺夫指数图,分析了激励频率和激励幅值对倍周期分岔的影响及变化规律,发现系统可出现周期、倍周期和混沌等复杂动力学响应。     

Several three-dimensional solutions for transversely isotropic functionally graded material plate welded with circular inclusion

The problem of a transversely isotropic functionally graded material (FGM) plate welded with a circular inclusion is considered. The analysis starts with the generalized England-Spencer plate theory for transversely isotropic FGM plates, which expresses a three-dimensional (3D) general solution in terms of four analytic functions. Several analytical solutions are then obtained for an infinite FGM plate welded with a circular inclusion and subjected to the loads at infinity. Three different cases are considered, i.e., a rigid circular inclusion fixed in the space, a rigid circular inclusion rotating about the x-, y-, and z-axes, and an elastic circular inclusion with different material constants from the plate itself. The static responses of the plate and/or the inclusion are investigated through numerical examples.     

Higher-order crack tip fields for functionally graded material plate with transverse shear deformation

The crack tip fields are investigated for a cracked functionally graded material (FGM) plate by Reissner’s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson’s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner’s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner’s effect when the in-homogeneity parameter approaches zero.     

Crack propagation in a functionally graded plate under thermal shock

Summary Thermal cracking in a ceramic/metal functionally graded plate is discussed. When a functionally graded plate is cooled from high temperature, curved or straight crack paths are experimentally observed on the ceramic surface. One of the reasons that make the crack paths to differ are the thermal or mechanical conditions. In order to clarify the influence of these conditions on the crack path, the crack propagation is simulated using finite element method. Received 29 September 1998; accepted for publication 2 August 1999     

Neural network optimization of material composition of a functionally graded material plate at arbitrary temperature range and temperature rise

Summary A neural network model is applied to optimization problems of material compositions for a functionally graded material plate with arbitrarily distributed and continuously varied material properties in the thickness direction. Unsteady temperature distribution is evaluated by taking into account the bounds of the number of the layers. Thermal stress components for an infinite functionally graded material plate are formulated under traction-free mechanical conditions. As a numerical example, a plate composed of zirconium oxide and titanium alloy is considered. In the optimization problem of minimizing the thermal stress distribution, the numerical calculations are carried out making use of the neural network. The optimum material composition is determined by taking into account the effect of temperature-dependence of material properties. The results obtained by neural network and ordinary nonlinear programming method are compared. Received 3 March 1998; accepted for publication 22 May 1998     

功能梯度双材料弱/微间断界面的冲击断裂分析

提出强间断、弱间断、微间断和全连续界面的概念与分类,建立功能梯度弹性双材料弱间断 界面冲击断裂问题的力学模型,采用积分变换法推导问题的Cauchy奇异积分方程,并用配 点法求得数值解. 分析表明,弱/微间断性对于FGMs界面裂纹应力强度因子有着重要影响, 而且微间断性是优于弱间断性的一种界面力学性能连接关系. 以FGMs界面某一侧 的力学性能函数在界面处的Taylor展开式的低阶项作为界面另一侧的力学性能函数,便可 以使FGMs界面成为``微间断'界面. 界面的一阶微间断对应力强度因子的减小作用较为明 显,而高阶(二阶及以上)微间断对应力强度因子的影响较小. 减小界面的弱间断程度或使 FGMs界面具备``微间断性',都将利于提高功能梯度双材料界面抗冲击断裂能力,在一定 程度上达到界面增韧的目的.     

Thermoelastic analysis of functionally graded annulus with arbitrary gradient

A thermoelastic problem of a circular annulus made of functionally graded materials with an arbitrary gradient is investigated. Different from previous works, our analysis neither requires a special form of the gradient of material properties nor demands partitioning the entire structure into a multilayered homogeneous structure. Instead, we propose a new method for solving the thermoelastic problem of a functionally graded circular annulus by transforming it to a Fredholm integral equation. The distribution of thermal stresses and radial displacement can be obtained by solving the resulting equation. Illustrative examples are given to show the effects of varying gradients on the thermal stresses and radial displacement for given temperature changes at the inner and outer surfaces. The results indicate that the thermal stresses can be relaxed for specified gradients, which is beneficial to design an inhomogeneous annulus to maintain structural integrity.     

梯度厚板热应力分层计算方法的改进

为了使功能梯度材料厚板的分层计算具有更快的收敛速度和更好的精度,而采用每层热力系数与梯度材料厚度变化为指数函数的分层方法,来平滑均匀分层法引起材料热力系数的跳跃,使热力系数更逼近实际,并将该方法命名为指数函数有限分层法.通过算例表明,在精度为0.1的情况下,均匀分层法的分层层数是指数函数分层法的分层层数的2倍多,并且随着精度的提高,分层层数相差的倍数就越大.     

Crack tip field in functionally gradient material with exponential variation of elastic constants in two directions

This paper presents an exact solution of the crack tip field in functionally gradient material with exponential variation of elastic constants. The dimensionless Poisson's ratios ν ₀ of the engineering materials (iron, glass . . .) are far less than one; therefore, neglecting them, one can simplify the basic equation and the exact solution is easy to obtain. Although the exact solution for the case ν ₀≠0 is also obtained, it is very complicated and the main result is the same with the case ν ₀=0 (it will be dealt with in Appendix VII). It has been found that the exponential term exp (ax+by) in the constitutive equations becomes exp (ax/2+by/2−kr/2) in the exact solution. The English text was polished by Keren Wang.     

粘贴压电层功能梯度材料Timoshenko梁的热过屈曲分析

研究了上下表面粘贴压电层的功能梯度材料Timoshenko梁在升温及电场作用下的过屈曲行为。在精确考虑轴线伸长和一阶横向剪切变形的基础上,建立了压电功能梯度Timoshenko层合梁在热-电-机械载荷作用下的几何非线性控制方程。其中,假设功能梯度的材料性质沿厚度方向按照幂函数连续变化,压电层为各向同性均匀材料。采用打靶法数值求解所得强非线性边值问题,获得了在均匀电场和横向非均匀升温场内两端固定Timoshenko梁的静态非线性屈曲和过屈曲数值解。并给出了梁的变形随热、电载荷及材料梯度参数变化的特性曲线。结果表明,通过施加电压在压电层产生拉应力可以有效地提高梁的热屈曲临界载荷,延缓热过屈曲发生。由于材料在横向的非均匀性,即使在均匀升温和均匀电场作用下,也会产生拉-弯耦合效应。但是对于两端固定的压电-功能梯度材料梁,在横向非均匀升温下过屈曲变形仍然是分叉形的。     

Bending of functionally graded piezoelectric rectangular plates

By introducing two displacement functions as well as two stress functions, two independent state equations with variable coefficients are derived from the three-dimensional theory equations of piezoelasticity for transverse isotropy. A laminated approximation is used to transform the state equations to those with constant coefficients in each sub-layer. The bending problem of a functionally graded rectangular plate is then analyzed based on the state equations. Numerical results are presented and the effect of material gradient index is discussed. Supported by the National Natural Sciences Foundation of China (No. 10002016).     

Free vibration of functionally graded material beams with surface-bonded piezoelectric layers in thermal environment

Free vibration of statically thermal postbuckled functionally graded material （FGM） beams with surface-bonded piezoelectric layers subject to both temperature rise and voltage is studied. By accurately considering the axial extension and based on the Euler-Bernoulli beam theory, geometrically nonlinear dynamic governing equations for FGM beams with surface-bonded piezoelectric layers subject to thermo-electro- mechanical loadings are formulated. It is assumed that the material properties of the middle FGM layer vary continuously as a power law function of the thickness coordinate, and the piezoelectric layers are isotropic and homogenous. By assuming that the amplitude of the beam vibration is small and its response is harmonic, the above mentioned non-linear partial differential equations are reduced to two sets of coupled ordinary differential equations. One is for the postbuckling, and the other is for the linear vibration of the beam superimposed upon the postbuckled configuration. Using a shooting method to solve the two sets of ordinary differential equations with fixed-fixed boundary conditions numerically, the response of postbuckling and free vibration in the vicinity of the postbuckled configuration of the beam with fixed-fixed ends and subject to transversely nonuniform heating and uniform electric field is obtained. Thermo-electric postbuckling equilibrium paths and characteristic curves of the first three natural frequencies versus the temperature, the electricity, and the material gradient parameters are plotted. It is found that the three lowest frequencies of the prebuckled beam decrease with the increase of the temperature, but those of a buckled beam increase monotonically with the temperature rise. The results also show that the tensional force produced in the piezoelectric layers by the voltage can efficiently increase the critical buckling temperature and the natural frequency.     

Homogenized and classical expressions for static bending solutions for functionally graded material Levinson beams

The exact relationship between the bending solutions of functionally graded material (FGM) beams based on the Levinson beam theory and those of the corresponding homogenous beams based on the classical beam theory is presented for the material properties of the FGM beams changing continuously in the thickness direction. The deflection, the rotational angle, the bending moment, and the shear force of FGM Levinson beams (FGMLBs) are given analytically in terms of the deflection of the reference homogenous Euler-Bernoulli beams (HEBBs) with the same loading, geometry, and end supports. Consequently, the solution of the bending of non-homogenous Levinson beams can be simplified to the calculation of transition coefficients, which can be easily determined by variation of the gradient of material properties and the geometry of beams. This is because the classical beam theory solutions of homogenous beams can be easily determined or are available in the textbook of material strength under a variety of boundary conditions. As examples, for different end constraints, particular solutions are given for the FGMLBs under specified loadings to illustrate validity of this approach. These analytical solutions can be used as benchmarks to check numerical results in the investigation of static bending of FGM beams based on higher-order shear deformation theories.     

Sound radiation of a functionally graded material cylindrical shell in water by mobility method

Based on the fundamental dynamic equations of functionally graded material (FGM) cylindrical shell, this paper investigates the sound radiation of vibrational FGM shell in water by mobility method. This model takes into account the exterior fluid loading due to the sound press radiated by the FGM shell. The FGM cylindrical shell was excited by a harmonic line radial force uniformly distributing along the generator. The FGM shell equations of motion, the Helmholtz equation in the exterior fluid medium and the continuity equation at fluid-shell interface are used in this vibroacoustic problem. The expressions of sound radiation efficiency and sound field of the FGM shell have been derived by mobility method. Radiation efficiency, modal mobility and the directivity pattern of the sound field are solved numerically. In particular, radiation efficiency and directivity pattern with various power law index are analyzed.     

Analytical solutions for plane problem of functionally graded magnetoelectric cantilever beam

In this paper, an exact analytical solution is presented for a transversely isotropic functionally graded magneto-electro-elastic (FGMEE) cantilever beam, which is subjected to a uniform load on its upper surface, as well as the concentrated force and moment at the free end. This solution can be applied for any form of gradient distribution. For the basic equations of plane problem, all the partial differential equations governing the stress field, electric, and magnetic potentials are derived. Then, the expressions of Airy stress, electric, and magnetic potential functions are assumed as quadratic polynomials of the longitudinal coordinate. Based on all the boundary conditions, the exact expressions of the three functions can be determined. As numerical examples, the material parameters are set as exponential and linear distributions in the thickness direction. The effects of the material parameters on the mechanical, electric, and magnetic fields of the cantilever beam are analyzed in detail.