FEM Implementation of the Coupled Elastoplastic/Damage Model: Failure Prediction of Fiber Reinforced Polymers (FRPs) Composites

Document Type : Research Paper


1 Laboratoire des Technologies Innovantes, LTI-EA 3899, Université de Picardie Jules Verne, Amiens, France

2 School of Engineering, Edith Cowan University, Joondalup, Perth, Australia



The coupled damage/plasticity model for meso-level which is ply-level in case of Uni-Directional (UD) Fiber Reinforced Polymers (FRPs) is implemented. The mathematical formulations, particularly the plasticity part, are discussed in a comprehensive manner. The plastic potential is defined in effective stress space and the damage evolution is based on the theory of irreversible thermodynamics. The model which is illustrated here has been implemented by different authors previously but, the complete pre-requisite algorithm ingredients used in the implicit scheme implementation are not available in the literature. This leads to the complexity in its implementation. Furthermore, this model is not available as a built-in material constitutive law in the commercial Finite Element Method (FEM) softwares. To facilitate the implementation and understanding, all the mathematical formulations are presented in great detail. In addition, the elastoplastic consistent operator needed for implementation in the implicit solution scheme is also derived. The model is formularized in incremental form to be used in the Return Mapping Algorithm (RMA). The quasi-static load carrying capability and non-linearity caused by the collaborative effect of damage and plasticity is predicted with User MATerial (UMAT) subroutine which solves the FEM problem with implicit techniques in ABAQUS.


[1] Ladevèze PLD E., 1992, Damage modelling of the elementary ply for laminated composites, Composites Science and Technology 43: 257-267.
[2] Liu PFZ J.Y., 2010, Recent developments on damage modeling and finite element analysis for composite laminates: A review, Materials & Design 31(8): 3825-3834.
[3] Wang Y., Tong M., Zhu S., 2009, Three dimensional continuum damage mechanics model of progressive failure analysis in fibre-reinforced composite laminates, Paper Presented at the Structures, Structural Dynamics, and Materials Conference, California.
[4] Edgren F., Mattsson D., Asp L.E., Varna J., 2004, Formation of damage and its effects on non-crimp fabric reinforced composites loaded in tension, Composites Science and Technology 64(5): 675-692.
[5] Herakovich C., Schroedter R., Gasser A., Guitard L., 2000, Damage evolution in [±45]s laminates with fiber rotation, Composite Science and Technology 60: 2781-2789.
[6] Abu Al-Rub R.K.V., George Z., 2003, On the coupling of anisotropic damage and plasticity models for ductile materials, International Journal of Solids and Structures 40(11): 2611-2643.
[7] Maire J., Chaboche J., 1997, A new formulation of continuum damage mechanics for composite materials, Aeropace Science and Technology 4: 247-257.
[8] O’Higgins R.M., McCarthy C.T., McCarthy M.A., 2011, Identification of damage and plasticity parameters for continuum damage mechanics modelling of carbon and glass fibre-reinforced composite materials, Strain 47(1): 105-115.
[9] Payan J., Hochard C., 2002, Damage modelling of laminated carbon/epoxy composites under static and fatigue loadings, International Journal of Fatigue 24: 299-306.
[10] Greve L., Pickett A.K., 2006, Modelling damage and failure in carbon/epoxy non-crimp fabric composites including effects of fabric pre-shear, Composites Part A: Applied Science and Manufacturing 37(11): 1983-2001.
[11] Boutaous A., Peseux B., Gornet L., Bélaidi A., 2006, A new modeling of plasticity coupled with the damage and identification for carbon fibre composite laminates, Composite Structures 74(1): 1-9.
[12] Kachanov L.M., 1958, Time of the rupture process under creep condition, Jzvestia Akad Nauk 8: 26-31.
[13] Chaboche J., 1988, Continuum damage mechanics: Part 1-General concepts, Journal of Applied Mechanics 55: 59-64.
[14] Ribeiro M.L., Tita V., Vandepitte D., 2012, A new damage model for composite laminates, Composite Structures 94(2): 635-642.
[15] Malgioglio F., Payan D., Magneville B., Farkas L., 2017, Material parameter identification challenge and procedure for intra-laminar damage prediction in unidirectional CFRP, Paper Presented at the 6th ECCOMAS Thematic Conference on the Mechanical Response of Composites.
[16] Chen J.F., Morozov E.V., Shankar K., 2012, A combined elastoplastic damage model for progressive failure analysis of composite materials and structures, Composite Structures 94(12): 3478-3489.
[17] Neto E.A.S., Peric D., Owen D.R.J., 2008, Computational Methods for Plasticity: Theory and Applications, John Wiley and Sons, West Sussex, United Kingdom.
[18] Ud Din I., Hao P., Franz G., Panier S., 2018, Elastoplastic CDM model based on Puck’s theory for the prediction of mechanical behavior of Fiber Reinforced Polymer (FRP) composites, Composite Structures 201: 291-302.
[19] Simo J., Taylor R., 1985, Consistent tangent operators for rate-independent elastoplasticity, Computer Methods in Applied Mechanics and Engineering 48: 101-118.
[20] O’Higgins R.M., McCarthy C.T., McCarthy M.A., 2012, Effects of shear-transverse coupling and plasticity in the formulation of an elementary Ply composites damage model, Part II: material characterization, Strain 48(1): 59-67.
[21] O’Higgins R.M., McCarthy C.T., McCarthy M.A., 2012, Effects of shear-transverse coupling and plasticity in the formulation of an elementary Ply composites damage model, Part I: model formulation and validation, Strain 48(1): 49-58.
[22] Truong T.C., Vettori M., Lomov S., Verpoest I., 2005, Carbon composites based on multi-axial multi-ply stitched preforms, Applied Science and Manufacturing 36(9): 1207-1221.