Next Article in Journal
Analysis of the Behavior of SMA Mixtures with Different Fillers Through the Semicircular Bend (SCB) Fracture Test
Next Article in Special Issue
Stress Concentration and Damage Factor Due to Central Elliptical Hole in Functionally Graded Panels Subjected to Uniform Tensile Traction
Previous Article in Journal
Optimal Design of Sandwich Composite Cradle for Computed Tomography Instrument by Analyzing the Structural Performance and X-ray Transmission Rate
Previous Article in Special Issue
Asymptotic Solution and Numerical Simulation of Lamb Waves in Functionally Graded Viscoelastic Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Finite Element Implementation of Functionally Graded Materials

by
Emilio Martínez-Pañeda
Department of Engineering, Cambridge University, Cambridge CB2 1PZ, UK
Materials 2019, 12(2), 287; https://doi.org/10.3390/ma12020287
Submission received: 24 December 2018 / Revised: 13 January 2019 / Accepted: 15 January 2019 / Published: 17 January 2019

Abstract

:
We investigate the numerical implementation of functionally graded properties in the context of the finite element method. The macroscopic variation of elastic properties inherent to functionally graded materials (FGMs) is introduced at the element level by means of the two most commonly used schemes: (i) nodal based gradation, often via an auxiliary (non-physical) temperature-dependence, and (ii) Gauss integration point based gradation. These formulations are extensively compared by solving a number of paradigmatic boundary value problems for which analytical solutions can be obtained. The nature of the notable differences revealed by the results is investigated in detail. We provide a user subroutine for the finite element package ABAQUS to overcome the limitations of the most popular approach for implementing FGMs in commercial software. The use of reliable, element-based formulations to define the material property variation could be key in fracture assessment of FGMs and other non-homogeneous materials.

Graphical Abstract

1. Introduction

There is an emerging interest in the analysis of the mechanical response of materials with spatially varying properties. New manufacturing technologies make it possible to engineer materials with functionally graded microstructures, so-called functionally graded materials (FGMs). The resulting spatial variation of material properties eliminates stress discontinuities at material interfaces and optimizes material performance under non-uniform service conditions. For example, the performance of coatings subjected to large thermal gradients can be significantly improved by using metal-ceramic FGMs [1], which combine the thermal and corrosive resistance of ceramics with the mechanical strength and high tenacity of metals. In addition, FGMs are now employed in a host of commercial applications, ranging from cutting tools to biomedical devices [2]. This widespread use of FGMs is largely due to their capacity to reduce residual stresses [3], increase the strength of joints [4], and tailor material microstructure to specific service requirements [5].
The complexity and cost associated with the manufacture and testing of functionally graded specimens has intensified the use of numerical tools to analyse their mechanical response. Although a variety of numerical techniques have been used, including mesh-free methods [6,7] and enriched formulations [8,9], the finite element method is by far the most popular approach [10,11,12,13,14,15]. Several formulations have been proposed to accurately capture a smooth material gradient by defining the material property variation at the element level [11,12,13]. While these formulations are been widely and indistinctly used, a performance assessment of the different types of graded elements has not been conducted yet. We investigate the performance of different types of functionally graded elements by comparing with analytical solutions of paradigmatic boundary value problems. We show that notable differences can be attained and that the most common approach in commercial software entails a number of limitations. An alternative implementation is presented in the context of the commercial finite element package ABAQUS.

2. Numerical Formulation

The assignment of material properties in the numerical model must reflect the material property distribution in the functionally graded specimen under consideration. However, an accurate characterization of the material gradient is not a straightforward task. Typically, the information available is the spatial variation of the volume fractions of constituent materials, which is provided as input to the production technique [16]. The macroscopic material property variation does not tend to mirror the volume fraction profile, but one can estimate the former from the latter by using homogenization laws [17]. However, the micromechanical assumptions upon which these theoretical mixing laws are built may hinder an accurate characterization of the macroscopic variation of material properties. An alternative approach is to determine the material property variation directly by experimentation. For example, by producing and testing individual homogeneous specimens with a range of volume fractions [18], by testing the graded material through indentation or ultrasonic techniques [19], or by cutting and testing small samples from a larger graded specimen [20]. Capturing this material gradation profile in the numerical model is key to designing optimal FGM specimens, as well as reproducing and gaining insight into experimental results.
From the numerical perspective, material properties can vary between elements or between nodes and integration points. Numerical works in the mechanics of functionally graded materials can be classified into two large groups depending on their approach to the implementation of the material gradient, see Figure 1, using either homogeneous elements (see, e.g., [10]) or graded elements (see [15] and references therein). The former is appropriate for layered functionally graded composites, but it constitutes a poor approximation otherwise. Assuming constant material properties within each element leads to a discontinuous step-type variation and requires uniform meshing along the material gradation direction. A material property variation at the element level is generally more appropriate and different graded elements formulations have been proposed.

2.1. Gauss Integration Point-Based Variation

Consider an isoparametric finite element with n number of nodes, the displacement field u ( x ) is interpolated from the nodal values u ^ i as
u = i = 1 n N i ( ξ , η , ζ ) u ^ i ,
where i is a given node and N i are the shape functions. For example, in an eight-node quadrilateral element, the shape functions read, for the corner nodes
N i = 1 4 1 + ξ ξ i 1 + η η i ξ ξ i + η η i 1
and for the mid-side nodes,
N i = 1 2 1 ξ 2 1 + η η i , for   ξ i = 0 ,
N i = 1 2 1 + ξ ξ i 1 η 2 , for   η i = 0 ,
with ξ , η denoting the intrinsic coordinates in the interval [ 1 , 1 ] , and ξ i , η i denoting the local coordinates of node i. Accordingly, the strain field ε ( x ) is computed from the nodal displacements by means of the strain-displacement matrix
ε = i = 1 n B i ξ , η , ζ u ^ i
with the matrix B i containing the appropriate derivatives of the shape functions N i . For example, in a plane strain element, the strain-displacement matrix for a node i reads
B i ξ , η = N i ξ , η / x 0 0 N i ξ , η / y 0 0 N i ξ , η / y N i ξ , η / x
so as to compute the strain components ε x x , ε y y , ε z z , γ x y from the nodal displacements.
Let us assume linear elastic behaviour, which is arguably appropriate for ceramic-based FGMs. The Cauchy stress field σ is related to the strain tensor ε through a spatially varying constitutive matrix C x as
σ = C x ε .
The principle of virtual work yields the relation between the deformation work given by the element and the elemental nodal force vector F e as
K e u ^ i = F e ,
where K e is the element strain-displacement matrix—see, for example, Ref. [21]. Accordingly, the element stiffness matrix over the volume of the element V e reads
K e = V e B e T C x B e d V ,
where B e is the element stiffness matrix, given by the assembly of B i over n nodes. Therefore, the linear elastic stiffness matrix is defined to match the material gradation profile at the Gauss integration points. This basic finite element formulation for functionally graded solids was presented by Santare and Lambros [12].

2.2. Nodal-Based Variation via Temperature Dependence

An alternative approach to develop a formulation for graded finite elements was proposed by Kim and Paulino [13]. They propose a generalized isoparametric finite element formulation where the same shape functions are employed to interpolate the nodal displacements, the geometry, and the material properties. Thus, consider a standard isoparametric formulation where the spatial coordinates ( x , y , z ) are interpolated as
x = i = 1 n N i x i , y = i = 1 n N i y i , z = i = 1 n N i z i .
The isoparametric concept can be generalized to interpolate the spatially varying Young’s modulus E x and Poisson’s ratio ν x as
E = i = 1 n N i E i , ν = i = 1 n N i ν i ,
where E i and ν i are the elastic properties defined at each node i. Hence, the material gradient is defined precisely at the nodes and subsequently interpolated to the Gauss integration points to compute the stresses through Equation (7).
A generalized isoparametric graded finite element can be easily implemented into a commercial finite element package by taking advantage of the possibility of defining temperature-dependent material properties [11,15]. For example, one can define E as a function of the temperature and provide the specimen with an initial temperature distribution that matches the Young’s modulus variation desired. Here, the temperature has no physical meaning and unwanted thermal strains are suppressed by assigning a zero thermal expansion coefficient. Since the temperature field is defined at the nodes and subsequently interpolated to the Gauss integration points, this technique constitutes a straightforward implementation of a generalized isoparametric graded element, enjoying great popularity. Evident drawbacks are the inability to (i) model thermomechanical problems, and (ii) define different profiles for Poisson’s ratio and Young’s modulus. Furthermore, to obtain a consistent variation of mechanical and thermal strains, many commercial codes interpolate nodal temperature values using shape functions one order lower than those used for the nodal displacements. Consequently, there is an inherent error in the presence of a nonlinear material gradation profile, as sketched in Figure 2. The implications of adopting this technique, relative to the Gauss-based approach defined in Section 2.1, are explored here.

