Method and system for comparative analysis of dam structure safety performance
The VTK-based method and system enable efficient and accurate comparison of finite element and monitoring data for dam safety analysis by interactive visualization and curve fitting, addressing inefficiencies in existing methods.
Patent Information
- Application Number
- CN202310038114.4
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-01-10
- Publication Date
- 2025-07-15
- Estimated Expiration
- 2043-01-10
AI Technical Summary
In the prior art, the comparison and analysis method of finite element calculation results and on-site monitoring data is time-consuming, labor-intensive, and error-prone, and cannot effectively improve the accuracy and efficiency of the comparison and analysis of the safety state of the dam.
Using three-dimensional visualization technology based on VTK library, through inverse distance weighted interpolation and quadratic B-spline fitting algorithm, interactive comparison and analysis of finite element calculation results and safety monitoring data is realized, including model construction, section cutting, data interpolation and curve fitting, improving the visualization and accuracy of data comparison.
It realizes efficient and accurate comparison and analysis of the safety status of the dam structure, reduces manual operation, improves comparison efficiency and accuracy, and enhances the visualization effect of the comparison results.
Smart Images

Figure CN116205099B_ABST
Abstract
Description
Technical Field
[0001] This application relates to the technical field of hydropower engineering, and particularly to a method and system for comparative analysis of the safety behavior of dam structures. Background Art
[0002] The safe operation of a hydropower station dam is related to the safe production and economic benefits of the power plant itself, and also to the safety of the lives and property of the people downstream of the dam. However, with the long-term operation of the dam, affected by external operating environments and the deterioration of its own materials and other factors, it is difficult for the dam to completely avoid local or partial safety hazards. How to effectively predict and analyze and give early warnings about the safe operation problems of the dam has become the key to effectively dealing with such problems.
[0003] Finite element calculation is a very useful means for analyzing and simulating the operating conditions of dams. To ensure the accuracy of the finite element calculation results, it is necessary to compare and analyze the finite element calculation results with on-site monitoring data. Since most finite element calculations are carried out using commercial software, while on-site monitoring data is stored in an independent acquisition system, and workers are used to extract data for the analysis and comparison of the two, the existing comparison methods are time-consuming, laborious, inefficient, the obtained results are not intuitive enough, and errors and omissions cannot be completely avoided.
[0004] To improve the accuracy and efficiency of the comparative analysis of the finite element calculation and the measured results of the dam, it is necessary to study a simple, fast and commercially software-unrestricted system and method for comparative analysis of the safety behavior of the dam. Summary of the Invention
[0005] This application provides a method and system for comparative analysis of the safety behavior of dam structures, which can conveniently, quickly and accurately compare and analyze the finite element calculation results of the dam with the safety monitoring data.
[0006] This application provides a method for comparative analysis of the safety behavior of dam structures, including the following steps:
[0007] Step S11: Construct using tetrahedral element vtkTetra, pentahedral element vtkHexahedron, and hexahedral element vtkHexahedron respectively according to the finite element calculation model element type, and use vtkSTLReader to read the three-dimensional model of the safety monitoring instrument;
[0008] Step S12: Cut the finite element calculation model using the vtkCutter shear function in VTK according to the safety monitoring section position, and then use the vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map, and obtain the measured data through interactive clicking operations on the three-dimensional model of the safety monitoring instrument on the section;
[0009] Step S13: According to the information of the instrument coordinate point P(x, y, z), search for the nearest calculation unit M. The types of M include tetrahedral elements, pentahedral elements, and hexahedral elements. A tetrahedral element has 4 nodes, a pentahedral element has 6 nodes, and a hexahedral element has 8 nodes. The nodes of the element M can be represented as P i (x i ,y i ,z i ). Use the inverse distance weighted interpolation algorithm to perform interpolation calculation on the calculated values of the element nodes to obtain the finite element calculated value s of the instrument position;
[0010] The inverse distance weighted interpolation function is
[0011]
[0012] wherein, is the distance from the point P(x, y, z) to the point P i (x i ,y i ,z i ); v represents the scalar value of the finite element calculation node, which can be the displacement, settlement, etc. of the dam; i = 1, 2,..., n, and k is a constant greater than 0, that is, the weighted power exponent, which is determined according to the actual interpolation effect;
[0013] Step S14: Interactively click on the finite element calculation section cloud map to extract the finite element analysis value s of the measuring point position of the dam under different time points in a certain period i and the measured value r i , i = 1, 2,..., n, and n is the number of time points in this period for comparison.
[0014] Preferably, it further includes the following steps:
[0015] Step S20: Use the curve fitting algorithm to respectively fit and display the multiple finite element analysis values s i of the measuring points and the measured value r i .
[0016] Preferably, step S20 includes the following steps:
[0017] Step S21: Use the quadratic B-spline curve fitting algorithm to respectively perform curve fitting on the finite element analysis value and the measured value of the measuring point position to obtain the finite element value curve as S and the measured value curve as R;
[0018] The quadratic B-spline curve function used is:
[0019]
[0020] wherein, A1 = x1 - x0, B1 = y1 - y0,
[0021] Step S22: Display and output the obtained finite element analysis value curve S and the measured value curve R on the same coordinate system graph.
[0022] Preferably, in step S21, every 3 adjacent points after the quadratic B-spline curve fitting are used to draw a parabola according to formula (2), and the front and back curves at the endpoints of each segment are tangent to each other.
[0023] On the other hand, the present application also provides a dam structure safety state comparison and analysis system, including:
[0024] Instruction receiving module: used to receive user operation instructions. The user selects the finite element calculation result file on the device touch display screen. The finite element calculation result file includes: model element node information, displacement information, etc., and touches the OK button to import.
[0025] Finite element calculation data reading module: used to detect whether there is an event of the user triggering the model data import device. If it is detected that there is an event of the user triggering the import device, read the model data path and save the model data in the device memory.
[0026] 3D model construction module: used to respectively establish a 3D visualization model of the finite element calculation result and the safety monitoring instrument in the VTK visualization scene. According to the finite element calculation model element type, tetrahedron element vtkTetra, pentahedron element vtkHexahedron, and hexahedron element vtkHexahedron are respectively used for construction, and the 3D model of the safety monitoring instrument is read using vtkSTLReader.
[0027] Monitoring section and measured data acquisition module: used to cut the finite element calculation model according to the safety monitoring section position using the vtkClipPolyData cutting function in VTK, and then use the vtkCutter shear function and vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map. Interactive clicking operations on the 3D model of the safety monitoring instrument on the section are used to obtain the measured data, and the data is stored in the device memory.
[0028] Calculated value search module: used to search for the finite element calculated value at the safety monitoring point position. According to the instrument coordinate information, search for the nearest calculation unit, and use the interpolation algorithm to interpolate and calculate the calculated value of the unit node to obtain the finite element calculated value at the instrument position, and the data is stored in the device memory.
[0029] Single - value comparison analysis module: It is used to obtain the finite - element calculation value and the monitoring value of the monitoring point while interactively clicking on the cross - section cloud map of the finite - element calculation, extract the finite - element analysis value and the measured value at the measuring point position under the required comparison time conditions within the comparison period of the dam, and store the data in the device memory.
[0030] Preferably, the interpolation algorithm used in the calculated - value search module is the inverse - distance weighted interpolation method.
[0031] Preferably, it further includes:
[0032] Curve fitting module: It is used to perform curve fitting on the finite - element analysis value and the measured value at the measuring point position respectively by using the curve fitting algorithm to obtain the fitting curve.
[0033] Curve comparison and plotting module: It is used to display the fitting curve and the measured - value curve.
[0034] Preferably, the curve fitting algorithm used in the curve fitting module is the quadratic B - spline curve.
[0035] The beneficial effects that this application can produce include:
[0036] 1) The method for comparing the safety state of the dam structure provided by this application, after performing three - dimensional visualization modeling on the finite - element calculation analysis results and the safety detection instrument setting position data by using the data of the visualization tool function library (VTK) class library configured on the Windows platform, obtaining the monitoring cross - section data and the measured data, and then searching for the finite - element calculation value at the instrument position by using the difference method, the operator only needs to click on each monitoring point on the three - dimensional visual model to compare the finite - element calculation value and the monitoring value of this point. The comparison result is directly visible, without manual comparison, with high comparison efficiency and accuracy, effectively solving the problems of the existing data comparison methods.
[0037] 2) The method for comparing the safety state of the dam structure provided by this application, after performing curve fitting on the finite - element analysis value and the measured value by using curve fitting, can also display the fitting comparison result on the model image, enhancing the visualization degree of the comparison effect and directly reading the comparison quantization result.
[0038] 3) The method for comparing the safety state of the dam structure provided by this application can directly load the finite - element calculation model and the three - dimensional safety monitoring model, arbitrarily cut the calculation model, extract the calculation results and the monitoring values of the monitoring cross - section, perform time - history curve comparison analysis, and display and output. Description of the Drawings
[0039] Figure 1 It is the schematic flow chart of the method for comparing the safety state of the dam structure provided by this application;
[0040] Figure 2Schematic diagram of the system module for comparative analysis of the safety behavior of the dam structure provided by this application;
[0041] Figure 3 VTK rendering flow chart used in the embodiment of this application;
[0042] Figure 4 Flow chart for reading the finite element calculation model and extracting profiles under VTK used in the embodiment of this application;
[0043] Figure 5 Flow chart for reading the 3D model of safety monitoring instruments under VTK used in the embodiment of this application;
[0044] Figure 6 Graphical schematic diagram for comparing before and after cutting the dam finite element calculation model under VTK in the embodiment of this application; where a) is before cutting the dam finite element calculation model; b) is after cutting the dam finite element calculation model;
[0045] Figure 7 Model and curve graph obtained after the system provided in the embodiment of this application processes the data of a certain measuring point, where a) is the comparison model graph; b) is the comparison curve graph. Detailed implementation manners
[0046] To make the purposes, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in 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. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. The components of the embodiments of the present invention usually described and illustrated herein can be arranged and designed in various different configurations.
[0047] Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed present invention, but merely represents the selected embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts fall within the scope of protection of the present invention.
[0048] The technical means not detailed in this application and not used to solve the technical problems of this application are set according to the common general knowledge in the art, and various ways of setting common general knowledge can be realized.
[0049] See Figure 1 , this application provides a method for comparative analysis of the safety behavior of the dam structure, including the following steps:
[0050] Step S11: Construct using tetrahedral elements vtkTetra, pentahedral elements vtkHexahedron, and hexahedral elements vtkHexahedron according to the finite element calculation model element types, and use vtkSTLReader to read the three-dimensional model of the safety monitoring instrument;
[0051] Step S12: Cut the finite element calculation model using the vtkCutter shear function in VTK according to the safety monitoring section position, and then use the vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map. Obtain the measured data through interactive clicking operations on the three-dimensional model of the safety monitoring instrument on the section;
[0052] Step S13: According to the information of the instrument coordinate point P(x, y, z), search for the nearest calculation unit M. The types of M include tetrahedral elements, pentahedral elements, and hexahedral elements. Tetrahedral elements have 4 nodes, pentahedral elements have 6 nodes, and hexahedral elements have 8 nodes. The nodes of the unit M can be represented as P i (x i ,y i ,z i ). Use the inverse distance weighted interpolation algorithm to interpolate the calculated values of the unit nodes to obtain the finite element calculated value s at the instrument position;
[0053] The inverse distance weighted interpolation function is
[0054]
[0055] where, is the distance from the point P(x, y, z) to the point P i (x i ,y i ,z i ); v represents the scalar value of the finite element calculation node, which can be the displacement, settlement, etc. of the dam; i = 1, 2,..., n, and k is a constant greater than 0, that is, the weighted power exponent, which is determined according to the actual interpolation effect;
[0056] Step S14: Interactively click on the finite element calculation section cloud map to extract the finite element analysis values s i at the measuring point positions under different time points in a certain period of the dam and compare with the measured values r i , i = 1, 2,..., n, where n is the number of time points in this period.
[0057] Through this method, a large amount of historical measured data and finite element analysis data can be visualized, which is convenient for comparison to find differences, reduces the manual comparison workload, improves efficiency, reduces the possibility of error generation, and improves the accuracy of the comparison results.
[0058] This method uses the VTK library for processing, which can accurately obtain the model and improve the comparison accuracy.
[0059] Preferably, it further includes the following steps: Step S20: Using a curve fitting algorithm to perform curve fitting on multiple finite element analysis values s of the measured points obtained i and the measured values r i respectively, and then display the results.
[0060] By displaying after curve fitting, it is convenient for users to compare historical measured data and finite element data, and improve the visibility of the comparison of measured point data.
[0061] Preferably, step S20 includes the following steps: Step S21: Using the quadratic B-spline curve fitting algorithm to perform curve fitting on the finite element analysis values and the measured values at the measured point positions respectively, obtaining the finite element value curve S and the measured value curve R;
[0062] The quadratic B-spline curve function used is:
[0063]
[0064] Among them, A1 = x1 - x0, B1 = y1 - y0,
[0065] Step S22: Display and output the obtained finite element analysis value curve S and the measured value curve R on the same coordinate system graph.
[0066] Using the quadratic B-spline curve for fitting, the accuracy of the obtained curve fitting result is relatively high.
[0067] Preferably, after the quadratic B-spline curve fitting in step S21, every 3 adjacent points are used to draw a parabola according to formula (2), and the front and back two curves are tangent at the endpoints of each segment.
[0068] See Figure 2 , on the other hand, this application also provides a dam structure safety state comparison analysis system, including:
[0069] Instruction receiving module 100: Used to receive user operation instructions. The user selects the finite element calculation result file on the device touch display screen. The finite element calculation result file includes: model unit node information, displacement information, etc., and touches the confirmation button to import;
[0070] Finite element calculation data reading module 101: Used to detect whether there is a device event for importing model data triggered by the user. If it detects that there is a user-triggered import device event, it reads the model data path and saves the model data in the device memory;
[0071] 3D Model Construction Module 102: It is used to establish 3D visualization models of finite element calculation results and safety monitoring instruments respectively in the VTK visualization scene. According to the element types of the finite element calculation model, it is constructed using tetrahedron elements vtkTetra, pentahedron elements vtkHexahedron, and hexahedron elements vtkHexahedron respectively. The 3D model of the safety monitoring instrument is read using vtkSTLReader;
[0072] Monitoring Section and Measured Data Acquisition Module 103: It is used to cut the finite element calculation model according to the position of the safety monitoring section using the vtkClipPolyData cutting function in VTK, and then use the vtkCutter shearing function and vtkPlane plane function to fill the cut area to generate a cloud map of the finite element calculation section. The measured data is obtained by interactive clicking operations on the 3D model of the safety monitoring instrument on the section, and the data is stored in the device memory;
[0073] Calculated Value Search Module 104: It is used to search for the finite element calculated values at the positions of the safety monitoring points. According to the instrument coordinate information, the nearest calculation element is searched, and the interpolation algorithm is used to interpolate the calculated values of the element nodes to obtain the finite element calculated values at the instrument positions, and the data is stored in the device memory;
[0074] Single-Value Comparison and Analysis Module 105: It is used to obtain and compare the finite element calculated values and measured values of the monitoring points while performing interactive clicks on the finite element calculation section cloud map, and extract the finite element analysis values and measured values at the measuring point positions under the required comparison time conditions during the comparison period of the dam. The data is stored in the device memory.
[0075] This device can effectively realize the visual display of historical measured data and finite element data according to the measuring points and their positions, realize efficient and accurate comparison, without manual operation, and has high processing efficiency.
[0076] Preferably, the interpolation algorithm used in the calculated value search module 104 is the inverse distance weighted interpolation method.
[0077] Preferably, it further includes:
[0078] Curve Fitting Module 106: It is used to perform curve fitting on the finite element analysis values and measured values at the measuring point positions respectively using the curve fitting algorithm to obtain the fitting curves;
[0079] Curve Comparison and Plotting Module 107: It is used to display the fitting curves and the measured value curves.
[0080] By adopting curve fitting, the data of the measuring points can be visually displayed, and the comparison efficiency and accuracy can be improved.
[0081] Preferably, the curve fitting algorithm used in the curve fitting module 106 is a quadratic B-spline curve.
[0082] Embodiment
[0083] The dam structure safety performance comparative analysis system provided by this application includes
[0084] Instruction receiving module 100: used to receive user operation instructions. The user selects a finite element calculation result file on the device touch display screen. The finite element calculation result file includes model unit node information, displacement information, etc., and touches the confirmation button to import.
[0085] Finite element calculation data reading module 101: used to detect whether there is a model data import device event triggered by the user. If it detects that there is a user-triggered import device event, it reads the model data path and saves the model data in the device memory.
[0086] 3D model construction module 102: used to respectively establish 3D visualization models of finite element calculation results and safety monitoring instruments in the VTK visualization scene. According to the finite element calculation model unit type, it is respectively constructed using tetrahedron unit vtkTetra, pentahedron unit vtkHexahedron, and hexahedron unit vtkHexahedron. The 3D model of the safety monitoring instrument is read using vtkSTLReader.
[0087] Monitoring section and measured data acquisition module 103: used to cut the finite element calculation model according to the safety monitoring section position using the vtkClipPolyData cutting function in VTK, and then use the vtkCutter shearing function and vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map. Through interactive click operations on the 3D model of the safety monitoring instrument on the section, the measured data is obtained and stored in the device memory.
[0088] Calculated value search module 104: used to search for the finite element calculated value at the safety monitoring point position. According to the instrument coordinate information, the nearest calculation unit is searched, and the interpolation algorithm is used to interpolate the calculated values of the unit nodes to obtain the finite element calculated value at the instrument position, which is stored in the device memory.
[0089] Single value comparative analysis module 105: used to compare the finite element calculated value and the measured value of the monitoring point while performing an interactive click on the finite element calculation section cloud map, and extract the finite element analysis value and the measured value at the measuring point position under the required comparison time conditions within the comparison period of the dam, which are stored in the device memory.
[0090] Curve fitting module 106: It is used to perform curve fitting on the finite element analysis values and measured values at the measuring point positions respectively by using a curve fitting algorithm to obtain fitting curves;
[0091] Curve comparison and plotting module 107: It is used to display the fitting curve and the measured value curve.
[0092] Preferably, the interpolation algorithm used in the calculated value search module 104 is the inverse distance weighted interpolation method.
[0093] Preferably, the curve fitting algorithm used in the curve fitting module 106 is the quadratic B-spline curve.
[0094] The module referred to in the present invention refers to a series of computer programs that can be executed by the processor 30 and can complete fixed functions, and are displayed in the display device 20. In this embodiment, for the specific functions of each module in the system, refer to the method flow chart.
[0095] Another aspect of the present application also provides a method for comparative analysis of the safety state of a dam structure, including the following steps:
[0096] Step S10: On the finite element calculation data reading module 101, after reading the finite element calculation analysis results of the dam under a certain reservoir water level condition, perform object rendering on the read results in the VTK environment;
[0097] In this embodiment, the object rendering of the read finite element calculation analysis results in the VTK environment is performed according to the existing method. For example, reference can be made to Figure 3 the VTK rendering flow chart in. Specifically refer to "The Visualization Toolkit" http: / / www.vtk.org.
[0098] Step S11: On the three-dimensional model construction module 102, construct respectively by using the tetrahedral element vtkTetra, pentahedral element vtkHexahedron, and hexahedral element vtkHexahedron according to the finite element calculation model element types, and use vtkSTLReader to read the three-dimensional model of the safety monitoring instrument; The method for reading the finite element calculation model in the VTK environment is performed according to the existing method. For details, refer to Figure 4 the flow chart for reading the finite element calculation model and extracting cross-sections in the VTK environment. Specifically refer to "TheVisualization Toolkit" http: / / www.vtk.org. The method for reading the three-dimensional model of the safety monitoring instrument is performed according to the existing method. For example, refer to Figure 5 the flow chart for reading the three-dimensional model of the safety monitoring instrument in the VTK environment; Specifically refer to "TheVisualization Toolkit" http: / / www.vtk.org.
[0099] Step S12: On the monitoring section and measured data acquisition module 103, according to the positions of the safety monitoring sections, use the vtkCutter shear function in VTK to cut the finite element calculation model, and then use the vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map. In this embodiment, the models before and after the cutting operation are as shown in Figure 6 a) and b). Interactive click operations on the three-dimensional model of the safety monitoring instrument on the section are used to obtain the measured data. The process of extracting the finite element calculation model profile under the VTK environment is as shown in Figure 5 .
[0100] Step S13: On the calculated value search module 104, according to the coordinate point P(x, y, z) information of the instrument in the model, search for the nearest calculation unit M. The types of M include tetrahedral elements, pentahedral elements, and hexahedral elements. Tetrahedral elements have 4 nodes, pentahedral elements have 6 nodes, and hexahedral elements have 8 nodes. The nodes of the unit M can be expressed as P i (x i , y i , z i ). The inverse distance weighted interpolation algorithm is used to interpolate the calculated values of the unit nodes. The inverse distance weighted interpolation function is as follows, and the finite element calculated value s at the instrument position is obtained;
[0101] The inverse distance weighted interpolation function is
[0102]
[0103] where is the distance from the point P(x, y, z) to the point P i (x i , y i , z i ); v represents the scalar value of the finite element calculation node, which can be the displacement, settlement, etc. of the dam; i = 1, 2,..., n, and k is a constant greater than 0, that is, the weighted power exponent, which is determined according to the actual interpolation effect.
[0104] Step S14: On the single value comparison and analysis module 105, while interactively clicking on the finite element calculation section cloud map, obtain the finite element calculation value and the measured value of the monitoring point for comparison. Using the above method, the finite element analysis value s of the measuring point position under different time points in a certain period of the dam can be extracted i and the measured value r i , i = 1, 2,..., n, where n is the number of time points in this period.
[0105] Step S21: On the curve fitting module 106, use the quadratic B-spline curve fitting algorithm to perform curve fitting on the finite element analysis values and measured values at the measurement point positions respectively. The finite element value curve is S, and the measured value curve is R.
[0106] The quadratic B-spline curve function is:
[0107]
[0108] Where, A1 = x1 - x0, B1 = y1 - y0,
[0109] When performing quadratic B-spline curve fitting, every adjacent 3 points are used to draw a parabola according to function (2), and the front and back two curves are tangent at the endpoints of each segment.
[0110] Step S22: On the curve comparison and plotting module 107, display and output the finite element analysis value curve S and the measured value curve R on the same coordinate system graph.
[0111] Using the method provided in this application, taking the deformation as shown in Figure 7 a as an example, take the monitoring data of a certain measurement point in a period. The specific data is shown in Table 1. Extract the calculated value of the finite element analysis at this monitoring point at the corresponding time point, and draw a curve on the same graph as shown in Figure 7 b for comparison
[0112] Table 1 Monitoring data and finite element calculation data of the surface deformation observation point TPDB4 - 8
[0113]
[0114] From Figure 7 it can be seen that the model obtained by using the method provided in this application is convenient for visually and visually obtaining the comparison results, with relatively high comparison efficiency and accuracy. The method provided in this application integrates the finite element analysis and the monitoring data, and the obtained results are more intuitive and convenient for comparison, effectively assisting engineering decision-making.
[0115] Although the present invention has been described in detail with reference to the foregoing embodiments, for those skilled in the art, they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements for some of the technical features. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.
Claims
1. A method for comparative analysis of the safety state of a dam structure, characterized in that, It includes the following steps: Step S11: Construct using tetrahedral element vtkTetra, pentahedral element vtkHexahedron, and hexahedral element vtkHexahedron respectively according to the finite element calculation model element type, and read the three-dimensional model of the safety monitoring instrument using vtkSTLReader; Step S12: Cut the finite element calculation model using the vtkCutter shear function in VTK according to the safety monitoring section position, then use the vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map, and obtain the measured data through interactive clicking operations on the three-dimensional model of the safety monitoring instrument on the section; Step S13: According to the information of the instrument coordinate point P(x, y, z), search for the nearest calculation unit M. The types of M include tetrahedral elements, pentahedral elements, and hexahedral elements. A tetrahedral element has 4 nodes, a pentahedral element has 6 nodes, and a hexahedral element has 8 nodes. The nodes of the unit M can be expressed as P i (x i , y i , z i ). Use the inverse distance weighted interpolation algorithm to perform interpolation calculations on the calculated values of the unit nodes to obtain the finite element calculated value s of the instrument position; The inverse distance weighted interpolation function is Among them, v i represents the finite element analysis value of the i-th node of element M, the displacement and settlement analysis values of the dam, where i = 1, 2, …, n; represents the power value of the distance from the i-th node of element M to the instrument coordinate point P. k is a constant greater than 0, and its value is determined according to the actual interpolation effect. When k = 1, is the distance from point P(x, y, z) to point P i (x i , y i , z i ); its value is determined according to the actual interpolation effect; Step S14: Interactively click and extract the finite element analysis value s of the measuring point position of the dam under different time points within a certain period in the finite element calculation section cloud diagram i and the measured value r i , where i = 1, 2, …, n, and n is the number of time points within this period, for comparison.
2. The method for comparative analysis of the safety behavior of the dam structure according to claim 1, wherein It also includes the following steps: Step S20: Use a curve fitting algorithm to fit and display multiple finite element analysis values s of the obtained measurement points i and the measured value r i after fitting respectively.
3. The method for comparative analysis of the safety state of a dam structure according to claim 2, wherein Step S20 includes the following steps: Step S21: Use the quadratic B-spline curve fitting algorithm to perform curve fitting on the finite element analysis values and measured values at the measuring point positions respectively to obtain the finite element value curve S and the measured value curve R; The quadratic B-spline curve function used is: wherein, A1 = x1 - x0, B1 = y1 - y0, Step S22: Display and output the obtained finite element analysis value curve S and measured value curve R on the same coordinate system graph.
4. The method for comparative analysis of the safety state of a dam structure according to claim 3, characterized in that, In Step S21, parabolas are drawn for every 3 adjacent points after quadratic B-spline curve fitting according to Equation (2), and the front and back curves are tangent at the endpoints of each segment.
5. A system for comparative analysis of the safety behavior of a dam structure, characterized in that, It includes: Instruction receiving module (100): Used to receive user operation instructions. The user selects the finite element calculation result file on the device touch display screen. The finite element calculation result file includes model element node information, displacement information, etc., and touches the OK button to import; Finite element calculation data reading module (101): Used to detect whether there is a device event for importing model data triggered by the user. If it detects that there is a user-triggered import device event, read the model data path and save the model data in the device memory; Three-dimensional model construction module (102): Used to establish the finite element calculation result and the three-dimensional visualization model of the safety monitoring instrument in the VTK visualization scene respectively. Construct using tetrahedral element vtkTetra, pentahedral element vtkHexahedron, and hexahedral element vtkHexahedron respectively according to the finite element calculation model element type, and read the three-dimensional model of the safety monitoring instrument using vtkSTLReader; Monitoring section and measured data acquisition module (103): Used to cut the finite element calculation model using the vtkClipPolyData cutting function in VTK according to the safety monitoring section position, then use the vtkCutter shear function and vtkPlane plane function to fill the cut area to generate a finite element calculation section cloud map, and obtain the measured data through interactive clicking operations on the three-dimensional model of the safety monitoring instrument on the section, and store the data in the device memory; Calculated value search module (104): It is used to search for the finite element calculated values at the safety monitoring point positions. According to the instrument coordinate information, it searches for the nearest calculation unit, and uses the interpolation algorithm to perform interpolation calculation on the calculated values of the unit nodes to obtain the finite element calculated values at the instrument positions, and stores the data in the device memory; Single value comparison and analysis module (105): It is used to compare the finite element calculated values and the measured values of the monitoring points while interactively clicking on the finite element calculation section contour map, and extract the finite element analysis values and the measured values at the measuring point positions under the required comparison time conditions during the period to be compared of the dam, and store the data in the device memory.
6. The comparative analysis system for the safety behavior of a dam structure according to claim 5, characterized in that, The interpolation algorithm used in the calculated value search module (104) is the inverse distance weighted interpolation method.
7. The dam structure safety performance comparison and analysis system according to claim 5, characterized in that It also includes: Curve fitting module (106): It is used to perform curve fitting on the finite element analysis values and the measured values at the measuring point positions respectively by using the curve fitting algorithm to obtain the fitting curves; Curve comparison and plotting module (107): It is used to display the fitting curves and the measured value curves.
8. The dam structure safety performance comparative analysis system according to claim 7, characterized in that The curve fitting algorithm used in the curve fitting module (106) is the quadratic B-spline curve.
Citation Information
Patent Citations
Large-span bridge monitoring data spatial interpolation and visualization method
CN104899349A
Marine structure observed behavior multi-dimensional analysis method based on virtual reality and system
CN106339820A