Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 8757132, 9 pages
https://doi.org/10.1155/2018/8757132

Research Article

The Displacement Analysis Investigation for Variable Gradient Parameters Functionally Graded Materials Based on Strip Element Method

J. H. Tian and L. L. Ma

CAE Analysis Room for Engineering Application, School of Mechatronic Engineering, Xi’an Technological University, Xi’an 710021, China

Correspondence should be addressed to J. H. Tian; 540673737@qq.com

Received 26 January 2018; Revised 22 March 2018; Accepted 5 April 2018; Published 24 May 2018

Academic Editor: Francesco Tornabene

Copyright © 2018 J. H. Tian and L. L. Ma. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper adopted the strip element method to research the displacement response of functionally graded material structure with variable gradient parameters under mechanical loading. Based on the dynamic equilibrium equation of the system, strip element method theory was established through the principle of virtual work to get the control equations of functionally graded materials plane problem, using the modal superposition method to complete the solution of the equation and obtain the displacement distribution of the FGM plates with variable gradient parameters. The accuracy of the calculation method was verified by comparison with the result of finite element method. The results show that the displacement response of the strip element method is consistent with the finite element method. With the increase of the gradient coefficient, the ceramic content is increased and the material stiffness is also increased; then the displacement response caused by external loading is decreased. With the increase of the gradient coefficient, the influence of the gradient coefficient on the displacement response is decreased gradually.

1. Introduction

Functionally graded materials (FGM) [1] are a newly developed material in recent years. Their distinguishing feature is that the composition of the material presents a continuous gradient change in the direction of thickness. The feature and function of the material also appear as a gradient change.

A large number of domestic and abroad scholars had done a lot of research on the mechanical problems of functionally graded materials [2, 3]. Tornabene [4] applied the Generalized Differential Quadrature (GDQ) Method to study four parameters functionally graded and laminated composite shells and panels of revolution. The free vibration of functional gradient materials under mechanical loading and temperature loading was analyzed using Peano-Baker series method by Liu [5]. The dynamic characteristic analytical solution of the quadrilateral and rhombic function gradient plate was obtained by using the beam function group [6]. Tornabene [7] analyzed the dynamic of doubly curved shell structures, and the natural frequencies and mode shapes of several structures are presented and compared to numerical solutions taken from the literature. Viola [8] studied the dynamic behavior of moderately thick functionally graded parabolic panels of revolution and was concerned with the influence of the three parameters of the power-law distribution on their mechanical behavior. Reddy [9] adopted the first-order shear deformation theory to study the structural response of the functional gradient plate structure under the thermomachine coupling. Based on the first-order shear deformation theory (FSDT), Tornabene [10] focused on the dynamic behavior of moderately thick functionally graded parabolic panels and shells of revolution and illustrated the influence of the parameters of the power-law distribution on the mechanical behavior of shell structures. And furthermore the study concerned six types of shell structures that illustrate the influence of the power-law exponent, of the power-law distribution and of the choice of the four parameters on the mechanical behavior of shell structures considered [11]. Using the first-order shear deformation theory, Thai [12] had studied the problem of bending and free vibration of functionally graded materials. Liew [13] studied the vibration response of functionally graded material plates under the coupling of mechanical load and electromagnetic load based on the first-order shear deformation theory using the finite element method. Zenkour [14] used first-order shear deformation theory and classical elastic theory to study the structural bending deformation of the sandwich structure under mechanical loading. Ferreira [15] used three-order shear deformation theory and mesh-less method to study the structure response of functional gradient materials under static load. Reddy [16] studied the deformation of functional gradient rectangular plate based on the three-order shear deformation theory. Brischetto [17] studied the analytic solution of the structural displacement response of functionally graded material plates under transverse loading. Ray [18] used the finite element method to study the structural response of the functional gradient fiber piezoelectric materials in exponential distribution. Kulikov [19] adopted surface sampling method to solve the problem of three-dimensional thermal stress of functional gradient materials.

