Rigid Material

Constitutive Law

This MPM material is actually two materials in one, depending on the value for its direction property.

Moving Boundary Condition

If direction is 0 through 7, the material is not really a material and points that are this type of rigid material will not really be part of the MPM analysis. These rigid material points are used to set boundary conditions on the grid that move with the particles. The rigid material points should be initialized to a velocity, temperature, or concentration when they are created in Region commands. As the rigid material points move through the grid, all nodes connected to those points will be assigned grid-based boundary conditions for the selected particle properties of velocity, temperature, and/or concentration. To have the boundary condition vary in time and position, use the settingFunction(2)(3) and valueFunction properties below to set the value rather than using the initial particle properties to pick a constant value.

Rigid Contact Material

If direction is 8, the material will correspond to actual material points that move through the analysis. Since they are rigid, they will move with their prescribed velocity. They will interact with other material points only through multimaterial contact or imperfect interface laws and thus this type of rigid material can only be used in multimaterial mode MPM calculations.

The advantage of this type of rigid material (unlike the other type) is that you can set frictional properties for contact between rigid materials and non-rigid materials or you can use multimaterial methods to model imperfect interfaces. These rigid materials will move in the direction defined by the initial velocity, or, by using settingFunction(2)(3) properties below, they can moved by any prescribed velocity. If two different rigid materials of this type interact at the same node, however, there is a conflict and the calculations will stop. Finally, because this type of rigid material does not set boundary conditions, it cannot be used to set temperature or concentration.

A direction=8 rigid material ignores cracks, but that actually makes it possible to use them for new crack problems. For example, these rigid particles can provide wedge opening of a crack. If the rigid particles run along an explicit crack plane, those particles may contact particles on both sides of the crack, which can move the top and bottom surfaces of the crack apart. It may be only way to wedge open a crack (i.e., the only material that can act like it is within a crack rather than on either side of the crack). All crack calculations (e.g., crack contact, traction laws, imperfect interfaces, updating of crack position, and J integral) will ignore the rigid particles. The interactions with the rigid particles arise only by contact laws between the rigid particles and the various crack velocity fields.

Material Properties

The rigid material properties are set with the following commands:

Property	Description	Units	Default
direction SetDirection	1 through 7: set grid velocities to various components of the rigid particle velocity (see table below). 8: create rigid material that interacts with non-rigid material by contact. 0: sets no velocity components, but can set temperature and/or concentration values. (note: the SetDirection property name can only be use in scripted files)	none	0
temperature	1 to set temperature boundary condition or 0 to not set it. This setting has no affect unless the analysis is doing thermal calculations. Can only be used if direction is 0 through 7. Use direction = 0 to having moving heat source only; use direction = 1 to 7 to combine moving heat source with velocity boundary conditions.	none	0
concentration	1 to set concentration boundary condition or 0 to not set it. This setting has no affect unless the analysis is doing diffusion calculations. Can only be used if direction is 0 through 7. Use direction = 0 to having moving concentration source only; use direction = 1 to 7 to combine moving concentration source with velocity boundary conditions.	none	0
settingFunction settingFunction2 settingFunction3	You can enter one to three user defined functions (settingFunction and settingFunction1 are synonyms) of time, time step, and position. The function should evaluate to the velocity in mm/s and it will set the rigid-particle velocity (any initial velocity setting will be ignored).	none	none
valueFunction	You can enter a user defined function of time, time step, and position. The function should evaluate to the temperature or concentration potential. Any initial temperature or concentration assigned to the rigid particles will be ignored. Because there is only one value function, you cannot independently set temperature and concentration with a function. Thus rigid particles with value functions should only set one of these quantities.	none	none

The direction options 1 through 7 set grid velocities to various components of the rigid particle velocity by bit-wise selection with first three bits being for x, y, and z directions (or R, Z, θ if axisymmetric). A value of 0 sets no direction:

Value	Binary	Set Directions
0	0x000	no direction
1	0x001	x only
2	0x010	y only
3	0x011	x and y
4	0x100	z only
5	0x101	x and z
6	0x110	y and x
7	0x111	x, y, and z

The (up to) three setting functions can set up to three velocity direction, and their meaning depends on the current setting for the direction property:

• For direction = 1 to 7 the functions must be given in order (settingFunction, settingFunction2, and settingFunction3) and they apply to the velocities that are controlled in order (e.g., when controlling y and z with direction=6, the settingFunction is for y velocity and settingFunction2 is for z velocity).
• For direction = 8, the three functions (regardless of order) are for setting x, y, and z velocities and zero to three can be used. Any not controlled will use the initial particle velocity.
• Note that in scripte files, you can use settingFunction1 (or x) as a synonym for settingFunction and you can use settingFunctiony (or z) as a synonym for settingFunction2 (or 3)).

History Data

Contact?

Notes

Note that <a href="multimaterialmode.html">multimaterial mode MPM</a> calculations can combine moving boundary condition rigid particles (direction=0 to 7) with contact rigid particles (direction=8) to provide different functions, but single-material mode MPM can only use boundary condition rigid particles (direction=0 to 7).

Note that creating rigid particles but not setting any velocity direction (direction=0) can be used to create moving temperature or concentration boundary conditions that move over the particles; i.e., the rigid particles will move through the real particles. Such moving temperature or concentration boundary conditions can only move at a constant velocity defined by the initial particle velocity, because the direction=0 setting implies no axis is controlled by a function.

If a rigid particle is added to a problem, but it does not set any velocity, temperature, or concentration (direction=temperature=concentration=0 or omitted), it will have no affect on the calculations. It may plot in visualization software and therefore appear to be part of the simulation, but it will be identical to the results that would be obtained if it had never been added to the problem.