MPM Archiving Options

These commands determine the way results of the calculations are archived for later analysis.

Archiving Particle Data

In MPM calculations, the calculation results are stored in a series of archive files. These files give periodic snapshots for the state of each material point. These files are read after the analysis to plot the results; in fact the results cannot be visualized without them. For more details, you can refer to the output file format definition.

Scripted Input Files

In scripted input files, the archiving options have two forms:

Archive (path)
ArchiveTime (time),<(first time)>
ToArchive (list on quantities)

or

ArchiveUnique (path)
ArchiveTime (time),<(first time)>
ToArchive (list on quantities)

where

(path) gives a relative path name from the saved output file to the archived files. The path name should be entered in Unix style or
```
folder/folder/folder/root
```
where there can be any number of folders followed by the root archive file name. The saved files will be in that location and have names
```
root.num
```
where num is the step number being archived. The folder names and root name cannot contain any colons (":"), to avoid conflict with MacOS path names) or spaces (to be consistent with most Unix systems).
If a series MPM calculations are saved to the same folder without changing the archive path name, the later calculations will most likely overwrite the earlier calculations causing you to lose results. The solution is to use the second form above that replaces the Archive command with an ArchiveUnique command. The ArchiveUnique command forces creation of a new folder within the last folder of the path. In the above path example, the ArchiveUnique option will write to the files:
```
folder/folder/folder/#/root.num
```
where # is chosen from 1, 2, 3, ... to create a unique folder. You can then safely write multiple MPM output files to the save folder without conflicts
(time) gives the time interval between archiving results in milliseconds (it can be an entity). The analysis will always archive the initial state (to define the problem). You only need to archive enough results to get good plots or movies after the analysis is done.
(first time) is an optional time to start archiving in milliseconds (it can be an entity). If this time is greater than zero, the analysis will archive the initial state (to define the problem) and not archive anymore more results until this time is reached. After this time, archiving will proceed at the time interval specified in (time).
(list on quantities) is a comma-separted list of quantities to be archived (see ToArchive command for the details.

ToArchive Command

The ToArchive command determines which calculation results to include in archived results by providing a comma-separated list of the following options (which can be in one or more ToArchive commands and are case sensitive):

velocity - particle velocities
stress - particle stress
strain - particle strain, which is elastic strain for plasticity materials
plasticstrain - particle plastic strain
strainenergy - particle strain energy, which is cumulative σ.d(ε-ε_res)
plasticenergy - cumulative dissipated energy
workenergy - particle work energy, which is cumulative σ.dε
temperature - particle temperature
concentration - particle concentration and concentration gradients
shearcomponents - separate components of total shear strain (namely u/y and v/x; these sum to total shear strain)
jintegral - crack tip J1 and J2 integrals
stressintensity - crack tip KI and KII
energybalance - energy balance results during crack growth by energy balance.
history# - where '#' is 1 to 4 to archive that history variable for a material. This history variable is only relevant for certain types of materials.
thermalenergy - approximate particle thermal energy (ρ C_p ΔT²/(2 T₀) where T₀ is the stress free temperature) - only meaningful when doing conduction calculations and only accurate when the temperature changes are small.
elementcrossings - number of times the particle has crossed an element boundary since the last archive time.
rotstrain - rotational strain. Archiving this strain is needed to plot transformed material points.

Certain particle results are always archived and thus there is no need to specify them in a ToArchive command. The standard archived data are:

mass - particle mass
material ID - material ID number for particle
material angle - material axis angle for anisotropic materials
thickness - thickness for 2D calculations
position - current particle position
original position - initial position for displacement calculations
crack position - crack particle positions
crack original position - initial crack particle positions
crack surface positions - location of crack surfaces

Finally, note that if no ToArchive command is used in an input file, the default archived quantities will be equivalent to the commands:

ToArchive velocity,stress,strain,plasticstrain
ToArchive plasticenergy,workenergy,temperature

But if one or more ToArchive commands are used, you have to specify all quantities you want (i.e., the default ones are not automatically included).

XML Input Files

In XML input files, the archiving commands, which must be within the <MPMHeader> element, are:

<ArchiveRoot unique='0'>(path)</ArchiveRoot>
<ArchiveTime>(time)</ArchiveTime>
<FirstArchiveTime>(first time)</FirstArchiveTime>
<MPMArchiveOrder>iYYYYNYYYNYYNNNNNY</MPMArchiveOrder>
<CrackArchiveOrder>iYYYN</CrackArchiveOrder>

where

(path) is the archive path as explain above.
(time) is the time interval between archiving. A units attribute determines the time units.
(time) is the first time results are archived. A units attribute determines the time units.

The calculations results included in the archive files are determined by a series of 'Y' and 'N' (or other) flags in the strings in the MPMArchiveOrder and CrackArchiveOrder commands.

<MPMArchiveOrder>

Byte 1: Byte order of archived files ('m' for old Macintosh order or 'i' for Intel chip order). This parameter is ignored on input but must be supplied for alignment.
Byte 2: Default particle properties (particle mass, position, material information, and thickness). This byte must always by 'Y'.
Byte 3: Particle velocity
Byte 4: Particle stress
Byte 5: Particle strain (which is elastic strain for plasticity materials)
Byte 6: Particle plastic strain
Byte 7: Should always be 'N'. Present to support reading old versions of archive files.
Byte 8: Particle work energy, which is cumulative σ.dε
Byte 9: Particle temperature
Byte 10: Particle plastic energy
Byte 11: Should always be 'N'. Present to support reading old versions of archive files.
Byte 12: Particle total shear strain components (du/dy and dv/dx)
Byte 13: Particle strain energy, which is cumulative σ.d(ε-ε_res)

Byte 14: Setting 'Y' or 'N' refers to particle history variable 1, but up to four history variables can be archived by using the following characters in this byte:

'1' or 'Y' - history 1
'2' - history 2
'3' - history 1 and 2
'4' - history 3
'5' - history 1 and 3
'6' - history 2 and 3
'7' - history 1, 2, and 3
'8' - history 4
'9' - history 1 and 4
':' - history 2 and 4
';' - history 1, 2, and 4
'<' - history 3 and 4
'=' - history 1, 3, and 4
'>' - history 2, 3, and 4
'?' - history 1, 2, 3, and 4

History variables are only meaningful for some material types.

Byte 15: Particle concentration and concentration gradients - only meaningful when doing diffusion calculations.
Byte 16: Approximate particle thermal energy (ρ C_p ΔT²/(2 T₀) where T₀ is the stress free temperature) - only meaningful when doing conduction calculations and only accurate when the temperature changes are small.
Byte 17: Number of times this particle crossed an element boundary since the last archive time.
Byte 18: Rotational strain; this strain must be archived for visualization tools to be able to plot "transformed" points (i.e., the stretch and shear plotted points to match the material point deformation state).

Default value: iYNN... - i.e. one Y and the rest N

<CrackArchiveOrder>

Byte 1: Byte order of archived files ('m' for Macintosh order or 'i' for Intel chip order). This parameter is ignored on input but must be supplied for alignment.
Byte 2: Default crack particle properties (location, surface positions). This byte must always by 'Y'.
Byte 3: J Integral at designated crack tips.
Byte 4: Stress intensity factors at designated crack tips.
Byte 5: Dynamic energy balancing results (in research, not meant for general use).

Default value: iYNNNN.

Archiving Mesh Data in VTK Files

See Custom Tasks.

Notes

There are no options for archiving kinetic energy because it is easily calculated in visualization software from mass and velocity. Thus, to do kinetic energy calculations, you must archive particle velocity.
Total energy is sum of strain energy and kinetic energy. To visualize total energy results, you must archive strain energy and velocity (for kinetic energy)
Potential energy is total energy minus external work. To visualize internal energy results, you must archive strain energy, velocity (for kinetic energy), and internal work. If the problem has not external work (i.e., force loaded particles), potential energy is equal to total energy.