Due to physical performance parameters of functionally gradient materials with the coordinate changing continuously, the constant coefficient differential equation turns into the variable coefficient differential equation while the theory of elastic body deformation is derived; it is difficult to solve such problem by analytical method. The numerical method is used to solve complex problem while the large computation is in process. The analytical solution is very difficult to obtain, while the precision of the solution will be declined for numerical method for FGM structure, and it also requires a lot of computational memories while using numerical method. In recent years, the researchers have proposed a seminumerical and semianalytic method, that is called the strip element method (SEM) [20], which combined the accuracy of the analytical method with the efficiency of the numerical method.

This study considered the boundary and initially applied strip element method to functionally graded materials structure. Meanwhile this paper firstly considered the displacement response of FGM structure with variable gradient parameters. The paper obtained the displacement response of functionally graded material structure under mechanical loading and further investigated the relationship between the gradient parameters and displacement response of functionally graded material structure. It provides guidance for the application of FGM.

2. Strip Element Method Theory

The strip element method is a series of strip elements separated by the solution domain as shown in Figure 1. The displacement function is the continuous function of direction, and the displacement of the element is obtained through the line displacement interpolation.

Figure 1: Strip element method for solving model.

Based on the finite element method in direction, it is assumed that the displacement functions of functionally graded materials are

where represents matrix of shape functions for displacement interpolation; represents continuous function of direction; represents a complex unit; represents circle frequency; represents the computing time.

The shape functions can be derived from the displacement function of a single element. The matrix of shape functions [20] is

where represents element thickness; represents the second order identity matrix.

According to Kausel’s equation, the dynamic equilibrium equation of the system is

where represents displacement vector; represent differential operator matrix; represents stress vector and ; superscript denotes transpose,

where the expression for is

It can be seen that the internal displacement of the element is obtained by the interpolation of the nodal displacement from (1). The residual value must be generated in the interpolation calculation, and the residual value is expressed as follows:

Since the internal displacement function of the element is interpolated by the nodal displacement, the residual value must exist in the element; then .

According to the principle of virtual work, for a deformed body in equilibrium state, the total virtual work done by the external force is equal to the sum of the work done by the internal force and the work done by the residual value, as follows:

where is tiny virtual displacement, represents external mechanical load, and represents nodal stress vector.

Given an external load as

where represents external load amplitude, the element stress is

where stress is

where ; is the elastic constant matrix; .

The displacement function is substituted into the virtual work equation based on the theory of thin plate deformation and we obtain

where , , and represent coefficient matrix. represents mass matrix.

It can be seen from (10) that the governing equation obtained by element discretization is the ordinary differential equation of , which reduces the difficulty of the solution; it also reflects the idea of direction discretization.

The coefficient matrix is

where

The expression for the mass matrix is

Assembling all elements, (10) becomes as

where , , and are the coefficients matrix for the element assembly in the solving domain; is the mass assembly matrix; is the external mechanical load; is the displacement function in the solving domain.

Assume that the displacement function form is

Assume that external mechanical load form is

Substituting (17) and (18) into (16), we obtain

represents the load acting on the nodal line, while (19) is transformed into a homogeneous equation and the eigenvalue equation for the k is obtained.

Equation (20) can be written

The 2M () eigenvalue is obtained by solving (21), is the element number, and the eigenvectors corresponding to the eigenvalues are represented by

By means of the modal superposition method, the solution of the equation can be expressed as

where

Equation (23) is a fundamental solution of the equation. There are 2M constants in this fundamental solution, so it is necessary to determine the unknown coefficients by the boundary conditions.

Since (23) is the fundamental solution in the problem domain, the displacement at any point satisfies (23) and the displacement is satisfied at the right boundary of the solution.

That can be written as

In the same way, on the left side

Equations (27) and (28) are assembled and sorted, and the constant expression can be obtained

where is the displacement at the boundary

By (24), can be expressed as

The superscript represents the right boundary.

Similarly, the left boundary condition has the same form as the displacement matrix at the right boundary. The difference between the two forms is that the coordinate value at the right boundary is changed into the coordinate value at the left boundary.

3. Application of the Boundary Conditions

The stress vector of the internal element can be expressed as

