Research Papers

Optimization of Composite Fracture Properties: Method, Validation, and Applications

[+] Author and Article Information
Grace X. Gu

Laboratory for Atomistic and
Molecular Mechanics (LAMM),
Department of Mechanical Engineering,
Massachusetts Institute of Technology,
77 Massachusetts Avenue,
Cambridge, MA 02139
e-mail: gracegu@mit.edu

Leon Dimas

Laboratory for Atomistic and
Molecular Mechanics (LAMM),
Department of Civil and
Environmental Engineering,
Massachusetts Institute of Technology,
77 Massachusetts Avenue,
Cambridge, MA 02139
e-mail: leon_dim@mit.edu

Zhao Qin

Laboratory for Atomistic and
Molecular Mechanics (LAMM),
Department of Civil and
Environmental Engineering,
Massachusetts Institute of Technology,
77 Massachusetts Avenue,
Cambridge, MA 02139
e-mail: qinzhao@mit.edu

Markus J. Buehler

Laboratory for Atomistic and
Molecular Mechanics (LAMM),
Department of Civil and
Environmental Engineering,
Massachusetts Institute of Technology,
77 Massachusetts Avenue,
Cambridge, MA 02139
e-mail: mbuehler@mit.edu

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received March 22, 2016; final manuscript received April 7, 2016; published online May 5, 2016. Editor: Yonggang Huang.

J. Appl. Mech 83(7), 071006 (May 05, 2016) (7 pages) Paper No: JAM-16-1149; doi: 10.1115/1.4033381 History: Received March 22, 2016; Revised April 07, 2016

