CN110472355B - 3D printing preview method based on multi-field coupling modeling and simulation solving - Google Patents
3D printing preview method based on multi-field coupling modeling and simulation solving Download PDFInfo
- Publication number
- CN110472355B CN110472355B CN201910771451.8A CN201910771451A CN110472355B CN 110472355 B CN110472355 B CN 110472355B CN 201910771451 A CN201910771451 A CN 201910771451A CN 110472355 B CN110472355 B CN 110472355B
- Authority
- CN
- China
- Prior art keywords
- field
- model
- printing
- structural part
- printing process
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 93
- 238000010146 3D printing Methods 0.000 title claims abstract description 68
- 238000010168 coupling process Methods 0.000 title claims abstract description 56
- 230000008878 coupling Effects 0.000 title claims abstract description 55
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 55
- 238000004088 simulation Methods 0.000 title claims abstract description 35
- 230000008569 process Effects 0.000 claims abstract description 56
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 230000007547 defect Effects 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 41
- 210000004027 cell Anatomy 0.000 claims description 24
- 238000009826 distribution Methods 0.000 claims description 23
- 238000010899 nucleation Methods 0.000 claims description 12
- 230000006911 nucleation Effects 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 9
- 238000009792 diffusion process Methods 0.000 claims description 6
- 230000001413 cellular effect Effects 0.000 claims description 4
- 239000000843 powder Substances 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000001276 controlling effect Effects 0.000 claims description 2
- 210000001787 dendrite Anatomy 0.000 claims description 2
- 238000007639 printing Methods 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims 1
- 238000012805 post-processing Methods 0.000 abstract 1
- 239000007790 solid phase Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 239000007791 liquid phase Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 239000013078 crystal Substances 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 239000012071 phase Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000007711 solidification Methods 0.000 description 4
- 230000008023 solidification Effects 0.000 description 4
- 238000004781 supercooling Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 238000003466 welding Methods 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000008018 melting Effects 0.000 description 2
- 238000002844 melting Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005184 irreversible process Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 229910001338 liquidmetal Inorganic materials 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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/20—Apparatus for additive manufacturing; Details thereof or accessories therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/12—Digital output to print unit, e.g. line printer, chain printer
- G06F3/1201—Dedicated interfaces to print systems
- G06F3/1202—Dedicated interfaces to print systems specifically adapted to achieve a particular effect
- G06F3/1203—Improving or facilitating administration, e.g. print management
- G06F3/1208—Improving or facilitating administration, e.g. print management resulting in improved quality of the output result, e.g. print layout, colours, workflows, print preview
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/12—Digital output to print unit, e.g. line printer, chain printer
- G06F3/1201—Dedicated interfaces to print systems
- G06F3/1223—Dedicated interfaces to print systems specifically adapted to use a particular technique
- G06F3/1237—Print job management
- G06F3/1253—Configuration of print job parameters, e.g. using UI at the client
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/12—Digital output to print unit, e.g. line printer, chain printer
- G06F3/1201—Dedicated interfaces to print systems
- G06F3/1278—Dedicated interfaces to print systems specifically adapted to adopt a particular infrastructure
- G06F3/1279—Controller construction, e.g. aspects of the interface hardware
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Computer Interaction (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Materials Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Mechanical Engineering (AREA)
- Optics & Photonics (AREA)
- Quality & Reliability (AREA)
Abstract
The invention discloses a 3D printing preview method based on multi-field coupling modeling and simulation solving, which aims at the 3D printing process of a structural part, utilizes a computer programming language and an interpolation algorithm to establish a multi-physical coupling model of a temperature field, a stress strain field, a flow field, a tissue field, a solute field and a defect field, solves and calculates the model by a finite element method, and carries out post-processing on the result to obtain the performance parameters of the 3D printing structural part under the mutual influence of different physical fields. The performance parameters of the 3D printing structure part in different physical fields are compared with the test results in real time and are judged and analyzed, so that the purposes of accurate preview and real-time regulation and control of the 3D printing process are achieved. According to the method, more reasonable 3D printing process parameters are worked out through judging and solving of the multi-field coupling model, the quality of the 3D printing structural part is improved, the 3D printing cost is effectively reduced, and the test efficiency is greatly improved.
Description
Technical Field
The invention relates to the technical field of 3D printing, in particular to a 3D printing preview method based on multi-field coupling modeling and simulation solving, and particularly relates to a digital preview and regulation and control technology of a 3D printing process.
Background
The 3D printing (also called Additive manufacturing, AM) technology is used as a rapid forming technology, and has the characteristics of integrated forming, high material utilization rate, high digitalization and the like, so that the method has a wide prospect in the field of manufacturing of precise and complex parts.
Whether a structural part formed using 3D printing techniques can be used in product manufacturing depends directly on the quality of the part. In order to obtain structural parts with required performance, parameter planning is required to be carried out before 3D printing so as to achieve the purpose of final preview, and real-time regulation and control are carried out in the 3D printing process so as to improve the production efficiency and the product quality. The traditional experimental method achieves the purpose of accumulating data by developing a large number of experiments, and has high cost, time and labor waste and difficulty in ensuring accuracy. The 3D printing process is an irreversible process, and the lack of a real-time regulation and control technology not only increases the time cost, but also has high requirements on the printing fault tolerance rate. With the improvement of computer simulation technology and related basic theory, there is a case of realizing preview of casting and welding process by simulating temperature field, stress strain field, microstructure field, etc. The 3D printing process is a process of rapid melting and rapid cooling, the existing time of a molten pool is very short, but the physical process of the molten pool is very complex, and the 3D printing process is simulated to realize 3D printing preview and regulation and control, which is a new challenge.
Disclosure of Invention
Aiming at the problems, the invention aims to design a 3D printing preview method based on multi-field coupling modeling and simulation solving so as to realize multi-field multi-scale visual preview of a 3D printing process. By means of carrying out multi-field coupling simulation of a temperature field, a flow field, a stress strain field, a microstructure, a defect field and the like in the 3D printing process of the structural part, the 3D printing technology and finite element analysis are comprehensively applied, the tissue, stress and performance prediction of the 3D printed structural part is achieved, and the purpose of 3D printing preview and regulation and control of the structural part is achieved.
In order to achieve the above object, the 3D print preview method based on multi-field coupling modeling and simulation solving in the present invention includes the following specific contents:
(1) constructing a multi-field coupling model of a 3D printing process of the structural part: establishing a three-dimensional geometric model and a grid model, defining material parameters, completing pretreatment of a solver, and completing establishment of a multi-field coupling model of a macroscopic temperature field, a stress strain field, a flow field, a microstructure field, a solute field and a defect field by adopting a computer programming language and finite element analysis software;
(2) establishing a solving algorithm of a multi-field coupling model: the macroscopic temperature field and the stress strain field adopt a finite element method, the macroscopic flow field adopts a finite volume method, the viscosity is realized by writing a User-defined function, the microstructure field, the solute field and the defect field adopt a specific Cellular Automata (CA) method, the conversion of a macroscopic physical field calculation result and microscopic physical field data is completed by an interpolation coupling algorithm, and the solving calculation of a multi-physical field coupling model is completed by applying an iterative algorithm;
(3) obtaining optimized 3D printing process parameters: and comparing and analyzing the calculation result of the multi-field coupling model with the preliminary test result in real time, judging that the multi-field coupling model returns to adjustment or continues to calculate according to the comparison result, obtaining the macroscopic and microscopic changes of the 3D printing process of the structural part when the performance parameters of the 3D printing structural part obtained by simulation completely accord with the test result, and outputting the 3D printing process parameters to achieve accurate preview and accurate guidance of the 3D printing process.
Preferably, in the step (1), the established three-dimensional geometric model close to reality is measured, analyzed and integrated based on partial experimental results, the geometric model is gridded by adopting a density-sparse combination gridding method, the material parameters comprise liquidus temperature, liquidus slope, solute distribution coefficient, liquid phase diffusion coefficient, solid phase diffusion coefficient, Gibbs-Thomson coefficient, initial concentration, cell size, time step, specific heat capacity, heat conductivity, young modulus, yield strength, density and the like, the thermosensitive material is defined as a function changing along with temperature, such as the specific heat capacity and the heat conductivity are defined as functions changing along with temperature, and the simulation software A, B, C, D and the programming software E are used for establishing the multi-field coupling model of the 3D printing structural part.
Preferably, the step (2) is to consider the calculation of the macroscopic flow field, approximate the flow variable to be solved by using a simple function in the finite element software C, substitute the approximate relation into the control equation of continuity to form a discrete equation system, and then solve the algebraic equation system. The CA method comprises a nucleation model, a dendrite growth model and a solute diffusion model, wherein the nucleation model is a uniform nucleation model and a non-uniform nucleation model which correspond to different algorithms and coexist. The coupled interpolation algorithm firstly discretizes a calculation region, divides the region into a limited number of grid nodes, each grid node is provided with a corresponding coordinate (i, j), then discretizes an equation by a finite difference method, approximately solves by using points around the coordinates, and changes continuous variables into discrete points; and finally, replacing the solution of the partial differential equation by an interpolation polynomial and the differential thereof to realize approximate solution. The iterative algorithm determines whether to continue the calculation by determining whether the state variables of all the cells are all "1".
Preferably, the macro-change in step (3) includes temperature field distribution, thermal cycle curve, residual stress distribution, residual deformation distribution, flow field distribution, etc. of the 3D printing process of the structural part, the micro-change includes microstructure distribution law, microstructure evolution, phase change, solute distribution evolution, defect distribution, formation law, etc., and 3D printing preview of the structural part can be realized according to the macro-and micro-changes and based on the 3D printed structural part performance parameters including porosity, crack distribution, deformation, microstructure size, etc. of the structural part. The judgment process of the comparison result on the multi-field coupling model starts from the judgment of the temperature field, the stress strain field and the flow field are in a mutual judgment relation, the calculation of the microscopic physical field can enter the next step only when the simulation results of the temperature field, the stress strain field and the flow field are matched with the preliminary test result, the microscopic tissue field, the solute field and the defect field adopt a sequential judgment sequence, when the final microscopic defect field is matched with the preliminary test result, 3D printing process parameters including laser power, powder feeding rate, laser scanning speed and layering thickness are output, the gradual realization of the judgment process of the multi-field coupling model through the comparison result can achieve the purpose of accurately regulating and controlling the 3D printing process.
The invention has the beneficial effects that:
aiming at the problem that the product quality and the production efficiency are difficult to meet the actual production requirements due to a plurality of factors such as low fault-tolerant rate, high time cost and the like in the current 3D printing process, the aims of accurate preview and accurate regulation and control of the 3D printing process of the structural part can be fulfilled by establishing a multi-physical-field coupling model and solving the model by adopting a specific algorithm. The model can well simulate the temperature field, can consider the influence of the temperature field change on the microscopic physical field, the macroscopic flow field and the stress strain field, establishes the coupling of the macroscopic physical field and the microscopic physical field, and realizes the multi-field coupling simulation of the powder melting and solidifying processes in the 3D printing process of the structural part. The simulation process can be compared with the preliminary experiment result in real time, so that model regulation and control and optimization of 3D printing process parameters are completed. The simulation result can accord with the macro and micro appearance of the 3D printing structure part in the laboratory.
Drawings
FIG. 1 is a flow diagram of an implementation of a three-dimensional finite element modeling and simulation method for a macroscopic temperature field and a stress-strain field of a 3D print preview model based on multi-field coupling solution;
FIG. 2 is a flow chart of an implementation of a three-dimensional finite element modeling and simulation method for a macroscopic flow field of a 3D print preview model based on multi-field coupling solution;
FIG. 3 is a flow diagram of an implementation of a microscopic simulation for a 3D print preview model based on multi-field coupling solution;
FIG. 4 is a flow diagram of a simulation algorithm implementation for a 3D print preview model based on multi-field coupling solution;
FIG. 5 is a flow diagram of an implementation of a 3D print preview model for multi-field coupling based solution.
Detailed Description
The following describes a 3D print preview method based on multi-field coupling modeling and simulation solving in detail with reference to the accompanying drawings.
The working flow of the method of the invention is shown in figures 1-3.
FIG. 1 is an implementation flow diagram of a three-dimensional finite element modeling and simulation method for macroscopic temperature and stress-strain fields of a 3D print preview model based on multi-field coupling solution.
and 2, establishing a grid model, and performing grid division on the geometric model of the 3D printing structure part in the finite element software B in a density transition grid division mode to ensure the efficiency and the accuracy of calculation.
Step 3-5 is to establish a 3D printing finite element model of the structural part in finite element software D finite element analysis software, which comprises the following steps:
(1) applying the material property. Defining the relation among the variables such as elastic modulus, yield strength, thermal expansion coefficient, thermal conductivity, specific heat capacity and temperature, defining Poisson's ratio and density of the material, and applying the defined material performance to the corresponding unit;
(2) a laser scan path is defined. The 3D printing process of the structural part is generally multilayer and multichannel laser scanning, so that a plurality of laser scanning paths need to be set, the paths are set by adopting a node method, the laser pointing method adopts a node method, and nodes are selected according to the relation between the previously set laser scanning paths and the laser irradiation direction;
(3) thermodynamic boundary conditions are imposed. The thermal boundary condition firstly selects a proper heat source model to characterize the action of the laser beam, and sets heat source parameters, wherein the heat source parameters comprise: welding energy, effective power coefficient, heat source width, heat source depth, heat source front half ellipsoid length, heat source rear half ellipsoid length, Gaussian heat source parameters: setting a welding speed after finishing setting heat source model parameters, selecting a welding path corresponding to the parameters, selecting a unit range which can be included by a heat source, applying convection boundary conditions of a workpiece and an external environment, setting the convection coefficient of the workpiece and the external environment to be 40, setting the temperature of the ambient environment to be 20 ℃, and selecting all surfaces capable of radiating to set radiating boundary conditions; the mechanical boundary condition setting is mainly to prevent the object from generating rigid displacement, and the constraint applying principle can ensure that the object does not generate rigid displacement, but can not add redundant rigid constraint. And (5) carrying out displacement constraint on the structural units in X, Y and Z directions to complete the setting of mechanical boundary conditions.
And 6, after the finite element model is established, calculating the temperature field and the stress field by adopting a finite element method.
FIG. 2 is a flow chart of an implementation of a three-dimensional finite element modeling and simulation method for a macroscopic flow field of a 3D print preview model based on multi-field coupling solution.
and step 12, outputting a numerical simulation result.
FIG. 3 is a flow diagram of an implementation of a microscopic simulation for a 3D print preview model based on multi-field coupling solution.
step 14-16 macroscopic microscopic temperature field coupling process. And (3) simulating the macroscopic temperature field and the microscopic structure by using two kinds of software, namely finite element software D and programming software E respectively, and adopting a weak coupling mode in the process of realizing the coupling of the macroscopic temperature field and the microscopic temperature field. In the weak coupling mode, the temperature of each unit cell is in a unique solidification path, the macroscopic time step is obtained by macroscopic node temperature interpolation, and the CA nucleation and growth model is circularly called to complete simulation.
FIG. 4 is a flow chart of a simulation algorithm implementation for a 3D print preview model based on multi-field coupling solution.
(1) calculating area discretization: the calculation region is divided into an M × N grid. Dividing the region into a finite number of mesh nodes, each mesh node having a corresponding coordinate (i, j);
(2) approximate substitution: discretizing an equation by a finite difference method, and according to a difference format, approximately solving by using points around coordinates, wherein continuous variables are changed into discrete points through the processing;
(3) and (3) approximation solving: this process can be seen as a process of replacing the solution of the partial differential equation with an interpolation polynomial and its differential, and the interpolation rule can be simply expressed as:
(1) dividing a grid: the calculation area is divided by adopting orthogonal grids, and the principle is mainly followed when the orthogonal grids are divided: 1) the grid adopts a square regular grid with the same size; 2) each grid is endowed with different characteristic variables and states, such as temperature, concentration, liquid phase, solid phase and the like; 3) the state change of the grid is determined according to the state of the neighbor grid; 4) the grid size is selected, and the calculation capability of the equipment and the instability caused by calculation are comprehensively considered;
(2) definition of thermophysical parameters: thermophysical parameters used include: density, liquidus temperature, solidus temperature, thermal conductivity, latent heat of solidification, specific heat capacity, and the like;
numerical value processing: the multi-field coupling numerical calculation model mainly comprises pretreatment, numerical calculation and post-treatment. 1) Preprocessing (initialization process of program): the thermal physical parameters, the temperature and the solute fraction are written into a program, one or more grain growth core points are arranged in a calculation area at the beginning of calculation, and the states of the rest grid points are respectively assigned to be solid-phase cells, liquid-phase cells or interface-layer cells according to the position of a growth core and a capture rule; 2) numerical value processing: and (4) calculating a control equation, and combining the initial condition and the interface judgment condition to iteratively solve the dispersed equation. After the initialization is completed, heat and mass transfer exists in each time step, and heat and mass transfer exists in the whole calculation area, so that heat and mass transfer calculation needs to be carried out on the grids in the whole area. According to the capture principle of the unit cell, when the unit cell meets the following two conditions, the unit cell becomes an interface layer unit cell, wherein whether at least one unit cell is a solid phase in the adjacent unit cell or not is judged, and whether the unit cell node is positioned in the interface layer or not is judged. In the multi-field coupling micro model, one or more grains which are arranged in advance are selected in the initial stage and are endowed with a preferential growth direction of 0-90 degrees. Then calculating the growth of the crystal grains, and using a new capture rule to make the interface continuously move towards the liquid phase and the crystal grains continuously grow up; 3) and (3) post-treatment: after the calculation is completed, a series of data files can be generated, and when the set conditions are met, the data stops being output. And (4) sorting and analyzing the simulation results, and performing certain processing on the results by using the visual software F to visually display the simulation results.
FIG. 5 is a flow diagram of an implementation of a 3D print preview model for multi-field coupling based solution.
And 24, calculating a temperature field, a stress strain field, a flow field, a solute field, a microstructure field and a microscopic defect field in sequence, wherein the steps 1 to 22 are a specific calculation process of the multi-physical field in the step 24. Firstly, calculating a temperature field to obtain the temperature field distribution characteristic of the structural part in the 3D printing process, comparing and judging the temperature field distribution characteristic with the physical property parameters of the actual structural part in the 3D printing process, if the temperature field distribution characteristic is consistent with the actual result, continuing the next step of calculation, if the temperature field distribution characteristic is not consistent with the actual result, redesigning the process parameters of the 3D printing structural part, and returning to the step 23; finishing temperature field judgment, performing stress-strain field calculation, comparing and judging a stress-strain characteristic result in the 3D printing process of the structural part with physical property parameters of the 3D printing process of the actual structural part, continuing the next step of calculation if the result is consistent with an actual result, redesigning the process parameters of the 3D printed structural part if the result is not consistent with the actual result, and returning to the step 23; and (5) finishing layer-by-layer calculation and judgment of the flow field, the solute field, the microstructure field and the microscopic defect field in sequence, and finally entering step 25 if all results are consistent with physical parameters of the actual structural part in the 3D printing process, outputting process parameters of the 3D printed structural part, and performing an actual 3D printing experiment on the structural part.
The foregoing is only a preferred embodiment of this invention and it should be noted that modifications can be made by those skilled in the art without departing from the principle of the invention and these modifications should also be considered as the protection scope of the invention.
Claims (3)
1. A3D printing preview method based on multi-field coupling modeling and simulation solving is characterized by comprising the following steps:
(1) constructing a multi-field coupling model of a 3D printing process of the structural part: establishing a three-dimensional geometric model and a grid model, defining material parameters, completing pretreatment of a solver, and completing establishment of a multi-field coupling model of a macroscopic temperature field, a stress strain field, a flow field, a microstructure field, a solute field and a defect field by adopting a computer programming language and finite element analysis software;
(2) establishing a solving algorithm of a multi-field coupling model: the macroscopic temperature field and the stress strain field adopt a finite element method, the macroscopic flow field adopts a finite volume method, the viscosity is realized by writing a User-defined function, the microstructure field, the solute field and the defect field adopt a specific Cellular Automata (CA) method, the conversion of a macroscopic physical field calculation result and microscopic physical field data is completed by an interpolation coupling algorithm, and the solving calculation of a multi-physical field coupling model is completed by applying an iterative algorithm;
(3) obtaining optimized 3D printing process parameters: comparing and analyzing the calculation result of the multi-field coupling model with the preliminary test result in real time, judging that the multi-field coupling model returns to adjustment or continues to calculate according to the comparison result, obtaining the macroscopic and microscopic changes of the 3D printing process of the structural part when the performance parameters of the 3D printing structural part obtained by simulation completely accord with the test result, and outputting the 3D printing process parameters to achieve accurate preview and accurate guidance of the 3D printing process;
the macroscopic change in the step (3) comprises temperature field distribution, a thermal cycle curve, residual stress distribution, residual deformation distribution and flow field distribution of the structural part in the 3D printing process, the microscopic change comprises a microstructure distribution rule, a microstructure evolution rule, phase change, solute distribution evolution, defect distribution and a formation rule, and 3D printing preview of the structural part can be realized according to the macroscopic change and the microscopic change and based on the performance parameters of the 3D printed structural part, including porosity, crack distribution, deformation and microstructure size of the structural part; the judgment process of the comparison result on the multi-field coupling model starts from the judgment of the temperature field, the stress strain field and the flow field are in a mutual judgment relation, the calculation of the microscopic physical field can enter the next step only when the simulation results of the temperature field, the stress strain field and the flow field are matched with the preliminary test result, the microscopic tissue field, the solute field and the defect field adopt a sequential judgment sequence, when the final microscopic defect field is matched with the preliminary test result, 3D printing process parameters including laser power, powder feeding rate, laser scanning speed and layering thickness are output, the gradual realization of the judgment process of the multi-field coupling model through the comparison result can achieve the purpose of accurately regulating and controlling the 3D printing process.
2. The 3D printing preview method based on multi-field coupling modeling and simulation solving of claim 1, wherein the step (1) is to integrate the established three-dimensional geometric model close to reality based on partial experimental result measurement and analysis, to perform grid division on the geometric model by adopting a density-sparse combination grid division method, to define the thermal sensitive material as a function varying with temperature, to establish the multi-field coupling model of the 3D printing structural part by using simulation software CATIA, Hypermesh, Fluent, MSC.
3. The 3D print preview method based on multi-field coupling modeling and simulation solution as claimed in claim 1, wherein the step (2) is to consider the calculation of the macroscopic flow field, approximate the flow variables to be solved in Fluent by using simple functions, substitute the approximate relationship into the control equation of continuity to form a discrete equation system, and then solve the algebraic equation system; the cellular automata method comprises a nucleation model, a dendrite growth model and a solute diffusion model, wherein the nucleation model is a uniform nucleation model and a non-uniform nucleation model which are corresponding to different algorithms and are coexisted; the interpolation coupling algorithm firstly discretizes a calculation region, divides the region into a limited number of grid nodes, each grid node is provided with a corresponding coordinate (i, j), then discretizes an equation by a finite difference method, approximately solves by using points around the coordinates, and changes continuous variables into discrete points; finally, an interpolation polynomial and the differential thereof are used for replacing the solution of a partial differential equation to realize approximate solution; the iterative algorithm determines whether to continue the calculation by determining whether the state variables of all the cells are all "1".
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910771451.8A CN110472355B (en) | 2019-08-20 | 2019-08-20 | 3D printing preview method based on multi-field coupling modeling and simulation solving |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910771451.8A CN110472355B (en) | 2019-08-20 | 2019-08-20 | 3D printing preview method based on multi-field coupling modeling and simulation solving |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110472355A CN110472355A (en) | 2019-11-19 |
CN110472355B true CN110472355B (en) | 2021-09-07 |
Family
ID=68513080
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910771451.8A Active CN110472355B (en) | 2019-08-20 | 2019-08-20 | 3D printing preview method based on multi-field coupling modeling and simulation solving |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110472355B (en) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111037143B (en) * | 2019-11-28 | 2022-09-06 | 新兴际华集团有限公司 | Method for regulating and controlling weld joint structure performance based on welding thermal cycle |
CN111168067B (en) * | 2020-01-22 | 2021-07-27 | 南京理工大学 | Pore prediction and control method based on laser directional energy deposition |
CN111203537B (en) * | 2020-04-22 | 2020-07-28 | 中国航发上海商用航空发动机制造有限责任公司 | Method for prefabricating fusion defect by controlling L MD process |
CN111203538B (en) * | 2020-04-22 | 2020-07-28 | 中国航发上海商用航空发动机制造有限责任公司 | Prefabricated crack defect, preparation method of built-in crack defect and prefabricated part |
CN111546636B (en) * | 2020-04-27 | 2021-06-29 | 中国水利水电科学研究院 | Simulation monitoring method for 3D printing of large structural body |
CN112182911B (en) * | 2020-10-26 | 2022-06-17 | 同济大学 | 3D printing filling design method based on force flow tube load path |
CN112517924B (en) * | 2020-11-11 | 2023-03-10 | 三峡大学 | Electric arc additive manufacturing method based on reverse compensation |
CN112364522A (en) * | 2020-11-20 | 2021-02-12 | 中国科学院长春光学精密机械与物理研究所 | Laser single-point curing modeling method applied to 3D printing lens |
CN112630046A (en) * | 2020-12-15 | 2021-04-09 | 中国科学院长春光学精密机械与物理研究所 | Method and system for measuring performance of high-temperature material |
CN112989626B (en) * | 2021-04-13 | 2022-12-23 | 清华大学 | Additive manufacturing organization simulation method, device, computer equipment and storage medium |
CN113158368B (en) * | 2021-04-16 | 2022-12-09 | 西安交通大学 | Method and equipment for predicting and controlling deformation of material-increasing and material-decreasing composite manufacturing part |
CN113118458B (en) * | 2021-04-20 | 2023-04-07 | 江西省科学院应用物理研究所 | Prediction method for tensile property of metal component formed by selective laser melting |
CN115310343A (en) * | 2021-05-07 | 2022-11-08 | 苏州奇流信息科技有限公司 | Sample database system, method for training and checking printing parameters and computer |
CN114872324B (en) * | 2022-04-15 | 2023-09-29 | 华中科技大学 | Laser additive manufacturing method based on multidimensional information coupling regulation performance |
CN115186557B (en) * | 2022-07-22 | 2023-09-26 | 南京航空航天大学 | Additive manufacturing component microstructure homogenization regulation and control method based on multi-physical field coupling solution heat accumulation effect |
CN116118189A (en) * | 2023-01-16 | 2023-05-16 | 石家庄铁道大学 | 3D printing technology-based rutting test block structure modulus targeting design method |
CN117648876B (en) * | 2024-01-30 | 2024-04-12 | 中南大学 | TPMS gradient hierarchical structure inverse design manufacturing method based on performance and BALANCE-CGAN |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021795A (en) * | 2016-06-03 | 2016-10-12 | 南昌航空大学 | Solidification temperature gradient controllable method for 3D printing of metal material |
WO2016187106A1 (en) * | 2015-05-19 | 2016-11-24 | President And Fellows Of Harvard College | Apparatus and method for high temperature 3d printing |
CN107368642A (en) * | 2017-07-13 | 2017-11-21 | 武汉大学 | The multiple dimensioned multiple physical field coupling simulation method of metal increasing material manufacturing |
WO2018052487A1 (en) * | 2016-09-19 | 2018-03-22 | Siemens Product Lifecycle Management Software Inc. | System and method for modeling characteristics of a melt pool that forms during an additive manufacturing process |
CN109261940A (en) * | 2018-09-28 | 2019-01-25 | 平顶山学院 | A kind of metal material increasing material manufacturing forming method and device |
CN109284524A (en) * | 2018-07-19 | 2019-01-29 | 西北工业大学 | A method of creation high-precision increasing material manufacturing finite element model |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108197384B (en) * | 2017-12-30 | 2021-11-02 | 大族激光科技产业集团股份有限公司 | Parameter optimization method of powder laser 3D printing forming process |
CN109858085B (en) * | 2018-12-26 | 2021-09-14 | 钢铁研究总院 | Austenitization determination method in heat treatment process of metal material |
-
2019
- 2019-08-20 CN CN201910771451.8A patent/CN110472355B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016187106A1 (en) * | 2015-05-19 | 2016-11-24 | President And Fellows Of Harvard College | Apparatus and method for high temperature 3d printing |
CN106021795A (en) * | 2016-06-03 | 2016-10-12 | 南昌航空大学 | Solidification temperature gradient controllable method for 3D printing of metal material |
WO2018052487A1 (en) * | 2016-09-19 | 2018-03-22 | Siemens Product Lifecycle Management Software Inc. | System and method for modeling characteristics of a melt pool that forms during an additive manufacturing process |
CN107368642A (en) * | 2017-07-13 | 2017-11-21 | 武汉大学 | The multiple dimensioned multiple physical field coupling simulation method of metal increasing material manufacturing |
CN109284524A (en) * | 2018-07-19 | 2019-01-29 | 西北工业大学 | A method of creation high-precision increasing material manufacturing finite element model |
CN109261940A (en) * | 2018-09-28 | 2019-01-25 | 平顶山学院 | A kind of metal material increasing material manufacturing forming method and device |
Non-Patent Citations (1)
Title |
---|
Invar合金激光熔化沉积过程微观组织建模与仿真研究;孟遥;《中国优秀硕士学位论文全文数据库(工程科技Ⅰ辑)》;20190215(第02期);B022-51 * |
Also Published As
Publication number | Publication date |
---|---|
CN110472355A (en) | 2019-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110472355B (en) | 3D printing preview method based on multi-field coupling modeling and simulation solving | |
Peng et al. | Fast prediction of thermal distortion in metal powder bed fusion additive manufacturing: Part 1, a thermal circuit network model | |
JP2017215957A (en) | Scalable finite element simulation of additive manufacturing | |
EP2612266B1 (en) | Interactive filling simulation on 3d injection molding models | |
CN109513924B (en) | Surface roughness control method in selective laser melting process | |
CN108062432A (en) | A kind of method for numerical simulation of selective laser fusion process | |
CN106373184B (en) | A kind of 3 D-printing model puts required amount of support Method of fast estimating | |
Park et al. | Design and simulation-based optimization of cooling channels for plastic injection mold | |
Kim et al. | Process-level modeling and simulation for HP's Multi Jet Fusion 3D printing technology | |
Neiva et al. | Numerical modelling of heat transfer and experimental validation in powder-bed fusion with the virtual domain approximation | |
Dalpadulo et al. | Integrated CAD platform approach for Design for Additive Manufacturing of high performance automotive components | |
US20220088881A1 (en) | System and method for simulation-assisted additive manufacturing | |
Moreira et al. | A multi-criteria h-adaptive finite-element framework for industrial part-scale thermal analysis in additive manufacturing processes | |
Korneev et al. | Fabricated shape estimation for additive manufacturing processes with uncertainty | |
Anand et al. | Additive manufacturing simulation tools in education | |
Chen | Understanding of the modeling method in additive manufacturing | |
Chen et al. | Data-driven prognostic model for temperature field in additive manufacturing based on the high-fidelity thermal-fluid flow simulation | |
Danylchenko | Comparative analysis of computer systems for casting processes simulation | |
US20200387652A1 (en) | Computational Modeling of Procedural Language Descriptors Prior to Manufacture | |
US11571740B2 (en) | Fabricated shape estimation for additive manufacturing processes | |
Garanger et al. | Foundations of intelligent additive manufacturing | |
US9087165B2 (en) | Automatic extremum detection on a surface mesh of a component | |
Mehmert | Residual stress analysis and geometrical tolerances in powder bed fusion and direct energy deposition processes | |
Patil | A novel numerical framework for simulation of multiscale spatio-temporally non-linear systems in additive manufacturing processes. | |
Bruggeman et al. | Residual Stress Generation in Additive Manufacturing of Complex Lattice Geometries |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |