Alloy cast hole loose prediction method and system, electronic device and storage medium

Abstract

The present application belongs to the technical field of casting process defect prediction, and discloses an alloy casting hole loose prediction method, system, electronic equipment and storage medium. The prediction method comprises the following steps: constructing a three-dimensional model of a pouring system and initializing the pouring system; searching for a connected liquid phase area; dividing the connected liquid phase area into cells, calculating the pressure loss of metal liquid flowing through any cell, and calculating the minimum pressure loss path and the corresponding minimum pressure loss value of any cell; taking the difference between the water head and the minimum pressure loss value as the boundary criterion for searching for an isolated liquid phase area; searching for the isolated liquid phase area; according to the evolution trend of the isolated liquid phase area, calculating the shrinkage increment of the isolated liquid phase area and distributing it to the isolated liquid phase area to form a hole loose; updating the liquid phase rate of each area and going to step two. In the process of predicting the hole loose distribution, the present application considers the influence of real-time pressure loss when the metal liquid flows, can more accurately realize the prediction of the hole loose distribution, and improves the prediction accuracy.

Description

Technical Field

[0001] This invention belongs to the technical field of casting process defect prediction, and more specifically, relates to a method, system, electronic device and storage medium for predicting porosity in alloy castings. Background Technology

[0002] Traditional numerical simulation models for solidification of aluminum alloy castings assume that the solidification process is layered solidification when predicting porosity, without considering the influence of the mushy region on feeding. For alloys with a wide solidification temperature range, this may lead to completely inaccurate simulation results. Therefore, in order to accurately simulate the solidification process of castings with a wide solidification temperature range and predict porosity defects, the influence of the mushy region must be considered.

[0003] Currently, some research has been conducted by those skilled in the art, but most algorithms divide compensable and uncompensable units by setting a fixed critical solidification fraction, resulting in low prediction accuracy. Accordingly, there is a technical need in this field to develop a more accurate method for predicting porosity in alloy castings over a wide solidification temperature range. Summary of the Invention

[0004] To address the aforementioned deficiencies or improvement needs of existing technologies, this invention provides a method, system, electronic device, and storage medium for predicting porosity in alloy castings. It solves the problem of low prediction accuracy in existing methods for predicting porosity defects during the solidification process of castings with a wide solidification temperature range. Based on the pressure loss generated when the molten metal flows through the mushy zone during feeding, it can more accurately predict the distribution of porosity, thus improving prediction accuracy.

[0005] To achieve the above objectives, according to a first aspect of the present invention, a method for predicting porosity in alloy castings is provided, comprising:

[0006] S1, Construct a three-dimensional model of the gating system and initialize the gating system at the beginning of solidification;

[0007] S2, search for connected liquid phase regions;

[0008] S3, the connected liquid phase region is divided into cells, the pressure loss of the molten metal flowing through any cell is calculated by Darcy's law, and the minimum pressure loss path and the corresponding minimum pressure loss value of any cell are calculated by Floyd's algorithm.

[0009] S4, use the difference between the water head and the minimum pressure loss value as the boundary criterion for searching the isolated liquid phase region, and search for the isolated liquid phase region;

[0010] S5, determine the evolution trend of the isolated liquid phase region;

[0011] S6, calculate the shrinkage increment of the isolated liquid phase region according to the evolution trend of the isolated liquid phase region, and allocate the shrinkage increment to the isolated liquid phase region to form porosity;

[0012] S7, update the liquid phase ratio of each region, go to S2, until there are no more connected liquid phase regions.

[0013] According to the alloy casting porosity prediction method provided by the present invention, S3 calculates the minimum pressure loss path and the corresponding minimum pressure loss value for any cell, specifically including:

[0014] The shortest path between any cell and the other cells is calculated using the Floyd algorithm. The pressure loss of the cells along the shortest path is then summed to obtain the pressure loss value corresponding to the shortest path.

[0015] The minimum pressure loss value corresponding to the shortest path between any cell and the other cells is taken as the minimum pressure loss value of that cell, and the corresponding shortest path is the minimum pressure loss path.

