124 NURBS-based geometric fracture growth representation

Numerical methods for fracture propagation model fracture growth as a geometric response to deformation. In contrast to the widely used faceted representations, a smooth Non-Uniform Rational B-Spline (NURBS) surface can be used to represent the fracture domain. Its benefits include low cost, resolutionindependent storage, and a parametric representation of a smooth domain. In the present work an interaction-free, deformation-informed, Gaussian-based modification algorithm of the fracture surface is presented, with localized stress intensity factor computations, and automatic resolution adjustment, which allow for geometric evolution without the need of appending or re-approximating the fracture surface. It is based on the movement of surface control points and on the systematic editing of weights and knots. It does not require trimming, and is able to shift fracture shape and capture its path evolution efficiently. Throughout growth, the number of points required for fracture representation remains fixed, and the discretization of the fracture surface is implicitly defined by the underlying parametric space. The proposed algorithm can be incorporated into any fracture propagation code that keeps track of fracture geometry and updates it as a function of deformation. The algorithm is demonstrated for a discrete finite element-based fracture propagation method.

[1]  Bharat K. Soni,et al.  A Concise Representation of Geometry Suitable for Mesh Generation , 2002, IMR.

[2]  John F. Hughes,et al.  Direct manipulation of free-form deformations , 1992, SIGGRAPH.

[3]  Stephan K. Matthäi,et al.  Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks , 2007 .

[4]  George Celniker,et al.  Linear constraints for deformable non-uniform B-spline surfaces , 1992, I3D '92.

[5]  Xiang Zhu,et al.  Modifying the shape of NURBS surfaces with geometric constraints , 2001, Comput. Aided Des..

[6]  M. Daniel,et al.  Controlled Local Deformation of NURBS Surfaces to Satisfy Several Punctual Constraints , 2006 .

[7]  Ole Sigmund,et al.  Isogeometric shape optimization of photonic crystals via Coons patches , 2011 .

[8]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[9]  T. Belytschko,et al.  X‐FEM in isogeometric analysis for linear fracture mechanics , 2011 .

[10]  A. Paluszny,et al.  Numerical modeling of discrete multi-crack growth applied to pattern formation in geological brittle media , 2009 .

[11]  Gerald E. Farin,et al.  Nurbs for Curve & Surface Design: From Projective Geometry to Practical Use , 1999 .

[12]  Guido Dhondt,et al.  Automatic 3-D crack propagation calculations: a pure hexahedral element approach versus a combined element approach , 2009 .

[13]  M. Pourazady,et al.  Direct manipulations of NURBS surfaces subjected to geometric constraints , 2006, Comput. Graph..

[14]  Paul A. Wawrzynek,et al.  Arbitrary crack representation using solid modeling , 1993, Engineering with Computers.

[15]  L. Piegl Modifying the shape of rational B-splines. part2: surfaces , 1989 .

[16]  Robert W. Zimmerman,et al.  Numerical simulation of multiple 3D fracture propagation using arbitrary meshes , 2011 .

[17]  W. Boehm,et al.  Bezier and B-Spline Techniques , 2002 .

[18]  Marc Daniel,et al.  Local Deformation of NURBS Curves , 2005 .

[19]  T. Rabczuk,et al.  A three-dimensional meshfree method for continuous multiple-crack initiation, propagation and junction in statics and dynamics , 2007 .