Substituting (1) into (34), we have

where and are element coefficient matrices.

is

is

Calculate and for all elements, and obtain the overall stress vector as follows:

where and are the assemble matrix of the element coefficient matrix.

The displacement is simplified by using the left and right boundaries

The stress vector of the inner nodal line is equal to the average value of the stress vector of the adjacent left and right elements. We obtain

In (40), is an external load acting on the left and right boundaries, where the stiffness matrix is

In (41), and have the same expressions as and . The difference is that and in and are replaced by and , respectively.

Equation (40) reflects the relationship between the stress boundary condition and displacement; the stress and displacement boundary conditions are transformed into displacement boundary conditions by (40); a constant can be obtained by substituting the boundary displacement into (29); from (23) and (1) the displacement can be obtained in the solving problem.

4. Establishment of Gradient Parameter Model

The functional gradient materials studied in this paper and the model is composite by metal and ceramic materials. The metal material is 1Cr₁₈Ni₁₁Nb and the ceramic material is Si₃N₄. The physical properties of the two materials are shown in Table 1.

Table 1: Physical properties of 1Cr₁₈Ni₁₁Nb and Si₃N₄ materials.

material
Modulus of elasticity/MPa
Poisson’s ratio

1Cr₁₈Ni₁₁Nb
2.38 × 10⁵
0.3177

Si₃N₄
3.22 × 10⁵
0.2400

Because functionally graded materials are continuously changed by different materials according to the design requirements, the physical performance parameters of the materials are expressed as a function of the volume fraction of the material, the physical properties of the materials, and the content of components; it denotes

where , , and are the elastic modulus, Poisson’s ratio, and density of FGM plates, respectively. , , and are the elastic modulus, Poisson’s ratio, and density of ceramic materials, respectively. , , and are the elastic modulus, Poisson’s ratio, and density of metal materials, respectively. is the volume fraction of the ceramic.

To ensure the continuity of the material, the sum of the metal volume fraction and the ceramic volume fraction at any location of the material is 1.

where is the volume fraction of metals.

The volume fraction change function of the functionally graded metal is

where is the vertical position and ; is FGM plate thickness; is the gradient parameters.

By (44), it can be seen that when , and . Substituting that into (42) we have , , and and functionally graded material degenerates into pure metal homogeneous material. When p tends to infinity, tends to 1 and (42) can be obtained , , and and the material is reduced to a pure ceramic homogeneous material.

When and are functionally graded materials, the volume fraction of metals varies as shown in Figure 2. Figure 2 can be shown: when , metal volume fraction is 0 and this position is pure ceramic material; when , metal volume fraction is equal to 1 and this position is the pure metal homogeneous material. So when and , the transition form of materials is a continuous transition of ceramic to metal. When the gradient parameter , the change of the volume fraction of the functionally graded material is continuous linear change. When the gradient parameter is lesser than 1, the change of the metal volume fraction decreases. When the gradient parameter is greater than 1, the metal volume fraction of the functional gradient material is slower in the early and then the volume fraction is accelerated with the thickness increased.

Figure 2: Law of change of volume fraction of metal with thickness under different gradient parameters at and .

5. Computational Model

5.1. Validation of Strip Element Method

Based on the finite element (FEM) theory, we used ANSYS to verify the results of SEM solution. The model of Figure 3 is solved by using the strip element method and finite element method. In Figure 3, = 100 mm, = 100 mm, the lower boundary is fixed, = 100 N/mm, and uniform load is applied on the upper boundary. The solution domain is divided into 20 elements and the material properties in each element are metal material. The metallic material property is shown in Table 1. The displacement of -direction and -direction at = 10 mm, = 50 mm, and = 90 mm is shown in Figures 4 and 5 by strip element method and finite element method (FEM).

Figure 3: Physical model of plane problem.

Figure 4: The displacement in -direction in different position by FEM and SEM.

Figure 5: The displacement in -direction in different position by FEM and SEM.

From Figures 4 and 5 we can see that the displacement obtained by strip element method and finite element method in -direction and -direction at different position is almost the same. From = 10 mm and = 50 mm to = 90 mm, the displacement in -direction is increased. The accuracy of the strip element method is verified.

