The preferential orientation for the initiation of a crack is the one which gives the maximum strain energy reduction for a given crack length. This proposed criterion is a logical extension of the maximum energy release rate criterion. It makes no assumptions on the configuration, the homogeneity, the stress condition on the crack faces, or the material response, consequently it is applicable to the usual engineering cases as well as to cases under compression and/or high confining pressures such as obtain inside the Earth. Numerical results for brittle materials (rocks) agree with laboratory and field data, and show that the criterion is an improvement over the empirical and approximate Coulomb-Mohr criterion which has been used for compressive fracture problems for more than 200 years. They also show that our method can be used in cases where it is not a priori evident whether the fracture will remain closed or will open.The mathematical formulation of the criterion is approached by way of constrained optimization, and the solution is proven to exist uniquely. The numerical implementation is based on a finite element scheme. An iterative method is employed to handle the material and geometric non-linearities.
NomenclatureA = stiffness matrix S 0, = minimum strain energy B = active set s = strain energy density B = boundary of the body 7" = prescribed traction on Bt Bc = contact fracture faces T = superscript for transpose Bd = detached fracture faces 0 = prescribed displacement on B. B: = fracture faces ff = displacement vector Bt = traction boundary V = domain of a body B. = displacement boundary ~? = position vector D = constitutive matrix x, y = spatial coordinates E = Young's modulus ~, 2 = indices /~ = body force bS = normalized strain energy reduction H = displacement solution space e = strain tensor i, j, k = indices 0 = crack orientation g = kinematically admissible set Os = preferential crack orientation m = total number of nodes # = Lagrange multiplier N~ = base function p~ = internal friction coefficient n = number of nodes in active set /~s = sliding friction coefficient P = potential energy v = Poisson's ratio p = number of nodes for prescribed displacement a = stress tensor q = m -p a. = normal stress R = real line cr t = tangential stress S = strain energy