3. Results

The variation in elastic properties inherent to FGMs is implemented at the element level by making use of user subroutines within the commercial finite element package ABAQUS. The graded elements described in Section 2 can be readily implemented by using a USDFLD subroutine, for a Gauss points-based implementation, or a UFIELD subroutine, for a nodal-based graded element. In addition, as discussed above, the temperature can be used as an auxiliary field to effectively implement a generalized isoparametric graded element. The user must provide the material properties as a function of a user defined field (or temperature). Then, a suitable field (or temperature distribution) is defined to match the material property variation desired. A direct comparison between the two approaches in terms of computational time is hindered by their different implementations; nevertheless, the sampling of material properties at integration points or nodes is achieved at a negligible computational cost. The performance of different types of graded elements will be benchmarked by considering a Gauss point-based implementation (Section 2.1) and, via temperature dependent properties, a generalized isoparametric approach (Section 2.2). For the sake of simplicity, we will consider bi-dimensional problems and four quadrilateral element types (see Figure 3): linear elements with reduced integration (Q4R), linear elements with full integration (Q4), quadratic elements with reduced integration (Q8R), and quadratic elements with full integration (Q8).
Three plane problems for which analytical solutions can be obtained will be addressed, as shown in Figure 4. Young’s modulus will be varied along the x-direction and Poisson’s ratio will be assumed to be constant. We assume plane stress conditions. As sketched in Figure 4, the three boundary value problems considered involve a functionally graded plate being subjected to (i) uniform displacement perpendicular to the material gradient direction, (ii) uniform traction perpendicular to the material gradient direction, and (iii) uniform traction in the direction parallel to material gradation.
In all cases, we assume that the Young’s modulus varies exponentially as
E ( x ) = E 0 exp β x
with E 0 and β being material constants. This choice is motivated by the existence of analytical solutions for functionally graded solids exhibiting an exponential variation of the elastic properties; see, for example, Refs. [13,22]. Many other functions have been employed in the literature (see Ref. [23] for a review), but our choice is appropriate for our aim: comparing on equal footing different graded finite element implementations. Consistent units are throughout the manuscript and, therefore, units will be omitted subsequently. We consider a width w = 1 , a total height of 2 h = 3 , and choose E 0 = 1 and β = ln 8 so as to vary E gradually in the x-direction from E ( 0 ) = 1 to E ( w ) = 8 .

3.1. Uniform Displacement Perpendicular to the Material Gradient Direction

Consider first the case of an FGM plate subjected to a remote strain ε 0 = u 0 / h , where u 0 denotes the displacement in the remote boundary and h denotes half the height of the plate. A constant Poisson’s ratio of ν = 0.3 throughout the plate is assumed. The relevant stress component is given by
σ y y x , y = E ( x ) u 0 / h .
This analytical solution is compared in Figure 5 with the finite element results obtained for the element types and graded element formulations described above. A remote displacement of u 0 = 0.2 is prescribed in all cases. As shown in the insets of Figure 5, a uniform mesh of 4 by 12 elements is employed.
Results reveal differences between the different types of graded elements. Consider first the case of the linear element with full integration Q4—see Figure 5a. The Gauss integration point-based approach accurately captures the material gradient depicted by the analytical solution. In fact, the numerical result is exact at the integration points since the displacement field is linear; a single Q4 element will suffice to capture the FGM response. However, the generalized isoparametric formulation (via temperature-dependent properties) exhibits a step-type variation with constant stress in each element. This behaviour is inherently related to how ABAQUS interpolates nodal temperature values. As many other finite element packages, ABAQUS interpolates nodal temperatures with shape functions that are one order lower than those used for the displacements, so as to obtain an equivalent distribution of mechanical and thermal strains. An average value of the temperature in the nodes is passed to the integration points when using linear elements and a linear variation is assumed in quadratic elements. In agreement with expectations, the use of linear elements with reduced integration (Q4R, see Figure 5b) exhibits the response inherent to homogeneous elements for both cases. However, the constant value of σ y y attained in each element depends on the implementation approach. The integration point-based scheme computes the exact σ y y at the element centroid, where E is sampled. On the other hand, the nodal-based approach averages nodal temperatures, introducing a source of error when the elastic properties vary in a nonlinear manner. The use of quadratic elements leads to a good agreement with the analytical solution for both schemes, although differences are observed (Figure 5c,d). We plot the error obtained with both schemes for the element Q8 in Figure 6. It is evident that the integration point-based approach reproduces the analytical result more accurately.
Further insight is gained by reproducing the analysis with a single element in the x-direction—see Figure 7. Inspection of Figure 7a reveals substantial differences between Gauss points-based and nodal-based implementations. The former reproduces precisely the analytical stress distribution, with the exact result being obtained at the integration points. Contrarily, the approximation is much poorer when the material gradient is implemented via the temperature. Differences are particularly significant at the left side of the specimen, where the error is on the order of 50%.
Remarkably, negative stresses are predicted for x = 0 . These non-physical compressive stresses arise as a consequence of the particularities of ABAQUS’ criterion for interpolating nodal temperatures, which does not correspond to the linear interpolation outlined in Figure 2. In ABAQUS, the nodal temperature values are multiplied by certain weights, such that the temperature T at an integration point i is given by
T i = j = 1 m T j W i j with   i = 1 , , n .
Here, T j is the temperature in node j, W i j the weight associated with the nodal temperature j and integration point i, and n and m respectively denote the total number of nodes and integration points. The specific values of W i j depend on a number of numerical considerations and can be easily obtained by means of a one-element model. This criterion is motivated by numerical convergence in thermomechanical problems, as it smoothens localized temperature peaks. The resulting variation in the elastic properties within the element is shown in Figure 7b. Differences with a direct linear interpolation are evident. The weighting procedure implemented in ABAQUS brings non-physical values of E, but it shows a better agreement with the material gradation profile at the Gauss integration points. In turn, this better approximation of E ( x ) reduces the error in the computation of the stresses, as quantified in Figure 7c as a function of the number of elements. The log–log plot of Figure 7c shows that the weighted interpolation of ABAQUS exhibits a smaller maximum error in the computation of σ y y at the Gauss points, as well as a faster convergence rate. Consequently, the error intrinsic to a temperature-based graded element is magnified if a standard linear interpolation is used and, therefore, the conclusions of the present study are even more relevant to finite element codes that employ non-weighted interpolations of nodal temperatures. Recall that the integration point-based scheme presented in Section 2.1 captures the analytical solution at the Gauss points with a single element.

3.2. Uniform Traction Perpendicular to the Material Gradient Direction

