3D architected isotropic materials with tunable stiffness and buckling strength
Introduction
Materials with extreme mechanical properties are highly attractive for many applications. Among them, stiffness and strength are fundamental for determining material load-bearing capability. Stiffness accounts for the ability to resist deformation while strength measures the ultimate load-carrying capability.
Many studies have been devoted to exploring materials with optimal stiffness (Hashin, 1962, Francfort and Murat, 1986, Milton, 1986, Sigmund, 1994, Sigmund, 2000, Berger et al., 2017, Tancogne-Dejean et al., 2018, Wang et al., 2019). It has been shown that isotropic stiffness-optimal plate materials can achieve the Hashin–Shtrikman upper bounds on isotropic elastic stiffness and show up to three times the stiffness of the isotropic stiffness-optimal truss materials in the low volume fraction limit (Berger et al., 2017, Tancogne-Dejean et al., 2018, Christensen, 1986, Sigmund et al., 2016). The stiffness advantage is attributed to the multiaxial stiffness offered by the constituent plates while bars in the truss materials only offer axial stiffness. However, a recent study has shown that isotropic stiffness-optimal truss material is superior from a bucking strength point of view compared to the isotropic plate counterparts for the same volume fraction when the volume fraction is below 31% (Andersen et al., 2021a) because of higher bending stiffness associated with the bars than the constituent plates. It is for example shown that the isotropic stiffness-optimal truss material possesses 48% higher uniaxial buckling strength and 52% lower Young’s modulus than the isotropic plate material for a volume fraction of 20%.
Material geometry strongly dictates material properties. Novel materials with exotic properties have been achieved through careful tailoring of material geometries via different design approaches. Among them, topology optimization methods (Bendsøe and Sigmund, 2003) have been proven powerful tools in designing novel materials ranging from mechanical properties, such as optimal stiffness (Sigmund, 1994, Andreassen et al., 2014), auxetic behavior (Sigmund, 1994, Andreassen et al., 2014, Wang, 2018), to acoustic and optical properties (Christiansen and Sigmund, 2016, Wang et al., 2018). Regarding elastic stability, material buckling failure may develop at different scales spanning from highly localized short wavelength modes to long wavelength modes. Previous numerical studies have employed homogenization methods assuming separation of scales (Guedes and Kikuchi, 1990) and Bloch–Floquet theory for detecting short and long wavelength buckling (Geymonat et al., 1993). A general methodology for characterizing material strength due to bifurcation failure was proposed in Triantafyllidis and Schnaidt (1993). Assuming small strains and ignoring material and geometric nonlinearities, topology optimization of material strength was first studied in Neves et al. (2002a), where only cell-periodic buckling modes were taken into account. Later, this work was extended to cover both local and global modes via the Bloch–Floquet theory (Neves et al., 2002b). More recently, 2D materials with enhanced buckling strength (Thomsen et al., 2018) have been systematically designed using topology optimization methods for different macro-level stress situations based on the homogenization theory and linear buckling analysis (LBA) with Bloch–Floquet boundary conditions. It was shown that optimized first-order hierarchical materials outperform their non-hierarchical counterparts in terms of buckling strength at the cost of slightly decreased stiffness. Further material evaluations considering both geometrical and material nonlinearity have proven that the superior buckling strength of the optimized hierarchical materials compared to reference materials also hold for finite structures (Wang and Sigmund, 2020) and geometrically nonlinear modeling (Bluhm et al., 2020). Hence, a buckling optimization approach assuming small strains has been demonstrated to be efficient and practically relevant, even for nonlinear structures.
This study extends the work in 2D material design with enhanced buckling strength in Thomsen et al. (2018) to 3D material design with tunable stiffness and buckling response utilizing a flexible framework for large scale topology optimization based on the Portable, Extensible Toolkit for Scientific Computation (PETSc) (Balay et al., 2016, Aage et al., 2015, Wang, 2018) and the Scalable Library for Eigenvalue Problem Computations (SLEPc) (Hernandez et al., 2005). Considered materials are constrained to be cubic symmetric and elastically isotropic. The homogenization method is employed to characterize the effective material properties, and LBA, together with Bloch–Floquet boundary conditions evaluated over the boundaries of the irreducible Brillouin zone (IBZ) (Brillouin, 1953), is employed to evaluate material buckling strength. The optimization problem for designing materials with tunable stiffness and buckling strength is formulated to minimize the weighted stiffness and buckling response. 3D materials are designed to achieve tunable stiffness and buckling response by assigning different weight factors for stiffness and buckling strength under uniaxial compression. Moreover, inspired by the topology optimized material configurations, a subsequent feature-based shape optimization approach is employed to simplify material geometries. In the feature-based approach, material architectures are parameterized using several hollow and one solid super-ellipsoids (Wang, 2018, Wein et al., 2020).
The paper is organized as follows. Section 2 presents finite element formulations to evaluate stiffness and buckling strength, and formulates the optimization problem for designing materials with tunable properties. Section 3 first validates the proposed optimization formulation by optimizing a material microstructure for maximum buckling strength under hydrostatic compression. Topology optimized single-length scale material microstructures with enhanced stiffness and strength under uniaxial compression are then systematically designed. Inspired by the optimized microstructure configurations, a shape optimization scheme is proposed to simplify the optimized microstructure geometries further. The paper ends with the conclusions in Section 4.
Section snippets
Optimization formulation of 3D material design with enhanced stiffness and buckling strength
This section presents the essential formulations for designing materials with enhanced stiffness and buckling strength using topology optimization. The finite element method is combined with homogenization theory used to evaluate material properties (Cook et al., 2002). It is well-known that linear elements, i.e., 8-node hexahedral element (), overestimate material stiffness and suffer from shear locking, which results in an inaccurate representation of stresses. To represent the stress
Results
The proposed optimization formulation is employed to design 3D isotropic microstructures with enhanced stiffness and buckling strength. The unit cell is discretized by 64 × 64 × 64 elements. The one-case robust formulation is employed with and a filter radius of , which corresponds to a relative minimal length scale around 0.02, i.e., the minimal feature size is around 2% of the microstructure size (da Silva et al., 2021). The volume fraction upper bound on the intermediate
Conclusion
3D isotropic microstructures with tunable stiffness and buckling strength have systematically been designed. The effective material properties are evaluated using homogenization, and linear buckling analysis is employed to predict the material buckling strength for a given macroscopic stress state, where both microscopic and macroscopic failure modes are captured using Bloch–Floquet boundary conditions. The optimization problem is formulated to improve the weighted stiffness and bucking
CRediT authorship contribution statement
Fengwen Wang: Conceptualization, Methodology, Software, Investigation, Formal analysis, Writing - original draft, Writing - review & editing. O. Sigmund: Conceptualization, Methodology, Formal analysis, Writing - review & editing, Funding acquisition.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgments
We acknowledge the financial support from the Villum Fonden, Denmark through the Villum Investigator Project InnoTop. We further acknowledge valuable discussions with Morten N. Andersen at the Department of Mechanical Engineering in Technical University of Denmark.
References (45)
- et al.
On the competition for ultimately stiffness and strong architected materials
Mater. Des.
(2021) - et al.
Design of manufacturable 3d extremal elastic microstructure
Mech. Mater.
(2014) - et al.
Nonlinear compressive stability of hyperelastic 2d lattices at finite volume fractions
J. Mech. Phys. Solids
(2020) Mechanics of low density materials
J. Mech. Phys. Solids
(1986)- et al.
Preprocessing and postprocessing for materials based on the homogenization method with adaptive finite element methods
Comput. Methods Appl. Mech. Engrg.
(1990) - et al.
A review of homogenization and topology optimization I - homogenization theory for media with periodic structure
Comput. Struct.
(1998) - et al.
Systematic control design by optimizing a vector performance index
Comput. Aided Design Control Syst.
(1980) Materials with prescribed constitutuve parameters: An inverse homogenization problem
Int. J. Solids Struct.
(1994)A new class of extremal composites
J. Mech. Phys. Solids
(2000)- et al.
Buckling strength topology optimization of 2d periodic materials based on linearized bifurcation analysis
Comput. Methods Appl. Mech. Engrg.
(2018)
Comparison of microscopic and macroscopic instabilities in a class of two-dimensional periodic composites
J. Mech. Phys. Solids
Systematic design of 3d auxetic lattice materials with programmable Poisson’s ratio for finite strains
J. Mech. Phys. Solids
Simple optimal lattice structures for arbitrary loadings
Extreme Mech. Lett.
Use of incompatible displacement modes for the calculation of element stiffnesses or stresses
Finite Elem. Anal. Des.
Incompatible Displacement Models
Topology optimization using petsc: An easy-to-use, fully parallel, open source topology optimization framework
Struct. Multidiscip. Optim.
Paraview: An end-user tool for large data visualization
Buckling and yield strength estimation of architected materials under arbitrary loads
PETSc Users ManualTechnical Report ANL-95/11 - Revision 3.7
Material interpolation schemes in topology optimization
Arch. Appl. Mech.
Topology Optimization: Theory, Methods and Applications
Mechanical metamaterials at the theoretical limit of isotropic elastic stiffness
Nature
Cited by (20)
Efficient inverse-designed structural infill for complex engineering structures
2024, Thin-Walled StructuresRational designs of mechanical metamaterials: Formulations, architectures, tessellations and prospects
2023, Materials Science and Engineering R: ReportsWang tiles enable combinatorial design and robot-assisted manufacturing of modular mechanical metamaterials
2023, Extreme Mechanics LettersLattice structure design optimization under localized linear buckling constraints
2023, Computers and StructuresProgramming and physical realization of extreme three-dimensional responses of metastructures under large deformations
2023, International Journal of Engineering Science