CN110472355A - A kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation - Google Patents
A kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation Download PDFInfo
- Publication number
- CN110472355A CN110472355A CN201910771451.8A CN201910771451A CN110472355A CN 110472355 A CN110472355 A CN 110472355A CN 201910771451 A CN201910771451 A CN 201910771451A CN 110472355 A CN110472355 A CN 110472355A
- Authority
- CN
- China
- Prior art keywords
- field
- printing
- model
- constitutional detail
- scenarios method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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 kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation, for the 3D printing process of constitutional detail, utilize computer programming language and interpolation algorithm, establish temperature field, stress-strain field, flow field, the coupling of multiple physics model for organizing field, solute field and defect field, solution calculating is carried out to model with FInite Element, and is post-processed to obtain the 3D printing constitutional detail performance parameter under influencing each other between different physical field to result.By by the 3D printing constitutional detail performance parameter under different physical field and test result real time contrast and carrying out decision analysis, to achieve the purpose that carry out accurate preview and real-time monitoring to 3D printing process.The invention makes more reasonable 3D printing technological parameter by the judgement solution to multi- scenarios method model, improves the quality of 3D printing constitutional detail, the cost of 3D printing is effectively reduced and greatly improves test efficiency.
Description
Technical field
The present invention relates to 3D printing technique fields, and in particular to it is a kind of based on multi- scenarios method modeling and simulation solve 3D beat
Print method for previewing, in particular to the digitlization preview of 3D printing process and control technique.
Background technique
3D printing (being also increasing material manufacturing, Additive manufacturing, AM) technology is as a kind of rapid shaping skill
Art has due to having the characteristics that integrated molding, stock utilization height, high digitlization in the manufacturing field of precision workpiece
Bright prospects.
The performance quality that product manufacturing directly depends on part can be used for using the constitutional detail that 3D printing technique shapes.
In order to obtain the constitutional detail of required performance, it is necessary to carry out parametric programming before 3D printing to achieve the purpose that final preview, In
Real-time monitoring is carried out during 3D printing to improve production efficiency and product quality.Traditional experimental method is a large amount of real by carrying out
It tests, to achieve the purpose that accumulate data, cost is quite high, time-consuming and laborious and accuracy is difficult to ensure.3D printing process is one
Irreversible procedure, lack real-time monitoring technology not only increase time cost also to printing serious forgiveness have high requirements.With calculating
Theoretical perfect of machine emulation technology and relevant rudimentary, has at present by moulds such as temperature field, stress-strain field, microstructure fields
Intend the case to realize casting and welding process preview.3D printing process is the process that a rapid melting is quickly cooled down again, is melted
Pond is extremely of short duration there are the time, but its physical process is sufficiently complex, the simulation to 3D printing process, realizes 3D printing preview
And regulation, and be a new challenge.
Summary of the invention
In view of the above problems, the present invention is directed to design it is a kind of based on multi- scenarios method modeling and simulation solve 3D beat
Method for previewing is printed, with more that realize 3D printing process multiple dimensioned intuitive previews.By to constitutional detail 3D printing process spread
The multi- scenarios methods such as temperature field, flow field, stress-strain field, microstructure and defect field simulation, by 3D printing technique and finite element fraction
Analysis carries out integrated application, realizes tissue, stress and the performance prediction of 3D printing constitutional detail, reaches constitutional detail 3D printing preview
With the purpose of regulation.
In order to achieve the above object, a kind of 3D printing solved based on multi- scenarios method modeling and simulation referred in the present invention
Method for previewing, particular content are as follows:
(1) it constructs the multi- scenarios method model of constitutional detail 3D printing process: establishing 3-D geometric model and grid model, it is fixed
Adopted material parameter completes solver pre-treatment, using computer programming language and finite element analysis software complete macro temperature field,
The foundation of stress-strain field, flow field, microstructure field, solute field and the multi- scenarios method model of defect field；
(2) establish the derivation algorithm of multi- scenarios method model: macro temperature field and stress-strain field are using FInite Element, macroscopic view
Flow field uses finite volume method, wherein viscosity is realized by writing User-defined function function, microstructure field,
Solute field and defect field use specific cellular automata (CA) method, interpolation coupling algorithm complete Macroscopic physical field computation result with
The conversion of microphysics field data, and calculated with the solution that iterative algorithm completes multiple physical field coupling model；
(3) the 3D printing technological parameter of optimization is obtained: the calculated result of multi- scenarios method model and Preliminary experiment results is real
When compare and analyze, determine that multi- scenarios method model returns according to comparing result and adjust or continue to calculate, when simulating resulting 3D
When print structure part performance parameter and test result comply fully with, the both macro and micro for obtaining constitutional detail 3D printing process becomes
Change, and export 3D printing technological parameter, reaches the accurate preview and accurate guidance to 3D printing process.
Preferably, step (1) is that the three-dimensional for the closing to reality established based on part of test results measurement analytical integration is several
What model carries out grid dividing to geometrical model using the Meshing Method that density combines, and the material parameter includes liquid phase
Line temperature, liquidous slopes, solute distribution coefficient, Liquid Diffusion Coefficient, solid phase diffusion welding, Gibbs-Thomson coefficient, just
Beginning concentration, cellular size, time step, specific heat capacity, thermal conductivity, Young's modulus, yield strength, density etc., thermo-sensitive material is defined
It is defined as the function varied with temperature at the function varied with temperature, such as specific heat capacity, thermal conductivity, with simulation softward A, B, C, D
The multi- scenarios method model of 3D printing constitutional detail is established with programming software E.
Preferably, step (2) is to consider the calculating in macroscopical flow field, in finite element software C using simple function approximation to
The flow variables asked substitute into the approximation relation in successional governing equation, form discrete equation group, solve algebra side later
Journey group.The CA method includes the homogeneous nucleation according to corresponding to algorithms of different and heterogeneous nucleation and the Nucleation Model deposited, branch
Crystals growth model and solutes accumulation model.The coupling interpolation algorithm is first by zoning discretization, region division at having
Limit mesh node, each grid has a corresponding coordinate (i, j), then by finite difference calculus to equation carry out from
It dissipates, using the point approximate solution around coordinate, continuous variable is become into discrete point；Finally with an interpolation polynomial and its
Differential replaces the solutions of partial differential equation, and realization approaches solution.The iterative algorithm is the shape by judging whether all cellulars
State variable is all " 1 " to determine whether continuing to calculate.
Preferably, the variation of macroscopic view described in step (3) includes the thermo parameters method of constitutional detail 3D printing process, thermal cycle
Curve, residual stress distribution, residual deformation distribution, Flow Field Distribution etc., the micro-variations include the microstructure regularity of distribution, micro-
Microstructure evolution, phase transformation, solute Distribution differentiation, defect distribution and formation rule etc. are seen, changes according to both macro and micro and is based on 3D
Print structure part performance parameter includes the porosity, crack distribution, deflection, microstructure size of structural member etc., Ke Yishi
Now to constitutional detail 3D printing preview.Comparing result to the decision process of multi- scenarios method model since determining temperature field, temperature
It is the relationship determined mutually between field and stress-strain field, flow field, when the simulation knot in temperature field, stress-strain field and flow field three
Fruit just can enter the calculating of next step microphysics field, microstructure field, solute field when completing and coincide with Preliminary experiment results
3D printing is exported when final microdefect field and identical Preliminary experiment results completion using successively judgement sequence with defect field
Technological parameter, including laser power, powder feeding rate, laser scanning speed, lift height, by comparing result to multi- scenarios method model
Decision process be done step-by-step, can achieve the purpose to 3D printing process accuracy controlling.
Beneficial effects of the present invention:
For product quality and production caused by many factors such as current 3D printing process fault tolerant rate is low, time cost is high
Efficiency is difficult to the problem of meeting actual production requirement, by establishing a multiple physical field coupling model and using special algorithm to it
It solves, can be realized the accurate preview of constitutional detail 3D printing process and the purpose accurately regulated and controled.Model can be well to temperature
Field is simulated, and can consider influence of the change of temperature field to microphysics field, macroscopical flow field and stress-strain field, it is established that macro
See physical field and microphysics field coupling, realize constitutional detail 3D printing during powder fusing, process of setting more couplings
Molding is quasi-.Simulation process can compare to complete model regulatory and 3D printing technological parameter in real time with preliminary experimental results
Optimization.Analog result can meet the both macro and micro pattern of 3D printing constitutional detail in laboratory.
Detailed description of the invention
Fig. 1 be 3D printing preview model for being solved based on multi- scenarios method macro temperature field and stress-strain field three
Tie up the implementation flow chart of finite element modeling and analogy method；
Fig. 2 be macroscopical flow field of 3D printing preview model for being solved based on multi- scenarios method Three-dimensional finite element modeling and
The implementation flow chart of analogy method；
Fig. 3 is the implementation flow chart of the microcosmic Simulation of the 3D printing preview model for being solved based on multi- scenarios method；
Fig. 4 is the simulation algorithm implementation flow chart of the 3D printing preview model for being solved based on multi- scenarios method；
Fig. 5 is the implementation flow chart of the 3D printing preview model for being solved based on multi- scenarios method.
Specific embodiment
The 3D printing method for previewing solved based on multi- scenarios method modeling and simulation a kind of to the present invention is made with reference to the accompanying drawing
It illustrates.
The method of the present invention workflow is as shown in Figure 1-Figure 3.
Fig. 1 be 3D printing preview model for being solved based on multi- scenarios method macro temperature field and stress-strain field three
Tie up the implementation flow chart of finite element modeling and analogy method.
Step 1 establishes 3D printing constitutional detail geometrical model, and specific steps include the reality according to 3D printing constitutional detail
Structure size carries out the building of 3D geometrical model in finite element modeling software A or L；
Step 2 establishes grid model, and the grid dividing mode of density transition is used to 3D printing constitutional detail geometrical model
Grid dividing is carried out, in finite element software B to guarantee the efficiency and accuracy that calculate.
Step 3-5 is that constitutional detail 3D printing finite element model is established in finite element software D finite element analysis software, is wrapped
It includes:
(1) apply material properties.Define elasticity modulus, yield strength, thermal expansion coefficient, thermal coefficient and specific heat capacity and temperature
The relationship of the variables such as degree, the material property of definition is applied on corresponding unit by the Poisson's ratio and density of definition material；
(2) laser beam scan path is defined.Constitutional detail 3D printing process is generally multilayer multiple tracks laser scanning, thus needs
Multi-stripe laser scan path is set, the setting in path uses nodal method, and the pointing method of laser uses nodal method, sets according to before
Relationship between fixed laser beam scan path and laser irradiation direction selects node；
(3) apply thermodynamic boundaries condition.Thermal boundary condition select first suitable heat source model to laser beam act on into
Row characterization, and Heat-Source Parameters are set, double-ellipsoid heat source parameter: welding energy, effective power coefficient, heat source width, heat
Semielliptical length after semielliptical length, heat source before Depth, heat source, Gauss Heat-Source Parameters: plane heat source power, plane heat source are effectively made
With radius, body heat source power, body heat source useful effect radius, body heat source effective interaction depth, heat source effective absorption coefficient, complete
Speed of welding is set after heat source model parameter setting, select the above parameter for path of welding, select heat source can
Including unit range, apply workpiece and external environment convection boundary condition, be arranged workpiece and external environment convection coefficient
It is 40, the temperature that ambient enviroment is arranged is 20 DEG C, selects all surface set heat dissipation boundary conditions that can be radiated；Mechanical boundary
Condition setting mainly prevents object from rigid displacement occurs, although the application principle constrained can guarantee that rigidity does not occur for object
Displacement, but extra rigid constraint can not be added.Displacement constraint is carried out to structure body unit in the x, y, and z directions, completes power
Learn the setting of boundary condition.
Step 6 completes temperature field and stress field calculation after the completion of finite element model is established, using finite element method.
Whether step 7 is completed by judging laser beam scan path to determine whether continuing to calculate, if necessary to calculating, entrance
Step 5 applies heat source model again, otherwise, exports analog result；
Fig. 2 be macroscopical flow field of 3D printing preview model for being solved based on multi- scenarios method Three-dimensional finite element modeling and
The implementation flow chart of analogy method.
Step 8 is solver setting.Before step 8, step 1-2 operation is completed, established grid model is read in
Finite element software C simultaneously checks grid.Using the flow variables that simple function approximation is to be asked in finite element software C, by the approximation
Relationship substitutes into successional governing equation, forms discrete equation group, later Solving Algebraic Equation group.
Step 9 is definition material attribute, comprising: density, specific heat capacity and the coefficient of heat conduction are equal to the related parameter of temperature, liquid
Liquidus temperature and solidus temperature, viscosity are realized by writing User-defined function function；
Step 10 is to define boundary condition, and the speed entrance side of the speed and scalar for defining flowing entrance boundary is arranged
Boundary's condition, comprising: define inflow velocity, setting temperature (when solving energy equation, needs the speed entrance boundary in temperature field
Set flowing static temperature), definition outflow normal pressure.The static pressure being arranged for defining flowing outlet (further includes it in reflux
His scalar) pressure export boundary condition, to need to input parameter include: static pressure, reflux condition to setting pressure export boundary condition
Part, total temperature, that is, stagnation temperature (be used for energy balane), turbulent parameters (being calculated for turbulent flow), the volume fraction of second level phase (for
The calculating of multiphase flow), discrete phase boundary condition (for discrete phase calculate)；
Step 11 solves the flow velocity of Equations of Turbulence each point and the Flow Field Distribution of entirety by iterative calculation；
Step 12 output numerical value analog result.
Fig. 3 is the implementation flow chart of the microcosmic Simulation of the 3D printing preview model for being solved based on multi- scenarios method.
The material property parameter for the 3D printing constitutional detail used in step 13 input temp, solute score and simulation, object
Property parameter includes liquidus temperature, liquidous slopes, solute distribution coefficient, Liquid Diffusion Coefficient, solid phase diffusion welding, Gibbs-
Thomson coefficient, initial concentration, cellular size, time step etc.；
Step 14-16 macro and micro temperature field coupling process.It is soft that finite element is respectively adopted in macro temperature field and microstructure
Two kinds of softwares of part D and programming software E are simulated, and use weak coupling mode in realizing macro microcosmic temperature field coupling.In weak coupling
In syntype, each cellular temperature is to obtain walking in macroscopic time by macroscopical node temperature interpolation under unique solidification path
In length, recursive call CA model of nucleation and growth completes simulation.
Step 17 is to establish forming core growth model, and model assigns each following information of cellular: temperature, solute concentration, crystal grain
Color variance, crystal grain are grown up orientation and crystalline form variable.Crystal grain color variance is used to indicate different solid-phase grains, each shape
The solid-phase grain of core has random orientation of growing up.Cellular state variation is indicated using continuous variable." 1 " represents solid phase, " 0 " generation
Table liquid phase represents solid/liquid interfaces between " 0 " and " 1 ", and the cellular in interfacial state can occur forming core and grow up.
Step 18 is the calculating process of forming core growth model, first determines whether degree of supercooling (Δ T) is greater than critical undercooling limits
(Δ Tn), when degree of supercooling (Δ T) is greater than critical undercooling limits (Δ Tn), generation solidifies forming core in liquid metal；Once core shape
Cheng Hou, nucleus continue to grow up and form crystal grain, therefore after forming core, and whether model, which grows up to nucleus, grow up progress
Judgement, when returning to program is to grow up signal, program is updated the state variable of the cellular, and the state of the cellular becomes
Amount becomes " 1 " from " 0 "；
Step 19 is by judging whether that the state variable of all cellulars is all " 1 " to determine whether continuing to calculate, by not
The disconnected return to cellular state variable, judges whether that the state variable of all cellulars is all " 1 ", when all cellulars become solid phase
When, process of setting is completed, and analog result is exported；If there are also cellulars to be in liquid phase, progress is secondary to be judged whether to reach time restriction,
If reaching time restriction, directly output analog result enters the simulation of subsequent time if not having, and has recycled until all
At output analog result；
Fig. 4 is the simulation algorithm implementation flow chart of the 3D printing preview model for being solved based on multi- scenarios method.
Step 20 is macro temperature field, stress-strain field and the calculating in flow field, and step 1-12 is the specific implementation of step 20
Process；
Step 21 is the calculating in microcosmic temperature field and flow field, this process is real by carrying out interpolation algorithm to macro temperature field
It is existing.Make (usually orthogonal equidistant) mesh generation of rule firstly the need of the region to calculating, basic step includes:
(1) zoning zoning discretization: is divided into the grid of M × N.Region division at limited grid knot
Point, each grid have a corresponding coordinate (i, j)；
(2) approximate substitution: it is discrete to equation progress by finite difference calculus, according to the format of difference, coordinate can be used
The point approximate solution of surrounding, by such processing, continuous variable becomes discrete point；
(3) approach solution: this process can be regarded as replacing partial differential side with an interpolation polynomial and its differential
The process of the solution of journey, interpolation rule can be represented simply as:
Step 22 is that the multi- scenarios method of microcosmic temperature field, microcosmic flow field, microstructure field, defect field and solute field calculates,
It is based on first three modeling to be calculated as a result, carrying out microcosmic coupling, it is necessary first to establish multi- scenarios method model.Include:
(1) division of grid: zoning is divided using orthogonal grid, mainly follows principle when dividing orthogonal grid: 1)
Grid uses the square regular grid of identical size；2) each grid can assign different characteristic variable and state, such as temperature
Degree, concentration, liquid phase and solid phase etc.；3) state change of grid is determined according to the state of neighbours' grid；4) selection of size of mesh opening,
Comprehensively considered from the computing capability of equipment and in terms of unstability two caused by calculating；
(2) definition of thermal physical property parameter: using to thermal physical property parameter include: density, liquidus temperature, solidus temperature
Degree, thermal coefficient, latent heat of solidification, specific heat capacity etc.；
Numerical value processing: multi- scenarios method mathematical calculation model mainly includes pre-treatment, numerical value calculates and post-processing.1) pre-treatment
(initialization procedure of program): it needs to calculate in thermophysical parameter, temperature and solute score write-in program and calculating when starting
The core point of one or more grain growths is set in region, according to the position of growth core and capture rule, respectively to remaining
Grid dotted state carries out assignment, is solid phase cellular, liquid phase cellular or boundary layer cellular；2) numerical value is handled: to the meter of governing equation
It calculates, and primary condition and interface Rule of judgment is combined to be iterated solution to the equation after discrete.After initialization is completed, every
In one time step, there is the transmission of heat and matter to need since heat, matter transmission all exist in entire zoning
Heat transfer and the calculating of mass transfer are all carried out to the grid in whole region.According to the capture principle of cellular, below cellular meets
When two conditions, become as boundary layer cellular, first is that in its neighbours' cellular, if having at least one cellular is solid phase, two
It is whether the cellular node is located in boundary layer.In multi- scenarios method micromodel, the incipient stage selects one or several shift to an earlier date
The crystal grain set, and it is endowed one 0 ° to 90 ° of selecting excellence evaluation.Then the growth of crystal grain is calculated, is used
New capture rule elapses interface constantly into liquid phase, and crystal grain is constantly grown up；3) post-process: calculating can after completing
Volume of data file is generated, when meeting the condition of setting, data just stop output.Analog result is arranged and is analyzed, is led to
It crosses and carries out certain processing to result using visual software F, intuitively show analog result.
Fig. 5 is the implementation flow chart of the 3D printing preview model for being solved based on multi- scenarios method.
The technological parameter for the 3D printing constitutional detail that step 23 input has designed, technological parameter include: that layer is high, printing is fast
Degree, print temperature, scan power and sweep span etc..
Step 24 successively carries out temperature field, stress-strain field, flow field, solute field, microstructure field and microdefect field meter
It calculates, step 1-22 is the specific calculating process of multiple physical field in step 24.Firstly, calculating temperature field, structure zero is obtained
Thermo parameters method characteristic during part 3D printing, by by thermo parameters method characteristic and practical structures part 3D printing process
Physical parameter compares judgement, if being consistent with actual result, continues to calculate in next step, if not being inconsistent together in actual result,
The technological parameter of 3D printing constitutional detail is redesigned, and return step 23；It completes temperature field to determine, carries out stress-strain field
It calculates, by the physical parameter of stress-strain characteristics result and border constitutional detail 3D printing process during constitutional detail 3D printing
Judgement is compared, if being consistent with actual result, continues to calculate in next step, if not being inconsistent together in actual result, to 3D printing
The technological parameter of constitutional detail redesigns, and return step 23；With the above operation, be successively performed flow field, flow field, solute field,
The calculating layer by layer of microstructure field and microdefect field and determine, if final all results with practical structures part 3D printing
The physical parameter of journey is consistent, then enters step 25, exports the technological parameter of 3D printing constitutional detail, carries out the reality of constitutional detail
3D printing experiment.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
It for member, can also make several improvements without departing from the principle of the present invention, these improvement also should be regarded as of the invention
Protection scope.
Claims (4)
1. a kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation, which is characterized in that comprise the following steps that
(1) it constructs the multi- scenarios method model of constitutional detail 3D printing process: establishing 3-D geometric model and grid model, define material
Expect parameter, complete solver pre-treatment, macro temperature field, stress are completed using computer programming language and finite element analysis software
The foundation of strain field, flow field, microstructure field, solute field and the multi- scenarios method model of defect field；
(2) establish the derivation algorithm of multi- scenarios method model: macro temperature field and stress-strain field are using FInite Element, macroscopical flow field
Using finite volume method, wherein viscosity is realized by writing User-defined function function, microstructure field, solute
And defect field use specific cellular automata (CA) method, interpolation coupling algorithm complete Macroscopic physical field computation result with it is microcosmic
The conversion of physical field data, and calculated with the solution that iterative algorithm completes multiple physical field coupling model；
(3) obtain the 3D printing technological parameter of optimization: by the calculated result of multi- scenarios method model and Preliminary experiment results in real time into
Row comparative analysis determines that multi- scenarios method model returns according to comparing result and adjusts or continue to calculate, when the resulting 3D printing of simulation
When constitutional detail performance parameter and test result comply fully with, the both macro and micro variation of constitutional detail 3D printing process is obtained,
And 3D printing technological parameter is exported, reach accurate preview and accurate guidance to 3D printing process.
2. a kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation according to claim 1, special
Sign is that above-mentioned steps (1) are the three-dimensional geometry moulds for the closing to reality established based on part of test results measurement analytical integration
Type carries out grid dividing to geometrical model using the Meshing Method that density combines, and the material parameter includes liquidus curve temperature
It is degree, liquidous slopes, solute distribution coefficient, Liquid Diffusion Coefficient, solid phase diffusion welding, Gibbs-Thomson coefficient, initial dense
Degree, cellular size, time step, specific heat capacity, thermal conductivity, Young's modulus, yield strength, density etc., by thermo-sensitive material be defined as with
The function of temperature change, if specific heat capacity, thermal conductivity are defined as the function varied with temperature, with simulation softward A, B, C, D and volume
Journey software E establishes the multi- scenarios method model of 3D printing constitutional detail.
3. a kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation according to claim 1, special
Sign is that above-mentioned steps (2) are to consider the calculating in macroscopical flow field, to be asked using simple function approximation in finite element software C
Flow variables substitute into the approximation relation in successional governing equation, form discrete equation group, later Solving Algebraic Equation
Group.The CA method includes the homogeneous nucleation according to corresponding to algorithms of different and heterogeneous nucleation and the Nucleation Model deposited, dendrite
Growth model and solutes accumulation model.The coupling interpolation algorithm is first by zoning discretization, region division at limited
A mesh node, each grid have a corresponding coordinate (i, j), then discrete to equation progress by finite difference calculus,
Using the point approximate solution around coordinate, continuous variable is become into discrete point；Finally with an interpolation polynomial and its micro-
Divide the solution to replace partial differential equation, realization approaches solution.The iterative algorithm is the state by judging whether all cellulars
Variable is all " 1 " to determine whether continuing to calculate.
4. a kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation according to claim 1, special
Sign is that the variation of macroscopic view described in above-mentioned steps (3) includes the thermo parameters method of constitutional detail 3D printing process, thermal cycle song
Line, residual stress distribution, residual deformation distribution, Flow Field Distribution etc., the micro-variations include the microstructure regularity of distribution, microcosmic
Microstructure evolution, phase transformation, solute Distribution differentiation, defect distribution and formation rule etc., are changed according to both macro and micro and are simultaneously beaten based on 3D
Printing constitutional detail performance parameter includes the porosity, crack distribution, deflection, microstructure size of structural member etc., be may be implemented
To constitutional detail 3D printing preview.Comparing result to the decision process of multi- scenarios method model since determining temperature field, temperature field
It is the relationship determined mutually between stress-strain field, flow field, when temperature field, the analog result of stress-strain field and flow field three
Just can enter the calculating of next step microphysics field when completing and coincide with Preliminary experiment results, microstructure field, solute field and
Defect field exports 3D printing work when final microdefect field and identical Preliminary experiment results completion using successively judgement sequence
Skill parameter, including laser power, powder feeding rate, laser scanning speed, lift height, by comparing result to multi- scenarios method model
Decision process is done step-by-step, and can achieve the purpose to 3D printing process accuracy controlling.
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 true CN110472355A (en) | 2019-11-19 |
| CN110472355B 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) |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111037143A (en) * | 2019-11-28 | 2020-04-21 | 新兴际华集团有限公司 | Method for regulating and controlling weld joint structure performance based on welding thermal cycle |
| CN111168067A (en) * | 2020-01-22 | 2020-05-19 | 南京理工大学 | Pore prediction and control method based on laser directional energy deposition |
| CN111203538A (en) * | 2020-04-22 | 2020-05-29 | 中国航发上海商用航空发动机制造有限责任公司 | Prefabricated crack defect, preparation method of built-in crack defect and prefabricated part |
| CN111203537A (en) * | 2020-04-22 | 2020-05-29 | 中国航发上海商用航空发动机制造有限责任公司 | Method for controlling defect of poor fusion of LMD process prefabrication |
| CN111546636A (en) * | 2020-04-27 | 2020-08-18 | 中国水利水电科学研究院 | Simulation monitoring method for 3D printing of large structural body |
| CN112182911A (en) * | 2020-10-26 | 2021-01-05 | 同济大学 | 3D printing filling design method based on force flow tube load path |
| CN112364522A (en) * | 2020-11-20 | 2021-02-12 | 中国科学院长春光学精密机械与物理研究所 | Laser single-point curing modeling method applied to 3D printing lens |
| CN112517924A (en) * | 2020-11-11 | 2021-03-19 | 三峡大学 | Electric arc additive manufacturing method based on reverse compensation |
| CN112630046A (en) * | 2020-12-15 | 2021-04-09 | 中国科学院长春光学精密机械与物理研究所 | Method and system for measuring performance of high-temperature material |
| CN112989626A (en) * | 2021-04-13 | 2021-06-18 | 清华大学 | Additive manufacturing organization simulation method, device, computer equipment and storage medium |
| CN113118458A (en) * | 2021-04-20 | 2021-07-16 | 江西省科学院应用物理研究所 | Prediction method for tensile property of metal component formed by selective laser melting |
| CN113158368A (en) * | 2021-04-16 | 2021-07-23 | 西安交通大学 | Prediction and control method and device for deformation of material-increasing and material-decreasing composite manufacturing part |
| CN114872324A (en) * | 2022-04-15 | 2022-08-09 | 华中科技大学 | Laser additive manufacturing method based on multidimensional information coupling regulation and control performance |
| CN115186557A (en) * | 2022-07-22 | 2022-10-14 | 南京航空航天大学 | Additive manufacturing component microstructure homogenization regulation and control method based on multi-physical-field coupling solution heat accumulation effect |
| WO2022233206A1 (en) * | 2021-05-07 | 2022-11-10 | 苏州奇流信息科技有限公司 | Sample database system, method for training and checking printing parameter, and computer |
| CN116118189A (en) * | 2023-01-16 | 2023-05-16 | 石家庄铁道大学 | 3D printing technology-based rutting test block structure modulus targeting design method |
| CN117648876A (en) * | 2024-01-30 | 2024-03-05 | 中南大学 | TPMS gradient hierarchical structure inverse design manufacturing method and device based on performance and BALANCE-CGAN, storage medium and electronic device |
Citations (8)
| 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 |
| CN108197384A (en) * | 2017-12-30 | 2018-06-22 | 大族激光科技产业集团股份有限公司 | The parameter optimization method of powder laser 3D printing moulding 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 |
| CN109858085A (en) * | 2018-12-26 | 2019-06-07 | 钢铁研究总院 | A kind of austenitizing measuring method in metal material heat treatment process |
-
2019
- 2019-08-20 CN CN201910771451.8A patent/CN110472355B/en active Active
Patent Citations (8)
| 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 |
| CN108197384A (en) * | 2017-12-30 | 2018-06-22 | 大族激光科技产业集团股份有限公司 | The parameter optimization method of powder laser 3D printing moulding process |
| 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 |
| CN109858085A (en) * | 2018-12-26 | 2019-06-07 | 钢铁研究总院 | A kind of austenitizing measuring method in metal material heat treatment process |
Non-Patent Citations (4)
| Title |
|---|
| CHANG LI等: "Numerical simulation and experimental study of cladding Fe60 on an ASTM 1045 substrate by laser cladding", 《SURFACE & COATINGS TECHNOLOGY》 * |
| LEI WANG等: "Simulation of dendrite growth in the laser welding pool of aluminumalloy 2024 under transient conditions", 《JOURNAL OF MATERIALS PROCESSING TECHNOLOGY》 * |
| 孟遥: "Invar合金激光熔化沉积过程微观组织建模与仿真研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅰ辑)》 * |
| 郝树青: "《空气反循环钻进流场理论及数值计算》", 30 September 2014, 徐州：中国矿业大学出版社 * |
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111037143A (en) * | 2019-11-28 | 2020-04-21 | 新兴际华集团有限公司 | Method for regulating and controlling weld joint structure performance based on welding thermal cycle |
| CN111037143B (en) * | 2019-11-28 | 2022-09-06 | 新兴际华集团有限公司 | Method for regulating and controlling weld joint structure performance based on welding thermal cycle |
| CN111168067A (en) * | 2020-01-22 | 2020-05-19 | 南京理工大学 | Pore prediction and control method based on laser directional energy deposition |
| CN111203538A (en) * | 2020-04-22 | 2020-05-29 | 中国航发上海商用航空发动机制造有限责任公司 | Prefabricated crack defect, preparation method of built-in crack defect and prefabricated part |
| CN111203537A (en) * | 2020-04-22 | 2020-05-29 | 中国航发上海商用航空发动机制造有限责任公司 | Method for controlling defect of poor fusion of LMD process prefabrication |
| CN111203538B (en) * | 2020-04-22 | 2020-07-28 | 中国航发上海商用航空发动机制造有限责任公司 | Prefabricated crack defect, preparation method of built-in crack defect and prefabricated part |
| CN111203537B (en) * | 2020-04-22 | 2020-07-28 | 中国航发上海商用航空发动机制造有限责任公司 | Method for prefabricating fusion defect by controlling L MD process |
| WO2021212847A1 (en) * | 2020-04-22 | 2021-10-28 | 中国航发上海商用航空发动机制造有限责任公司 | Methods for preparing prefabricated crack defect and built-in crack defect, and prefabricated member |
| CN111546636B (en) * | 2020-04-27 | 2021-06-29 | 中国水利水电科学研究院 | Simulation monitoring method for 3D printing of large structural body |
| CN111546636A (en) * | 2020-04-27 | 2020-08-18 | 中国水利水电科学研究院 | Simulation monitoring method for 3D printing of large structural body |
| CN112182911A (en) * | 2020-10-26 | 2021-01-05 | 同济大学 | 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 |
| CN112517924A (en) * | 2020-11-11 | 2021-03-19 | 三峡大学 | 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 |
| CN112989626A (en) * | 2021-04-13 | 2021-06-18 | 清华大学 | Additive manufacturing organization simulation method, device, computer equipment and storage medium |
| CN112989626B (en) * | 2021-04-13 | 2022-12-23 | 清华大学 | Additive manufacturing organization simulation method, device, computer equipment and storage medium |
| CN113158368A (en) * | 2021-04-16 | 2021-07-23 | 西安交通大学 | Prediction and control method and device for deformation of material-increasing and material-decreasing composite manufacturing part |
| CN113118458A (en) * | 2021-04-20 | 2021-07-16 | 江西省科学院应用物理研究所 | Prediction method for tensile property of metal component formed by selective laser melting |
| WO2022233206A1 (en) * | 2021-05-07 | 2022-11-10 | 苏州奇流信息科技有限公司 | Sample database system, method for training and checking printing parameter, and computer |
| CN114872324A (en) * | 2022-04-15 | 2022-08-09 | 华中科技大学 | Laser additive manufacturing method based on multidimensional information coupling regulation and control performance |
| CN114872324B (en) * | 2022-04-15 | 2023-09-29 | 华中科技大学 | Laser additive manufacturing method based on multidimensional information coupling regulation performance |
| CN115186557A (en) * | 2022-07-22 | 2022-10-14 | 南京航空航天大学 | Additive manufacturing component microstructure homogenization regulation and control method based on multi-physical-field coupling solution heat accumulation effect |
| 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 |
| CN117648876A (en) * | 2024-01-30 | 2024-03-05 | 中南大学 | TPMS gradient hierarchical structure inverse design manufacturing method and device based on performance and BALANCE-CGAN, storage medium and electronic device |
| CN117648876B (en) * | 2024-01-30 | 2024-04-12 | 中南大学 | TPMS gradient hierarchical structure inverse design manufacturing method based on performance and BALANCE-CGAN |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110472355B (en) | 2021-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110472355A (en) | A kind of 3D printing method for previewing solved based on multi- scenarios method modeling and simulation | |
| Xinhua et al. | An investigation on distortion of PLA thin-plate part in the FDM process | |
| DE60006701T2 (en) | METHOD AND COMPUTER DEVICE FOR MODELING THE SPLASHING OF A FLUIDUM INTO A FORM CAVITY | |
| US9939394B2 (en) | Process mapping of cooling rates and thermal gradients | |
| Di Angelo et al. | A reliable build orientation optimization method in additive manufacturing: The application to FDM technology | |
| KR101996933B1 (en) | A method for determining stresses and shape deviations in a stacked structure, a computer-readable data carrier, a computer program, and a simulator | |
| CN108399307A (en) | A kind of laser 3D printing Finite Element Method | |
| CN108062432A (en) | A kind of method for numerical simulation of selective laser fusion process | |
| CN108292325A (en) | Based on heat/structural simulation of the component to being produced via 3D printer come the system and method in optimization tool path | |
| CN109513924B (en) | Surface roughness control method in selective laser melting process | |
| Park et al. | Design and simulation-based optimization of cooling channels for plastic injection mold | |
| CN109284524A (en) | A method of creation high-precision increasing material manufacturing finite element model | |
| Boissier et al. | Additive manufacturing scanning paths optimization using shape optimization tools | |
| US11073824B1 (en) | System and method of simulating and optimizing surface quality based on location and orientation of additively manufactured build parts | |
| US20150286757A1 (en) | Method for Efficiently Predicting the Quality of Additively Manufactured Metal Products | |
| CN115186557B (en) | Additive manufacturing component microstructure homogenization regulation and control method based on multi-physical field coupling solution heat accumulation effect | |
| CN110245449A (en) | A kind of magnesium alloy cast part uneven components Numerical Predicting Method | |
| Danylchenko | Comparative analysis of computer systems for casting processes simulation | |
| CN111036909B (en) | Grain structure numerical value prediction method in metal rapid forming process | |
| CN106407547A (en) | Numerical simulation method for aiming at casting residual stress of anisotropic material | |
| Ravi | Metal casting 4.0: closing the loop between design and manufacturing | |
| CN115408842A (en) | Method for analyzing similarity of characteristics of injection molding simulation process | |
| Garanger et al. | Foundations of intelligent additive manufacturing | |
| Sideris et al. | End-to-end path planning for homogeneous temperature fields in additive manufacturing | |
| CN112084723A (en) | Method and device for simulating one-step stretching process of optical fiber preform |
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 |