Image-based goal-oriented adaptive isogeometric analysis with application to the micro-mechanical modeling of trabecular bone

• Published in: Computer methods in applied mechanics and engineering 2015-02, Vol.284, p.138-164
• Main Authors: Verhoosel, C.V., van Zwieten, G.J., van Rietbergen, B., de Borst, R.
• Format: Article
• Language: English
• Subjects: Bones; Computation; Computer simulation; Finite cell method; Goal-oriented error estimation and adaptivity; Hierarchical splines; Image-based analysis; Isogeometric analysis; Mathematical analysis; Mathematical models; Methodology; Modulus of elasticity; Splines; Three dimensional; Trabecular bone micro-structures
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2014.07.009

Isogeometric analysis (IGA) of geometrically complex three-dimensional objects is possible when used in combination with the Finite Cell method (FCM). In this contribution we propose a computational methodology to automatically analyze the effective elastic properties of scan-based volumetric objects of arbitrary geometric and topological complexity. The first step is the reconstruction of a smooth geometry from scan-based voxel data using a B-spline level set function. The second step is a goal-oriented adaptive isogeometric linear elastic analysis. Elements are selected for refinement using dual-weighted residual shape function indicators, and hierarchical splines are employed to construct locally refined spline spaces. The proposed methodology is studied in detail for various numerical test cases, including the computation of the effective Young’s modulus of a trabecular bone micro-structure reproduced from μCT-scan data. •A smooth geometry is reconstructed from voxel data using a B-spline approximation.•A goal-adaptive isogeometric finite cell method is used to compute elastic properties.•A tessellation-based procedure is proposed to accurately integrate trimmed elements.•The proposed methodology is applied to μCT-scan data of a trabecular bone specimen.