Numerical Investigation on the Fracture Behaviors of Three-Dimensional Functionally Graded Materials

Abstract:

Functionally graded materials (FGMs) with continuous varying properties have absorbed great attention for the purpose of eliminating the mismatch of material properties which may result in cracking. In this paper, three-dimensional finite element method (3D FEM) based on nonhomogeneous elements is used to study the fracture behaviors of a 3D FGM plate. Since real material properties at Gaussian integration points are adopted during forming the element stiffness matrix, the nonhomogeneous material properties can be applied in each element. Moreover, 20-node singular elements are used around the crack front to deal with the singularity of stress fields at the crack front. By this way, the stress intensity factors (SIFs) can be calculated with high efficiency and accuracy. Therefore, compared with the general FEM using homogeneouos elements, the calculating efficiency and accuracy can be increased. Finally, parameter analysis is conducted. It is found that the material nonhomogeneity constant and the crack parameter have significant influences on the SIFs.

• Periodical: Key Engineering Materials (Volumes 353-358)
• Pages: 1098-1101
• DOI: https://doi.org/10.4028/www.scientific.net/KEM.353-358.1098
• Online since: September 2007
• Authors: Hong Jun Yu, Li Cheng Guo, Lin Zhi Wu
• Keywords: Crack, Functionally Graded Material (FGM), Nonhomogeneous Elements, Stress Intensity Factor (SIF), Three-Dimensional

References

[1] F. Delale, F. Erdogan: J. Appl. Mech. Vol. 50 (1983), p.609.

[2] F. Erdogan: J. Appl. Mech. Vol. 52 (1985), p.823.

[3] J. H. Kim, G. H. Paulino: Int. J. Numer. Meth. Engng. Vol. 53 (2002), p. (1903).

[4] J. H. Kim, G. H. Paulino: Eng. Fracture Mech. Vol. 69 (2002), p.1557.

[5] M. C. Walters, G. H. Paulino, Robert H: Int. J. Solids Struct. Vol. 41 (2004), p.1081.

[6] B. Yildirim, S. Dag and F. Erdogan: Int. J. Fracture. Vol. 132 (2005), p.369.

[7] Z. Z. Liang, C. A. Tang and D. S. Zhao: Key Eng. Mater. Vol. 297-300 (2005), p.1113.

[8] R. W. J. Koers: Int. J. Fracture. Vol. 40 (1989), p.79.

[9] S. Kadioglu, S. Dag, S. Yahsi. Int. J. Fracture. Vol. 94 (1998), p.63.

[10] L. C. Guo, L. Z. Wu, T. Zeng, L. Ma: Eur. J. Mech. A/Solids. Vol. 23 (2004), p.219.

[11] T. Nakamura, D.M. Parks: J. Appl. Mech. Vol. 55 (1988), p.805.

0 0. 2 0. 4 0. 6 0. 8 1. 0.

[1] 0.

[1] 1.

[1] 2.

[1] 3 a0=0. 05 Normalized SIFs at crack front-a z/(t/2) E2/E1=.

1.

2.

5 1 2 5 10.

0 0. 2 0. 4 0. 6 0. 8 1. 0.

[1] 5.

[1] 8.

[2] 1.

[2] 4.

[2] 7.

[3] 0 a0=0. 4 E2/E1=.

1.

2.

5 1 2 5 10 Normalized SIFs at crack front-a z/(t/2) Fig. 4 the SIFs along crack front-a in a 3D plate (a) 0 0. 05a = ; (b) 0 0. 4a = (a) (b).