Galvanic corrosion simulation method and system for aircraft dissimilar metal lap joint structures

Abstract

The present invention relates to the technical field of corrosion simulation, and specifically provides a method and system for simulating galvanic corrosion of dissimilar metal overlap structures of aircraft, comprising: constructing a Heaviside function for characterizing boundary discontinuities based on defined domains and boundaries; converting the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and constructing an interface description model based on the basis function; constructing a level set mechanism model based on the interface description model and level set basic functions; setting the same particle concentration in the electrolyte domain, solving the level mechanism model, and obtaining dynamic changes in the corrosion interface. The present invention combines the Heaviside function and the level set mechanism, and combines the Nernst-Planck equation to construct a level set mechanism model, realizes simulation of the dynamic interface evolution of corrosion, and has a higher degree of consistency with the actual galvanic corrosion process of the dissimilar metal overlap structure of the aircraft.

Description

Technical Field

[0001] The invention belongs to the technical field of corrosion simulation, and in particular relates to a method and system for simulating galvanic corrosion of an aircraft dissimilar metal overlap structure. Background Art

[0002] Aircraft operate in complex environments, often exposed to humid, marine, and industrial atmospheres. Corrosion under these conditions is a significant form of structural damage. Because aircraft structures often contain numerous dissimilar materials, galvanic corrosion caused by these materials is one of the most common forms of corrosion in aircraft structures.

[0003] Currently, the most widely used method for galvanic corrosion structural deformation is the Arbitrary Lagrangian-Eulerian Formulation (ALE). This method not only absorbs the advantages of both the Lagrangian and Eulerian algorithms for finite element simulation, but also overcomes some of the shortcomings of these two traditional methods. Specifically, ALE discretizes the finite element model and writes the governing equations as functions of mesh coordinates. These mesh coordinates have a one-to-one mapping relationship with the spatial coordinates of the current computational domain. If the mesh displacement is too large or the mesh quality is too low during the calculation process, the computational domain will be re-meshed, and all variables to be solved will be mapped from the old mesh to the new mesh. This redefines the governing equations on the freely varying mesh. The ALE method allows the mesh coordinates to be separated from the material coordinates, so it can describe physical problems with large deformations with high accuracy. However, when the structure undergoes large topological deformation, this method will encounter serious numerical singularity problems, resulting in non-convergence of the results.

[0004] Due to the high computational cost, the grid being prone to singularities, and the computation being prone to non-convergence, existing galvanic corrosion models are unable to calculate structural deformation, and therefore cannot accurately simulate the galvanic corrosion of aircraft dissimilar metal lap joints. Summary of the Invention

[0005] In view of the above-mentioned deficiencies in the prior art, the present invention provides a method and system for simulating galvanic corrosion of aircraft dissimilar metal overlap structures to solve the above-mentioned technical problems.

[0006] In a first aspect, the present invention provides a method for simulating galvanic corrosion of an aircraft dissimilar metal overlap structure, comprising:

[0007] Based on the defined domain and boundary, a Heaviside function is constructed to characterize boundary discontinuities.

[0008] Converting the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and constructing an interface description model based on the basis function;

[0009] Based on the interface description model and the level set basic function, a level set mechanism model is constructed;

[0010] The particle concentration in the electrolyte domain is assumed to be the same, and the horizontal mechanism model is solved to obtain the dynamic changes of the corrosion interface.

[0011] In an optional embodiment, based on the defined domain and boundary, a Heaviside function for characterizing boundary discontinuities is constructed, including:

[0012] Assume that the domain Ω is divided into two subdomains Ω1 and Ω2, and Ω2=Ω\Ω1, and assume that Г is the boundary between subdomains Ω1 and Ω2, and transform the function The bound Г is implicitly defined as the regularized characteristic function, Phase transition is achieved through the interface Г;

[0013] Introduce the Heaviside function: .

[0014] In an optional embodiment, the Heaviside function is converted into a basis function that introduces a first-order fuzzy approximate interface, and an interface description model is constructed based on the basis function, including:

[0015] Introduce the Heaviside function of the first-order fuzzy approximation interface: ;

[0016] in Represents the thickness of the interface Г, which usually depends on the mesh size.

[0017] In an optional embodiment, a level set mechanism model is constructed based on the interface description model and the level set basic function, including:

[0018] Based on the non-fixedness of the interface profile and thickness, the smoothness of the discontinuous interface is maintained and a compression term is added after each time step: , represents the compression flux;

[0019] By setting , so that the compression flux acts on region and along the normal direction of the interface, where is the set time value;