[0016] According to the porosity prediction method for alloy castings provided by this invention, the pressure loss ΔP of the molten metal flowing through any cell in S3 is specifically as follows:

[0017]

[0018] Among them, C T ρ is the secondary dendrite arm spacing coefficient; s ρ is the density of the solid phase; l μ is the liquid density; μ is the viscosity. G is the rate of change of the temperature field; T1 and T2 are the temperatures on either side of the cell, respectively; T l T s These are the liquidus line and the solidus line, respectively.

[0019] According to the alloy casting porosity prediction method provided by the present invention, S6 specifically includes:

[0020] S61, when there is annihilation in the isolated liquid phase region, the shrinkage increment caused by the complete solidification of the parent isolated liquid phase region is calculated by the liquid phase volume of the isolated liquid phase region, i.e., the parent isolated liquid phase region, at the previous moment, and the shrinkage increment is allocated to the parent isolated liquid phase region to form porosity.

[0021] S62, when the isolated liquid phase region is in a state of contraction or splitting, the contraction increment is calculated by the difference between the liquid phase volume of the isolated liquid phase region (i.e., the parent isolated liquid phase region) at the previous moment and the liquid phase volume of at least one isolated liquid phase region (i.e., the child isolated liquid phase region) that still exists in the parent isolated liquid phase region at this moment, and the contraction increment is distributed proportionally to the child isolated liquid phase region according to the volume of the child isolated liquid phase region to form porosity.

[0022] According to the method for predicting porosity in alloy castings provided by the present invention, the shrinkage increment ΔV generated by the complete solidification of the parent isolated liquid phase region in S61 父 Specifically:

[0023]

[0024] Where, ρ s ρ is the density of the solid phase; l f is the liquid phase density; l Ω represents the liquid phase ratio of the cell. 父1 The volume of the cell representing the parent isolated liquid phase region;

[0025] The shrinkage increment in S62 is calculated by the following formula:

[0026]

[0027] Where, ΔV i Ω is the shrinkage increment allocated to the i-th sub-isolated liquid phase region. 子i Let be the volume of the cell in the i-th sub-isolated liquid phase region; n is the number of sub-isolated liquid phase regions.

[0028] According to the alloy casting porosity prediction method provided by the present invention, S6, which involves allocating the shrinkage increment to the isolated liquid phase region to form porosity, specifically includes:

[0029] The priority of the cell's shrinkage increment allocation is determined based on the difference between the cell's head and the minimum pressure loss value, with a larger difference indicating a higher priority.

[0030] The shrinkage increment within the isolated liquid phase region is allocated according to the priority of the shrinkage increment allocation of the cell.

[0031] The method for predicting porosity in alloy castings provided by this invention, which determines the priority of shrinkage increment allocation for cells based on the difference between the water head and the minimum pressure loss value of a cell, further includes:

[0032] Set pressure parameters. For any two cells, namely the first cell and the second cell, the difference between the water head and the minimum pressure loss value in the first cell is the first difference value, and the difference between the water head and the minimum pressure loss value in the second cell is the second difference value. When the first difference value is greater than the second difference value, and the difference between the first difference value and the second difference value is less than or equal to the pressure parameter, the priority of the second cell is set to be the same as the priority of the first cell.

[0033] According to a second aspect of the present invention, a system for predicting porosity in alloy castings is provided, comprising:

[0034] The initialization module is used to build a three-dimensional model of the gating system and initializes the gating system at the beginning of solidification.

[0035] The first search module is used to search for connected liquid phase regions;

[0036] The pressure loss calculation module is used to divide the connected liquid phase region into cells, calculate the pressure loss of the molten metal flowing through any cell using Darcy's law, and calculate the minimum pressure loss path and the corresponding minimum pressure loss value for any cell using the Floyd algorithm.

[0037] The second search module is used to use the difference between the water head and the minimum pressure loss value as the boundary criterion for searching the isolated liquid phase region;

[0038] The judgment module is used to judge the evolution trend of the isolated liquid phase region;

[0039] The determination module is used to calculate the shrinkage increment of the isolated liquid phase region based on the evolution trend of the isolated liquid phase region, and allocate the shrinkage increment to the formation of porosity in the isolated liquid phase region.

[0040] The update module is used to update the liquid phase ratio of each region, and then proceeds to the first search module until no connected liquid phase regions are found.

[0041] According to a third aspect of the invention, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the steps of the alloy casting porosity prediction method described above.

[0042] According to a fourth aspect of the present invention, a non-transitory computer-readable storage medium is provided, on which a computer program is stored, which, when executed by a processor, implements the steps of the alloy casting porosity prediction method described in any of the preceding claims.

[0043] In summary, compared with the prior art, the alloy casting porosity prediction method, system, electronic device, and storage medium provided by this invention offer the following advantages:

[0044] 1. Based on the pressure loss generated when molten metal flows through the mushy zone during the feeding process, a method for predicting porosity in alloy castings with a wide solidification temperature range was designed. The method considers the influence of real-time pressure loss during molten metal flow when predicting porosity distribution. Compared with existing technologies that usually divide compensable and non-compensable units based on a fixed value, this prediction method can more accurately determine the compensating ability of a unit, thereby enabling more accurate prediction of porosity distribution and improving prediction accuracy.

[0045] 2. A method for predicting porosity in alloy castings over a wide solidification temperature range has been developed. Within a limited time, based on the three-dimensional model of the gating system, corresponding material properties, and casting process parameters, the method can accurately simulate the location, size, and morphology of shrinkage cavities under this casting process, providing valuable reference for actual production. Attached Figure Description

[0046] Figure 1 This is a flowchart illustrating the porosity prediction method for alloy castings provided by the present invention.

[0047] Figure 2 This is a schematic diagram of cell priority sorting provided by the present invention;

[0048] Figure 3 This is a flowchart of the method for predicting porosity in alloy castings provided by the present invention in a specific embodiment;

[0049] Figure 4 This is a schematic diagram of the electronic device provided by the present invention. Detailed Implementation

[0050] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention. Furthermore, the technical features involved in the various embodiments of this invention described below can be combined with each other as long as they do not conflict with each other.

[0051] Please see Figure 1 This invention provides a method for predicting porosity in alloy castings, the method comprising:

[0052] S1, Construct a three-dimensional model of the casting system. At the beginning of solidification, initialize the casting system; at the initial moment, treat the entire casting system as an isolated liquid phase region.

[0053] S2, Search for connected liquid phase regions; in one specific implementation, a breadth-first search algorithm can be used to search for connected liquid phase regions;

[0054] S3, the connected liquid phase region is divided into cells, the pressure loss of the molten metal flowing through any cell is calculated by Darcy's law, and the minimum pressure loss path and the corresponding minimum pressure loss value of any cell are calculated by Floyd's algorithm.

[0055] S4, use the difference between the water head and the minimum pressure loss value as the boundary criterion for searching isolated liquid phase regions, and search for isolated liquid phase regions; the difference between the water head and the minimum pressure loss value of each cell is the condition for determining whether the isolated liquid phase region is connected or not.

[0056] S5, determine the evolution trend of the isolated liquid phase region;

[0057] S6, calculate the shrinkage increment of the isolated liquid phase region according to the evolution trend of the isolated liquid phase region, and allocate the shrinkage increment to the isolated liquid phase region to form porosity;

[0058] S7, update the liquid phase ratio of each region, go to S2, until there are no more connected liquid phase regions.

[0059] This embodiment derives the pressure loss of molten metal flowing through a cell based on Darcy's law, constructs the entire connected liquid phase region as a graph, uses the Floyd algorithm to calculate the shortest compensation pressure loss value between cells, and finally obtains the minimum pressure loss value of the cell. The difference between the water head of the cell and the minimum pressure loss value is used as a critical value to determine whether the cell can be compensated, and based on this, isolated liquid phase regions are searched, and the distribution of shrinkage cavities and porosity is obtained according to the evolution of isolated liquid phase regions.

[0060] The porosity prediction method for alloy castings provided by this invention is based on the pressure loss generated when the molten metal flows through the mushy zone during the feeding process. It is designed for predicting porosity in alloy castings with a wide solidification temperature range. In the process of predicting porosity distribution, the influence of real-time pressure loss during molten metal flow is considered. Compared with the existing technology, which usually divides compensable and non-compensable units based on a fixed value, this prediction method can more accurately determine the compensating ability of a unit, thereby achieving more accurate prediction of porosity distribution and improving prediction accuracy.

[0061] This invention realizes a method for predicting porosity in alloy castings with a wide solidification temperature range. Within a limited time, based on the three-dimensional model of the gating system, the corresponding material properties, and the casting process parameters, it can accurately simulate the location, size, and morphology of shrinkage cavities under this casting process, which has good reference value for actual production.

[0062] Furthermore, the solidified paste region is treated as a porous medium when deriving the boundary criterion.

[0063] Furthermore, S3 calculates the minimum pressure loss path and the corresponding minimum pressure loss value for any cell, specifically including:

[0064] The shortest path between any cell and the other cells is calculated using the Floyd algorithm. The pressure loss of the cells along the shortest path is then summed to obtain the pressure loss value corresponding to the shortest path.

[0065] The minimum pressure loss value corresponding to the shortest path between any cell and the other cells is taken as the minimum pressure loss value of that cell, and the corresponding shortest path is the minimum pressure loss path.

[0066] In other words, the Floyd algorithm can be used to find the shortest path between any cell and all other cells. For each shortest path, the pressure loss of the cell is calculated, and then the pressure loss values ​​on each shortest path are accumulated. Finally, the minimum pressure loss value is taken as the minimum pressure loss value.

[0067] Furthermore, the pressure loss ΔP of the molten metal flowing through any cell in S3 is specifically as follows:

[0068]

[0069] Among them, C T ρ is the secondary dendrite arm spacing coefficient; s ρ is the density of the solid phase; l μ is the liquid density; μ is the viscosity. G is the rate of change of the temperature field; T1 and T2 are the temperatures on either side of the cell, respectively; T l T s These are the liquidus line and the solidus line, respectively.

[0070] Furthermore, S6 specifically includes:

[0071] S61, when there is annihilation in the isolated liquid phase region, the shrinkage increment caused by the complete solidification of the parent isolated liquid phase region is calculated by the liquid phase volume of the isolated liquid phase region, i.e., the parent isolated liquid phase region, at the previous moment, and the shrinkage increment is allocated to the parent isolated liquid phase region to form porosity.

[0072] S62, when the isolated liquid phase region is in a state of contraction or splitting, the contraction increment is calculated by the difference between the liquid phase volume of the isolated liquid phase region (i.e., the parent isolated liquid phase region) at the previous moment and the liquid phase volume of at least one isolated liquid phase region (i.e., the child isolated liquid phase region) that still exists in the parent isolated liquid phase region at this moment, and the contraction increment is distributed proportionally to the child isolated liquid phase region according to the volume of the child isolated liquid phase region to form porosity.

[0073] The evolution trend of isolated liquid phase regions is as follows: they gradually shrink as cooling progresses, encompassing three evolutionary pathways: ① Contraction of isolated liquid phase regions, where a large parent isolated liquid phase region evolves into a small child isolated liquid phase region; ② Splitting of isolated liquid phase regions, where a large parent isolated liquid phase region evolves into multiple small child isolated liquid phase regions; ③ Annihilation of isolated liquid phase regions, where all cells of the parent isolated liquid phase region are transformed into solid phase cells. Essentially, since contraction and splitting are merely differences in the number of child isolated liquid phase regions, they can be treated in one manner.

[0074] Furthermore, the shrinkage increment ΔV resulting from the complete solidification of the isolated liquid phase region in S61 父 Specifically:

[0075]

[0076] Where, ρ s ρ is the density of the solid phase; l f is the liquid phase density; l Ω represents the liquid phase ratio of the cell. 父1 This represents the volume of the parent isolated liquid phase region cell. In the annihilation scenario, the shrinkage increment is allocated to the parent isolated liquid phase region cell from the previous time step. Within this parent isolated liquid phase region, the shrinkage increment creates porosity at the corresponding location. The above formula is the formula for calculating the shrinkage amount.

[0077] Furthermore, the shrinkage increment in S62 is calculated using the following formula:

[0078]

[0079] Where, ΔV i Ω is the shrinkage increment allocated to the i-th sub-isolated liquid phase region. 父1 The volume of the isolated liquid phase region of the parent cell; Ω 子i ρ is the volume of the cell in the i-th sub-isolated liquid phase region; n is the number of sub-isolated liquid phase regions; s ρ is the density of the solid phase; l f is the liquid phase density; l This represents the liquid phase ratio of a cell. In the case of contraction or splitting, the contraction increment is allocated to the sub-isolated liquid phase regions formed after the evolution at that moment. Within each sub-isolated liquid phase region, the contraction increment creates porosity at the corresponding location. The above formula is for calculating the contraction increment allocated to any sub-isolated liquid phase region. Specifically, in the case of contraction, the number of sub-isolated liquid phase regions, n, is 1, and all the contraction increment is allocated to this sub-isolated liquid phase region. In the case of splitting, the contraction increment is allocated according to the volume ratio of each sub-isolated liquid phase region; that is, the proportion of the contraction increment allocated to each sub-isolated liquid phase region is consistent with the volume ratio of each sub-isolated liquid phase region.

[0080] Furthermore, in S6, allocating the shrinkage increment to the isolated liquid phase region to form porous structures specifically includes:

[0081] The priority of the cell's shrinkage increment allocation is determined based on the difference between the cell's head and the minimum pressure loss value, with a larger difference indicating a higher priority.

[0082] The shrinkage increment within the isolated liquid phase region is allocated according to the priority of the shrinkage increment allocation of the cell.

[0083] In this embodiment, when shrinkage increments form porosity within an isolated liquid phase region, the location of porosity is allocated according to the priority order of the cells. Specifically, shrinkage increments are allocated first to cells with higher priority, forming porosity. This strategy determines the priority of shrinkage increment allocation based on the difference between the cell's head and minimum pressure loss value. A larger difference indicates a smaller minimum pressure loss value, making it easier for the flow to compensate for the shrinkage in other cells, thus facilitating porosity formation. This allocation strategy considers the real-time pressure loss during molten metal flow, enabling a more accurate assessment of the cell's compensability and resulting in more precise prediction of porosity distribution.

[0084] Furthermore, determining the priority of cell shrinkage increment allocation based on the difference between the cell's head and the minimum pressure loss value also includes:

[0085] Set pressure parameters. For any two cells, namely the first cell and the second cell, the difference between the water head and the minimum pressure loss value in the first cell is the first difference value, and the difference between the water head and the minimum pressure loss value in the second cell is the second difference value. When the first difference value is greater than the second difference value, and the difference between the first difference value and the second difference value is less than or equal to the pressure parameter, the priority of the second cell is set to be the same as the priority of the first cell.

[0086] Specifically, when allocating shrinkage increments within an isolated liquid phase region, the priority of shrinkage increment allocation is determined by the difference between the water head and the minimum pressure loss value of a cell; the larger the difference, the higher the priority. This method also introduces a pressure parameter, Critical_DeltaP, which represents the pressure increment caused by the volume force over a cell length. Under gravity casting, this pressure parameter is ρgh, where h represents the cell height. When the difference between the value corresponding to a cell of the previous priority and the value corresponding to a cell of the next priority is less than or equal to this parameter, that cell can be merged into the previous priority for shrinkage increment allocation calculation. This design facilitates faster allocation of shrinkage increments, thereby improving the efficiency of porosity prediction.

[0087] refer to Figure 2 A diagram illustrating the division of cells into different priorities, such as... Figure 2 As shown, in this embodiment, A1-A4 are priority 1 cells, B1-B6 are priority 2 cells, and then the priorities of CDEF decrease sequentially. Thus, when allocating the shrinkage increment, it is preferentially allocated to A to form a loose hole, and then allocated to B in sequence according to the priority order.

[0088] Furthermore, the present invention also provides an alloy casting porosity prediction system, which is used to implement the alloy casting porosity prediction method described in any of the above embodiments. This prediction system and the above prediction method can be understood in correspondence. The prediction system includes:

[0089] The initialization module is used to build a three-dimensional model of the gating system and initializes the gating system at the beginning of solidification.

[0090] The first search module is used to search for connected liquid phase regions;

[0091] The pressure loss calculation module is used to divide the connected liquid phase region into cells, calculate the pressure loss of the molten metal flowing through any cell using Darcy's law, and calculate the minimum pressure loss path and the corresponding minimum pressure loss value for any cell using the Floyd algorithm.

[0092] The second search module is used to use the difference between the water head and the minimum pressure loss value as the boundary criterion for searching the isolated liquid phase region;

[0093] The judgment module is used to judge the evolution trend of the isolated liquid phase region;

[0094] The determination module is used to calculate the shrinkage increment of the isolated liquid phase region based on the evolution trend of the isolated liquid phase region, and allocate the shrinkage increment to the formation of porosity in the isolated liquid phase region.

[0095] The update module is used to update the liquid phase ratio of each region, and then proceeds to the first search module until no connected liquid phase regions are found.

[0096] Furthermore, considering that existing methods for predicting porosity in casting processes do not take into account the impact of pressure loss during feeding, resulting in low prediction accuracy, this embodiment provides a method for predicting porosity in alloy castings over a wide solidification temperature range based on pressure loss. (Refer to...) Figure 3 The gravity casting shrinkage cavity defect prediction method mainly includes the following steps:

[0097] S1, construct a three-dimensional model of the casting system. At the actual moment, initialize the casting system, set the minimum pressure loss value of the entire casting system to 0, and regard the entire casting system as an isolated liquid phase region.

[0098] S2, search for connected liquid phase regions; if a connected liquid phase region exists, proceed to S3; otherwise, end; in this implementation scheme, the search for connected liquid phase regions is implemented using a breadth-first search algorithm. Specifically, the temperature field of the casting system can be solved, the solid-liquid ratio of each cell can be calculated based on the temperature field, and then the connected liquid phase regions can be determined.

[0099] S3. Search for the minimum pressure loss path within each connected liquid phase region and calculate the minimum pressure loss value at each location within all connected liquid phase regions. Specifically, the connected liquid phase region is considered as a whole composed of multiple cells. The pressure loss of the molten metal flowing through a cell is calculated using Darcy's law. The connected liquid phase regions are constructed into a graph, with cells as nodes of the graph. The pressure loss of the molten metal flowing through a cell is the weight between nodes. The shortest path between each cell is calculated using the Floyd algorithm, and the minimum pressure loss value of each cell is calculated based on this.

[0100] S4 uses the difference between the head and the minimum pressure loss value being less than zero as the boundary criterion for searching isolated liquid phase regions. Then, a breadth-first search algorithm is used to search for the distribution of isolated liquid phase regions at a certain moment. Specifically, it can first be determined whether each cell inside the connected liquid phase region can be shrunk. Once the non-shrunk cells are determined (i.e., the difference between the head and the minimum pressure loss value is less than zero), the connected liquid phase regions are searched after excluding the non-shrunk cells as isolated liquid phase regions. That is, the boundary of the isolated liquid phase region is the non-shrunk cell, and the interior is the connected liquid phase region.

[0101] S5, determine whether there is an isolated liquid phase region annihilation, that is, a certain region was an isolated liquid phase region in the previous moment and the region has been completely solidified at this moment. If it exists, go to S61; otherwise, go to S62.

[0102] S61, the shrinkage amount generated at one time step is allocated to the isolated liquid phase region at the previous time step to form shrinkage cavities; specifically, the shrinkage amount generated by the complete solidification of a certain isolated liquid phase region (called the parent isolated liquid phase region) is calculated based on the liquid phase volume at the previous time step, and the shrinkage amount is allocated to the corresponding position of the parent isolated liquid phase region to produce porosity.

[0103] S62, the shrinkage amount generated at a given time step is proportionally distributed to the volume of the isolated liquid phase region at the current time step to form shrinkage cavities; specifically, the shrinkage amount generated at a given time step is calculated by the difference between the liquid phase volume of the parent isolated liquid phase region and the liquid phase volume of the isolated liquid phase region (called the child isolated liquid phase region) that still exists in this region at the current time step, and then distributed to these child isolated liquid phase regions proportionally according to their volumes to form loose cavities;

[0104] Specifically, when porosity forms at the corresponding position in the isolated liquid phase region, the shrinkage amount is allocated to the corresponding cell according to the shrinkage rate allocation criterion, i.e., the cell allocation priority sorting, to form porosity.

[0105] S7, then update the liquid phase ratio of each region, switch to S2, until there are no more connected liquid phase regions.

[0106] Specifically, the gas phase fraction can be calculated based on the allocated shrinkage amount, and then the liquid phase fraction can be calculated based on the equilibrium equation "solid phase fraction + liquid phase fraction + shrinkage fraction = 100%". After updating the liquid phase fraction of all regions, the time step is increased to the next moment and the cycle is repeated until all isolated liquid phase regions disappear, so as to obtain the final shrinkage distribution.

[0107] Furthermore, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of the method for predicting porosity in alloy castings as described in any of the above embodiments.

[0108] Furthermore, a non-transitory computer-readable storage medium stores a computer program thereon, which, when executed by a processor, implements the steps of the alloy casting porosity prediction method as described in any of the above embodiments.

[0109] Figure 4 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 4 As shown, the electronic device may include a processor, a communications interface, a memory, and a communication bus, wherein the processor, communications interface, and memory communicate with each other via the communication bus. The processor can call logical instructions in the memory to execute the method for predicting porosity in alloy castings.

[0110] Furthermore, the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, and can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0111] On the other hand, the present invention also provides a computer program product, the computer program product comprising a computer program stored on a non-transitory computer-readable storage medium, the computer program comprising program instructions, and when the program instructions are executed by a computer, the computer is able to execute the alloy casting porosity prediction method provided by the above methods.

[0112] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to perform the alloy casting porosity prediction method provided by the methods described above.

[0113] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0114] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0115] Those skilled in the art will readily understand that the above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for predicting porosity in alloy castings, characterized in that, include: S1, Construct a three-dimensional model of the gating system and initialize the gating system at the beginning of solidification; S2, search for connected liquid phase regions; S3, the connected liquid phase region is divided into cells, the pressure loss of the molten metal flowing through any cell is calculated by Darcy's law, and the minimum pressure loss path and the corresponding minimum pressure loss value of any cell are calculated by Floyd's algorithm. S4, use the difference between the water head and the minimum pressure loss value as the boundary criterion for searching the isolated liquid phase region, and search for the isolated liquid phase region; S5, determine the evolution trend of the isolated liquid phase region; S6, calculate the shrinkage increment of the isolated liquid phase region according to the evolution trend of the isolated liquid phase region, and allocate the shrinkage increment to the isolated liquid phase region to form porosity; S7, update the liquid phase fraction of each region, go to S2, until there are no more connected liquid phase regions; S3 calculates the minimum pressure loss path and the corresponding minimum pressure loss value for any cell, specifically including: The shortest path between any cell and the other cells is calculated using the Floyd algorithm. The pressure loss of the cells along the shortest path is then summed to obtain the pressure loss value corresponding to the shortest path. The minimum pressure loss value corresponding to the shortest path between any cell and the other cells is taken as the minimum pressure loss value of that cell, and the corresponding shortest path is the minimum pressure loss path. Pressure loss of molten metal flowing through any cell in S3 Specifically: ； Among them, C T ρ is the secondary dendrite arm spacing coefficient; s ρ is the density of the solid phase; l μ is the liquid density; μ is the viscosity. G is the rate of change of the temperature field; T1 and T2 are the temperatures on either side of the cell, respectively; T l T s These are the liquidus line and the solidus line, respectively.