A paradigm in nature is to architect composites with excellent material properties compared to its constituents, which themselves often have contrasting mechanical behavior. Most engineering materials sacrifice strength for toughness, whereas natural materials do not face this tradeoff. However, biology's designs, adapted for organism survival, may have features not needed for some engineering applications. Here, we postulate that mimicking nature's elegant use of multimaterial phases can lead to better optimization of engineered materials. We employ an optimization algorithm to explore and design composites using soft and stiff building blocks to study the underlying mechanisms of nature's tough materials. For different applications, optimization parameters may vary. Validation of the algorithm is carried out using a test suite of cases without cracks to optimize for stiffness and compliance individually. A test case with a crack is also performed to optimize for toughness. The validation shows excellent agreement between geometries obtained from the optimization algorithm and the brute force method. This study uses different objective functions to optimize toughness, stiffness and toughness, and compliance and toughness. The algorithm presented here can provide researchers a way to tune material properties for a vast number of engineering problems by adjusting the distribution of soft and stiff materials.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Sih, G. , 1973, “ Some Basic Problems in Fracture Mechanics and New Concepts,” Eng. Fract. Mech., 5(2), pp. 365–377. [CrossRef]
Gao, H. , Ji, B. , Jäger, I. L. , Arzt, E. , and Fratzl, P. , 2003, “ Materials Become Insensitive to Flaws at Nanoscale: Lessons From Nature,” Proc. Natl. Acad. Sci., 100(10), pp. 5597–5600. [CrossRef]
Bonderer, L. J. , Studart, A. R. , and Gauckler, L. J. , 2008, “ Bioinspired Design and Assembly of Platelet Reinforced Polymer Films,” Science, 319(5866), pp. 1069–1073. [CrossRef] [PubMed]
Sen, D. , and Buehler, M. J. , 2011, “ Structural Hierarchies Define Toughness and Defect-Tolerance Despite Simple and Mechanically Inferior Brittle Building Blocks,” Sci. Rep., 1, p. 35.
Meyers, M. A. , Chen, P.-Y. , Lin, A. Y.-M. , and Seki, Y. , 2008, “ Biological Materials: Structure and Mechanical Properties,” Prog. Mater. Sci., 53(1), pp. 1–206. [CrossRef]
Yang, W. , Sherman, V. R. , Gludovatz, B. , Mackey, M. , Zimmermann, E. A. , Chang, E. H. , Schaible, E. , Qin, Z. , Buehler, M. J. , and Ritchie, R. O. , 2014, “ Protective Role of Arapaima Gigas Fish Scales: Structure and Mechanical Behavior,” Acta Biomater., 10(8), pp. 3599–3614. [CrossRef] [PubMed]
Launey, M. E. , Buehler, M. J. , and Ritchie, R. O. , 2010, “ On the Mechanistic Origins of Toughness in Bone,” Annu. Rev. Mater. Res., 40(1), pp. 25–53. [CrossRef]
Dunlop, J. W. , and Fratzl, P. , 2010, “ Biological Composites,” Annu. Rev. Mater. Res., 40(1), pp. 1–24. [CrossRef]
Mayer, G. , 2005, “ Rigid Biological Systems as Models for Synthetic Composites,” Science, 310(5751), pp. 1144–1147. [CrossRef] [PubMed]
Fratzl, P. , 2007, “ Biomimetic Materials Research: What Can We Really Learn From Nature's Structural Materials?” J. R. Soc. Interface, 4(15), pp. 637–642. [CrossRef] [PubMed]
Wegst, U. , and Ashby, M. , 2004, “ The Mechanical Efficiency of Natural Materials,” Philos. Mag., 84(21), pp. 2167–2186. [CrossRef]
Sanchez, C. , Arribart, H. , and Guille, M. M. G. , 2005, “ Biomimetism and Bioinspiration as Tools for the Design of Innovative Materials and Systems,” Nat. Mater., 4(4), pp. 277–288. [CrossRef] [PubMed]
Almeida, F. S. , and Awruch, A. M. , 2009, “ Design Optimization of Composite Laminated Structures Using Genetic Algorithms and Finite Element Analysis,” Compos. Struct., 88(3), pp. 443–454. [CrossRef]
Naik, G. N. , Gopalakrishnan, S. , and Ganguli, R. , 2008, “ Design Optimization of Composites Using Genetic Algorithms and Failure Mechanism Based Failure Criterion,” Compos. Struct., 83(4), pp. 354–367. [CrossRef]
Paluch, B. , Grediac, M. , and Faye, A. , 2008, “ Combining a Finite Element Programme and a Genetic Algorithm to Optimize Composite Structures With Variable Thickness,” Compos. Struct., 83(3), pp. 284–294. [CrossRef]
Muc, A. , and Gurba, W. , 2001, “ Genetic Algorithms and Finite Element Analysis in Optimization of Composite Structures,” Compos. Struct., 54(2–3), pp. 275–281. [CrossRef]
Sigmund, O. , and Torquato, S. , 1997, “ Design of Materials With Extreme Thermal Expansion Using a Three-Phase Topology Optimization Method,” J. Mech. Phys. Solids, 45(6), pp. 1037–1067. [CrossRef]
Wang, M. Y. , Wang, X. , and Guo, D. , 2003, “ A Level Set Method for Structural Topology Optimization,” Comput. Methods Appl. Mech. Eng., 192(1), pp. 227–246. [CrossRef]
Rozvany, G. I. , 2009, “ A Critical Review of Established Methods of Structural Topology Optimization,” Struct. Multidiscip. Optim., 37(3), pp. 217–237. [CrossRef]
Hassani, B. , and Hinton, E. , 2012, Homogenization and Structural Topology Optimization: Theory, Practice and Software, Springer Science & Business Media, London.
Sigmund, O. , 2001, “ A 99 Line Topology Optimization Code Written in matlab,” Struct. Multidiscip. Optim., 21(2), pp. 120–127. [CrossRef]
Hajela, P. , Lee, E. , and Lin, C.-Y. , 1993, “ Genetic Algorithms in Structural Topology Optimization,” Topology Design of Structures, Springer, Dordrecht, The Netherlands, pp. 117–133.
Sigmund, O. , and Petersson, J. , 1998, “ Numerical Instabilities in Topology Optimization: A Survey on Procedures Dealing With Checkerboards, Mesh-Dependencies and Local Minima,” Struct. Optim., 16(1), pp. 68–75. [CrossRef]
Guo, X. , and Gao, H. , 2006, “ Bio-Inspired Material Design and Optimization,” Proceedings of the IUTAM Symposium on Topological Design Optimization of Structures, Machines and Materials, Springer, Dordrecht, The Netherlands, pp. 439–453.
Challis, V. J. , Roberts, A. P. , and Wilkins, A. H. , 2008, “ Fracture Resistance Via Topology Optimization,” Struct. Multidiscip. Optim., 36(3), pp. 263–271. [CrossRef]
Zhang, Z. , Schwartz, S. , Wagner, L. , and Miller, W. , 2000, “ A Greedy Algorithm for Aligning DNA Sequences,” J. Comput. Biol., 7(1–2), pp. 203–214. [CrossRef] [PubMed]
Ruiz, R. , and Stützle, T. , 2007, “ A Simple and Effective Iterated Greedy Algorithm for the Permutation Flowshop Scheduling Problem,” Eur. J. Oper. Res., 177(3), pp. 2033–2049. [CrossRef]
Dunstan, F. , and Welsh, D. , 1973, “ A Greedy Algorithm for Solving a Certain Class of Linear Programmes,” Math. Program., 5(1), pp. 338–353. [CrossRef]
Irwin, G. R. , 1997, “ Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate,” SPIE Milestone Series MS, 137, pp. 167–170.
Lawn, B. , 1993, Fracture of Brittle Solids, Cambridge University Press, Melbourne, Australia.
Anderson, T. L. , and Anderson, T. , 2005, Fracture Mechanics: Fundamentals and Applications, CRC Press, Boca Raton, FL.
Dimas, L. S. , Bratzel, G. H. , Eylon, I. , and Buehler, M. J. , 2013, “ Tough Composites Inspired by Mineralized Natural Materials: Computation, 3D Printing, and Testing,” Adv. Funct. Mater., 23(36), pp. 4629–4638. [CrossRef]
Dimas, L. S. , and Buehler, M. J. , 2014, “ Modeling and Additive Manufacturing of Bio-Inspired Composites With Tunable Fracture Mechanical Properties,” Soft Matter, 10(25), pp. 4436–4442. [CrossRef] [PubMed]
Finnemore, A. , Cunha, P. , Shean, T. , Vignolini, S. , Guldin, S. , Oyen, M. , and Steiner, U. , 2012, “ Biomimetic Layer-by-Layer Assembly of Artificial Nacre,” Nat. Commun., 3, p. 966. [CrossRef] [PubMed]
Tang, Z. Y. , Kotov, N. A. , Magonov, S. , and Ozturk, B. , 2003, “ Nanostructured Artificial Nacre,” Nat. Mater., 2(6), pp. U413–U418. [CrossRef]
Munch, E. , Launey, M. E. , Alsem, D. H. , Saiz, E. , Tomsia, A. P. , and Ritchie, R. O. , 2008, “ Tough, Bio-Inspired Hybrid Materials,” Science, 322(5907), pp. 1516–1520. [CrossRef] [PubMed]
Launey, M. E. , Munch, E. , Alsem, D. H. , Saiz, E. , Tomsia, A. P. , and Ritchie, R. O. , 2010, “ A Novel Biomimetic Approach to the Design of High-Performance Ceramic–Metal Composites,” J. R. Soc. Interface, 7(46), pp. 741–753. [CrossRef] [PubMed]
Gu, G. X. , Su, I. , Sharma, S. , Voros, J. L. , Qin, Z. , and Buehler, M. J. , 2016, “ Three-Dimensional-Printing of Bio-Inspired Composites,” ASME J. Biomech. Eng., 138(2), p. 021006. [CrossRef]