Consider now the case where the remote load is prescribed as a traction perpendicular to the elastic gradient—see Figure 4c. The Dirichlet boundary conditions of the problem read
u x ( 0 , 0 ) = 0 ,
u y ( x , 0 ) = 0 .
In the case of a plate with infinite height, the only non-zero component of the Cauchy stress tensor is σ y y . Following Refs. [13,22], a membrane resultant N along the x = w / 2 line can be defined as a function of the remote stress σ 0 and the width,
N = σ 0 w .
The compatibility condition 2 ε y y / x 2 = 0 requires the strain component to be of the form
ε y y ( x ) = A x + B ,
and, consequently, one can readily obtain the stress field by considering the exponential elastic modulus variation assumed (12) and making use of Hooke’s law as
σ y y ( x ) = E 0 e β x ( A x + B ) .
The coefficients A and B are obtained by solving
0 w σ y y ( x ) d x = N ,
0 w σ y y ( x ) x d x = N w 2 ,
such that
A = β N 2 E 0 w β 2 e β w 2 β e β w + w β 2 + 2 β e β w β 2 w 2 e 2 β w + 2 e β w 1 ,
B = β N 2 E 0 · e β w [ e β w ( w 2 β 2 + 3 β w 4 ) + w 2 β 2 2 β w + 8 ] β w 4 ( e β w 1 ) ( e β w β 2 w 2 e 2 β w + 2 e β w 1 ) .
The displacement solution can be readily obtained by making use of the strain-displacement relations and applying the boundary conditions (15) and (16),
u x ( x , y ) = ν A 2 x 2 + B x A 2 y 2 ,
u y ( x , y ) = A x + B y .
The analytical solution in the middle line y = h / 2 is compared with the numerical predictions for σ 0 = 2 . Results are shown in Figure 8 for a uniform mesh of 48 plane stress elements.
Differences between the integration point-based scheme and the nodal/temperature-based implementation are particularly significant for the case of linear elements with full integration (Q4, Figure 8a). Sampling the material gradient directly at the Gauss points leads to a good agreement with the analytical solution. However, using temperature-based properties renders the homogeneous element solution. Both the analytical and Gauss point-based solutions show an increasing σ y y along x within those elements close to the left edge. Contrarily, the inverse response is observed when using temperature-dependent properties, as the strain field decreases with x and E is constant element-wise. The Q4 element predicts in all cases a linear variation of σ y y within each element for the quadratic displacement solution under consideration. On the other hand, identical predictions between graded element schemes are obtained when using linear elements with reduced integration (Q4R, Figure 8b). This is unlike the case of a prescribed displacement (see Figure 5b), as the error in the approximation of the strain field ε y y (exact at the integration points only for the Gauss points-based scheme) is compensated. Quadratic elements (Q8 and Q8R) introduce an element-wise variation of E in both approaches and, therefore, differences appear to be smaller than in linear elements—see Figure 8c,d. The error in the approximation is shown in Figure 9 for the Q8 element case. When using full integration, the Gauss point-based approach outperforms the temperature-based, generalized isoparameteric graded element.
Further insight is gained by analysing a very coarse mesh with a single element in the x-direction; results are shown in Figure 10. Regarding the stress (Figure 10a), the prediction obtained from the integration point-based implementation of graded elements agrees well with the analytical solution, being exact at the Gauss points (symbols). On the other hand, the approximation via a nodal/temperature-based approach introduces a significant source of error (larger than 20% in all the integration points). Moreover, when the load is prescribed as a traction, the approximation of the material gradient also influences the strain field, see Equations (8) and (9). As shown in Figure 10b, a better approximation is attained if the material properties are sampled directly at the Gauss points. Both schemes differ at the edges with the analytical solution for a plate of infinite height.

3.3. Uniform Traction Parallel to the Material Gradient Direction

The last case study involves a functionally graded plate subjected to traction in the x-direction, parallel to the elastic gradient, see Figure 4d. Under those conditions, the normal stress component equals the applied stress
σ x x ( x , y ) = σ 0
if Poisson’s ratio is made equal to zero, ν = 0 . The results obtained for a uniform mesh of 75 plane stress elements are shown in Figure 11.
Contrarily to what has been observed so far, the Q4 results (Figure 11a) show that the nodal/temperature-based approach outperforms the integration point-based counterpart. Hooke’s law requires the strains to vary according to an inverse exponential distribution to obtain a constant stress for an exponentially varying E. However, linear elements predict a constant strain field and, consequently, an effectively homogeneous element will predict a constant stress. This trend is inverted for the case of a quadratic element with full integration (Q8). As shown in Figure 11b, again the use of a Gauss point-based graded element formulation approximates the analytical solution better. The results pertaining to reduced integration elements (Q4R and Q8R), not shown for brevity, reveal a perfect agreement with the analytical solution in all cases. Thus, reduced integration improves precision in this specific case study as the resulting stress is constant—either because there is a single integration point (Q4R), leading to a constant ε x x and E, or because both ε x x and E are element-wise linear (Q8R).

