Isogeometric Boundary Element Methods

figure1

Legato-team members: E. Atroshchenko, H. Lian and X. Peng

During the recent years one of the most promising directions of    research in computational mechanics has been towards the synthesis of the CAD geometries and methods of stress analysis. Since the boundary element method only requires the discretization of the boundary and its accuracy is highly dependent on the accuracy in the boundary representation, the CAD-parameterization of the exact geometry can be directly embedded into the governing equations, yielding highly accurate solutions.

Approach which utilizes the same basis functions for both – description of the geometry and approximation of the solution fields – is known as isogeometric BEM (IGABEM). The Legato team has successfully implemented IGABEM in two applications:

IGABEM for shape optimisation

Shape optimisation is a process to find the optimal shape of a component or structure under given constraints and objectives. In finite element method (FEM), a mesh must be created to approximate the geometry and discretize the governing partial differential equation (PDE) to allow analysis to be performed. Shape optimization is an iterative procedure and geometries vary at each step, thus leading to cumbersome remeshing procedures, which occupy around 80% of the total problem solving time for linear problems at each iterative step and accumulate to an unaffordable computational burden. Moreover, it may happen in industrial practices that the geometry is so complex that available mesh generators fail, or require significant human intervention, which precludes the automated algorithm. So the mesh burden has become the bottleneck to achieving an automated and fast shape optimization. To solve this problem, the isogeometric boundary element method is applied to shape optimization in linear elasticity. Based on the same boundary representation as CAD, IGABEM can achieve a seamless CAD and analysis integration, which renders IGABEM immediate advantages in the application of shape optimization: 1) the mesh generation/regeneration can be avoided, 2) the automation of the optimization process is more realistic because no human intervention is introduced for constructing analysis-suitable model and meshing, 3) the free-form representation for shape optimization can be naturally reached, and 4) the optimal geometrical model output can be used directly in CAD without any postprocessing or “smoothing” procedure.

figure2optimalChairSide

 

IGABEM for fatigue fracture

In conventional fatigue simulations as performed industrially using the finite element based methods, the key difficulty is the accurate computation of the crack driving force, namely the stress intensity factors (SIFs). The second dificulty is that the domain mesh used for stress analysis and hence for the detection of “sensitive” regions in the component, where initial flaws are introduced, is typically at least one order of magnitude too coarse to provide quality SIFs. The third difficulty lies in the geometrical complexity of the domain which, if the predicted fatigue life is deemed inadequate must be redesigned. For each new design, and for each crack configuration, a new mesh typically needs to be generated, not only to conform to the new chosen geometry, but also to properly resolve stresses in the vicinity of the crack tip. Even when enriched finite element methods are used, some level of remeshing is required.

IGABEM is an strong contender to attack fracture mechanics problems, because it requires only boundary discretization, which simplifies the insertion of new crack segments during growth and offers superior accuracy for the computation of the SIFs for the same number of degrees of freedom compared to other methods.

figure3

figure4

 

Future work

Future may focus on the following aspects:

1) 3D fatigue crack growth and cohesive zone modeling

2) IGABEM for contact problem

3) IGABEM–FEM coupling with model reduction

Reference

H. Lian, R. N. Simpson and S. P. A. Bordas, Sensitivity analysis and shape optimization using isogeomgetric boundary element methods. 11th World Congress on Computational Mechanics. 2014. http://orbilu.uni.lu/handle/10993/16044

X. Peng, E. Atroshchenko, R. N. Simpon, S. Kulasegaram and S. P. A. Bordas, Crack growth analysis by a NURBS-based isogeometric boundary element method. 1th World Congress on Computational Mechanics. 2014. http://orbilu.uni.lu/handle/10993/14135

R. N. Simpson, S. P. A. Bordas, J. Trevelyan and T. Rabczuk, A two-dimensional isogeometric boundary element method for elastostatic analysis. Comp. Meth. in App. Mech. and Engng., Vol. 209-212(0), 87-100, 2012. http://orbilu.uni.lu/handle/10993/13849

H. Lian, R. N. Simpson and S. P. A. Bordas, Stress analysis without meshing: Isogeometric boundary-element method. Proceedings of the Institution of Civil Engineers: Engineering and Computational Mechanics, Vol. 166(2), 88-99, 2013. http://orbilu.uni.lu/handle/10993/10039

R. N. Simpson, S. P. A. Bordas, H. Lian and J. Travelyan, An isogeometric boundary element method for elastostatic analysis: 2D implementation aspects. Computers & Structures, Vol. 118, 2-12, 2013. http://orbilu.uni.lu/handle/10993/12157

M. A. Scott, R. N. Simpson, J. A. Evans, S. Lipton, S. P. A. Bordas, T. J. R. Hughes and T.W. Sederberg, Isogeometric boundary element analysis using unstructured T-splines. Comp. Meth. in App. Mech. and Engng., Vol. 254(0), 197-221, 2013. http://orbilu.uni.lu/handle/10993/11850