Numerical study of composite laminates subjected to low-velocity impact using a localized damage algorithm of Puck's 3D IFF criterion. (1st April 2020)

• Record Type: Journal Article
• Authors: Gu, Fuwei; Yuan, Xinyi; Zhu, Xiaolei; Lu, Xiaofeng; Fang, Daining; Li, Li
• Is Part Of: Engineering fracture mechanics. Volume 228(2020)
• Publication Date: 2020-04-01
• Subjects: Localized damage -- Low-velocity impact -- Failure criteria -- Computational cost -- User subroutine; Fracture mechanics -- Periodicals; Rupture, Mécanique de la -- Périodiques; Fracture mechanics; Periodicals; 620.112605
• DOI: 10.1016/j.engfracmech.2020.106901
• Languages: English
• ISSNs: 0013-7944

Abstract

Highlights: The localized damage characteristic of FRP under impact was summarized. A localized damage algorithm was proposed. A simplified calculation method of delamination area was established. Three-dimensional damage envelope was obtained by developing user subroutine. Abstract: The damage area of the composite laminates subjected to low-velocity impact has a localized damage characteristic. Based on this characteristic, a localized damage algorithm applied in a finite element model of composite materials under low-velocity impact was proposed. Firstly, the laminated plate subjected to low-velocity impact was simplified as a simply supported isotropic thin rectangular plate subjected to a central concentrated force. And the localized damage area defined as a reference value was conservatively estimated by the Navier deflection solution. Then, the numerical analysis was performed on a composite laminate specimen under the impact energy of 4J and 8J. By developing user subroutine of ABAQUS VUMAT, the progressive damage models based on Puck's 3D IFF criterion was defined in this localized zone and the linear elastic constitutive model was defined in the other zones. Compared with Traversal Search Method and Semi-Analytical Approximation Method applied in Puck's 3D IFF criterion, the localized damage algorithm was in good agreement with experimental results, and had little effect on the calculation accuracy. According to the 3D damage envelope inside the localized impact zone, the matrix tension damage developed perpendicularly to the ply orientation and presented as a peanut-like shape. Finally, compared with Traversal Search Method and Semi-analytical Approximation Method, the localized damage algorithm could greatly improve the computational efficiency.