4. Conclusions

We have explored the influence of element order, integration scheme and graded element formulation in the finite element analysis of functionally graded materials (FGMs). Two graded element formulations are presented to account for the variation in space of material properties: nodal and integration point based gradations. The nodal based variation is implemented by defining temperature-dependent properties with a zero thermal expansion coefficient, a simple approach that enables the use of this scheme in commercial finite element packages. Important insight is gained by solving, analytically and numerically, three boundary value problems involving remote tractions and displacements, applied parallel and perpendicular to the material gradation direction.
Results reveal that integration point-based graded elements generally outperform a nodal-based implementation through temperature-dependent properties. The former approximates better the analytical solution in all the boundary value problems considered if quadratic shape functions are used. A much finer mesh is needed to attain a similar degree of precision with a nodal-based approach. However, the temperature-based generalized isoparametric graded element is more accurate when linear elements are employed and the traction is applied parallel to the direction of material gradation. These observations are inherent to the interpolation of nodal temperatures with shape functions that are one order lower than those employed for the displacement field, as done in most finite element codes to ensure an equivalent variation of thermal and mechanical strains. In addition, in the case of the commercial package ABAQUS, the nodal temperature averaging criterion employed can lead to non-physical results.
The results presented have implications in the analysis of functionally graded structures, both in terms of computation time and local precision, particularly relevant for fracture studies. A user subroutine for ABAQUS is presented to overcome the number of drawbacks identified with the most popular graded finite element implementation. The user subroutine can be downloaded from www.empaneda.com/codes.

Funding

