IP Library Granted Patent US 12,229,476
Granted Patent B2
US 12,229,476 · App. 17/236,917 · Granted Feb 18, 2025

Computer aided generative design with modal analysis driven shape modification process

Inventors: Ravi Kumar Burla (Novi, MI); Jesus Rodriguez (Farmington, MI); Siavash Navadeh Meshkat (Carlsbad, CA); Jaesung Eom (Wexford, PA)
Assignee: Autodesk, Inc.
G06F30/12G06F17/16G06F30/10G06F30/20G06N3/0475
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,229,476
App. No.
17/236,917
Granted
Feb 18, 2025
Kind
B2
Abstract

Methods, systems, and apparatus, including medium-encoded computer program products, for computer aided design of physical structures using generative design processes include: obtaining one or more design criteria, a model of an object for which a physical structure is to be manufactured, and a set of eigenmodes from a modal analysis of the model of the object; extracting a proper subset of non-zero eigenmodes from the set of eigenmodes, wherein the proper subset of non-zero eigenmodes include at least three lowest valued, non-zero eigenmodes; combining data of the proper subset of non-zero eigenmodes to form a strain energy field for the model of the object; iteratively modifying a generatively designed shape of the model of the object using the strain energy field to drive changes to the generatively designed shape of the model; and providing the generatively designed shape of the model of the object.

Claims (72)

1. A method comprising:

obtaining one or more design criteria and a three dimensional model of an object for which a corresponding physical structure is to be manufactured, wherein the three dimensional model is a computer-implemented model having a design space in which the physical structure is to be designed;

selecting non-rigid body eigenmodes of a finite element stiffness matrix from a modal analysis of the three dimensional model of the object, wherein the non-rigid body eigenmodes that are selected are less than ten in number;

superposing strain energy density fields extracted from the non-rigid body eigenmodes to form a synthetic strain energy field for the three dimensional model of the object using weighted strain energies from the strain energy density fields;

iteratively modifying a shape of the three dimensional model of the object using the synthetic strain energy field to drive changes to the shape of the three dimensional model, wherein the iteratively modifying comprises

calculating a shape change velocity field for the shape of the three dimensional model of the object using the synthetic strain energy field, and

updating computer-implemented geometry of the shape within the design space using the shape change velocity field, wherein the geometry of the shape comprises computer representation of individual shape elements; and

providing the shape of the three dimensional model of the object, which is useable in manufacturing the physical structure using one or more computer-controlled manufacturing systems.

2. The method of claim 1 , wherein the selecting comprises selecting a number of the non-rigid body eigenmodes based on computer processing speed and a desired upper limit to processing time for the iteratively modifying.

3. The method of claim 2 , comprising receiving input from a user specifying the desired upper limit to processing time for the iteratively modifying.

4. The method of claim 1 , wherein the selecting comprises:

obtaining one or more frequencies of vibration for which structural performance of the physical structure is to be enhanced; and

extracting the strain energy density fields from the non-rigid body eigenmodes that correspond to the one or more frequencies of vibration.

5. The method of claim 4 , wherein obtaining one or more frequencies of vibration comprises receiving user specification of a specified frequency or a range of frequencies for which structural performance of the physical structure is to be enhanced, and the selecting comprises identifying the non-rigid body eigenmodes that correspond to the one or more frequencies of vibration as those eigenmodes with resonate frequencies that fall within the range of frequencies or that are within a tolerance band around the specified frequency.

6. The method of claim 1 , wherein the iteratively modifying comprises:

checking for completion of a percentage reduction in volume for the shape that is less than a target volume reduction of an original volume of the three dimensional model of the object; and

re-computing the synthetic strain energy field in response to the completion of the percentage reduction in volume for the shape that is less than the target volume reduction of the original volume, the re-computing comprising:

performing an updated modal analysis on a modified version of the shape of the three dimensional model of the object after the completion of the percentage reduction in volume for the shape that is less than the target volume reduction of the original volume;

selecting updated non-rigid body eigenmodes of an updated finite element stiffness matrix from the updated modal analysis; and

