Table of Contents
List of computational modules and programs
PMD is divided into several computational modules that provide the necessary programs for the computation of a given class of problems.
The programs marked with an asterisk are available in two variants: the standard variant, which employs the frontal solver only, and the alternative variant, which also employs the parallel sparse direct solver. Alternative programs have the letter A appended to their name (e.g., FEFS → FEFSA, HPLS → HPLSA, etc.).
Within the computation the same linear solver must be used; it is selected by the key $\mathtt{KSOL}$ in the input file name.i4 or name.iB.
The digits 2 and 3 in some programs’ names (e.g., RMD2 and RMD3, HPP2 and HPP3, etc.) distinguish programs designated for 2D problems from programs designated for 3D problems. If the incorrect program in terms of problem dimension is run within the computation, older releases of PMD (up to version 12) will abort with an error.
In newer releases of PMD (version 13 and above) the corresponding pairs of programs are still available for backward compatibility, but they are identical (merged) internally and, therefore, it does not matter which one is eventually used within the computation.
The problem dimension is determined automatically at the start of the computation; if the parameter $\mathtt{ITED}\in\{4,5,6,7\}$ or the key $\mathtt{KSS}\ne0$ is specified in the input file name.i1, a 2D problem is assumed, otherwise a 3D problem is assumed.
In this way the majority of problems are identified correctly without any need to modify the input file.
The solution of heat transfer problems is divided between four programs (XT2S, XT2T, XT3S, and XT3T), each of which is designated for one combination of problem dimension and type (2D / 3D, steady-state / transient). If the incorrect program in terms of problem dimension/type is run within the computation, older releases of PMD (up to version 12) will abort with an error.
In newer releases of PMD (version 13 and above) all four programs are still available for backward compatibility, but they are identical (merged) internally and, therefore, it does not matter which one is eventually used within the computation.
The problem dimension is determined by the parameters in the input file name.i1 (see above) and the problem’s steady-state nature by the key $\mathtt{KSU}$ in the input file name.iB.
Heat transfer module
| Program | Function | Input | Log |
|---|---|---|---|
| XRM2* | Process the description of geometric properties, 2D problem | name.i1 | name.o1 |
| XRM3* | Process the description of geometric properties, 3D problem | name.i1 | name.o1 |
| XRPD | Process the description of physical properties and solver control parameters | name.iB | name.oB |
| XT2S* | Solve a steady-state 2D heat transfer problem (modified Newton-Raphson method) | name.oT |
|
| XT3S* | Solve a steady-state 3D heat transfer problem (modified Newton-Raphson method) | name.oT |
|
| XT2T* | Solve a transient 2D heat transfer problem (modified Newton-Raphson method) | name.oT |
|
| XT3T* | Solve a transient 3D heat transfer problem (modified Newton-Raphson method) | name.oT |
Elastostatics module
| Program | Function | Input | Log |
|---|---|---|---|
| RMD2* | Process the description of geometric properties, 2D problem | name.i1 | name.o1 |
| RMD3* | Process the description of geometric properties, 3D problem | name.i1 | name.o1 |
| RPD2 | Process the description of physical properties, 2D problem | name.i2 | name.o2 |
| RPD3 | Process the description of physical properties, 3D problem | name.i2 | name.o2 |
| SRH2 | Calculate element stiffness matrices and right-hand side vectors, 2D problem | name.i3 | name.o3 |
| SRH3 | Calculate element stiffness matrices and right-hand side vectors, 3D problem | name.i3 | name.o3 |
| FEFS* | Solve a linear elastostatic problem (frontal method) | name.i4 | name.o4 |
| STR2 | Calculate strains and stresses, 2D problem | name.i5 | name.o5 |
| STR3 | Calculate strains and stresses, 3D problem | name.i5 | name.o5 |
Dynamics module
This module also requires the programs of the Elastostatics module.
| Program | Function | Input | Log |
|---|---|---|---|
| HMOT | Calculate element mass matrices | name.iM | name.oM |
| HCRE | Calculate element damping matrices | name.iC | name.oC |
| HFRO | Factorize the effective global stiffness matrix for dynamic solvers (Frontal method) | name.iR | name.oR |
| HEIG* | Solve a generalized eigenvalue problem (subspace iteration method) | name.iE | name.oE |
| HFRQ | Print normalized eigenmodes and eigenfrequencies | name.iF | name.oF |
| HMOD* | Solve a dynamic problem (modal superposition method) | name.iD | name.oD |
| HNEW | Solve a dynamic problem (direct integration method) | name.iW | name.oW |
Nonlinear module
This module also requires the programs of the Elastostatics module.
| Program | Function | Input | Log |
|---|---|---|---|
| HPP2 | Prepare data for nonlinear solvers, 2D problem | name.iP | name.oP |
| HPP3 | Prepare data for nonlinear solvers, 3D problem | name.iP | name.oP |
| HPLS* | Solve a nonlinear static problem (BFGS method) | name.iL | name.oL |
| HDYN | Solve a nonlinear dynamic problem (central difference method) | name.iN | name.oN |
Stability module
This module also requires the programs of the Elastostatics module and Dynamics module.