Complete constitutive model of CFRP including continuous damage in low strain rate compression and temperature generation in high strain rate impact

Abstract This paper reports a method to construct the complete constitutive model of carbon fiber‐reinforced plastic (CFRP) by distinguishing the strain rate state and modifying the material stiffness. In quasistatic (low strain rate) compression, the initial nonlinear effect caused by the material defects was described by introducing an nonlinear increasing factor. In dynamic (high strain rate) impact, the nonlinear strengthening effect was described by introducing the dynamic amplification factors of stress and strain based on the reference state. Considering the instantaneous temperature rise obtained from the impact work–temperature generation equation and the influence of temperature on modulus and strength, the coupled thermomechanical constitutive models combining the strain rate effect and the temperature effect were established. The user‐defined material subroutines (VUMAT) were developed, and then these constitutive models were verified by finite element analysis (FEA). Finally, the developed material subroutine was applied to predict the impact penetration of CFRP laminates, and the results show that the opening holes, damage and energy dissipation are in good agreement with the reference experiment. This work comprehensively analyzed the construction method of CFRP constitutive models, which would provide a guidance for the coupled thermomechanical behavior under dynamic impact. Highlights An anisotropic continuum mechanical constitutive behavior considering strain rate effect and damage evolution behavior was established. A macroscopic anisotropic constitutive model including the coupled impact temperature generation‐mechanical behavior was established. Combined with the incremental finite element method, the three‐dimensional user‐defined material subroutines (VUMAT) were developed and verified. The thermomechanical impact penetration behaviors of CFRP laminates was predicted, and the stress field, temperature field and failure characteristics were visualized.