superposing updated strain energy density fields extracted from the updated non-rigid body eigenmodes selected from the updated modal analysis to update the synthetic strain energy field used in the iteratively modifying after the completion of the percentage reduction in volume for the shape that is less than the target volume reduction of the original volume.

7. The method of claim 1 , wherein the one or more design criteria comprise a manufacturing constraint, and the iteratively modifying comprises:

checking whether the manufacturing constraint is violated in a given iteration of the iteratively modifying; and

re-computing the synthetic strain energy field in response to the manufacturing constraint being violated in the given iteration, the re-computing comprising:

performing an updated modal analysis on a modified version of the shape of the three dimensional model of the object after the manufacturing constraint is violated in the given iteration;

selecting updated non-rigid body eigenmodes of an updated finite element stiffness matrix from the updated modal analysis; and

superposing updated strain energy density fields extracted from the updated non-rigid body eigenmodes selected from the updated modal analysis to update the synthetic strain energy field used in the iteratively modifying after the given iteration.

8. The method of claim 1 , wherein the calculating, in at least one iteration of the iteratively modifying, comprises:

performing numerical simulation based on a boundary condition specification for the three dimensional model of the object to form a strain energy field for the three dimensional model of the object; and

modifying the strain energy field formed by the numerical simulation using the synthetic strain energy field or modifying the synthetic strain energy field using the strain energy field formed by the numerical simulation.

9. The method of claim 1 , wherein the iteratively modifying is performed without any input regarding boundary conditions for the three dimensional model of the object, and the method comprises providing a user interface that provides user control over weighting of the strain energies from the strain energy density fields to favor some eigenmodes over other eigenmodes.

10. The method of claim 1 , wherein the providing comprises:

generating toolpath specifications for a manufacturing machine using the shape of the three dimensional model of the object; and

manufacturing at least a portion of the physical structure, or a mold for the physical structure, with the manufacturing machine using the toolpath specifications.

11. A system comprising:

a non-transitory storage medium having instructions of a computer aided design program stored thereon; and

one or more data processing apparatus configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to

obtain one or more design criteria and a three dimensional model of an object for which a corresponding physical structure is to be manufactured, wherein the three dimensional model is a computer-implemented model having a design space in which the physical structure is to be designed;

select non-rigid body eigenmodes of a finite element stiffness matrix from a modal analysis of the three dimensional model of the object, wherein the non-rigid body eigenmodes that are selected are less than ten in number;

superpose strain energy density fields extracted from the non-rigid body eigenmodes to form a synthetic strain energy field for the three dimensional model of the object using weighted strain energies from the strain energy density fields;

iteratively modify a shape of the three dimensional model of the object using the synthetic strain energy field to drive changes to the shape of the three dimensional model, wherein the iterative modification comprises

calculating a shape change velocity field for the shape of the three dimensional model of the object using the synthetic strain energy field, and

updating computer-implemented geometry of the shape within the design space using the shape change velocity field, wherein the geometry of the shape comprises computer representation of individual shape elements; and

provide the shape of the three dimensional model of the object, which is useable in manufacturing the physical structure using one or more computer-controlled manufacturing systems.

12. The system of claim 11 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to select a number of the non-rigid body eigenmodes based on computer processing speed and a desired upper limit to processing time for the iterative modification.

13. The system of claim 12 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to receive input from a user specifying the desired upper limit to processing time for the iterative modification.

14. The system of claim 11 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to:

obtain one or more frequencies of vibration for which structural performance of the physical structure is to be enhanced; and

extract the strain energy density fields from the non-rigid body eigenmodes that correspond to the one or more frequencies of vibration.

15. The system of claim 14 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to iteratively modify the shape of the three dimensional model of the object by being configured to:

receive user specification of a specified frequency or a range of frequencies for which structural performance of the physical structure is to be enhanced; and

identify the non-rigid body eigenmodes that correspond to the one or more frequencies of vibration as those eigenmodes with resonate frequencies that fall within the range of frequencies or that are within a tolerance band around the specified frequency.

16. The system of claim 11 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to iteratively modify the shape of the three dimensional model of the object by being configured to:

check for completion of a percentage reduction in volume for the shape that is less than a target volume reduction of an original volume of the three dimensional model of the object; and