2. The method for predicting porosity in alloy castings as described in claim 1, characterized in that, S6 specifically includes: S61, when there is annihilation in the isolated liquid phase region, the shrinkage increment caused by the complete solidification of the parent isolated liquid phase region is calculated by the liquid phase volume of the isolated liquid phase region, i.e., the parent isolated liquid phase region, at the previous moment, and the shrinkage increment is allocated to the parent isolated liquid phase region to form porosity. S62, when the isolated liquid phase region is in a state of contraction or splitting, the contraction increment is calculated by the difference between the liquid phase volume of the isolated liquid phase region (i.e., the parent isolated liquid phase region) at the previous moment and the liquid phase volume of at least one isolated liquid phase region (i.e., the child isolated liquid phase region) that still exists in the parent isolated liquid phase region at this moment, and the contraction increment is distributed proportionally to the child isolated liquid phase region according to the volume of the child isolated liquid phase region to form porosity.

3. The method for predicting porosity in alloy castings as described in claim 2, characterized in that, The shrinkage increment caused by the complete solidification of the isolated liquid phase region in S61 Specifically: Where, ρ s ρ is the density of the solid phase; l The density of the liquid phase; The liquid phase ratio of the cell; The volume of the cell representing the parent isolated liquid phase region; The shrinkage increment in S62 is calculated by the following formula: in, The shrinkage increment allocated to the i-th sub-isolated liquid phase region, Let be the volume of the cell in the i-th sub-isolated liquid phase region; n is the number of sub-isolated liquid phase regions.

4. The method for predicting porosity in alloy castings as described in claim 1, characterized in that, S6 specifically includes the following: Distributing the shrinkage increment to the isolated liquid phase region to form porous structures. The priority of the cell's shrinkage increment allocation is determined based on the difference between the cell's head and the minimum pressure loss value, with a larger difference indicating a higher priority. The shrinkage increment within the isolated liquid phase region is allocated according to the priority of the shrinkage increment allocation of the cell.

5. The method for predicting porosity in alloy castings as described in claim 4, characterized in that, The priority of allocating cell shrinkage increments is determined based on the difference between the cell's head and the minimum pressure loss value, and also includes: Set pressure parameters. For any two cells, namely the first cell and the second cell, the difference between the water head and the minimum pressure loss value in the first cell is the first difference value, and the difference between the water head and the minimum pressure loss value in the second cell is the second difference value. When the first difference value is greater than the second difference value, and the difference between the first difference value and the second difference value is less than or equal to the pressure parameter, the priority of the second cell is set to be the same as the priority of the first cell.

6. A system for predicting porosity in alloy castings, characterized in that, The method for predicting porosity in alloy castings according to any one of claims 1-5 includes: The initialization module is used to build a three-dimensional model of the gating system and initializes the gating system at the beginning of solidification. The first search module is used to search for connected liquid phase regions; The pressure loss calculation module is used to divide the connected liquid phase region into cells, calculate the pressure loss of the molten metal flowing through any cell using Darcy's law, and calculate the minimum pressure loss path and the corresponding minimum pressure loss value for any cell using the Floyd algorithm. The second search module is used to use the difference between the water head and the minimum pressure loss value as the boundary criterion for searching the isolated liquid phase region; The judgment module is used to judge the evolution trend of the isolated liquid phase region; The determination module is used to calculate the shrinkage increment of the isolated liquid phase region based on the evolution trend of the isolated liquid phase region, and allocate the shrinkage increment to the formation of porosity in the isolated liquid phase region. The update module is used to update the liquid phase ratio of each region, and then proceeds to the first search module until no connected liquid phase regions are found.

7. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps of the method for predicting porosity in alloy castings as described in any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for predicting porosity in alloy castings as described in any one of claims 1 to 5.

Images