5.2. The Influence of Element Number for the Results

In the calculation process, because of the discrete of the solving domain, the size of the element has an influence on the accuracy of the displacement solution. Figures 6 and 7 are the -direction displacement and the -direction displacement at = 10 mm using 10 elements, 20 elements, and 50 elements, respectively. It can be known that displacement results have small difference when the element numbers are 10 and 20, while when the element numbers are 20 and 50, the result difference is not obvious, so in this calculation the element number is set to 20 for ensuring that precision and calculation speed.

Figure 6: Calculation results of direction displacement at = 10 mm.

Figure 7: Calculation results of direction displacement at = 10 mm.

5.3. Structural Response Analysis of Functionally Graded Materials

The displacements with different gradient parameters are shown in Figure 8 at = 10 mm in the direction. It can be seen from Figure 8, with the increase of gradient coefficient, that the displacement response decreases, and this is due to the fact that metal volume fraction decreases in the functionally graded materials and ceramic volume fraction increases, causing the structure stiffness to increase. In addition, the displacement in the direction at = 10 mm has nonlinear characteristics.

Figure 8: direction displacement of different gradient parameters at = 10 mm.

The displacements with different gradient parameters are shown in Figure 9 at = 10 mm in the direction. It can be seen from Figure 9 that the displacement in direction at = 10 mm with the increase of coordinate has the characteristics of first increasing, then decreasing, and then increasing. As the gradient parameter increases, the displacement of the same position decreases. As the gradient coefficient increases, the influence of the gradient coefficient for the displacement response decreases gradually.

Figure 9: directional displacement of different gradient parameters at = 10 mm.

direction of displacement is shown in Figures 10 and 11 at = 50 mm and = 90 mm. They have the same law at = 50 mm and = 90 mm. The displacement absolute value of direction with respect to = 50 mm is symmetrical and, with the increase of gradient coefficient, direction displacement is reduced; this is due to the increase of the gradient coefficient; the ceramics content of the material keeps increasing causing the structure stiffness to increase.

Figure 10: directional displacement of different gradient parameters at = 50 mm.

Figure 11: direction displacement for different gradient parameters at = 90 mm.

direction of displacement is shown in Figures 12 and 13 at = 50 mm and = 90 mm. It can be seen that the displacement in direction has class parabolic distribution, and as the change of spatial coordinates, the curve is approximatively changed to straight with the increase of gradient coefficient, structural stiffness is gradually increased and the displacement response under external loading is reduced.

Figure 12: direction displacement of different gradient parameters at = 50 mm.

Figure 13: direction displacement for different gradient parameters at = 90 mm.

6. Conclusion

This paper presents the variable gradient parameters model of the functional gradient materials. The displacement response of considering different gradient parameters model is analyzed in mechanics loading using the strip element method; there are some useful results.

(1) The displacement response result is almost the same by using finite element method and strip element method, so as to verify the validation of the strip element method.

(2) With the ceramic volume fraction increasing, the metal volume fraction decreases in the functionally graded materials, the structure stiffness increases, and displacement response decreases.

(3) As the gradient coefficient increases, the influence of the gradient coefficient on the structure displacement response decreases gradually.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

This work was presented at “the 8th International Conference on Computational Methods (ICCM2017)” for communication.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The present paper was financially supported by the National Natural Science Foundation of China (no. 11302159), Young Talent Fund of University Association for Science and Technology in Shaanxi (no. 20170519), and the Principal Foundation of Xi’an Technological University (no. XAGDXJJ16009).

References

M. Koizumi, “FGM activities in Japan,” Composites Part B: Engineering, vol. 28, no. 1-2, pp. 1–4, 1997.

Z. Zhen, W. Linzhi, and C. Weiqiu, “Progress in the Study on Mechanics Problems of Functionally Graded Materials and Structures,” Advances in Mechanics, vol. 40, no. 5, pp. 528–541, 2010.

