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.     

Aerothermoelastic behavior of supersonic rotating thin-walled beams made of functionally graded materials

In this study, a thin-walled beam made of functionally graded material (FGM) which is used as rotating blades in turbomachinery under aerothermoelastic loading is investigated. The governing equations, which are based on first-order shear deformation theory, include the effects of the presetting angle, the secondary warping, temperature gradient through the wall thickness of the beam and also the rotational speed. Moreover, quasi-steady aerodynamic pressure loadings are determined using first-order piston theory, and steady beam surface temperature is obtained from gas dynamics theory. Then, the blade partial differential equations are transformed into a set of ordinary differential equations using the extended Galerkin method. Finally, having solved the resulting structural–fluid–thermal eigenvalue system of equations, the effects of Mach number and geometric parameters on natural frequencies are presented. The results demonstrate that the natural frequencies decrease under aerothermoelastic loading at high Mach numbers.     

Elastoplastic deformations of rotating parabolic solid disks using Tresca''s yield criterion

Analytical solutions for the stress distribution in rotating parabolic solid disks are obtained. The analysis is based on Tresca's yield criterion, its associated flow rule and linear strain hardening. It is shown that, the deformation behavior of the convex parabolic disk is similar to that of the uniform thickness disk, but in the case of concave parabolic solid disk, it is different. In the latter, the plastic core consists of three different plastic regions with different mathematical forms of the yield criteria. Accordingly, three different stages of elastic–plastic deformation occur. All these stages of elastic–plastic deformation are studied in detail. It is also shown mathematically that in the limiting case the parabolic disk solution reduces to the solution of rotating uniform thickness solid disk.     

Analytical solution for the electromechanical behavior of a rotating functionally graded piezoelectric hollow shaft

In the present paper, the analytical solution for a radially piezoelectric functionally graded rotating hollow shaft is presented. The variation of material properties is assumed to follow a power law along the radial direction of the shaft. Two resulting fully coupled differential equations in terms of the displacement and electric potential are solved directly. Numerical results for different shaft geometries with different profiles of inhomogeneity are also graphically displayed.     

Von mises yield criterion and nonlinearly hardening variable thickness rotating annular disks with rigid inclusion

A computational model is developed to investigate inelastic deformations of variable thickness rotating annular disks mounted on rigid shafts. The von Mises yield condition and its flow rule are combined with Swift’s hardening law to simulate nonlinear hardening material behavior. An efficient numerical solution procedure is designed and used throughout to handle the nonlinearities associated with the von Mises yield condition and the boundary condition at the shaft–annular disk interface. The results of the computations are verified by comparison with an analytical solution employing Tresca’s criterion available in the literature. Inelastic stresses and deformations are calculated for rotating variable thickness disks described by two different commonly used disk profile functions i.e. power and exponential forms. Plastic limit angular velocities for these disks are calculated for different values of the geometric and hardening parameters. These critical angular velocities are found to increase as the edge thickness of the disk reduces. Lower plastic limit angular velocities are obtained for disks made of nonlinearly hardening materials.     

On the screw dislocation in a functionally graded material

This paper presents the stress field of a screw dislocation in a medium graded in y-direction. The medium is exponentially graded. For such a graded material theories of elasticity as well as gradient elasticity are applied. By means of the stress function technique we found exact analytical solutions of the corresponding master equations. Using the stress field, the Peach–Koehler force is given. The axial symmetry of a screw dislocation is lost due to the gradation in the y-direction.     

Coupled thermoelasticity of functionally graded cylindrical shells

The coupled thermoelastic response of a functionally graded circular cylindrical shell is studied. The coupled thermoelastic and the energy equations are simultaneously solved for a functionally graded axisymmetric cylindrical shell subjected to thermal shock load. A second-order shear deformation shell theory that accounts for the transverse shear strains and rotations is considered. Including the thermo-mechanical coupling and rotary inertia, a Galerkin finite element formulation in space domain and the Laplace transform in time domain are used to formulate the problem. The inverse Laplace transform is obtained using a numerical algorithm. The shell is graded through the thickness assuming a volume fraction of metal and ceramic, using a power law distribution. The results are validated with the known data in the literature.     

Natural frequencies of rotating functionally graded cylindrical shells

Love’s first approximation theory is used to analyze the natural frequencies of rotating functionally graded cylindrical shells.To verify the validity of the present method,the natural frequencies of the simply supported non-rotating isotropic cylindrical shell and the functionally graded cylindrical shell are compared with available published results.Good agreement is obtained.The effects of the power law index,the wave numbers along the x-and θ-directions,and the thickness-to-radius ratio on the natural frequencies of the simply supported rotating functionally graded cylindrical shell are investigated by several numerical examples.It is found that the fundamental frequencies of the backward waves increase with the increasing rotating speed,the fundamental frequencies of the forward waves decrease with the increasing rotating speed,and the forward and backward waves frequencies increase with the increasing thickness-to-radius ratio.     

Analysis of cracked functionally graded piezoelectric strip