re-compute the synthetic strain energy field in response to the completion of the percentage reduction in volume for the shape that is less than the target volume reduction of the original volume by being configured to:

perform an updated modal analysis on a modified version of the shape of the three dimensional model of the object after the completion of the percentage reduction in volume for the shape that is less than the target volume reduction of the original volume;

select updated non-rigid body eigenmodes of an updated finite element stiffness matrix from the updated modal analysis; and

superpose updated strain energy density fields extracted from the updated non-rigid body eigenmodes selected from the updated modal analysis to update the synthetic strain energy field used in the iterative modification after the completion of the percentage reduction in volume for the shape that is less than the target volume reduction of the original volume.

17. The system of claim 11 , wherein the one or more design criteria comprise a manufacturing constraint, and the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to:

check whether the manufacturing constraint is violated in a given iteration of the iterative modification; and

re-compute the synthetic strain energy field in response to the manufacturing constraint being violated in the given iteration by performing operations comprising:

performing an updated modal analysis on a modified version of the shape of the three dimensional model of the object after the manufacturing constraint is violated in the given iteration;

selecting updated non-rigid body eigenmodes of an updated finite element stiffness matrix from the updated modal analysis; and

superposing updated strain energy density fields extracted from the updated non-rigid body eigenmodes selected from the updated modal analysis to update the synthetic strain energy field used in the iterative modification after the given iteration.

18. The system of claim 11 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to:

perform numerical simulation based on a boundary condition specification for the three dimensional model of the object to form a strain energy field for the three dimensional model of the object; and

modify the strain energy field formed by the numerical simulation using the synthetic strain energy field or modify the synthetic strain energy field using the strain energy field formed by the numerical simulation.

19. The system of claim 11 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to:

iteratively modify the shape of the three dimensional model of the object without any input regarding boundary conditions for the three dimensional model of the object; and

provide a user interface that provides user control over weighting of the strain energies from the strain energy density fields to favor some eigenmodes over other eigenmodes.

20. The system of claim 11 , wherein the one or more data processing apparatus are configured to run the instructions of the computer aided design program to cause the one or more data processing apparatus to:

generate toolpath specifications for a manufacturing machine using the shape of the three dimensional model of the object; and

manufacture at least a portion of the physical structure, or a mold for the physical structure, with the manufacturing machine using the toolpath specifications.

