**j_integral_lorenzi#

Description#

This option calculates the \(J\) integral using the method of Horst and DeLorenzi [U7]. The implementation is analogous to the **parks (see parks ) in its use of a virtual crack extension. This method is more reliable (accurate, insensitive to mesh density, and path independent) for small deformation analysis than the output method *J described in *J . It is also valid for axisymmetric or 3D analysis where *J is not. However, unlike that method the **j_integral_lorenzi is limited to small deformation, and can not calculate \(\Delta J\).

Remember that is the case of linear elasticity:

(100)#\[\begin{split}J = \begin{matrix} \frac{K^2}{ E}(1-\nu^2) & \hbox{ plane strain} \\ \frac{K^2}{ E} & \hbox{ plane stress} \\ \end{matrix}\end{split}\]

Syntax#

**j_integral_lorenzi \(~\,~\,\) perturb (elset | tip) (name1  node_num) \(~\,~\,~\,~\,\) % continuing with as many lines as \(G\) calculations \(~\,\) [ tip_radius val ] \(~\,~\,\) da (dax day)

Refer to the **parks option for a complete description of these options.