Constitutive Relation Error Method (CRE) for linear problems using Isogeometric Analysis

Dr. Hoang Phuong Thai

Feb. 4, 2020, 9:15 a.m. S2 416-1

The paper focuses on the application of the Constitutive Relation Error method (CRE) [1] for a posteriori error estimation in Isogeometric Analysis (IGA). IGA enables us to represent exactly the geometry of mechanical problems based on basis functions generated from Non-uniform Rational B-Splines (NURBS) [2]. The exact geometry in IGA can be maintained at the coarsest level of discretization. This helps to directly eliminate the error sources due to the violation of geometry, approximation of applied loads and boundary conditions without additional assumption.

In the present work, the CRE concept is used to assess the discretization error, and potentially other error sources, from the approximate IGA solutions not satisfying the equilibrium conditions. Essentially, the principles for the construction of so-called admissible (equilibrated) fields construction in the framework of IGA are identical to these in FEM. However, since B-splines or NURBS shape functions of IGA are established for control points instead of physical nodes, a relevant nomination for components of force densities was proposed to facilitate their determination. The construction is furthermore performed at the knot span level. Two first numerical examples are presented: a linear elasticity problem of L-shape console beam with a singularity point, and a heat transfer problem in an annulus. They are typical test cases where IGA was proven to be superior to FEM. The obtained results show the potential power of CRE in estimating errors of IGA solutions for linear problems.

[1] P. Ladevèze and J. P. Pelle Mastering Calculations in Linear and Nonlinear Mechanics. Springer, 2005.
[2] J.A. Cottrell, T.J.R. Hughes, and Y. Bazilevs Isogeometric analysis: toward integration of CAD and FEA. John Wiley & Sons,2009