Z. Cao, “Unified expression of natural frequency solutions for functionally graded composite rectangular plates under various boundary conditions,” Fuhe Cailiao Xuebao/Acta Materiae Compositae Sinica, vol. 22, no. 5, pp. 172–177, 2005.

F. Tornabene and E. Viola, “Static analysis of functionally graded doubly-curved shells and panels of revolution,” Meccanica, vol. 48, no. 4, pp. 901–930, 2013.

W. X. Liu, “Three-dimensional thermoelastic analysis on functionally graded plates with arbitrary graded distribution functions,” Tongji Daxue Xuebao. Ziran Kexue Ban. Journal of Tongji University. Natural Science, vol. 40, no. 2, pp. 205–210, 2012.

C. A. O. Zhiyuan, “Analytical Solution of Dynamic Character for Functionally Graded Material Parallelogrammic Plates with Various Boundary Conditions,” Chinese Quarterly of Mechanics, vol. 27, no. 2, pp. 255–261, 2006.

F. Tornabene, N. Fantuzzi, and M. Bacciocchi, “Free vibrations of free-form doubly-curved shells made of functionally graded materials using higher-order equivalent single layer theories,” Composites Part B: Engineering, vol. 67, pp. 490–509, 2014.

E. Viola and F. Tornabene, “Free vibrations of three parameter functionally graded parabolic panels of revolution,” Mechanics Research Communications, vol. 36, no. 5, pp. 587–594, 2009.

J. N. Reddy, “Analysis of functionally graded plates,” International Journal for Numerical Methods in Engineering, vol. 47, no. 1–3, pp. 663–684, 2000.

F. Tornabene and E. Viola, “Free vibrations of four-parameter functionally graded parabolic panels and shells of revolution,” European Journal of Mechanics - A/Solids, vol. 28, no. 5, pp. 991–1013, 2009.

F. Tornabene, “Free vibration analysis of functionally graded conical, cylindrical shell and annular plate structures with a four-parameter power-law distribution,” Computer Methods Applied Mechanics and Engineering, vol. 198, no. 37-40, pp. 2911–2935, 2009.

H.-T. Thai and D.-H. Choi, “A simple first-order shear deformation theory for the bending and free vibration analysis of functionally graded plates,” Composite Structures, vol. 101, pp. 332–340, 2013.

K. M. Liew, X. Q. He, and S. Kitipornchai, “Finite element method for the feedback control of FGM shells in the frequency domain via piezoelectric sensors and actuators,” Computer Methods Applied Mechanics and Engineering, vol. 193, no. 3-5, pp. 257–273, 2004.

A. M. Zenkour, “A comprehensive analysis of functionally graded sandwich plates: part 1-deflection and stresses,” International Journal of Solids and Structures, vol. 42, no. 18-19, pp. 5224–5242, 2005.

A. J. M. Ferreira, R. C. Batra, C. M. C. Roque, L. F. Qian, and P. A. L. S. Martins, “Static analysis of functionally graded plates using third-order shear deformation theory and a meshless method,” Composite Structures, vol. 69, no. 4, pp. 449–457, 2005.

J. N. Reddy and C. D. Chin, “Thermomechanical analysis of functionally graded cylinders and plates,” Journal of Thermal Stresses, vol. 21, no. 6, pp. 593–626, 1998.

S. Brischetto and E. Carrera, “Advanced mixed theories for bending analysis of functionally graded plates,” Computers & Structures, vol. 88, no. 23-24, pp. 1474–1483, 2010.

M. C. Ray and H. M. Sachade, “Finite element analysis of smart functionally graded plates,” International Journal of Solids and Structures, vol. 43, no. 18-19, pp. 5468–5484, 2006.

G. M. Kulikov and S. V. Plotnikova, “A sampling surfaces method and its implementation for 3D thermal stress analysis of functionally graded plates,” Composite Structures, vol. 120, pp. 315–325, 2015.

G. R. Liu and J. D. Achenbach, “A strip element method for stress analysis of anisotropic linearly elastic solids,” American Society of Mechanical Engineers. Transactions of the ASME. Journal of Applied Mechanics, vol. 61, no. 2, pp. 270–277, 1994.