[Bug]: Fix PD Solid Plastic

[JTHesse]

PeriLab Version

1.1.5

Julia Version

1.10.3

MPI Version

No response

What is the expected behavior?

PD Solid Plastic material model is not tested yet!

What is the actual behavior?

Infinite forces

Reproduction steps

julia --project=. src/main.jl test/fullscale_tests/test_PD_solid_plastic/strain_xy.yaml

Relevant log output

[Info] 
 PeriLab version: 1.1.6
 Copyright: Dr.-Ing. Christian Willberg, M.Sc. Jan-Timo Hesse
 Contact: [email protected], [email protected]
 GitHub: https://github.com/PeriHub/PeriLab.jl
 DOI: 10.1016/j.softx.2024.101700
 License: BSD-3-Clause
 ---------------------------------------------------------------

[Info] 2024-05-15 14:08:37
[Warn] Commit: c3180f4912544a763b887f986ed8b5f193f005be, Local changes detected
[Info] Read input file test/fullscale_tests/test_PD_solid_plastic/strain_xy.yaml
[Debug] Read node sets
[Info] Read mesh file test/fullscale_tests/test_PD_solid_plastic/test_mesh.txt
[Info] 2d problem assumed, please define plane stress or plane strain if needed in Physics
[Debug] Init Neighborhoodlist
[Debug] Finished init Neighborhoodlist
[Debug] Start distribution
[Debug] Finished distribution
[Debug] Create Overlap
[Debug] Finished Overlap
  Mesh input overview
─────────────────────────────────────────────────────
                  Number of nodes   ..........   12
  Geometrical degrees of freedoms   ..........    2
[Debug] Get node sets
[Debug] Finish init data
[Info] Init solver
[Debug] Read properties
[Debug] Init models
[Info] Init material models
[Debug] Init pre calculation models
[Debug] Init Boundary Conditions
[Warn] Missing boundary condition type for BC_1. Assuming Dirichlet.
[Warn] Missing boundary condition type for BC_2. Assuming Dirichlet.
[Debug] Init Verlet
  Verlet Solver
────────────────────────────────────────────────────────
            Initial time   ..........          0.0   s
              Final time   ..........       5.0e-6   s
  Minimal time increment   ..........   2.67322e-6   s
           Safety factor   ..........          1.0
          Time increment   ..........       2.5e-6   s
         Number of steps   ..........            2
       Numerical Damping   ..........          0.0
[Debug] Finished Init Solver
┌─────────┬──────────┬────────┬─────────┬───────────┬──────────┬─────────┬─────────┬─────────────────┐
│   Block │ Material │ Damage │ Thermal │ Corrosion │ Additive │ Density │ Horizon │ Number of Nodes │
├─────────┼──────────┼────────┼─────────┼───────────┼──────────┼─────────┼─────────┼─────────────────┤
│ block_1 │    Mat_1 │        │         │           │          │  2.7e-9 │     2.1 │               8 │
│ block_2 │    Mat_2 │        │         │           │          │  2.7e-9 │     2.1 │               4 │
└─────────┴──────────┴────────┴─────────┴───────────┴──────────┴─────────┴─────────┴─────────────────┘
[Debug] Init write results
[Warn] No Output File Type defined for Output1, defaulting to Exodus
[Warn] No Output File Type defined for Output1, defaulting to Exodus
[Info] Create output test/fullscale_tests/test_PD_solid_plastic/strain_xy.e
[Info] Run Verlet Solver
[Error] Forces is infinite.
[Warn] File already closed
[Debug] Close result files

Which release of PeriLab are you seeing the problem on?

• [ ] Julia Package
• [ ] Linux Release
• [ ] Docker Image