1. A method for correcting simulation data in civil engineering, characterized in that it comprises the following steps:

s1: acquiring geometric parameters and material parameters of a measured object, establishing a finite element model of the measured object, and dividing a grid;

s2: selecting a plurality of measuring points in the finite element model, and respectively constructing a discrete weight function of each measuring point;

s3: applying real boundary conditions to the finite element model, and calculating to obtain ideal data of each measuring point;

s4: correcting the ideal data of each measuring point by using the discrete weight function and the measured data of each measuring point, so as to obtain the corrected data of each node in the finite element model of the measuring object;

the order of the above steps S2 and S3 is not fixed.

2. The method of claim 1, wherein the step of constructing the discrete weight function of each measurement point in S2 specifically comprises:

s21: applying loads of any physical quantity type of a unit I at a measuring point i in the finite element model, and applying constraints of the same physical quantity type to other measuring points; i is a positive integer and is used for representing the serial number of the measuring point;

s32: obtaining simulation values of the physical quantity types of all grid nodes in the finite element model in a numerical simulation mode, namely obtaining a discrete weight function of the measuring point i;

s33: and repeating the steps S31-S32 until the discrete weight functions of all the measuring points are obtained.

3. The method of correcting simulation data in civil engineering work according to claim 2, wherein the discrete weighting function is specifically:

N_i＝{u₁,u₂,...,u_n}；

wherein N is_iIs a discrete weight function of the measuring point i; u. of₁,u₂,...,u_nIdeal data of the 1 st node, the 2 nd node to the nth node in the finite element model are obtained; and n is a positive integer, and the total number of the grid nodes in the finite element model is taken.

4. The method of correcting simulation data in civil engineering work according to claim 1, wherein said step S4 specifically includes:

the measured data u of each measuring point_i ^rSubstitution formulaThe correction data of the type of the measured physical quantity can be calculated;

wherein m is the total number of the measuring points; i is the serial number of the measuring point; n is a radical of_iIs a discrete weight function of the measuring point i;representing a hadamard product; u. of^eA corrected data matrix for the type of physical quantity measured; u. of^sAn ideal data matrix of the measured physical quantity type; u. of_i ^rMeasured data of a measuring point i is obtained; u. of_i ^sIs the ideal data for point i.

5. The method of correcting simulation data in civil engineering work of claim 4, wherein the ideal data u is associated with a certain measure point_i ^sWhen it is equal to 0, then let u_i ^r/u_i ^sEqual to 1.

6. The method of correcting simulation data in civil engineering work according to claim 1, wherein the plurality of measurement points are located on mesh nodes of the finite element model.

7. The method of correcting simulation data in civil engineering work of claim 1, wherein the discrete weight function of the same sample point is changed in accordance with a change in the type of the physical quantity of the correction data.

8. The simulation data correcting method in civil engineering according to claim 1, wherein the type of the physical quantity of the correction data includes any one of temperature, speed, concentration, stress, strain, displacement.

9. A simulation data correction system in civil engineering, characterized by comprising:

a processor and a memory for storing executable instructions;

wherein the processor is configured to execute the executable instructions to perform the simulation data correction method in civil engineering as claimed in any one of claims 1 to 8.

10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements a simulation data correcting method in civil engineering as claimed in any one of claims 1 to 8.

Background

Structural data in civil engineering are often difficult to determine, the result of numerical calculation is greatly influenced, the measurement of the structural data in a laboratory has the problem of scale effect, the number of field samples is often not large in consideration of economic cost, and the real structural data of the whole engineering area cannot be obtained, so that a method of space difference is necessary to be adopted to correct laboratory simulation data.

The estimation accuracy of the spatial interpolation method depends mainly on whether all sample information is fully and reasonably utilized. The general spatial interpolation method only focuses on the spatial distribution and the numerical size of the sample and ignores hidden information behind the sample, resulting in insufficient interpolation precision. For example, conventional methods do not distinguish between processing the sample data for distortion and temperature. In addition, common spatial interpolation methods usually assume some mathematical models to supplement spatial correlation information, but their determination is empirical and fuzzy, and interpolation accuracy is difficult to guarantee.