The fracture behavior of a cracked strip under antiplane mechanical and inplane electrical loading is studied. A functionally graded piezoelectric strip with exponential material gradation is under consideration. The mechanical and electrical loading is combined via loading coupling factor. The problem of a graded piezoelectric strip containing a screw dislocation is solved. This solution results in stress and electric displacement components with Cauchy singularity. Based on the solution achieved for the dislocation, the distributed dislocation technique (DDT) is utilized to form any geometry of multiple cracks and analyze the behavior of a cracked strip under antiplane mechanical and inplane electrical loading. This technique is capable of the analysis of a strip with a system of interacting cracks. Several examples including strips with single crack, two straight cracks and two curved cracks are presented.     

Impact response of a functionally graded piezoelectric plate with a vertical crack

The dynamic fracture problem for a functionally graded piezoelectric plate containing a crack perpendicular to the free boundaries is considered in this study. It is assumed that the electroelastic properties of the medium vary continuously in the thickness direction. Integral transform techniques and dislocation density function are employed to reduce the problem to the solution of a singular integral equation. Mode I dynamic energy density factors are presented for an internal crack as well as an edge crack for various values of dimensionless parameters representing the size and location of the crack and the material nonhomogeneity.     

A preparation method of functionally graded materials with phase transition under shock loading

The propagation of phase boundary in a material undergoing shock induced irreversible phase transition is studied theoretically using a model based on simple-mixture rule. It is found that along with the decay of the phase boundary, a functionally graded material (FGM) forms in the mixed-phase region. Such FGMs are composed of parent phase and product phase, and the composition and physical properties are changing continuously without apparent macro-interfaces. The effect of stress boundary conditions on formation of the FGM is investigated in detail with a numerical method. The possibility of producing FGMs with impact method is proposed and the limit of this method is discussed.     

Thermomagnetic viscoelastic responses in a functionally graded hollow structure

This paper presents an analytical solution for the interaction of electric potentials,electric displacements,elastic deformations,and thermoelasticity,and describes electromagnetoelastic responses and perturbation of the magnetic field vector in hollow structures(cylinder or sphere),subjected to mechanical load and electric potential.The material properties,thermal expansion coefficient and magnetic permeability of the structure are assumed to be graded in the radial direction by a power law distribution.In the present model we consider the solution for the case of a hollow structure made of viscoelastic isotropic material,reinforced by elastic isotropic fibers,this material is considered as structurally anisotropic material.The exact solutions for stresses and perturbations of the magnetic field vector in FGM hollow structures are determined using the infinitesimal theory of magnetothermoelasticity,and then the hollow structure model with viscoelastic material is solved using the correspondence principle and Illyushin’s approximation method.Finally,numerical results are carried out and discussed.     

Micromechanics-based thermoelastic model for functionally graded particulate materials with particle interactions

Thermoelastic behavior of functionally graded particulate materials is investigated with a micromechanical approach. Based on a special representative volume element constructed to represent the graded microstructure of a macroscopic material point, the relation between the averaged strains of the particle and matrix phases is derived with pair-wise particle interactions, and a set of governing equations for the thermoelastic behavior of functionally graded materials is presented. The effective coefficient of thermal expansion at a material point is solved through the overall averaged strain of two phases induced by temperature change under the stress-free condition, and is shown to exhibit a weak anisotropy due to the particle interactions within the graded microstructure. When the material gradient is eliminated, the proposed model predicts the effective coefficient of thermal expansion for uniform composites as expected. If the particle interactions are disregarded, the proposed model recovers the Kerner model. The proposed semi-analytical scheme is consistent and general, and can handle any thermal loading variation. As examples, the thermal stress distributions of graded thermal barrier coatings are solved for two types of thermal loading: uniform temperature change and steady-state heat conduction in the gradation direction.     

Crack in functionally graded piezoelectric strip bonded to elastic surface layers under electromechanical loading

Solved is the problem of a crack in a functionally graded piezoelectric material (FGPM) bonded to two elastic surface layers. It is assumed that the elastic stiffness, piezoelectric constant, and dielectric permittivity of the FGPM vary continuously along the thickness of the strip. The outside layers are under antiplane mechanical loading and in-plane electric loading. The solution involves solving singular integral equations by application of the Gauss–Jacobi integration formula. Numerical calculations are carried out to obtain the energy density factors. Their variations with the geometric, loading and material parameters are shown graphically.     

Periodic cracks in a functionally graded piezoelectric layer bonded to a piezoelectric half-plane

Studied is the problem of a periodic array of cracks in a functionally graded piezoelectric strip bonded to a homogeneous piezoelectric material. The properties of the functionally graded piezoelectric strip, such as elastic modulus, piezoelectric constant and dielectric constant, are assumed in exponential forms and vary along the crack direction. The crack surface condition is assumed to be electrically impermeable or permeable. Integral transform and dislocation density functions are employed to reduce the problem to the solution of a system of singular integral equations. The effects of the periodic crack spacing, material constants and the geometry parameters on the stress intensity factor, the energy release ratio and the energy density factor are studied.     