Grahic Jump Location
Fig. 4

Larger grid size solutions for minimizing stiffness (maximizing compliance) and maximizing stiffness. Convergence of solutions: (a) For a 16 × 16 system, the optimal solution obtained from the algorithm remains the same. The graph shows that the initial effective stiffness starts out high and decreases as iteration increases. (b) The optimal solution remains the same for maximizing stiffness and the graph shows an increase in effective stiffness with iteration. These two graphs show the solutions converging to the theoretical optimal values of compliance and stiffness.

Grahic Jump Location
Fig. 3

Optimized solutions for maximizing effective compliance and effective stiffness for 8 × 8 and 12 × 12 grid size systems. (a) The maximum compliance solution is a geometry in which the soft (black) and stiff (gray) materials are in series with the loading conditions. From different initial random geometries, the algorithm leads to the optimized design of in-series materials. (b) The maximum stiffness solution is a geometry in which the soft and stiff materials are in parallel. Similarly, the algorithm leads to the optimized design starting from random initial geometries. These final geometries show that there is no single optimum solution for these problems, but many optimal solutions solving the same objective function.

Grahic Jump Location
Fig. 2

Algorithm organization. The algorithm takes in as an input a random initial population of soft and stiff elements. At every iteration, an objective function value of the system is calculated for every element when only it is switched. Elements switch from soft to stiff and vice versa. The switch that generates the highest objective value increase is kept for the next iteration. This process is repeated until there is no switch that generates a higher objective value compared to the previous iteration at which point the algorithm exits outputting the final geometry.

Grahic Jump Location
Fig. 1

Model formulation. Targeted material property optimization is performed by algorithmic assignment of stiff and soft elements in a multiphase building block. The prescribed binary distribution of element stiffness defines stiffness and volume ratios for each initial geometry. The material contains an edge crack and undergoes tensile loading under mode I failure with displacement controlled boundary conditions (“dx”). a is the length of the sample (square), and b is the length of the crack.

Grahic Jump Location
Fig. 5

Solutions from brute force method to validate algorithm when optimizing for toughness. The chart describes the differences between systems A and B. System A has a grid size of 8 × 8 and an edge crack that is 50% of the system length, with a being the length of the sample. The algorithm begins with a random geometry and generates the solution that exactly matches the solution obtained from brute force (i.e., checking every possible solution and selecting the best). The algorithm, however, is orders of magnitude faster than the brute force method. System B has a grid size of 10 × 10 and an edge crack that is 20% of the system length. The algorithm also obtains the same solution as the brute force method. Tbrute is the toughness obtained from the brute force method and Talg is the toughness obtained from our algorithm, and the ratio shows unity. This performance confirms that our algorithm is effective and robust.

Grahic Jump Location
Fig. 6

Various case studies. (a) The geometry for the objective function that optimizes toughness modulus only. (b) The geometry for the objective function that optimizes toughness along with stiffness. More elements need to be in parallel to maintain high stiffness and also toughness. (c) The geometry for the objective function that optimizes toughness along with compliance. More elements are spread out to be in series with each other. Variables: f is the objective function, T is the toughness modulus, T0 is the initial geometry toughness modulus, Eeff is the effective stiffness of the system, and E0 is the effective stiffness of the initial population.

Grahic Jump Location
Fig. 7

Strain field for homogenous design compared to optimized design for 40 × 40 grid size. White space represents the stiff material and black space represents the soft material. Delocalization of strain seen in optimized design from crack tip area to soft material locations.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In