Computational Mechanics

In this paper we present a fully-coupled, two-scale homogenization method for dynamic loading in the spirit of FE $$^2$$ methods. The framework considers the balance of linear momentum including inertia at the microscale to capture possible dynamic effects arising from micro heterogeneities. A finite-strain formulation is adapted to account for geometrical nonlinearities enabling the study of e.g. plasticity or fiber pullout, which may be associated with large deformations. A consistent kinematic scale link is established as displacement constraint on the whole representative volume element. The consistent macroscopic material tangent moduli are derived including micro inertia in closed form. These can easily be calculated with a loop over all microscopic finite elements, only applying existing assembly and solving procedures. Thus, making it suitable for standard finite element program architectures. Numerical examples of a layered periodic material are presented and compared to direct numerical simulations to demonstrate the capability of the proposed framework. In addition, a simulation of a split Hopkinson tension test showcases the applicability of the framework to engineering problems.

The contact between bodies is a complex phenomenon that involves mechanical interaction, frictional sliding and heat transfer, among others. A common (and convenient) approach for the mechanical interaction in a finite element framework is to directly use the geometry of the elements to formulate the contact. The main drawback lies in the sharp corners that occur when straight finite elements are connected leading eventually to contact singularities. To circumvent this issue, particularly in the context of beam-to-beam contact, the present work proposes a pointwise contact formulation based on smooth C1 continuous spline contact elements. The proposed spline-based formulation, which can be directly attached to any quadratic beam finite element formulation, guarantees a smooth description for the whole set of elements, where contact takes place. A specific nonlinear normal contact interaction law and a rheological model for friction, both with elastic and viscous damping contributions, are developed increasing robustness in practical applications. To demonstrate this robustness, specific examples are considered including comparisons with a similar surface-to-surface formulation and an alternative smooth contact scheme, smooth contact with finite elements having sharp corners, modeling of a knot tightening with self-contact, and a simulation involving multiple pointwise contacts.

A “quasi-conforming” triangular laminated shell element based on a refined first-order shear deformation theory is presented. The Hu-Washizu variational principle, involving strain and displacement fields as variables, with stresses being considered as Lagrange multipliers, is used to develop the laminate composite shell element. Both strains and displacements are discretized in the element, while displacements alone are discretized at the boundary. The inter-element C 1 continuity is satisfied a posteriori in a weak form. Due to the importance of rotations and shear deformation in the geometrically non-linear analyses of shells, 7 degrees of freedom per node are chosen, viz. three displacements, two first-derivatives in the in-plane directions of the out-of-plane displacement, and two transverse shear strains at each node. To consider the effect of transverse shear deformation on the global behavior of the laminated composite shell, the Reissner-Mindlin first-order theory, with shear correction factors of Chow and Whitney, is adopted. The transverse shear stresses are obtained through the integration of the 3-D equilibrium equations; and the warping induced by transverse shear is considered in the calculation of the in-plane stresses to improve their accuracy. Numerical examples show that the element has good convergence properties and leads to highly accurate stresses.

In patient-specific arterial fluid–structure interaction (FSI) computations the image-based arterial geometry comes from a configuration that is not stress-free. We present a method for estimation of element-based zero-stress (ZS) state. The method has three main components. (1) An iterative method, which starts with an initial guess for the ZS state, is used for computing the element-based ZS state such that when a given pressure load is applied, the image-based target shape is matched. (2) A method for straight-tube geometries with single and multiple layers is used for computing the element-based ZS state so that we match the given diameter and longitudinal stretch in the target configuration and the “opening angle.” (3) An element-based mapping between the arterial and straight-tube configurations is used for mapping from the arterial configuration to the straight-tube configuration, and for mapping the estimated ZS state of the straight tube back to the arterial configuration, to be used as the initial guess for the iterative method that matches the image-based target shape. We present a set of test computations to show how the method works.

It is still a challenge to model swelling phenomena occurring in charged hydrated porous media. This is not only due to the overall complexity of the model but also to the fact that boundary conditions occur, which depend on internal variables. In the present contribution, a multi-component model based on the Theory of Porous Media (TPM) is presented. The advantage of this model is that it is thermodynamically consistent and it consists of only three primary variables. As a result of the boundary conditions depending on internal variables, the numerical treatment within the finite element method (FEM) by use of the mixed finite element scheme reveals artificial oscillations in the numerical results. To overcome these oscillations, we propose to fulfil boundary conditions weakly.

Strain gradient plasticity theories are being widely used for fracture assessment, as they provide a richer description of crack tip fields by incorporating the influence of geometrically necessary dislocations. Characterizing the behavior at the small scales involved in crack tip deformation requires, however, the use of a very refined mesh within microns to the crack. In this work a novel and efficient gradient-enhanced numerical framework is developed by means of the extended finite element method (X-FEM). A mechanism-based gradient plasticity model is employed and the approximation of the displacement field is enriched with the stress singularity of the gradient-dominated solution. Results reveal that the proposed numerical methodology largely outperforms the standard finite element approach. The present work could have important implications on the use of microstructurally-motivated models in large scale applications. The non-linear X-FEM code developed in MATLAB can be downloaded from www.empaneda.com/codes .

Computational homogenization is the gold standard for concurrent multi-scale simulations (e.g., FE2) in scale-bridging applications. Often the simulations are based on experimental and synthetic material microstructures represented by high-resolution 3D image data. The computational complexity of simulations operating on such voxel data is distinct. The inability of voxelized 3D geometries to capture smooth material interfaces accurately, along with the necessity for complexity reduction, has motivated a special local coarse-graining technique called composite voxels (Kabel et al. Comput Methods Appl Mech Eng 294: 168–188, 2015). They condense multiple fine-scale voxels into a single voxel, whose constitutive model is derived from the laminate theory. Our contribution generalizes composite voxels towards composite boxels (ComBo) that are non-equiaxed, a feature that can pay off for materials with a preferred direction such as pseudo-uni-directional fiber composites. A novel image-based normal detection algorithm is devised which (i) allows for boxels in the firsts place and (ii) reduces the error in the phase-averaged stresses by around 30% against the orientation cf. Kabel et al. (Comput Methods Appl Mech Eng 294: 168–188, 2015) even for equiaxed voxels. Further, the use of ComBo for finite strain simulations is studied in detail. An efficient and robust implementation is proposed, featuring an essential selective back-projection algorithm preventing physically inadmissible states. Various examples show the efficiency of ComBo against the original proposal by Kabel et al. (Comput Methods Appl Mech Eng 294: 168–188, 2015) and the proposed algorithmic enhancements for nonlinear mechanical problems. The general usability is emphasized by examining various Fast Fourier Transform (FFT) based solvers, including a detailed description of the Doubly-Fine Material Grid (DFMG) for finite strains. All of the studied schemes benefit from the ComBo discretization.