[0020] To avoid breakpoints in the interface, Introducing the viscosity term into ,in stands for Laplace operator;

[0021] Will and level set basic functions Perform the combined operation to obtain the level set control function ,in , γ is the thickness change threshold, The interface thickness parameter represented by ;

[0022] Integrate over a domain: ;

[0023] Solving for the conductivity in the domain , the formula is: ,in is the electrolyte conductivity of the electrolyte, Electrolyte conductivity of the electrode;

[0024] Constructing the Heaviside function H The derivative in the normal direction, that is, the Dirac function: , where is the normal vector on the level set boundary;

[0025] Will Substituting the Dirac function, we get: ,in, ;

[0026] The surface integral on the interface Г is defined as: ;

[0027] By the Dirac function Domain Except for the boundary near Г Delete all areas outside the interface and track the interface Г;

[0028] right The function is regularized and approximated to ensure , obtained by regularization The function is: ;

[0029] After regularization The function is: ;

[0030] Introducing the electrolyte current source term Q l Convert the boundary reaction integral into a domain integral for solution, and the electrolyte current source term Q l The expression function is: , i locis the local current density on the electrode surface.

[0031] In an optional embodiment, the particle concentration in the electrolyte domain is set to be the same, and the horizontal mechanism model is solved to obtain the dynamic changes of the corrosion interface, including:

[0032] Assume that the electrolyte solution is a static dilute solution, that is, the particle concentration in the electrolyte domain is the same everywhere and there is no concentration gradient. Ignore the convection and diffusion terms of the Nernst-Planck equation. The Nernst-Planck equation includes:

[0033] ,in is the diffusion term, is the convection term, is the electromigration term, is the net flux;

[0034] Charged particles in the electrolyte solution only form current through electromigration. The current density vector in the electrolyte solution is i l The representation is:

[0035] ;

[0036] In dilute solutions, the conductivity of an electrolyte can be expressed as the sum of the contributions of each ion, i.e.:

[0037] ;

[0038] The current density vector i l The representation is converted to:

[0039] ;

[0040] Based on the preset net charge in the electrolyte domain being 0, the current conservation equation ,right The solution is:

[0041] ;

[0042] By solving the electrolyte potential φ l The Laplace equation is used to obtain the potential distribution law in the corrosion electric field, and then substitute it into the formula , and obtain the current density distribution;

[0043] The current density value is calculated according to Faraday's law and the current density value equation is:

[0044] ;

[0045] The metal corrosion rate is obtained from the current density equation:

[0046] ;

[0047] Where: It reflects the change of corrosion depth over time, that is, the corrosion rate v ; i loc is the local current density; n is the chemical equivalence coefficient; M is the relative molecular mass; is the density;

[0048] Local current density at the interface within the domain i loc for:

[0049] ;

[0050] Where: overpotential , is the electrode potential, which is taken as 0; is the electrolyte potential;

[0051] Corrosion rate v It is the driving force that controls the change of the level set function. In the two-dimensional solution domain, the level set function The velocity vector is x and y Directional component and Can be characterized as:

[0052] ;

[0053] ;

[0054] The amount and Substituting the level set governing equation for solution, the dynamic changes of the corrosion interface are obtained.

[0055] In a second aspect, the present invention provides a galvanic corrosion simulation system for aircraft dissimilar metal overlap structures, comprising:

[0056] Function building module, used to build Heaviside functions for characterizing boundary discontinuities based on defined domains and boundaries;

[0057] A function transformation module, used to convert the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and to construct an interface description model based on the basis function;

[0058] A model construction module, used for constructing a level set mechanism model based on the interface description model and the level set basic function;

[0059] The model solving module is used to set the particle concentration in the electrolyte domain to be the same, solve the horizontal mechanism model, and obtain the dynamic changes of the corrosion interface.

[0060] In an optional embodiment, the function building module includes:

[0061] Assume that the domain Ω is divided into two subdomains Ω1 and Ω2, and Ω2=Ω\Ω1, and assume that Г is the boundary between subdomains Ω1 and Ω2, and transform the function The bound Г is implicitly defined as the regularized characteristic function, Phase transition is achieved through the interface Г;

[0062] Introduce the Heaviside function: .

[0063] In an optional embodiment, the function transformation module includes:

[0064] Introduce the Heaviside function of the first-order fuzzy approximation interface: ;

[0065] in Represents the thickness of the interface Г, which usually depends on the mesh size.

[0066] In an optional embodiment, the model building module includes:

[0067] Based on the non-fixedness of the interface profile and thickness, the smoothness of the discontinuous interface is maintained and a compression term is added after each time step: , represents the compression flux;

[0068] By setting , so that the compression flux acts on region and along the normal direction of the interface, where is the set time value;

[0069] To avoid breakpoints in the interface, Introducing the viscosity term into ,in stands for Laplace operator;

[0070] Will and level set basic functions Perform the combined operation to obtain the level set control function ,in , γ is the thickness change threshold, The interface thickness parameter represented by ;

[0071] Integrate over a domain: ;

[0072] Solving for the conductivity in the domain , the formula is: ,in is the electrolyte conductivity of the electrolyte, Electrolyte conductivity of the electrode;

[0073] Constructing the Heaviside function H The derivative in the normal direction, that is, the Dirac function: , where is the normal vector on the level set boundary;

[0074] Will Substituting the Dirac function, we get: ,in, ;

[0075] The surface integral on the interface Г is defined as: ;

[0076] By the Dirac function Domain Except for the boundary near Г Delete all areas outside the interface and track the interface Г;

[0077] right The function is regularized and approximated to ensure , obtained by regularization The function is: ;

[0078] After regularization The function is: ;

[0079] Introducing the electrolyte current source term Q l Convert the boundary reaction integral into a domain integral for solution, and the electrolyte current source term Q l The expression function is: , i loc is the local current density on the electrode surface.

[0080] In an optional embodiment, the model solving module includes:

[0081] Assume that the electrolyte solution is a static dilute solution, that is, the particle concentration in the electrolyte domain is the same everywhere and there is no concentration gradient. Ignore the convection and diffusion terms of the Nernst-Planck equation. The Nernst-Planck equation includes:

[0082] ,in is the diffusion term, is the convection term, is the electromigration term, is the net flux;

[0083] Charged particles in the electrolyte solution only form current through electromigration. The current density vector in the electrolyte solution is i l The representation is:

[0084] ;

[0085] In dilute solutions, the conductivity of an electrolyte can be expressed as the sum of the contributions of each ion, i.e.:

[0086] ;

[0087] The current density vector i l The representation is converted to:

[0088] ;

[0089] Based on the preset net charge in the electrolyte domain being 0, the current conservation equation ,right The solution is:

[0090] ;

[0091] By solving the electrolyte potential φ l The Laplace equation is used to obtain the potential distribution law in the corrosion electric field, and then substitute it into the formula , and obtain the current density distribution;

[0092] The current density value is calculated according to Faraday's law and the current density value equation is:

[0093] ;

[0094] The metal corrosion rate is obtained from the current density equation:

[0095] ;

[0096] Where: It reflects the change of corrosion depth over time, that is, the corrosion rate v ; i loc is the local current density; n is the chemical equivalence coefficient; M is the relative molecular mass; is the density;

[0097] Local current density at the interface within the domain i loc for:

[0098] ;

[0099] Where: overpotential , is the electrode potential, which is taken as 0; is the electrolyte potential;

[0100] Corrosion rate v It is the driving force that controls the change of the level set function. In the two-dimensional solution domain, the level set function The velocity vector is x and y Directional component and Can be characterized as:

[0101] ;

[0102] ;

[0103] The amount and Substituting the level set governing equation for solution, the dynamic changes of the corrosion interface are obtained.

[0104] The beneficial effect of the present invention is that the method and system for simulating galvanic corrosion of aircraft dissimilar metal lap structures provided by the present invention combine the Heaviside function and the level set mechanism, and combine with the Nernst-Planck equation to construct a level set mechanism model to realize the simulation of the dynamic interface evolution of corrosion, which is more consistent with the actual galvanic corrosion process of aircraft dissimilar metal lap structures.

[0105] In addition, the present invention has a reliable design principle, a simple structure and a very broad application prospect. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0107] Figure 1 is a schematic flow chart of a method according to an embodiment of the present invention.

[0108] Figure 2 It is a schematic diagram of a model of a method according to an embodiment of the present invention.

[0109] Figure 3 FIG. 4 is a schematic block diagram of a system according to an embodiment of the present invention. DETAILED DESCRIPTION

[0110] In order to enable those skilled in the art to better understand the technical solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of the present invention.

[0111] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art of the present invention. The terms used in this specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention.

[0112] The method for simulating galvanic corrosion of an aircraft dissimilar metal lap joint structure provided by an embodiment of the present invention is executed by a computer device. Accordingly, the system for simulating galvanic corrosion of an aircraft dissimilar metal lap joint structure runs in the computer device.

[0113] Figure 1 is a schematic flow chart of a method according to an embodiment of the present invention. Figure 1 The execution subject can be a galvanic corrosion simulation system for aircraft dissimilar metal overlap structures. According to different requirements, the order of the steps in the flow chart can be changed, and some steps can be omitted.

[0114] like Figure 1 As shown, the method includes:

[0115] Step 110 , constructing a Heaviside function for characterizing boundary discontinuities based on the defined domain and boundary;

[0116] Step 120, converting the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and constructing an interface description model based on the basis function;

[0117] Step 130: constructing a level set mechanism model based on the interface description model and the level set basic function;

[0118] Step 140 , setting the particle concentration in the electrolyte domain to be the same, solving the horizontal mechanism model, and obtaining the dynamic changes of the corrosion interface.

[0119] To facilitate understanding of the present invention, the following further describes the method for simulating galvanic corrosion of dissimilar metal lap joints of aircraft provided by the present invention based on the principle of the method for simulating galvanic corrosion of dissimilar metal lap joints of aircraft provided by the present invention and the process of simulating galvanic corrosion of dissimilar metal lap joints of aircraft provided by the present invention in combination with the embodiment.

[0120] For details, please refer to Figure 2 The simulation method of galvanic corrosion of aircraft dissimilar metal joint structures includes:

[0121] S1. Based on the defined domain and boundary, a Heaviside function is constructed to characterize boundary discontinuities.

[0122] The level set method is a conservative level set tracking dynamic interface model under the premise of a given velocity vector. It has been widely used in the tracking of phase interfaces in two-phase flow. Assume that the domain Ω is divided into two subdomains Ω1 and Ω2, and Ω2=Ω\Ω1. At the same time, assume that Г is the boundary between subdomains Ω1 and Ω2. The bound Г is implicitly defined as the regularized characteristic function, Phase transition is achieved through the interface Г;

[0123] Introduce the Heaviside function: .

[0124] S2. Convert the Heaviside function into a basis function for introducing a first-order fuzzy approximation interface, and construct an interface description model based on the basis function.

[0125] Usually choose The contour line of is taken as the interface Г. Since this method is conservative, is completely conserved. If the Heaviside function with a clear interface is used, the evolution of the closed curve (surface) can be accurately tracked. However, on a discrete grid, it is difficult to achieve a sudden change from 0 to 1 in the position of the interface. To enhance numerical robustness, a first-order fuzzy approximation of the interface Heaviside function is introduced: ;

[0126] in Represents the thickness of the interface Г, which usually depends on the mesh size.

[0127] S3. Constructing a level set mechanism model based on the interface description model and the level set basic function.

[0128] Since the contour and thickness of the interface are difficult to remain constant during the calculation process, based on the instability of the contour and thickness of the interface and maintaining the smoothness of the discontinuous interface, a compression term is added after each time step. This can be regarded as an intermediate step in solving the conservation law: , represents the compression flux;

[0129] By setting , so that the compression flux acts on region and along the normal direction of the interface, where It is the set time value, not the actual time t .

[0130] It can be seen that the formula is a hyperbolic differential equation. Assuming time As the interface increases, the stability of the interface is difficult to ensure. Introducing the viscosity term into ,in stands for Laplace operator.

[0131] Will and level set basic functions Perform the combined operation to obtain the level set control function ;

[0132] in , γ is the thickness change threshold, reinitialize the parameters γ It ensures that the gradient in the level set function gradually concentrates on the free surface, that is, the phase interface, over time. It does not set the surface thickness, but ensures that the changes in the level set function are within the thickness range. If the parameter value is set too low, the changes in the level set function may be trapped in the same domain, and then a two-phase interface is generated out of thin air. Too high a value often results in too small a time step and too long a calculation time. The interface thickness parameter is more intuitive to understand, and its only function is to control the thickness of the region where the level set function changes. Too thick will blur the shape of the free surface. Although this makes the problem easier to solve, it loses the shape accuracy of the free surface. A value smaller than the mesh cell size is good for generating clear surfaces, but this also requires the mesh to have a corresponding resolution. Setting a value smaller than the mesh cell size is meaningless, as this will make the interface irregular and the surface tension that cannot be resolved may cause velocity distortion. The setting is generally 1 / 4 of the largest unit near the surface.

[0133] The integral over the domain is defined using the Heaviside function: ;

[0134] To solve the conductivity in the domain For example, the formula is: ,in is the electrolyte conductivity of the electrolyte, Electrolyte conductivity of the electrode; electrolyte conductivity in the electrolyte domain describes actual chemical problems. The definition of electrolyte conductivity in the electrode domain is only to help numerical convergence and has no practical significance.

[0135] Since the electrode reactions all occur on the electrode surface, that is, the interface between the electrolyte and the electrode, according to the definition, the Heaviside function H The derivative in the normal direction is the Dirac function: , where is the normal vector on the level set boundary, Substituting the Dirac function, we get:

[0136] ,

[0137] in, .

[0138] The surface integral on the interface Г is defined as: ;

[0139] Dirac function Deletable domain Except for the boundary near Г All areas outside the interface are tracked.

[0140] right The function is regularized and approximated to ensure , obtained by regularization The function is: ;

[0141] After regularization The function is: ;

[0142] The local current density on the electrode surface i loc For example, when defining the electrode surface reaction, unlike the ALE method, simply defining the cathode and anode reactions on the boundary cannot meet the requirements. The electrolyte current source term is introduced. Q l Convert the boundary reaction integral into a domain integral for solution, and the electrolyte current source term Q l The expression function is: , i locis the local current density on the electrode surface.

[0143] S4. Assuming the same particle concentration in the electrolyte domain, solving the horizontal mechanism model to obtain the dynamic changes of the corrosion interface.

[0144] Establish assumptions upfront, including:

[0145] (1) Assume that the oxygen reduction reaction only occurs on the surface of TA15 titanium alloy, and the active dissolution of metal only occurs on the surface of 2A12 aluminum alloy;

[0146] (2) Assume that the electrolyte is an infinitely dilute solution, that is, all physical parameters are independent of the substance concentration;

[0147] (3) An acidic electrolyte solution with a pH of 4 is used. It is assumed that the corrosion products do not undergo dense deposition and have no effect on the corrosion reaction.

[0148] Assume that the electrolyte solution is a static dilute solution, that is, the particle concentration in the electrolyte domain is the same everywhere and there is no concentration gradient. Ignore the convection and diffusion terms of the Nernst-Planck equation. The Nernst-Planck equation includes:

[0149] ,in is the diffusion term, is the convection term, is the electromigration term, is the net flux;

[0150] Charged particles in the electrolyte solution only form current through electromigration. The current density vector in the electrolyte solution is i l The representation is:

[0151] ;

[0152] In dilute solutions, the conductivity of an electrolyte can be expressed as the sum of the contributions of each ion, i.e.:

[0153] ;

[0154] The current density vector i l The representation is converted to:

[0155] ;

[0156] This equation has the same form as Ohm's law, based on the assumption that the net charge in the electrolyte domain is zero, and the current conservation equation ,right The solution is:

[0157] ;

[0158] By solving the electrolyte potential φ l The Laplace equation is used to obtain the potential distribution law in the corrosion electric field, and then substitute it into the formula , and obtain the current density distribution.

[0159] The current density value is calculated according to Faraday's law and the current density value equation is:

[0160] ;

[0161] The metal corrosion rate is obtained from the current density equation:

[0162] ;

[0163] Where: It reflects the change of corrosion depth over time, that is, the corrosion rate v ; i loc is the local current density; n is the chemical equivalence coefficient; M is the relative molecular mass; is the density.

[0164] Local current density at the interface within the domain i loc for:

[0165] ;

[0166] Where: overpotential , is the electrode potential, which is taken as 0; is the electrolyte potential.

[0167] Corrosion rate v It is the driving force that controls the change of the level set function. In the two-dimensional solution domain, the level set function The velocity vector is x and y Directional component and Can be characterized as:

[0168] ;

[0169] ;

[0170] The amount and Substitute the level set governing equation in step S3 and solve it to obtain the dynamic changes of the corrosion interface.

[0171] In some embodiments, the aircraft dissimilar metal lap joint structure galvanic corrosion simulation system may include multiple functional modules composed of computer program segments. The computer program of each program segment in the aircraft dissimilar metal lap joint structure galvanic corrosion simulation system may be stored in a memory of a computer device and executed by at least one processor to perform (see Figure 1 Description) Functionality for galvanic corrosion simulation of aircraft dissimilar metal joints.

[0172] In this embodiment, the aircraft dissimilar metal overlap structure galvanic corrosion simulation system can be divided into multiple functional modules according to the functions it performs, such as Figure 3 As shown. The functional modules of system 300 may include: a function construction module 310, a function transformation module 320, a model construction module 330, and a model solution module 340. A module as referred to in the present invention refers to a series of computer program segments that can be executed by at least one processor and can perform fixed functions, and is stored in a memory. In this embodiment, the functions of each module will be described in detail in subsequent embodiments.

