Techniques to accelerate boundary element contributions in elasticity

The problem of rapid re-analysis of small problems in elasticity is investigated. The aim is to enable updated stress contours to be displayed in real-time as a design geometry is dynamically modified. The focus of this work is small to medium sized problems; as a result it cannot be assumed that the solution phase dominates, and so the evaluation of boundary integrals is considered as well as the equation solution. Two strategies are employed for acceleration of boundary element integrals: the use of Look-Up Tables (LUTs) containing precomputed integrals and the use of approximate analytical expressions derived from surface fits. These may be used in the matrix assembly and internal point calculations. LUTs are derived for both flat and circular arc elements for both the displacement and stress boundary integral equation. Details are provided on suitable LUT refinements and the approach is benchmarked against conventional Gauss-Legendre quadrature. The surface fit approach is presented as an alternative to LUTs that does not incur the considerable memory cost associated with LUTs. This approach has been limited to flat elements. The equation solution is cast in a re-solution framework, in which we use a GM-RES iterative solver. Convergence is greatly accelerated by using an approximate but complete LU preconditioner updated periodically using multi-threading. Consideration of the period of update is investigated with reference to the spread of eigenvalues in the preconditioned system. The resulting system achieves the aim of providing real time update of contours for small to medium size problems on a PC. This development is expected to allow a qualitative change in the way engineers might use computer aided engineering tools, in which design ideas may rapidly be assessed immediately as a change is made

Vademecum-based GFEM (V-GFEM): optimal enrichment for transient problems

This is the accepted version of the following article: [Canales, D., Leygue, A., Chinesta, F., González, D., Cueto, E., Feulvarch, E., Bergheau, J. -M., and Huerta, A. (2016) Vademecum-based GFEM (V-GFEM): optimal enrichment for transient problems. Int. J. Numer. Meth. Engng, 108: 971–989. doi: 10.1002/nme.5240.], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5240/fullThis paper proposes a generalized finite element method based on the use of parametric solutions as enrichment functions. These parametric solutions are precomputed off-line and stored in memory in the form of a computational vademecum so that they can be used on-line with negligible cost. This renders a more efficient computational method than traditional finite element methods at performing simulations of processes. One key issue of the proposed method is the efficient computation of the parametric enrichments. These are computed and efficiently stored in memory by employing proper generalized decompositions. Although the presented method can be broadly applied, it is particularly well suited in manufacturing processes involving localized physics that depend on many parameters, such as welding. After introducing the vademecum-generalized finite element method formulation, we present some numerical examples related to the simulation of thermal models encountered in welding processes.Peer ReviewedPostprint (author's final draft

Topology optimization for human proximal femur considering bi-modulus behavior of cortical bones

© Springer International Publishing Switzerland 2015. The material in the human proximal femur is considered as bi-modulus material and the density distribution is predicted by topology optimization method. To reduce the computational cost, the bi-modulus material is replaced with two isotropic materials in simulation. The selection of local material modulus is determined by the previous local stress state. Compared with density prediction results by traditional isotropic material in proximal femur, the bi-modulus material layouts are different obviously. The results also demonstrate that the bi-modulus material model is better than the isotropic material model in simulation of density prediction in femur bone

A high-performance boundary element method and its applications in engineering

As a semi-numerical and semi-analytical method, owing to the inherent advantage, of boundary-only discretisation, the boundary element method (BEM) has been widely applied to problems with complicated geometries, stress concentration problems, infinite domain problems, and many others. However, domain integrals and non-symmetrical and dense matrix systems are two obstacles for BEM which have hindered the its further development and application. This thesis is aimed at proposing a high-performance BEM to tackle the above two drawbacks and broaden the application scope of BEM. In this thesis, a detailed introduction to the traditional BEM is given and several popular algorithms are introduced or proposed to enhance the performance of BEM. Numerical examples in heat conduction analysis, thermoelastic analysis and thermoelastic fracture problems are performed to assess the efficiency and correction of the algorithms. In addition, necessary theoretical derivations are embraced for establishing novel boundary integral equations (BIEs) for specific engineering problems. The following three parts are the main content of this thesis. (1) The first part (Part II consisting of two chapters) is aimed at heat conduction analysis by BEM. The coefficient matrix of equations formed by BEM in solving problems is fully-populated which occupy large computer memory. To deal with that, the fast multipole method (FMM) is introduced to energize the line integration boundary element method (LIBEM) to performs better in efficiency. In addition, to compute domain integrals with known or unknown integrand functions which are caused by heat sources or heterogeneity, a novel BEM, the adaptive orthogonal interpolation moving least squares (AOIMLS) method enhanced LIBEM, which also inherits the advantage of boundary-only discretisation, is proposed. Unlike LIBEM, which is an accurate and stable method for computing domain integrals, but only works when the mathematical expression of integral function in domain integrals is known, the AOIMLS enhanced LIBEM can compute domain integrals with known or unknown integral functions, which ensures all the nonlinear and nonhomogeneous problems can be solved without domain discretisation. In addition, the AOIMLS can adaptively avoid singular or ill-conditioned moment matrices, thus ensuring the stability of the calculation results. (2) In the second part (Part III consisting of four chapters), the thermoelastic problems and fracture problems are the main objectives. Due to considering thermal loads, domain integrals appear in the BIEs of the thermoelastic problems, and the expression of integrand functions is known or not depending on the temperature distribution given or not, the AOIMLS enhanced LIBEM is introduced to conduct thermoelasticity analysis thereby. Besides, a series of novel unified boundary integral equations based on BEM and DDM are derived for solving fracture problems and thermoelastic fracture problems in finite and infinite domains. Two sets of unified BIEs are derived for fracture problems in finite and infinite domains based on the direct BEM and DDM respectively, which can provide accurate and stable results. Another two sets of BIEs are addressed by employing indirect BEM and DDM, which cannot ensure a stable result, thereby a modified indirect BEM is proposed which performs much more stable. Moreover, a set of novel BIEs based on the direct BEM and DDM for cracked domains under thermal stress is proposed. (3) In the third part (Part IV consisting of one chapter), a high-efficiency combined BEM and discrete element method (DEM) is proposed to compute the inner stress distribution and particle breakage of particle assemblies based on the solution mapping scheme. For the stress field computation of particles with similar geometry, a template particle is used as the representative particle, so that only the related coefficient matrices of one template particle in the local coordinate system are needed to be calculated, while the coefficient matrices of the other particles, can be obtained by mapping between the local and global coordinate systems. Thus, the combined BEM and DEM is much more effective when modelling a large-scale particle system with a small number of distinct possible particle shapes. Furthermore, with the help of the Hoek-Brown criterion, the possible cracks or breakage paths of a particle can be obtained

ADVANCES IN GEOMETRY INDEPENDENT APPROXIMATIONS

We present recent advances in geometry independent field approximations. The GIFT approach is a generalisation of isogeometric analysis where the approximation used to describe the field variables no-longer has to be identical to the approximation used to describe the geometry of the domain. As such, the geometry can be described using usual CAD representations, e.g. NURBS, which are the most common in the CAD area, whilst local refinement and meshes approximations can be used to describe the field variables, enabling local adaptivity. We show in which cases the approach passes the patch test and present applications to various mechanics, fracture and multi-physics problems

Vademecum-based GFEM (V-GFEM): optimal enrichment for transient problems

This paper proposes a generalized finite element method based on the use of parametric solutions as enrichment functions. These parametric solutions are precomputed off-line and stored in memory in the form of a computational vademecum so that they can be used on-line with negligible cost. This renders a more efficient computational method than traditional finite element methods at performing simulations of processes. One key issue of the proposed method is the efficient computation of the parametric enrichments. These are computed and efficiently stored in memory by employing proper generalized decompositions. Although the presented method can be broadly applied, it is particularly well suited in manufacturing processes involving localized physics that depend on many parameters, such as welding. After introducing the vademecum-generalized finite element method formulation, we present some numerical examples related to the simulation of thermal models encountered in welding processes