Disclosure of Invention

In view of the above, the invention provides a method and a system for correcting simulation data in civil engineering, which can reasonably expand data of any physical quantity type of a plurality of measuring points to the whole structure, realize a high-precision spatial interpolation effect, and provide direct structure working state judgment basis for scientific research personnel and engineering experts.

In order to achieve the purpose, the invention adopts the technical scheme that:

according to a first aspect of the present invention, there is provided a simulation data correcting method in civil engineering, the method comprising the steps of:

s1: acquiring geometric parameters and material parameters of a measured object, establishing a finite element model of the measured object, and dividing a grid;

s2: selecting a plurality of measuring points in the finite element model, and respectively constructing a discrete weight function of each measuring point;

s3: applying real boundary conditions to the finite element model, and calculating to obtain ideal data of each measuring point;

s4: correcting the ideal data of each measuring point by using the discrete weight function and the measured data of each measuring point, so as to obtain the corrected data of each node in the finite element model of the measuring object;

the order of the above steps S2 and S3 is not fixed.

Further, the constructing of the discrete weight function of each measurement point in S2 specifically includes:

s21: applying loads of any physical quantity type of a unit I at a measuring point i in the finite element model, and applying constraints of the same physical quantity type to other measuring points; i is a positive integer and is used for representing the serial number of the measuring point;

s32: obtaining simulation values of the physical quantity types of all grid nodes in the finite element model in a numerical simulation mode, namely obtaining a discrete weight function of the measuring point i;

s33: and repeating the steps S31-S32 until the discrete weight functions of all the measuring points are obtained.

Further, the discrete weight function is specifically:

N_i＝{u₁,u₂,...,u_n}；

wherein N is_iIs a discrete weight function of the measuring point i; u. of₁,u₂,...,u_nIdeal data of the 1 st node, the 2 nd node to the nth node in the finite element model are obtained; and n is a positive integer, and the total number of the grid nodes in the finite element model is taken.

Further, the step S4 specifically includes:

the measured data u of each measuring point_i ^rSubstitution formulaThe correction data of the type of the measured physical quantity can be calculated;

wherein m is the total number of the measuring points; i is the serial number of the measuring point; n is a radical of_iIs a discrete weight function of the measuring point i;representing a hadamard product; u. of^eA corrected data matrix for the type of physical quantity measured; u. of^sAn ideal data matrix of the measured physical quantity type; u. of_i ^rMeasured data of a measuring point i is obtained; u. of_i ^sIs the ideal data for point i.

Furthermore, if the ideal data u corresponding to a certain measuring point_i ^sWhen it is equal to 0, then let u_i ^r/u_i ^sEqual to 1.

Furthermore, the measuring points are all located on the grid nodes of the finite element model.

Further, the discrete weight function of the same sample point is changed with the change of the type of the physical quantity of the correction data.

Further, the type of the physical quantity of the correction data includes any one of temperature, velocity, concentration, stress, strain, and displacement.

According to a second aspect of the present invention, there is provided a simulation data correcting system in civil engineering, comprising:

a processor and a memory for storing executable instructions;

wherein the processor is configured to execute the executable instructions to perform the simulation data correction method in civil engineering described above.

According to a third aspect of the invention, there is provided a computer-readable storage medium on which a computer program is stored, characterized in that the computer program, when executed by a processor, implements the above-described simulation data correction method in civil engineering.

Compared with the prior art, the civil engineering measurement data correction method and system provided by the invention have the following advantages:

(1) the data correction method provided by the invention effectively combines experimental data and simulation data, so that a comprehensive and real working state of an engineering structure is obtained, the implementation process is simple, the number and distribution of actual measuring points in the engineering are not limited, equations and optimization problems are not solved, and the method is convenient to popularize and apply;

(2) the invention solves the problem that the simulation result is inconsistent with the experimental result due to the unknown defect of the structure by correcting the simulation result obtained by simulation through part of the measured data in the engineering, so that the simulation result can more truly reflect the actual stress state of the structure.