[0173] Function building module, used to build Heaviside functions for characterizing boundary discontinuities based on defined domains and boundaries;

[0174] A function transformation module, used to convert the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and to construct an interface description model based on the basis function;

[0175] A model construction module, used for constructing a level set mechanism model based on the interface description model and the level set basic function;

[0176] The model solving module is used to set the particle concentration in the electrolyte domain to be the same, solve the horizontal mechanism model, and obtain the dynamic changes of the corrosion interface.

[0177] Optionally, as an embodiment of the present invention, the function building module includes:

[0178] Assume that the domain Ω is divided into two subdomains Ω1 and Ω2, and Ω2=Ω\Ω1, and assume that Г is the boundary between subdomains Ω1 and Ω2, and transform the function The bound Г is implicitly defined as the regularized characteristic function, Phase transition is achieved through the interface Г;

[0179] Introduce the Heaviside function: .

[0180] Optionally, as an embodiment of the present invention, the function transformation module includes:

[0181] Introduce the Heaviside function of the first-order fuzzy approximation interface: ;

[0182] in Represents the thickness of the interface Г, which usually depends on the mesh size.

[0183] Optionally, as an embodiment of the present invention, the model building module includes:

[0184] Based on the non-fixedness of the interface profile and thickness, the smoothness of the discontinuous interface is maintained and a compression term is added after each time step: , represents the compression flux;

[0185] By setting , so that the compression flux acts on region and along the normal direction of the interface, where is the set time value;

[0186] To avoid breakpoints in the interface, Introducing the viscosity term into ,in stands for Laplace operator;

[0187] Will and level set basic functions Perform the combined operation to obtain the level set control function ,in , γ is the thickness change threshold, The interface thickness parameter represented by ;

[0188] Integrate over a domain: ;

[0189] Solving for the conductivity in the domain , the formula is: ,in is the electrolyte conductivity of the electrolyte, Electrolyte conductivity of the electrode;

[0190] Constructing the Heaviside function H The derivative in the normal direction, that is, the Dirac function: , where is the normal vector on the level set boundary;

[0191] Will Substituting the Dirac function, we get: ,in, ;

[0192] The surface integral on the interface Г is defined as: ;

[0193] By the Dirac function Domain Except for the boundary near Г Delete all areas outside the interface and track the interface Г;

[0194] right The function is regularized and approximated to ensure , obtained by regularization The function is: ;

[0195] After regularization The function is: ;

[0196] Introducing the electrolyte current source term Q l Convert the boundary reaction integral into a domain integral for solution, and the electrolyte current source term Q l The expression function is: , i loc is the local current density on the electrode surface.

[0197] Optionally, as an embodiment of the present invention, the model solving module includes:

[0198] Assume that the electrolyte solution is a static dilute solution, that is, the particle concentration in the electrolyte domain is the same everywhere and there is no concentration gradient. Ignore the convection and diffusion terms of the Nernst-Planck equation. The Nernst-Planck equation includes:

[0199] ,in is the diffusion term, is the convection term, is the electromigration term, is the net flux;

[0200] Charged particles in the electrolyte solution only form current through electromigration. The current density vector in the electrolyte solution is i l The representation is:

[0201] ;

[0202] In dilute solutions, the conductivity of an electrolyte can be expressed as the sum of the contributions of each ion, i.e.:

[0203] ;

[0204] The current density vector il The representation is converted to:

[0205] ;

[0206] Based on the preset net charge in the electrolyte domain being 0, the current conservation equation ,right The solution is:

[0207] ;

[0208] By solving the electrolyte potential φ l The Laplace equation is used to obtain the potential distribution law in the corrosion electric field, and then substitute it into the formula , and obtain the current density distribution;

[0209] The current density value is calculated according to Faraday's law and the current density value equation is:

[0210] ;

[0211] The metal corrosion rate is obtained from the current density equation:

[0212] ;

[0213] Where: It reflects the change of corrosion depth over time, that is, the corrosion rate v ; i loc is the local current density; n is the chemical equivalence coefficient; M is the relative molecular mass; is the density;

[0214] Local current density at the interface within the domain i loc for:

[0215] ;

[0216] Where: overpotential , is the electrode potential, which is taken as 0; is the electrolyte potential;

