US20180319087A1 - Estimating physical property of 3d printed parts - Google Patents
Estimating physical property of 3d printed parts Download PDFInfo
- Publication number
- US20180319087A1 US20180319087A1 US15/970,291 US201815970291A US2018319087A1 US 20180319087 A1 US20180319087 A1 US 20180319087A1 US 201815970291 A US201815970291 A US 201815970291A US 2018319087 A1 US2018319087 A1 US 2018319087A1
- Authority
- US
- United States
- Prior art keywords
- unit cell
- properties
- build
- extrusion printer
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B29—WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
- B29C—SHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
- B29C64/00—Additive manufacturing, i.e. manufacturing of three-dimensional [3D] objects by additive deposition, additive agglomeration or additive layering, e.g. by 3D printing, stereolithography or selective laser sintering
- B29C64/30—Auxiliary operations or equipment
- B29C64/386—Data acquisition or data processing for additive manufacturing
- B29C64/393—Data acquisition or data processing for additive manufacturing for controlling or regulating additive manufacturing processes
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B29—WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
- B29C—SHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
- B29C64/00—Additive manufacturing, i.e. manufacturing of three-dimensional [3D] objects by additive deposition, additive agglomeration or additive layering, e.g. by 3D printing, stereolithography or selective laser sintering
- B29C64/10—Processes of additive manufacturing
- B29C64/106—Processes of additive manufacturing using only liquids or viscous materials, e.g. depositing a continuous bead of viscous material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y10/00—Processes of additive manufacturing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y50/00—Data acquisition or data processing for additive manufacturing
- B33Y50/02—Data acquisition or data processing for additive manufacturing for controlling or regulating additive manufacturing processes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Materials Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Manufacturing & Machinery (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Optics & Photonics (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
Abstract
Description
- This application claims priority to U.S. Provisional Application No. 62/502,953, filed on May 8, 2017. The disclosure of the prior application is considered part of and is incorporated by reference in the disclosure of this application.
- This specification relates to estimating the physical properties of 3D printed objects.
- Additive manufacturing (AM), also known as solid free form fabrication or three-dimensional (3D) printing, refers to any manufacturing process where 3D objects are built up from raw material (generally powders, liquids, suspensions, or molten solids) in a series of two-dimensional layers or cross-sections. An example of additive manufacturing is Fused Filament Fabrication (FFF). A 3D extrusion printer typically uses FFF to lay down material, such as a plastic filament or metal wire that is unwound from a coil, in layers to produce a 3D printed object.
- Traditionally, seasoned engineers make an educated guess about the physical properties of a 3D printed object produced by FFF, based on the properties of the filaments used for FFF and guided by their experience. If a more accurate estimation is needed, a test object is fabricated using a 3D extrusion printer for experimental determination of the physical properties of interest. Additional test objects are fabricated to determine the impact of different 3D printing parameters on the physical properties of the printed object.
- Finite element analysis (FEA) is a numerical method for simulating, for example, structural mechanics, heat transfer, fluid flow, mass transport, and electromagnetism. It is typically a two-step process. First, the domain of the problem is divided into a collection of subdomains, with each subdomain represented by a set of element equations to the original problem. Second, all sets of element equations are systematically recombined into a global system of equations for the final calculation. Mesh generation techniques are used for dividing a complex problem into a collection of subdomains. The subdomains generated through mesh generation are called mesh elements. FEA can be used to simulate various physical properties of an object contained within the domain.
- Unit cell homogenization is a technique that enables replacement of a complex medium by an ‘equivalent’ homogeneous medium to solve a global problem. A representative volume element (RVE) is often used for performing unit cell homogenization. A RVE is understood to be a sub-volume of an inhomogeneous medium that is of sufficient size for providing geometrical information necessary for obtaining an appropriate homogenized behavior.
- This specification relates to techniques and methods for estimating the physical properties of 3D printed objects.
- In general, one or more aspects of the subject matter described in this specification can be embodied in one or more methods that include: obtaining a three dimensional (3D) model of a 3D object to be manufactured by a 3D extrusion printer; obtaining at least a portion of tool path data for the 3D extrusion printer to build the 3D object in accordance with the 3D model; creating unit cell geometry from the at least a portion of tool path data using at least one limitation on one or more unit cells of the unit cell geometry with respect to an aspect of the 3D extrusion printer's build of the 3D object; generating a numerical simulation of one or more microstructural properties of the 3D object to be built using the unit cell geometry; and estimating one or more macroscale properties, of the 3D object to be built, from the numerical simulation.
- In various implementations, creating the unit cell geometry can be based on various criteria. For example, creating the unit cell geometry can include limiting a minimum volume of each unit cell to ensure the unit cell spans at least three layers of the 3D extrusion printer's build of the 3D object, as determined from a layer thickness aspect of the 3D extrusion printer's build of the 3D object. As another example, creating the unit cell geometry can include checking, for each of the one or more unit cells, a percentage of the unit cell falling inside a containing envelope of the 3D extrusion printer's build of the 3D object; and moving any of the one or more unit cells, which have less than fifty percent falling inside the containing envelope, to have a higher percentage of the unit cell within the containing envelope. As another example, creating the unit cell geometry can include limiting a maximum volume of each unit cell to ensure that locality is maintained. The maximum volume of each unit cell, for example, can be a volume of a mesh element containing the unit cell. These and other criteria can be used in various implementations, either alone or in combination with each other.
- In various implementations, the one or more macroscale properties can include one or more mechanical properties, one or more thermal properties, one or more electromagnetic properties, or a combination thereof.
- Particular embodiments of the subject matter described in this specification can be implemented to realize one or more of the following advantages. Detailed microstructural properties of 3D printed objects can be estimated as early as at the 3D modeling phase. Design cycle of 3D printed objects can be shortened by predicting possible failures of the 3D printed objects during design and modeling phase. Physical properties can be quickly estimated by simulating limited unit cell volumes. The estimation can be robust and applicable to various 3D printing techniques including filling or mixing multiple filament types. The estimated microstructural properties can be extended to macroscale FEA simulation of the entire 3D object. A smooth workflow is provided by automating estimation of microstructural or macroscale physical properties through setting up and running FEA simulations from a selected region in the 3D model of an object. The macroscale properties of a 3D printed object can be determined with reduced computational resource and time. The macroscale properties of a 3D printed object can be automatically checked against a set of requirements for compliance. A toolpath meeting requirements for physical properties when fabricated can be generated using design of experiment or iterative parametric optimization. Multiple physical properties governed by different domains of physics (e.g., thermal and mechanical) can be simultaneously simulated to generate a toolpath, when fabricated, meeting requirements for physical properties of more than one physics domain (e.g., safe strength and favorable thermal conductivity.)
- The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
-
FIG. 1 shows an example of a system to design and estimate the macroscale properties of a 3D printed object. -
FIGS. 2A and 2B show a schematic diagram of a cross-section of a 3D printed object. -
FIG. 3 is a flowchart showing an example of a process to perform macroscale simulation of a 3D model of a 3D printed object using unit cell homogenization. -
FIG. 4 is a schematic diagram of an example mesh element used in FEA. -
FIG. 5 is a flowchart showing an example of a process to simulate the macroscale properties of a 3D printed object. - Like reference symbols in the various drawings indicate like elements.
- Typically, a designer of a 3D object for fabrication by a 3D extrusion printer does not accurately know what the macroscale properties of the manufactured 3D printed part will be during the design phase. For example, the macroscale mechanical properties (e.g. elasticity and strength) of a 3D printed object can be highly dependent, among many others, on the printing orientations and sequences of a 3D extrusion printer's fabrication flow. For example, a 3D printed object produced by FFF can have different printing orientation, paths, and infill density, potentially at different layers or portions of the 3D printed object. The combination of different printing orientations, paths, and infill densities lead to different local ‘microstructures’, with different void distribution and fused filament orientations. These differences in the local microstructures affect the local microstructural properties of a 3D printed object, which in turn affects the object's macroscale properties as a whole. Furthermore, 3D printed objects can have different macroscale properties along different axes, a phenomenon called anisotropy. For example, some mechanical properties of 3D objects printed using FFF can be anisotropic. The mechanical anisotropy of a 3D printed object can be attributed in part to the microstructures and filament-to-filament bond strength. Due to these factors, macroscale properties of a 3D printed object cannot be simply determined from the material properties of the individual filaments that constitute the object.
- In fact, some macroscale properties may still be unknown to the designers even after the 3D printed object has been fabricated, as experimental determination of macroscale properties can be complicated, time consuming, expensive, and destructive. However, the designers often desire to know the macroscale properties such as the mechanical, thermal, or electromagnetic properties as these properties may affect the function or behavior of the fabricated object. Therefore, a method and system for estimating the macroscale properties of a 3D printed object from a 3D model is desired.
-
FIG. 1 shows an example of asystem 100 to design and estimate the macroscale properties of a 3D printed object. Acomputer 110 includes aprocessor 112 and amemory 114, and thecomputer 110 can be connected to anetwork 140, which can be a private network, a public network, a virtual private network, etc. One or moreremote computer systems 150 and a 3D extrusion printer 160 can be connected to thecomputer 110 through thenetwork 140. In some implementations, the 3D extrusion printer 160 is connected directly to thecomputer 110. Theprocessor 112 can be one or more hardware processors, which can each include multiple processor cores. Thememory 114 can include both volatile and non-volatile memory, such as Random Access Memory (RAM) and Flash RAM. Thecomputer 110 can include various types of computer storage media and devices, which can include thememory 114, to store instructions of programs that run on theprocessor 112. - Such programs can include a
3D modeling program 116 a, a toolpath generation program 116 b, and anumerical simulation program 116 c, which can each run locally oncomputer 110, remotely on a computer of one or moreremote computer systems 150, or both. The3D modeling program 116 a presents a user interface (UI) 122 on adisplay device 120 of thecomputer 110, which can be operated using one ormore input devices 118 of the computer 110 (e.g., keyboard and mouse). Note that while shown as separate devices inFIG. 1 , thedisplay device 120 and/orinput devices 118 can also be integrated with each other and/or with thecomputer 110, such as in a tablet computer. The toolpath generation program 116 b and thenumerical simulation program 116 c can be separate from the3D modeling program 116 a or integrated into the3D modeling program 116 a or into each other. Thus, a “program” as referenced herein can be a procedure within another program. - A
user 190 can interact with the3D modeling program 116 a to create a3D model 132 of a 3D object to be fabricated by the 3D extrusion printer 160. This can be done using known graphical user interface tools. In the example shown, the3D model 132 is a dragon figurine composed of multiple load bearing regions including the legs and the neck. In general, various suitable forms of3D model 132 can be created to represent a wide variety of 3D objects that can be manufactured using additive manufacturing systems and techniques. The3D model 132 can also contain mesh elements for performing FEA simulations. - During the design process, the
user 190 can specify a region of interest (ROI) 124 inUI 122 to mark a 3D region for estimating the physical properties of a portion of the 3D model contained within theROI 124. The size of theROI 124 can be changed to cover a small portion of the 3D model 132 (e.g., 1 mm3) to probe the microstructural properties of the modeled 3D object at theROI 124, or cover a large portion of the 3D model 132 (e.g., the entire model) to estimate the macroscale properties of the modeled 3D object. - Estimated
values 125 of the physical properties pertaining to theROI 124 are displayed in aphysical properties window 126. In some implementations, the estimatedvalues 125 can be updated at the request of theuser 190. In some implementations, the estimatedvalues 125 can be updated automatically based on various user inputs, such as moving or resizing theROI 124, or modifying the3D model 132. - The
toolpath generation program 116 b generatestoolpath data 130 based on the3D model 132. Thetoolpath data 130 is usable by the 3D extrusion printer 160 to fabricate a 3D printedobject 170 corresponding to the3D model 132. An example of toolpath data is G-code. Commercial and open source examples of thetoolpath generation program 116 b include Slic3r from Slic3r project (slic3r.org) and Cura from Ultimaker (Geldermalsen, Netherlands). - The inputs of the
toolpath generation program 116 b include the3D model 132, printer-hardware parameters, and microstructure-related parameters. The printer-hardware parameters can include nozzle diameter, filament diameter, extrusion temperature, and printer bed temperature. The microstructure-related parameters can include the layer height, number of perimeters, infill pattern, infill ratio, raster angle, raster to raster air gap, and perimeter to raster air gap. As a result, the generatedtoolpath data 130 contains the information and instructions necessary for the 3D extrusion printer 160 to fabricate the 3D printedobject 170. Examples of instructions contained within thetoolpath data 130 include printer head movement, stage movement, material feed rate, etc. for each layer. The 3D extrusion printer 160 follows the sets of instructions for each layer until all layers have been deposited to complete the fabrication of the 3D printedobject 170. - The
numerical simulation program 116 c is responsible for simulation of the physical properties of both the microstructures resulting from 3D printing, as well as the3D model 132 as a whole. Thenumerical simulation program 116 c, for example, can be a FEA simulation program for simulating mechanical, electrical, or thermal properties. Some other examples of thenumerical simulation program 116 c include a meshless method (element free Galerkin), XFEM (eXtended Finite element Method), BEM (Boundary Element Method), FVM (Finite Volume Method), DGFEM (Discontinuous Galerkin Finite Element Method), IGA (Isogeometric Analysis), and analytical methods based on volume fractions. Examples of mechanical properties to be simulated include elasticity, stiffness, flexibility, shear strength, tensile strength, yield strength, and specific weight. Examples of electrical properties to be simulated include the permittivity, permeability, resistivity, dielectric strength, inductance, and capacitance. Examples of thermal properties to be simulated include thermal conductivity and coefficient of thermal expansion. Commercial examples of thenumerical simulation program 116 c implementing FEA include NASTRAN from Autodesk Inc. (San Rafael, Calif.) and COMSOL Multiphysics from COMSOL Inc. (Burlington, Mass.). -
FIGS. 2A and 2B show a schematic diagram of a cross-section of a 3D printed object. A 3D printedobject 200 is a stack of individual 3D printed layers, including afirst layer 210, and asecond layer 220. For example, the 3D printedobject 200 shown is a half-sphere, in which individual layers are shaped to form an approximate surface of an ideal sphere. Referring toFIG. 2B , thefirst layer 210 has afirst microstructure 212, and thesecond layer 220 has asecond microstructure 222. Thefirst microstructure 212 is determined by the instructions contained within thetoolpath data 130. In some implementations, thefirst microstructure 212 can include one ormore perimeters 214 and aninfill pattern 216. For example, themicrostructure 212 shown has 2 perimeter filaments, and theinfill 216 is a raster pattern with its characteristics raster to raster air gap, perimeter to raster air gap, filament width, and raster angle. More generally, theinfill 216 can have various geometries including hexagonal lattices, concentric rectangles, Hilbert curves, and Archimedean chords. - The
second microstructure 222 is also determined by the instructions contained withintoolpath data 130. In some implementations, thesecond microstructure 222 can also include one ormore perimeters 224 and aninfill pattern 226. For example, themicrostructure 222 as shown has characteristics that are similar to themicrostructure 212, but theinfill 226 starts at a different corner than doesinfill 216, resulting in different filament orientations. In this example, theinfills object 200. In some cases, the 3D printedobject 200 can be formed by alternating betweenmicrostructures -
FIG. 3 is a flowchart showing an example of aprocess 300 to perform macroscale simulation of a 3D model of a 3D printed object using unit cell homogenization. Initially, at 310, theuser 190 uses the3D modeling program 116 a to create the3D model 132 as described in connection withFIG. 1 , or load a previously created3D model 132. - At 320, during any phase of the modeling process, the
user 190 may want to estimate the physical properties of a region of the3D model 132 by specifying theROI 124. For example, the specifiedROI 124 can be the load bearing neck region of the dragon figurine as shown, and theuser 190 may want to estimate the macroscale stiffness of the region. In another example, the specifiedROI 124 can be a metamaterial region containing dielectric materials and metallic traces, and theuser 190 may want to estimate the macroscale effective permittivity and effective permeability of the region. In another example, the specifiedROI 124 can be a critical heat conduction path to estimate the macroscale thermal resistance of the path. - At 330, once the
ROI 124 have been specified by theuser 190, thetoolpath generation program 116 b generates thetoolpath data 130 containing the instructions for 3D extrusion printer 160. In some implementations, thetoolpath data 130 is generated for theentire 3D model 132. In some implementations, thetoolpath data 130 is generated locally for the part of the3D model 132 covered by theROI 124. In some implementations, the toolpath data is generated hierarchically to avoid re-generation of toolpath data covering the regions that have not been modified since the last generation of the toolpath data. In some implementations, thetoolpath data 130 is generated only for the 3D printed layers that overlap with theROI 124. Such partial generation of toolpath data may be beneficial for improving the speed of the toolpath generation process, and hence improving the speed of the entire process for determining the macroscale properties. - At 340, one or
more unit cells 342 are defined for creating one or moreunit cell geometries 344, based on theROI 124 and thetoolpath data 130. Theunit cell 342 is a 3D region within theROI 124. Theunit cell 342 delineates a portion of the3D model 132 to be simulated using thenumerical simulation program 116 c to determine microstructural properties. In some implementations, theunit cell 342 can be defined by the3D modeling program 116 a. Theunit cell 342 is used to construct theunit cell geometry 344, which contains the microstructural information necessary for running a numerical simulation to determine theunit cell 342's microstructural properties. Theunit cell geometry 344 can be constructed by taking the union of the 3D region occupied by theunit cell 342 and thetoolpath data 130. For example, in a conformal method, filament geometries derived from thetoolpath data 130 can be implicitly added in the3D model 132 and converted into to finite element mesh through automatic 3D mesh generation. Another example is a voxel-based method. In the voxel-based method, a voxel domain can be created over the3D model 132. The voxel domain, for example, can be of a volumetric array of discrete voxels subdividing the3D model 132. The voxel domain can then be overlaid onto thetoolpath data 130. For each voxel, a check can be performed to determine whether a spatial overlap of a voxel with the filaments defined by thetoolpath data 130 exceeds a threshold (e.g., more than 1%, 5%, 10%, 20%, 50%, 75%, 99%). If the threshold is exceeded, the entire voxel is considered to be filled with the filament material, and vice versa. The resulting voxel domain consisting of filled and empty voxels can be directly translated to a finite element mesh. For example, each voxel can be translated 1-to-1 to a mesh element. - The
unit cell 342 can be defined in various ways. In some examples, theunit cell 342 defined by the3D modeling program 116 a is theROI 124. In some examples, theunit cell 342 is smaller than theROI 124. For example, theROI 124 may contain periodic or quasi-periodic microstructures. In such scenarios, according to the theory of unit cell homogenization, one or more representative volume elements (RVEs) can be determined, which can be used as the unit cells. The resulting unit cells can accurately represent the microstructural properties across theROI 124 without having to simulate theentire ROI 124. An illustrative example is a brick wall. A brick wall is made of individual bricks glued together by mortar. Because the brick and mortar stacking is repeated quasi-periodically throughout the entire wall, a small subsection of the wall containing just a few bricks glued by a mortar layer can be used as the unit cell, or a RVE, of the entire wall. Therefore, the entire wall made of bricks and mortar can be replaced by a wall of homogeneous material that has the microstructural properties of the unit cell (the ‘average’ property of a brick wall). - In some examples, the
ROI 124 may contain two ormore unit cells 342. This may be the case when the portion of the3D model 132 that falls within theROI 124 is not sufficiently uniform as informed by thetoolpath data 130. For example, theROI 124 may encompass both the infill region and the perimeter region of a3D model 132 as determined from thetoolpath data 130. A perimeter is typically one or more parallel filaments defining the outer surface of a 3D printedobject 170, and this filament configuration of the perimeter is typically substantially different from that of the infill region, where the filaments may be sparse and not parallel to each other. In this case, more than one unit cells should be created and positioned to cover both the perimeter region as well as the infill region to sufficiently describe, or sample, theROI 124. - In general, when two or more unit cells are needed to adequately sample the
ROI 124, several considerations exist for the number, locations, and distribution of the unit cells. One consideration is the mesh density of the3D model 132. In some implementations, mesh density of the3D model 132 can be used as a weighing criterion to distribute more unit cells in regions of high mesh density, and less unit cells in regions of low mesh density. Another consideration is the variation of microstructures from layer to layer. In some toolpath data, little variations occur from layer to layer (e.g., uniform hexagonal lattice over multiple layers). Given such toolpath data, more unit cells should be distributed along the plane of the layers to sample variations in microstructure, and less unit cells should be distributed to sample variations, if any, along the out-of-plane direction. - In some implementations, one unit cell can be placed at each of the mesh elements of the
3D model 132. In this case, a unit cell is placed at a representative location of each mesh element as determined using theory of unit cell homogenization. The microstructural properties of the respective unit cells can then be used to homogenize the volume covered by the respective mesh elements. In some examples, the representative location can be the geometrical center of the mesh element. - At 350, the
numerical simulation program 116 c receives one or moreunit cell geometries 344 and generates amicrostructural simulation 354, a numerical simulation for generatingmicrostructural properties 356, for each received unit cell geometry. In some implementations, the generation of themicrostructural simulation 354 can include recreating theunit cell geometry 344 within thenumerical simulation program 116 c, assigning of material properties to created geometries, generating mesh, and setting up an appropriate system of equations for the microstructural properties to be simulated. - The
numerical simulation program 116 c runs the generatedmicrostructural simulation 354 to calculate themicrostructural properties 356. In some implementations, quick completion of themicrostructural simulation 354 is desired to provide theuser 190 with the estimatedvalues 125 of the microstructural properties of theROI 124 in a quasi-real-time fashion. In such implementations, the number of mesh elements generated by thenumerical simulation program 116 b can to be reduced to decrease simulation run time. However, it should be noted that a tradeoff exists between the number of mesh elements and the accuracy of a FEA simulation. Therefore, thesystem 100 may askuser 190 to provide, for example, a desired level of accuracy or a maximum runtime of themicrostructural simulation 354 for improved user experience. - At 360, a
macroscale simulation 364 encompassing a part of, or theentire 3D model 132 as defined by theROI 124 is performed using themicrostructural properties 356 to determine themacroscale properties 366. Based on the unit cell homogenization technique, the estimatedmicrostructural properties 356 can be used to replace the microstructural details of the3D model 132 as provided by thetoolpath data 130 with one or more homogeneous materials having the estimatedmicrostructural properties 356. If the homogeneous materials are sufficiently representative of the microstructures of the 3D printedobject 170 over theROI 124 simulated in themacroscale simulation 364, the simulatedmacroscale properties 366 can be a good approximation of the true macroscale properties of the 3D printedobject 170 at reduced computational complexity, resource requirement, and time. Examples of macroscale properties that can be simulated are described above in connection with thenumerical simulation program 116 c andFIG. 1 . Additional examples of mechanical macroscale properties that can be simulated using unit cell homogenization simulation include Young's modulus, Poisson's ratio of elasticity, natural frequency, principal strain, principal stress, and measures of stress including Von Mises Stress, Tsai-Wu Failure index, and safe factor. In some implementations, weighted sampling point integration scheme (e.g., Gauss numerical integration, Cubature method) is used to integrate governing equations within a unit cell as part of themacroscale simulation 364. -
FIG. 4 is a schematic diagram of an example mesh element used in FEA. Amesh element 400 includes Gaussian integration points 410 andunit cells 420. In some implementations, theunit cells 420 are each centered on the Gaussian integration points 410. In the example shown, themesh element 400 is a hexahedron, but other types are possible, including a pyramid, a triangular prism, and a tetrahedron. - Some physical properties, such as the elasticity matrix, may require more than one unit cell geometry to be created and simulated within a mesh element for accurate determination. An approximation of a definite integral of a mathematical function can be stated as a weighted sum of the evaluated function values at the Gaussian integration points 410 within the
mesh element 400. The definite integral of a mathematical function is needed for evaluation of certain physical properties, such as elasticity. For example, in evaluating the elasticity, the local stress and strain are evaluated at the Gaussian integration points 410 to provide the most accurate estimation of the elasticity matrix. To evaluate local stress and strain at the Gaussian integration points 410, the microstructural properties of the Gaussian integration points 410 are needed. Therefore, a unit cell can be placed at each of the Gaussian integration points 410 to determine the microstructural properties of the respective Gaussian integration points 410 using the unit cell homogenization technique. -
FIG. 5 is a flowchart showing an example of aprocess 500 to simulate the macroscale properties of a 3D printed object. Initially, a 3D model of a 3D object to be manufactured by a 3D extrusion printer is obtained 510. The3D model 132 of a 3D object can be obtained from the3D modeling program 116 a in various ways. For example, the3D modeling program 116 a can export the3D model 132 into a portable digital representation. Examples of the portable digital representation of the 3D model include stereo lithography ‘.STL’ files, wavefront ‘.OBJ’ files, 3D manufacturing format ‘0.3MF’ files, and additive manufacturing file format ‘.AMF’ files. These digital representations can contain information describing the surface geometry of the3D model 132 and the physical properties of the material constituting the modeled 3D object. - A toolpath data for 3D extrusion printer is then obtained 520. The entire toolpath data or a partial toolpath data can be obtained as described above in connection with
FIG. 3 at 330. In some implementations, the generation of thetoolpath data 130 can be automated. For example, the3D modeling program 116 a can export a suitable output file (e.g., .STL, .OBJ, or .AMF file), which can be passed onto thetoolpath generation program 116 b using an application programming interface of thetoolpath generation program 116 b. Based on the received commands and data files, thetoolpath data 130 can be generated without intervention of theuser 190. - A unit cell geometry is then created 530. Creating 530 is performed using at least one limitation on one or
more unit cells 342 of theunit cell geometry 344 with respect to an aspect of the 3D extrusion printer's build of the 3D object. In some implementations, a unit cell encloses at least three layers of the 3D extrusion printer's build of the 3D object. A unit cell encompassing at least three layers contains both the upper interface and the lower interface of the middle layer as sandwiched by the top and the bottom layer. The inclusion of both interfaces improves the accuracy of the microstructural properties simulated from the unit cell geometry. - In some implementations, the layer height parameter of the 3D extrusion printer, for example, is used to determine the minimum dimensions of a unit cell in corresponding design units (e.g., mm or mil) of the
3D modeling program 116 a. For example, a typical layer height ranges from 0.2-0.3 mm, leading to a minimum unit cell height of 0.6-0.9 mm. In some implementations, the minimum width and length of the unit cell are equal to the minimum height. - In some implementations, a unit cell has an upper limit in size. For example, the upper limit can come from consideration of locality, a condition used in micromechanics. For the locality condition to hold, the unit cell needs to be smaller than the mesh element that it represents. If the unit cell is larger than the mesh element itself, the unit cell is no longer local to the mesh element and ceases to be representative of the mesh element, violating the locality condition. In some implementations, when multiple Gaussian integration points 410 of the
mesh element 400 are used for evaluation of physical properties, the maximum volume of each unit cell is the volume of themesh element 400 divided by the number of Gaussian integration points of the mesh element to avoid overlapping of unit cells and violation of the locality condition. For example, if an 8-point Gauss integration scheme is used for a macroscale simulation, the unit cell volume should be less than ⅛ of the volume of the macroscale mesh element. - The location of a unit cell can be determined using various criteria. In some implementations, a unit cell is placed so that a substantial fraction of the unit cell is within the containing envelope of the 3D model. For example, if a unit cell is placed near the containing envelope (e.g., a surface) of the 3D model and a large portion of the unit cell falls outside of the 3D model, the unit cell should be moved towards the inside of the 3D model to have more overlap with the 3D model (e.g., more than 1%, 5%, 10%, 20%, 50%, 75%, 99% in volume overlap). In some examples, a unit cell may fall within a sparse region of the 3D model as informed by the toolpath data 130 (e.g., region of low infill) where little or no microstructures exist within the unit cell. In such cases, the unit cell can be moved or enlarged so that the unit cell has greater overlap with the microstructures (e.g., more than 1%, 5%, 10%, 20%, 50%, 75%, 99% in volume overlap).
- A numerical simulation is generated 540. The numerical simulation can be generated as described above in connection with
FIG. 3 at 350. In some cases, two or more numerical simulations can be run in parallel based on the generated numerical simulation. For example, when calculating linear elasticity, the effective elasticity tensor can be identified from the strain energies of the microstructure under specific boundary conditions. The elasticity tensor can have 6 basic strain modes, which are independent from one another. Therefore, numerical simulations for the 6 basic strain mode can be run in parallel. - Macroscale properties are estimated 550. The macroscale properties can be calculated as described above in connection with
FIG. 3 at 360. - A check is performed to determine 560 whether the estimated macroscale properties meet the user's requirements. In some implementations, the
user 190 can provide thesystem 100 with one or more requirements on the microstructural properties or the macroscale properties of the 3D printed object at theROI 124. Based on estimating 550, if the user requirements are met, success is reported 565 to theuser 190. When the user requirements are not met, failure is reported 570, and a retry of the toolpath generation may be initiated by returning to obtaining 520 different toolpath data for 3D extrusion printer generated using different parameters. - In some implementations, the retry of the toolpath generation can be automated. For example, printer-hardware parameters and microstructure-related parameters can be automatically modified during the retry of obtaining 520. In some implementations, the parameters to be changed can be determined automatically by the
3D modeling program 116 a. The automatic determination can be based on various rules derived from first principles and designer experience to make the automatic iteration more efficient. For example, theprogram 116 a may increase the infill density to increase stiffness, or change infill pattern to change the anisotropy as derived from first principles. Theprocess 500 can be repeated until the requirements on the macroscale properties are satisfied. - A computer program (also known as a program, software, software application, script, or code) can be written in any form of programming language, including compiled or interpreted languages, declarative or procedural languages, and it can be deployed in any form, including as a stand alone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program does not necessarily correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub programs, or portions of code). A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.
- The processes and logic flows described in this specification can be performed by one or more programmable processors executing one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).
- Processors suitable for the execution of a computer program include, by way of example, both general and special purpose microprocessors, and any one or more processors of any kind of digital computer. Generally, a processor will receive instructions and data from a read only memory or a random access memory or both. The essential elements of a computer are a processor for performing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto optical disks, or optical disks. However, a computer need not have such devices. Moreover, a computer can be embedded in another device, e.g., a mobile telephone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a Global Positioning System (GPS) receiver, or a portable storage device (e.g., a universal serial bus (USB) flash drive), to name just a few. Devices suitable for storing computer program instructions and data include all forms of non volatile memory, media and memory devices, including by way of example semiconductor memory devices, e.g., EPROM (Erasable Programmable Read-Only Memory), EEPROM (Electrically Erasable Programmable Read-Only Memory), and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto optical disks; and CD ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
- To provide for interaction with a user, embodiments of the subject matter described in this specification can be implemented on a computer having a display device, e.g., LCD (liquid crystal display), OLED (organic light emitting diode) or other monitor, for displaying information to the user and a keyboard and a pointing device, e.g., a mouse or a trackball, by which the user can provide input to the computer. Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback, e.g., visual feedback, auditory feedback, or tactile feedback; and input from the user can be received in any form, including acoustic, speech, or tactile input.
- The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other. Embodiments of the subject matter described in this specification can be implemented in a computing system that includes a back end component, e.g., as a data server, or that includes a middleware component, e.g., an application server, or that includes a front end component, e.g., a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the subject matter described is this specification, or any combination of one or more such back end, middleware, or front end components. The components of the system can be interconnected by any form or medium of digital data communication, e.g., a communication network. Examples of communication networks include a local area network (“LAN”) and a wide area network (“WAN”), an inter-network (e.g., the Internet), and peer-to-peer networks (e.g., ad hoc peer-to-peer networks).
- While this specification contains many implementation details, these should not be construed as limitations on the scope of the invention or of what may be claimed, but rather as descriptions of features specific to particular embodiments of the invention. Certain features that are described in this specification in the context of separate embodiments can also be implemented in combination in a single embodiment. Conversely, various features that are described in the context of a single embodiment can also be implemented in multiple embodiments separately or in any suitable subcombination. Moreover, although features may be described above as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a subcombination or variation of a subcombination.
- Similarly, while operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results. In certain circumstances, multitasking and parallel processing may be advantageous. Moreover, the separation of various system components in the embodiments described above should not be understood as requiring such separation in all embodiments, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.
- Thus, particular embodiments of the invention have been described. Other embodiments are within the scope of the following claims. In addition, the actions recited in the claims can be performed in a different order and still achieve desirable results.
Claims (18)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US15/970,291 US11130292B2 (en) | 2017-05-08 | 2018-05-03 | Estimating physical property of 3D printed parts |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201762502953P | 2017-05-08 | 2017-05-08 | |
US15/970,291 US11130292B2 (en) | 2017-05-08 | 2018-05-03 | Estimating physical property of 3D printed parts |
Publications (2)
Publication Number | Publication Date |
---|---|
US20180319087A1 true US20180319087A1 (en) | 2018-11-08 |
US11130292B2 US11130292B2 (en) | 2021-09-28 |
Family
ID=64014041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/970,291 Active 2040-03-21 US11130292B2 (en) | 2017-05-08 | 2018-05-03 | Estimating physical property of 3D printed parts |
Country Status (1)
Country | Link |
---|---|
US (1) | US11130292B2 (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190026953A1 (en) * | 2016-05-16 | 2019-01-24 | Hewlett-Packard Development Company, L.P. | Generating a shape profile for a 3d object |
US20190121329A1 (en) * | 2017-10-24 | 2019-04-25 | General Electric Company | Apparatus and method for integrating multiple products in a single process |
US10489525B2 (en) * | 2016-02-17 | 2019-11-26 | Toyota Motor Engineering & Manufacturing North America, Inc. | Systems and methods for composite thermal interface material microstructure property prediction |
US20200130281A1 (en) * | 2017-06-20 | 2020-04-30 | Toray Engineering Co., Ltd. | Method for predicting strength of structure, method for modeling structure, support method for additive manufacturing of structure, and recording medium |
US20200387652A1 (en) * | 2019-06-07 | 2020-12-10 | Intact Solutions, Inc. | Computational Modeling of Procedural Language Descriptors Prior to Manufacture |
US20210046545A1 (en) * | 2019-08-14 | 2021-02-18 | Fuji Xerox Co., Ltd. | Three-dimensional shape data generation apparatus, three-dimensional modeling apparatus, and non-transitory computer readable medium storing three-dimensional shape data generation program |
US10928805B2 (en) * | 2018-11-15 | 2021-02-23 | Lawrence Livermore National Security, Llc | Additive manufacturing system and method having toolpath analysis |
US10981375B2 (en) * | 2016-10-11 | 2021-04-20 | Lawrence Livermore National Security, Llc | In situ monitoring of direct ink write process using electromagnetic resonant coupling |
US20210370584A1 (en) * | 2020-05-06 | 2021-12-02 | Eaton Intelligent Power Limited | Method of additively manufacturing transparent lenses for luminaries |
US11250186B2 (en) * | 2019-10-04 | 2022-02-15 | Autodesk, Inc. | Machine learning approach to piecewise linear interface construction |
US20220288861A1 (en) * | 2019-10-30 | 2022-09-15 | Hewlett-Packard Development Company, L.P. | Print operation using plural sets of process values |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160096318A1 (en) * | 2014-10-03 | 2016-04-07 | Disney Enterprises, Inc. | Three dimensional (3d) printer system and method for printing 3d objects with user-defined material parameters |
US9895845B2 (en) * | 2015-02-16 | 2018-02-20 | Arevo Inc. | Method and a system to optimize printing parameters in additive manufacturing process |
US10353378B2 (en) * | 2016-08-18 | 2019-07-16 | Wisconsin Alumni Research Foundation | Homogenization of material properties in additively manufactured structures |
-
2018
- 2018-05-03 US US15/970,291 patent/US11130292B2/en active Active
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10489525B2 (en) * | 2016-02-17 | 2019-11-26 | Toyota Motor Engineering & Manufacturing North America, Inc. | Systems and methods for composite thermal interface material microstructure property prediction |
US11043042B2 (en) * | 2016-05-16 | 2021-06-22 | Hewlett-Packard Development Company, L.P. | Generating a shape profile for a 3D object |
US20190026953A1 (en) * | 2016-05-16 | 2019-01-24 | Hewlett-Packard Development Company, L.P. | Generating a shape profile for a 3d object |
US10981375B2 (en) * | 2016-10-11 | 2021-04-20 | Lawrence Livermore National Security, Llc | In situ monitoring of direct ink write process using electromagnetic resonant coupling |
US20200130281A1 (en) * | 2017-06-20 | 2020-04-30 | Toray Engineering Co., Ltd. | Method for predicting strength of structure, method for modeling structure, support method for additive manufacturing of structure, and recording medium |
US11465356B2 (en) * | 2017-06-20 | 2022-10-11 | Toray Engineering Co., Ltd. | Method for predicting strength of structure, method for modeling structure, support method for additive manufacturing of structure, and recording medium |
US20190121329A1 (en) * | 2017-10-24 | 2019-04-25 | General Electric Company | Apparatus and method for integrating multiple products in a single process |
US10928805B2 (en) * | 2018-11-15 | 2021-02-23 | Lawrence Livermore National Security, Llc | Additive manufacturing system and method having toolpath analysis |
US20200387652A1 (en) * | 2019-06-07 | 2020-12-10 | Intact Solutions, Inc. | Computational Modeling of Procedural Language Descriptors Prior to Manufacture |
US20210046545A1 (en) * | 2019-08-14 | 2021-02-18 | Fuji Xerox Co., Ltd. | Three-dimensional shape data generation apparatus, three-dimensional modeling apparatus, and non-transitory computer readable medium storing three-dimensional shape data generation program |
US11886781B2 (en) * | 2019-08-14 | 2024-01-30 | Fujifilm Business Innovation Corp. | Three-dimensional shape data generation apparatus, three-dimensional modeling apparatus, and non-transitory computer readable medium storing three-dimensional shape data generation program |
US11250186B2 (en) * | 2019-10-04 | 2022-02-15 | Autodesk, Inc. | Machine learning approach to piecewise linear interface construction |
US20220288861A1 (en) * | 2019-10-30 | 2022-09-15 | Hewlett-Packard Development Company, L.P. | Print operation using plural sets of process values |
US20210370584A1 (en) * | 2020-05-06 | 2021-12-02 | Eaton Intelligent Power Limited | Method of additively manufacturing transparent lenses for luminaries |
Also Published As
Publication number | Publication date |
---|---|
US11130292B2 (en) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11130292B2 (en) | Estimating physical property of 3D printed parts | |
CN109313670B (en) | Method and system for generating lattice suggestions in computer aided design applications | |
US10353378B2 (en) | Homogenization of material properties in additively manufactured structures | |
US8463584B2 (en) | Interactive filling simulation on 3D injection molding models | |
CN108351906B (en) | System and method for modeling of a component having a lattice structure | |
JP7417815B2 (en) | Computer-aided generative design with global thickness control to accelerate manufacturing and structural performance | |
US10599791B2 (en) | Method and device for digital reconstruction of an elementary volume representing a microstructure of composite material | |
Bacciaglia et al. | A systematic review of voxelization method in additive manufacturing | |
Liu et al. | A survey of modeling and optimization methods for multi-scale heterogeneous lattice structures | |
US11947333B2 (en) | Dual lattice representation for crash simulation and manufacturing | |
US10422925B2 (en) | 2.5D stadia meshing | |
Kilingar et al. | Computational generation of open-foam representative volume elements with morphological control using distance fields | |
GB2529957B (en) | 2.75D meshing algorithm | |
CN115114706A (en) | Slope design method, device, equipment and computer readable storage medium | |
Liang et al. | Numerical quantification of the impact of microstructure on the mechanical behavior of particulate Al/SiC composites in 2D | |
US10436939B2 (en) | Lofting algorithm for discrete network meshing | |
Di Paola et al. | Algorithmic-Generative Architecture | |
Galicia | 3D printing speed optimization by minimizing void paths | |
JP2019111779A (en) | Three dimension printer control method, three dimension printer control program, and information processing apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
AS | Assignment |
Owner name: AUTODESK, INC., CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:EOM, JAESUNG;BURLA, RAVI KUMAR;REEL/FRAME:046104/0758 Effective date: 20170505 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: AUTODESK, INC., CALIFORNIA Free format text: CHANGE OF ADDRESS FOR ASSIGNEE;ASSIGNOR:AUTODESK, INC.;REEL/FRAME:061572/0061 Effective date: 20220711 |