Document Type : Original Article


1 Faculty of Aerospace Engineering, Amirkabir University of Technology -

2 Aerospace Research Institute, Ministry of Science, Research and Technology

3 - Faculty of Aerospace Engineering, Amirkabir University of Technology


This paper investigates the micromechanical damage behavior of composite materials on the transverse tensile response of representative volume elements (RVEs) with various fiber volume fractions. A new algorithm is developed to generate the random distribution of fibers in the RVE, and it is possible to create a fiber distribution with high fiber volume fraction. The fibers and matrix are considered linear elastic and elastoplastic with Drucker–Prager's criterion, respectively. The fiber-matrix debonding is modeled by cohesive elements which lead to matrix cracking. To achieve a more realistic microstructure, the normal cohesive properties distribution is applied for interfaces between fibers and the matrix. To investigate the effect of the position of fibers with the weakest cohesive strength, sensitivity analysis concerning the different arrangements of specific normal cohesive properties on the RVE’s strength are performed. Moreover, the effects of different damage parameters such as fiber random distributions and various cohesive strengths on the overall damage behavior of the RVE are described in detail. It is revealed that the application of a normal cohesive distribution strongly reduces the maximum strength of the RVE and shifts the strain of damage initiation point and crack propagation path.


[1] Gamstedt, E. K., and B. A. Sjögren. "Micromechanisms in tension-compression fatigue of composite laminates containing transverse plies." Composites Science and Technology 59.2 (1999): 167-178. [2] Hobbiebrunken T, Hojo M, Adachi T, De Jong C, Fiedler B. Evaluation of interfacial strength in CF/epoxies using FEM and in-situ experiments. Compos Part A – Appl Sci 2006;37(12):2248–56.
[3] Totry E, Gonzlez C, Llorca J. Failure locus of fiber[1]reinforced composites under transverse compression and out-of-plane shear. Compos Sci Technol 2008;68(3–4):829–39.
[4] Hojo M, Mizuno M, Hobbiebrunken T, Adachi T, Tanaka M, Ha SK. Effect of fiber array irregularities on microscopic interfacial normal stress states of transversely loaded UD-CFRP from viewpoint of failure initiation. Compos Sci Technol 2009;69(11– 12):1726–34.
[5] Asp LE, Berglund LA, Talreja R. Prediction of matrix[1]initiated transverse failure in polymer composites. Compos Sci Technol 1996;56(9):1089–97.
[6] Maligno AR, Warrior NA, Long AC. Finite element investigations on the microstructure of fibre-reinforced composites. Express Polym Lett 2008;2(9):665–76.
[7] Trias D, Costa J, Mayugo JA, Hurtado JE. Random models versus periodic models for fibre reinforced composites. Comput Mater Sci 2006;38(2):316–24.
[8] Vaughan, T. J., & McCarthy, C. T. (2010). A combined experimental–numerical approach for generating statistically equivalent fibre distributions for high strength laminated composite materials. Composites Science and Technology, 70(2), 291-297. doi:
[9] Ismail, Yaser, Dongmin Yang, and Jianqiao Ye. "Discrete element method for generating random fibre distributions in micromechanical models of fibre reinforced composite laminates." Composites Part B: Engineering 90 (2016): 485-492.
[10] Gonzlez C, Llorca J. Mechanical behavior of unidirectional fiber-reinforced polymers under transverse compression: microscopic mechanisms and modeling. Compos Sci Technol 2007;67(13):2795–806.
[11] Alfaro, MV Cid, et al. "Numerical homogenization of cracking processes in thin fibre-epoxy layers." European Journal of Mechanics-A/Solids 29.2 (2010): 119-131.
[12] Naderi, M., N. Apetre, and N. Iyyer. "Effect of interface properties on transverse tensile response of fiber[1]reinforced composites: Three-dimensional micromechanical modeling." Journal of Composite Materials 51.21 (2017): 2963-2977.
[13] Trias D, Costa J, Turon A, Hurtado J. Determination of the critical size of a statistical representative volume element (SRVE) for carbon reinforced polymers. Acta Mater 2006;54(13):3471e84.
[14] Yang L, Yan Y, Ran Z, Liu Y. A new method for generating random fibre distributions for fibre reinforced composites. Compos Sci Technol 2013;76: 14e20.
[15] Melro A, Camanho P, Pinho S. Generation of random distribution of fibres in long-fibre reinforced composites. Compos Sci Technol 2008;68(9):2092e102.
[16] Xia, Z., Zhang, Y., & Ellyin, F. (2003). A unified periodical boundary conditions for representative volume elements of composites and applications. International Journal of Solids and Structures, 40(8), 1907-1921.
[17] Nguyen, V. D., Béchet, E., Geuzaine, C., & Noels, L. (2012). Imposing periodic boundary condition on arbitrary meshes by polynomial interpolation. Computational Materials Science, 55, 390-406.
[18] Barbero, E. J. (2013). Finite element analysis of composite materials using AbaqusTM. CRC press.
[19] Melro, A. R. D. O. S. (2011). Analytical and numerical modelling of damage and fracture of advanced composites.
[20] Soden, P. D., M. J. Hinton, and A. S. Kaddour. "Lamina properties, lay-up configurations and loading conditions for a range of fibre-reinforced composite laminates." Composites Science and Technology 58.7 (1998): 1011-1022.
[21] Vaughan TJ and McCarthy CT. Micromechanical modelling of the transverse damage behaviour in fibre reinforced composites. Compos Sci Technol 2011; 71: 388–396.
[22] Asp LE, Berglund LA and Talreja R. A criterion for crack initiation in glassy polymers subjected to a composite-like stress state. Compos Sci Technol 1996; 56: 1291–1301.
[23] Asp LE, Berglund LA and Gudmundson P. Effects of a composite-like stress state on the fracture of epoxies. Compos Sci Technol 1995; 53: 27–37.
[24] Yang, Lei, et al. "Effects of triangle-shape fiber on the transverse mechanical properties of unidirectional carbon fiber reinforced plastics." Composite Structures 152 (2016): 617-625