[0217] Corrosion rate v It is the driving force that controls the change of the level set function. In the two-dimensional solution domain, the level set function The velocity vector is x and y Directional component and Can be characterized as:

[0218] ;

[0219] ;

[0220] The amount and Substituting the level set governing equation for solution, the dynamic changes of the corrosion interface are obtained.

[0221] Although the present invention has been described in detail with reference to the accompanying drawings and in conjunction with preferred embodiments, the present invention is not limited thereto. Without departing from the spirit and essence of the present invention, persons of ordinary skill in the art may make various equivalent modifications or substitutions to the embodiments of the present invention, and such modifications or substitutions shall be within the scope of the present invention. Any changes or substitutions that can be easily conceived by persons skilled in the art within the technical scope disclosed in the present invention shall be within the scope of protection of the present invention.

Claims

1. A method for simulating galvanic corrosion of aircraft dissimilar metal lap joint structures, characterized in that: include: Based on the defined domain and boundary, a Heaviside function is constructed to characterize boundary discontinuities. Converting the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and constructing an interface description model based on the basis function; Based on the interface description model and the level set basic function, a level set mechanism model is constructed; Assuming the same particle concentration in the electrolyte domain, the horizontal mechanism model is solved to obtain the dynamic changes of the corrosion interface; Based on the defined domain and boundary, the Heaviside function is constructed to characterize the boundary discontinuity, including: Assume that the domain Ω is divided into two subdomains Ω1 and Ω2, and Ω2=Ω\Ω1, and assume that Г is the boundary between subdomains Ω1 and Ω2, and transform the function The bound Г is implicitly defined as the regularized characteristic function, Phase transition is achieved through the interface Г; Introduce the Heaviside function: ; The Heaviside function is converted into a basis function for introducing a first-order fuzzy approximate interface, and an interface description model is constructed based on the basis function, including: Introduce the Heaviside function of the first-order fuzzy approximation interface: ; in represents the thickness of the interface Г, which depends on the mesh size; Based on the interface description model and the level set basic function, a level set mechanism model is constructed, including: Based on the non-fixedness of the interface profile and thickness, the smoothness of the discontinuous interface is maintained and a compression term is added after each time step: , represents the compression flux; By setting , so that the compression flux acts on region and along the interface normal direction, where is the set time value; To avoid breakpoints in the interface, Introducing the viscosity term into ,in stands for Laplace operator; Will and level set basic functions Perform the combined operation to obtain the level set control function ,in , γ is the thickness change threshold, The interface thickness parameter represented by ; Integrate over a domain: ; Solving for the conductivity in the domain , the formula is: ,in is the electrolyte conductivity of the electrolyte, Electrolyte conductivity of the electrode; Constructing the Heaviside function H The derivative in the normal direction, that is, the Dirac function: , where is the normal vector on the level set boundary; Will Substituting the Dirac function, we get: ,in, ; The surface integral on the interface Г is defined as: ; By the Dirac function Domain Except for the boundary near Г Delete all areas outside the interface and track the interface Г; right The function is regularized and approximated to ensure , obtained by regularization The function is: ; After regularization The function is: ; Introducing the electrolyte current source term Q l Convert the boundary reaction integral into a domain integral for solution, and the electrolyte current source term Q l The expression function is: , i loc is the local current density on the electrode surface.

2. The method according to claim 1, characterized in that Assuming the same particle concentration in the electrolyte domain, the horizontal mechanism model is solved to obtain the dynamic changes of the corrosion interface, including: Assume that the electrolyte solution is a static dilute solution, that is, the particle concentration in the electrolyte domain is the same everywhere and there is no concentration gradient. Ignore the convection and diffusion terms of the Nernst-Planck equation. The Nernst-Planck equation includes: ,in is the diffusion term, is the convection term, is the electromigration term, is the net flux; Charged particles in the electrolyte solution only form current through electromigration. The current density vector in the electrolyte solution is i l The representation is: ； In dilute solutions, the conductivity of an electrolyte can be expressed as the sum of the contributions of each ion, i.e.: ； The current density vector i l The representation is converted to: ； Based on the preset net charge in the electrolyte domain being 0, the current conservation equation ,right The solution is: ； By solving the electrolyte potential φ l The Laplace equation is used to obtain the potential distribution law in the corrosion electric field, and then substitute it into the formula , and obtain the current density distribution; The current density value is calculated according to Faraday's law and the current density value equation is: ； The metal corrosion rate is obtained from the current density equation: ； Where: It reflects the change of corrosion depth over time, that is, the corrosion rate v ; i loc is the local current density; n is the chemical equivalence coefficient; M is the relative molecular mass; is the density; Local current density at the interface within the domain i loc for: ； Where: overpotential , is the electrode potential, which is taken as 0; is the electrolyte potential; Corrosion rate v It is the driving force that controls the change of the level set function. In the two-dimensional solution domain, the level set function The velocity vector is x and y Directional component and Can be characterized as: ； ； The amount and Substituting the level set governing equation for solution, the dynamic changes of the corrosion interface are obtained.

3. A galvanic corrosion simulation system for aircraft dissimilar metal lap joint structures, characterized in that: include: Function building module, used to build Heaviside functions for characterizing boundary discontinuities based on defined domains and boundaries; A function transformation module, used to convert the Heaviside function into a basis function for introducing a first-order fuzzy approximate interface, and to construct an interface description model based on the basis function; A model construction module, used for constructing a level set mechanism model based on the interface description model and the level set basic function; A model solving module is used to set the particle concentration in the electrolyte domain to be the same, solve the horizontal mechanism model, and obtain the dynamic changes of the corrosion interface; The function building block includes: Assume that the domain Ω is divided into two subdomains Ω1 and Ω2, and Ω2=Ω\Ω1, and assume that Г is the boundary between subdomains Ω1 and Ω2, and transform the function The bound Г is implicitly defined as the regularized characteristic function, Phase transition is achieved through the interface Г; Introduce the Heaviside function: ; The function transformation module includes: Introduce the Heaviside function of the first-order fuzzy approximation interface: ; in represents the thickness of the interface Г, which depends on the mesh size; The model building module includes: Based on the non-fixedness of the interface profile and thickness, the smoothness of the discontinuous interface is maintained and a compression term is added after each time step: , represents the compression flux; By setting , so that the compression flux acts on region and along the interface normal direction, where is the set time value; To avoid breakpoints in the interface, Introducing the viscosity term into ,in stands for Laplace operator; Will and level set basic functions Perform the combined operation to obtain the level set control function ,in , γ is the thickness change threshold, The interface thickness parameter represented by ; Integrate over a domain: ; Solving for the conductivity in the domain , the formula is: ,in is the electrolyte conductivity of the electrolyte, Electrolyte conductivity of the electrode; Constructing the Heaviside function H The derivative in the normal direction, that is, the Dirac function: , where is the normal vector on the level set boundary; Will Substituting the Dirac function, we get: ,in, ; The surface integral on the interface Г is defined as: ; By the Dirac function Domain Except for the boundary near Г Delete all areas outside the interface and track the interface Г; right The function is regularized and approximated to ensure , obtained by regularization The function is: ; After regularization The function is: ; Introducing the electrolyte current source term Q l Convert the boundary reaction integral into a domain integral for solution, and the electrolyte current source term Q l The expression function is: , i loc is the local current density on the electrode surface.

4. The system according to claim 3, characterized in that The model solving module includes: Assume that the electrolyte solution is a static dilute solution, that is, the particle concentration in the electrolyte domain is the same everywhere and there is no concentration gradient. Ignore the convection and diffusion terms of the Nernst-Planck equation. The Nernst-Planck equation includes: ,in is the diffusion term, is the convection term, is the electromigration term, is the net flux; Charged particles in the electrolyte solution only form current through electromigration. The current density vector in the electrolyte solution is i l The representation is: ； In dilute solutions, the conductivity of an electrolyte can be expressed as the sum of the contributions of each ion, i.e.: ； The current density vector i l The representation is converted to: ； Based on the preset net charge in the electrolyte domain being 0, the current conservation equation ,right The solution is: ； By solving the electrolyte potential φ l The Laplace equation is used to obtain the potential distribution law in the corrosion electric field, and then substitute it into the formula , and obtain the current density distribution; The current density value is calculated according to Faraday's law and the current density value equation is: ； The metal corrosion rate is obtained from the current density equation: ； Where: It reflects the change of corrosion depth over time, that is, the corrosion rate v ; i loc is the local current density; n is the chemical equivalence coefficient; M is the relative molecular mass; is the density; Local current density at the interface within the domain i loc for: ； Where: overpotential , is the electrode potential, which is taken as 0; is the electrolyte potential; Corrosion rate v It is the driving force that controls the change of the level set function. In the two-dimensional solution domain, the level set function The velocity vector is x and y Directional component and Can be characterized as: ； ； The amount and Substituting the level set governing equation for solution, the dynamic changes of the corrosion interface are obtained.

Images