Assignments (2)
CHANGE OF ADDRESS FOR ASSIGNEE Recorded Aug 19, 2022
From: AUTODESK, INC.
To: AUTODESK, INC.
Reel/Frame 061572/0061 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Apr 29, 2021
From: BURLA, RAVI KUMAR; RODRIGUEZ, JESUS; MESHKAT, SIAVASH NAVADEH; EOM, JAESUNG
To: AUTODESK, INC.
Reel/Frame 056089/0791 →
Continuity (1)
Related Publication 20220358255A1 · Nov 10, 2022
References Cited (68)
US 5841892A · McGrath · 1998 [cited by examiner]
US 7552032B2 · Krishnapillai · 2009 [cited by examiner]
US 9978177B2 · Mehr · 2018 [cited by examiner]
US 11321508B2 · Harris · 2022 [cited by examiner]
US 11663379B2 · Harris · 2023 [cited by examiner]
US 11675333B2 · Bandara · 2023 [cited by examiner]
US 20040194051A1 · Croft · 2004 [cited by examiner]
US 20060139347A1 · Choi · 2006 [cited by examiner]
US 20060259283A1 · Brughmans · 2006 [cited by examiner]
US 20120179430A1 · Aish · 2012 [cited by examiner]
US 20140214370A1 · Olhofer · 2014 [cited by examiner]
US 20160009368A1 · Kismarton · 2016 [cited by examiner]
US 20160253444A1 · Rodrigues · 2016 [cited by examiner]
US 20160321384A1 · Pal · 2016 [cited by examiner]
US 20170193699A1 · Mehr · 2017 [cited by examiner]
US 20180345647A1 · Morris et al. · 2018 [cited by applicant]
US 20180349531A1 · Morris et al. · 2018 [cited by applicant]
US 20190236220A1 · Nomura · 2019 [cited by examiner]
US 20190272356A1 · Trivedi · 2019 [cited by examiner]
US 20200065429A1 · Nakhjavani · 2020 [cited by examiner]
US 20200143009A1 · Schmidt · 2020 [cited by examiner]
US 20200150623A1 · Bandara et al. · 2020 [cited by applicant]
US 20200151286A1 · Willis et al. · 2020 [cited by applicant]
US 20200156323A1 · Woytowitz · 2020 [cited by examiner]
US 20200191766A1 · Kim · 2020 [cited by examiner]
US 20200356638A1 · Nomura · 2020 [cited by examiner]
US 20210182456A1 · Couret · 2021 [cited by examiner]
US 20210279384A1 · Bahamonde Jacome · 2021 [cited by examiner]
US 20210350036A1 · Burla · 2021 [cited by examiner]
US 20210406412A1 · Cramer · 2021 [cited by examiner]
US 20220250342A1 · Metzner · 2022 [cited by examiner]
US 20220405442A1 · Schmidt · 2022 [cited by examiner]
US 20240012956A1 · Rauh · 2024 [cited by examiner]
Allaire, Grégoire, and Georgios Michailidis. “Modal basis approaches in shape and topology optimization of frequency response problems.” International Journal for Numerical Methods in Engineering 113.8 (2018): 1258-1299… [cited by examiner]
Dewulf, W.; De Roeck, G.. Modal strain energy method critically revised Proceedings of the SPIE—The International Society for Optical Engineering 3989: 99-110. SPIE—Int. Soc. Opt. Eng. (2000) (Year: 2000). [cited by examiner]
Sigmund, Ole, and Joakim Petersson. “Numerical instabilities in topology optimization: a survey on procedures dealing with checkerboards, mesh-dependencies and local minima.” Structural optimization 16 (1998): 68-75. (Y… [cited by examiner]
Dienemann, R., A. Schumacher, and S. Fiebig. “Topology optimization for finding shell structures manufactured by deep drawing.” Structural and Multidisciplinary Optimization 56 (2017): 473-485. (Year: 2017). [cited by examiner]
Völkl, Harald, et al. “An efficient bionic topology optimization method for transversely isotropic materials.” Composite Structures 204 (2018): 359-367. (Year: 2018). [cited by examiner]
Xia, L., Xia, Q., Huang, X et al. Bi-directional Evolutionary Structural Optimization on Advanced Structures and Materials: A Comprehensive Review. Arch Computat Methods Eng 25, 437-478, Published Nov. 19, 2016 (Year: 2… [cited by examiner]
Allaire et al., “Modal basis approaches in shape and topology optimization of frequency response problems,” International Journal for Numerical Methods in Engineering, Feb. 2017, 113(8):1258-1299. [cited by applicant]
Allaire et al., “Structural optimization under overhang constraints imposed by additive manufacturing technologies,” Journal of Computational Physics, Sep. 2017, 351:295-328. [cited by applicant]
Dalklint et al., “Eigenfrequency constrained topology optimization of finite strain hyperelastic structures,” Struct Multidisc Optim, May 2020, 61(6):2577-2594. [cited by applicant]
Extended European Search Report in European Appln. No. 22167040.9, dated Sep. 2, 2022, 16 pages. [cited by applicant]
Liu et al., “Eigenvalue topology optimization of structures using a parameterized level set method,” Struct Multidisc Optim, Apr. 2014, 50(4):573-591. [cited by applicant]
Yamada et al., “A topology optimization method based on the level set method incorporating a fictitious interface energy,” Comput. Methods Appl. Mech. Engrg., Nov. 2010, 199(45-48):2876-2891. [cited by applicant]
Arisoy et al., “Design and Topology Optimization of Lattice Structures Using Deformable Implicit Surfaces for Additive Manufacturing,” in ASME 2015 International Design Engineering Technical Conferences and Computers an… [cited by applicant]
Autodesk.com [online], “Section 21: Modal or Natural Frequency Analysis,” published on Nov. 8, 2018, retrieved Apr. 12, 2021, retrieved from URL<https://knowledge.autodesk.com/support/inventor-nastran>, 11 pages. [cited by applicant]
Autodesk.com [online], “Section 22: Modal Analysis Exercise,” published on Nov. 8, 2018, retrieved on Apr. 12, 2021, retrieved from URL<https://knowledge.autodesk.com/>, 21 pages. [cited by applicant]
Brackett et al., “Topology Optimization for Additive Manufacturing,” in Proceedings of the Solid Freeform Fabrication Symposium (2011) 15 pages. [cited by applicant]
Dapogny et al., “Geometric Constraints for Shape and Topology Optimization in Architectural Design,” Computational Mechanics, Springer-Verlag, 2017, 59(6):933-965.10.1007/s00466-017-1383-6. hal-01354004v3. [cited by applicant]
Gibou et al., “A Review of Level-set Methods and Some Recent Applications,” Journal of Computational Physics (2018) 353:82-109. [cited by applicant]
Guest and Zhu, “Casting and Milling Restrictions in Topology Optimization via Projection-based Algorithms,” in ASME 2012 International Design Engineering Technical Conferences and Computers and Information in Engineerin… [cited by applicant]
Ikeya and Shimoda, “Multi-objective Free-form Optimization for Shape and Thickness of Shell Structures with Composite Materials,” in 11th World Congress on Structural and Multidisciplinary Optimisation (2015) 6 pages. [cited by applicant]
Joshi et al., “CAD-integrated Topology Optimization (BGCE Honours Project),” Department of Informatics, Technical University of Munich (2016) 77 pages. [cited by applicant]
Langelaar, “Topology Optimization for Multi-axis Machining,” Comput. Methods Appl. Mech. Engrg. (2019) 351:226-252. [cited by applicant]
Liu and Ma, “A Survey of Manufacturing Oriented Topology Optimization Methods,” Advances in Engineering Software (2016) 100:161-175. [cited by applicant]
Liu et al., “Current and Future Trends in Topology Optimization for Additive Manufacturing,” Structural and Multidisciplinary Optimization (2018) 57:2457-2483. [cited by applicant]
Lu et al., “Build-to-Last: Strength to Weight 3D Printed Objects,” ACM Trans. Graph. (2014) 33(4):97:1-97:10. [cited by applicant]
Nakayama and Shimoda, “Shape-topology Optimization for Designing Shell Structures,” in VII European Congress on Computational Methods in Applied Sciences and Engineering (2016) 10 pages. [cited by applicant]
PTC.com [online], “Model Analysis,” published on or before Apr. 12, 2021, retrieved on Feb. 12, 2021, retrieved from URL<https://www.ptc.com/en/technologies/cad/simulation-and-analysis/model-analysis>, 3 pages. [cited by applicant]
Sigmund and Maute, “Topology Optimization Approaches—a Comparative Review,” Struct. Multidisc. Optim. (2013) 48:1031-1055. [cited by applicant]
Unknown author, “Topology Optimization R18.0 Feature and Usage Highlights,” © 2016 ANSYS, Inc., Mar. 12, 2017, 29 pages. [cited by applicant]
Van Dijk et al., “Level-set Methods for Structural Topology Optimization: a Review,” Struct. Multidisc. Optim. (2013) 48:437-472. [cited by applicant]
Vatanabe et al., “Topology Optimization with Manufacturing Constraints: a Unified Projection-based Approach,” Advances in Engineering Software (2016) 100:97-112. [cited by applicant]
Visweswaran, “Modal Analysis—What is it and How is it Useful?,” Graphler Technology Solutions, Dec. 25, 2020, 4 pages. [cited by applicant]
Xia et al., “A Level Set Based Method for the Optimization of Cast Part,” Struct. Multidisc. Optim. (2010) 41:735-747. [cited by applicant]
Liu et al., “Concurrent deposition path planning and structural topology optimization for additive manufacturing,” Rapid Prototyping Journal, Aug. 2017, 23(5):930-942. [cited by applicant]
Office Action in European Appln. No. 22167040.9, dated Dec. 16, 2024, 14 pages. [cited by applicant]