On tailoring fracture resistance of brittle structures

A level set interface-enriched topology optimization approach

Journal Article (2022)
Author(s)

J. Zhang (TU Delft - Computational Design and Mechanics)

F. Van Keulen (TU Delft - Computational Design and Mechanics)

Alejandro M. Aragón (TU Delft - Computational Design and Mechanics)

Research Group
Computational Design and Mechanics
Copyright
© 2022 J. Zhang, A. van Keulen, A.M. Aragon
DOI related publication
https://doi.org/10.1016/j.cma.2021.114189
More Info
expand_more
Publication Year
2022
Language
English
Copyright
© 2022 J. Zhang, A. van Keulen, A.M. Aragon
Research Group
Computational Design and Mechanics
Volume number
388
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

We propose a fully immersed topology optimization procedure to design structures with tailored fracture resistance under linear elastic fracture mechanics assumptions for brittle materials. We use a level set function discretized by radial basis functions to represent the topology and the Interface-enriched Generalized Finite Element Method (IGFEM) to obtain an accurate structural response. The technique assumes that cracks can nucleate at right angles from the boundary, at the location of enriched nodes that are added to enhance the finite element approximation. Instead of performing multiple finite element analyses to evaluate the energy release rates (ERRs) of all potential cracks—a procedure that would be computationally intractable—we approximate them by means of topological derivatives after a single enriched finite element analysis of the uncracked domain. ERRs are then aggregated to construct the objective function, and the corresponding sensitivity formulation is derived analytically by means of an adjoint formulation. Several numerical examples demonstrate the technique's ability to tailor fracture resistance, including the well-known benchmark L-shaped bracket and a multiple-loading optimization problem for obtaining a structure with fracture resistance anisotropy.