The jointed material model:
is intended to provide a simple continuum model for a material containing a high density of parallel joint surfaces where each system of parallel joints is associated with a particular orientation, such as sedimentary rock;
assumes that the spacing of the joints of a particular orientation is sufficiently close compared to characteristic dimensions in the domain of the model such that the joints can be smeared into a continuum of slip systems;
provides for opening or frictional sliding of the joints in each of these systems (a “system” in this context is a joint orientation in a particular direction at a material calculation point); and
assumes that the elastic behavior of the material is isotropic and linear when all joints at a point are closed (isotropic linear elastic behavior must be included in the material definition; see Defining isotropic elasticity” in “Linear elastic behavior, Section 17.2.1).
The jointed material model is intended primarily for applications where the stresses are mainly compressive. The model provides a joint opening capability when the stress normal to the joint tries to become tensile. In this case the stiffness of the material normal to the joint plane becomes zero instantaneously. ABAQUS/Standard uses a stress-based joint opening criterion, whereas joint closing is monitored based on strain. Joint system a opens when the estimated pressure stress across the joint (normal to the joint surface) is no longer positive:
The shear response of open joints is governed by the shear retention parameter, , which represents the fraction of the elastic shear modulus retained when the joints are open (=0 means no shear stiffness associated with open joints, while =1 corresponds to elastic shear stiffness in open joints; any value between these two extremes can be used). When a joint opens, the shear behavior may be brittle, depending on the shear retention factor used for open joints. In addition, the stiffness of the material normal to the joint plane suddenly goes to zero. For these reasons, in situations where the confining stresses are low or significant regions experience tensile behavior, the joint systems may experience a sequence of alternate opening and closing states from iteration to iteration. Typically such behavior manifests itself as oscillating global residual forces. The convergence rate associated with such discontinuous behavior may be very slow and, thus, prohibit obtaining a solution. This type of failure is more probable in cases where more than one joint system is modeled.
When the repeated opening and closing of joints makes convergence difficult, you can improve convergence by preventing a joint from opening. In this case an elastic stiffness is always associated with the joint. It is most useful when the opening and closing of joints is limited to small regions of the model. You can prevent a joint from opening only when the joint direction is specified, as described below.
|Input File Usage:|
*JOINTED MATERIAL, NO SEPARATION, JOINT DIRECTION
You must specify nonzero shear retention in open joints directly. The parameter can be defined as a tabular function of temperature and predefined field variables.
|Input File Usage:|
*JOINTED MATERIAL, SHEAR RETENTION
The failure surface for sliding on joint system a is defined by
is the rate of inelastic shear strain in direction on the joint surface ( are orthogonal directions on the joint surface),
is the magnitude of the inelastic strain rate,
is a component of the shear stress on the joint surface,
is the dilation angle for this joint system (choosing provides pure shear flow on the joint, while causes dilation of the joint as it slips), and
is the inelastic strain normal to the joint surface.
Up to three joint directions can be included in the material description. The orientations of the joint directions are given by referring to the names of user-defined local orientations (Orientations, Section 2.2.5) that define the joint orientations in the original configuration. Output of stress and strain components is in the global directions unless a local orientation is also used in the material's section definition.
The parameters , , and can be specified as tabular functions of temperature and/or predefined field variables for each joint direction.
In addition to the joint systems, the jointed material model includes a bulk material failure mechanism, which is based on the Drucker-Prager failure criterion:
is the Mises equivalent deviatoric stress,
is the deviatoric stress,
is the equivalent pressure stress,
is the friction angle for the bulk material, and
is the cohesion for the bulk material.
If bulk material failure is to be modeled, a jointed material behavior must be specified to define the parameters associated with bulk material failure behavior. Thus, up to five jointed material behaviors can appear in the same material definition: three joint directions, shear retention in open joints, and bulk material failure.
The parameters , , and can be specified as a tabular function of temperature and/or predefined field variables.
|Input File Usage:|
*JOINTED MATERIAL (the JOINT DIRECTION parameter must be omitted)
If in any joint system, whether it be associated with the joint surfaces or the bulk material, the flow in that system is “nonassociated.” The implication is that the material stiffness matrix is not symmetric. Therefore, the unsymmetric matrix solution scheme should be used for the analysis step (Procedures: overview, Section 6.1.1), especially when large regions of the model are expected to flow plastically and when the difference between and is large. If the difference between and is not large, a symmetric approximation to the matrix can provide an acceptable rate of convergence of the equilibrium equations and, hence, a lower overall solution cost. Therefore, the unsymmetric matrix solution scheme is not invoked automatically when jointed material behavior is defined.
The jointed material model can be used with plane strain, generalized plane strain, axisymmetric, and three-dimensional solid (continuum) elements in ABAQUS/Standard. This model cannot be used with elements for which the assumed stress state is plane stress (plane stress, shell, and membrane elements).