A peridynamics material model based on Gaussian function for composite laminate

Purpose The numerical modeling of damage in composite laminates remains a challenge due to their heterogeneous structure and complex failure mechanisms. This study proposes a new material model for fiber-reinforced composite laminates and transversely isotropic materials within the framework of bond-based peridynamics (PD). Design/methodology/approach Gaussian functions are introduced to characterize the micromodulus and critical bond stretch, capturing anisotropic behavior. A FORTRAN-based implementation was developed on the Visual Studio platform to simulate damage evolution and crack propagation. Model accuracy and convergence were assessed through tensile simulations by comparing the off-axis modulus with analytical solutions. Horizon size and grid density were optimized to ensure computational reliability. Findings Validation was performed using unidirectional multilayer laminates with a central orifice, showing strong agreement with experimental damage modes and crack paths. Additionally, the predicted stiffness, strength and failure patterns of multidirectional laminates closely matched experimental results, demonstrating the model’s effectiveness for complex composite structures. Originality/value This study proposes a new material model for fiber-reinforced composite laminates and transversely isotropic materials within the framework of bond-based peridynamics (PD). Gaussian functions are introduced to characterize the micromodulus and critical bond stretch, capturing anisotropic behavior.