The objective of this contribution is to present a unifying review on strain-driven computational homogenization at finite strains, thereby elaborating on computational aspects of the finite element method. The underlying assumption of computational homogenization is separation of length scales, and hence, computing the material response at the macroscopic scale from averaging the microscopic behavior. In doing so, the energetic equivalence between the two scales, the Hill-Mandel condition, is guaranteed via imposing proper boundary conditions such as linear displacement, periodic displacement and antiperiodic traction, and constant traction boundary conditions. Focus is given on the finite element implementation of these boundary conditions and their influence on the overall response of the material. Computational frameworks for all canonical boundary conditions are briefly formulated in order to demonstrate similarities and differences among the various boundary conditions. Furthermore, we detail on the computational aspects of the classical Reuss' and Voigt's bounds and their extensions to finite strains. A concise and clear formulation for computing the macroscopic tangent necessary for FE 2 calculations is presented. The performances of the proposed schemes are illustrated via a series of two-and three-dimensional numerical examples. The numerical examples provide enough details to serve as benchmarks.
Although both computational and analytical homogenization are well-established today, a thorough and systematic study to compare them is missing in the literature. This manuscript aims to provide an exhaustive comparison of numerical computations and analytical estimates, such as Voigt, Reuss, Hashin–Shtrikman, and composite cylinder assemblage. The numerical computations are associated with canonical boundary conditions imposed on either tetragonal, hexagonal, or circular representative volume elements using the finite-element method. The circular representative volume element is employed to capture an effective isotropic material response suitable for comparison with associated analytical estimates. The analytical results from composite cylinder assemblage are in excellent agreement with the numerical results obtained from a circular representative volume element. We observe that the circular representative volume element renders identical responses for both linear displacement and periodic boundary conditions. In addition, the behaviors of periodic and random microstructures with different inclusion distributions are examined under various boundary conditions. Strikingly, for some specific microstructures, the effective shear modulus does not lie within the Hashin–Shtrikman bounds. Finally, numerical simulations are carried out at finite deformations to compare different representative volume element types in the nonlinear regime. Unlike other canonical boundary conditions, the uniform traction boundary conditions result in nearly identical effective responses for all types of representative volume element, indicating that they are less sensitive with respect to the underlying microstructure. The numerical examples furnish adequate information to serve as benchmarks.
In the past decades computational homogenization has proven to be a powerful strategy to compute the overall response of continua. Central to computational homogenization is the Hill-Mandel condition. The Hill-Mandel condition is fulfilled via imposing displacement boundary conditions (DBC), periodic boundary conditions (PBC) or traction boundary conditions (TBC) collectively referred to as canonical boundary conditions. While DBC and PBC are widely implemented, TBC remains poorly understood, with a few exceptions. The main issue with TBC is the singularity of the stiffness matrix due to rigid body motions. The objective of this manuscript is to propose a generic strategy to implement TBC in the context of computational homogenization at finite strains. To eliminate rigid body motions, we introduce the concept of semi-Dirichlet boundary conditions. Semi-Dirichlet boundary conditions are non-homogeneous Dirichlet-type constraints that simultaneously satisfy the Neumann-type conditions. A key feature of the proposed methodology is its applicability for both strain-driven as well as stress-driven homogenization. The The term semi-Dirichlet exists in mathematics but, in the completely different context of semi-Dirichlet forms and should not be confused with our non-homogeneous Dirichlet-type constraints.
Finite-thickness interphases between different constituents in heterogeneous materials are often replaced by a zerothickness interface model. Commonly accepted interface models intuitively assume that the interface layer is situated exactly in the middle of its associated interphase. Furthermore, it has been reported in the literature that this assumption is necessary to guarantee the balance of angular momentum on the interface. While the interface coincides with the mid-layer of a uniform interphase, we argue that this assumption fails to sufficiently capture the behavior of graded or inhomogeneous interphases. This contribution extends the formulation of the general interface model to account for arbitrary interface positions. The issue of angular momentum balance on general interfaces is critically revisited. It is proven that the interface position does not necessarily have to coincide with the mid-layer in order to satisfy the angular momentum balance. The analysis here leads to a unique definition of the controversially discussed interface configuration. The presented general interface model is essentially based upon the weighted average operator instead of the commonly accepted classical average operator. The framework is geometrically exact and suitable for finite deformations. The significance of the interface position is demonstrated via a series of examples where the interface position is identified based on a full resolution interphase.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.