Constitutive and numerical framework for modeling joints and faults in rock

A three-dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre-existing gaps or the formation of gaps. Representative stress-strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright (c) 2015 John Wiley & Sons, Ltd.