Peter Wriggers

Vien Minh Nguyen-Thanh, Xiaoying Zhuang, Hung Nguyen-Xuan et al.

This paper presents the Virtual Element Method (VEM) for the modeling of crack propagation in 2D within the context of linear elastic fracture mechanics (LEFM). By exploiting the advantage of mesh flexibility in the VEM, we establish an adaptive mesh refinement strategy based on the superconvergent patch recovery for triangular, quadrilateral as well as for arbitrary polygonal meshes. For the local stiffness matrix in VEM, we adopt a stabilization term which is stable for both isotropic scaling and ratio. Stress intensity factors (SIFs) of a polygonal mesh are discussed and solved by using the interaction domain integral. The present VEM formulations are finally tested and validated by studying its convergence rate for both continuous and discontinuous problems, and are compared with the optimal convergence rate in the conventional Finite Element Method (FEM). Furthermore, the adaptive mesh refinement strategies used to effectively predict the crack growth with the existence of hanging nodes in nonconforming elements are examined.

Nima Noii, Fadi Aldakheel, Thomas Wick et al.

This work addresses an efficient Global-Local approach supplemented with predictor-corrector adaptivity applied to anisotropic phase-field brittle fracture. The phase-field formulation is used to resolve the sharp crack surface topology on the anisotropic/non-uniform local state in the regularized concept. To resolve the crack phase-field by a given single preferred direction, second-order structural tensors are imposed to both the bulk and crack surface density functions. Accordingly, a split in tension and compression modes in anisotropic materials is considered. A Global-Local formulation is proposed, in which the full displacement/phase-field problem is solved on a lower (local) scale, while dealing with a purely linear elastic problem on an upper (global) scale. Robin-type boundary conditions are introduced to relax the stiff local response at the global scale and enhancing its stabilization. Another important aspect of this contribution is the development of an adaptive Global-Local approach, where a predictor-corrector scheme is designed in which the local domains are dynamically updated during the computation. To cope with different finite element discretizations at the interface between the two nested scales, a non-matching dual mortar method is formulated. Hence, more regularity is achieved on the interface. Several numerical results substantiate our developments.

Carsten Carstensen, Philipp Bringmann, Friederike Hellwig et al.

The discontinuous Petrov-Galerkin method is a minimal residual method with broken test spaces and is introduced for a nonlinear model problem in this paper. Its lowest-order version applies to a nonlinear uniformly convex model example and is equivalently characterized as a mixed formulation, a reduced formulation, and a weighted nonlinear least-squares method. Quasi-optimal a priori and reliable and efficient a posteriori estimates are obtained for the abstract nonlinear dPG framework for the approximation of a regular solution. The variational model example allows for a built-in guaranteed error control despite inexact solve. The subtle uniqueness of discrete minimizers is monitored in numerical examples.

Jan Niklas Fuhg, Christoph Boehm, Nikolaos Bouklas et al.

Computational multiscale methods for analyzing and deriving constitutive responses have been used as a tool in engineering problems because of their ability to combine information at different length scales. However, their application in a nonlinear framework can be limited by high computational costs, numerical difficulties, and/or inaccuracies. In this paper, a hybrid methodology is presented which combines classical constitutive laws (model-based), a data-driven correction component, and computational multiscale approaches. A model-based material representation is locally improved with data from lower scales obtained by means of a nonlinear numerical homogenization procedure leading to a model-data-driven approach. Therefore, macroscale simulations explicitly incorporate the true microscale response, maintaining the same level of accuracy that would be obtained with online micro-macro simulations but with a computational cost comparable to classical model-driven approaches. In the proposed approach, both model and data play a fundamental role allowing for the synergistic integration between a physics-based response and a machine learning black-box. Numerical applications are implemented in two dimensions for different tests investigating both material and structural responses in large deformation.

Dengpeng Huang, Jan Niklas Fuhg, Christian Weißenfels et al.

Data-driven material models have many advantages over classical numerical approaches, such as the direct utilization of experimental data and the possibility to improve performance of predictions when additional data is available. One approach to develop a data-driven material model is to use machine learning tools. These can be trained offline to fit an observed material behaviour and then be applied in online applications. However, learning and predicting history dependent material models, such as plasticity, is still challenging. In this work, a machine learning based material modelling framework is proposed for both elasticity and plasticity. The machine learning based hyperelasticity model is developed with the Feed forward Neural Network (FNN) directly whereas the machine learning based plasticity model is developed by using of a novel method called Proper Orthogonal Decomposition Feed forward Neural Network (PODFNN). In order to account for the loading history, the accumulated absolute strain is proposed to be the history variable of the plasticity model. Additionally, the strain-stress sequence data for plasticity is collected from different loading-unloading paths based on the concept of sequence for plasticity. By means of the POD, the multi-dimensional stress sequence is decoupled leading to independent one dimensional coefficient sequences. In this case, the neural network with multiple output is replaced by multiple independent neural networks each possessing a one-dimensional output, which leads to less training time and better training performance. To apply the machine learning based material model in finite element analysis, the tangent matrix is derived by the automatic symbolic differentiation tool AceGen. The effectiveness and generalization of the presented models are investigated by a series of numerical examples using both 2D and 3D finite element analysis.