CN110598358A - Additive manufacturing stress deformation simulation method, device, equipment and storage medium - Google Patents
Additive manufacturing stress deformation simulation method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN110598358A CN110598358A CN201910919933.3A CN201910919933A CN110598358A CN 110598358 A CN110598358 A CN 110598358A CN 201910919933 A CN201910919933 A CN 201910919933A CN 110598358 A CN110598358 A CN 110598358A
- Authority
- CN
- China
- Prior art keywords
- heat source
- parameters
- source model
- molten pool
- acquiring
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The invention discloses a simulation calculation method, a device, equipment and a computer storage medium for additive manufacturing stress deformation, wherein the method comprises the following steps: acquiring the current appearance of a molten pool in the forming process of a structural part; matching parameters of a proper point-like heat source model according to the current appearance; according to the parameters of the point-like heat source model, constructing a linear heat source model based on energy loading distribution; acquiring parameters and time step length of a linear heat source model, and acquiring energy temperature field distribution of a structural member according to the parameters and the time step length of the linear heat source model; and obtaining the stress deformation distribution of the structural member according to the energy temperature field distribution. The invention carries out simulation calculation by means of energy loading distribution, and can greatly reduce the calculation amount on the basis of not influencing the accuracy of the calculation result.
Description
Technical Field
The invention relates to the technical field of additive manufacturing numerical simulation, in particular to a method, a device, equipment and a storage medium for additive manufacturing stress deformation simulation.
Background
The so-called additive manufacturing stress deformation simulation may be considered to be a process in which the stress deformation distribution is obtained by simulating the entire additive manufacturing forming process by a computer. Because the stress distribution of the formed workpiece is closely related to the temperature field distribution in the forming process, the calculation of the temperature field is particularly important in the actual numerical simulation process. However, in the additive manufacturing process, the heat source has the characteristics of concentration and movement, and an uneven temperature field with a large gradient in space and time can be formed. Due to the high concentration of the heat source, if the traditional moving point-like heat source is directly adopted for calculation, when a finite element model is established, grid units need to be finely divided, and a large amount of time steps are needed for iterative calculation in the calculation, so that the calculation efficiency is extremely low, and the forming process of some components with complicated structures and large sizes cannot be simulated.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a method, an apparatus, a device, and a storage medium for simulating additive manufacturing stress deformation, where the method performs simulation calculation in an energy loading distribution manner, so that the calculation amount can be greatly reduced without affecting the accuracy of the calculation result, and the problems of excessive calculation amount and high difficulty in the additive manufacturing numerical simulation are solved.
The embodiment of the invention provides a simulation method for stress deformation in additive manufacturing, which comprises the following steps:
acquiring the current appearance of a molten pool in the forming process of a structural part;
matching parameters of a proper point-like heat source model according to the current appearance;
according to the parameters of the point-like heat source model, constructing a linear heat source model based on energy loading distribution;
acquiring parameters and time step length of a linear heat source model, and acquiring energy temperature field distribution of a structural member according to the parameters and the time step length of the linear heat source model;
and obtaining the stress deformation distribution of the structural member according to the energy temperature field distribution.
Preferably, the step of acquiring the current appearance of the molten pool is preceded by the step of acquiring the current appearance of the molten pool
Acquiring a finite element model generated after a pre-established three-dimensional geometric model of a structural part is subjected to network division;
acquiring forming process parameters of the corresponding structural part according to the finite element analysis model;
and initializing parameters of a molten pool in the forming process of the structural member according to the forming process parameters of the structural member so as to obtain the current appearance of the molten pool.
Preferably, the forming process parameters include the width of the forming path and the height of each layer of the forming path, and the parameters of the molten pool in the forming process of the structural member are initialized according to the forming process parameters of the structural member to obtain the current shape of the molten pool, specifically:
according to the width of the forming path, obtaining the width of the molten pool;
according to the height of each layer of the forming path, obtaining the depth of the molten pool;
and according to the width of the molten pool and the depth of the molten pool, acquiring the current appearance of the molten pool through a process-shaped energy input mode.
Preferably, according to the parameters of the point-like heat source model, a linear heat source model is constructed based on energy loading distribution, specifically:
obtaining the maximum value of the heat flux density of the linear heat source according to the heat input power, the heat efficiency coefficient, the welding voltage, the welding current and the parameters of the point heat source model;
obtaining the heat flux density of each position in the internal space of the linear heat source according to the maximum value of the heat flux density of the linear heat source and the parameters of the point heat source model;
according to the motion speed and the parameters of the point-shaped heat source model, the heating time of the linear heat source is obtained;
and constructing a linear heat source model according to the heating time of the linear heat source and the heat flux density of each position in the internal space of the linear heat source.
Preferably, the linear heat source model tableThe expression is as follows: q. q.ss(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；Wherein q issIs the heat flux density of each position in the internal space of the heat source, a and b are the width and height of the point heat source model, QmIs the heat input power, QmEta UI, eta is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, qsmIs the maximum value of the heat flux density, t, of the linear heat sourcesThe heating time of the linear heat source model is shown.
Preferably, the time step of the linear heat source model is 10 times of the time step of the punctiform heat source model.
The embodiment of the invention also provides a material increase manufacturing stress deformation simulation device, which comprises:
the current appearance obtaining unit is used for obtaining the current appearance of a molten pool in the forming process of the structural part;
the parameter matching unit is used for matching the parameters of the proper point-like heat source model according to the current appearance;
the linear heat source model building unit is used for building a linear heat source model based on energy loading distribution according to the parameters of the point heat source model;
the energy temperature field distribution acquisition unit is used for acquiring parameters and time step length of the linear heat source model and acquiring the energy temperature field distribution of the structural member according to the parameters and the time step length of the linear heat source model;
and the stress deformation distribution acquisition unit is used for acquiring the stress deformation distribution of the structural member according to the energy temperature field distribution.
Preferably, also comprises
The finite element model acquisition unit is used for acquiring a finite element model generated after a pre-established three-dimensional geometric model of a structural part is subjected to network division;
the forming process parameter acquiring unit is used for acquiring forming process parameters of the corresponding structural part according to the finite element analysis model;
and the initialization unit is used for initializing the parameters of the molten pool in the forming process of the structural member according to the forming process parameters of the structural member so as to obtain the current appearance of the molten pool.
Preferably, the forming process parameters include the width of the forming path and the height of each layer of the forming path, and the parameters of the molten pool in the forming process of the structural member are initialized according to the forming process parameters of the structural member to obtain the current shape of the molten pool, specifically:
according to the width of the forming path, obtaining the width of the molten pool;
according to the height of each layer of the forming path, obtaining the depth of the molten pool;
and according to the width of the molten pool and the depth of the molten pool, acquiring the current appearance of the molten pool through a process-shaped energy input mode.
Preferably, according to the parameters of the point-like heat source model, a linear heat source model is constructed based on energy loading distribution, specifically:
obtaining the maximum value of the heat flux density of the linear heat source according to the heat input power, the heat efficiency coefficient, the welding voltage, the welding current and the parameters of the point heat source model;
obtaining the heat flux density of each position in the internal space of the linear heat source according to the maximum value of the heat flux density of the linear heat source and the parameters of the point heat source model;
according to the motion speed and the parameters of the point-shaped heat source model, the heating time of the linear heat source is obtained;
and constructing a linear heat source model according to the heating time of the linear heat source and the heat flux density of each position in the internal space of the linear heat source.
Preferably, the linear heat source model expression is: q. q.ss(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；Wherein q issIs the heat flux density at each position of the internal space of the heat source, a and b are respectively point heat sourcesWidth and height of the model, QmIs the heat input power, QmEta UI, eta is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, qsmIs the maximum value of the heat flux density, t, of the linear heat sourcesThe heating time of the linear heat source model is shown.
Preferably, the time step of the linear heat source model is 10 times of the time step of the punctiform heat source model.
An embodiment of the present invention further provides an additive manufacturing stress deformation simulation method apparatus, including a processor, a memory, and a computer program stored in the memory, where the computer program is executable by the processor to implement the additive manufacturing stress deformation simulation method according to the first aspect.
The embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, where when the computer program runs, a device on which the computer-readable storage medium is located is controlled to execute the additive manufacturing stress deformation simulation method described above.
The embodiment of the invention has the following beneficial effects:
1. according to the method, the linear heat source model is constructed based on the energy loading distribution according to the parameters of the point-like heat source model, and the simulation calculation is carried out in the energy loading distribution mode, so that the calculation amount is greatly reduced on the basis of not influencing the accuracy of the calculation result, and the problems of overlarge calculation amount and high difficulty in the numerical simulation of the material increase manufacturing are solved.
2. The invention reduces the requirement on the size of the additive manufacturing stress simulation numerical simulation grid, does not need to divide a large number of dense grids to ensure the accuracy of the calculation result, and can divide the grid size along the forming direction greatly because the heat flux density of the linear distribution heat source along the forming direction is uniformly distributed, thereby greatly reducing the total grid number.
3. The method reduces the requirement on the numerical simulation time step length of the material increase manufacturing stress simulation, does not need to approximately fit the whole forming process by dividing dense time step lengths, and greatly reduces the calculated amount.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a simulation method of additive manufacturing stress deformation according to a first embodiment of the present invention.
Fig. 2 is a schematic structural diagram of a linear heat source model according to a second embodiment of the present invention.
Fig. 3 is a schematic structural diagram of an additive manufacturing stress deformation simulation method and apparatus according to a second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
For better understanding of the technical solutions of the present invention, the following detailed descriptions of the embodiments of the present invention are provided with reference to the accompanying drawings.
It should be understood that the described embodiments are only some embodiments of the invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The terminology used in the embodiments of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the examples of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It should be understood that the term "and/or" as used herein is merely one type of association that describes an associated object, meaning that three relationships may exist, e.g., a and/or B may mean: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
The word "if" as used herein may be interpreted as "at … …" or "when … …" or "in response to a determination" or "in response to a detection", depending on the context. Similarly, the phrases "if determined" or "if detected (a stated condition or event)" may be interpreted as "when determined" or "in response to a determination" or "when detected (a stated condition or event)" or "in response to a detection (a stated condition or event)", depending on the context.
In the embodiments, the references to "first \ second" are merely to distinguish similar objects and do not represent a specific ordering for the objects, and it is to be understood that "first \ second" may be interchanged with a specific order or sequence, where permitted. It should be understood that "first \ second" distinct objects may be interchanged under appropriate circumstances such that the embodiments described herein may be practiced in sequences other than those illustrated or described herein.
The first embodiment is as follows:
referring to fig. 1 and fig. 2, a first embodiment of the present invention provides an additive manufacturing stress-deformation simulation method, which can be executed by an additive manufacturing stress-deformation simulation apparatus, and in particular, by one or more processors in the additive manufacturing stress-deformation simulation apparatus, and at least includes the following steps:
s101, acquiring the current appearance of a molten pool in the forming process of the structural part.
In this embodiment, the additive manufacturing stress deformation simulation equipment initializes parameters of a molten pool in a forming process of the structural member according to the forming process parameters of the structural member, so as to obtain a current shape of the molten pool. Wherein the current appearance of the molten pool comprises the penetration depth, the width and the shape of the molten pool. The forming process parameters include wire feed speed, heat source power, heat source radius, scanning speed, forming path width, and forming path height per layer, etc., which enable a simulated environment to approach an actual environment indefinitely. Specifically, according to the width of the forming path, to obtain the width of the molten pool; and according to the height of each layer of the forming path, obtaining the depth of the molten pool; and then, acquiring the current appearance of the molten pool through a process-shaped energy input mode according to the width of the molten pool and the depth of the molten pool.
It should be noted that the forming process parameters are obtained according to the finite element analysis model, specifically, a finite element model generated by network division of a pre-established three-dimensional geometric model of a structural member is obtained by the additive manufacturing stress deformation simulation equipment; and acquiring the forming process parameters of the corresponding structural part according to the finite element analysis model, and initializing the parameters of a molten pool in the forming process of the structural part according to the forming process parameters of the structural part so as to obtain the current appearance of the molten pool. The geometric model of the structural part is built in three-dimensional modeling software, then grid division software is introduced into the three-dimensional modeling software to divide grids, and it needs to be explained that in order to greatly enlarge the grid size required by the finite element model and reduce the grid number of the grids, the grid type is hexahedral grid, the grid density is dense at the position close to the center of a welding seam, and sparse at the position far away from the welding seam, so that the calculation amount can be effectively reduced compared with the grid division required by the traditional heat source model.
And S102, matching parameters of a proper point-like heat source model according to the current appearance.
In this embodiment, since the radius of the point-like heat source model is smaller than half of the width of the molten pool during the material addition process, in order to match an appropriate point-like heat source model parameter, the initial value of the radius of the point-like heat source model is always half of the width of the molten pool, and based on this, the appropriate radius of the point-like heat source model is obtained according to the width of the current appearance of the molten pool; and obtaining the height of a proper point-like heat source model according to the height of the current appearance of the molten pool.
S103, constructing a linear heat source model based on energy loading distribution according to the parameters of the point heat source model.
In this embodiment, when a conventional heat loading manner is used for simulation, since the dot heat source has a small volume, a relatively dense grid needs to be divided at a forming path part, and a fine time step is adopted to divide the whole forming process, a large amount of time steps are required for calculation, resulting in a huge calculation amount, and therefore, according to the parameters of the dot heat source model including the radius and height of the dot heat source model, the heat input power, the movement speed and the like, the maximum value of the heat flow density of the linear heat source is obtained according to the radius of the dot heat source model, the height of the dot heat source model, the heat input power, the heat efficiency coefficient, the welding voltage and the welding current; obtaining the heat flux density of each position in the internal space of the linear heat source according to the maximum value of the heat flux density of the linear heat source, the radius of the point heat source model and the height of the point heat source model; according to the movement speed and the radius of the point-shaped heat source model, the heating time of the linear heat source is obtained; and constructing a linear heat source model according to the heating time of the linear heat source and the heat flux density of each part in the internal space of the linear heat source, and establishing the linear heat source model in an energy loading mode, so that the large calculation amount brought by the original point heat source model is reduced.
Wherein, the energy loading distribution equation is as follows: q. q.ss(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2) (ii) a The linear heat source model expression is as follows: q. q.ss(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；qsIs the heat flux density of each position in the internal space of the heat source, a and b are the radius and height of the point heat source model, QmIs the heat input power, QmEta UI, eta is the thermal efficiency coefficient, U is the welding voltage, I is the welding currentStream, qsmIs the maximum value of the heat flux density, t, of the linear heat sourcesThe heating time of the linear heat source model is shown.
And S104, acquiring parameters and time step length of the linear heat source model, and acquiring energy temperature field distribution of the structural member according to the parameters and the time step length of the linear heat source model.
In this embodiment, the time step of the linear heat source model is 10 times that of the point-like heat source model; the time step length of the point-like heat source only needs to meet the simulation calculation conditions of the finite element model, namely:where Δ x is the size of the divided finite element mesh, cpIs the specific heat capacity. Of course, the time step of the linear heat source model may also be 8 times or 7 times that of the point heat source model, and the like, and the description of the present invention is omitted here.
In this embodiment, since the linear heat source model parameters and the time step are the conditions that must be set in advance for the solver to calculate the temperature field, the linear heat source model parameters and the time step are input to the solver to calculate the energy temperature field distribution of the structural member, wherein the solver is software for calculating the temperature field, and ABAQUS is used in the present invention.
And S105, obtaining stress deformation distribution of the structural member according to the energy temperature field distribution.
In this embodiment, the final stress-strain distribution of the structural member is obtained by thermal-stress coupling based on the energy-temperature field distribution.
In conclusion, according to the parameters of the point-like heat source model, the linear heat source model is constructed based on the energy loading distribution, and the simulation calculation is carried out in the energy loading distribution mode, so that the calculation amount can be greatly reduced on the basis of not influencing the accuracy of the calculation result, and the problems of overlarge calculation amount and high difficulty in the numerical simulation of the additive manufacturing are solved. Meanwhile, the method reduces the requirement on the size of the additive manufacturing stress simulation numerical simulation grid, does not need to divide a large number of dense grids to ensure the accuracy of a calculation result, and can divide the size of the grid in the forming direction greatly because the heat flow density of the linear distribution heat source in the forming direction is uniformly distributed, thereby greatly reducing the total grid number. The invention reduces the requirement on the numerical simulation time step length of the material increase manufacturing stress simulation, does not need to approximately fit the whole forming process by dividing dense time step lengths, and greatly reduces the calculated amount.
For the convenience of understanding of the present invention, the following description will be made of an application of the present embodiment in a practical application scenario.
Assuming that the formed piece is a short beam piece, the whole short beam forming is divided into two stages: a base portion and a rib portion, wherein the size of the base is 563X 80X 20mm, and the base portion is stacked in 8 layers. The stacking part on the substrate is ribs, the transverse width in the x direction is 12mm, the longitudinal width in the y direction is 8mm, and the height is 50 mm. The layers 33 are stacked according to different process paths.
The base portion and rib portion forming process parameters are shown in tables 1 and 2 below.
Table 1: basic parameters of matrix forming (CX) process
| Focusing current | Scanning amplitude | Platform lift | Beam current | Speed of movement | Wire feed speed | Scanning direction | Scanning pitch |
| 900/490 | 300/300 | 30mm | 67mA | 200mm/min | 21mm/s | X or Y | 7mm |
Table 2: fundamental parameters of rib forming (JT) process
| Focusing current | Scanning amplitude | Platform lift | Beam current | Speed of movement | Wire feed speed | Scanning direction | Scanning pitch |
| 900/490 | 300/300 | 30mm | 67mA | 200mm/min | 21mm/s | ——— | ——— |
Specifically, in step S1, a solid three-dimensional finite element model is established, a geometric model of the workpiece is established in the three-dimensional modeling software, then a mesh is divided by introducing the three-dimensional finite element model into mesh division software, linear hexahedral elements are adopted, the element size is sparsely divided, and the total number of the elements is 38962. In step S2, the width of the molten pool is determined to be about 7mm, the depth is determined to be 8mm, and the cross section is determined to be a double ellipsoid according to the process parameters in the table. In step S3, the radius width of the heat source model is determined to be 5mm and the depth is determined to be 7mm based on the molten pool width depth parameter and the shape. In step S4, the linear heat source model is determined to have a of 5mm and b of 7mm based on the parameters of the point-like heat source in step S3, and the heat input power Q is determinedm1900W, speed of movement vm3.33 mm/s. Then, constructing a linear heat source model based on energy loading distribution, wherein an energy loading distribution equation is as follows:x, y, z are coordinates (mm). In step S5, the time step is set to 0.5S, and the time step is input into a solver, which is a temperature field calculation module of the ABAQUS software in this embodiment. In step S6, the temperature field distribution of the component is finally solved in the solver based on the linear heat source model parameters, the time step, and various conditions necessary for calculating the temperature field, including various boundary conditions. In step S7, the final stress deformation distribution of the member is obtained by thermal-stress coupling based on the temperature field distribution obtained in step S6, and the maximum deformation is 10.9mm and the warpage at both ends is 5.1 mm.
And comparing the simulation calculation result with the workpiece obtained after actual forming, wherein the maximum deformation is about 1.5mm after the base body is formed (incompletely cooled, the temperature is about 200 ℃), and the whole body is in a wavy deformation with a middle tilted state. In the rib forming process, the influence of the single-layer rib forming process on the integral deformation of the workpiece is small. Wherein, in the process of forming the transverse short ribs, the deformation is gradually increased, and the change interval of the integral deformation of the single layer is less than 0.3 mm. The deformation trend of the whole workpiece after forming shows a mode that two ends are tilted and the middle is concave, and the tilting deformation amount of the two ends is 5.7 mm. The deformation of the base plate is large, the maximum deformation is 11.7mm, and the overall deformation trend is wavy. The maximum deformation measured finally in the experiment is 10.9mm, and the tilting deformation of the two ends is 5.1mm, so that the deformation trend is the same as that obtained by simulation, the establishment of a linear heat source model for accelerating the simulation efficiency of the material additive manufacturing stress can be well realized, and the time of simulation calculation can be greatly shortened.
Second embodiment of the invention:
referring to fig. 3, an embodiment of the present invention further includes an additive manufacturing stress deformation simulation method and apparatus, including:
the current appearance obtaining unit 100 is used for obtaining the current appearance of a molten pool in the forming process of a structural part;
the parameter matching unit 200 is used for matching the parameters of the appropriate point-like heat source model according to the current appearance;
a linear heat source model building unit 300, configured to build a linear heat source model based on energy loading distribution according to the parameters of the dotted heat source model;
the energy temperature field distribution acquisition unit 400 is used for acquiring parameters and time step length of the linear heat source model and acquiring the energy temperature field distribution of the structural member according to the parameters and the time step length of the linear heat source model;
and the stress deformation distribution acquisition unit 500 is used for acquiring the stress deformation distribution of the structural member according to the energy temperature field distribution.
On the basis of the above embodiments, a preferred embodiment of the present invention further includes
The finite element model acquisition unit is used for acquiring a finite element model generated after a pre-established three-dimensional geometric model of a structural part is subjected to network division;
the forming process parameter acquiring unit is used for acquiring forming process parameters of the corresponding structural part according to the finite element analysis model;
and the initialization unit is used for initializing the parameters of the molten pool in the forming process of the structural member according to the forming process parameters of the structural member so as to obtain the current appearance of the molten pool.
On the basis of the above embodiment, in a preferred embodiment of the present invention, the forming process parameters include a width of the forming path and a height of each layer of the forming path, and the parameters of the molten pool in the forming process of the structural member are initialized according to the forming process parameters of the structural member to obtain a current shape of the molten pool, specifically:
according to the width of the forming path, obtaining the width of the molten pool;
according to the height of each layer of the forming path, obtaining the depth of the molten pool;
and according to the width of the molten pool and the depth of the molten pool, acquiring the current appearance of the molten pool through a process-shaped energy input mode.
On the basis of the foregoing embodiment, in a preferred embodiment of the present invention, a linear heat source model is constructed based on energy loading distribution according to parameters of the point-like heat source model, specifically:
obtaining the maximum value of the heat flux density of the linear heat source according to the heat input power, the heat efficiency coefficient, the welding voltage, the welding current and the parameters of the point heat source model;
obtaining the heat flux density of each position in the internal space of the linear heat source according to the maximum value of the heat flux density of the linear heat source and the parameters of the point heat source model;
according to the motion speed and the parameters of the point-shaped heat source model, the heating time of the linear heat source is obtained;
and constructing a linear heat source model according to the heating time of the linear heat source and the heat flux density of each position in the internal space of the linear heat source.
On the basis of the foregoing embodiment, in a preferred embodiment of the present invention, the linear heat source model expression is: q. q.ss(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；Wherein q issIs the heat flux density of each position in the internal space of the heat source, a and b are the width and height of the point heat source model, QmIs the heat input power, QmEta UI, eta is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, qsmIs the maximum value of the heat flux density, t, of the linear heat sourcesThe heating time of the linear heat source model is shown.
On the basis of the foregoing embodiments, in a preferred embodiment of the present invention, the time step of the linear heat source model is 10 times that of the point-like heat source model.
Third embodiment of the invention:
a third embodiment of the present invention provides an additive manufacturing stress-deformation simulation method apparatus, including a processor, a memory, and a computer program stored in the memory, the computer program being executable by the processor to implement the additive manufacturing stress-deformation simulation method as described above.
The fourth embodiment of the present invention:
a fourth embodiment of the present invention provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program, where the computer program, when running, controls a device on which the computer-readable storage medium is located to execute the additive manufacturing stress deformation simulation method as described above.
Illustratively, the computer program may be divided into one or more units, which are stored in the memory and executed by the processor to accomplish the present invention. The one or more units may be a series of computer program instruction segments capable of performing specific functions, and the instruction segments are used for describing the execution process of the computer program in the additive manufacturing stress deformation simulation method device.
The additive manufacturing stress deformation simulation method and device can comprise but not limited to a processor and a memory. It will be understood by those skilled in the art that the schematic diagrams are merely examples of an additive manufacturing stress deformation simulation method apparatus and do not constitute a limitation of an additive manufacturing stress deformation simulation method apparatus, and may include more or fewer components than those shown, or combine certain components, or different components, for example, the additive manufacturing stress deformation simulation method apparatus may further include an input-output device, a network access device, a bus, etc.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may be any conventional processor or the like, and the control center of the additive manufacturing stress deformation simulation method apparatus connects the various parts of the entire additive manufacturing stress deformation simulation method apparatus by using various interfaces and lines.
The memory may be configured to store the computer program and/or the module, and the processor may implement various functions of the additive manufacturing stress deformation simulation method apparatus by executing or executing the computer program and/or the module stored in the memory and calling up data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The unit integrated by the additive manufacturing stress deformation simulation method and the device can be stored in a computer readable storage medium if the unit is realized in the form of a software functional unit and sold or used as an independent product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, etc. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
It should be noted that the above-described device embodiments are merely illustrative, where the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. In addition, in the drawings of the embodiment of the apparatus provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, and may be specifically implemented as one or more communication buses or signal lines. One of ordinary skill in the art can understand and implement it without inventive effort.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.
Claims (10)
1. An additive manufacturing stress deformation simulation method is characterized by comprising the following steps:
acquiring the current appearance of a molten pool in the forming process of a structural part;
matching parameters of a proper point-like heat source model according to the current appearance;
according to the parameters of the point-like heat source model, constructing a linear heat source model based on energy loading distribution;
acquiring parameters and time step length of a linear heat source model, and acquiring energy temperature field distribution of a structural member according to the parameters and the time step length of the linear heat source model;
and obtaining the stress deformation distribution of the structural member according to the energy temperature field distribution.
2. The method of additive manufacturing stress-deformation simulation of claim 1, wherein the step of obtaining a current topography of the melt pool is preceded by the step of obtaining a current topography of the melt pool further comprising
Acquiring a finite element model generated after a pre-established three-dimensional geometric model of a structural part is subjected to network division;
acquiring forming process parameters of the corresponding structural part according to the finite element analysis model;
and initializing parameters of a molten pool in the forming process of the structural member according to the forming process parameters of the structural member so as to obtain the current appearance of the molten pool.
3. The method according to claim 2, wherein the forming process parameters include a width of a forming path and a height of each layer of the forming path, and parameters of a molten pool in a forming process of the structural member are initialized according to the forming process parameters of the structural member to obtain a current shape of the molten pool, specifically:
according to the width of the forming path, obtaining the width of the molten pool;
according to the height of each layer of the forming path, obtaining the depth of the molten pool;
and according to the width of the molten pool and the depth of the molten pool, acquiring the current appearance of the molten pool through a process-shaped energy input mode.
4. The additive manufacturing stress-deformation simulation method of claim 3,
according to the parameters of the point-like heat source model, constructing a linear heat source model based on energy loading distribution, specifically:
obtaining the maximum value of the heat flux density of the linear heat source according to the heat input power, the heat efficiency coefficient, the welding voltage, the welding current and the parameters of the point heat source model;
obtaining the heat flux density of each position in the internal space of the linear heat source according to the maximum value of the heat flux density of the linear heat source and the parameters of the point heat source model;
according to the motion speed and the parameters of the point-shaped heat source model, the heating time of the linear heat source is obtained;
and constructing a linear heat source model according to the heating time of the linear heat source and the heat flux density of each position in the internal space of the linear heat source.
5. The additive manufacturing stress-deformation simulation method of claim 4, wherein the linear heat source model expression is: q. q.ss(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；Wherein q issThe heat flux density at each position of the internal space of the heat source, a and b are respectively a point heat source dieWidth and height of pattern, QmIs the heat input power, QmEta UI, eta is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, qsmIs the maximum value of the heat flux density, t, of the linear heat sourcesThe heating time of the linear heat source model is shown.
6. The additive manufacturing stress-deformation simulation method of claim 4,
the time step of the linear heat source model is 10 times of the time step of the punctiform heat source model.
7. An additive manufacturing stress deformation simulation device, comprising:
the current appearance obtaining unit is used for obtaining the current appearance of a molten pool in the forming process of the structural part;
the parameter matching unit is used for matching the parameters of the proper point-like heat source model according to the current appearance;
the linear heat source model building unit is used for building a linear heat source model based on energy loading distribution according to the parameters of the point heat source model;
the energy temperature field distribution acquisition unit is used for acquiring parameters and time step length of the linear heat source model and acquiring the energy temperature field distribution of the structural member according to the parameters and the time step length of the linear heat source model;
and the stress deformation distribution acquisition unit is used for acquiring the stress deformation distribution of the structural member according to the energy temperature field distribution.
8. The additive manufacturing stress-deformation simulation method of claim 7, further comprising
The finite element model acquisition unit is used for acquiring a finite element model generated after a pre-established three-dimensional geometric model of a structural part is subjected to network division;
the forming process parameter acquiring unit is used for acquiring forming process parameters of the corresponding structural part according to the finite element analysis model;
and the initialization unit is used for initializing the parameters of the molten pool in the forming process of the structural member according to the forming process parameters of the structural member so as to obtain the current appearance of the molten pool.
9. An additive manufacturing stress-deformation simulation method apparatus comprising a processor, a memory, and a computer program stored in the memory, the computer program being executable by the processor to implement the additive manufacturing stress-deformation simulation method of any one of claims 1 to 6.
10. A computer-readable storage medium comprising a stored computer program, wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform a method of computing an additive manufacturing stress-deformation simulation according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910919933.3A CN110598358A (en) | 2019-09-26 | 2019-09-26 | Additive manufacturing stress deformation simulation method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910919933.3A CN110598358A (en) | 2019-09-26 | 2019-09-26 | Additive manufacturing stress deformation simulation method, device, equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110598358A true CN110598358A (en) | 2019-12-20 |
Family
ID=68863768
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910919933.3A Pending CN110598358A (en) | 2019-09-26 | 2019-09-26 | Additive manufacturing stress deformation simulation method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110598358A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111291512A (en) * | 2020-01-22 | 2020-06-16 | 南京理工大学 | Design method of variable-density Gyroid lattice structure with heat insulation function |
| CN111687553A (en) * | 2020-05-28 | 2020-09-22 | 西安交通大学 | Method for improving residual stress distribution of electric arc additive manufacturing structural part through ultrasonic impact |
| CN112287489A (en) * | 2020-11-06 | 2021-01-29 | 中国石油大学(北京) | Method for determining parameters of arc fuse additive manufacturing process based on multi-physical-field simulation |
| CN112507583A (en) * | 2020-11-24 | 2021-03-16 | 深圳技术大学 | SF (sulfur hexafluoride)6Finite element simulation analysis method for arc plasma multi-physical field in opening process of breaker contact |
| CN113139314A (en) * | 2021-04-29 | 2021-07-20 | 四川大学 | Heat source numerical simulation method for laser additive manufacturing process |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101334807A (en) * | 2008-07-28 | 2008-12-31 | 中国航空工业第一集团公司北京航空制造工程研究所 | Electro-beam welding joint melting-coagulation area shape factor modeling and solving method |
| JP2009250829A (en) * | 2008-04-08 | 2009-10-29 | Radioactive Waste Management Funding & Research Center | Method for simple three-dimensional analysis of welding deformation and residual stress |
| JP2010175093A (en) * | 2009-01-27 | 2010-08-12 | Taikisha Ltd | Heat source equipment control system |
| US20110213594A1 (en) * | 2010-02-03 | 2011-09-01 | Hitachi, Ltd. | Method for simulation of welding distortion |
| CN102728932A (en) * | 2012-06-11 | 2012-10-17 | 江苏新瑞峰信息科技有限公司 | Method and device for automatically optimizing electric-arc welding process parameters |
| CN103049623A (en) * | 2013-01-18 | 2013-04-17 | 哈尔滨工业大学 | Building method for laser welding heat source model |
| CN104268343A (en) * | 2014-09-28 | 2015-01-07 | 北京理工大学 | Method of cutting force prediction and temperature prediction for end-milling cutting |
| CN106909714A (en) * | 2017-01-19 | 2017-06-30 | 西南交通大学 | A kind of thin-wall member electric arc silk filling increasing material manufacturing temperature field prediction method |
| CN107368642A (en) * | 2017-07-13 | 2017-11-21 | 武汉大学 | The multiple dimensioned multiple physical field coupling simulation method of metal increasing material manufacturing |
| CN108062432A (en) * | 2017-11-17 | 2018-05-22 | 西安铂力特增材技术股份有限公司 | A kind of method for numerical simulation of selective laser fusion process |
| CN108334692A (en) * | 2018-01-30 | 2018-07-27 | 清华大学 | The method for predicting increasing material manufacturing part deformation |
-
2019
- 2019-09-26 CN CN201910919933.3A patent/CN110598358A/en active Pending
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009250829A (en) * | 2008-04-08 | 2009-10-29 | Radioactive Waste Management Funding & Research Center | Method for simple three-dimensional analysis of welding deformation and residual stress |
| CN101334807A (en) * | 2008-07-28 | 2008-12-31 | 中国航空工业第一集团公司北京航空制造工程研究所 | Electro-beam welding joint melting-coagulation area shape factor modeling and solving method |
| JP2010175093A (en) * | 2009-01-27 | 2010-08-12 | Taikisha Ltd | Heat source equipment control system |
| US20110213594A1 (en) * | 2010-02-03 | 2011-09-01 | Hitachi, Ltd. | Method for simulation of welding distortion |
| CN102728932A (en) * | 2012-06-11 | 2012-10-17 | 江苏新瑞峰信息科技有限公司 | Method and device for automatically optimizing electric-arc welding process parameters |
| CN103049623A (en) * | 2013-01-18 | 2013-04-17 | 哈尔滨工业大学 | Building method for laser welding heat source model |
| CN104268343A (en) * | 2014-09-28 | 2015-01-07 | 北京理工大学 | Method of cutting force prediction and temperature prediction for end-milling cutting |
| CN106909714A (en) * | 2017-01-19 | 2017-06-30 | 西南交通大学 | A kind of thin-wall member electric arc silk filling increasing material manufacturing temperature field prediction method |
| CN107368642A (en) * | 2017-07-13 | 2017-11-21 | 武汉大学 | The multiple dimensioned multiple physical field coupling simulation method of metal increasing material manufacturing |
| CN108062432A (en) * | 2017-11-17 | 2018-05-22 | 西安铂力特增材技术股份有限公司 | A kind of method for numerical simulation of selective laser fusion process |
| CN108334692A (en) * | 2018-01-30 | 2018-07-27 | 清华大学 | The method for predicting increasing material manufacturing part deformation |
Non-Patent Citations (5)
| Title |
|---|
| LVJIE LIANG 等: "A CFD-FEM Model of Residual Stress for Electron Beam Welding Including the Weld Imperfection Effect", 《METALL MATER TRANS A》 * |
| 李良碧 编: "高斯热源模型", 《船舶结构疲劳强度与焊接残余应力》 * |
| 董明晔 等: "铝合金筒壁电弧增材制造数值模拟中分段弧形体热源模型的建立", 《HTTP://KNS.CNKI.NET/KCMS/DETAIL/11.2223.N.20190610.1654.004.HTML》 * |
| 贾金龙 等: "基于温度函数法的铝合金电弧增材制造残余应力与变形数值模拟", 《焊接学报》 * |
| 高金岭: "FDM快速成型机温度场及应力场的数值模拟仿真", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111291512A (en) * | 2020-01-22 | 2020-06-16 | 南京理工大学 | Design method of variable-density Gyroid lattice structure with heat insulation function |
| CN111291512B (en) * | 2020-01-22 | 2021-08-31 | 南京理工大学 | Design method of variable-density Gyroid lattice structure with heat insulation function |
| CN111687553A (en) * | 2020-05-28 | 2020-09-22 | 西安交通大学 | Method for improving residual stress distribution of electric arc additive manufacturing structural part through ultrasonic impact |
| CN112287489A (en) * | 2020-11-06 | 2021-01-29 | 中国石油大学(北京) | Method for determining parameters of arc fuse additive manufacturing process based on multi-physical-field simulation |
| CN112287489B (en) * | 2020-11-06 | 2022-10-21 | 中国石油大学(北京) | Method for determining parameters of arc fuse additive manufacturing process based on multi-physical-field simulation |
| CN112507583A (en) * | 2020-11-24 | 2021-03-16 | 深圳技术大学 | SF (sulfur hexafluoride)6Finite element simulation analysis method for arc plasma multi-physical field in opening process of breaker contact |
| CN113139314A (en) * | 2021-04-29 | 2021-07-20 | 四川大学 | Heat source numerical simulation method for laser additive manufacturing process |
| CN113139314B (en) * | 2021-04-29 | 2022-09-27 | 四川大学 | Heat source numerical simulation method for laser additive manufacturing process |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110598358A (en) | Additive manufacturing stress deformation simulation method, device, equipment and storage medium | |
| US11335074B2 (en) | Arrangement determination for 3D fabricated parts | |
| CN107633295B (en) | Method and device for adapting parameters of a neural network | |
| CN112669463A (en) | Method for reconstructing curved surface of three-dimensional point cloud, computer device and computer-readable storage medium | |
| CN107153724B (en) | Chip temperature analysis method based on iterative algorithm | |
| CN115510798B (en) | Chip typesetting method and device, computer readable storage medium and terminal equipment | |
| CN112749450A (en) | Simulation analysis method and system based on Hyperworks | |
| CN110532670B (en) | Method and system for constructing spline entity model suitable for analysis | |
| CN110705158B (en) | Prediction method, device and equipment for gas-liquid-solid three-phase behavior of selective laser melting | |
| CN114580249B (en) | Multi-loop FDTD electromagnetic field simulation analysis method, system, equipment and medium | |
| CN108733886B (en) | Lofting method and device for hemispherical-square ship structure and computer-readable storage medium | |
| CN113352619B (en) | Skirt edge generation method and device for starting position, computer equipment and storage medium | |
| CN115114777A (en) | Injection mold pouring simulation method, device, equipment and storage medium | |
| CN111737787B (en) | Calculation method, calculation system, storage medium and electronic equipment for bent materials | |
| CN112711231A (en) | Method and device for planning additive manufacturing processing path | |
| Adamoudis et al. | Heuristic repairing operators for 3D tetrahedral mesh generation using the advancing-front technique | |
| CN115034081B (en) | Foam concrete pore structure construction method, device, equipment and storage medium | |
| US7412359B2 (en) | Mutual immittance calculation program, mutual immittance calculation apparatus and method, and electromagnetic-field strength calculation program | |
| JP2006139723A (en) | Fdtd computation device and fdtd computation method | |
| CN108671417A (en) | Pencil beam Response characteristics based on self-consistency | |
| CN112426716B (en) | Three-dimensional hair model processing method, device, equipment and storage medium | |
| Chen et al. | Surface triangular mesh and volume tetrahedral mesh generations for biomolecular modeling | |
| CN118171539B (en) | Quasi-periodic array electromagnetic calculation method based on finite element boundary element region decomposition | |
| CN117934753A (en) | Grid model building method and device, electronic equipment and storage medium | |
| CN116305681A (en) | Envelope generation method, device, electronic equipment and storage medium |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Pang Shengyong Inventor after: Wang Jingsheng Inventor after: Li Yifan Inventor after: Liang Lvjie Inventor after: Huang Anguo Inventor before: Pang Shengyong Inventor before: Li Yifan Inventor before: Liang Lvjie Inventor before: Wang Jingsheng Inventor before: Huang Anguo |
|
| CB03 | Change of inventor or designer information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191220 |
|
| RJ01 | Rejection of invention patent application after publication |