Highlight

PDMATLAB2D: A peridynamics MATLAB two-dimensional code

Crack branching simulation
Crack branching simulation of a pre-notched soda-lime glass thin plate under traction loading. Left: strain energy density. Right: damage.
2023: Q3

Achievement

This paper describes PDMATLAB2D, a newly released meshfree peridynamics code for simulation of two-dimensional fracture problems. To view PDMATLAB2D click here.

Significance and Impact

Peridynamics has become the method of choice for fracture computations for many researchers worldwide, resulting in an increased demand for accessible peridynamics computational tools in recent years. While several peridynamics codes are currently available, few of these codes focus on an entry-level, easily modifiable tool that can suit the wide variety of needs across disciplines for both research and education. PDMATLAB2D provides an entry-level peridynamics computational tool for educational and training purposes as well as an accessible and easily modifiable research tool for researchers who would like to adapt the code for a multitude of peridynamics simulations. 

Research Details

This paper provides an overview of PDMATLAB2D, a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. PDMATLAB2D provides an entry-level peridynamics computational tool, which can serve educational and training goals, and an accessible and easily modifiable computational environment that peridynamics researchers can adapt to simulate different peridynamics problems. The code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking, and it is designed to be extendable for other peridynamic models and computational features. We reviewed the code structure and functions at a high level, while providing confirmations of the correctness of the computations via figures or tables. The MATLAB functions themselves are carefully commented; thus, we focused their descriptions on the rationale and implementation. Numerical examples for wave propagation and crack branching are presented to demonstrate application of PDMATLAB2D. The corresponding input decks should serve users as a reference for running alternative simulations. While the present version of the code is quite flexible and incorporates various peridynamics modeling and computational features, the addition of other components will be considered in the future, such as state-based peridynamic models, displacement and velocity boundary conditions, and contact algorithms.

Overview

PDMATLAB2D is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The purpose of this code is twofold. First, it provides an entry-level peridynamics computational tool for educational and training purposes. Second, it serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulation scenarios. The current version of the code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking. However, the code is designed to be extendable for other peridynamic models and computational features. In this paper, we provide an overview of the code structure and functions with illustrative examples. Due to the integrated computation and postprocessing MATLAB capabilities, PDMATLAB2D can serve as an effective testbed for testing new constitutive models and advanced numerical features for peridynamics computations.

To view the publication click here

 

Last Updated: February 13, 2024 - 8:33 am