Exact solutions for radial deformations of a functionally graded isotropic and incompressible second-order elastic cylinder

We find closed-form solutions for axisymmetric plane strain deformations of a functionally graded circular cylinder comprised of an isotropic and incompressible second-order elastic material with moduli varying only in the radial direction. Cylinder's inner and outer surfaces are loaded by hydrostatic pressures. These solutions are specialized to cases where only one of the two surfaces is loaded. It is found that for a linear through-the-thickness variation of the elastic moduli, the hoop stress for the first-order solution (or in a cylinder comprised of a linear elastic material) is a constant but that for the second-order solution varies through the thickness. The radial displacement, the radial stress and the hoop stress do not depend upon the second-order elastic constant but the hydrostatic pressure and hence the axial stress depends upon it. When the two elastic moduli vary as the radius raised to the power two or four, the radial and the hoop stresses in an infinite space with a pressurized cylindrical cavity equal the pressure in the cavity. For an affine variation of the elastic moduli, the hoop stress in an internally loaded cylinder made of a linear elastic isotropic and incompressible material at the point is the same as that in a homogeneous cylinder. Here R_in and R_ou equal, respectively, the inner and the outer radius of the undeformed cylinder and R the radial coordinate of a point in the unstressed reference configuration.     

High-order shear theory for static analysis of functionally graded plates with porosities

The bending responses of porous functionally graded (FG) thick rectangular plates are investigated according to a high-order shear deformation theory. Both the effect of shear strain and normal deformation are included in the present theory and so it does not need any shear correction factor. The equilibrium equations according to the porous FG plates are derived. The solution to the problem is derived by using Navier's technique. Numerical results have been reported and compared with those available in the open literature for non-porous plates. The effects of the exponent graded and porosity factors are investigated.     

Effects of SH waves in a functionally graded plate

A computational method is presented to investigate SH waves in functionally graded material (FGM) plates. The FGM plate is first divided into quadratic layer elements (QLEs), in which the material properties are assumed as a quadratic function in the thickness direction. A general solution for the equation of motion governing the QLE has been derived. The general solution is then used together with the boundary and continuity conditions to obtain the displacement and stress in the wave number domain for an arbitrary FGM plate. The displacements and stresses in the frequency domain and time domain are obtained using inverse Fourier integration. Furthermore, a simple integral technique is also proposed for evaluating modified Bessel functions with complex valued order. Numerical examples are presented to demonstrate this numerical technique for SH waves propagating in FGM plates.     

Reissner-Sagoci problem for a homogeneous coating on a functionally graded half-space

This paper is concerned with the axisymmetric elastostatic problem related to the rotation of a rigid punch which is bonded to the surface of a nonhomogeneous half-space. The half-space is composed of an isotropic homogeneous coating in the form of layer, which is attached to the functionally graded half-space. The shear modulus of the FGM is assumed to vary in the direction of axis Oz normal to the boundary as μ₁(z) = μ₀(1 + αz)^β, where μ₀, α, β are positive constants. The punch undergoes rotation due to the action of the internal loads. By using Hankel's integral transforms, the mixed boundary value problem is reduced to dual integral equations, and next, to a Fredholm's integral equation of the second kind, which is solved numerically for the case of β = 2. The final results show the effect of non-homogeneity on the shear stresses and an unknown moment of punch rotation.     

Recursive geometric integrators for wave propagation in a functionally graded multilayered elastic medium

The differential equations governing transfer and stiffness matrices and acoustic impedance for a functionally graded generally anisotropic magneto-electro-elastic medium have been obtained. It is shown that the transfer matrix satisfies a linear 1st order matrix differential equation, while the stiffness matrix satisfies a nonlinear Riccati equation. For a thin nonhomogeneous layer, approximate solutions with different levels of accuracy have been formulated in the form of a transfer matrix using a geometrical integration in the form of a Magnus expansion. This integration method preserves qualitative features of the exact solution of the differential equation, in particular energy conservation. The wave propagation solution for a thick layer or a multilayered structure of inhomogeneous layers is obtained recursively from the thin layer solutions. Since the transfer matrix solution becomes computationally unstable with increase of frequency or layer thickness, we reformulate the solution in the form of a stable stiffness-matrix solution which is obtained from the relation of the stiffness matrices to the transfer matrices. Using an efficient recursive algorithm, the stiffness matrices of the thin nonhomogeneous layer are combined to obtain the total stiffness matrix for an arbitrary functionally graded multilayered system. It is shown that the round-off error for the stiffness-matrix recursive algorithm is higher than that for the transfer matrices. To optimize the recursive procedure, a computationally stable hybrid method is proposed which first starts the recursive computation with the transfer matrices and then, as the thickness increases, transits to the stiffness matrix recursive algorithm. Numerical results show this solution to be stable and efficient. As an application example, we calculate the surface wave velocity dispersion for a functionally graded coating on a semispace.