alt

ea_back


FCM for Geometrically and Materially Nonlinear Analysis of Metal Foams

 

Project team leader Dr. rer. nat. Ralf-Peter Mundani
Stefan Kollmannsberger, M.Sc.
Prof. Dr. rer. nat. Ernst Rank

PhD student Dipl.-Ing. Dominik Schillinger, M.Sc.

 

Project Description

The conventional finite element method requires the discretization of the domain of interest into a finite element mesh, whose boundaries have to coincide with the physical boundaries of the problem. While this constraint can be easily achieved for most applications of solid mecha- nics, it becomes a huge problem for structures with very complex boundaries, whose geometry is only available in voxel format (see Fig. 1). Mesh generation algorithms that are able to transfer voxel based geometry discriptions into finite element meshes are error prone, often yield largely distorted elements and are computationally very expensive, so that usually more than 90% of the total FE analysis time is consumed by mesh generation alone.

The recently introduced Finite Cell Method (FCM) is targeted at circumventing the complete mesh generation. Its main feature is the combination of high-order Legendre based finite elements by a fictitious domain approach, which allows the extension of the physical domain beyond its boundaries. It complements the weak formulation of the mechanical problem by an additional function α, which is 1 inside the physical domain and 0 outside. Thus, the boundary information is implicitly incorporated into the functional and does not need to be represented by the mesh. Since the extension of solution fields are smooth beyond the physical domain, FCM can preserve exponential convergence of conventionl p-FEM.

 

alt


 

 

 

 

 

 

 

 

 

 

The FCM strategy is briefly illustrated in Figs. 2 through 4 by the simple benchmark problem of a perforated plate. Due to its fictitious domain capability, the plate can be discretized irrespective of the circular hole in the middle by just 4 Legendre elements (p=10). The accuracy of the method significantly depends on capturing the discontinuity of α along the physical boundary during numerical integration (α=0 inside the circular hole, α=1 in the plate). This is achieved by an adaptive partitioning with sub-cells, which each contain a grid of 10x10 Gauss points.

Up to now, FCM has been successfully applied in the linear elastic analysis of complex voxel-based structures, such as the proximal femur bone in Fig. 5. However, porous structures at higher load levels respond strongly nonlinear and can therefore not be described accurately by linearized theories. In the context of metal foams, these nonlinearities include most importantly local buckling of cell walls, plastification of parts of the physical domain and self-contact of cell walls during the loading process. The main target of the present project is the integration of nonlinear capabilities into the existing linear FCM approach, so that an elasto-plastic deformation process at finite strains can be realistically simulated for the aluminium foam example.



alt

 

 

 

 

 

 

 

 

 

 

Fig. 5: Linear elastic FCM analysis of a femur (voxel geometry, high-order 3D hex elements, von Mises stress).

 

Publications

  1. Düster, A.; Parvizian, J.; Yang, Z,; Rank, E.;
    The FInite Cell Method for Three-Dimensional Problems of Solid Mechanics. Computer Methods of Applied Mechanics and Engineering 197, 3768 - 3782, 2008.
  2. Parvizian, J.; Düster, A.; Rank, E.;
    Finite Cell Method: h- and p- extension for embedded domain methods in solid mechanics. Computational Mechanics 41, 122-133, 2007.

Contact: Dipl.-Ing. Dominik Schillinger, M.Sc.

Fig. 1: CT based voxel model of a section of an aluminium foam used in lightweight panels of aerospace structures (courtesy of Prof. Banhart, Helmholtz Zentrum Berlin)
Fig. 2: Perforated plate benchmark (FCM mesh and boundary conditions)


Fig 3: Adaptive sub-cell structure for numerical integration
Fig. 4: Normal stress results in y-direction