model domain command

Syntax

model domain keyword

Specify the domain boundaries and conditions. DFN fractures only exist within the domain extent.

Note

In FLAC3D, zones and structural elements are not aware of the model domain.

In PFC, walls clumps and balls must be inside of the domain and are aware of the domain boundary conditions.

condition keyword <keyword keyword>

Set the domain boundary conditions. Each domain condition applies to the two box sides that are perpendicular to the direction.

The first keyword specifies the condition for the x direction. The second for the y, and the third for the z. If only one keyword is specified, then y and z are set to the same condition as x. The default condition is the stop condition.

The following keywords can be used to specify domain conditions.

stop

The stop condition imposes the constraint that if the body centroid falls outside the model domain, the velocity and spin of the body are nulled. For PFC walls, once a facet centroid falls outside the model domain, the velocity and spin are nulled, including any vertex velocities.

reflect

The reflect condition imparts the opposite sign of the velocity, and the spin is unchanged if the body centroid falls outside the model domain. The use of this condition as a contact rule is not recommended since it may induce instabilities. The reflect condition is not compatible with walls intersecting the domain.

destroy

The destroy condition deletes balls, clumps, and wall facets whose individual extents (i.e., an axis-aligned bounding box that tightly encompasses the object) fall outside the model domain.

periodic

The periodic condition applies periodic boundary conditions. When the ball or clump centroid falls outside of the model domain, they are translated back to the opposite side of the model. To ensure that contacts are created as if the model was continuous, “ghost” balls and clumps are introduced. The periodic condition is not compatible with walls intersecting the domain.

extent fxl fxu <fyl fyu fzl fzu >

Specify the locations of the sides of the domain. Each set of values applies to the lower and upper bound of that dimension of the domain. If fxl and fxu are specified alone, then the domain is a cube.

strain-rate fxx fyy fzz fxy fxz fyz fyx fzx fyx

Specify the full strain rate tensor. All components must be specified in order, otherwise subsequent components after the last specified component are zeroed.

strain-rate-xx f

Specify the xx-component of the strain rate tensor. Other components are not modified.

strain-rate-xy f

Specify the xy-component of the strain rate tensor. Other components are not modified.

strain-rate-xz f (3D only)

Specify the xz-component of the strain rate tensor. Other components are not modified.

strain-rate-yx f

Specify the yx-component of the strain rate tensor. Other components are not modified.

strain-rate-yy f

Specify the yy-component of the strain rate tensor. Other components are not modified.

strain-rate-yz f (3D only)

Specify the yz-component of the strain rate tensor. Other components are not modified.

strain-rate-zx f (3D only)

Specify the zx-component of the strain rate tensor. Other components are not modified.

strain-rate-zy f (3D only)

Specify the zy-component of the strain rate tensor. Other components are not modified.

strain-rate-zz f (3D only)

Specify the zz-component of the strain rate tensor. Other components are not modified.

tolerance f

Specify the domain remapping tolerance. By default, this value is set to zero and all pieces are remapped when the domain extent is modified. When a non-zero domain strain-rate is used, it is recommended to use a non-zero remapping tolerance to avoid this and achieve better runtimes.