Computing Upper and Lower Bounds for the J-Integral in Two-Dimensional Linear Elasticity
Name
HPCES004.pdf
Size
1.89 MB
Format
Adobe PDF
Checksum (MD5)
107893f23681ba422681713a3435a45f
Author(s) • • •
Xuan, Z.C.
Lee, Kwok Hong
Patera, Anthony T.
Peraire, Jaime
Date Issued
January 2004
Series/Report no.
High Performance Computation for Engineered Systems (HPCES);
Abstract
We present an a-posteriori method for computing rigorous upper and lower bounds of the J-integral in two dimensional linear elasticity. The J-integral, which is typically expressed as a contour integral, is recast as a surface integral which yields a quadratic continuous functional of the displacement. By expanding the quadratic output about an approximate finite element solution, the output is expressed as a known computable quantity plus linear and quadratic functionals of the solution error. The quadratic component is bounded by the energy norm of the error scaled by a continuity constant, which is determined explicitly. The linear component is expressed as an inner product of the errors in the displacement and in a computed adjoint solution, and bounded using standard a-posteriori error estimation techniques. The method is illustrated with two fracture problems in plane strain elasticity.
Subjects
J-integral
fracture mechanics
linear elasticity
Persistent DSpace Link