PERI2D
Peridynamics for 2D Time-Dependent Problem

PERI2D is a C program which sets up and solves a 2D time-dependent peridynamics problem, by Miroslav Stoyanov;

Usage:

Typical usage:

        ./pdsolver -Nx 10 -Ny 10 -delta 0.2

Options:

--help, print this information;
-s, this is a steady state problem;
-i, this is a time integration problem;
-Nx (integer), specify the number of cells in X direction;
-Ny (integer), specify the number of cells in Y direction;
-delta (double), specify the peridynamics delta (the "horizon");
-dt (double), specify the time discretization parameter dt;
-Tend (double), specify the end time;
-maxi (integer), specify the maximum number of refinement iterations;
-refine, use iterative refinement for the steady state (need a metric);
-mesh (string), specify an input file with the first mesh to use;
-fmesh (string), filename into which the final mesh will be written;
-o (string), filename into which the final solution will be written;

Libraries:

PERI2D requires the SuiteSparse package of sparse solvers, written by Timothy Davis. This package is available at http://www.cise.ufl.edu/research/sparse/SuiteSparse/.

You may need to make some local modifications. Suppose, for instance, that you install SuiteSparse within the directory:

        /panfs/panasas1/users/jburkardt/software/SuiteSparse

Then you need to go to the directory

        /panfs/panasas1/users/jburkardt/software/SuiteSparse/SuiteSparse_config

and modify the file SuiteSparse_config.mk so that the library and include directories are defined as

        INSTALL_LIB = 
        INSTALL_INCLUDE = /panfs/panasas1/users/jburkardt/software/SuiteSparse/Inc

Now, you need to create these two directories, by moving to the SuiteSparse directory and typing:

        mkdir Lib
        mkdir Inc

Now type

        make

and then

        make install

Assuming you don't see any error messages, SuiteSparse has been set up. Now you need to tell PERI2D where this is. In the main PERI2D directory, you should see a file called Make.conf. Edit this file, and set the values of LSS and ISS to correspond to the locations of the SuiteSparse Lib and Inc directories:

        LSS = /panfs/panasas1/users/jburkardt/software/SuiteSparse/Lib
        ISS = /panfs/panasas1/users/jburkardt/software/SuiteSparse/Inc

Now, in the main PERI2D directory, you should type

        make

If this command executes properly, then in the PDSolver directory, you should find the executable program pdsolver, which you can now run.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Languages:

PERI2D is available in a C version.

Related Data and Programs:

PERI1D, a C program which sets up and solves a 1D time-dependent peridynamics problem, by Miroslav Stoyanov;

PERIDYNAMICS_1D_STEADY, a MATLAB library which solves a 1D steady version of the Poisson equation, using the non-local peridynamics model, by Marta D'Elia.

Author:

Miroslav Stoyanov, Oak Ridge National Laboratory.

Reference:

Qiang Du, Max Gunzburger, Rich Lehoucq, Kun Zhou,
Analysis and approximation of nonlocal diffusion problems with volume constraints,
SIAM Review,
Volume 54, Number 4, pages 667-696, 2012.
Max Gunzburger, Rich Lehoucq,
A nonlocal vector calculus with application to nonlocal boundary value problems,
Multiscale Modeling and Simulation,
Volume 8, Number 5, 2010, pages 1581-1598.
Miroslav Stoyanov,
Quadrature for Peridynamics Problems.

Source Code:

peri2d.tar.gz, a GZIP'd TAR file containing the source code.

Examples and Tests:

RUN1: ./peri2d -Nx 10 -Ny 10 -delta 0.1

run1_output.txt, the output file.

RUN2: ./peri2d -Nx 10 -Ny 10 -delta 0.1 -fmesh run2_fmesh.txt

run2_output.txt, the output file.
run2_fmesh.txt, the final mesh.

RUN1: ./peri2d -Nx 10 -Ny 10 -delta 0.1 -o run3_solution.txt

run3_output.txt, the output file.
run3_solution.txt, the solution file.

You can go up one level to the C source codes.

Last revised on 29 November 2012.

PERI2D Peridynamics for 2D Time-Dependent Problem