The author gratefully acknowledges financial support from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007-2013) under REA grant agreement No. 609405 (COFUNDPostdocDTU).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Kawasaki, A.; Watanabe, R. Finite element analysis of thermal stress of the metal/ceramic multi-layer composites with controlled compositional gradients. J. Jpn. Inst. Metals 1987, 51, 525–529. [Google Scholar] [CrossRef]
  2. Uemura, S. The Activities of FGM on New Application. Mater. Sci. Forum 2003, 423–425, 1–10. [Google Scholar] [CrossRef]
  3. Lee, Y.D.; Erdogan, F. Residual/thermal stresses in FGM and laminated thermal barrier coatings. Int. J. Fract. 1994, 69, 145–165. [Google Scholar] [CrossRef]
  4. Ramaswamy, P.; Seetharamu, S.; Varma, K.; Rao, K. Al2O3-ZrO2 Composite coatings for thermal-barrier applications. Compos. Sci. Technol. 1997, 57, 81–89. [Google Scholar] [CrossRef]
  5. Marur, P.R.; Tippur, H.H.V. Evaluation of mechanical properties of functionally graded materials. J. Test. Eval. 1998, 26, 539–545. [Google Scholar]
  6. Rao, B.N.; Rahman, S. Mesh-free analysis of cracks in isotropic functionally graded materials. Eng. Fract. Mech. 2003, 70, 1–27. [Google Scholar] [CrossRef]
  7. Liu, P.; Yu, T.; Bui, T.Q.; Zhang, C. Transient dynamic crack analysis in non-homogeneous functionally graded piezoelectric materials by the X-FEM. Comput. Mater. Sci. 2013, 69, 542–558. [Google Scholar] [CrossRef]
  8. Comi, C.; Mariani, S. Extended finite element simulation of quasi-brittle fracture in functionally graded materials. Comput. Methods Appl. Mech. Eng. 2007, 196, 4013–4026. [Google Scholar] [CrossRef]
  9. Natarajan, S.; Baiz, P.M.; Bordas, S.; Rabczuk, T.; Kerfriden, P. Natural frequencies of cracked functionally graded material plates by the extended finite element method. Compos. Struct. 2011, 93, 3082–3092. [Google Scholar] [CrossRef]
  10. Bao, G.; Wang, L. Multiple cracking in functionally graded ceramic/metal coatings. Int. J. Solids Struct. 1995, 32, 2853–2871. [Google Scholar] [CrossRef]
  11. Rousseau, C.E.; Tippur, H.V. Compositionally graded materials with cracks normal to the elastic gradient. Acta Mater. 2000, 48, 4021–4033. [Google Scholar] [CrossRef]
  12. Santare, M.H.; Lambros, J. Use of Graded Finite Elements to Model the Behavior of Nonhomogeneous Materials. J. Appl. Mech. 2000, 67, 819–822. [Google Scholar] [CrossRef]
  13. Kim, J.H.; Paulino, G.H. Isoparametric Graded Finite Elements for Nonhomogeneous Isotropic and Orthotropic Materials. J. Appl. Mech. 2002, 69, 502. [Google Scholar] [CrossRef]
  14. Valizadeh, N.; Natarajan, S.; Gonzalez-Estrada, O.A.; Rabczuk, T.; Bui, T.Q.; Bordas, S.P.A. NURBS-based finite element analysis of functionally graded plates: Static bending, vibration, buckling and flutter. Compos. Struct. 2013, 99, 309–326. [Google Scholar] [CrossRef] [Green Version]
  15. Martínez-Pañeda, E.; Gallego, R. Numerical analysis of quasi-static fracture in functionally graded materials. Int. J. Mech. Mater. Design 2015, 11, 405–424. [Google Scholar] [CrossRef]
  16. Butcher, R.J.; Rousseau, C.E.; Tippur, H.V. A functionally graded particulate composite: Preparation, measurements and failure analysis. Acta Mater. 1998, 47, 259–268. [Google Scholar] [CrossRef]
  17. Mathew, T.V.; Natarajan, S.; Martínez-Pañeda, E. Size effects in elastic-plastic functionally graded materials. Compos. Struct. 2018, 204, 43–51. [Google Scholar] [CrossRef]
  18. Carrillo-Heian, E.M.; Carpenter, R.D.; Paulino, G.H.; Gibeling, J.C.; Munir, Z.A. Dense Layered Molybdenum Disilicide-Silicon Carbide Functionally Graded Composites Formed by Field-Activated Synthesis. J. Am. Ceram. Soc. 2001, 84, 962–968. [Google Scholar] [CrossRef] [Green Version]
  19. Krumova, M.; Klingshirn, C.; Haupert, F.; Friedrich, K. Microhardness studies on functionally graded polymer composites. Compos. Sci. Technol. 2001, 61, 557–563. [Google Scholar] [CrossRef]
  20. Abanto-Bueno, J.; Lambros, J. Parameters controlling fracture resistance in functionally graded materials under mode I loading. Int. J. Solids Struct. 2006, 43, 3920–3939. [Google Scholar] [CrossRef] [Green Version]
  21. Bower, A.F. Applied Mechanics of Solids; CRC Press, Taylor & Francis: Boca Raton, FL, USA, 2009. [Google Scholar]
  22. Erdogan, F.; Wu, B.H. The surface crack problem for a plate with functionally graded properties. J. Appl. Mech. 1997, 64, 449. [Google Scholar] [CrossRef]
  23. Sayyad, A.S.; Ghugal, Y.M. Modeling and analysis of functionally graded sandwich beams: A review. Mech. Adv. Mater. Struct. 2018. [Google Scholar] [CrossRef]
