A noniterative method for robustly computing the intersections between a line and a curve or surface

The need to compute the intersections between a line and a high-order curve or surface arises in a large number of finite element applications. Such intersection problems are easy to formulate but hard to solve robustly. We introduce a non-iterative method for computing intersections by solving a matrix singular value decomposition (SVD) and an eigenvalue problem. That is, all intersection points and their parametric coordinates are determined in one-shot using only standard linear algebra techniques available in most software libraries. As a result, the introduced technique is far more robust than the widely used Newton-Raphson iteration or its variants. The maximum size of the considered matrices depends on the polynomial degree $q$ of the shape functions and is $2q \times 3q$ for curves and $6 q^2 \times 8 q^2$ for surfaces. The method has its origin in algebraic geometry and has here been considerably simplified with a view to widely used high-order finite elements. In addition, the method is derived from a purely linear algebra perspective without resorting to algebraic geometry terminology. A complete implementation is available from this http URL.

[1]  P. Dooren,et al.  An improved algorithm for the computation of Kronecker's canonical form of a singular pencil , 1988 .

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

[3]  Laurent Busé,et al.  Implicit matrix representations of rational Bézier curves and surfaces , 2014, Comput. Aided Des..

[4]  Nicholas M. Patrikalakis,et al.  Shape Interrogation for Computer Aided Design and Manufacturing , 2002, Springer Berlin Heidelberg.

[5]  G. Karniadakis,et al.  Spectral/hp Element Methods for Computational Fluid Dynamics , 2005 .

[6]  Falai Chen,et al.  Implicitization using moving curves and surfaces , 1995, SIGGRAPH.

[7]  Laurent Busé,et al.  Matrix-based implicit representations of rational algebraic curves and applications , 2010, Comput. Aided Geom. Des..

[8]  B. Laurent,et al.  Implicit matrix representations of rational Bézier curves and surfaces , 2014 .

[9]  Peter Wriggers,et al.  Computational Contact Mechanics , 2002 .

[10]  G. Strang.  Introduction to Linear Algebra , 1993 .

[11]  David Moxey,et al.  Curvilinear mesh generation using a variational framework , 2017, Comput. Aided Des..

[12]  Fehmi Cirak,et al.  Subdivision-stabilised immersed b-spline finite elements for moving boundary flows , 2012 .

[13]  P. Wriggers.  Computational contact mechanics , 2012 .

[14]  K. Morgan,et al.  The generation of arbitrary order curved meshes for 3D finite element analysis , 2013 .

[15]  Malcolm A. Sabin,et al.  Interrogation of spline surfaces with application to isogeometric design and analysis of lattice-skin structures , 2018, Computer Methods in Applied Mechanics and Engineering.

[16]  Thomas W. Sederberg,et al.  COMPUTER AIDED GEOMETRIC DESIGN , 2012 .