MPM Global Archiving Options

From OSUPDOCS
Revision as of 21:58, 9 September 2013 by Nairnj (talk | contribs)
Jump to navigation Jump to search

Global archiving is an option to save selected calculated results in a single file. This section explains how to use this feature.

Introduction

In MPM calculations, the solution is stored in a series of archive files. You can optionally create a global results file which will store selected calculated result such as average stress, total energy, etc., in a single text file. This file can plotted latter in NairnFEAMPM or in NairnFEAMPMViz or can be opened later any spread-sheet software. The file is tab delimited with time (in msec) in the first column and all requested results in subsequent columns.

Input Commands

In script input files, global archiving is set up with the following commands:

GlobalArchive (quantity),<(material)>>
GlobalArchiveTIme (time)

In XML input files, global archive is set up with the following commands, which must be within the <MPMHeader> element:

<GlobalArchive type='(quantity)' mat='(material)'/>
<GlobalArchiveTime>1</GlobalArchiveTime>

where

  • (quantity) is the global quantity to be save in the global results file. The case-sensitive options are:
    • sxx, syy, szz, sxy, sxz, or syz - average element of the stress tensor in MPa
    • sRR, sZZ, sTT, or sRZ - average element of the axisymmetric stress tensor in MPa (synonyms for sxx, syy, szz, and sxy)
    • exx, eyy, ezz, exy, exz, or eyz - average element of the total strain tensor in %, which is sum of elastic and plastic strain for plasticity materials
    • eRR, eZZ, eTT, or eRZ - average element of the axisymmetric total strain tensor % (synonyms for exx, eyy, ezz, and exy)
    • exxe, eyye, ezze, exye, exze, or eyze - average element of the strain tensor in %, which is elastic strain only for plasticity materials
    • eRRe, eZZe, eTTe, or eRZe - average element of the axisymmetric elastic strain tensor % (synonyms for exxe, eyye, ezze, and exye)
    • exxp, eyyp, ezzp, exyp, exzp, or eyzp - average element of the plastic strain tensor in %
    • eRRp, eZZp, eTTp, or eRZp - average element of the axisymmetric plastic strain tensor % (synonyms for exxe, eyye, ezze, and exye)
    • velx, vely, or velz - average component of velocity in mm/sec
    • velR or velZ - average component of axisymmetric velocity in mm/sec (synonyms for velx or vely)
    • dispx, dispy, or dispz - average component of displacement in mm
    • dispR or dispZ - average component of axisymmetric displacement in mm (synonyms for dispx or dispy)
    • contactx, contacty, or contactz - component of the total contact force on the grid for <a href="multimaterialmode.html">multimaterial mode</a> simulations when they include <a href="material.html#rigid">rigid materials</a> that have direction=8. It is a sum of all contact forces for the rigid material on the object in Newtons. If the simulation has activated the <a href="customtask.html#VTKArchive">VTKArchive custom task</a>, that task must archive contactforce or else these global archive results might get invalid zero forces.
    • contactR or contactZ - component of the total contact force per radian on the axisymmetric grid in N (synonyms for contactx or contacty)
    • reactionx, reactiony, or reactionz - component of the reaction force at nodes with velocity boundary conditions in N. If #2 specifies a material, the force will be for all velocity conditions created by that rigid material; if #2 is omitted the force will sum all velocity boundary conditions (specified or created by rigid particles); if #2 &lt 0, the force will be only for <a href="velocity.html">velocity</a> boundary conditions with that <a href="boundaryid.html">boundary condition ID</a>.
    • reactionR or reactionZ - component of the reaction force per radian at axisymmetric nodes with velocity boundary conditions in N (synonyms for reactionx or reactiony and see them for meaning of #2)
    • temp - average temperature
    • concentration - weight fraction concentration (when doing <a href="diffusion.html">diffusion</a> calculations)
    • Strain Energy - total strain energy in J
    • Heat Energy - total heat energy in J
    • Kinetic Energy - total kinetic energy in J on the particles
    • Grid Kinetic Energy - total kinetic energy in J on the grid
    • Energy - sum of total strain and kinetic energies in J
    • Interface Energy - total energy associated with cracks having imperfect interfaces in J. Sum "Strain Energy" and "Interface Energy" to get total strain energy in an elastic material with imperfect interfaces.
    • External Work - total external work in J
    • Potential Energy - total energy minus total external work in J
    • Plastic Energy - total plastic energy in J
    • Thermal Energy - approximate total thermal energy (ρ Cp ΔT2/(2 T0) where T0 is the <a href="stressfreetemp.html">stress free temperature</a>) - only meaningful when doing <a href="conduction.html">conduction</a> calculations and only accurate when temperature changes are small.
    • alpha - the evolving global damping coefficient, α, when using <a href="feedbackdamping.html">feedback damping</a> (in 1/sec)
    • Step number - the current MPM step number
    • Elapsed time - elapsed clock time for the current calculation (in secs)
    • CPU time - total CPU time for the current calculation (in secs)
  • #2 is the material ID for a previously defined <a href="material.html">material</a>. Most of the global quantities are averaged over particles. You can specify this optional second parameter with a material ID and the average will be only for particles of that material type. Omit argument #2 to indicate an average over all particles. Most quantites on the grid (e.g., "Grid Kinetic Energy" but not "reactionx(y,z,R,Z)") will ignore the material ID.
  • #3 is the time interval for storing global results in milliseconds (it can be an <a href="entity.html">entity</a>). If the time interval is never specified, it will be set to archiving interval specified in the <a href="archive.html">Archive</a> command. Specifying another time lets these two archiving times be different if desired.
  • Notes

    1. You can store as many quantities as you want in the global results file by having any number of GlobalArchive commands, but only one option can be in each command.