Figure 1. Sketch outlining the (a) gradual variation of Young’s modulus E along the y-coordinate, as captured by (b) homogeneous elements and (c) graded elements.
Figure 1. Sketch outlining the (a) gradual variation of Young’s modulus E along the y-coordinate, as captured by (b) homogeneous elements and (c) graded elements.
Materials 12 00287 g001
Figure 2. Sketch outlining the gradual variation of Young’s modulus E and its associated interpolation by means of temperature-based generalized isoparametric graded element for an equivalent interpolation of thermal and mechanical strains.
Figure 2. Sketch outlining the gradual variation of Young’s modulus E and its associated interpolation by means of temperature-based generalized isoparametric graded element for an equivalent interpolation of thermal and mechanical strains.
Materials 12 00287 g002
Figure 3. Graded finite elements employed and notation used.
Figure 3. Graded finite elements employed and notation used.
Materials 12 00287 g003
Figure 4. Boundary value problems under consideration: (a) functionally graded plate with spatially varying Young’s modulus subjected to (b) uniform displacement perpendicular to the material gradient direction, (c) uniform traction perpendicular to the material gradient direction, and (d) uniform traction in the direction parallel to material gradation—consistent units.
Figure 4. Boundary value problems under consideration: (a) functionally graded plate with spatially varying Young’s modulus subjected to (b) uniform displacement perpendicular to the material gradient direction, (c) uniform traction perpendicular to the material gradient direction, and (d) uniform traction in the direction parallel to material gradation—consistent units.
Materials 12 00287 g004
Figure 5. Uniform displacement perpendicular to the material gradient for different kinds of elements: (a) Q4; (b) Q4R; (c) Q8; and (d) Q8R—consistent units.
Figure 5. Uniform displacement perpendicular to the material gradient for different kinds of elements: (a) Q4; (b) Q4R; (c) Q8; and (d) Q8R—consistent units.
Materials 12 00287 g005
Figure 6. Uniform displacement perpendicular to the material gradient. Error analysis for the Q8 element. Consistent units.
Figure 6. Uniform displacement perpendicular to the material gradient. Error analysis for the Q8 element. Consistent units.
Materials 12 00287 g006
Figure 7. Uniform displacement perpendicular to the material gradient, (a) tensile stress for one Q8 element in the x-direction; (b) Young’s modulus interpolation through different schemes; and (c) mesh-sensitivity error analysis—consistent units.
Figure 7. Uniform displacement perpendicular to the material gradient, (a) tensile stress for one Q8 element in the x-direction; (b) Young’s modulus interpolation through different schemes; and (c) mesh-sensitivity error analysis—consistent units.
Materials 12 00287 g007
Figure 8. Uniform traction perpendicular to the material gradient for different kinds of elements: (a) Q4; (b) Q4R; (c) Q8; and (d) Q8R—consistent units.
Figure 8. Uniform traction perpendicular to the material gradient for different kinds of elements: (a) Q4; (b) Q4R; (c) Q8; and (d) Q8R—consistent units.
Materials 12 00287 g008
Figure 9. Uniform traction perpendicular to the material gradient. Error analysis for the Q8 element—consistent units.
Figure 9. Uniform traction perpendicular to the material gradient. Error analysis for the Q8 element—consistent units.
Materials 12 00287 g009
Figure 10. Uniform traction perpendicular to the material gradient, tensile (a) stress and (b) strain for one Q8 element in the x-direction—consistent units.
Figure 10. Uniform traction perpendicular to the material gradient, tensile (a) stress and (b) strain for one Q8 element in the x-direction—consistent units.
Materials 12 00287 g010
Figure 11. Uniform traction parallel to the material gradient for different kinds of elements: (a) Q4 and (b) Q8—consistent units.
Figure 11. Uniform traction parallel to the material gradient for different kinds of elements: (a) Q4 and (b) Q8—consistent units.
Materials 12 00287 g011

Share and Cite

MDPI and ACS Style

Martínez-Pañeda, E. On the Finite Element Implementation of Functionally Graded Materials. Materials 2019, 12, 287. https://doi.org/10.3390/ma12020287

AMA Style

Martínez-Pañeda E. On the Finite Element Implementation of Functionally Graded Materials. Materials. 2019; 12(2):287. https://doi.org/10.3390/ma12020287

Chicago/Turabian Style

Martínez-Pañeda, Emilio. 2019. "On the Finite Element Implementation of Functionally Graded Materials" Materials 12, no. 2: 287. https://doi.org/10.3390/ma12020287

APA Style

Martínez-Pañeda, E. (2019). On the Finite Element Implementation of Functionally Graded Materials. Materials, 12(2), 287. https://doi.org/10.3390/ma12020287

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop