Supervisor:
Project members:
- M.Sc. Meysam Joulaian
Project description:
The numerical modeling of multi-scale structural problems with heterogeneities is a numerically demanding task. The difficulty lies in the existence of size-dependent phenomena requiring an appropriate solution. The fact that problems of this nature usually feature a complex geometry, which calls for the generation of complex mesh models, is an additional complication. There are various techniques for circumventing the difficulties posed by multi-scale simulation, such as the Generalized Finite Element Method (GFEM), the eXtended Finite Element Method (XFEM) or the Meshless Methods (MM). For this particular project we chose to employ the recently proposed Finite Cell Method (FCM), which is a combination of the high-order finite element method (p-FEM) coupled with a fictitious domain approach.
The FCM has an inherent capacity for facilitating the handling of linear and non-linear problems connected with complex geometry and it is a promising method for many intricate applications. Despite the significant success of the method, however, the FCM is still in the developing stages and needs to be improved on several levels. In the case of multi-material modeling, for instance, where the solution features a loss of regularity along the material interface, the FCM exhibits oscillatory results particularly in the derivatives. Besides, the loss of regularity hinders optimal convergence. A simple, yet effective remedy is to enrich the finite element space with an appropriate enrichment term designed to represent the weak discontinuity at the material interface. This enrichment can be provided either by choosing the pre-known functions that can describe the special phenomenon in the domain of interest or by a local h-, p- or hp-extension. Irrespective to the type of enrichment, it can be added to the FE space by means of the hp-d method which places the enrichment terms on a superimposed mesh, the so- called overlay mesh. Figure 1 schematically demonstrates the main idea behind this approach. This method of enrichment complies effectively with both the FCM framework and the multi-scale simulation.
Figure 1: FCM with the local enrichment using the hp-d method
Publications:
- A. Düster, E. Rank. “Die Finite Cell Methode – Eine Fictitious Domain Methode mit Finite-Element Ansätzen hoher Ordnung.” GAMM Rundbrief, 2/2011:6-13, 2011. pp. 1755-1762.
- D. Schillinger, A. Düster, E. Rank. “The hp-d-adaptive finite cell method for geometrically nonlinear problems of solid mechanics.” International Journal for Numerical Methods in Engineering, 89:1171-1202, 2012.
- M. Joulaian, A. Düster. “Local enrichment of the finite cell method for problems with material interfaces.” Computational Mechanics, 52:741-762, 2013.
- S. Duczek, M. Joulaian, A. Düster, U. Gabbert. “Simulation of Lamb waves using the spectral cell method.”Health Monitoring of Structural and Biological Systems, 86951U, DOI: 10.1117/12.2009983 , 2013.
- M. Joulaian, A. Düster. “The hp-d version of the finite cell method with local enrichment for multiscale problems.” PAMM, 13: 259-260, 2013.
- S. Duczek, M. Joulaian, A. Düster, U. Gabbert. ” Numerical analysis of Lamb waves using the finite and spectral cell method.”;International Journal for Numerical Methods in Engineering, 99:26-53, 2014.
- M. Joulaian, S. Duczek, U. Gabbert, A. Düster. ” Finite and spectral cell method for wave propagation in heterogeneous materials”; Computational Mechanics, 54:661-675, 2014.
- M. Joulaian, S. Duczek, U. Gabbert, A. Düster. “Efficient simulation of wave propagation in heterogeneous materials”; PAMM, 14:715-716, 2014.
- S. Heinze, M. Joulaian, A. Düster. “Numerical homogenization of hybrid metal foams using the finite cell method”; Computers & Mathematics with Applications, 70:1501-1517, 2015.
- M. Joulaian, N. Zander, T. Bog, S. Kollmannsberger, E. Rank, A Düster. “A high-order enrichment strategy for the finite cell method”; PAMM, 15:207-208, 2015.
- M. Joulaian, S. Hubrich, A. Düster. “Numerical integration of discontinuities on arbitrary domains based on moment fitting”; Computational Mechanics, 57:979-999, 2016.