CN111460709A - Method for predicting temperature distribution and deformation of part in fused deposition manufacturing process - Google Patents
Method for predicting temperature distribution and deformation of part in fused deposition manufacturing process Download PDFInfo
- Publication number
- CN111460709A CN111460709A CN202010235780.3A CN202010235780A CN111460709A CN 111460709 A CN111460709 A CN 111460709A CN 202010235780 A CN202010235780 A CN 202010235780A CN 111460709 A CN111460709 A CN 111460709A
- Authority
- CN
- China
- Prior art keywords
- unit
- thermal analysis
- printed
- temperature
- finite element
- 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.)
- Withdrawn
Links
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a method for predicting temperature distribution and deformation of a part in a fused deposition manufacturing process. The method for predicting the temperature distribution and deformation of the part in the fused deposition manufacturing process comprises the following steps: defining the type of a thermal analysis unit and thermophysical parameters of a material used when the part is printed, establishing a simulation model according to the shape and size of the part to be printed, carrying out grid division on the simulation model according to the set type of the thermal analysis unit to obtain a finite element model comprising nodes and units, activating the units of the finite element model, and applying initial conditions and boundary conditions on the unit nodes to obtain the temperature distribution of the part to be printed in the fused deposition molding process and a corresponding thermal analysis result file, converting the thermal analysis unit into a structural unit, and defining the mechanical properties of the material to be used to obtain the mechanical model of the part to be printed. The invention can predict the warping deformation of the part.
Description
Technical Field
The invention relates to the technical field of fused deposition manufacturing, in particular to a method for predicting temperature distribution and deformation of a part in the fused deposition manufacturing process.
Background
In recent years, Fused Deposition Modeling (FDM) has been widely used due to its characteristics of low cost, easy maintenance, wide range of molding materials, low price, small occupied space, and the like. When the model is formed by the FDM technology, the adopted material is basically thermoplastic material, the material is heated at a nozzle, the molten material is extruded from the nozzle and is stacked layer by layer along with the movement of the nozzle, and then the molten material is cooled gradually, and finally the required model is formed. In the molding process, the material undergoes multiple uneven heating and cooling cycles, which causes different curing times at various positions of the model and uneven distribution of temperature field and stress field, and further causes buckling deformation, even failure phenomena such as interlayer debonding and cracking, and the molding precision is affected.
In order to avoid part failure as much as possible and reduce the buckling deformation of parts, a test block is generally required to be printed before formal production for trial and error test, so that more reasonable printing parameters are obtained. However, the process is limited by capital and time costs, and the experiment is limited by measuring means, so that the distribution of temperature and stress of the part in the forming process cannot be obtained well, and the warping deformation of the whole part cannot be predicted, and therefore, improvement is urgently needed.
Disclosure of Invention
In view of the above-mentioned drawbacks of the prior art, an object of the present invention is to provide a method for predicting the temperature distribution and deformation of a part in a fused deposition manufacturing process, which is used to solve the problems of the prior art that the distribution of the temperature and stress of the part in the molding process cannot be obtained well, and the warpage deformation of the whole part cannot be predicted.
To achieve the above and other related objects, the present invention provides a method for predicting a temperature distribution and deformation of a part in a fused deposition manufacturing process, the method comprising:
defining the type of a thermal analysis unit and thermophysical parameters of a material used when printing a part;
establishing a simulation model according to the shape and size of a part to be printed, and carrying out meshing on the simulation model according to the set type of a thermal analysis unit to obtain a finite element model comprising nodes and units;
activating the units of the finite element model, and applying initial conditions and boundary conditions on unit nodes to obtain the temperature distribution of the part to be printed in the fused deposition molding process and a corresponding thermal analysis result file;
converting the thermal analysis unit into a structural unit, and defining mechanical properties of a used material to obtain a mechanical model of the part to be printed;
and adding displacement constraint at the bottom of the mechanical model, loading the thermal analysis result file, and loading the temperature of the node of the finite element model as a body load onto the corresponding node of the finite element model to obtain the stress distribution and deformation of the part to be printed in the fused deposition molding process.
In an embodiment of the present invention, the creating a simulation model according to the shape and size of the part to be printed, and the meshing the simulation model according to the set thermal analysis unit type includes:
establishing a simulation model according to the shape and the size of the part to be printed;
obtaining the printing line width and the layer thickness by using a finite volume method and combining the layering thickness and the printing line width adopted by the part to be printed;
and taking the printing line width as the length and the width of the thermal analysis unit, taking the layer thickness as the height of the thermal analysis unit, and performing grid division on the simulation model according to the set type of the thermal analysis unit.
In an embodiment of the invention, the step of activating the elements of the finite element model comprises:
killing all the elements of the finite element model by using a life-death element technology, and activating the Nth element in the finite element model according to a printing path, wherein N represents any one element in the finite element model;
applying initial conditions and boundary conditions on the Nth unit of the finite element model, and solving the temperature of the temperature field according to a set load loading mode and load steps to delete the initial temperature load of the Nth unit of the finite element model;
activating the (N + 1) th unit of the finite element model, and solving the temperature of the temperature field to delete the initial temperature load of the (N + 1) th unit of the finite element model; until all elements of the finite-element model are activated.
In an embodiment of the present invention, the step of obtaining the temperature distribution of the part to be printed in the fused deposition modeling process and the corresponding thermal analysis result file according to the initial condition and the boundary condition applied to the unit node includes:
and solving the temperature field according to a set load loading mode and load steps, and applying initial conditions and boundary conditions on the activated nodes of each unit to obtain the temperature distribution of the part to be printed and a corresponding thermal analysis result file.
In an embodiment of the present invention, the thermophysical parameters of the used material include density, specific heat capacity, thermal conductivity, and overall heat transfer coefficient.
In an embodiment of the present invention, the set load loading manner is a step-wise load loading, and the step size of the load step is set according to the printing speed and the length of each unit.
In an embodiment of the invention, the initial conditions comprise an initial temperature comprising one or more of a nozzle temperature, a forming chamber temperature, a sole plate temperature.
In one embodiment of the present invention, the boundary condition includes heat transfer between the molded material, the molded material and air in the molding chamber transfer heat by thermal convection and thermal radiation, and the molded material gradually cools from a molten state to a solid state with latent heat of phase change of the phase.
The invention provides a system for predicting the temperature distribution and deformation of a part in a fused deposition manufacturing process, which comprises the following steps:
the processor unit is used for defining the type of the thermal analysis unit and thermophysical parameters of materials used when the parts are printed;
the system comprises a simulation model establishing unit, a thermal analysis unit and a thermal analysis unit, wherein the simulation model establishing unit is used for establishing a simulation model according to the shape and the size of a part to be printed and carrying out grid division on the simulation model according to the set type of the thermal analysis unit so as to obtain a finite element model comprising nodes and units;
the activator unit is used for activating the units of the finite element model and applying initial conditions and boundary conditions on unit nodes to obtain the temperature distribution of the part to be printed in the fused deposition molding process and a corresponding thermal analysis result file;
the mechanical model establishing unit is used for converting the thermal analysis unit into a structural unit and defining the mechanical properties of the used materials so as to obtain a mechanical model of the part to be printed;
and the displacement constraint unit is used for adding displacement constraint at the bottom of the mechanical model, loading the thermal analysis result file, and loading the temperature of the node of the finite element model as a body load onto the corresponding node of the finite element model to obtain the stress distribution and deformation of the part to be printed in the fused deposition molding process.
The invention provides electronic equipment which comprises a processor and a memory, wherein the memory stores program instructions, and the processor runs the program instructions to realize the prediction method of the temperature distribution and deformation of the part in the fused deposition manufacturing process.
As described above, the method for predicting the temperature distribution and deformation of a part in the fused deposition manufacturing process according to the present invention has the following advantages:
the method for predicting the temperature distribution and deformation of the part in the fused deposition manufacturing process can obtain the temperature field and the stress field of the part in the forming process by combining thermal analysis and mechanical analysis.
The invention establishes a temperature and deformation prediction system, and provides technical support for optimizing printing parameters, thereby reducing the warping deformation generated by printing parts by using a process fused deposition manufacturing technology.
The method can calculate the temperature distribution, the stress distribution and the deformation of the part in the fused deposition manufacturing process at one time, predict the temperature distribution and the final warping deformation in the part forming process, and avoid the defects of time and labor waste and error proneness in the manual modeling and solving by using a user interface for multiple times.
Drawings
Fig. 1 is a flowchart illustrating a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure.
Fig. 2 is a flowchart illustrating a step S2 of a method for predicting a temperature distribution and deformation of a part in a fused deposition manufacturing process in fig. 1 according to an embodiment of the present disclosure.
Fig. 3 is a flowchart illustrating a step S3 of a method for predicting a temperature distribution and deformation of a part in a fused deposition manufacturing process in fig. 1 according to an embodiment of the present disclosure.
Fig. 4 is a flowchart illustrating a step S3 of a method for predicting a temperature distribution and deformation of a part in a fused deposition manufacturing process in fig. 1 according to an embodiment of the present disclosure.
Fig. 5 is a flowchart illustrating a step S5 of a method for predicting a temperature distribution and deformation of a part in a fused deposition manufacturing process in fig. 1 according to an embodiment of the present disclosure.
FIG. 6 is a schematic block diagram illustrating a system for predicting temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure.
Fig. 7 is a schematic block diagram of a structure of an electronic device according to an embodiment of the present disclosure.
FIG. 8 is a schematic diagram of an application of a system for predicting temperature distribution and distortion of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure.
Fig. 9 is a schematic model diagram of a method for predicting temperature distribution and deformation of a part in a fused deposition manufacturing process according to an embodiment of the present disclosure.
Fig. 10 is a schematic scanning diagram illustrating a method for predicting temperature distribution and deformation of a part in a fused deposition manufacturing process according to an embodiment of the present disclosure.
Fig. 11(a) and (b) are temperature distributions and stress distribution diagrams of a part to be printed, which are predicted by a method for predicting the temperature distribution and deformation of the part in a fused deposition manufacturing process according to an embodiment of the present application, just before printing and forming.
Fig. 12(a) and (b) are temperature distributions and stress distribution diagrams of a part to be printed when the part to be printed is cooled for 20s, which are predicted by the prediction method for the temperature distribution and deformation of the part in the fused deposition manufacturing process provided by the embodiment of the present application.
Description of the element reference numerals
1 prediction system for temperature distribution and deformation of parts in fused deposition manufacturing process
2 Server
3 third layer mesh model
4 second layer grid model
5 first layer mesh model
10 processor unit
20 simulation model establishing unit
30 activator unit
40 mechanical model building unit
50 displacement restraint unit
60 processor
70 memory
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It is to be noted that the features in the following embodiments and examples may be combined with each other without conflict.
It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention, and the drawings only show the components related to the present invention rather than the number, shape and size of the components in actual implementation, and the type, quantity and proportion of the components in actual implementation may be changed freely, and the layout of the components may be more complicated.
Referring to fig. 1, 2, 3, 4, and 5, fig. 1 is a flow chart of a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure, fig. 2 is a flow chart of a step S2 of a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure, fig. 3 is a flow chart of a step S3 of a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure, fig. 4 is a flow chart of a step S3 of a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure, fig. 5 is a flow chart of a step S5 of a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure, fig. 5 is a flow chart of a method for predicting a temperature distribution and deformation of a part during a fused deposition manufacturing process according to a embodiment of the present disclosure, a temperature distribution and deformation of a portion of a fused deposition process, a heat transfer coefficient of a heat transfer material L, a heat transfer coefficient of a heat sink L, a heat transfer material, a heat transfer coefficient of a heat transfer heat sink, a heat transfer heat analysis unit 3580, a heat transfer heat analysis unit 3580, a heat analysis heat transfer heat analysis heat unit 3580, a heat analysis heat transfer heatThe material is temperature dependent in real time. The comprehensive heat transfer coefficient is obtained by comprehensively considering heat convection and heat radiation, and the comprehensive heat transfer coefficient h is as follows: h ═ σ (T)p 2+Tc 2)(Tp+Tc)+hcWherein σ represents Stefan-Boltzmann constant, and has a value of 5.67 × 10-8W/m2·K4Is the blackness of the material, TpIndicating the temperature, T, of the partcDenotes the ambient temperature, hcRepresenting the natural convective heat transfer coefficient. S2, establishing a simulation model according to the shape and size of the part to be printed, and carrying out grid division on the simulation model according to the set thermal analysis unit type to obtain a finite element model comprising nodes and units. Specifically, the step of establishing a simulation model according to the shape and size of the part to be printed in step S2, and the step of meshing the thermal analysis unit includes: and S21, establishing a simulation model according to the shape and the size of the part to be printed. And S22, obtaining the printing line width and the layer thickness by using a finite volume method and combining the layering thickness and the printing line width adopted by the part to be printed. And S23, taking the printing line width as the length and width of the thermal analysis unit, taking the layer thickness as the height of the thermal analysis unit, and performing grid division on the simulation model according to the set type of the thermal analysis unit. And S3, activating the elements of the finite element model, and applying initial conditions and boundary conditions on element nodes to obtain the temperature distribution of the part to be printed in the fused deposition modeling process and a corresponding thermal analysis result file. Specifically, the step of activating the elements of the finite element model in step S3 includes: and S31, killing all the elements of the finite element model by using a life-death element technology, and activating the Nth element in the finite element model according to the printing path, wherein N represents any one element in the finite element model. And S32, applying initial conditions and boundary conditions on the Nth unit of the finite element model, and solving the temperature of the temperature field according to the set load loading mode and the load step to delete the initial temperature load of the Nth unit of the finite element model. S33, activating the (N + 1) th unit of the finite element model, and solving the temperature of the temperature fieldSolving to delete the initial temperature load of the (N + 1) th unit of the finite element model; until all elements of the finite-element model are activated. S4, converting the thermal analysis unit into a structural unit, and defining the mechanical properties of the used material to obtain a mechanical model of the part to be printed. And S5, adding displacement constraint at the bottom of the mechanical model, loading the thermal analysis result file, and loading the temperature of the node of the finite element model as a body load onto the corresponding node of the finite element model to obtain the stress distribution and deformation of the part to be printed in the fused deposition molding process.
Referring to fig. 1, fig. 2, fig. 3, fig. 4, and fig. 5, before activating the elements of the finite element model in step S3, it is necessary to sequentially activate each divided element by using a living and dead element technique according to a selected scanning mode, and the temperature distribution of the part to be printed and a corresponding thermal analysis result file in the fused deposition molding process can be obtained through, but not limited to, transient thermal analysis. The stress distribution and deformation of the part to be printed in the fused deposition molding process can be obtained by solving through a finite element method. The method can be used for dividing the grids of the built simulation model by combining the layered thickness and the printing line width adopted when the part is actually printed by using a finite volume method, taking the line width as the length and the width of the cell and taking the layer thickness as the height of the cell. The life and death unit technology means that if materials are added or deleted in the model, corresponding units in the model exist or disappear. The "death" of a cell is not removed from the model, but rather its stiffness matrix or conduction matrix or other characteristic matrix is multiplied by a small factor (ESTIF) which defaults to 1.0e-6Setting the unit load, the mass, the rigidity, the specific heat and the like to be 0, excluding the mass and the energy of a dead unit from the solving result of the model, and setting the strain of the unit to be 0 while killing; the 'birth' of the unit is not to add a new unit in the model, but to reactivate the unit, and the rigidity, the mass, the unit load and the like of the unit can restore the original values after the unit is activated. The load loading mode can be but is not limited to loading the load in a step mode, and the step length of the load stepThe initial condition is the temperature at the beginning of printing, the initial temperature is known because the machine preheats according to the temperature set by the user before printing with FDM (fused deposition modeling), the initial temperature includes but is not limited to one or more of nozzle temperature, molding chamber temperature and bottom plate temperature, the boundary condition is that heat is transferred between the molded material by heat conduction, the molded material and the air in the molding chamber transfer heat by heat convection and heat radiation, and phase change exists during the gradual cooling of the material from the molten state to the solid state, namely phase change exists, the molded material is processed by using a specific heat capacity mutation method, namely the effect of replacing latent heat by an abrupt change of latent heat capacity of the material, the latent heat is calculated by the enthalpy of the material in YS software, the unit of the specific heat capacity is KJ/Kg, the calculation formula is as H ═ c (T) dT, wherein H represents the value of the material ρ c, T represents the density of the material, T represents the temperature of the initial temperature, and the density of the molding material, and the unit of the initial temperature and the initial condition is set as a temperature after the initial load command, the activation of the activation unit, the activation of the molding process, and the activation of the heating unit, the activation of the molding process, and the unit, the activation of the molding process, the activation of the molding process of the activation of the molding process, the activation of the molding process, the molding process of the activation of the molding process, the activation ofThe method comprises the steps of S317, entering a cooling stage after printing is completed, selecting surface units of a model, reapplying boundary conditions, S318, setting cooling time of the part to be printed, solving according to set boundary conditions, same load step size and load sub steps, completing thermal analysis of the whole part from molding to cooling, S319, entering a post-processing stage after solving is completed, obtaining relevant information such as temperature distribution and temperature gradient of the part to be printed at each moment in the fused deposition molding process, converting the thermal analysis unit into a structural unit in step S4 to convert an SO L ID70 unit into an SO L ID185 unit, adding displacement constraint to the bottom of the mechanical model through a displacement constraint unit 50 in step S5 to obtain a displacement constraint to obtain a distribution of the part to be printed in the fused deposition molding process, setting the displacement constraint to be a three-step load step S511, reading a three-step load distribution, reading a three-step load step S82, reading a three-step load step S, a three-step load step S82, and a three-step load step S82.
The thermal deformation characteristics of the parts during the fused deposition manufacturing process are defined by a thermal expansion coefficient and a thermal deformation coefficient of the parts during the fused deposition manufacturing process, and are calculated by a thermal expansion coefficient calculation unit, a thermal expansion coefficient calculation unit, a thermal calculation unit, a thermal calculation unit, a thermal calculation unit, a thermal calculation unit, a thermal calculation unit.
Referring to fig. 6, 7 and 8, fig. 6 is a schematic structural diagram of a system for predicting temperature distribution and deformation of a part in a fused deposition manufacturing process according to an embodiment of the present disclosure. Fig. 7 is a schematic block diagram of a structure of an electronic device according to an embodiment of the present disclosure. FIG. 8 is a schematic diagram of an application of a system for predicting temperature distribution and distortion of a part during a fused deposition manufacturing process according to an embodiment of the present disclosure. Similar to the principle of the method for predicting the temperature distribution and deformation of a part in a fused deposition manufacturing process of the present invention, the present invention also provides a system for predicting the temperature distribution and deformation of a part in a fused deposition manufacturing process, which includes, but is not limited to, a processor unit 10, a simulation model establishing unit 20, an activator unit 30, a mechanical model establishing unit 40, and a displacement constraining unit 50. The processor unit 10, the simulation model building unit 20, the activator unit 30, the mechanical model building unit 40 and the displacement constraining unit 50 are part of a central processor. The processor unit 10 is used to define the type of thermal analysis unit and the thermophysical parameters of the material used when printing the part. The simulation model establishing unit 20 is configured to establish a simulation model according to the shape and size of the part to be printed, and perform mesh division on the simulation model according to the set type of the thermal analysis unit to obtain a finite element model including nodes and units. The activator unit 30 is configured to activate the elements of the finite element model, and apply initial conditions and boundary conditions to the element nodes to obtain a temperature distribution of the part to be printed during the fused deposition modeling process and a corresponding thermal analysis result file. The mechanical model building unit 40 is configured to convert the thermal analysis unit into a structural unit, and define mechanical properties of a used material to obtain a mechanical model of the part to be printed. The displacement constraint unit 50 is configured to add displacement constraint to the bottom of the mechanical model, load the thermal analysis result file, and load the temperature of the node of the finite element model as a body load onto the corresponding node of the finite element model, so as to obtain stress distribution and deformation of the part to be printed in the fused deposition molding process. The invention also provides an electronic device, which comprises a processor 60 and a memory 70, wherein the memory 70 stores program instructions, and the processor 60 executes the program instructions to realize the prediction method of the temperature distribution and deformation of the part in the fused deposition manufacturing process. It should be noted that the Processor 60 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; or a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, or a discrete hardware component; the Memory 70 may include a Random Access Memory (RAM), and may further include a Non-volatile Memory (Non-volatile Memory), such as at least one disk Memory. The Memory 70 may also be an internal Memory of Random Access Memory (RAM) type, and the processor 60 and the Memory 70 may be integrated into one or more independent circuits or hardware, such as: application Specific Integrated Circuit (ASIC). It should be noted that the computer program in the memory 70 may be implemented in the form of software functional units and stored in a computer readable storage medium when the computer program is sold or used as a stand-alone product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, an electronic device, or a network device) to perform all or part of the steps of the method according to the embodiments of the present invention.
In summary, the method for predicting the temperature distribution and deformation of the part in the fused deposition manufacturing process can obtain the temperature field and the stress field of the part in the molding process by combining thermal analysis and mechanical analysis, can predict the warping deformation of the part, and provides technical support for optimizing printing parameters and reducing the defects of a printed product.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.
Claims (10)
1. A method for predicting the temperature distribution and deformation of a part in a fused deposition manufacturing process is characterized by comprising the following steps:
defining the type of a thermal analysis unit and thermophysical parameters of a material used when printing a part;
establishing a simulation model according to the shape and size of a part to be printed, and carrying out meshing on the simulation model according to the set type of a thermal analysis unit to obtain a finite element model comprising nodes and units;
activating the units of the finite element model, and applying initial conditions and boundary conditions on unit nodes to obtain the temperature distribution of the part to be printed in the fused deposition molding process and a corresponding thermal analysis result file;
converting the thermal analysis unit into a structural unit, and defining mechanical properties of a used material to obtain a mechanical model of the part to be printed;
and adding displacement constraint at the bottom of the mechanical model, loading the thermal analysis result file, and loading the temperature of the node of the finite element model as a body load onto the corresponding node of the finite element model to obtain the stress distribution and deformation of the part to be printed in the fused deposition molding process.
2. The method as claimed in claim 1, wherein the step of creating a simulation model according to the shape and size of the part to be printed, and the step of meshing the simulation model according to the set thermal analysis unit type comprises:
establishing a simulation model according to the shape and the size of the part to be printed;
obtaining the printing line width and the layer thickness by using a finite volume method and combining the layering thickness and the printing line width adopted by the part to be printed;
and taking the printing line width as the length and the width of the thermal analysis unit, taking the layer thickness as the height of the thermal analysis unit, and performing grid division on the simulation model according to the set type of the thermal analysis unit.
3. The method of claim 1, wherein the step of activating the elements of the finite element model comprises:
killing all the elements of the finite element model by using a life-death element technology, and activating the Nth element in the finite element model according to a printing path, wherein N represents any one element in the finite element model;
applying initial conditions and boundary conditions on the Nth unit of the finite element model, and solving the temperature of the temperature field according to a set load loading mode and load steps to delete the initial temperature load of the Nth unit of the finite element model;
activating the (N + 1) th unit of the finite element model, and solving the temperature of the temperature field to delete the initial temperature load of the (N + 1) th unit of the finite element model; until all elements of the finite-element model are activated.
4. The method as claimed in claim 3, wherein the step of obtaining the temperature distribution of the part to be printed and the corresponding thermal analysis result file from the initial conditions and the boundary conditions applied to the unit nodes comprises:
and solving the temperature field according to a set load loading mode and load steps, and applying initial conditions and boundary conditions on the activated nodes of each unit to obtain the temperature distribution of the part to be printed and a corresponding thermal analysis result file.
5. The method of claim 1, wherein the step of predicting the temperature distribution and deformation of the part during the fused deposition fabrication process comprises: the thermophysical parameters of the used material comprise density, specific heat capacity, heat conductivity coefficient and comprehensive heat exchange coefficient.
6. The method of claim 4, wherein the step of predicting the temperature distribution and deformation of the part during the fused deposition fabrication process comprises: the set load loading mode is a step mode load loading mode, and the step length of the load step is set according to the printing speed and the length of each unit.
7. The method of claim 4, wherein the step of predicting the temperature distribution and deformation of the part during the fused deposition fabrication process comprises: the initial conditions include an initial temperature, and the initial temperature includes one or more of a nozzle temperature, a forming chamber temperature, and a soleplate temperature.
8. The method of claim 4, wherein the step of predicting the temperature distribution and deformation of the part during the fused deposition fabrication process comprises: the boundary conditions include heat transfer between the molded material, heat transfer between the molded material and the air in the molding chamber by thermal convection and thermal radiation, and latent heat of phase change of the phase change during gradual cooling of the molded material from a molten state to a solid state.
9. A system for predicting temperature distribution and distortion of a part during a fused deposition manufacturing process, the system comprising:
the processor unit is used for defining the type of the thermal analysis unit and thermophysical parameters of materials used when the parts are printed;
the system comprises a simulation model establishing unit, a thermal analysis unit and a thermal analysis unit, wherein the simulation model establishing unit is used for establishing a simulation model according to the shape and the size of a part to be printed and carrying out grid division on the simulation model according to the set type of the thermal analysis unit so as to obtain a finite element model comprising nodes and units;
the activator unit is used for activating the units of the finite element model and applying initial conditions and boundary conditions on unit nodes to obtain the temperature distribution of the part to be printed in the fused deposition molding process and a corresponding thermal analysis result file;
the mechanical model establishing unit is used for converting the thermal analysis unit into a structural unit and defining the mechanical properties of the used materials so as to obtain a mechanical model of the part to be printed;
and the displacement constraint unit is used for adding displacement constraint at the bottom of the mechanical model, loading the thermal analysis result file, and loading the temperature of the node of the finite element model as a body load onto the corresponding node of the finite element model to obtain the stress distribution and deformation of the part to be printed in the fused deposition molding process.
10. An electronic device comprising a processor and a memory, the memory storing program instructions, characterized in that: the processor executes the program instructions to implement the method for predicting the temperature distribution and deformation of the part in the fused deposition manufacturing process according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010235780.3A CN111460709A (en) | 2020-03-30 | 2020-03-30 | Method for predicting temperature distribution and deformation of part in fused deposition manufacturing process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010235780.3A CN111460709A (en) | 2020-03-30 | 2020-03-30 | Method for predicting temperature distribution and deformation of part in fused deposition manufacturing process |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111460709A true CN111460709A (en) | 2020-07-28 |
Family
ID=71681516
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010235780.3A Withdrawn CN111460709A (en) | 2020-03-30 | 2020-03-30 | Method for predicting temperature distribution and deformation of part in fused deposition manufacturing process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111460709A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113128099A (en) * | 2021-05-08 | 2021-07-16 | 江苏师范大学 | Turning workpiece frequency prediction method |
| CN114536769A (en) * | 2022-02-15 | 2022-05-27 | 浙江工业大学 | Deformation regulation and control method of deployable antenna back frame structure based on 4D printing |
| CN114919181A (en) * | 2022-05-30 | 2022-08-19 | 北京航空航天大学 | Continuous fiber 3D printing process dynamic simulation and printed part buckling deformation prediction method |
-
2020
- 2020-03-30 CN CN202010235780.3A patent/CN111460709A/en not_active Withdrawn
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113128099A (en) * | 2021-05-08 | 2021-07-16 | 江苏师范大学 | Turning workpiece frequency prediction method |
| CN114536769A (en) * | 2022-02-15 | 2022-05-27 | 浙江工业大学 | Deformation regulation and control method of deployable antenna back frame structure based on 4D printing |
| CN114536769B (en) * | 2022-02-15 | 2024-04-02 | 浙江工业大学 | Deformation regulation and control method of deployable antenna back frame structure based on 4D printing |
| CN114919181A (en) * | 2022-05-30 | 2022-08-19 | 北京航空航天大学 | Continuous fiber 3D printing process dynamic simulation and printed part buckling deformation prediction method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liang et al. | Modified inherent strain method for efficient prediction of residual deformation in direct metal laser sintered components | |
| CN111460709A (en) | Method for predicting temperature distribution and deformation of part in fused deposition manufacturing process | |
| Zhou et al. | Topology optimization of thermal conductive support structures for laser additive manufacturing | |
| Li et al. | Efficient predictive model of part distortion and residual stress in selective laser melting | |
| Nagesha et al. | Review on characterization and impacts of the lattice structure in additive manufacturing | |
| Ramos et al. | New scanning strategy to reduce warpage in additive manufacturing | |
| Yavari et al. | Thermal modeling in metal additive manufacturing using graph theory | |
| You et al. | Design and additive manufacturing of thermal metamaterial with high thermal resistance and cooling capability | |
| Xu et al. | Residual stress constrained self-support topology optimization for metal additive manufacturing | |
| Van Belle et al. | Comparisons of numerical modelling of the Selective Laser Melting | |
| Baran et al. | Thermo-chemical modelling strategies for the pultrusion process | |
| JP7125764B2 (en) | Laminate-molded article analysis method, laminate-molded article analysis apparatus, laminate-molded article manufacturing method, and laminate-molded article manufacturing apparatus | |
| Cao et al. | A novel high-efficient finite element analysis method of powder bed fusion additive manufacturing | |
| Stavropoulos et al. | Warping in SLM additive manufacturing processes: estimation through thermo-mechanical analysis | |
| Özcan et al. | Residual stresses in metal deposition modeling: discretizations of higher order | |
| Lee et al. | Design optimization of thermally conductive support structure for laser powder-bed fusion process with part-scale thermal history | |
| Takezawa et al. | Optimally variable density lattice to reduce warping thermal distortion of laser powder bed fusion | |
| Ferreira et al. | A numerical and experimental study of fracture in RP stereolithography patterns and ceramic shells for investment casting | |
| Ernot et al. | Analytical and numerical predictions of the thermal performance of multi-layered lattice structures | |
| Subedi et al. | Towards the optimal design of support structures for laser powder bed fusion-based metal additive manufacturing via thermal equivalent static loads | |
| Carraturo et al. | Two-level method part-scale thermal analysis of laser powder bed fusion additive manufacturing | |
| Da et al. | Design of heterogeneous mesostructures for nonseparated scales and analysis of size effects | |
| Gan et al. | Prediction of residual deformation and stress of laser powder bed fusion manufactured Ti-6Al-4V lattice structures based on inherent strain method | |
| Daviran et al. | Differential quadrature method for thermal shock analysis of CNT reinforced metal-ceramic functionally graded disc | |
| US20110264425A1 (en) | Molten alloy solidification analyzing method and solidification analyzing program for performing the same |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200728 |