CN110598358A - Stress-deformation simulation method, device, equipment and storage medium for additive manufacturing - Google Patents
Stress-deformation simulation method, device, equipment and storage medium for additive manufacturing 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
- linear heat
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 104
- 238000004088 simulation Methods 0.000 title claims abstract description 61
- 239000000654 additive Substances 0.000 title claims abstract description 50
- 230000000996 additive effect Effects 0.000 title claims abstract description 50
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 50
- 238000003860 storage Methods 0.000 title claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 66
- 238000009826 distribution Methods 0.000 claims abstract description 56
- 238000011068 loading method Methods 0.000 claims abstract description 22
- 230000004907 flux Effects 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 23
- 238000003466 welding Methods 0.000 claims description 20
- 238000010438 heat treatment Methods 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000000155 melt Substances 0.000 claims 2
- 238000012876 topography Methods 0.000 claims 2
- 238000004364 calculation method Methods 0.000 abstract description 34
- 230000035882 stress Effects 0.000 description 18
- 239000010410 layer Substances 0.000 description 10
- 230000006870 function Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000007493 shaping process Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 230000008646 thermal stress Effects 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 230000003416 augmentation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
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)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种增材制造应力变形仿真算方法、装置、设备以及计算机存储介质，方法包括：获取结构件成形过程中熔池的当前形貌；根据所述当前形貌，匹配合适的点状热源模型的参数；根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型；获取线状热源模型的参数以及时间步长，并根据所述线状热源模型的参数以及时间步长，获得结构件的能量温度场分布；根据所述能量温度场分布，以获得结构件的应力变形分布。本发明通过能量加载分布的方式进行仿真计算，能够在不影响计算结果准确性的基础上，大幅度降低计算量。
The invention discloses a stress-deformation simulation calculation method, device, equipment and computer storage medium for additive manufacturing. The method includes: obtaining the current shape of the molten pool during the forming process of a structural part; matching suitable points according to the current shape The parameters of the linear heat source model; according to the parameters of the point heat source model, the linear heat source model is constructed based on the energy loading distribution; the parameters and the time step of the linear heat source model are obtained, and according to the parameters and time of the linear heat source model The step size is to obtain the energy temperature field distribution of the structural member; according to the energy temperature field distribution, the stress deformation distribution of the structural member is obtained. The present invention performs simulation calculations in the manner of energy loading distribution, and can greatly reduce the amount of calculations without affecting the accuracy of calculation results.
Description
技术领域technical field
本发明涉及增材制造数值模拟仿真技术领域，尤其涉及一种增材制造应力变形仿真方法、装置、设备和存储介质。The present invention relates to the technical field of numerical simulation of additive manufacturing, in particular to a method, device, equipment and storage medium for stress deformation simulation of additive manufacturing.
背景技术Background technique
所谓的增材制造应力变形仿真，可以认为是通过计算机模拟整个增材制造成形过程，得到其过程的应力变形分布情况。由于成形工件的应力分布与成形过程的温度场分布联系紧密，因而在实际数值模拟过程中，温度场的计算显得格外重要。然而在增材制造过程中，热源具有集中、移动的特点，会形成在空间和时间上梯度都很大的不均匀温度场。由于其热源的高度集中性，如果直接采用传统的移动点状热源进行计算，在建立有限元模型时，需要将网格单元划分得很细，计算中也需要大量的时间步进行迭代计算，从而使得计算效率极为低下，也无法对一些结构较为复杂体型较大的构件成形过程进行模拟。The so-called stress and deformation simulation of additive manufacturing can be considered as the simulation of the entire additive manufacturing forming process by computer to obtain the stress and deformation distribution of the process. Since the stress distribution of the formed workpiece is closely related to the temperature field distribution during 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, which will form a non-uniform temperature field with a large gradient in space and time. Due to the high concentration of its heat source, if the traditional moving point heat source is directly used for calculation, the mesh unit needs to be divided very finely when establishing the finite element model, and a large number of time steps are required for iterative calculation in the calculation, so that The calculation efficiency is extremely low, and it is impossible to simulate the forming process of some components with relatively complex structures and large sizes.
发明内容Contents of the invention
针对上述问题，本发明的目的在于提供一种增材制造应力变形仿真方法、装置、设备和存储介质，本发明通过能量加载分布的方式进行仿真计算，能够在不影响计算结果准确性的基础上，大幅度的降低了计算量，解决增材制造数值模拟中计算量过大，难度高的问题。In view of the above problems, the purpose of the present invention is to provide a stress deformation simulation method, device, equipment and storage medium for additive manufacturing. The present invention performs simulation calculations in the form of energy loading distribution, which can be achieved without affecting the accuracy of calculation results. , which greatly reduces the amount of calculation, and solves the problem of excessive calculation and high difficulty in the numerical simulation of additive manufacturing.
本发明实施例提供了一种增材制造应力变形仿真方法，包括：An embodiment of the present invention provides a stress-deformation simulation method for additive manufacturing, including:
获取结构件成形过程中熔池的当前形貌；Obtain the current shape of the molten pool during the forming process of the structural part;
根据所述当前形貌，匹配合适的点状热源模型的参数；Matching parameters of a suitable point heat source model according to the current shape;
根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型；Constructing a linear heat source model based on energy loading distribution according to the parameters of the point heat source model;
获取线状热源模型的参数以及时间步长，并根据所述线状热源模型的参数以及时间步长，获得结构件的能量温度场分布；Obtaining the parameters and the time step of the linear heat source model, and obtaining the energy temperature field distribution of the structural member according to the parameters and the time step of the linear heat source model;
根据所述能量温度场分布，以获得结构件的应力变形分布。According to the energy temperature field distribution, the stress deformation distribution of the structural member is obtained.
优选地，获取熔池的当前形貌步骤之前，还包括Preferably, before the step of obtaining the current morphology of the molten pool, it also includes
获取预先建立的结构件的三维几何模型进行网络划分后的生成的有限元模型；Obtain the finite element model generated after the network division of the pre-established three-dimensional geometric model of the structural part;
根据所述有限元分析模型，获取对应的结构件的成形工艺参数；According to the finite element analysis model, the forming process parameters of the corresponding structural parts are obtained;
根据所述结构件的成形工艺参数，对结构件成形过程中熔池的参数进行初始化，以获得熔池的当前形貌。According to the forming process parameters of the structural part, the parameters of the molten pool during the forming process of the structural part are initialized to obtain the current shape 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, then according to the forming process parameters of the structural part, the parameters of the molten pool during the forming process of the structural part are initialized to obtain the The current appearance, specifically:
根据成形路径的宽度，以获得熔池的宽度；According to the width of the forming path, the width of the molten pool is obtained;
根据成形路径每层的高度，以获得熔池的深度；According to the height of each layer of the forming path, the depth of the molten pool is obtained;
根据所述熔池的宽度以及熔池的深度，经工艺成形的能量输入方式，以获取熔池的当前形貌。According to the width and depth of the molten pool, the current shape of the molten pool is obtained through the energy input mode of process shaping.
优选地，根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型，具体为：Preferably, according to the parameters of the point heat source model, a linear heat source model is constructed based on energy loading distribution, specifically:
根据热输入功率、热效率系数、焊接电压、焊接电流以及点状热源模型的参数，以获得线状热源的热流密度最大值；According to the parameters of heat input power, thermal efficiency coefficient, welding voltage, welding current and point heat source model, the maximum heat flux density of linear heat source can be obtained;
根据线状热源的热流密度最大值以及点状热源模型的参数，获得线状热源内部空间各处热流密度；According to the maximum value of the heat flux of the linear heat source and the parameters of the point heat source model, the heat flux of each place in the inner space of the linear heat source is obtained;
根据运动速度以及点状热源模型的参数，以获得线状热源的加热时间；According to the speed of movement and the parameters of the point heat source model, the heating time of the linear heat source is obtained;
根据所述线状热源的加热时间以及线状热源内部空间各处热流密度，构建线状热源模型。A model of the linear heat source is constructed according to the heating time of the linear heat source and the heat flux in each space inside the linear heat source.
优选地，所述线状热源模型表达式为：qs(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；其中，qs是热源内部空间各处热流密度，a、b分别是点状热源模型的宽度和高度，Qm是热输入功率，Qm＝ηUI，η是热效率系数，U是焊接电压，I是焊接电流，qsm是线状热源的热流密度最大值，ts是线状热源模型的加热时间。Preferably, the linear heat source model expression is: q s (x,y,z)=q sm exp(-3x 2 /a 2 )exp(-3z 2 /b 2 ); Among them, q s is the heat flux density everywhere in the inner space of the heat source, a and b are the width and height of the point heat source model respectively, Q m is the heat input power, Q m =ηUI, η is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, q sm is the maximum heat flux density of the linear heat source, and t s is the heating time of the linear heat source model.
优选地，所述线状热源模型的时间步长为点状热源模型的时间步长的10倍。Preferably, the time step of the linear heat source model is 10 times the time step of the point heat source model.
本发明实施例还提供了一种增材制造应力变形仿真装置，包括：The embodiment of the present invention also provides a stress-deformation simulation device for additive manufacturing, including:
当前形貌获取单元，用于获取结构件成形过程中熔池的当前形貌；The current shape acquisition unit is used to obtain the current shape of the molten pool during the forming process of the structural part;
参数匹配单元，用于根据所述当前形貌，匹配合适的点状热源模型的参数；A parameter matching unit, configured to match parameters of a suitable point heat source model according to the current shape;
线状热源模型构建单元，用于根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型；A linear heat source model construction unit, configured to construct a linear heat source model based on energy loading distribution according to the parameters of the point heat source model;
能量温度场分布获取单元，用于获取线状热源模型的参数以及时间步长，并根据所述线状热源模型的参数以及时间步长，获得结构件的能量温度场分布；An energy temperature field distribution acquisition unit configured to acquire parameters and a time step of the linear heat source model, and obtain the energy temperature field distribution of the structural member according to the parameters and the time step of the linear heat source model;
应力变形分布获取单元，用于根据所述能量温度场分布，以获得结构件的应力变形分布。The stress and deformation distribution acquisition unit is configured to obtain the stress and deformation distribution of the structural member according to the energy and temperature field distribution.
优选地，还包括Preferably, also include
有限元模型获取单元，用于获取预先建立的结构件的三维几何模型进行网络划分后的生成的有限元模型；The finite element model obtaining unit is used to obtain the finite element model generated after the three-dimensional geometric model of the pre-established structural part is divided into networks;
成形工艺参数获取单元，用于根据所述有限元分析模型，获取对应的结构件的成形工艺参数；The forming process parameter acquisition unit is used to acquire the forming process parameters of the corresponding structural parts according to the finite element analysis model;
初始化单元，用于根据所述结构件的成形工艺参数，对结构件成形过程中熔池的参数进行初始化，以获得熔池的当前形貌。The initialization unit is configured to initialize the parameters of the molten pool during the forming process of the structural part according to the forming process parameters of the structural part, so as to obtain the current shape 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, then according to the forming process parameters of the structural part, the parameters of the molten pool during the forming process of the structural part are initialized to obtain the The current appearance, specifically:
根据成形路径的宽度，以获得熔池的宽度；According to the width of the forming path, the width of the molten pool is obtained;
根据成形路径每层的高度，以获得熔池的深度；According to the height of each layer of the forming path, the depth of the molten pool is obtained;
根据所述熔池的宽度以及熔池的深度，经工艺成形的能量输入方式，以获取熔池的当前形貌。According to the width and depth of the molten pool, the current shape of the molten pool is obtained through the energy input mode of process shaping.
优选地，根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型，具体为：Preferably, according to the parameters of the point heat source model, a linear heat source model is constructed based on energy loading distribution, specifically:
根据热输入功率、热效率系数、焊接电压、焊接电流以及点状热源模型的参数，以获得线状热源的热流密度最大值；According to the parameters of heat input power, thermal efficiency coefficient, welding voltage, welding current and point heat source model, the maximum heat flux density of linear heat source can be obtained;
根据线状热源的热流密度最大值以及点状热源模型的参数，获得线状热源内部空间各处热流密度；According to the maximum value of the heat flux of the linear heat source and the parameters of the point heat source model, the heat flux of each place in the inner space of the linear heat source is obtained;
根据运动速度以及点状热源模型的参数，以获得线状热源的加热时间；According to the speed of movement and the parameters of the point heat source model, the heating time of the linear heat source is obtained;
根据所述线状热源的加热时间以及线状热源内部空间各处热流密度，构建线状热源模型。A model of the linear heat source is constructed according to the heating time of the linear heat source and the heat flux in each space inside the linear heat source.
优选地，所述线状热源模型表达式为：qs(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；其中，qs是热源内部空间各处热流密度，a、b分别是点状热源模型的宽度和高度，Qm是热输入功率，Qm＝ηUI，η是热效率系数，U是焊接电压，I是焊接电流，qsm是线状热源的热流密度最大值，ts是线状热源模型的加热时间。Preferably, the linear heat source model expression is: q s (x,y,z)=q sm exp(-3x 2 /a 2 )exp(-3z 2 /b 2 ); Among them, q s is the heat flux density everywhere in the inner space of the heat source, a and b are the width and height of the point heat source model respectively, Q m is the heat input power, Q m =ηUI, η is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, q sm is the maximum heat flux density of the linear heat source, and t s is the heating time of the linear heat source model.
优选地，所述线状热源模型的时间步长为点状热源模型的时间步长的10倍。Preferably, the time step of the linear heat source model is 10 times the time step of the point heat source model.
本发明实施例还提供了一种增材制造应力变形仿真方法设备，包括处理器、存储器以及存储在所述存储器内的计算机程序，所述计算机程序能够被所述处理器执行以实现如第一方面所述的增材制造应力变形仿真方法。The embodiment of the present invention also provides a stress-deformation simulation method device for additive manufacturing, including a processor, a memory, and a computer program stored in the memory, the computer program can be executed by the processor to achieve the first The stress-deformation simulation method for additive manufacturing described in the aspect.
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如上述所述的增材制造应力变形仿真方法。An embodiment of the present invention also provides a computer-readable storage medium, the computer-readable storage medium includes a stored computer program, wherein, when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the above-mentioned The stress-deformation simulation method for additive manufacturing.
实施本发明实施例，具有如下有益效果：Implementing the embodiment of the present invention has the following beneficial effects:
1、本发明根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型，通过能量加载分布的方式进行仿真计算，能够在不影响计算结果准确性的基础上，大幅度的降低了计算量，解决增材制造数值模拟中计算量过大，难度高的问题。1. According to the parameters of the point heat source model, the present invention builds a linear heat source model based on the energy loading distribution, and performs simulation calculations in the way of energy loading distribution, which can greatly reduce the temperature without affecting the accuracy of the calculation results. The amount of calculation is reduced, and the problem of excessive calculation and high difficulty in the numerical simulation of additive manufacturing is solved.
2、本发明降低了对增材制造应力仿真数值模拟网格尺寸的要求，不需要通过划分大量的致密网格来保证计算结果的准确性，由于线状分布热源在沿成形方向上的热流密度为均匀分布，所以可以将沿成形方向上的网格尺寸划分较大，大大降低了总体的网格数量。2. The present invention reduces the requirement for the numerical simulation grid size of additive manufacturing stress simulation, and does not need to divide a large number of dense grids to ensure the accuracy of the calculation results. Due to the heat flux density of the linearly distributed heat source along the forming direction For uniform distribution, the grid size along the forming direction can be divided into a larger size, which greatly reduces the overall number of grids.
3、本发明降低了对增材制造应力仿真数值模拟时间步长的要求，不需要通过划分密集的时间步长来近似拟合整个成形过程，大大降低了计算量。3. The present invention reduces the requirement on the numerical simulation time step of additive manufacturing stress simulation, and does not need to divide dense time steps to approximately fit the entire forming process, greatly reducing the amount of calculation.
附图说明Description of drawings
为了更清楚地说明本发明的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。In order to illustrate the technical solution of the present invention more clearly, the accompanying drawings used in the implementation will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some implementations of the present invention. As far as the skilled person is concerned, other drawings can also be obtained based on these drawings on the premise of not paying creative work.
图1是本发明第一实施例提供的增材制造应力变形仿真方法的流程示意图。Fig. 1 is a schematic flow chart of a stress-deformation simulation method for additive manufacturing provided by a first embodiment of the present invention.
图2是本发明第二实施例提供的线性热源模型的结构示意图。Fig. 2 is a schematic structural diagram of a linear heat source model provided by the second embodiment of the present invention.
图3是本发明第二实施例提供的增材制造应力变形仿真方法装置的结构示意图。Fig. 3 is a schematic structural diagram of an apparatus for a stress-deformation simulation method for additive manufacturing provided by a second embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.
为了更好的理解本发明的技术方案，下面结合附图对本发明实施例进行详细描述。In order to better understand the technical solutions of the present invention, the embodiments of the present invention will be described in detail below in conjunction with the accompanying drawings.
应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。It should be clear that the described embodiments are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present invention.
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。Terms used in the embodiments of the present invention are only for the purpose of describing specific embodiments, and are not intended to limit the present invention. As used in the embodiments of the present invention and the appended claims, the singular forms "a", "said" and "the" are also intended to include the plural forms unless the context clearly indicates otherwise.
应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。It should be understood that the term "and/or" used herein is only an association relationship describing associated objects, which means that there may be three relationships, for example, A and/or B, which may mean that A exists alone, and A and B exist simultaneously. B, there are three situations of B alone. In addition, the character "/" in this article generally indicates that the contextual objects are an "or" relationship.
取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。类似地，取决于语境，短语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测(陈述的条件或事件)”。Depending on the context, the word "if" as used herein may be interpreted as "at" or "when" or "in response to determining" or "in response to detecting". Similarly, depending on the context, the phrases "if determined" or "if detected (the stated condition or event)" could be interpreted as "when determined" or "in response to the determination" or "when detected (the stated condition or event) )" or "in response to detection of (a stated condition or event)".
实施例中提及的“第一\第二”仅仅是是区别类似的对象，不代表针对对象的特定排序，可以理解地，“第一\第二”在允许的情况下可以互换特定的顺序或先后次序。应该理解“第一\第二”区分的对象在适当情况下可以互换，以使这里描述的实施例能够以除了在这里图示或描述的那些以外的顺序实施。The "first\second" mentioned in the embodiment is only to distinguish similar objects, and does not represent a specific ordering of objects. It is understandable that "first\second" can be interchanged with specific sequence or sequence. It should be understood that the terms "first\second" are interchangeable under appropriate circumstances such that the embodiments described herein can be practiced in sequences other than those illustrated or described herein.
实施例一：Embodiment one:
请参阅图1以及图2，本发明第一实施例提供了一种增材制造应力变形仿真方法，其可由增材制造应力变形仿真设备来执行，特别的，由增材制造应力变形仿真设备内的一个或多个处理器来执行，并至少包括如下步骤：Please refer to Fig. 1 and Fig. 2, the first embodiment of the present invention provides a stress-deformation simulation method for additive manufacturing, which can be performed by stress-deformation simulation equipment for additive manufacturing. executed by one or more processors, and include at least the following steps:
S101，获取结构件成形过程中熔池的当前形貌。S101. Obtain the current shape of the molten pool during the forming process of the structural part.
在本实施例中，增材制造应力变形仿真设备根据所述结构件的成形工艺参数，对结构件成形过程中熔池的参数进行初始化，以获得熔池的当前形貌。其中，所述熔池的当前形貌包括熔池的熔深、熔宽以及熔池的形状等。所述成形工艺参数包括送丝速度、热源功率、热源半径、扫描速度、成形路径的宽度以及成形路径每层的高度等，这些参数能够使得模拟的环境无限接近于实际的环境。具体地，根据成形路径的宽度，以获得熔池的宽度；以及根据成形路径每层的高度，以获得熔池的深度；然后根据所述熔池的宽度以及熔池的深度，经工艺成形的能量输入方式，以获取熔池的当前形貌。In this embodiment, the additive manufacturing stress-deformation simulation device initializes the parameters of the molten pool during the forming process of the structural part according to the forming process parameters of the structural part, so as to obtain the current shape of the molten pool. Wherein, the current morphology of the molten pool includes the molten pool depth, molten width, shape of the molten pool, and the like. The forming process parameters include wire feeding speed, heat source power, heat source radius, scanning speed, width of the forming path and height of each layer of the forming path, etc. These parameters can make the simulated environment infinitely close to the actual environment. Specifically, according to the width of the forming path, the width of the molten pool is obtained; and according to the height of each layer of the forming path, the depth of the molten pool is obtained; then according to the width of the molten pool and the depth of the molten pool, the process-shaped Energy input method to obtain the current shape of the molten pool.
需要说明的是，所述成形工艺参数，根据所述有限元分析模型获取的，具体地，增材制造应力变形仿真设备获取预先建立的结构件的三维几何模型进行网络划分后的生成的有限元模型；根据所述有限元分析模型，获取对应的结构件的成形工艺参数，根据所述结构件的成形工艺参数，对结构件成形过程中熔池的参数进行初始化，以获得熔池的当前形貌。其中，所述结构件的几何模型是在三维建模软件中建立的，然后在三维建模软件中导入网格划分软件进行划分网格，需要说明的是，为了使得所述有限元模型所需网格尺寸可大幅度扩大，降低网格格数，本发明采用网格类型为六面体网格，网格密度为靠近焊缝中心处网格密集，远离焊缝处网格稀疏，相对于传统热源模型需要的网格划分来说，能够有效降低计算量。It should be noted that the forming process parameters are obtained according to the finite element analysis model, specifically, the stress-deformation simulation equipment for additive manufacturing obtains the pre-established three-dimensional geometric model of the structural part and performs network division to generate the finite element model; according to the finite element analysis model, the forming process parameters of the corresponding structural parts are obtained, and according to the forming process parameters of the structural parts, the parameters of the molten pool in the forming process of the structural parts are initialized to obtain the current shape of the molten pool appearance. Wherein, the geometric model of the structural part is established in the three-dimensional modeling software, and then the grid division software is imported into the three-dimensional modeling software to divide the grid. It should be noted that, in order to make the required finite element model The grid size can be greatly expanded to reduce the number of grid cells. The grid type used in the present invention is a hexahedral grid. The grid density is dense near the center of the weld and sparse at the distance from the weld. Compared with the traditional heat source model In terms of the required grid division, the calculation amount can be effectively reduced.
S102，根据所述当前形貌，匹配合适的点状热源模型的参数。S102. Match the parameters of a suitable point heat source model according to the current shape.
在本实施例中，由于在增材过程中，所述点状热源模型的半径比熔池的宽度的一半还要小，所以为了匹配合适的点状热源模型参数，常令点状热源模型半径的初值为熔池的宽度的一半，基于此，则根据熔池当前形貌的宽度，以获得合适的点状热源模型的半径；根据熔池当前形貌的高度，以获得合适的点状热源模型的高度。In this embodiment, since the radius of the point heat source model is smaller than half the width of the molten pool during the additive process, in order to match the appropriate point heat source model parameters, the radius of the point heat source model is often set to The initial value of is half of the width of the molten pool. Based on this, the radius of a suitable point heat source model can be obtained according to the width of the current shape of the molten pool; the appropriate point shape can be obtained according to the height of the current shape of the molten pool. The height of the heat source model.
S103，根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型。S103. Construct a linear heat source model based on energy loading distribution according to the parameters of the point heat source model.
在本实施例中，传统的热量加载方式进行模拟时，由于点状热源体积小，在成形路径部位需要划分相当致密的网格，并采用细小时间步长来对整个成形过程进行划分，因而需要大量的时间步进行计算，导致巨大的计算量，因此，本发明根据点状热源模型的参数包括点状热源模型的半径以及高度、热输入功率以及运动速度等，则根据点状热源模型的半径、点状热源模型的高度、热输入功率、热效率系数、焊接电压、焊接电流，以获得线状热源的热流密度最大值；根据线状热源的热流密度最大值、根据点状热源模型的半径以及点状热源模型的高度，获得线状热源内部空间各处热流密度；根据运动速度以及根据点状热源模型的半径，以获得线状热源的加热时间；根据所述线状热源的加热时间以及线状热源内部空间各处热流密度，构建线状热源模型，通过能量加载方式去建立线状热源模型，降低了原来点状热源模型带来的大计算量。In this embodiment, when simulating with the traditional heat loading method, due to the small size of the point heat source, it is necessary to divide a rather dense mesh in the forming path, and use a small time step to divide the entire forming process, so it is necessary A large number of time steps are calculated, resulting in a huge amount of calculation. Therefore, according to the parameters of the point heat source model in the present invention, including the radius and height of the point heat source model, heat input power and movement speed, etc., then according to the radius of the point heat source model , the height of the point heat source model, heat input power, thermal efficiency coefficient, welding voltage, and welding current to obtain the maximum heat flux density of the linear heat source; according to the maximum heat flux density of the linear heat source, according to the radius of the point heat source model and According to the height of the point heat source model, the heat flux density in the internal space of the linear heat source is obtained; according to the moving speed and the radius of the point heat source model, the heating time of the linear heat source is obtained; according to the heating time of the linear heat source and the line The heat flux density in the internal space of the heat source is used to construct a linear heat source model, and the linear heat source model is established by energy loading, which reduces the large amount of calculation brought by the original point heat source model.
其中，能量加载分布方程为：qs(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；所述线状热源模型表达式为：qs(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；qs是热源内部空间各处热流密度，a、b分别是点状热源模型的半径和高度，Qm是热输入功率，Qm＝ηUI，η是热效率系数，U是焊接电压，I是焊接电流，qsm是线状热源的热流密度最大值，ts是线状热源模型的加热时间。Wherein, the energy loading distribution equation is: q s (x,y,z)=q sm exp(-3x 2 /a 2 )exp(-3z 2 /b 2 ); the linear heat source model expression is: q s (x,y,z)=q sm exp(-3x 2 /a 2 )exp(-3z 2 /b 2 ); q s is the heat flux density everywhere in the inner space of the heat source, a and b are the radius and height of the point heat source model respectively, Q m is the heat input power, Q m =ηUI, η is the thermal efficiency coefficient, U is the welding voltage, I is the welding Current, q sm is the maximum heat flux density of the linear heat source, t s is the heating time of the linear heat source model.
S104，获取线状热源模型的参数以及时间步长，并根据所述线状热源模型的参数以及时间步长，获得结构件的能量温度场分布。S104. Obtain the parameters and time step of the linear heat source model, and obtain the energy temperature field distribution of the structural part according to the parameters and time step of the linear heat source model.
在本实施例中，所述线状热源模型的时间步长为点状热源模型的时间步长的10倍；其中，点状热源时间步长只要满足有限元模型的模拟计算条件即可，即：式中Δx为划分的有限元网格大小，cp为比热容。当然，所述线状热源模型的时间步长也可以为点状热源模型的时间步长的8倍或者7倍等，在此，本发明不再赘述。In this embodiment, the time step of the linear heat source model is 10 times the time step of the point heat source model; wherein, the time step of the point 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 grid, and c p is the specific heat capacity. Certainly, the time step of the linear heat source model may also be 8 times or 7 times of the time step of the point heat source model, etc., which will not be repeated in the present invention.
在本实施例中，由于线状热源模型参数和时间步长为求解器计算温度场所必需先设定的条件，因此，通过将所述线状热源模型的参数和时间步长输入至求解器，以计算结构件的能量温度场分布，其中，求解器为温度场计算的软件，本发明使用的是ABAQUS。In this embodiment, since the parameters of the linear heat source model and the time step are the conditions that must be set in advance for the solver to calculate the temperature location, by inputting the parameters and the time step of the linear heat source model into the solver, To calculate the energy temperature field distribution of the structural parts, the solver is the software for temperature field calculation, and the present invention uses ABAQUS.
S105，根据所述能量温度场分布，以获得结构件的应力变形分布。S105. Obtain the stress and deformation distribution of the structural member according to the energy temperature field distribution.
在本实施例中，根据所述能量温度场分布，通过热-应力耦合得到结构件的最终的应力变形分布。In this embodiment, according to the energy-temperature field distribution, the final stress-deformation distribution of the structural member is obtained through thermal-stress coupling.
综上，本发明根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型，通过能量加载分布的方式进行仿真计算，能够在不影响计算结果准确性的基础上，大幅度的降低了计算量，解决增材制造数值模拟中计算量过大，难度高的问题。同时本发明降低了对增材制造应力仿真数值模拟网格尺寸的要求，不需要通过划分大量的致密网格来保证计算结果的准确性，由于线状分布热源在沿成形方向上的热流密度为均匀分布，所以可以将沿成形方向上的网格尺寸划分较大，大大降低了总体的网格数量。以及本发明降低了对增材制造应力仿真数值模拟时间步长的要求，不需要通过划分密集的时间步长来近似拟合整个成形过程，大大降低了计算量。In summary, the present invention constructs a linear heat source model based on the energy loading distribution based on the parameters of the point heat source model, and performs simulation calculations in the manner of energy loading distribution, which can greatly increase the accuracy of the calculation results without affecting the accuracy of the calculation results. The amount of calculation is reduced, and the problem of excessive calculation and high difficulty in the numerical simulation of additive manufacturing is solved. At the same time, the present invention reduces the requirement for the grid size of the stress simulation numerical simulation of additive manufacturing, and does not need to divide a large number of dense grids to ensure the accuracy of the calculation results. Since the heat flux of the linearly distributed heat source along the forming direction is Uniform distribution, so the mesh size along the forming direction can be divided into larger, which greatly reduces the overall number of meshes. And the present invention reduces the requirement on the numerical simulation time step of additive manufacturing stress simulation, does not need to approximate the entire forming process by dividing dense time steps, and greatly reduces the amount of calculation.
为便于对本发明的理解，下面以实际的应用场景来说明本实施例的应用。In order to facilitate the understanding of the present invention, the application of this embodiment will be described below using actual application scenarios.
假设成形件为短梁件，整个短梁成型分为两个阶段：基体部分和筋条部分，其中基体的大小为563×80×20mm，分8层进行堆积。基体上面堆积部分为筋条，x方向横向宽度为12mm，y方向纵向宽度为8mm，高度为50mm。分33层按照不同的工艺路径进行堆叠。Assuming that the formed part is a short beam, the entire short beam forming is divided into two stages: the base part and the rib part, where the size of the base is 563 × 80 × 20mm, and it is stacked in 8 layers. The stacked part on the substrate is a rib, the horizontal width in the x direction is 12mm, the longitudinal width in the y direction is 8mm, and the height is 50mm. It is divided into 33 layers and stacked according to different process paths.
基体部分和筋条部分成形工艺参数如下表1和表2所示。The forming process parameters of the matrix part and the rib part are shown in Table 1 and Table 2 below.
表1：基体成形(CX)工艺基本参数Table 1: Basic parameters of matrix forming (CX) process
表2：筋条成形(JT)工艺基本参数Table 2: Basic parameters of rib forming (JT) process
具体地，在步骤S1中，建立实体三维有限元模型，首先在三维建模软件中建立工件的几何模型，然后导入网格划分软件中划分网格，采用线性六面体单元，单元尺寸可划较疏，单元总数为38962。在步骤S2中，通过上表中工艺参数判断熔池宽度约为7mm，深度为8mm，截面呈双椭球形。在步骤S3中，根据熔池宽度深度参数以及形状确定热源模型中半径宽度为5mm，深度为7mm。在步骤S4中，根据S3步骤中点状热源的参数确定线状热源模型中a＝5mm，b＝7mm，热输入功率Qm＝1900W，运动速度vm＝3.33mm/s。其中，然后基于能量加载分布构建线性热源模型，其中，能量加载分布方程为：x，y，z为坐标(mm)。在步骤S5中，将时间步长定为0.5s，输入求解器中，本项实施案例中求解器为ABAQUS软件的温度场计算模块。在步骤S6中，根据线状热源模型参数、时间步长以及各种包括边界条件等各项计算温度场必须的条件，来最终在求解器中求解出构件的温度场分布。在步骤S7中，根据S6所求的温度场分布，通过热-应力耦合得到构件最终的应力变形分布，得到最大变形为10.9mm，两端翘起变形量为5.1mm。Specifically, in step S1, a solid three-dimensional finite element model is established. First, the geometric model of the workpiece is established in the three-dimensional modeling software, and then imported into the mesh division software to divide the mesh. Linear hexahedral elements are used, and the element size can be drawn relatively sparsely. , the total number of units is 38962. In step S2, it is judged from the process parameters in the above table that the width of the molten pool is about 7mm, the depth is 8mm, and the cross section is double ellipsoidal. In step S3, according to the parameters of the width and depth of the molten pool and the shape, it is determined that the radius width in the heat source model is 5 mm, and the depth is 7 mm. In step S4, a=5mm, b=7mm, heat input power Q m =1900W, and moving speed v m =3.33mm/s are determined in the linear heat source model according to the parameters of the point heat source in step S3. Among them, a linear heat source model is then constructed based on the energy loading distribution, where the energy loading distribution equation is: x, y, z are coordinates (mm). In step S5, set the time step to 0.5s, and input it into the solver. In this implementation case, the solver is the temperature field calculation module of ABAQUS software. In step S6, according to the parameters of the linear heat source model, the time step and various conditions necessary for calculating the temperature field including boundary conditions, the temperature field distribution of the component is finally solved in the solver. In step S7, according to the temperature field distribution obtained in S6, the final stress and deformation distribution of the component is obtained through thermal-stress coupling, and the maximum deformation is 10.9 mm, and the warping deformation at both ends is 5.1 mm.
将模拟计算结果与实际成形后得到的工件进行比较，基体成形完成后(未完全冷却，温度约200℃)，最大变形约1.5mm，整体呈中间翘起的波浪状变形。筋条成形过程中，单层筋条成形工序对于工件整体变形影响较小。其中在成形横向短筋条的过程中，变形逐渐增加，单层整体变形量变化区间小于0.3mm。整体工件成形后变形趋势呈现两端翘起，中间下凹形式，两端翘起变形量为5.7mm。基板变形量较大，最大变形量为11.7mm，整体变形趋势为波浪状。实验最终测量的最大变形为10.9mm，两端翘起变形量为5.1mm，因此，变形趋势与仿真得到的相同，可以很好地实现一种加速增材制造应力仿真效率的线状热源模型的建立，能够大幅度缩短模拟计算的时间。Comparing the simulated calculation results with the workpiece obtained after actual forming, after the matrix forming is completed (not completely cooled, the temperature is about 200°C), the maximum deformation is about 1.5mm, and the overall deformation is wavy in the middle. During the rib forming process, the single-layer rib forming process has little effect on the overall deformation of the workpiece. Among them, in the process of forming short transverse ribs, the deformation gradually increases, and the variation range of the overall deformation of a single layer is less than 0.3mm. After the overall workpiece is formed, the deformation trend shows that both ends are upturned and the middle is concave, and the deformation amount of both ends is 5.7mm. The deformation of the substrate is large, the maximum deformation is 11.7mm, and the overall deformation trend is wavy. The maximum deformation finally measured in the experiment is 10.9 mm, and the warping deformation at both ends is 5.1 mm. Therefore, the deformation trend is the same as that obtained by the simulation, which can well realize a linear heat source model that accelerates the stress simulation efficiency of additive manufacturing. The establishment can greatly shorten the time of simulation calculation.
本发明第二实施例：Second embodiment of the present invention:
参见图3，本发明实施例还包括一种增材制造应力变形仿真方法装置，包括：Referring to Fig. 3, the embodiment of the present invention also includes a stress-deformation simulation method device for additive manufacturing, including:
当前形貌获取单元100，用于获取结构件成形过程中熔池的当前形貌；The current shape acquisition unit 100 is used to acquire the current shape of the molten pool during the forming process of the structural part;
参数匹配单元200，用于根据所述当前形貌，匹配合适的点状热源模型的参数；A parameter matching unit 200, configured to match parameters of an appropriate point heat source model according to the current shape;
线状热源模型构建单元300，用于根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型；A linear heat source model construction unit 300, configured to construct a linear heat source model based on energy loading distribution according to the parameters of the point heat source model;
能量温度场分布获取单元400，用于获取线状热源模型的参数以及时间步长，并根据所述线状热源模型的参数以及时间步长，获得结构件的能量温度场分布；The energy temperature field distribution acquisition unit 400 is used to acquire the parameters and time step of the linear heat source model, and obtain the energy temperature field distribution of the structural member according to the parameters and time step of the linear heat source model;
应力变形分布获取单元500，用于根据所述能量温度场分布，以获得结构件的应力变形分布。The stress and deformation distribution acquisition unit 500 is configured to obtain the stress and deformation distribution of the structural member according to the energy temperature field distribution.
在上述实施例的基础上，本发明一优选实施例中，还包括On the basis of the above embodiments, in a preferred embodiment of the present invention, it also includes
有限元模型获取单元，用于获取预先建立的结构件的三维几何模型进行网络划分后的生成的有限元模型；The finite element model obtaining unit is used to obtain the finite element model generated after the three-dimensional geometric model of the pre-established structural part is divided into networks;
成形工艺参数获取单元，用于根据所述有限元分析模型，获取对应的结构件的成形工艺参数；The forming process parameter acquisition unit is used to acquire the forming process parameters of the corresponding structural parts according to the finite element analysis model;
初始化单元，用于根据所述结构件的成形工艺参数，对结构件成形过程中熔池的参数进行初始化，以获得熔池的当前形貌。The initialization unit is configured to initialize the parameters of the molten pool during the forming process of the structural part according to the forming process parameters of the structural part, so as to obtain the current shape of the molten pool.
在上述实施例的基础上，本发明一优选实施例中，所述成形工艺参数包括成形路径的宽度以及成形路径每层的高度，则根据所述结构件的成形工艺参数，对结构件成形过程中熔池的参数进行初始化，以获得熔池的当前形貌，具体为：On the basis of the above embodiments, in a preferred embodiment of the present invention, the forming process parameters include the width of the forming path and the height of each layer of the forming path, then according to the forming process parameters of the structural part, the forming process of the structural part Initialize the parameters of the molten pool in order to obtain the current morphology of the molten pool, specifically:
根据成形路径的宽度，以获得熔池的宽度；According to the width of the forming path, the width of the molten pool is obtained;
根据成形路径每层的高度，以获得熔池的深度；According to the height of each layer of the forming path, the depth of the molten pool is obtained;
根据所述熔池的宽度以及熔池的深度，经工艺成形的能量输入方式，以获取熔池的当前形貌。According to the width and depth of the molten pool, the current shape of the molten pool is obtained through the energy input mode of process shaping.
在上述实施例的基础上，本发明一优选实施例中，根据所述点状热源模型的参数，基于能量加载分布构建线状热源模型，具体为：On the basis of the above embodiments, in a preferred embodiment of the present invention, according to the parameters of the point heat source model, a linear heat source model is constructed based on the energy loading distribution, specifically:
根据热输入功率、热效率系数、焊接电压、焊接电流以及点状热源模型的参数，以获得线状热源的热流密度最大值；According to the parameters of heat input power, thermal efficiency coefficient, welding voltage, welding current and point heat source model, the maximum heat flux density of linear heat source can be obtained;
根据线状热源的热流密度最大值以及点状热源模型的参数，获得线状热源内部空间各处热流密度；According to the maximum value of the heat flux of the linear heat source and the parameters of the point heat source model, the heat flux of each place in the inner space of the linear heat source is obtained;
根据运动速度以及点状热源模型的参数，以获得线状热源的加热时间；According to the speed of movement and the parameters of the point heat source model, the heating time of the linear heat source is obtained;
根据所述线状热源的加热时间以及线状热源内部空间各处热流密度，构建线状热源模型。A model of the linear heat source is constructed according to the heating time of the linear heat source and the heat flux in each space inside the linear heat source.
在上述实施例的基础上，本发明一优选实施例中，所述线状热源模型表达式为：qs(x,y,z)＝qsmexp(-3x2/a2)exp(-3z2/b2)；其中，qs是热源内部空间各处热流密度，a、b分别是点状热源模型的宽度和高度，Qm是热输入功率，Qm＝ηUI，η是热效率系数，U是焊接电压，I是焊接电流，qsm是线状热源的热流密度最大值，ts是线状热源模型的加热时间。On the basis of the above embodiments, in a preferred embodiment of the present invention, the linear heat source model expression is: q s (x, y, z) = q sm exp(-3x 2 /a 2 ) exp(- 3z 2 /b 2 ); Among them, q s is the heat flux density everywhere in the inner space of the heat source, a and b are the width and height of the point heat source model respectively, Q m is the heat input power, Q m =ηUI, η is the thermal efficiency coefficient, U is the welding voltage, I is the welding current, q sm is the maximum heat flux density of the linear heat source, and t s is the heating time of the linear heat source model.
在上述实施例的基础上，本发明一优选实施例中，所述线状热源模型的时间步长为点状热源模型的时间步长的10倍。On the basis of the above embodiments, in a preferred embodiment of the present invention, the time step of the linear heat source model is 10 times the time step of the point heat source model.
本发明第三实施例：The third embodiment of the present invention:
本发明第三实施例提供了一种增材制造应力变形仿真方法设备，包括处理器、存储器以及存储在所述存储器内的计算机程序，所述计算机程序能够被所述处理器执行以实现如上述所述的增材制造应力变形仿真方法。The third embodiment of the present invention provides a stress-deformation simulation method device for additive manufacturing, including a processor, a memory, and a computer program stored in the memory, and the computer program can be executed by the processor to achieve the above-mentioned The stress-deformation simulation method for additive manufacturing.
本发明第四实施例：Fourth embodiment of the present invention:
本发明第四实施例提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如上述所述的增材制造应力变形仿真方法。A fourth embodiment of the present invention provides a computer-readable storage medium, the computer-readable storage medium includes a stored computer program, wherein, when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the following steps: The stress-deformation simulation method for additive manufacturing described above.
示例性的，所述计算机程序可以被分割成一个或多个单元，所述一个或者多个单元被存储在所述存储器中，并由所述处理器执行，以完成本发明。所述一个或多个单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序在增材制造应力变形仿真方法设备中的执行过程。Exemplarily, the computer program may be divided into one or more units, and the one or more units are stored in the memory and executed by the processor to implement the present invention. The one or more units may be a series of computer program instruction segments capable of completing specific functions, and the instruction segments are used to describe the execution process of the computer program in the stress-deformation simulation method equipment for additive manufacturing.
所述增材制造应力变形仿真方法设备可包括但不仅限于处理器、存储器。本领域技术人员可以理解，所述示意图仅仅是增材制造应力变形仿真方法设备的示例，并不构成对增材制造应力变形仿真方法设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述增材制造应力变形仿真方法设备还可以包括输入输出设备、网络接入设备、总线等。The equipment for the stress-deformation simulation method of additive manufacturing may include, but not limited to, a processor and a memory. Those skilled in the art can understand that the schematic diagram is only an example of the equipment for the stress-deformation simulation method of additive manufacturing, and does not constitute a limitation on the equipment for the stress-deformation simulation method of additive manufacturing, and may include more or less components than those shown in the illustration , or combine certain components, or different components, for example, the stress-deformation simulation method for additive manufacturing. The device may also include input and output devices, network access devices, buses, and the like.
所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述增材制造应力变形仿真方法设备的控制中心，利用各种接口和线路连接整个增材制造应力变形仿真方法设备的各个部分。The so-called processor can be a central processing unit (Central Processing Unit, CPU), and can also be other general-purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf Programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc., and the control center of the stress-deformation simulation method device of the additive manufacturing uses various interfaces and lines to connect the entire stress-deformation simulation of the additive manufacturing Various parts of the method equipment.
所述存储器可用于存储所述计算机程序和/或模块，所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现所述增材制造应力变形仿真方法设备的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(SecureDigital,SD)卡，闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。The memory can be used to store the computer programs and/or modules, and the processor realizes the augmentation by running or executing the computer programs and/or modules stored in the memory and calling the data stored in the memory. Various functions of equipment for stress deformation simulation method of material manufacturing. The memory may mainly include a program storage area and a data storage area, wherein the program storage area may store an operating system, at least one application program required by a function (such as a sound playback function, an image playback function, etc.) and the like; the storage data area may store Data created based on the use of the mobile phone (such as audio data, phonebook, etc.), etc. In addition, the memory can include high-speed random access memory, and can also include non-volatile memory, such as hard disk, memory, plug-in hard disk, smart memory card (Smart Media Card, SMC), secure digital (SecureDigital, SD) card, A flash memory card (Flash Card), at least one magnetic disk storage device, flash memory device, or other volatile solid state storage devices.
其中，所述增材制造应力变形仿真方法设备集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。Wherein, if the integrated unit of the stress-deformation simulation method for additive manufacturing is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present invention realizes all or part of the processes in the methods of the above embodiments, and can also be completed by instructing related hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer When the program is executed by the processor, the steps in the above-mentioned various method embodiments can be realized. Wherein, the computer program includes computer program code, and the computer program code may be in the form of source code, object code, executable file or some intermediate form. The computer-readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer memory, a read-only memory (ROM, Read-OnlyMemory), Random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content contained in the computer-readable medium may be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable Excludes electrical carrier signals and telecommunication signals.
需说明的是，以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外，本发明提供的装置实施例附图中，模块之间的连接关系表示它们之间具有通信连接，具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。It should be noted that the device embodiments described above are only illustrative, and the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physically separated. A unit can be located in one place, or it can be distributed to multiple network units. Part or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. In addition, in the drawings of the device embodiments provided by the present invention, the connection relationship between the modules indicates that they have a communication connection, which can be specifically implemented as one or more communication buses or signal lines. It can be understood and implemented by those skilled in the art without creative effort.
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。The above description is a preferred embodiment of the present invention, it should be pointed out that for those skilled in the art, without departing from the principle of the present invention, some improvements and modifications can also be made, and these improvements and modifications are also considered Be the protection scope of the present invention.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910919933.3A CN110598358A (en) | 2019-09-26 | 2019-09-26 | Stress-deformation simulation method, device, equipment and storage medium for additive manufacturing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910919933.3A CN110598358A (en) | 2019-09-26 | 2019-09-26 | Stress-deformation simulation method, device, equipment and storage medium for additive manufacturing |
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 | Stress-deformation simulation method, device, equipment and storage medium for additive manufacturing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110598358A (en) |
Cited By (6)
| 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 |
| CN115048830A (en) * | 2022-05-17 | 2022-09-13 | 华中科技大学 | Additive manufacturing partition process simulation method |
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 (9)
| 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 | 南京理工大学 | A Design Method for Variable Density Gyroid Lattice Structure with Thermal 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 | 四川大学 | A Numerical Simulation Method of Heat Source for Laser Additive Manufacturing Process |
| CN115048830A (en) * | 2022-05-17 | 2022-09-13 | 华中科技大学 | Additive manufacturing partition process simulation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110598358A (en) | Stress-deformation simulation method, device, equipment and storage medium for additive manufacturing | |
| US10274935B2 (en) | System, method, and computer program for creating geometry-compliant lattice structures | |
| CN104102919A (en) | Image classification method capable of effectively preventing convolutional neural network from being overfit | |
| CN111444573B (en) | Ship symmetric segment model generation method and device, storage medium and terminal | |
| CN114175092A (en) | Image-Based Defect Detection in Additive Manufacturing | |
| US20230325244A1 (en) | Object processing method and apparatus, computer device, storage medium, and program product | |
| US20170323034A1 (en) | Particle etching or depositing evolutionary simulation method and device and computer readable medium | |
| CN117672437B (en) | A thermal simulation method based on phonon Monte Carlo simulation and related components | |
| CN116522660A (en) | Numerical simulation method and system for explosion problem | |
| CN115495968A (en) | Mesh division method, terminal and storage medium | |
| CN104484527A (en) | Automatic and dynamic modification method of uniformly-distributed load in topological optimization process of discrete structure | |
| CN111291450A (en) | A method, device and equipment for determining a three-dimensional model of hood springback | |
| CN115114777A (en) | Injection mold pouring simulation method, device, equipment and storage medium | |
| CN105447222B (en) | The method that technique change for integrated circuit is analyzed | |
| CN112487674B (en) | Method and device for generating models of lattice structures | |
| CN114444416A (en) | Grid aggregation method and device for fluid dynamics simulation | |
| CN118839572A (en) | Finite element grid dividing method and system suitable for calculating thermal stress of turbine component | |
| CN116911236B (en) | Numerical simulation method, device and equipment of semiconductor device and storage medium | |
| CN115034081B (en) | Foam concrete pore structure construction method, device, equipment and storage medium | |
| CN117725756A (en) | Parameter determination method and device, electronic equipment and storage medium | |
| CN111787608B (en) | Anchor node layout method and device | |
| CN115983144A (en) | Method and device for calculating equivalent roughness of fully developed wind power plant and storage medium | |
| CN115470624A (en) | Injection mold simulation parameter design method, device, equipment and storage medium | |
| CN114820991A (en) | Cross processing method and device for unstructured boundary layer grids | |
| CN109324612B (en) | Robot walking parameter optimization method, computer equipment and readable 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 | ||
| 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 |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191220 |