Post-processing and visualization techniques for isogeometric analysis results

Abstract Isogeometric Analysis (IGA) introduced in 2005 by Hughes et al. (2005)  [1] exploits one mathematical basis representation for computer aided design (CAD), geometry and analysis during the entire engineering process. In this paper we extend this concept also for visualization. The presented post-processing and visualization techniques thereby strengthen the relation between geometry, analysis and visualization. This is achieved by facilitating the same mathematical function space used for geometry and analysis also for post-processing and visualization purposes. During non-linear analysis derivatives are incrementally computed and stored with different basis function representations. We introduce and investigate projection methods to be able to use the same function space for both displacements and stresses without loss of accuracy. To obtain a common representation for structured and unstructured meshes like hierarchical spline, locally refined B-spline (LR B-spline) and T-spline techniques we exploit Bezier decomposition in a post-processing step resulting in a Bezier element representation and constitute it as generalized representation. The typically used unrelated (fictitious) finite element mesh representation for visualization purposes are easily replaced without changing the underlying geometry as well as the algorithmic data structure. One further benefit of the used Bezier decomposition lies in the fact that it facilitates a natural parallel implementation on Graphics Processor Units (GPUs) exploiting shader programming. In this paper we have developed and investigated an accurate, efficient and practical post-processing pipeline for visualization of isogeometric analysis results . The proposed IGA visualization pipeline consists of three steps: (1) Projection , (2) Bezier decomposition and (3) Pixel-accurate rendering . We have tested four different projection methods. A description on how to perform Bezier decomposition of LR B-splines are given (whereas for hierarchical and T-splines this has been done before). Furthermore, the use of GPU shader programming to enable efficient and pixel-accurate visualization is detailed. The performance of the four different projection techniques has been tested on manufactured problems as well as on realistic benchmark problems. Furthermore, the IGA visualization pipeline has been demonstrated on a number of real-world applications.

[1]  Ingrid Scholl,et al.  Comparing GPU-based multi-volume ray casting techniques , 2010, Computer Science - Research and Development.

[2]  Tore Børvik,et al.  Impact against X65 steel pipes - An experimental investigation , 2013 .

[3]  Tony DeRose,et al.  Feature-adaptive GPU rendering of Catmull-Clark subdivision surfaces , 2012, TOGS.

[4]  Trond Kvamsdal,et al.  Adaptive isogeometric finite element analysis of steady‐state groundwater flow , 2016 .

[5]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[6]  Leila De Floriani,et al.  Multiresolution modeling and visualization of volume data based on simplicial complexes , 1994, VVS '94.

[7]  M. Levoy Volume Rendering Display of Surfaces from Volume Data , 1988 .

[8]  H.-F. Pabst,et al.  Ray Casting of Trimmed NURBS Surfaces on the GPU , 2006, 2006 IEEE Symposium on Interactive Ray Tracing.

[9]  Jim X. Chen,et al.  OpenGL Shading Language , 2009 .

[10]  Adrian J. Lew,et al.  Variational time integrators for finite‐dimensional thermo‐elasto‐dynamics without heat conduction , 2011 .

[11]  Trond Kvamsdal,et al.  On the similarities and differences between Classical Hierarchical, Truncated Hierarchical and LR B-splines , 2015 .

[12]  John A. Evans,et al.  An Isogeometric design-through-analysis methodology based on adaptive hierarchical refinement of NURBS, immersed boundary methods, and T-spline CAD surfaces , 2012 .

[13]  Elaine Cohen,et al.  Direct Isosurface Visualization of Hex-Based High-Order Geometry and Attribute Representations , 2012, IEEE Transactions on Visualization and Computer Graphics.

[14]  P. Gould Introduction to Linear Elasticity , 1983 .

[15]  Trond Kvamsdal,et al.  Isogeometric analysis using LR B-splines , 2014 .

[16]  T. Hughes,et al.  Converting an unstructured quadrilateral mesh to a standard T-spline surface , 2011 .

[17]  T. R. Hughes,et al.  Mathematical foundations of elasticity , 1982 .

[18]  Trond Kvamsdal,et al.  Simple a posteriori error estimators in adaptive isogeometric analysis , 2015, Comput. Math. Appl..

[19]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

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

[21]  Kjell Magne Mathisen,et al.  Isogeometric analysis of finite deformation nearly incompressible solids , 2011 .

[22]  Trond Kvamsdal,et al.  Superconvergent patch recovery and a posteriori error estimation technique in adaptive isogeometric analysis , 2017 .

[23]  Heinrich Müller,et al.  Adaptive generation of surfaces in volume data , 2005, The Visual Computer.

[24]  L. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communications.

[25]  A. Stroud,et al.  Gaussian quadrature formulas , 1966 .

[26]  Takashi Kanai,et al.  Surface quality assessment of subdivision surfaces on programmable graphics hardware , 2004, Proceedings Shape Modeling Applications, 2004..

[27]  Markus Geimer,et al.  Interactive Ray Tracing of Trimmed Bicubic Bézier Surfaces without Triangulation , 2005, WSCG.

[28]  I. Wald,et al.  Interactive Isosurface Ray Tracing of Large Octree Volumes , 2006, 2006 IEEE Symposium on Interactive Ray Tracing.

[29]  W. Boehm Inserting New Knots into B-spline Curves , 1980 .

[30]  John A. Evans,et al.  Robustness of isogeometric structural discretizations under severe mesh distortion , 2010 .

[31]  Graham Sellers,et al.  OpenGL SuperBible: Comprehensive Tutorial and Reference , 2007 .

[32]  Jon Mikkelsen Hjelmervik,et al.  Heterogeneous Computing with Focus on Mechanical Engineering , 2009 .

[33]  Marc Levoy,et al.  Display of surfaces from volume data , 1988, IEEE Computer Graphics and Applications.

[34]  Nicholas S. North,et al.  T-spline simplification and local refinement , 2004, SIGGRAPH 2004.

[35]  Odd Sture Hopperstad,et al.  NUMERICAL SIMULATIONS OF SUBMERGED AND PRESSURISED X65 STEEL PIPES , 2013 .

[36]  John A. Evans,et al.  Isogeometric finite element data structures based on Bézier extraction of NURBS , 2011 .

[37]  Robert L. Taylor,et al.  Convergence of an efficient local least-squares fitting method for bases with compact support , 2012 .

[38]  Ahmad H. Nasri,et al.  T-splines and T-NURCCs , 2003, ACM Trans. Graph..

[39]  Tore Børvik,et al.  Damage and Failure in an X65 Steel Pipeline Caused by Trawl Gear Impact , 2013 .

[40]  Tom Lyche,et al.  Polynomial splines over locally refined box-partitions , 2013, Comput. Aided Geom. Des..

[41]  T. N. E. Greville,et al.  Theory and applications of spline functions , 1969 .

[42]  B. Simeon,et al.  A hierarchical approach to adaptive local refinement in isogeometric analysis , 2011 .

[43]  Trond Kvamsdal,et al.  Divergence-conforming discretization for Stokes problem on locally refined meshes using LR B-splines , 2015 .