{"id":{"repo_id":"uiuc","oai_identifier":"oai:www.ideals.illinois.edu:2142/95442"},"canonical_url":"https://search.dev.ndltd.org/etd/uiuc/oai:www.ideals.illinois.edu:2142/95442","repository":{"repo_id":"uiuc","name":"University of Illinois - Urbana-Champaign","base_url":"https://www.ideals.illinois.edu/oai-pmh"},"display":{"title":"Multiscale design of nonlinear materials using a shape optimization scheme based on an interface-enriched GFEM","abstract":"Motivated by key advances in manufacturing techniques, the tailoring of materials with specific macroscopic properties has been the focus of active research in mechanical engineering and materials science over the past decade. The key challenge in this line of work is how to optimize the material microstructure to achieve a desired macroscopic constitutive response. The overwhelming majority of this type of inverse design work relies on topology optimization based, primarily, on linear theory. In this work, we develop and implement a method to design particulate composites at the mesoscale using a shape optimization scheme to minimize or maximize a nonlinear cost function at the macroscale while satisfying a set of constraints associated, for example, with the volume fraction of inclusions or with the manufacturing technique. The optimization method relies on three key ‘modules’: multiscale modeling, sensitivity analysis, and optimization. The multiscale modeling is based on a nonlinear finite element solver, which combines a classical homogenization scheme with a NURBS-based Interface-enriched Generalized Finite Element Method (NIGFEM) used to capture accurately and efficiently the displacement field in a heterogeneous material with a finite element discretization that does not conform to the material interfaces. Damage evolution is captured using a three-parameter isotropic damage model able to simulate a wide range of failure responses. The proposed gradient-based shape optimization scheme relies on the stationary nature of the non-conforming meshes used to discretize the periodic unit cell, thereby avoiding mesh distortion issues that plague conventional finite-element-based shape optimization studies. In the current approach, the finite element approximation space used in the NIGFEM is augmented with NURBS to allow for the accurate capture of the weak discontinuity present along complex, curvilinear material interfaces. NURBS are also used to parameterize the design geometry precisely and compactly by a small number of design variables. To compute the derivatives of the cost and constraint functions with respect to the design variables, we also formulate an analytic nonlinear sensitivity, which is simplified by the fact that only the enrichment control points on material interfaces move, appear or disappear during the shape optimization process. The derivations uncover subtle but important new terms involved in the sensitivity of shape functions and their spatial derivatives. Our analytic nonlinear shape sensitivity avoids the technical difficulties encountered in the finite difference or semi-analytical schemes when the boundary intersects an element very close to a node in a non-conforming mesh. In these situations, the boundary may move to another element during the design perturbation step, resulting in changes of the mesh topology, making the differentiation of the stiffness matrix and load vector problematic. We apply the NIGFEM shape optimization scheme to several 2D and 3D structural problems including some benchmark and application examples to demonstrate the performance and accuracy of the method. Based on the multiscale approach, we also design the microstructure of a periodic particulate composite to optimize the volume fraction and distribution of the inclusions for a desired macroscopic nonlinear stress-strain curve.","abstract_html":"Motivated by key advances in manufacturing techniques, the tailoring of materials with specific macroscopic properties has been the focus of active research in mechanical engineering and materials science over the past decade. The key challenge in this line of work is how to optimize the material microstructure to achieve a desired macroscopic constitutive response. The overwhelming majority of this type of inverse design work relies on topology optimization based, primarily, on linear theory. In this work, we develop and implement a method to design particulate composites at the mesoscale using a shape optimization scheme to minimize or maximize a nonlinear cost function at the macroscale while satisfying a set of constraints associated, for example, with the volume fraction of inclusions or with the manufacturing technique. The optimization method relies on three key ‘modules’: multiscale modeling, sensitivity analysis, and optimization. The multiscale modeling is based on a nonlinear finite element solver, which combines a classical homogenization scheme with a NURBS-based Interface-enriched Generalized Finite Element Method (NIGFEM) used to capture accurately and efficiently the displacement field in a heterogeneous material with a finite element discretization that does not conform to the material interfaces. Damage evolution is captured using a three-parameter isotropic damage model able to simulate a wide range of failure responses. The proposed gradient-based shape optimization scheme relies on the stationary nature of the non-conforming meshes used to discretize the periodic unit cell, thereby avoiding mesh distortion issues that plague conventional finite-element-based shape optimization studies. In the current approach, the finite element approximation space used in the NIGFEM is augmented with NURBS to allow for the accurate capture of the weak discontinuity present along complex, curvilinear material interfaces. NURBS are also used to parameterize the design geometry precisely and compactly by a small number of design variables. To compute the derivatives of the cost and constraint functions with respect to the design variables, we also formulate an analytic nonlinear sensitivity, which is simplified by the fact that only the enrichment control points on material interfaces move, appear or disappear during the shape optimization process. The derivations uncover subtle but important new terms involved in the sensitivity of shape functions and their spatial derivatives. Our analytic nonlinear shape sensitivity avoids the technical difficulties encountered in the finite difference or semi-analytical schemes when the boundary intersects an element very close to a node in a non-conforming mesh. In these situations, the boundary may move to another element during the design perturbation step, resulting in changes of the mesh topology, making the differentiation of the stiffness matrix and load vector problematic. We apply the NIGFEM shape optimization scheme to several 2D and 3D structural problems including some benchmark and application examples to demonstrate the performance and accuracy of the method. Based on the multiscale approach, we also design the microstructure of a periodic particulate composite to optimize the volume fraction and distribution of the inclusions for a desired macroscopic nonlinear stress-strain curve.","abstract_has_math":false,"creators":["Raeisi Najafi, Ahmad"],"institution":"University of Illinois at Urbana-Champaign","degree_name":"Ph.D.","degree_level":"Dissertation","degree_discipline":"Theoretical & Applied Mechans","degree_department":null,"school":null,"contributors":["Geubelle, Philippe H.","Tortorelli, Daniel A.","Duarte, Armando C.","James, Kai"],"advisors":[],"committee_chairs":[],"committee_members":[],"year":2017,"date_issued":"2017-03-01T16:36:34Z","date_published":"2017-03-01T16:36:34Z","updated_at":"2026-07-22T22:26:37Z","subjects":["Gradient-based shape optimization scheme","Generalized Finite Element Method (GFEM)"],"languages":["en"],"rights":["Copyright 2016 Ahmad Raeisi Najafi"],"rights_urls":[],"identifier_entries":[]},"links":{"outbound_url":"http://hdl.handle.net/2142/95442","outbound_label":"Handle","outbound_source":"dc:identifier"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor","label":"Contributor","values":["Geubelle, Philippe H.","Tortorelli, Daniel A.","Duarte, Armando C.","James, Kai"]},{"key":"dc:creator","label":"Author","values":["Raeisi Najafi, Ahmad"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date","label":"Dc Date","values":["2017-03-01T16:36:34Z","2019-03-02T10:15:30Z","2016-07-28","2016-12"]},{"key":"dc:type","label":"Dc Type","values":["text"]},{"key":"thesis:degree_discipline","label":"Discipline","values":["Theoretical & Applied Mechans"]},{"key":"thesis:degree_level","label":"Degree Level","values":["Dissertation"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Ph.D."]},{"key":"thesis:institution_name","label":"Thesis Institution Name","values":["University of Illinois at Urbana-Champaign"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Gradient-based shape optimization scheme","Generalized Finite Element Method (GFEM)"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language","label":"Dc Language","values":["en"]},{"key":"dc:rights","label":"Dc Rights","values":["Copyright 2016 Ahmad Raeisi Najafi"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier","label":"Identifier","values":["http://hdl.handle.net/2142/95442"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["Motivated by key advances in manufacturing techniques, the tailoring of materials with specific macroscopic properties has been the focus of active research in mechanical engineering and materials science over the past decade. The key challenge in this line of work is how to optimize the material microstructure to achieve a desired macroscopic constitutive response. The overwhelming majority of this type of inverse design work relies on topology optimization based, primarily, on linear theory. In this work, we develop and implement a method to design particulate composites at the mesoscale using a shape optimization scheme to minimize or maximize a nonlinear cost function at the macroscale while satisfying a set of constraints associated, for example, with the volume fraction of inclusions or with the manufacturing technique. The optimization method relies on three key ‘modules’: multiscale modeling, sensitivity analysis, and optimization. The multiscale modeling is based on a nonlinear finite element solver, which combines a classical homogenization scheme with a NURBS-based Interface-enriched Generalized Finite Element Method (NIGFEM) used to capture accurately and efficiently the displacement field in a heterogeneous material with a finite element discretization that does not conform to the material interfaces. Damage evolution is captured using a three-parameter isotropic damage model able to simulate a wide range of failure responses. The proposed gradient-based shape optimization scheme relies on the stationary nature of the non-conforming meshes used to discretize the periodic unit cell, thereby avoiding mesh distortion issues that plague conventional finite-element-based shape optimization studies. In the current approach, the finite element approximation space used in the NIGFEM is augmented with NURBS to allow for the accurate capture of the weak discontinuity present along complex, curvilinear material interfaces. NURBS are also used to parameterize the design geometry precisely and compactly by a small number of design variables. To compute the derivatives of the cost and constraint functions with respect to the design variables, we also formulate an analytic nonlinear sensitivity, which is simplified by the fact that only the enrichment control points on material interfaces move, appear or disappear during the shape optimization process. The derivations uncover subtle but important new terms involved in the sensitivity of shape functions and their spatial derivatives. Our analytic nonlinear shape sensitivity avoids the technical difficulties encountered in the finite difference or semi-analytical schemes when the boundary intersects an element very close to a node in a non-conforming mesh. In these situations, the boundary may move to another element during the design perturbation step, resulting in changes of the mesh topology, making the differentiation of the stiffness matrix and load vector problematic. We apply the NIGFEM shape optimization scheme to several 2D and 3D structural problems including some benchmark and application examples to demonstrate the performance and accuracy of the method. Based on the multiscale approach, we also design the microstructure of a periodic particulate composite to optimize the volume fraction and distribution of the inclusions for a desired macroscopic nonlinear stress-strain curve.","Submission published under a 24 month embargo labeled 'U of I Access', the embargo will last until 2018-12-01","The student, Ahmad Raeisi Najafi, accepted the attached license on 2016-07-26 at 15:58.","The student, Ahmad Raeisi Najafi, submitted this Dissertation for approval on 2016-07-26 at 16:00.","This Dissertation was approved for publication on 2016-07-28 at 08:47.","DSpace SAF Submission Ingestion Package generated from Vireo submission #10089 on 2017-02-28 at 14:35:36","Made available in DSpace on 2017-03-01T16:36:34Z (GMT). No. of bitstreams: 2 RAEISINAJAFI-DISSERTATION-2016.pdf: 77314353 bytes, checksum: 63df37dc95ee209b26426508e1c32aae (MD5) LICENSE.txt: 4216 bytes, checksum: c1c45c7b4457a2a2a003365fe0e3e87b (MD5) Previous issue date: 2016-07-28","Embargo set by: Seth Robbins for item 98558 Lift date: 2019-03-01T16:37:19Z Reason: Author requested U of Illinois access only (OA after 2yrs) in Vireo ETD system","U of I Only Restriction Lifted for Item 98558 on 2019-03-02T10:15:30Z."]},{"key":"dc:format","label":"Dc Format","values":["application/pdf"]},{"key":"dc:title","label":"Title","values":["Multiscale design of nonlinear materials using a shape optimization scheme based on an interface-enriched GFEM"]}]}],"canonical_facts":{"dc:contributor":["Geubelle, Philippe H.","Tortorelli, Daniel A.","Duarte, Armando C.","James, Kai"],"dc:creator":["Raeisi Najafi, Ahmad"],"dc:date":["2017-03-01T16:36:34Z","2019-03-02T10:15:30Z","2016-07-28","2016-12"],"dc:description":["Motivated by key advances in manufacturing techniques, the tailoring of materials with specific macroscopic properties has been the focus of active research in mechanical engineering and materials science over the past decade. The key challenge in this line of work is how to optimize the material microstructure to achieve a desired macroscopic constitutive response. The overwhelming majority of this type of inverse design work relies on topology optimization based, primarily, on linear theory. In this work, we develop and implement a method to design particulate composites at the mesoscale using a shape optimization scheme to minimize or maximize a nonlinear cost function at the macroscale while satisfying a set of constraints associated, for example, with the volume fraction of inclusions or with the manufacturing technique. The optimization method relies on three key ‘modules’: multiscale modeling, sensitivity analysis, and optimization. The multiscale modeling is based on a nonlinear finite element solver, which combines a classical homogenization scheme with a NURBS-based Interface-enriched Generalized Finite Element Method (NIGFEM) used to capture accurately and efficiently the displacement field in a heterogeneous material with a finite element discretization that does not conform to the material interfaces. Damage evolution is captured using a three-parameter isotropic damage model able to simulate a wide range of failure responses. The proposed gradient-based shape optimization scheme relies on the stationary nature of the non-conforming meshes used to discretize the periodic unit cell, thereby avoiding mesh distortion issues that plague conventional finite-element-based shape optimization studies. In the current approach, the finite element approximation space used in the NIGFEM is augmented with NURBS to allow for the accurate capture of the weak discontinuity present along complex, curvilinear material interfaces. NURBS are also used to parameterize the design geometry precisely and compactly by a small number of design variables. To compute the derivatives of the cost and constraint functions with respect to the design variables, we also formulate an analytic nonlinear sensitivity, which is simplified by the fact that only the enrichment control points on material interfaces move, appear or disappear during the shape optimization process. The derivations uncover subtle but important new terms involved in the sensitivity of shape functions and their spatial derivatives. Our analytic nonlinear shape sensitivity avoids the technical difficulties encountered in the finite difference or semi-analytical schemes when the boundary intersects an element very close to a node in a non-conforming mesh. In these situations, the boundary may move to another element during the design perturbation step, resulting in changes of the mesh topology, making the differentiation of the stiffness matrix and load vector problematic. We apply the NIGFEM shape optimization scheme to several 2D and 3D structural problems including some benchmark and application examples to demonstrate the performance and accuracy of the method. Based on the multiscale approach, we also design the microstructure of a periodic particulate composite to optimize the volume fraction and distribution of the inclusions for a desired macroscopic nonlinear stress-strain curve.","Submission published under a 24 month embargo labeled 'U of I Access', the embargo will last until 2018-12-01","The student, Ahmad Raeisi Najafi, accepted the attached license on 2016-07-26 at 15:58.","The student, Ahmad Raeisi Najafi, submitted this Dissertation for approval on 2016-07-26 at 16:00.","This Dissertation was approved for publication on 2016-07-28 at 08:47.","DSpace SAF Submission Ingestion Package generated from Vireo submission #10089 on 2017-02-28 at 14:35:36","Made available in DSpace on 2017-03-01T16:36:34Z (GMT). No. of bitstreams: 2 RAEISINAJAFI-DISSERTATION-2016.pdf: 77314353 bytes, checksum: 63df37dc95ee209b26426508e1c32aae (MD5) LICENSE.txt: 4216 bytes, checksum: c1c45c7b4457a2a2a003365fe0e3e87b (MD5) Previous issue date: 2016-07-28","Embargo set by: Seth Robbins for item 98558 Lift date: 2019-03-01T16:37:19Z Reason: Author requested U of Illinois access only (OA after 2yrs) in Vireo ETD system","U of I Only Restriction Lifted for Item 98558 on 2019-03-02T10:15:30Z."],"dc:format":["application/pdf"],"dc:identifier":["http://hdl.handle.net/2142/95442"],"dc:language":["en"],"dc:rights":["Copyright 2016 Ahmad Raeisi Najafi"],"dc:subject":["Gradient-based shape optimization scheme","Generalized Finite Element Method (GFEM)"],"dc:title":["Multiscale design of nonlinear materials using a shape optimization scheme based on an interface-enriched GFEM"],"dc:type":["text"],"thesis:degree_discipline":["Theoretical & Applied Mechans"],"thesis:degree_level":["Dissertation"],"thesis:degree_name":["Ph.D."],"thesis:institution_name":["University of Illinois at Urbana-Champaign"]},"updated_at":"2026-07-22T22:26:37Z"}