CN117150867B - Simulation method, device and equipment for chloride ion erosion process in marine building - Google Patents
Simulation method, device and equipment for chloride ion erosion process in marine building Download PDFInfo
- Publication number
- CN117150867B CN117150867B CN202311415094.4A CN202311415094A CN117150867B CN 117150867 B CN117150867 B CN 117150867B CN 202311415094 A CN202311415094 A CN 202311415094A CN 117150867 B CN117150867 B CN 117150867B
- Authority
- CN
- China
- Prior art keywords
- order
- chloride ion
- time
- derivative
- preset
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 title claims abstract description 125
- 238000000034 method Methods 0.000 title claims abstract description 101
- 230000003628 erosive effect Effects 0.000 title claims abstract description 70
- 230000008569 process Effects 0.000 title claims abstract description 66
- 238000004088 simulation Methods 0.000 title claims abstract description 65
- 238000009792 diffusion process Methods 0.000 claims abstract description 75
- 238000009826 distribution Methods 0.000 claims abstract description 39
- 230000015654 memory Effects 0.000 claims description 24
- 230000006870 function Effects 0.000 claims description 22
- 230000002123 temporal effect Effects 0.000 claims description 18
- 238000000992 sputter etching Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 abstract description 6
- 238000005536 corrosion prevention Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 description 14
- 238000005260 corrosion Methods 0.000 description 5
- 230000007797 corrosion Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000013535 sea water Substances 0.000 description 3
- 239000004568 cement Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002161 passivation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of engineering simulation and numerical simulation, and discloses a simulation method, a device and equipment for a chloride ion erosion process in a marine structure, wherein the method utilizes historical erosion data of a second sea area to invert model parameters of a one-dimensional variable order time fractal derivative flow diffusion model under a finite difference discrete format to obtain an inversion parameter set; then, inputting the inversion parameter set and the chloride ion concentration of the first sea area into a model, and solving the model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area in a preset simulation time length; the one-dimensional variable-order time fractal derivative flow diffusion model considers the change of the erosion rate of chloride ions to concrete in a long-time process, so that the erosion process is accurately simulated, the solving efficiency and the accuracy of the model are effectively improved on the premise of not increasing space distribution points, and powerful technical support is provided for the corrosion prevention work of marine structures.
Description
Technical Field
The invention relates to the technical field of engineering simulation and numerical simulation, in particular to a simulation method, a simulation device and simulation equipment for a chloride ion erosion process in a marine building.
Background
With the gradual development of ocean resources, marine structures are gradually increased, and corrosion problems become important influencing factors of the service lives of the marine structures. The chloride ion corrosion is the most main seawater corrosion effect, simulates the diffusion of chloride ions in the marine structure for a long time, and has the advantages of low cost, intuitionism and convenience in evaluating the seawater corrosion effect and analyzing the service life of the marine structure.
Because the chloride ion erosion process in concrete is an abnormal diffusion phenomenon, namely the square of the mean square displacement of chloride ion erosion is not in direct proportion to time, and the fractal derivative model has the characteristics of simple structure and capability of describing complex diffusion behavior of solute, and is often used for describing the solute migration process in complex media. However, the erosion rate of chloride ions on concrete is not always consistent because the erosion of chloride ions over a long period of time can reduce the alkaline environment inside the concrete, destroy the passivation layer therein and reduce the barrier of the concrete to the migration of chloride ions.
Therefore, if the erosion process of chloride ions to concrete is simulated by using the fractal derivative model, the change of the erosion rate is not considered, so that the simulation result is inaccurate.
Disclosure of Invention
In view of the above, the invention provides a simulation method, a device and equipment for a chloride ion erosion process in a marine structure, so as to solve the problem of inaccurate simulation of the chloride ion erosion process on concrete.
In a first aspect, the invention provides a method for simulating a chloride ion erosion process in a marine structure, the method comprising:
acquiring sea state data of a first sea area and historical erosion data of a second sea area, wherein the sea state data comprises chloride ion concentration, and the sea state data of the second sea area and the sea state data of the first sea area belong to the same interval; inverting model parameters of a flow diffusion model by utilizing historical erosion data and carrying out one-dimensional variable order time fractal derivative under a finite difference discrete format to obtain an inversion parameter set; inputting the chloride ion concentration and inversion parameter set into a one-dimensional variable-order time fractal derivative convective diffusion model, and solving the one-dimensional variable-order time fractal derivative convective diffusion model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area within a preset simulation duration.
According to the simulation method of the chloride ion erosion process in the marine structure, firstly, the model parameters of a one-dimensional variable order time fractal derivative flow diffusion model under a finite difference discrete format are inverted by using the historical erosion data of a second sea area, so that an inversion parameter set is obtained; then, inputting the inversion parameter set and the chloride ion concentration of the first sea area into a one-dimensional variable-order time fractal derivative convection diffusion model, and solving the model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area in a preset simulation duration; the one-dimensional variable-order time fractal derivative flow diffusion model considers the change of the erosion rate of chloride ions to concrete in a long-time process, so that the erosion process is accurately simulated, the solving efficiency and the accuracy of the model are effectively improved on the premise of not increasing space distribution points, and powerful technical support is provided for the corrosion prevention work of marine structures.
In an alternative embodiment, the inversion parameter set includes a variable order fractal derivative order function, a diffusion coefficient, and a convection coefficient.
In an alternative embodiment, inverting model parameters of a flow diffusion model by one-dimensional variable order time fractal derivative under a finite difference discrete format to obtain an inversion parameter set includes:
and inverting model parameters of the flow diffusion model by adopting a least square fitting method to obtain an inversion parameter set.
In an alternative embodiment, the one-dimensional variable order time fractal derivative model is constructed using the following formula:
wherein,fractal derivative for time variant order, +.>For the variable order fractal derivative order function,，/>for convection parameters->For diffusion parameter->Is->Temporal spatial position->Chloride ion concentration at the site.
In an alternative embodiment, inputting the chloride ion concentration and inversion parameter set into a one-dimensional variable-order time fractal derivative convective diffusion model, and solving the one-dimensional variable-order time fractal derivative convective diffusion model by using a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of a first sea area within a preset simulation duration, wherein the method comprises the following steps:
determining a time step based on preset time distribution information and a variable order fractal derivative order function; solving a one-dimensional variable-order time fractal derivative flow diffusion model by using a finite difference discrete format of a time step to obtain chloride ion concentrations respectively corresponding to each moment at different positions; and obtaining the chloride ion erosion process of the first sea area in the preset simulation time based on the chloride ion concentrations respectively corresponding to the moments at different positions.
Because the numerical accuracy of the uniform distribution mode is low, the time step determined by the embodiment is continuously changed along with the fractal derivative order of the variable order. Because the lower the order is, the lower the numerical accuracy is, the embodiment improves the numerical accuracy by arranging more monitoring moments at the lower order, thereby improving the accuracy of the simulation process.
In an alternative embodiment, the finite difference discrete format is:
wherein,is->Fractional derivative order of corresponding order at each moment,/->Is the firstIndividual temporal spatial position->Concentration of chloride ion at->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->Time steps corresponding to the respective moments.
In an alternative embodiment, the preset time distribution information includes a preset analog duration and a preset total number of time distribution, and the time step is determined by the following formula:
wherein,is->Time step corresponding to each moment +.>Is->At all times->Is->At all times->For a preset simulation time length, < >>Is the total number of preset time points->Is->The fractional derivative order of the corresponding order at each moment.
In a second aspect, the present invention provides an apparatus for simulating a chloride ion etching process in a marine structure, the apparatus comprising:
the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring sea state data of a first sea area and historical erosion data of a second sea area, the sea state data comprise chloride ion concentration, and the sea state data of the second sea area and the sea state data of the first sea area belong to the same interval; the inversion module is used for inverting model parameters of the flow diffusion model by utilizing the historical erosion data and carrying out inversion on one-dimensional variable order time fractal derivative under the finite difference discrete format to obtain an inversion parameter set; the generation module is used for inputting the chloride ion concentration and the inversion parameter set into the one-dimensional variable-order time fractal derivative convection diffusion model, and solving the one-dimensional variable-order time fractal derivative convection diffusion model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area within a preset simulation duration.
In an alternative embodiment, the inversion parameter set generated by the inversion module includes a variable order fractal derivative order function, a diffusion coefficient, and a convection coefficient.
In an alternative embodiment, an inversion module includes:
and the inversion sub-module is used for inverting the model parameters of the flow diffusion model by adopting a least square fitting method to obtain an inversion parameter set.
In an alternative embodiment, the one-dimensional variable order time fractal derivative in the inversion module builds the convective diffusion model using the following formula:
wherein,fractal derivative for time variant order, +.>For the variable order fractal derivative order function,，/>for convection parameters->For diffusion parameter->Is->Temporal spatial position->Chloride ion concentration at the site.
In an alternative embodiment, the generating module includes:
the determining submodule is used for determining a time step based on preset time distribution information and a variable order fractal derivative order function; the solving submodule is used for solving the one-dimensional variable order time fractal derivative flow diffusion model by utilizing a finite difference discrete format of taking time step length to obtain chloride ion concentrations respectively corresponding to each moment at different positions; and the determination submodule is used for obtaining the chloride ion erosion process of the first sea area in the preset simulation time based on the chloride ion concentrations respectively corresponding to the moments at different positions.
In an alternative embodiment, the finite difference discrete format in the solution sub-module is:
wherein,is->Fractional derivative order of corresponding order at each moment,/->Is the firstIndividual temporal spatial position->Concentration of chloride ion at->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->Time steps corresponding to the respective moments.
In a third aspect, the present invention provides a computer device comprising: the processor executes the computer instructions, thereby executing the simulation method of the chloride ion erosion process in the marine engineering building according to the first aspect or any corresponding implementation mode thereof.
In a fourth aspect, the present invention provides a computer readable storage medium having stored thereon computer instructions for causing a computer to perform the simulation method of the chloride ion etching process in the marine structure of the first aspect or any one of its corresponding embodiments.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a simulation method of a chloride ion etching process in a marine structure according to an embodiment of the present invention;
FIG. 2 is a flow chart of a simulation method of a chloride ion etching process in another marine structure according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a distribution of monitoring moments within a preset analog duration according to an embodiment of the present invention;
FIG. 4 is a block diagram of a simulation apparatus of a chloride ion etching process in a marine structure according to an embodiment of the present invention;
fig. 5 is a schematic diagram of a hardware structure of a computer device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In order to solve the problem that the erosion process of chloride ions to concrete simulated by utilizing a fractal derivative model is inaccurate, the embodiment of the invention provides a simulation method of the erosion process of chloride ions in a marine structure, and a inversion parameter set which can enable a simulation result to be more accurate is obtained by inverting model parameters of a one-dimensional variable-order time fractal derivative convection diffusion model under a finite difference discrete format, so that the erosion process of chloride ions to concrete is simulated by utilizing the one-dimensional variable-order time fractal derivative convection diffusion model comprising the inversion parameter set, and the one-dimensional variable-order time fractal derivative convection diffusion model is solved by utilizing the finite difference discrete format, so that the effect of more accurate simulation result is achieved.
According to an embodiment of the present invention there is provided an embodiment of a method of simulating a chloride ion attack process in a marine structure, it being noted that the steps shown in the flow chart of the drawings may be performed in a computer system such as a set of computer executable instructions and, although a logical sequence is shown in the flow chart, in some cases, the steps shown or described may be performed in a different order than that shown herein.
In this embodiment, a method for simulating a chloride ion etching process in a marine structure is provided, which may be used in a computer device, and fig. 1 is a flowchart of a method for simulating a chloride ion etching process in a marine structure according to an embodiment of the present invention, and as shown in fig. 1, the flowchart includes the following steps:
step S101, sea state data of a first sea area and historical erosion data of a second sea area are acquired.
Specifically, the first sea area is a sea area to be simulated of the process of erosion of concrete by chloride ions, and the second sea area and the first sea area have any area with the same or similar sea state data, so that the sea state data of the second sea area and the sea state data of the first sea area belong to the same section. The interval in the present embodiment may be determined by those skilled in the art based on the intrinsic characteristics of sea state data, and is not particularly limited herein.
In particular, the sea state data includes chloride ion concentration, and in addition to the chloride ion concentration, the sea state data includes water cement ratio and PH of the marine structure concrete in the first sea area. The concrete to cement ratio and PH are related to the erosion rate of chloride ions in the concrete.
In particular, historical erosion data, i.e. erosion process of chloride ions in the second sea area over a historical period of time, is data of changes in chloride ion concentration over time at various locations of the marine structure.
And S102, inverting model parameters of a flow diffusion model by utilizing historical erosion data and carrying out one-dimensional variable order time fractal derivative under a finite difference discrete format to obtain an inversion parameter set.
Specifically, a one-dimensional variable-order time fractal derivative convection diffusion model is constructed by adopting the following formula (I):
wherein,fractal derivative for time variant order, +.>Fractal derivative order as a function of time>，/>For convection parameters->For diffusion parameter->Is->Temporal spatial position->Chloride ion concentration at the site.
The finite difference discrete format is constructed using the following equation (two):
wherein,is->Fractional derivative order of corresponding order at each moment,/->Is the firstIndividual temporal spatial position->Concentration of chloride ion at->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->Time steps corresponding to the respective moments. It should be emphasized that the time step in the finite difference discrete format may be obtained according to a uniform point distribution manner, or may be obtained by a preset point distribution manner, and the determining manner of the time step in this embodiment is not particularly limited.
Based on the finite difference discrete format described above, the following equation (three) can be determined:
wherein,is->Personal space position->In->Time->Is used for the concentration of chloride ions,，/>for the length of space->Is the total space segment.
Substituting the formula (II) and the formula (III) into the formula (I) to obtain the one-dimensional variable order time fractal derivative convection diffusion model under the finite difference discrete format.
Specifically, the inversion parameter set includes a variable order fractal derivative order function, a diffusion coefficient, and a convection coefficient.
By using historical erosion data, a least square fitting method is adopted to invert model parameters of a flow diffusion model by one-dimensional variable-order time fractal derivative under a finite difference discrete format, and a variable-order fractal derivative order function, a diffusion coefficient and a convection coefficient are obtained.
Step S103, inputting the chloride ion concentration and inversion parameter set into a one-dimensional variable-order time fractal derivative convection diffusion model, and solving the one-dimensional variable-order time fractal derivative convection diffusion model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area within a preset simulation duration.
Specifically, the chloride ion concentration and inversion parameters of the first sea area, namely model parameters, are input into a one-dimensional variable-order time fractal derivative convection diffusion model, and according to preset time distribution information, the one-dimensional variable-order time fractal derivative convection diffusion model is solved by utilizing a finite difference discrete format, so that the chloride ion concentrations of different positions at all moments are obtained.
According to the simulation method of the chloride ion erosion process in the marine structure, firstly, the model parameters of a one-dimensional variable order time fractal derivative flow diffusion model under a finite difference discrete format are inverted by using the historical erosion data of a second sea area, so that an inversion parameter set is obtained; then, inputting the inversion parameter set and the chloride ion concentration of the first sea area into a one-dimensional variable-order time fractal derivative convection diffusion model, and solving the model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area in a preset simulation duration; the one-dimensional variable-order time fractal derivative flow diffusion model considers the change of the erosion rate of chloride ions to concrete in a long-time process, so that the erosion process is accurately simulated, the solving efficiency and the accuracy of the model are effectively improved on the premise of not increasing space distribution points, and powerful technical support is provided for the corrosion prevention work of marine structures.
In this embodiment, a method for simulating a chloride ion etching process in a marine structure is provided, which may be used in a computer device, and fig. 2 is a flowchart of a method for simulating a chloride ion etching process in a marine structure according to an embodiment of the present invention, and as shown in fig. 2, the flowchart includes the following steps:
step S201, sea state data of a first sea area and historical erosion data of a second sea area are acquired. Please refer to step S101 in the embodiment shown in fig. 1 in detail, which is not described herein.
Step S202, inversion is carried out on model parameters of a flow diffusion model by utilizing historical erosion data and using one-dimensional variable order time fractal derivative under a finite difference discrete format, so as to obtain an inversion parameter set. Please refer to step S102 in the embodiment shown in fig. 1 in detail, which is not described herein.
Step S203, inputting the chloride ion concentration and inversion parameter set into a one-dimensional variable-order time fractal derivative convection diffusion model, and solving the one-dimensional variable-order time fractal derivative convection diffusion model by utilizing a finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area within a preset simulation duration.
Specifically, the step S203 includes:
step S2031, determining a time step based on the preset time distribution information and the variable order fractal derivative order function.
Specifically, the preset time distribution information includes a preset simulation duration and a preset time distribution total number, wherein the preset simulation duration is the total duration of the corrosion process of the chloride ions to be simulated on the concrete, and the preset simulation duration can be determined by a worker according to the characteristics of the marine construction concrete and the concentration of the chloride ions in sea water, and is not particularly limited herein. The total number of preset time distribution points is the total number of monitoring moments, if the preset simulation time length is 3 months, the total number of the preset time distribution points is 180, namely, in the simulation time length of 3 months, the chloride ion concentration at each position of the marine engineering and the building is respectively obtained at the 180 monitoring moments.
Specifically, the time step is the duration between two adjacent monitoring moments, for example, the preset time distribution number isEach monitoring time is ∈ ->、/>、/>、…、/>，/>And->The time length between them, i.e. 1 st time corresponding to time step 1, < >>And->The duration between which is the time step 2 corresponding to the second moment, in this embodiment co +.>A time step. It is emphasized that in this embodiment the time steps are not identical.
Illustratively, in an alternative embodiment, the step of time is determined by equation (four):
wherein,is->Time step corresponding to each moment +.>Is->At all times->Is->At all times->For a preset simulation time length, < >>Is the total number of preset time points->Is->The fractional derivative order of the corresponding order at each moment.
Exemplary, if the preset simulation time is 4 days, the total time distribution is 20, and the initial variable fractal derivative order is 0.6 (initial fractal order is the first time within the preset simulation timeCorresponding variant order fractal derivative order +.>) And (3) obtaining a time step sequence in a preset simulation duration according to the formula (IV), wherein the variable order fractal derivative order of the final moment is 1, and determining each monitoring moment based on the initial moment and the time non-long sequence. The final time step sequence obtained in this embodiment is as follows: i.e. 0.0318877968410468, 0.0398594629012216, 0.0492347099879916, 0.0601503405303135, 0.0727383102717604, 0.0871209676429161, 0.103405047409503, 0.121674164070866, 0.141979498135607, 0.164328289447756, 0.188669623448730, 0.214876773663112, 0.242724953735531, 0.271862531005949, 0.301772078792860, 0.331713840659476, 0.360634527989132, 0.386995720859243, 0.408367913975951, 0.420003448631034, the distribution of monitoring moments in a preset analog period according to the time step sequence is shown in fig. 3, and as can be seen from fig. 3, the dense condition of the monitoring moment distribution gradually decreases with time.
Because the numerical accuracy of the uniform distribution mode is low, the time step determined by the embodiment is continuously changed along with the fractal derivative order of the variable order. Because the lower the order is, the lower the numerical accuracy is, the embodiment improves the numerical accuracy by arranging more monitoring moments at the lower order, thereby improving the accuracy of the simulation process.
Step S2032, solving the one-dimensional variable-order time fractal derivative flow diffusion model by using a finite difference discrete format of taking a time step to obtain chloride ion concentrations respectively corresponding to each moment at different positions.
Specifically, after the time step is obtained, each monitoring time within the preset analog duration is determined. And constructing a finite difference discrete format of the time step by utilizing the time step corresponding to the last moment of the current monitoring moment and the chloride ion concentration corresponding to the last moment of the current monitoring moment at each position, and solving a one-dimensional variable order time fractal derivative flow diffusion model by utilizing the discrete format to obtain the chloride ion concentration of the current monitoring moment at each position, thereby obtaining the chloride ion concentration of the marine construction at different positions at different moments.
Step S2033, obtaining a chloride ion erosion process of the first sea area within a preset simulation duration based on the chloride ion concentrations respectively corresponding to each moment at different positions.
Specifically, the chloride ion concentration corresponding to each position at different time points is displayed in a function or graphic mode, so that the change of the chloride ion concentration at each position along with time can be obtained, the change of the chloride ion concentration at all positions along with time is synthesized, and the corresponding chloride ion erosion process of the first sea area in the preset simulation duration can be obtained.
In this embodiment, a simulation device for a chloride ion erosion process in a marine structure is further provided, and the device is used to implement the foregoing embodiments and preferred embodiments, which have been described and will not be repeated. As used below, the term "module" may be a combination of software and/or hardware that implements a predetermined function. While the means described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware, is also possible and contemplated.
The embodiment provides a simulation device for a chloride ion erosion process in a marine structure, as shown in fig. 4, including:
the acquiring module 401 is configured to acquire sea state data of a first sea area and historical erosion data of a second sea area, where the sea state data of the second sea area and the sea state data of the first sea area belong to the same interval.
And the inversion module 402 is configured to invert model parameters of the flow diffusion model by using historical erosion data and using one-dimensional variable order time fractal derivative under a finite difference discrete format to obtain an inversion parameter set.
The generating module 403 is configured to input the chloride ion concentration and the inversion parameter set into a one-dimensional variable-order time fractal derivative convective diffusion model, and solve the one-dimensional variable-order time fractal derivative convective diffusion model by using a finite difference discrete format based on preset time distribution information, so as to obtain a chloride ion erosion process of the first sea area within a preset simulation duration.
In some alternative embodiments, the inversion parameter set generated by the inversion module 402 includes a variable order fractal derivative order function, a diffusion coefficient, and a convection coefficient.
In some alternative embodiments, the inversion module 402 includes:
and the inversion sub-module is used for inverting the model parameters of the flow diffusion model by adopting a least square fitting method to obtain an inversion parameter set.
In some alternative embodiments, the one-dimensional variable order time fractal derivative in inversion module 402 builds a convective diffusion model using the following formula:
wherein,fractal derivative for time variant order, +.>For the variable order fractal derivative order function,，/>for convection parameters->For diffusion parameter->Is->Temporal spatial position->Chloride ion concentration at the site.
In some alternative embodiments, the generating module 403 includes:
the determining submodule is used for determining the time step based on preset time distribution information and the variable order fractal derivative order function.
And the solving sub-module is used for solving the one-dimensional variable-order time fractal derivative flow diffusion model by utilizing a finite difference discrete format of taking time step length to obtain chloride ion concentrations respectively corresponding to each moment at different positions.
And the determination submodule is used for obtaining the chloride ion erosion process of the first sea area in the preset simulation time based on the chloride ion concentrations respectively corresponding to the moments at different positions.
In some alternative embodiments, the finite difference discrete format in the solution sub-module is:
wherein,is->Fractional derivative order of corresponding order at each moment,/->Is the firstIndividual temporal spatial position->Concentration of chloride ion at->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->Time steps corresponding to the respective moments.
In some alternative embodiments, determining the preset time-point information in the sub-module includes determining a preset analog duration and a preset total number of time-points, and determining the time step in the sub-module is determined by the following formula:
wherein,is->Time step corresponding to each moment +.>Is->At all times->Is->At all times->For a preset simulation time length, < >>Is the total number of preset time points->Is->The fractional derivative order of the corresponding order at each moment.
Further functional descriptions of the above respective modules and units are the same as those of the above corresponding embodiments, and are not repeated here.
The simulation of the chloride ion attack process in the marine structure of the present embodiment is presented in the form of functional units, here referred to as ASIC (Application Specific Integrated Circuit ) circuits, processors and memories executing one or more software or fixed programs, and/or other devices that can provide the above described functionality.
The embodiment of the invention also provides computer equipment, which is provided with the simulation device for the chloride ion erosion process in the marine structure shown in the figure 4.
Referring to fig. 5, fig. 5 is a schematic structural diagram of a computer device according to an alternative embodiment of the present invention, as shown in fig. 5, the computer device includes: one or more processors 10, memory 20, and interfaces for connecting the various components, including high-speed interfaces and low-speed interfaces. The various components are communicatively coupled to each other using different buses and may be mounted on a common motherboard or in other manners as desired. The processor may process instructions executing within the computer device, including instructions stored in or on memory to display graphical information of the GUI on an external input/output device, such as a display device coupled to the interface. In some alternative embodiments, multiple processors and/or multiple buses may be used, if desired, along with multiple memories and multiple memories. Also, multiple computer devices may be connected, each providing a portion of the necessary operations (e.g., as a server array, a set of blade servers, or a multiprocessor system). One processor 10 is illustrated in fig. 5.
The processor 10 may be a central processor, a network processor, or a combination thereof. The processor 10 may further include a hardware chip, among others. The hardware chip may be an application specific integrated circuit, a programmable logic device, or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general-purpose array logic, or any combination thereof.
Wherein the memory 20 stores instructions executable by the at least one processor 10 to cause the at least one processor 10 to perform the methods shown in implementing the above embodiments.
The memory 20 may include a storage program area that may store an operating system, at least one application program required for functions, and a storage data area; the storage data area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include high-speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid-state storage device. In some alternative embodiments, memory 20 may optionally include memory located remotely from processor 10, which may be connected to the computer device via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
Memory 20 may include volatile memory, such as random access memory; the memory may also include non-volatile memory, such as flash memory, hard disk, or solid state disk; the memory 20 may also comprise a combination of the above types of memories.
The computer device also includes a communication interface 30 for the computer device to communicate with other devices or communication networks.
The embodiments of the present invention also provide a computer readable storage medium, and the method according to the embodiments of the present invention described above may be implemented in hardware, firmware, or as a computer code which may be recorded on a storage medium, or as original stored in a remote storage medium or a non-transitory machine readable storage medium downloaded through a network and to be stored in a local storage medium, so that the method described herein may be stored on such software process on a storage medium using a general purpose computer, a special purpose processor, or programmable or special purpose hardware. The storage medium can be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, a solid state disk or the like; further, the storage medium may also comprise a combination of memories of the kind described above. It will be appreciated that a computer, processor, microprocessor controller or programmable hardware includes a storage element that can store or receive software or computer code that, when accessed and executed by the computer, processor or hardware, implements the methods illustrated by the above embodiments.
Although embodiments of the present invention have been described in connection with the accompanying drawings, various modifications and variations may be made by those skilled in the art without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope of the invention as defined by the appended claims.
Claims (6)
1. A method of simulating a chloride ion erosion process in a marine structure, the method comprising:
acquiring sea state data of a first sea area and historical erosion data of a second sea area, wherein the sea state data comprises chloride ion concentration, and the sea state data of the second sea area and the sea state data of the first sea area belong to the same interval;
inverting model parameters of a flow diffusion model by utilizing the historical erosion data and carrying out one-dimensional variable order time fractal derivative under a finite difference discrete format to obtain an inversion parameter set;
inputting the chloride ion concentration and the inversion parameter set into the one-dimensional variable-order time fractal derivative convective diffusion model, and solving the one-dimensional variable-order time fractal derivative convective diffusion model by utilizing the finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area in a preset simulation duration; the inversion parameter set comprises a variable-order fractal derivative order function, a diffusion coefficient and a convection coefficient;
inputting the chloride ion concentration and the inversion parameter set into the one-dimensional variable-order time fractal derivative convective diffusion model, and solving the one-dimensional variable-order time fractal derivative convective diffusion model by utilizing the finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area within a preset simulation duration, wherein the method comprises the following steps:
determining a time step based on the preset time distribution point information and a variable order fractal derivative order function;
solving a flow diffusion model by using the finite difference discrete format of the time step, and obtaining chloride ion concentrations respectively corresponding to each moment at different positions;
obtaining a chloride ion erosion process of the first sea area in a preset simulation time based on the chloride ion concentrations respectively corresponding to each moment at different positions;
the one-dimensional variable-order time fractal derivative convection diffusion model is constructed by adopting the following formula:
wherein,fractal derivative for time variant order, +.>Fractal derivative order function for variable order, +.>，/>For convection parameters->For diffusion parameter->Is->Temporal spatial position->Chloride ion concentration at the site;
the finite difference discrete format is:
wherein,is->Fractional derivative order of corresponding order at each moment,/->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->A time step corresponding to each moment;
the preset time distribution information comprises preset simulation duration and preset time distribution total number, and the time step is determined by the following formula:
wherein,is->Time step corresponding to each moment +.>Is->At all times->Is->At all times->For a preset simulation time length, < >>Is the total number of preset time points->Is->The fractional derivative order of the corresponding order at each moment.
2. The method of claim 1, wherein inverting model parameters of the flow diffusion model with one-dimensional variable order time fractal derivative in finite difference discrete format to obtain an inversion parameter set, comprises:
and inverting model parameters of the flow diffusion model by adopting a least square fitting method to obtain the inversion parameter set.
3. A simulation device for a chloride ion etching process in a marine structure, the device comprising:
the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is used for acquiring sea state data of a first sea area and historical erosion data of a second sea area, the sea state data comprise chloride ion concentration, and the sea state data of the second sea area and the sea state data of the first sea area belong to the same interval;
the inversion module is used for inverting model parameters of the flow diffusion model by utilizing the historical erosion data and carrying out inversion on one-dimensional variable order time fractal derivative under a finite difference discrete format to obtain an inversion parameter set;
the generation module is used for inputting the chloride ion concentration and the inversion parameter set into the one-dimensional variable-order time fractal derivative convection diffusion model, and solving the one-dimensional variable-order time fractal derivative convection diffusion model by utilizing the finite difference discrete format based on preset time distribution information to obtain a chloride ion erosion process of the first sea area within a preset simulation duration;
the inversion parameter set generated by the inversion module comprises a variable order fractal derivative order function, a diffusion coefficient and a convection coefficient;
the generating module comprises:
the determining submodule is used for determining a time step based on the preset time distribution information and the variable-order fractal derivative order function;
the solving submodule is used for solving the one-dimensional variable-order time fractal derivative flow diffusion model by utilizing the finite difference discrete format of the time step to obtain chloride ion concentrations respectively corresponding to each moment at different positions;
the determining submodule is used for obtaining a chloride ion erosion process of the first sea area in a preset simulation duration based on the chloride ion concentrations respectively corresponding to each moment at different positions;
the one-dimensional variable-order time fractal derivative convection diffusion model is constructed by adopting the following formula:
wherein,fractal derivative for time variant order, +.>Fractal derivative order function for variable order, +.>，/>For convection parameters->For diffusion parameter->Is->Temporal spatial position->Chloride ion concentration at the site;
the finite difference discrete format is:
wherein,is->Fractional derivative order of corresponding order at each moment,/->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->Individual temporal spatial position->Concentration of chloride ion at->Is->A time step corresponding to each moment;
the preset time distribution information comprises preset simulation duration and preset time distribution total number, and the time step is determined by the following formula:
wherein,is->Time step corresponding to each moment +.>Is->At all times->Is->At all times->For a preset simulation time length, < >>Is the total number of preset time points->Is->The fractional derivative order of the corresponding order at each moment.
4. The apparatus of claim 3, wherein the inversion module comprises:
and the inversion sub-module is used for inverting the model parameters of the flow diffusion model by adopting a least square fitting method to obtain the inversion parameter set.
5. A computer device, comprising:
a memory and a processor, said memory and said processor being communicatively coupled to each other, said memory having stored therein computer instructions, said processor executing said computer instructions to perform the method of simulating a chloride ion etching process in a marine structure according to claim 1 or 2.
6. A computer readable storage medium having stored thereon computer instructions for causing a computer to perform the simulation of a chloride ion etching process in a marine structure according to claim 1 or 2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311415094.4A CN117150867B (en) | 2023-10-30 | 2023-10-30 | Simulation method, device and equipment for chloride ion erosion process in marine building |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311415094.4A CN117150867B (en) | 2023-10-30 | 2023-10-30 | Simulation method, device and equipment for chloride ion erosion process in marine building |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117150867A CN117150867A (en) | 2023-12-01 |
| CN117150867B true CN117150867B (en) | 2024-01-26 |
Family
ID=88908405
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311415094.4A Active CN117150867B (en) | 2023-10-30 | 2023-10-30 | Simulation method, device and equipment for chloride ion erosion process in marine building |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117150867B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108088770A (en) * | 2017-11-23 | 2018-05-29 | 河海大学 | Chlorion change exponent number of unusual diffusion in concrete divides shape derivative analogue method |
| CN112033881A (en) * | 2020-08-11 | 2020-12-04 | 南京理工大学 | Method for calculating chloride ion concentration in concrete under action of chloride salt-corrosion |
| CN113919115A (en) * | 2020-07-08 | 2022-01-11 | 中核武汉核电运行技术股份有限公司 | Method for establishing existing concrete chloride ion erosion life prediction model |
| CN115712993A (en) * | 2022-11-08 | 2023-02-24 | 中国长江三峡集团有限公司 | Chloride ion erosion model construction method and chloride ion erosion simulation method |
| CN116312903A (en) * | 2023-05-25 | 2023-06-23 | 长江三峡集团实业发展(北京)有限公司 | Sea sand concrete chloride ion erosion analysis method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220383990A1 (en) * | 2021-05-28 | 2022-12-01 | DigiM Solution LLC | System and Method to Predict Mass Transport from Complex Release Systems Using Experimental Data-based Modeling |
-
2023
- 2023-10-30 CN CN202311415094.4A patent/CN117150867B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108088770A (en) * | 2017-11-23 | 2018-05-29 | 河海大学 | Chlorion change exponent number of unusual diffusion in concrete divides shape derivative analogue method |
| CN113919115A (en) * | 2020-07-08 | 2022-01-11 | 中核武汉核电运行技术股份有限公司 | Method for establishing existing concrete chloride ion erosion life prediction model |
| CN112033881A (en) * | 2020-08-11 | 2020-12-04 | 南京理工大学 | Method for calculating chloride ion concentration in concrete under action of chloride salt-corrosion |
| CN115712993A (en) * | 2022-11-08 | 2023-02-24 | 中国长江三峡集团有限公司 | Chloride ion erosion model construction method and chloride ion erosion simulation method |
| CN116312903A (en) * | 2023-05-25 | 2023-06-23 | 长江三峡集团实业发展(北京)有限公司 | Sea sand concrete chloride ion erosion analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117150867A (en) | 2023-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777302B (en) | Method and device for converting space geographic coordinates | |
| CN113792642B (en) | River and lake ecological management data processing method and system based on intelligent technology | |
| Smith et al. | Towards a hydrodynamic modelling framework appropriate for applications in urban flood assessment and mitigation using heterogeneous computing | |
| WO2024098951A1 (en) | Method and apparatus for simulating rapid transport of sediment during reservoir discharging | |
| CN116167289B (en) | Power grid operation scene generation method and device, computer equipment and storage medium | |
| CN114861473B (en) | Effective storage capacity analysis method and device for river channel type reservoir and electronic equipment | |
| CN117150867B (en) | Simulation method, device and equipment for chloride ion erosion process in marine building | |
| CN114911788A (en) | Data interpolation method and device and storage medium | |
| CN114065593A (en) | Method and device for underground water coupling simulation in multi-scale fracture matrix system | |
| Defne et al. | Toward a comprehensive water-quality modeling of Barnegat Bay: development of ROMS to WASP coupler | |
| CN116720448B (en) | Wind power generation random simulation method, device, equipment and medium | |
| CN115712993A (en) | Chloride ion erosion model construction method and chloride ion erosion simulation method | |
| CN114841077A (en) | Wind power prediction method, device and medium | |
| CN114242181A (en) | Desert sand concrete strength prediction model training method, device, equipment and medium | |
| CN117150975B (en) | Hydrodynamic model parameter optimization and hydrodynamic process simulation method and device | |
| Branisavljević et al. | Uncertainty reduction in water distribution network modelling using system inflow data | |
| CN117172162B (en) | Simulation method and device for saline solution migration process in seawater migration process | |
| CN110795791A (en) | Engineering model slicing method and device, computer equipment and storage medium | |
| CN110288256B (en) | Scene analysis method and device | |
| Raschke | Opportunities of the minimum Anderson–Darling estimator as a variant of the maximum likelihood method | |
| CN117933729B (en) | Cultural heritage vulnerability assessment method and device under flood influence | |
| CN110569404A (en) | data comparison method and device, electronic equipment and storage medium | |
| CN116757003B (en) | Plain river network simulation method, plain river network simulation device, computer equipment and medium | |
| CN116029151B (en) | Water flow resistance prediction method, training method, flow prediction method and device | |
| CN114923544B (en) | Method and device for predicting water level of outlet section and electronic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |