MPM Region and Hole Commands
In MPM, the most common way to discretize an object is to add material points to the background grid using Region commands and to designate regions of the grid has having no material points using Hole commands. Enclosed in these commands are various shape commands to define the object. An alternative and powerful option for discretizing an object is to create material points automatically from an image of the object.
Region Commands
Once the background grid is created, material points are added to the grid using a series of Region commands. In scripted files, the Region command defines material type, initial velocity, thickness (in 2D), and optionally an initial angle (for some anisotropic materials), temperature, and concentration. The command is
Region (matid),(velx),(vely),(velz) or (thick),<(property),(value)>... (interspersed initial velocity commands) (any number of shape commands) (interspersed rotation commands) EndRegion
In XML files, Region command are replaced by <Body> commands and optionally one <PointList> block. All these elements must be within the single <MaterialPoints> element in the input file:
<MaterialPoints>
<PointList>
(see help on PointList block)
</PointList>
<Body mat='(matid)' vx='(velx)' vy='(vely)' vz='(velz)' thick='(thick)'
angle='(angle)' temp='(temp)' conc='(conc)' wtconc='(wtconc)'>
(interspersed initial velocity commands)
(any number of shape commands)
(interspersed rotation commands)
</Body>
....
</MaterialPoints>
where
- (matid) is the material ID for a previously defined material.
- (velx), (vely), and (velz) are the components of the initial velocity for all particles in the region (velz is only used for 3D simulations). The units are velocity units.
- (thick) is the thickness of each material point in length units. The thickness is only used for planar 2D calculations. In scripted files for axisymmetric calculations, the thickness parameter is still required, but it is only used for parameter alignment (its value is ignored).
Besides the above required parameters, you can provide optional parameters to set initial conditions for all particles in the region. In scripted files, the optional parameters appear in pairs. The first item in the pair (i.e., (property)) is the name of the property to set and the second item (i.e., (value)) is the initial value for that property. In XML files, the optional properties are set using optional attributes:
- angle,(angle) to set the initial rotation angle for the material axes about the z axis. The units are degrees. This parameter is sufficient for 2D analyses with constant material angle. See help on setting material angles for methods to set variable angles or to set more than one rotation angle in 3D calculations.
- temp,(temp) to set the initial particle temperature. The units are degrees.
- conc,(conc) to set the initial particle concentration potential. The dimensionless concentration potential must be from 0 to 1. Only used when doing diffusion calculations.
- (wtconc)to set concentration by weight fraction concentration instead of by concentration potential. The resulting potential will be (wtconc) divided by the material's saturation concentration. This option is only available in XML input files. The units are dimensionless and must be from 0 to the material's saturation concentration. Only used when doing diffusion calculations.
The optional PointList Block can be used to define individual material points, each with individual initial conditions. The option is only available to XML input files.
Initial Velocity Commands
Between the Region command and the subsequent EndRegion command (or within each <Body> element), the shape commands can be interspersed with initial velocity commands to set initial velocity or angular velocity of the particles. The details on these shape command are given elsewhere:
Shape Commands
Between the Region command and the subsequent EndRegion command (or within each <Body> element), there can be any number of shape commands to define material point positions and the shapes. The details on these shape command are given elsewhere:
Rotation Commands
Between the Region command and the subsequent EndRegion command (or within each <Body> element), the shape commands can be interspersed with rotation commands to define anisotropic material orientation. The details on these shape command are given elsewhere:
Hole Commands
The Hole command is used to designate regions of the background grid that should not have material points. In scripted files, the Hole command is:
Hole (any number of shape commands) EndHole
In XML files, the <Hole> block is:
<Hole> (any number of shape commands) </Hole>
Between the Hole command and the subsequent EndHole command (or within the <Hole> block), there can be any number of shape commands to define hole locations. The details on these shape command are given elsewhere:
An object with holes (or voids) is created by a sequence of Region and Hole commands. Each Region command fills every location within its shapes with a material point unless that location was previously designated as a hole. Each Hole command designates every location within its shapes as being a hole, but, if any location already has a material point, that point is left intact. In other words, the order of these commands is crucial. By selecting the correct order, you have greatly flexibility in creating objects with solid regions and void (or hole) regions.