Numerical simulation method of steel plate corrosion based on coupling of three-dimensional cellular automata and finite element method

Abstract

The present invention relates to the field of numerical simulation technology, and specifically to a numerical simulation method for steel plate corrosion coupled with a three-dimensional cellular automation and a finite element. The present invention analyzes the initial corrosion rate and initial corrosion degree of each cellular grid area at different times through initial simulation data, sets the unit element size of the finite element grid decomposition, divides the cellular grid area into multiple unit element areas, and obtains fine simulation data of each unit element area, thereby obtaining adjusted simulation data of each cellular grid area at different times. Then, by analyzing the difference between the adjusted simulation data and the initial simulation data, the adjusted corrosion rate and adjusted corrosion degree of each cellular grid area at different times are obtained, and finally the steel plate corrosion process is dynamically displayed. By establishing a "CA-FRA bidirectional dynamic coupling" mechanism, the present invention can achieve accurate steel plate corrosion performance at different times during the steel plate corrosion simulation process, and obtain an accurate and excellent steel plate corrosion simulation process.

Description

Technical Field

[0001] The present invention relates to the technical field of numerical simulation, and in particular to a numerical simulation method for steel plate corrosion by coupling a three-dimensional cellular automaton with a finite element. Background Art

[0002] After heat treatment and pre-stretching, the mechanical properties and corrosion resistance of steel plates are significantly improved, and therefore they are widely used in coastal construction, shipbuilding, marine engineering and other fields. However, the service environment of such structures is extremely complex, and they often face harsh working conditions such as high salt, high humidity, and high temperature. They are not only subjected to continuous erosion by electrochemical corrosion, but also need to withstand the multi-field coupling of external loads. Among them, pitting is a common form of localized corrosion, which often forms deep and narrow pits on the surface of the structure. These pits not only directly weaken the effective load-bearing cross-section of the material, but also act as a source of stress concentration, inducing and accelerating the initiation and propagation of cracks under mechanical loads, significantly reducing the fatigue life and load-bearing capacity of the structure, and even causing catastrophic failure.

[0003] Traditional steel plate corrosion simulation research methods often use finite element decomposition to help obtain temperature and stress data for local areas of the steel plate. However, during the steel plate corrosion process, the temperature and stress data at different locations on the steel plate are different. When using units of the same size to obtain temperature and other data at different steel plate locations, it is difficult to reflect the local temperature changes caused by the corrosion reaction. If smaller finite element decomposition units are used, more accurate regional temperature, stress and other data can be obtained, but at the same time, the efficiency of system analysis will be affected. Therefore, traditional methods have limitations when dealing with steel plate corrosion simulations and are unable to fully and efficiently capture the complex changes in the corrosion process. Summary of the Invention

[0004] In order to solve the above technical problems, the purpose of the present invention is to provide a numerical simulation method for steel plate corrosion by coupling three-dimensional cellular automata with finite elements.

[0005] According to a first aspect of an embodiment of the present invention, a numerical simulation method for steel plate corrosion is provided by coupling a three-dimensional cellular automation with a finite element method. The technical solution adopted is as follows:

[0006] Based on the 3D cellular automation module and finite element analysis module, the initial simulation data of each cellular grid area at different times during the steel plate corrosion simulation process are collected;

[0007] Based on the initial simulation data, analyzing the initial corrosion rate and initial corrosion degree of each cellular grid area at different times, setting the unit element size of the finite element mesh decomposition, and dividing each cellular grid area at different times into a number of unit element areas;

[0008] Based on the three-dimensional cellular automaton module and the finite element analysis module, the detailed simulation data of each unit element area at different times are obtained;

[0009] According to the refined simulation data, combined with the position distribution of the unit cell area, the adjusted simulation data of each cellular grid area at different times are obtained;

[0010] Analyzing the difference between the adjusted simulation data and the initial simulation data, and combining the initial corrosion rate of the cellular grid area, obtaining the adjusted corrosion rate and the adjusted corrosion degree of each cellular grid area at different times;

[0011] Based on the adjusted corrosion rate and the adjusted corrosion degree during the steel plate corrosion simulation process, the steel plate corrosion process is dynamically displayed.

[0012] In some embodiments of the present invention, the initial simulation data includes initial temperature data and initial stress data.

[0013] In some embodiments of the present invention, analyzing the initial corrosion rate of each cellular grid region at different times based on the initial simulation data includes:

[0014] Count the number of eroded cellular grids in the upward or downward direction of each cellular grid area at different times;

[0015] The difference in the number of corroded cellular grids in each cellular grid area at different times and the previous time is analyzed, and the initial corrosion rate of each cellular grid area at different times is obtained by combining the initial temperature data of each cellular grid area at different times.

[0016] In some embodiments of the present invention, analyzing the initial corrosion degree of each cellular grid region at different times based on the initial simulation data includes:

[0017] Analyzing the corrosion performance of each cellular grid area at different times according to the initial corrosion rate and the number of corroded cellular grids;

[0018] According to the corrosion performance, combined with the initial stress data of each cellular grid area at different times, the initial corrosion degree of each cellular grid area at different times is obtained.

[0019] In some embodiments of the present invention, the unit element size of the finite element mesh decomposition is set to divide each cellular mesh region into a number of unit element regions, including:

[0020] According to the initial corrosion degree, the unit element size of the finite element mesh decomposition is set;

[0021] Each cellular grid region is divided into a number of unit cell regions according to the unit cell size, and the number of unit cell regions in each cellular grid region is obtained.

[0022] In some embodiments of the present invention, obtaining adjusted simulation data of each cellular grid region at different times based on the refined simulation data and in combination with the position distribution of the unit cell region includes:

[0023] Obtain the distance between each unit cell area and the corrosion center in each cellular grid area at different times;

[0024] According to the distance, the distance weights of each unit cell area in each cellular grid area at different times are obtained;

[0025] According to the refined simulation data and in combination with the distance weight, the adjusted simulation data of each cellular grid area at different times are obtained.

[0026] In some embodiments of the present invention, adjusting simulation data includes adjusting temperature data and adjusting stress data.

[0027] In some embodiments of the present invention, analyzing the difference between the adjusted simulation data and the initial simulation data, and combining the initial corrosion rate of the cellular grid area, obtaining the adjusted corrosion rate and the adjusted corrosion degree of each cellular grid area at different times, includes:

[0028] Calculating a ratio of the adjusted temperature data to the initial temperature data, and combining the initial corrosion rate of the cellular grid area to obtain an adjusted corrosion rate of each cellular grid area at different times;

[0029] According to the adjusted corrosion speed, the adjusted corrosion degree of each cellular grid area at different times is obtained.

[0030] In some embodiments of the present invention, obtaining the adjusted corrosion degree of each cellular grid region at different times according to the adjusted corrosion rate includes:

[0031] Count the number of eroded cellular grids in the upward or downward direction of each cellular grid area at different times;

[0032] According to the adjusted corrosion rate and the number of corroded cellular grids, analyzing the adjusted corrosion performance of each cellular grid area at different times;

[0033] According to the adjusted corrosion performance, combined with the adjusted stress data of each cellular grid area at different times, the adjusted corrosion degree of each cellular grid area at different times is obtained.

[0034] In some embodiments of the present invention, the number of eroded cellular grids in the upward or downward direction of each cellular grid region at different times is counted, which also includes:

[0035] Obtaining the cellular state of the cellular grid region based on the cellular state transition rule, wherein the cellular state includes: corrosion completed, corrosion in progress, and uncorroded;

[0036] Then, the number of corroded cellular grids in the upper or lower direction of each corroded cellular grid area at different times is counted.

[0037] Compared with the existing technology, the numerical simulation method of steel plate corrosion coupled with three-dimensional cellular automation and finite element method provided by the present invention has the following beneficial effects:

[0038] The present invention simulates steel plate corrosion by coupling a three-dimensional cellular automaton with a finite element method, solving the problem that traditional experimental research methods cannot reflect the dynamic development of corrosion damage and changes in environmental conditions during the corrosion process.

[0039] The present invention obtains initial simulation data, analyzes the initial corrosion rate and initial corrosion degree of each cellular grid area at different times, sets the unit element size of the finite element grid decomposition, divides each cellular grid area at different times into a number of unit element areas, obtains detailed simulation data of each unit element area at different times, and then combines the position distribution of the unit element areas to obtain adjusted simulation data of each cellular grid area at different times. Then, the adjusted corrosion rate and adjusted corrosion degree of each cellular grid area are analyzed to dynamically display the steel plate corrosion process. That is, the present invention establishes a "three-dimensional cellular automaton and finite element bidirectional dynamic coupling" mechanism, obtains the corrosion performance of the steel plate corrosion simulation process through the three-dimensional cellular automaton, intelligently controls the unit cloud size of the finite element decomposition to obtain more refined temperature and other data performance, and then reacts to analyze the corrosion performance of the steel plate based on the obtained temperature and other data, thereby achieving precise steel plate corrosion performance at different times during the steel plate corrosion simulation process, and obtaining a more accurate and excellent steel plate corrosion simulation process. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the technical solutions and advantages of the embodiments of the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0041] Figure 1 A schematic diagram of the basic process of a numerical simulation method for steel plate corrosion coupled with a three-dimensional cellular automation and finite element method provided by one embodiment of the present invention;

[0042] Figure 2 A schematic diagram of a cellular grid area division provided by an embodiment of the present invention;

[0043] Figure 3 A schematic longitudinal cross-sectional view of a corrosion position of a steel plate provided by one embodiment of the present invention;

[0044] Figure 4 A schematic diagram of the corrosion results of a steel plate corrosion numerical simulation result on a vertical cross section at 0.06s provided by one embodiment of the present invention;

[0045] Figure 5 A schematic diagram of the corrosion results of a steel plate corrosion numerical simulation result on a vertical cross section at 0.08s provided by one embodiment of the present invention;

[0046] Figure 6 A schematic diagram of the corrosion results of a steel plate corrosion numerical simulation result on a vertical cross section at 0.10s provided by one embodiment of the present invention;

[0047] Figure 7 A schematic diagram of the corrosion results of a steel plate corrosion numerical simulation result on a vertical cross section at 0.12s provided by one embodiment of the present invention;

[0048] Figure 8 A schematic diagram of the corrosion results of a steel plate corrosion numerical simulation result on a vertical cross section at 0.14s provided by one embodiment of the present invention;

[0049] Figure 9 A schematic diagram of the corrosion results of a steel plate corrosion numerical simulation result on a vertical cross section at 0.16s provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0050] To further illustrate the technical means and effectiveness of the present invention in achieving its intended purpose, the following, in conjunction with the accompanying drawings and preferred embodiments, describes in detail the specific implementation, structure, features, and effectiveness of the numerical simulation method for steel plate corrosion using a three-dimensional cellular automation coupled with finite element methods proposed in the present invention. In the following description, references to different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, specific features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.

[0051] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the art to which this invention pertains. Terms such as "comprises," "comprising," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a circuit structure, article, or device comprising a series of elements includes not only those elements but also other elements not explicitly listed, or elements inherent to such article or device. In the absence of further limitations, the phrase "comprising a ..." to define an element does not preclude the presence of other identical elements in the article or device comprising the element.

[0052] The specific scheme of the steel plate corrosion numerical simulation method coupled with three-dimensional cellular automation and finite element method provided by the present invention is described in detail below with reference to the accompanying drawings.

[0053] See also Figure 1 , which shows the basic process of a numerical simulation method for steel plate corrosion coupled with a three-dimensional cellular automation and finite element method provided by an embodiment of the present invention.

[0054] like Figure 1 As shown, a numerical simulation method for steel plate corrosion coupled with a three-dimensional cellular automation and a finite element method according to an embodiment of the present invention specifically includes:

[0055] S100: Establish a three-dimensional cellular automation and finite element coupling mechanism to collect initial simulation data of each cellular grid area at different times during the steel plate corrosion simulation process.

[0056] This invention primarily simulates the corrosion of steel plates in extremely complex environments. The simulation and analysis of steel plate corrosion primarily involves two modules: a 3D cellular automaton module and a finite element analysis module. The 3D cellular automaton module simulates the formation and growth of corrosion pits, while the finite element analysis module analyzes the stress and temperature distribution in the area surrounding the pits.

[0057] Therefore, in the embodiment of the present invention, by establishing a three-dimensional cellular automaton and finite element coupling mechanism, the volume of the steel plate is discretized into a three-dimensional cubic cellular grid (with a size of micrometers) based on the three-dimensional cellular automaton, and the cellular grid area of ​​the steel plate (such as Figure 2 As shown in the figure, each cellular mesh region discretized from the steel plate volume is of the same size. Finite element analysis is used to collect initial simulation data for each cellular mesh region at different times during the steel plate corrosion simulation. The initial simulation data includes initial temperature data and initial stress data.

[0058] S200: Based on the initial simulation data, the initial corrosion rate and initial corrosion degree of each cellular grid area at different times are analyzed, the unit element size of the finite element grid decomposition is set, and each cellular grid area at different times is divided into a number of unit element areas.

[0059] Traditional methods for simulating steel corrosion suffer from a disconnect between the dynamic development of corrosion damage and material property analysis. However, as the corrosion progresses, the corrosion level within the corroded area gradually reaches saturation, while the corrosion level outside the area continues to increase. As the corrosion reaction progresses, the temperature in the local area gradually rises and changes, further promoting the corrosion reaction. Therefore, a 3D cellular automaton module simulates the electrochemical corrosion growth process of the steel plate based on the iterative performance of the cellular meshes divided into sections, thereby assisting the finite element analysis module in performing unit element division. The finite element analysis module intelligently adjusts the corrosion growth performance of the steel plate during the corrosion simulation by analyzing the stress-stability distribution of different unit elements. This establishes a "bidirectional dynamic coupling mechanism between 3D cellular automaton and finite element analysis (CA-FRA)" to facilitate the simulation of the steel plate corrosion and enhance the accuracy of the data generated during the simulation. It is important to note that the interaction between the 3D cellular automaton module and the finite element analysis module occurs simultaneously during the simulated corrosion process.

[0060] To establish the "CA-FRA bidirectional dynamic coupling" mechanism, we first need to use three-dimensional cellular automata to analyze the corrosion process of steel plates during the simulation. In the process of steel plate corrosion simulation, the corrosion phenomenon often starts from a local position on the steel plate and spreads from point to surface (such as Figure 3 The steel plate's ambient temperature and other conditions also have a certain impact on its corrosion resistance, and the extent of the steel plate's corrosion reaction also affects local temperature changes. For areas with more intense corrosion, a more detailed analysis is required to obtain data such as temperature.

[0061] Therefore, it is necessary to first analyze the degree of corrosion in each cellular grid area at different times during the steel plate corrosion simulation process, so as to adjust the accuracy of the finite element analysis, thereby ensuring the accuracy of the simulation process while reducing the amount of calculation of the system to a certain extent.

[0062] Based on the above analysis, in an embodiment of the present invention, based on the initial simulation data, the initial corrosion rate and initial corrosion degree of each cellular grid area at different times are analyzed, the unit element size of the finite element mesh decomposition is set, and each cellular grid area at different times is divided into several unit element areas.

[0063] Since the volume of the steel plate is discretized into a three-dimensional solid, a three-dimensional space system is constructed, and the centroid position of each cellular grid area is used to represent the spatial position of the cellular grid area. In the three-dimensional cube cellular network, the cellular grid area is randomly selected as the initial corrosion core according to the environmental parameters (temperature, pH value, etc.) to perform steel plate corrosion simulation. During the steel plate corrosion simulation process, the corrosion performance of the steel plate is different at different times. For the location area with strong corrosion reaction, the temperature change caused by the corrosion reaction is more drastic, and a more detailed analysis of its temperature and other parameter distribution is required. Therefore, the first analysis is conducted. The corrosion rate and corrosion degree of the steel plate during the corrosion simulation at the time. Further:

[0064] Based on the initial simulation data, the initial corrosion rate of each cellular grid area at different times is analyzed, including:

[0065] First, the cell state transition rule is defined based on the principle of electrochemical kinetics, and then the cell state of the cell grid area is obtained based on the cell state transition rule. The cell state includes: corrosion end, corrosion, and non-corrosion. Then, the three-dimensional space system is judged. The cellular state of the cellular grid area at the position. For the cellular grid area at the position where the corrosion is complete and the corrosion is not complete, there is no need to analyze the corrosion rate of the cellular grid area at that position. This is because when the steel plate at a certain position is completely corroded, the steel plate at this position is empty and has no meaning to be analyzed; and the uncorroded position also has no meaning to be analyzed.

[0066] Then, the number of eroded cell grids in the up or down direction of each eroded cell grid area at different times is counted. Specifically, in the three-dimensional space system, we can obtain At the moment, the coordinate Location The number of corrosion cell grid areas on the axis, but it should be noted that due to the three-dimensional shape of the steel plate, corrosion problems may occur on the upper and lower sides of the same position of the steel plate at the same time. Therefore, the At the moment, the coordinate Location The number of cellular grid regions representing corrosion in the positive direction of the axis ,as well as Coordinates at the moment Location The number of cellular grid areas that exist in the negative direction of the axis, indicating corrosion It should be noted that the number of cellular grid areas that are corroded here includes the number of cellular grid areas in the corrosion end state and the number of cellular grid areas in the corrosion process.

[0067] Finally, the difference in the number of corroded cell grids in each cell grid area at different times and the previous time is analyzed, and the initial temperature data of each cell grid area at different times are combined to obtain the initial corrosion rate of each cell grid area at different times.

[0068] The change of the corrosion state of the cellular grid area at the same position at adjacent moments can reflect the corrosion preview of the cellular grid area. At the same time, factors such as ambient temperature determine the upper limit of corrosion. Therefore, according to the collected corrosion conditions of the discrete cellular grid area of ​​the steel plate, the corrosion rate is calculated. The erosion rate of the cellular grid area at time:

[0069] Since the steel plate is in the coordinate Location There are two surfaces on the shaft that may be corroded, so it is necessary to determine The cellular grid area at the position is located on which surface of the steel plate, and it is judged whether there is a cellular grid area in the positive direction of the position with the corrosion end state. If the cellular states in the positive direction of the position are all in the corrosion end state, it means The cellular grid area at the position is in the negative direction; similarly, if the cellular states in the negative direction of the position are all in the corrosion end state, it means The cell grid area at position is in the positive direction.

[0070] To facilitate subsequent analysis, we take the cell grid area in the negative direction as an example. The cellular grid area at the position and in the erosion state is calculated. The moment and the previous moment ( The difference in the number of eroded cell grids at the moment At this moment The initial temperature data of the steel plate at position is obtained by finite element analysis, and the Moment The initial corrosion rate of the cellular grid area at position and in the corroding state is:

[0071]

[0072] Where, express Moment The initial corrosion rate of the cellular grid area at position and in the corroding state; express Coordinates at the moment Location The number of cellular grid areas that exist in the negative direction of the axis, indicating corrosion (corrosion end and corrosion in progress); express Coordinates at the moment Location The number of cellular grid areas that exist in the negative direction of the axis, indicating corrosion (corrosion end and corrosion in progress); Indicates the time difference between adjacent moments ( Moment and Time difference between the moments); express At this moment Initial temperature data of the cellular mesh area at location obtained by finite element analysis.

[0073] The above initial corrosion rate is analyzed based on the corrosion change of the cellular grid area at each position at adjacent moments. When the number of cellular grid areas with corrosion change in the corrosion change direction at a certain position at adjacent moments is greater, that is, The larger the value, the faster the corrosion rate of the cellular grid area at that location during the simulation, and the greater the initial corrosion rate of the cellular grid area. At the same time, temperature also has a certain degree of influence on the corrosion rate of the cellular grid area. The increase in temperature enhances the kinetic energy of ions and molecules in the solution, increases the collision frequency and energy, and thus promotes the corrosion process, that is, The larger the value, the faster the initial erosion rate of the cellular grid area.

[0074] Based on the initial simulation data, the initial corrosion degree of each cellular grid area at different times is analyzed, including:

[0075] First, based on the initial corrosion rate and the number of corroded cell grids, the corrosion performance of each cell grid area at different times is analyzed. Then, based on the corrosion performance, combined with the initial stress data of each cell grid area at different times, the initial corrosion degree of each cell grid area at different times is obtained. Specifically, when the corrosion change rate of the cell grid area at a certain position at a certain time is faster, it means that the corrosion reaction at this position is getting more and more severe, and the corrosion at this position is deeper, so the corrosion rate is calculated. At this moment The degree of corrosion of the cellular grid area at the location for:

[0076]

[0077]

[0078] Where, express At this moment Corrosion performance of the cellular grid area at the location; express Moment The initial corrosion rate of the cellular grid area at position and in the corroding state; express Coordinates at the moment Location The number of cellular grid areas that exist in the negative direction of the axis, indicating corrosion (corrosion end and corrosion in progress); express At this moment The degree of corrosion of the cellular grid area at the location; express At this moment Initial stress data of the cell grid area at location; represents the linear normalization function.

[0079] The above formula is mainly based on Moment The corrosion performance of the cellular grid area at a certain location is used to calculate the corrosion degree of the cellular grid area at that location. The faster the corrosion rate of the steel plate at a certain location, and the larger the corrosion scale at that location, the greater the corrosion degree of the steel plate at that location. Stress data directly affects the failure behavior of the material. Larger stresses will destroy the passivation film and induce crack initiation. The larger the stress data obtained at that location, the greater the corrosion degree at that location.

[0080] The cellular mesh regions decomposed from the steel plate exhibit varying degrees of corrosion at different locations. This means the severity of the corrosion reaction at different locations varies. Under varying degrees of corrosion, the measured ambient temperature and stress data vary. Therefore, finite element analysis at different locations should be performed with varying densities to obtain temperature, stress, and other data at varying levels of detail. Therefore, the corrosion degree of the cellular mesh regions at each location, as analyzed by a three-dimensional cellular automation, helps intelligently control the finite element to divide the unit element regions into different densities.

[0081] Therefore, in an embodiment of the present invention, after obtaining the initial corrosion degree, the unit element size of the finite element mesh decomposition is set, and each cellular mesh area is divided into a number of unit element areas, including: setting the unit element size of the finite element mesh decomposition according to the initial corrosion degree; dividing each cellular mesh area into a number of unit element areas according to the unit element size, and obtaining the number of unit element areas in each cellular mesh area. The specific implementation method is as follows:

[0082] The above process obtained At this moment The degree of corrosion of the cellular grid area at the location When the degree of corrosion of the steel plate at a certain position is large, it means that the corrosion reaction at this position is more severe at the current moment, and the changes in temperature, stress and other data at this position are greater, and a more precise finite element network is needed for division.

[0083] The unit type of the finite element mesh is set to C3D8R (hexahedron), and the unit size is set to 30mm. At this moment The corrosion degree of the steel plate at the position is used to intelligently control the division size of the unit element of the cellular grid area at that position. for:

[0084]

[0085] Where, express At this moment The unit cell size within the cellular grid area at position; express At this moment The degree of corrosion of the cellular grid area at the location; Expressed as a natural constant An exponential function with base .

[0086] The above formula is mainly calculated based on the degree of corrosion of the cellular grid area. The unit element size within the cellular grid area at the current moment. When the degree of corrosion of the steel plate in a certain cellular grid area is high, it means that the corrosion reaction in this cellular grid area at the current moment is more severe, and a finer grid division is required to obtain more accurate temperature, stress, etc., so the unit element size divided at this location should be smaller.

[0087] According to the unit cell size, each cellular grid area is divided into several unit cell areas, and the number of unit cell areas in each cellular grid area is:

[0088]

[0089] Where, express At this moment The number of unit cell regions within the cellular grid region at position; express At this moment The volume of the cell grid area at the location; express At this moment The volume of the unit cell area divided within the cellular grid area at the position is calculated by the unit cell size, for example, a size value closest to the unit size is determined, and At this moment Each side length of the cellular grid area at the position can be divided by the size value closest to the unit size, and the volume of the unit cell area is calculated according to the size value closest to the unit size.

[0090] So far, through The corrosion degree of different cellular grid areas at different times was calculated, and the unit element size of the finite element mesh decomposition of each cellular grid area of ​​the steel plate was obtained.

[0091] S300: Based on the coupling mechanism of three-dimensional cellular automata and finite elements, the detailed simulation data of each unit element area at different times are obtained.

[0092] Based on the coupling mechanism of three-dimensional cellular automaton and finite element, the fine simulation data of each unit element area at different times are obtained. Among them, the fine simulation data includes fine temperature data and fine stress data. Specifically, the finite element acquisition is used to obtain At this moment The first cell grid area at the position Fine temperature data at the unit cell area , and fine stress data .

[0093] S400: According to the refined simulation data and in combination with the position distribution of the unit cell area, the adjusted simulation data of each cellular grid area at different times are obtained.

[0094] Different unit cell sizes yield different degrees of refinement in the regional temperature and stress data obtained. Smaller unit cell sizes yield more refined regional temperature and stress data, which differ from the initially measured temperature and stress data. To achieve a more refined steel plate corrosion simulation, the refined temperature and stress data analyzed must be used to provide feedback and adjust the corrosion level at different steel plate locations during the simulation.

[0095] Therefore, the temperature and stress data of the cellular grid area can be analyzed based on the temperature and stress data obtained at the unit element area in the cellular grid area of ​​the steel plate during the simulation process to obtain more accurate simulation data of the cellular grid area.

[0096] Based on the above analysis, in some embodiments of the present invention, based on the refined simulation data and the position distribution of the unit cell regions, adjusted simulation data for each cellular grid region at different times is obtained. A specific implementation method is as follows: the distance between each unit cell region and the corrosion center in each cellular grid region at different times is obtained, where the corrosion center is the centroid position of the initial corrosion core randomly selected based on environmental parameters (temperature, pH value, etc.); based on the distance, the distance weight of each unit cell region in each cellular grid region at different times is obtained as follows:

[0097]

[0098] Where, express At this moment The first cell grid area at the position The distance weight of the unit cell area; express At this moment The first cell grid area at the position The distance between the unit area and the corrosion center; express At this moment The number of unit cell regions within the cell grid region at position.

[0099] According to the detailed simulation data and combined with the distance weight, the adjusted simulation data of each cell grid area at different times are obtained, wherein the adjusted simulation data includes the adjusted temperature data and the adjusted stress data. At this moment The adjusted temperature data and adjusted stress data of the cellular grid area at position are:

[0100]

[0101]

[0102] Where, express At this moment Adjusted temperature data of the cell grid area at the location; express At this moment The first cell grid area at the position Fine temperature data at each unit cell area; express At this moment The first cell grid area at the position The distance weight of the unit cell area; express At this moment Adjusted stress data for the cellular mesh region at location; At this moment The first cell grid area at the position Detailed stress data at each unit area;

[0103] The above formula primarily utilizes the precise temperature and stress data at the unit element, combined with the location of the divided unit element from the corrosion center, for calculation. Since corrosion in steel plates spreads from a central point, locations closer to the corrosion center exhibit longer corrosion times and more stable reactions, and are therefore more representative of the relevant data for that cellular grid area. Therefore, unit element areas closer to the corrosion center are given larger weights to obtain the temperature and stress data for that cellular grid area.

[0104] At this point, the adjusted temperature data and adjusted stress data are obtained after precise analysis at each cellular grid area position of the steel plate during the simulation process.

[0105] S500: Analyze the difference between the adjusted simulation data and the initial simulation data, and combine the initial corrosion rate of the cellular grid area to obtain the adjusted corrosion rate and the adjusted corrosion degree of each cellular grid area at different times.

[0106] After obtaining the adjusted simulation data, the difference between the simulation data and the initial simulation data can be analyzed to help adjust the corrosion degree performance in different cellular grid areas during the simulation process.

[0107] Therefore, in the embodiment of the present invention, by analyzing the difference between the adjusted simulation data and the initial simulation data and combining the initial corrosion rate of the cellular grid area, the adjusted corrosion rate and adjusted corrosion degree of each cellular grid area at different times are obtained.

[0108] First, calculate the ratio of the adjusted temperature data to the initial temperature data, and combine it with the initial corrosion rate of the cellular grid area to obtain the adjusted corrosion rate of each cellular grid area at different times:

[0109]

[0110] Where, express At this moment Adjusted erosion rate of the cellular grid area at the location; express At this moment The initial corrosion rate of the cellular grid area at the location; express At this moment Adjusted temperature data of the cell grid area at the location; express At this moment Initial temperature data for the cell grid region at location.

[0111] Intelligently adjust the corrosion rate at different locations during the steel plate corrosion simulation based on temperature change data. Accurately calculate the corrosion rate during the corrosion simulation using accurate temperature data, thereby displaying a more accurate corrosion rate curve.

[0112] Then, based on the adjusted corrosion rate, the adjusted corrosion degree of each cellular grid area at different times is obtained. The specific steps are: counting the number of corroded cellular grids in the upward or downward direction of each cellular grid area at different times; analyzing the adjusted corrosion performance of each cellular grid area at different times based on the adjusted corrosion rate and the number of corroded cellular grids; based on the adjusted corrosion performance and the adjusted stress data of each cellular grid area at different times, the adjusted corrosion degree of each cellular grid area at different times is obtained as follows:

[0113]

[0114]

[0115] Where, express At this moment The cell grid area at position Adjustment of corrosion performance in direction; express At this moment Adjusted erosion rate of the cellular grid area at the location; express Coordinates at the moment Location The number of cellular grid areas that exist in the negative direction of the axis, indicating corrosion (corrosion end and corrosion in progress); express At this moment The degree of erosion of the cell grid area at the position; express At this moment Adjusted stress data for the cellular mesh region at location; represents the linear normalization function.

[0116] At this point, the feedback adjustment of corrosion performance by finite element analysis during the steel plate corrosion simulation process has been completed.

[0117] S600: Based on the adjustment of corrosion rate and corrosion degree during the steel plate corrosion simulation process, the steel plate corrosion process is dynamically displayed.

[0118] Based on the adjustment of the corrosion rate and corrosion severity during the steel plate corrosion simulation, the corrosion performance data (adjusted corrosion rate and corrosion severity) of the steel plate at different times during the simulation is accurately analyzed and displayed. At the same time, relevant data during the steel plate corrosion simulation is adjusted to make the data representation more accurate and realistic, thereby improving the accuracy of the steel plate corrosion simulation. Then, based on the relevant parameters during the steel plate corrosion simulation, the steel plate corrosion process is dynamically displayed. Figure 4-Figure 9 The figure shows the change of the numerical simulation results of steel plate corrosion over time on a vertical section (the vertical section is parallel to the plane formed by the Y-axis and the Z-axis), where Figure 4 is the corrosion result on the vertical section at 0.06s, Figure 5 is the corrosion result on the vertical section at 0.08s, Figure 6 is the corrosion result on the vertical section at 0.10s, Figure 7 is the corrosion result on the vertical section at 0.12s, Figure 8 is the corrosion result on the vertical section at 0.14s, Figure 7 This is the corrosion result on the vertical section at 0.16s.

[0119] It should be noted that the order in which the embodiments of the present invention are described above is for illustrative purposes only and does not necessarily represent the superiority or inferiority of the embodiments. The processes depicted in the accompanying drawings do not necessarily require the specific order or sequential order shown to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0120] The various embodiments in this specification are described in a progressive manner. The same or similar parts between the various embodiments can be referred to each other. Each embodiment focuses on the differences from other embodiments.

Claims

1. A numerical simulation method for steel plate corrosion by coupling three-dimensional cellular automata with finite elements, characterized in that: The method comprises: Based on the 3D cellular automation module and finite element analysis module, the initial simulation data of each cellular grid area at different times during the steel plate corrosion simulation process are collected; Based on the initial simulation data, analyzing the initial corrosion rate and initial corrosion degree of each cellular grid area at different times, setting the unit element size of the finite element mesh decomposition, and dividing each cellular grid area at different times into a number of unit element areas; Based on the three-dimensional cellular automaton module and the finite element analysis module, the detailed simulation data of each unit element area at different times are obtained; According to the refined simulation data, combined with the position distribution of the unit cell area, the adjusted simulation data of each cellular grid area at different times are obtained; Analyzing the difference between the adjusted simulation data and the initial simulation data, and combining the initial corrosion rate of the cellular grid area, obtaining the adjusted corrosion rate and the adjusted corrosion degree of each cellular grid area at different times; Based on the adjusted corrosion rate and the adjusted corrosion degree during the steel plate corrosion simulation process, the steel plate corrosion process is dynamically displayed.

2. The steel plate corrosion numerical simulation method coupled with three-dimensional cellular automation and finite element method according to claim 1 is characterized in that: The initial simulation data includes initial temperature data and initial stress data.

3. The method for numerical simulation of steel plate corrosion by coupling three-dimensional cellular automata with finite elements according to claim 2, characterized in that: Based on the initial simulation data, the initial corrosion rate of each cellular grid area at different times is analyzed, including: Count the number of eroded cellular grids in the upward or downward direction of each cellular grid area at different times; The difference in the number of corroded cellular grids in each cellular grid area at different times and the previous time is analyzed, and the initial corrosion rate of each cellular grid area at different times is obtained by combining the initial temperature data of each cellular grid area at different times.

4. The method for numerical simulation of steel plate corrosion by coupling three-dimensional cellular automata with finite elements according to claim 3, characterized in that: Based on the initial simulation data, the initial corrosion degree of each cellular grid area at different times is analyzed, including: Analyzing the corrosion performance of each cellular grid area at different times according to the initial corrosion rate and the number of corroded cellular grids; According to the corrosion performance, combined with the initial stress data of each cellular grid area at different times, the initial corrosion degree of each cellular grid area at different times is obtained.

5. The method for numerical simulation of steel plate corrosion by coupling three-dimensional cellular automata with finite elements according to claim 4, characterized in that: Set the unit element size for finite element mesh decomposition, dividing each cellular mesh region into several unit element regions, including: According to the initial corrosion degree, the unit element size of the finite element mesh decomposition is set; Each cellular grid region is divided into a number of unit cell regions according to the unit cell size, and the number of unit cell regions in each cellular grid region is obtained.

6. The method for numerical simulation of steel plate corrosion by coupling three-dimensional cellular automata with finite elements according to claim 2, characterized in that: According to the refined simulation data and in combination with the position distribution of the unit cell area, the adjusted simulation data of each cellular grid area at different times is obtained, including: Obtain the distance between each unit cell area and the corrosion center in each cellular grid area at different times; According to the distance, the distance weights of each unit cell area in each cellular grid area at different times are obtained; According to the refined simulation data and in combination with the distance weight, the adjusted simulation data of each cellular grid area at different times are obtained.

7. The method for numerical simulation of steel plate corrosion by coupling three-dimensional cellular automata with finite elements according to claim 6, characterized in that: The adjusting simulation data includes adjusting temperature data and adjusting stress data.

8. The method for numerical simulation of steel plate corrosion by coupling three-dimensional cellular automata with finite elements according to claim 7, characterized in that: Analyze the difference between the adjusted simulation data and the initial simulation data, and combine the initial corrosion rate of the cellular grid area to obtain the adjusted corrosion rate and adjusted corrosion degree of each cellular grid area at different times, including: Calculating a ratio of the adjusted temperature data to the initial temperature data, and combining the initial corrosion rate of the cellular grid area to obtain an adjusted corrosion rate of each cellular grid area at different times; According to the adjusted corrosion speed, the adjusted corrosion degree of each cellular grid area at different times is obtained.

9. The steel plate corrosion numerical simulation method coupled with three-dimensional cellular automation and finite element method according to claim 8, characterized in that: According to the adjusted corrosion rate, the adjusted corrosion degree of each cellular grid area at different times is obtained, including: Count the number of eroded cellular grids in the upward or downward direction of each cellular grid area at different times; According to the adjusted corrosion rate and the number of corroded cellular grids, analyzing the adjusted corrosion performance of each cellular grid area at different times; According to the adjusted corrosion performance, combined with the adjusted stress data of each cellular grid area at different times, the adjusted corrosion degree of each cellular grid area at different times is obtained.

10. The steel plate corrosion numerical simulation method coupled with three-dimensional cellular automation and finite element method according to claim 3, characterized in that: Count the number of eroded cellular grids in the upward or downward direction of each cellular grid area at different times. Previously, it also included: Obtaining the cellular state of the cellular grid region based on the cellular state transition rule, wherein the cellular state includes: corrosion completed, corrosion in progress, and uncorroded; Then, the number of corroded cellular grids in the upper or lower direction of each corroded cellular grid area at different times is counted.

Images