(3) The method reasonably expands the data of any physical quantity type of a plurality of measuring points to the whole structure, realizes a high-precision spatial interpolation effect, and provides direct structural working state judgment basis for scientific research personnel and engineering experts.

Drawings

The accompanying drawings are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention.

FIG. 1 is a finite element model of a three-span bridge according to embodiment 1 of the present invention;

FIG. 2 is a real and ideal simulation of the deflection of a three-span bridge according to example 1 of the present invention;

FIG. 3 is a schematic diagram of a discrete weight function of the deflection of a three-span bridge according to embodiment 1 of the present invention;

FIG. 4 is a comparison of the simulation adjustment results of the deflection of the three-span bridge according to the embodiment 1 of the present invention;

FIG. 5 is a finite element model of the exterior wall of the building according to embodiment 2 of the present invention;

FIG. 6 is a real and ideal simulation of the temperature of the exterior wall of the building as described in example 2 of the present invention;

FIG. 7 is a graph illustrating a discrete weighting function of the temperature of the exterior wall of the building according to example 2 of the present invention;

FIG. 8 is a comparison of the simulation results of adjusting the temperature of the exterior wall of the building according to example 2 of the present invention;

Detailed Description

Reference will now be made in detail to the exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present invention. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the invention, as detailed in the appended claims.

The terms first, second and the like in the description and in the claims of the present invention are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are, for example, capable of operation in sequences other than those illustrated or otherwise described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.

A plurality, including two or more.

And/or, it should be understood that, as used herein, the term "and/or" is merely one type of association that describes an associated object, meaning that three types of relationships may exist. For example, a and/or B, may represent: a exists alone, A and B exist simultaneously, and B exists alone.

A simulation data correction method in civil engineering comprises the following specific technical scheme:

a. establishing a finite element model of a measuring object, dividing the finite element model into a plurality of units, applying specific boundary conditions, and obtaining a discrete weight function N of each measuring point by a numerical simulation method_i。

b. Performing conventional simulation, applying real boundary conditions, inputting material attributes, and obtaining a conventional simulation result u^sAnd extracting a conventional analog value u at the measurement point_i ^s。

c. The measured data u of each measuring point_i ^rSubstitution formulaI.e. an accurate estimate u of the measured physical quantity field can be calculated^eWherein m is the total number of the measuring points,representing the hadamard product.

Preferably, all the measurement points in step a should be exactly on the mesh nodes of the model.

Preferably, the specific boundary condition in step a refers to a weight function N in calculating the ith measurement point_iAt point i, the physical quantity in unity (which is identical or close to the data type) is added, and zero is set at the remaining points.

Preferably, the numerical simulation in step a only considers the linear elasticity of the material, and the value of each node physical quantity obtained by the simulation is the discrete weight function N forming the measuring point i_i＝{u₁,u₂,...,u_j,...,u_nIn which u_jIs the analog value at the jth node, and n is the total number of nodes.

Preferably, u in step c^e、u^sAnd N_iAll are vectors related to the coordinates of the nodes, and the jth element is the exact correction value, ideal simulation value and weight function value at the jth node.

Preferably, the analog value u corresponding to the measured point in step c_i ^sWhen it is 0, the coefficient u is set_i ^r/u_i ^sIs 1.

Preferably, the automatic loop-constrained program is programmed in ANSYS-APDL language and the discrete weight function vectors are output in bulk.

Preferably, the present invention is particularly suitable for physical quantities having transfer characteristics such as deformation, velocity, temperature, concentration, etc., which, when the system deviates from the equilibrium state, a transfer phenomenon occurs to bring the system toward the equilibrium state. For physical quantities that are not transitive, a weighting function can be calculated with similar physical quantities, for example, a weighting function constructed with displacements to adjust the strain data.

Example 1: simulation result adjustment of three-span beam

The three-span beam is a common structure in civil engineering, and is generally characterized in that two ends are fixedly connected, and two sliding supports are arranged in the middle. The geometrical parameters of the three-span beam structure at a certain position are measured as follows, the length of the beam is 3m, the length of each span is 1m, and the section of the beam is 100 multiplied by 100mm²Simultaneously measuring the elastic modulus E of the three-span beam material₀200 GPa. From this a finite element model as shown in fig. 1 was built.

The three-span beam is applied with uniform load of 1kN/m, if a damaged part of the real beam is found, the elastic modulus of the real beam is 1/20 which is normal, and the ideal simulation does not consider the weak part, and the simulation results of the two parts are shown in figure 2.

Extracting real simulation displacement values of the three spans as actual measurement samples u_i ^rAnd extracting the corresponding ideal analog value u_i ^s. In fact, the displacement of the four supports is also known information and can be used as a measuring point for supplement, and the precision of analog adjustment is improved.

Then, discrete weight functions of seven measuring points are calculated, and the weight function at the measuring point 4 is constructed as an example. According to the constructive nature of the weighting function, add y-direction unit one at the measuring point 4And (4) displacement, namely adding fixed constraint in the y direction to other measuring points, and limiting rigid body displacement in other directions of the plate. Then, static analysis is carried out, and the obtained y-direction displacement field is the discrete weight function N₄＝{u₁,u₂,...,u_k,...,u_nAs shown in fig. 3. Similarly, other discrete weighting functions N may be obtained_i(i-1, 2, …,7), wherein the weighting function N is₇A schematic diagram is also given in fig. 3.

After obtaining the discrete weight function of each measuring point, the ideal simulation result u is obtained^sMeasured value u of the measuring point_i ^r(from a real simulation) and a simulated value u_i ^sSubstitution formulaThe adjustment result of the beam deflection simulation can be obtained, wherein u is at four branch points_i ^r/u_i ^sIs set to 1. FIG. 4 compares the results of the real simulation, the ideal simulation, and the simulation adjustment, and it can be seen that the ideal simulation is adjusted to be closer to the real simulation, illustrating the effectiveness of the method.

Example 2: simulation result adjustment of wall temperature

Measuring the wall size of a section of wall of a certain building to be 300 multiplied by 200mm²The temperature outside the wall is 50 ℃, and the convection heat transfer coefficient inside the wall is 0.5. The thermal conductivity of the normal insulation board is 0.4W/m.K (material 2), and the thermal conductivity of the concrete is 1.28W/m.K (material 2). A finite element model as shown in fig. 5 was built and divided into 60 x 40 Shell 132 cells.

If the heat insulation board of the real outer wall is seeped water at this time, the heat conductivity coefficient of the real outer wall becomes 4W/m.K, and the damage of the heat insulation board is not considered in the ideal simulation, and the simulation results of the heat insulation board and the ideal simulation are respectively shown as a and b in fig. 6.

Three real simulation temperature values are extracted as actual measurement samples u_i ^rAnd extracting the corresponding ideal analog value u_i ^s. In fact, the temperature of the outer wall boundary is also known information, and the whole is regarded as a measuring point.

Discrete weights for four measurement points are then calculatedThe function is exemplified by constructing a weight function at measurement point 1. According to the constructive nature of the weighting function, a temperature of 1 ℃ is added at point 1, and 0 ℃ is added at the other points. Then carrying out steady-state thermodynamic analysis to obtain a temperature field which is a discrete weight function N₁＝{u₁,u₂,...,u_k,...,u_nAs shown in fig. 7 a. Similarly, other discrete weighting functions N may be obtained_i(i-1, 2, …,4), wherein the weighting function N is₂The schematic diagram is given as shown in fig. 7 b.

After obtaining the discrete weight function of each measuring point, the ideal simulation result u is obtained^sMeasured value u of the measuring point_i ^r(from a real simulation) and a simulated value u_i ^sSubstitution formulaThe adjustment results for beam deflection simulation can be obtained as shown in figure 8 a. Fig. 8b compares the results of the real simulation, the ideal simulation, and the simulation adjustment at the central axis of the wall, and it can be seen that the ideal simulation is adjusted to be closer to the real simulation.

The above-mentioned serial numbers of the embodiments of the present invention are merely for description and do not represent the merits of the embodiments.

While the present invention has been described with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, which are illustrative and not restrictive, and it will be apparent to those skilled in the art that various changes and modifications can be made therein without departing from the spirit and scope of the invention as defined in the appended claims.