Flood discharge and sand flushing service gate vibration and water flow excitation coupling system and method

Abstract

The invention relates to the technical field of gate engineering, in particular to a vibration and water flow excitation coupling system and method for a flood discharge and sand flushing service gate. The method comprises the following steps that a multi-cylinder hydraulic jacking structure drawing is obtained, and a multi-cylinder hydraulic jacking three-dimensional model is constructed; the multi-cylinder hydraulic jacking three-dimensional model is used for gate opening simulation, and gate opening parameters are obtained; the instantaneous discharge amount is calculated through gate opening parameters, and the vortex-induced vibration degree is analyzed; performing gate leaf vibration detection based on the vortex-induced vibration degree to obtain gate leaf vibration data; performing structural fatigue analysis according to the gate leaf vibration data of the gate to obtain gate leaf structural fatigue data; predicting the gate leaf fatigue life based on the gate leaf structure fatigue data; and calculating the instantaneous water pressure according to the gate opening parameter. The gate structure state real-time analysis is realized based on the gate engineering technology, and the flood discharge safety and the operation reliability are improved.

Description

Technical Field

[0001] This invention relates to the field of gate engineering technology, and in particular to a coupling system and method for vibration and water flow excitation of a flood discharge and sand flushing gate. Background Technology

[0002] Currently, flood discharge and sediment flushing gates are widely used in hydraulic engineering projects, ship locks, and large-scale hydraulic machinery. However, under high flow, frequent opening and closing, and complex hydraulic conditions, the gate leaves are susceptible to vortex-induced vibration, fluid impact, and instantaneous water pressure, leading to increased vibration amplitude, intensified structural fatigue, and localized material damage, thus affecting the gate's operational safety and service life. Existing technologies primarily focus on increasing the overall thickness of the gate plates or using high-strength steel to improve rigidity. However, this approach not only increases the weight and manufacturing cost of the gate leaves but also places higher demands on the lifting drive system. Furthermore, traditional vibration analysis and fatigue assessment often rely on empirical formulas or localized testing, making it difficult to comprehensively reflect the dynamic load distribution under different opening and flow conditions. This results in inaccurate prediction of localized fatigue risks, easily leading to fatigue cracks and material failure. In addition, existing gate lifting schemes are mostly single-cylinder or double-cylinder lifting, which suffers from poor synchronization, insufficient consideration of vibration superposition effects, and inadequate utilization of localized reinforcement structures under multi-cylinder synchronization and complex hydraulic excitation, failing to effectively guarantee long-term reliable operation. Summary of the Invention

[0003] Based on this, it is necessary for the present invention to provide a coupling system and method for vibration and water flow excitation of flood discharge and sand flushing gate, so as to solve at least one of the above-mentioned technical problems.

[0004] To achieve the above objectives, a coupling method for the vibration and water flow excitation of a flood discharge and sediment flushing gate includes the following steps: Step S1: Obtain the structural drawings of the multi-cylinder hydraulic jacking structure and construct a three-dimensional model of the multi-cylinder hydraulic jacking; use the three-dimensional model of the multi-cylinder hydraulic jacking to simulate the gate opening and obtain the gate opening parameters; Step S2: Calculate the instantaneous discharge flow using the gate opening parameters and analyze the degree of vortex-induced vibration; perform gate leaf vibration detection based on the degree of vortex-induced vibration to obtain gate leaf vibration data; Step S3: Perform structural fatigue analysis based on the gate leaf vibration data to obtain gate leaf structural fatigue data; predict the gate leaf fatigue life based on the gate leaf structural fatigue data; Step S4: Calculate the instantaneous water pressure based on the gate opening parameters; calculate the stress on the gate leaf plate using the instantaneous water pressure; correct the gate leaf fatigue life based on the stress on the gate leaf plate to obtain the gate life data.

[0005] Preferably, this specification also provides a coupling system for vibration and water flow excitation of a flood discharge and sediment flushing gate, used to execute the coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate as described above. The coupling system for vibration and water flow excitation of the flood discharge and sediment flushing gate includes: The gate opening simulation module is used to acquire drawings of the multi-cylinder hydraulic jacking structure and construct a three-dimensional model of the multi-cylinder hydraulic jacking; the gate opening is simulated using the three-dimensional model of the multi-cylinder hydraulic jacking to obtain the gate opening parameters. Gate leaf vibration detection is used to calculate instantaneous discharge flow using gate opening parameters and analyze the degree of vortex-induced vibration; based on the degree of vortex-induced vibration, gate leaf vibration detection is performed to obtain gate leaf vibration data; The gate leaf fatigue life prediction is used to perform structural fatigue analysis based on gate leaf vibration data to obtain gate leaf structural fatigue data; and to predict the gate leaf fatigue life based on the gate leaf structural fatigue data. The instantaneous water pressure calculation module is used to calculate the instantaneous water pressure based on the gate opening parameters; calculate the stress on the gate leaf plate using the instantaneous water pressure; and correct the fatigue life of the gate leaf based on the stress on the gate leaf plate to obtain the gate life data.

[0006] The present invention has the following beneficial effects: Firstly, it enables accurate understanding of the gate opening variation pattern during the design and simulation stages, providing reliable data input for discharge flow and vibration analysis, thereby improving the accuracy and reliability of the design.

[0007] Secondly, it can quantify the impact of water flow excitation on the vortex-induced vibration of the door leaf, obtain door leaf vibration characteristic data, provide a scientific basis for structural health monitoring and anomaly early warning, thereby improving operational safety.

[0008] Thirdly, it can assess the fatigue accumulation and structural damage of the door leaf during long-term operation, predict the fatigue life of the door leaf, provide a reference for scientifically scheduling maintenance cycles, reduce the risk of fatigue failure, and improve reliability.

[0009] Fourthly, by analyzing the stress state of the gate leaf plate, the fatigue life assessment can be corrected, making the life prediction closer to the actual stress situation, thereby extending the service life of the gate and reducing operation and maintenance costs. Attached Figure Description

[0010] Other features, objects, and advantages of the invention will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings: Figure 1 This is a schematic flowchart illustrating the steps of a coupling method for vibration and water flow excitation of a flood discharge and sand flushing gate according to the present invention. Figure 2 This is a detailed flowchart of step S1 in the present invention; Figure 3 This is a detailed flowchart of step S2 in the present invention; The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0011] The technical method of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. Based on the embodiments of this invention, all other embodiments obtained by those skilled in the art without inventive effort are within the scope of protection of this invention.

[0012] Furthermore, the accompanying drawings are merely illustrative of the invention and are not necessarily drawn to scale. The same reference numerals in the drawings denote the same or similar parts, and therefore repeated descriptions of them will be omitted. Some block diagrams shown in the drawings are functional entities and do not necessarily correspond to physically or logically independent entities. These functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor methods and / or microcontroller methods.

[0013] It should be understood that although the terms "first," "second," etc., may be used herein to describe various units, these units should not be limited by these terms. These terms are used merely to distinguish one unit from another. For example, without departing from the scope of the exemplary embodiments, a first unit may be referred to as a second unit, and similarly, a second unit may be referred to as a first unit. The term "and / or" as used herein includes any and all combinations of one or more of the associated listed items.

[0014] To achieve the above objectives, please refer to Figures 1 to 3 This invention provides a coupling method for the vibration and water flow excitation of a flood discharge and sediment flushing gate, the method comprising the following steps: Step S1: Obtain the structural drawings of the multi-cylinder hydraulic jacking structure and construct a three-dimensional model of the multi-cylinder hydraulic jacking; use the three-dimensional model of the multi-cylinder hydraulic jacking to simulate the gate opening and obtain the gate opening parameters; In this embodiment, information such as the hydraulic cylinder layout, cylinder diameter, piston rod length, cylinder end fixed support position, and overall support frame dimensions are obtained from the multi-cylinder hydraulic lifting structure drawings. Specifically, the starting and ending coordinates of the axis of each hydraulic cylinder in the two-dimensional drawings are recorded and converted into three-dimensional spatial coordinates, with accuracy controlled to the millimeter level. Based on the cylinder diameter and piston rod diameter indicated in the drawings, a solid model of each hydraulic cylinder is constructed in three-dimensional space. The wall thickness and piston rod diameter are directly assigned according to the drawings, for example, a cylinder wall thickness of three centimeters and a piston rod diameter of fifteen centimeters. The hydraulic cylinder model is assembled onto the gate support frame, ensuring that the hydraulic cylinder mounting holes are perfectly aligned with the gate hinge points, with an assembly error not exceeding one millimeter. Then, based on the relationship between the hydraulic cylinder piston extension and contraction and the gate rotation, the gate's change from fully closed to fully open is simulated, with each step incrementing by half a degree. The gate leaf angle, gate leaf edge height, and the relative position of the gate leaf to the water surface are recorded at each moment, forming a complete dataset of gate opening parameters. The dimensions, thickness, and hinge point positions of the gate leaf must be strictly entered according to the construction drawings. For example, the gate height is five meters, the gate width is three meters, the thickness is twelve centimeters, and the hinge point is thirty centimeters from the bottom.

[0015] Step S2: Calculate the instantaneous discharge flow using the gate opening parameters and analyze the degree of vortex-induced vibration; perform gate leaf vibration detection based on the degree of vortex-induced vibration to obtain gate leaf vibration data; In this embodiment, the effective flow area of ​​the gate is determined. This area is the opening area from the bottom of the gate leaf to the water surface, calculated based on the gate leaf height and the upstream water head height collected by the water level sensor. The water level is collected every second, with a height range from ten to fifty meters. The flow velocity through the gate gap is determined based on the opening parameters and water level height. Instantaneous discharge data is obtained by analyzing the flow velocity and opening area. Then, based on the changes in water velocity and depth, the velocity gradient information of the water flow is extracted. The velocity gradient is used to analyze the water flow rotation and the vortex regions formed, while recording the core region position and rotation intensity of each vortex. By analyzing the size of the vortex core region, rotation intensity, and the location of the vortex, the point of action of the vortex is correlated with the gate leaf surface, and the water flow excitation force at each point of action is measured along the gate leaf normal. The instantaneous excitation forces at all points of action are spatially superimposed to obtain the vibration load distribution at each position of the gate leaf, ultimately forming complete gate leaf vibration data.

[0016] Step S3: Perform structural fatigue analysis based on the gate leaf vibration data to obtain gate leaf structural fatigue data; predict the gate leaf fatigue life based on the gate leaf structural fatigue data; In this embodiment, vibration amplitude and period information are extracted to assess the stiffness of the door leaf plate. The plate stiffness is determined by the material type and thickness; the material is steel plate with a thickness of 12 cm. Local stress analysis is performed on the door leaf plate to identify instability areas. If the stress caused by local vibration exceeds the local bearing capacity of the steel plate, the material thinning value in that area is recorded. The thinning range is typically between 0 and 3 mm. The degree of local collapse is assessed based on the material thinning data, generating plate damage data. Deformation analysis is performed on the door leaf frame. The frame displacement is extracted from the door leaf vibration signal, and the maximum deflection of the frame is obtained through vibration signal processing. The deformation area of ​​the frame is then defined based on the deflection. The plate damage data and frame deformation data are integrated to form door leaf structural fatigue data. Based on this, using steel fatigue performance data and cyclic loading information, the number of cycles the door leaf can withstand is calculated according to the cumulative damage principle, thereby predicting the fatigue life of the door leaf.

[0017] Step S4: Calculate the instantaneous water pressure based on the gate opening parameters; calculate the stress on the gate leaf plate using the instantaneous water pressure; correct the gate leaf fatigue life based on the stress on the gate leaf plate to obtain the gate life data.

[0018] In this embodiment, water pressure is measured along the gate leaf height according to water depth distribution, with water pressure values ​​sampled every five centimeters. Upstream water level is collected via sensors. The water pressure is then converted into stress distribution on the gate leaf surface. Surface stress is evaluated by measuring water pressure, the surface area subjected to force, and the plate thickness. The surface stress at each sampling point is recorded individually. Subsequently, the surface stress data is combined with the gate leaf structural fatigue data obtained in step S3, and local damage and remaining lifespan are corrected according to the principle of cumulative fatigue. After fatigue life correction at each sampling point, a gate leaf lifespan dataset is generated. Continuous calculations are performed on the entire gate leaf from bottom to top to ensure a complete opening and closing cycle of the gate. The data update cycle is once every 0.1 seconds to obtain continuous lifespan change information.

[0019] Preferably, step S1 specifically includes: Step S11: Obtain the drawings of the multi-cylinder hydraulic jacking structure; In this embodiment, complete design drawings of the multi-cylinder hydraulic jacking structure are obtained from the construction or design unit. These drawings include a hydraulic cylinder layout plan, longitudinal section, gate support frame structure, hinge point locations, and hydraulic system piping diagrams. The scale of the drawings must be clearly defined, requiring a scale of 1:50 or 1:100 to ensure measurement accuracy at the millimeter level. The hydraulic cylinder dimensions recorded in the drawings include cylinder diameter, piston rod diameter, cylinder length, and cylinder wall thickness. All parameters must be consistent with the design documents; for example, a hydraulic cylinder diameter of 200 mm, a piston rod diameter of 150 mm, a cylinder wall thickness of 30 mm, and a maximum piston stroke of 1.5 meters. By reading the coordinates marked on the drawings, the center points of the mounting holes of each hydraulic cylinder, the start and end points of the cylinder shaft are recorded in a two-dimensional coordinate table, with the coordinate unit being millimeters. The three-dimensional dimensions of the gate leaf, the hinge point locations, and the thickness and width of the gate leaf also need to be obtained from the drawings; for example, a gate height of 5 meters, a gate width of 3 meters, a thickness of 12 centimeters, and a hinge point distance of 30 centimeters from the bottom. The data acquisition process uses digital measurement tools, such as electronic rulers or coordinate measuring machines, to convert the annotation points on the two-dimensional drawings into precise values, providing basic data for subsequent three-dimensional modeling.

[0020] Step S12: Use the multi-cylinder hydraulic jacking structure drawing to identify the two-dimensional contour and determine the coordinates of the two-dimensional contour; In this embodiment, after acquiring the drawing data, a two-dimensional drawing tool is used to depict the outlines of the main components of the hydraulic jacking device as two-dimensional planar graphics, including the outer outline of the hydraulic cylinder, the outline of the piston rod, the outline of the support frame, and the position of the gate hinge point. First, the boundary lines of each component in the drawing are vectorized. The outline dimensions are obtained by measuring the length of the line segments between the marked points on the drawing; for example, the hydraulic cylinder is 1.8 meters long and 0.2 meters wide, and the piston rod is 1.2 meters long and 0.15 meters in diameter. Then, the coordinates of the endpoints of each boundary line are recorded as two-dimensional coordinates in millimeters. The origin is selected at the lower left corner of the gate mounting foundation to ensure that the two-dimensional outline of the entire hydraulic jacking device is in the same coordinate system. For the beam-column nodes of the support frame, the two-dimensional coordinates of the nodes are recorded, including the X and Y directions, with an accuracy controlled within ±1 millimeter. For the hinge point and the center point of the mounting hole, the two-dimensional coordinates of the center point are calculated using the radius or hole diameter marked on the drawing, and the hole diameter data is recorded; for example, the hole diameter is 40 millimeters. Once the two-dimensional contour coordinate table is formed, it contains the coordinate information of all hydraulic cylinders, pistons, support frame nodes, and gate hinge points, providing a complete reference for subsequent three-dimensional patch construction.

[0021] Step S13: Construct a three-dimensional surface patch for multi-cylinder hydraulic jacking based on two-dimensional contour coordinates; In this embodiment, each component is stretched vertically using two-dimensional contour coordinates to generate a three-dimensional patch. The hydraulic cylinder is stretched vertically along its two-dimensional circular contour to form a three-dimensional patch. The cylinder length is based on the design drawings, for example, 1.8 meters, with a wall thickness of 30 millimeters. During the stretching process, the patch thickness is ensured to match the cylinder wall thickness. The piston rod is stretched vertically along its two-dimensional circular contour, resulting in a length of 1.2 meters and a diameter of 0.15 meters. The beams and columns of the support frame are stretched vertically along two-dimensional rectangular contours, resulting in beams 0.4 meters high and 0.3 meters wide, and columns 2 meters high with a cross-sectional width and thickness of 50 millimeters. The patch node accuracy is controlled within 1 millimeter. The gate leaf is stretched vertically along its two-dimensional rectangular contour to generate a three-dimensional patch with a thickness of 12 centimeters. During the construction of the three-dimensional patches, the three-dimensional vertex coordinates, boundary points, and normal information of each component are recorded completely to ensure the accuracy of the component's shape and size in three-dimensional space. The gaps between the patches are consistent with the design gaps; for example, the gap between the hydraulic cylinder and the support frame is 5 millimeters. All the patch data form a 3D component library, providing a complete geometric basis for subsequent assembly.

[0022] Step S14: Assemble the components according to the three-dimensional surface of the multi-cylinder hydraulic jacking to obtain the three-dimensional model of the multi-cylinder hydraulic jacking. In this embodiment, the hydraulic cylinder faceplate is positioned in three-dimensional space according to two-dimensional contour coordinates. The cylinder mounting hole is aligned with the support frame hole, and the installation error must not exceed 1 mm. The piston rod is inserted into the hydraulic cylinder, and its end is connected to the gate hinge point, ensuring that the extension and retraction direction is consistent with the rotation axis of the gate leaf. The piston rod has a total stroke of 1.2 meters, and the relative position of a node is recorded every 10 centimeters. The support frame beams and columns are assembled according to the design coordinates, with the node positions precisely matching the beam and column faceplates. The connection surfaces of the beams and columns fit together, with no gap deviation greater than 2 mm. The gate leaf is connected to the support frame through the hinge point, with the center of the hinge point aligned with the center of the support frame hole, and the hinge gap controlled to 3 mm. During the assembly process, the relative position and angle of each component are recorded. Each hinge point angle is divided into 10 sampling points from 0° to 60°, forming a complete three-dimensional assembly coordinate table to provide basic data for opening simulation.

[0023] Step S15: Use a multi-cylinder hydraulic lifting three-dimensional model to simulate the gate opening and obtain the gate opening parameters.

[0024] In this embodiment, the extension length of each hydraulic cylinder piston rod is evenly divided into 30 sampling points within the range of 0 to 1.5 meters, and the gate leaf rotation angle corresponding to the extension length is recorded. The gate leaf rotation angle is calculated through geometric relationships, and the gate leaf edge height and the distance between the gate leaf and the water surface are recorded every 0.5 degrees. For example, when fully closed, the distance between the gate leaf and the water surface is 0 meters, and when fully open, the upper edge of the gate leaf is 2.5 meters from the water surface. During the simulation, the extension and retraction of each hydraulic cylinder, the position of the hinge point, and the constraint conditions of the support frame are combined and calculated to form a complete three-dimensional attitude dataset of the gate leaf for each sampling point. All sampling points form a continuous sequence of opening parameters, including the gate leaf rotation angle, edge height, and relative height to the water surface. The data is updated every 0.1 seconds to ensure that the simulation results cover the dynamic changes of the gate from fully closed to fully open. The data is saved as a three-dimensional coordinate table, containing the vertex coordinates and angle information of all key components, providing basic data for subsequent discharge flow calculation and vibration analysis.

[0025] Preferably, step S2 specifically includes: Step S21: Extract the effective flow area using the gate opening parameters; In this embodiment, the gate leaf rotation angle, the relative height of the lower edge of the gate leaf to the water surface, and the position of the gate leaf edge are acquired. Based on the gate leaf thickness and width, the opening area formed by the gate leaf and the water surface is divided into several horizontally equal segments, each segment being 5 centimeters high, with each segment corresponding to a rectangular opening unit. The width of each unit is determined by the gate leaf width and rotation angle, and the thickness is consistent with the gate leaf thickness. The area of ​​each rectangular unit is calculated, and the total effective flow area is obtained by summing the areas of all rectangular units, in square meters. For example, in the fully open state, the lower edge of the gate leaf is 2.5 meters from the water surface; the gate opening is divided into 50 rectangular units, each with an area of ​​0.15 square meters, for a total effective flow area of ​​7.5 square meters. This operation uses numerical integration to accumulate the area of ​​each segment, ensuring that the accuracy of the height and width data for each segment is controlled within 1 millimeter. The effective flow area at each time sampling point is recorded at a time interval of 0.1 seconds, forming a continuous effective flow area curve, providing a basis for subsequent flow velocity calculations.

[0026] Step S22: Obtain the upstream head height; calculate the gate velocity data based on the upstream head height and the effective flow area; In this embodiment, the upstream water head height is obtained by a water level sensor deployed upstream of the gate, collecting data once per second. The water head range is between 10 meters and 50 meters, and the sensor accuracy is ±5 millimeters. The collected water head data is correlated with the effective flow area obtained in step S21 to form time series data. Based on the opening area and water head height, the instantaneous velocity of the water flowing through the gate gap is calculated. The water flow velocity is segmented vertically, with each segment consisting of 5 centimeters, corresponding to the ratio of the water head difference to the effective flow area in each segment. A water head driving force is applied to each segment, considering the water flow friction coefficient and local resistance coefficient. The friction coefficient is set to 0.02, and the local resistance coefficient is determined to be 0.05 based on the shape of the gate leaf edge and hinge point. The water flow velocity data of each segment is integrated to obtain the instantaneous flow velocity distribution of the entire gate gap, in meters per second. The sampling frequency of the flow velocity data is consistent with the water head sampling frequency, forming a continuous time series for subsequent discharge calculation and vortex analysis.

[0027] Step S23: Perform flow rate integration calculation based on gate flow velocity data to obtain instantaneous discharge flow rate; In this embodiment, based on the flow velocity distribution, the gate's flow cross-section is divided into several small units, each corresponding to a rectangular area 5 cm high and equally divided in the width direction of the gate leaf. For each rectangular unit, the instantaneous flow rate of that unit is obtained by multiplying the water flow velocity of that unit by the unit area, in cubic meters per second. The instantaneous flow rates of all units are accumulated to obtain the instantaneous discharge rate of the entire gate leaf at that time point. For example, when the opening is 45°, the instantaneous flow rates of 50 units are calculated separately, and the accumulated instantaneous discharge rate is 8.3 cubic meters per second. The time step is 0.1 seconds, forming a complete instantaneous discharge rate sequence. Throughout the calculation process, the unit area error is controlled within 0.5 square centimeters, the flow rate data is recorded in cubic meters per second, and the time series corresponds completely with the opening parameter sequence, ensuring that the discharge rate calculation is synchronized with the gate leaf opening.

[0028] Step S24: Analyze the degree of vortex-induced vibration using instantaneous discharge flow rate; In this embodiment, instantaneous discharge data is used as the water flow excitation input. A vortex analysis grid is divided along the leaf blade surface, with each grid cell measuring 0.1 meters along the width and 0.05 meters along the height. For each grid cell, local velocity information is extracted based on the flow rate and velocity gradient, and the location and rotation intensity of the local vortex generation region are calculated. The vortex core region is defined as the area where the velocity gradient exceeds a threshold of 0.8 meters per second per meter. The core region area is obtained by measuring the coordinates of the vertices of the rotating region, and the rotation intensity is obtained by integrating the local flow velocity along the closed path within the core region. Based on the rotation intensity and the vortex core area, the direction and magnitude of the excitation force at each vortex point are determined. The direction is along the leaf blade normal, and the magnitude is calculated proportionally to the vortex intensity, in kilonewtons. The excitation forces at all vortex points are spatially superimposed to obtain the vortex-induced vibration degree of the entire leaf blade under the instantaneous discharge.

[0029] Step S25: Detect the gate leaf vibration based on the degree of vortex-induced vibration to obtain the gate leaf vibration data.

[0030] In this embodiment, vibration measurement nodes are installed along the surface of the door leaf according to the vortex excitation force distribution, with each node spaced 0.2 meters apart. Each node acquires normal vibration displacement and vibration velocity data, with displacement ranging from 0 to 5 millimeters and vibration velocity ranging from 0 to 0.5 meters per second. The node vibration signals are recorded using accelerometers at a sampling frequency of 100 Hz. The node vibration data is integrated according to the door leaf surface and frame structure partitions, recording the maximum vibration amplitude, vibration frequency, and vibration direction of each partition. A complete door leaf vibration data table is formed by mapping the node data to the vortex excitation force, including node location coordinates, vibration amplitude, vibration frequency, and vibration direction. The door leaf vibration data covers the entire height and width of the door leaf, with each sampling point synchronously recording a timestamp. Data accuracy is controlled to within 0.1 millimeters, forming continuous dynamic vibration monitoring data, providing an input basis for structural fatigue analysis.

[0031] Preferably, step S24 specifically includes: Step S241: Calculate the local Reynolds number using the instantaneous discharge flow rate; In this embodiment, the flow passage section of the gate is divided vertically into several equal-height units, each unit being 5 cm high, with the total height covering the flow area from the lower edge of the gate leaf to the water surface. Within each unit, the water flow velocity is measured in meters per second (m / s), obtained by dividing the flow rate by the unit area. Water flow characteristic parameters include water density (1000 kg / m³) and dynamic viscosity. The local Reynolds number is calculated based on a 0.1-second time interval, with the characteristic length of the flow cross-section taken as a unit height of 5 cm. The local Reynolds number is formed by combining local flow velocity, characteristic length, water density, and dynamic viscosity. Reynolds number is calculated for each unit, creating a complete spatial and temporal distribution matrix. For example, at an instantaneous discharge of 8 m / s, the unit flow velocity is 1.5 m / s, corresponding to a characteristic length of 5 cm, resulting in a local Reynolds number of 7500. The Reynolds number at each calculation point is recorded with integer precision to ensure accurate differentiation between laminar and turbulent regions in subsequent vortex analysis. After recording, a two-dimensional local Reynolds number matrix is ​​generated, with a sampling point every 0.1 meters horizontally along the door leaf width and every 0.05 meters vertically along the height, providing the basic input for vortex analysis.

[0032] Step S242: Perform vortex analysis based on the local Reynolds number to obtain vortex data; In this embodiment, the flow passage is divided into vortex analysis grids with a grid size of 0.1 meters along the width direction and 0.05 meters along the height direction. Each grid cell corresponds to a Reynolds number. Grid cells with a Reynolds number greater than 2000 are defined as turbulent regions, and local rotational motion is calculated by extracting velocity gradients. The vorticity field of each cell is measured using velocity gradients, with the direction of the vorticity field rotating perpendicular to the main flow direction. The vortex core region is determined by a continuous grid region with a local velocity gradient exceeding 0.8 m / s. The core area is obtained by statistically analyzing the areas of continuous grid cells, in square meters. Velocity integration is performed along a closed path in the vortex core region to obtain core rotation intensity data, in meters per second. The rotation center coordinates, core area, and rotation intensity of each core region are recorded to form a complete vortex data table, containing the positions and intensities of all core vortices. The time sampling interval is 0.1 seconds to ensure that the vortex data is synchronized with the instantaneous discharge flow, providing an input basis for vortex frequency calculation.

[0033] Step S243: Calculate the vortex frequency based on the vortex data, and evaluate the door leaf excitation force based on the vortex frequency; In this embodiment, the rotational intensity and core area of ​​the vortex core region are extracted. The vortex frequency is calculated by the ratio of the square root of the vortex rotational intensity to the core area, and the unit is Hertz. The vortex frequency is distributed across the entire leaf surface to form a frequency field, with each vortex frequency corresponding to its core region coordinates and direction of action. Based on the location of the vortex's point of action, the direction of the excitation force is determined along the leaf normal. The magnitude of the excitation force is proportional to the vortex rotational intensity and the core area, and the unit is kilonewtons. Within the entire height and width of the leaf, the excitation force is recorded at each vortex core point, with a lateral sampling interval of 0.1 meters, a longitudinal sampling interval of 0.05 meters, and a time step of 0.1 seconds. The excitation force at each vortex point of action is mapped to the leaf structure boundary to form a complete instantaneous excitation force matrix. The matrix records the direction, magnitude, and vortex frequency of each excitation force point on the leaf surface, providing input for subsequent vortex-induced vibration quantification.

[0034] Step S244: Optimize the degree of vortex-induced vibration based on the door leaf excitation force.

[0035] In this embodiment, vibration measuring points are arranged on the leaf surface at a node spacing of 0.2 meters. The normal displacement and vibration velocity are recorded at each node, with the displacement ranging from 0 to 5 millimeters and the vibration velocity ranging from 0 to 0.5 meters per second. The excitation force data of each node are superimposed to calculate the node vibration amplitude and frequency response, generating quantitative data of vortex-induced vibration across the entire leaf surface. Vibration amplitude and vibration velocity are maintained with an accuracy of 0.01 millimeters and 0.01 meters per second, respectively. Spatially weighted averages are performed on the node data from different regions of the leaf, weighted according to the relative distance between the node position and the vortex action point, to obtain the quantitative vortex-induced vibration degree of each region. The data recording time step is 0.1 seconds, forming a continuous dynamic sequence of vortex-induced vibration, providing fundamental data for leaf fatigue analysis and structural life correction.

[0036] Preferably, step S242 specifically includes: Extract velocity gradient information based on local Reynolds number; In this embodiment, based on the local Reynolds number two-dimensional matrix, uniform grids are established along the transverse and longitudinal directions of the gate flow section. The grid spacing is 0.1 meters along the width direction and 0.05 meters along the height direction, with the total grid covering the lower edge of the gate leaf to the full height of the water surface. For each grid cell, the velocity gradient vector is calculated using the Reynolds number and local flow velocity. The velocity gradient vector is sampled along both the horizontal and vertical directions with an accuracy of 0.01 meters per second per meter. The velocity gradient is calculated along both the main flow direction and the normal direction. The normal direction gradient is used to determine the intensity of rotational motion, and the main flow direction gradient is used to determine the shear rate. The data record for each grid cell includes the magnitude, direction, and corresponding coordinates of the velocity gradient. By continuously sampling at a time step of 0.1 seconds, the instantaneous velocity gradient matrix is ​​obtained, ensuring that the velocity change of each grid cell is completely described throughout the entire discharge process, providing the basic input for vorticity field calculation.

[0037] Calculate the vorticity field using velocity gradient information; In this embodiment, local rotation (vorticity) is calculated for each grid cell based on the instantaneous velocity gradient matrix. The vorticity vector direction is perpendicular to the main flow direction, and the unit is per second. The calculation method utilizes the rotation matrix formed by the velocity difference between adjacent grid cells and the grid spacing. The rotation must be measured in both the horizontal and vertical directions. The vorticity field covers the entire flow cross-section of the door leaf, and the magnitude, direction, and position of the vorticity for each grid cell are recorded. The calculation accuracy of vorticity is maintained at 0.01 seconds. The vorticity data is spatially smoothed along the flow direction and vertical direction with a smoothing radius of 0.05 meters to ensure that local numerical fluctuations do not lead to misjudgment of the vortex core region. The time sampling interval is 0.1 seconds to form a complete vorticity field matrix, providing the input basis for vortex region identification.

[0038] Identify vortex regions based on vorticity fields and determine the rotational intensity of vortex regions; In this embodiment, vorticity fields are used to cluster grid cells with continuous vorticity values ​​exceeding 0.8 m / s, forming vortex candidate regions. By detecting the number of continuous grid cells and the total area within each candidate region, vortex cells smaller than 0.0025 m² are excluded, and the remaining regions are defined as vortex core regions. The center coordinates, number of grid cells, total area, and average vorticity intensity are recorded for each vortex region. The average vorticity intensity is calculated along the core region, in units of m / s. The vortex core area and rotation intensity are used to characterize the vortex scale and dynamic level. The rotation direction of each vortex core region is recorded horizontally and labeled as an angle ranging from 0 to 360 degrees. Further analysis of the vorticity distribution gradient yields the vortex boundary. Boundary grid data records the boundary coordinates and vorticity change rate, ensuring the vortex region is spatially continuous, complete, and accurate.

[0039] Based on the analysis of vortex disturbance dynamics by rotational intensity, vortex data are obtained.

[0040] In this embodiment, a disturbance force distribution matrix is ​​formed based on the vortex core region and rotation intensity, using the rotation intensity and vortex action area. The disturbance force acts along the door leaf normal, measured in kilonewtons. The disturbance force is uniformly distributed within each vortex core region in a grid pattern with a grid spacing of 0.05 meters. The direction of the disturbance force is determined by the vortex rotation direction, with clockwise or counterclockwise corresponding to the force direction. The total disturbance force is calculated for each vortex core region, and its dynamic changes are recorded at a time step of 0.1 seconds. A complete vortex data table is generated, including the center coordinates, core area, rotation intensity, disturbance force magnitude, disturbance force direction, and time series for each vortex core, providing direct input for door leaf vibration detection and vortex-induced vibration quantification.

[0041] Preferably, the calculation of vortex frequency based on vortex data is specifically as follows: ; in, The vortex frequency; The vortex circulation is obtained by integrating the velocity in the core region of the vortex along a closed path. This is the projection of the vortex core area, used to describe the vortex scale; This indicates the effect of core rotation intensity on frequency, which is proportional to the square root of the rotation intensity.

[0042] In this embodiment, velocity field information and core area projection of each vortex core region are extracted from the vortex data. The core region is divided by a spatial grid with a grid spacing of 0.05 meters horizontally and 0.02 meters vertically. Each grid cell records the instantaneous velocity value and its coordinate position in the core region. The vortex circulation is obtained by integrating the velocity of each vortex core region along a closed path. The closed path is defined as an equidistant grid path around the vortex center with a path spacing of 0.05 meters. The instantaneous velocity magnitude and direction are collected at each node on the path, and the tangential component of the velocity vector is accumulated along the path direction during integration. The circulation unit is meters per second, and the integration time step is 0.1 seconds to ensure that the dynamic process is completely captured. The core area projection is calculated through the vortex boundary grid. The number of grids is counted in both the horizontal and vertical directions, and multiplied by the grid cell area to obtain the projected area, in square meters. This area is used to describe the spatial scale of the vortex and provides necessary parameters for vortex frequency calculation. The core rotation intensity is calculated by the average vortex in the core region, in seconds, and recorded in the vortex data table at each time step. The rotational intensity value is square-rooted to obtain the influence factor of the vortex frequency. The vortex frequency is dynamically calculated along the time series, with the frequency corresponding to each vortex core recorded independently, in Hertz. The frequency calculation steps include: combining the circulation value with the core area projection ratio, multiplying by the square root of the core rotational intensity to obtain the instantaneous vortex frequency, with a frequency accuracy of 0.01 Hertz. Each vortex core region is processed iteratively, reading the vortex center coordinates, boundary grid, instantaneous velocity, and vorticity information one by one from the vortex data table, calculating the circulation using the closed-path integration method; obtaining the core area projection using grid statistical methods; reading the rotational intensity and performing square-root processing; finally, calculating the vortex frequency by correlating circulation, area projection, and rotational intensity to form a complete vortex frequency data table, including vortex center coordinates, core area projection, instantaneous circulation, rotational intensity, and corresponding frequency. This data table provides direct input for subsequent door blade excitation force calculations, and complete vortex frequency information is recorded at each time step to ensure that the vortex dynamic process is continuously tracked throughout the entire flood discharge cycle.

[0043] Preferably, the evaluation of the door leaf excitation force based on the vortex frequency is as follows: Match the vortex point of action according to the vortex frequency; In this embodiment, the center coordinates, core area projection, and instantaneous vortex frequency information of each vortex core region are read from the vortex frequency data table. A spatial grid is established along the leaf blade surface, with a grid spacing of 0.05 meters along the width direction and 0.02 meters along the height direction, covering the entire leaf blade surface. Based on the comparison between the vortex frequency and the natural vibration frequency of the leaf blade, grid regions with a frequency difference of less than 0.2 Hz are identified, and these grid nodes are marked as vortex action points. For each action point, its corresponding grid coordinates on the leaf blade surface, the corresponding vortex center coordinates, and the vortex core area projection information are recorded. Through this matching method, the precise spatial mapping of the vortex is achieved, and the state of all action points is recorded at each time step of 0.1 seconds.

[0044] Identify the normal force along the door leaf at the point of action of the vortex; identify the excitation force along the door leaf at the point of action of the vortex; In this embodiment, for each identified vortex point of action, the normal force is calculated along the normal direction of the door leaf surface. The magnitude of the normal force is calculated as a proportion of the square root of the vortex core rotation intensity to the core area projection, in kilonewtons (kN). The rotation intensity in the calculation formula is taken from the vortex data table, and the core area projection is obtained through grid statistical methods. The direction of the normal force is parallel to the door leaf normal. The excitation force is calculated tangentially. The magnitude of the excitation force is calculated from the vortex frequency and the local water flow velocity, fluid density, and local door leaf area corresponding to the point of action, in kilonewtons (kN). The direction of the excitation force is tangential and consistent with the flow direction. The normal force and excitation force at each point of action are recorded with a time step of 0.1 seconds, forming a complete force time series table.

[0045] The instantaneous excitation force at the point of action of the vortex is synthesized from the normal force and the excitation force. In this embodiment, the normal force vector and the tangential excitation force vector at each vortex point of action are vector-superimposed to obtain the instantaneous excitation force vector, including its magnitude and direction, in kilonewtons. Vector superposition is performed according to the corresponding node coordinates, with each point of action calculated independently, and a time step of 0.1 seconds. The instantaneous excitation force is recorded in a data table, including the coordinates of the point of action, the time series, the normal force, the tangential force, and the magnitude and direction of the resultant force. The synthesis process ensures that the dynamic force applied by the vortex is completely described in both space and time at each point of action.

[0046] The instantaneous excitation forces at the vortex action points of each point are spatially superimposed to determine the door leaf excitation force.

[0047] In this embodiment, the instantaneous excitation forces at all vortex action points on the leaf blade surface are spatially superimposed. Vectors are summed along the corresponding positions of the grid nodes to obtain the total excitation force for each grid node at each time step. During the superposition process, a three-dimensional vector summation method is used to accumulate the normal and tangential force components separately, and then calculate the magnitude and direction of the vector composite force. The excitation forces across the entire leaf blade surface form a complete spatial distribution map. With a time step of 0.1 seconds, the instantaneous resultant force values ​​and direction information of all grid nodes are recorded, forming a leaf blade excitation force data table. This data table includes node coordinates, time, instantaneous resultant force magnitude and direction, covering the entire leaf blade surface, achieving a continuous spatial and temporal description of vortex forces, and providing direct input for subsequent structural fatigue analysis and leaf blade life prediction.

[0048] Preferably, step S3 specifically includes: Step S31: Determine the vibration amplitude of the gate leaf based on the gate leaf vibration data; evaluate the plate stiffness using the gate leaf vibration amplitude; identify the unstable area based on the plate stiffness; perform material thinning detection on the unstable area to obtain material thinning data; In this embodiment, vibration acceleration, velocity, and displacement signals are read from the door leaf vibration data table with a time step of 0.01 seconds and a spatial resolution of 0.05 m × 0.05 m grid spacing on the door leaf surface. Peak values ​​of the vibration displacement signals for each grid node are extracted to obtain the node vibration amplitude in millimeters. The local stiffness of the plate surface is calculated using the vibration amplitude data. The plate surface stiffness is defined using bending stiffness in Newton-meters (N·m). The local stiffness of each grid node is equal to the ratio of the node's normal force to the vibration amplitude, with the force taken as the average normal force over a 0.1-second time step. All node local stiffnesses are compared to a set stiffness threshold of 1000 N / mm. Nodes below the threshold are marked as potentially unstable areas. For the marked unstable areas, ultrasonic thickness measurement is used to collect material thickness data every 5 mm along the thickness direction to obtain the actual plate thickness in millimeters. Material thinning data is obtained by subtracting the measured plate thickness from the original designed plate thickness. Nodes with a thickness reduction exceeding 5 mm are recorded as critical material thinning locations. All data is stored according to a spatial grid for easy subsequent processing.

[0049] Step S32: Determine the local collapse degree based on the material thinning data, and determine the degree of damage to the plate surface based on the local collapse degree to obtain the plate surface damage data; In this embodiment, the material thinning data is processed. The local collapse degree is defined according to the ratio of the plate thickness reduction, and the calculation formula is the thinned thickness divided by the original design thickness. The result is a dimensionless value, and each grid node is calculated individually. Grid nodes with a local collapse degree exceeding 0.15 are marked as high-damage nodes. Based on the relationship between local collapse degree and plate instability, the local collapse degree is mapped to the plate damage level, with a damage level range of 0-1, where 0 represents no damage and 1 represents complete failure. The damage levels of all plate nodes are stored in the plate damage data table, including node coordinates, local collapse degree, and plate damage level.

[0050] Step S33: Analyze the frame deformation based on the gate leaf vibration data to obtain the frame deformation data; In this embodiment, vibration displacement signals are extracted from the frame nodes along three-dimensional spatial coordinates with a time step of 0.01 seconds and a spatial resolution of one node every 0.05 meters along the frame. The vibration displacement of the frame nodes is integrated twice to obtain the frame deformation, including deflection values ​​along the X, Y, and Z directions, in millimeters. The maximum deflection of each node is calculated along the length of the frame to obtain the frame deformation amplitude. Stiffness analysis is performed on the frame structure based on the deflection values. Nodes with frame stiffness below a set threshold of 500 N / mm are identified as locally unstable nodes, and the frame deflection and stiffness information is stored in a frame deformation data table.

[0051] Step S34: Integrate the plate surface damage data and plate surface damage data to obtain the door leaf structure fatigue data; In this embodiment, the plate surface damage data table from step S32 and the border deformation data table from step S33 are matched according to the spatial grid node coordinates. The plate surface damage level and border deflection data of the same grid node are merged to form a complete door leaf structure data table. Based on the proportional relationship between plate surface damage level and border deflection, plate surface damage and border deformation are uniformly quantified into fatigue damage indices with a value range of 0-1. The node fatigue damage indices are calculated by linear weighted superposition, with a plate surface weight of 0.6 and a border weight of 0.4, to obtain the door leaf structure fatigue data for each node. All node data are stored according to time series and spatial coordinates to form a complete three-dimensional door leaf structure fatigue distribution dataset.

[0052] Step S35: Predict the fatigue life of the door leaf based on the fatigue data of the door leaf structure.

[0053] In this embodiment, fatigue damage indices for each grid node are extracted from the fatigue data table of the door leaf structure obtained in step S34, with units ranging from 0 to 1. Combined with the door leaf cycle count records, fatigue cycles are accumulated at 0.1-second time steps to calculate the cumulative fatigue factor for each node. According to the fatigue life standard, nodes with a cumulative fatigue factor reaching 1 are marked as fatigue failure nodes. The failure times of each node are statistically analyzed along the entire door leaf surface, and the earliest failure time is taken as the overall fatigue life of the door leaf, with units of hours or flood discharge cycles. The door leaf fatigue life data table includes node coordinates, cumulative fatigue factor, expected failure time, and corresponding time step records, providing a complete data foundation for door leaf life analysis.

[0054] Preferably, step S33 specifically includes: Step S331: Extract vibration signals based on the gate leaf vibration data; In this embodiment, acceleration signals along the gate leaf edge nodes are extracted from the previous gate leaf vibration data. The edge nodes are arranged along the gate leaf edge with a grid spacing of 0.05 meters, and the time step is set to 0.01 seconds. The data unit is meters per second squared (m / s²). The acceleration signal of each node is filtered by a high-pass filter to remove drift components below 0.5 Hz and by a low-pass filter to remove noise above 50 Hz. A fourth-order Butterworth filter is used to ensure clear edge vibration signals. The filtered acceleration signals are stored as a vibration signal data table, containing time series, node numbers, acceleration values, and node three-dimensional coordinate information. The vibration signals of all edge nodes are time-aligned to ensure consistent time steps in subsequent integration calculations.

[0055] Step S332: Perform a second integration based on the vibration signal to obtain the frame displacement data; In this embodiment, the velocity signal is initially integrated using the trapezoidal integration method, with each time step being 0.01 seconds and the initial integration condition set to zero initial velocity. A second integration of the velocity signal yields the displacement signal, and the displacement increments at each time step are accumulated to obtain the total node displacement. To reduce cumulative integration errors, a zero-drift correction technique is employed, adjusting the average displacement value over the entire time period to zero initial displacement. The displacement unit is meters, and the displacement of each node is calculated along the X, Y, and Z directions to obtain complete three-dimensional displacement data for the bounding box. The data is stored in the form of node numbers, time series, and displacement value tables for easy subsequent deflection analysis.

[0056] Step S333: Calculate the deflection of the border based on the border displacement data; In this embodiment, the relative displacement difference of each node is calculated along the length of the border using the border displacement data obtained in step S332, and is used as the local deflection value. Border deflection is defined as the displacement difference between adjacent nodes in the normal direction, in millimeters. The maximum deflection value is extracted along the entire length of the border as the border extreme deflection. The local deflection calculated for each node is smoothed using the three-point difference method, with a smoothing window set to 0.15 meters to eliminate local measurement errors. The border deflection data table includes node coordinates, normal displacement, local deflection, and maximum deflection value, ensuring that it reflects the overall deformation trend of the border.

[0057] Step S334: Evaluate the degree of border deformation based on border deflection to obtain border deformation data.

[0058] In this embodiment, a frame deformation threshold of 5 mm is set, and nodes with local deflection exceeding the threshold are marked as deformed nodes. The degree of frame deformation is quantified according to the ratio of deflection to the threshold, with a dimensionless range of 0 to 1, where the deformation degree is 1 when the deflection equals the threshold. The frame deformation degrees of all nodes are summarized to form a frame deformation data table, including node number, three-dimensional coordinates, local deflection, maximum deflection, and deformation degree value. The data table provides spatial visualization of the overall frame deformation trend, ensuring that the deformation information of each node can be accurately used in subsequent structural fatigue analysis.

[0059] Therefore, the embodiments should be considered as exemplary and non-limiting in all respects, and the scope of the invention is defined by the appended claims rather than the foregoing description. Thus, all variations falling within the meaning and scope of the equivalents of the application are intended to be included within the invention.

[0060] The above description is merely a specific embodiment of the present invention, enabling those skilled in the art to understand or implement the invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the invention. Therefore, the present invention is not to be limited to the embodiments shown herein, but is to be accorded the widest scope consistent with the principles and novel features of the invention herein.

Claims

1. A coupling method for vibration and water flow excitation of a flood discharge and sediment flushing gate, characterized in that, Includes the following steps: Step S1: Obtain the structural drawings of the multi-cylinder hydraulic jacking system and construct a three-dimensional model of the multi-cylinder hydraulic jacking system; use the three-dimensional model of the multi-cylinder hydraulic jacking system to simulate the gate opening and obtain the gate opening parameters; Step S2: Calculate the instantaneous discharge flow using the gate opening parameters and analyze the degree of vortex-induced vibration; perform gate leaf vibration detection based on the degree of vortex-induced vibration to obtain gate leaf vibration data; Step S3: Perform structural fatigue analysis based on the gate leaf vibration data to obtain gate leaf structural fatigue data; predict the gate leaf fatigue life based on the gate leaf structural fatigue data; Step S4: Calculate the instantaneous water pressure based on the gate opening parameters; calculate the stress on the gate leaf plate using the instantaneous water pressure; correct the gate leaf fatigue life based on the stress on the gate leaf plate to obtain the gate life data.

2. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 1, characterized in that, Step S1 is as follows: Step S11: Obtain the drawings of the multi-cylinder hydraulic jacking structure; Step S12: Use the multi-cylinder hydraulic jacking structure drawing to identify the two-dimensional contour and determine the coordinates of the two-dimensional contour; Step S13: Construct a three-dimensional surface patch for multi-cylinder hydraulic jacking based on two-dimensional contour coordinates; Step S14: Assemble the components according to the three-dimensional surface of the multi-cylinder hydraulic jacking to obtain the three-dimensional model of the multi-cylinder hydraulic jacking. Step S15: Use a multi-cylinder hydraulic lifting three-dimensional model to simulate the gate opening and obtain the gate opening parameters.

3. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 1, characterized in that, Step S2 is as follows: Step S21: Extract the effective flow area using the gate opening parameters; Step S22: Obtain the upstream head height; calculate the gate velocity data based on the upstream head height and the effective flow area; Step S23: Perform flow rate integration calculation based on gate flow velocity data to obtain instantaneous discharge flow rate; Step S24: Analyze the degree of vortex-induced vibration using instantaneous discharge flow rate; Step S25: Detect the gate leaf vibration based on the degree of vortex-induced vibration to obtain the gate leaf vibration data.

4. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 2, characterized in that, Step S24 is as follows: Step S241: Calculate the local Reynolds number using the instantaneous discharge flow rate; Step S242: Perform vortex analysis based on the local Reynolds number to obtain vortex data; Step S243: Calculate the vortex frequency based on the vortex data, and evaluate the door leaf excitation force based on the vortex frequency; Step S244: Optimize the degree of vortex-induced vibration based on the door leaf excitation force.

5. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 4, characterized in that, Step S242 is as follows: Extract velocity gradient information based on local Reynolds number; Calculate the vorticity field using velocity gradient information; Identify vortex regions based on vorticity fields and determine the rotational intensity of vortex regions. Based on the analysis of vortex disturbance dynamics by rotational intensity, vortex data are obtained.

6. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 4, characterized in that, The vortex frequency is calculated based on vortex data as follows: ； in, The vortex frequency; The vortex circulation is obtained by integrating the velocity in the core region of the vortex along a closed path. This is the projection of the vortex core area, used to describe the vortex scale; This indicates the effect of core rotation intensity on frequency, which is proportional to the square root of the rotation intensity.

7. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 4, characterized in that, The specific method for evaluating the excitation force of the door leaf based on the vortex frequency is as follows: Match the vortex point of action according to the vortex frequency; Identify the normal force along the door leaf at the point of action of the vortex; identify the excitation force along the door leaf at the point of action of the vortex; The instantaneous excitation force at the point of action of the vortex is synthesized from the normal force and the excitation force. The instantaneous excitation forces at the vortex action points of each point are spatially superimposed to determine the door leaf excitation force.

8. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 1, characterized in that, Step S3 is as follows: Step S31: Determine the vibration amplitude of the gate leaf based on the gate leaf vibration data; evaluate the plate stiffness using the gate leaf vibration amplitude; identify the unstable area based on the plate stiffness; perform material thinning detection on the unstable area to obtain material thinning data; Step S32: Determine the local collapse degree based on the material thinning data, and determine the degree of damage to the plate surface based on the local collapse degree to obtain the plate surface damage data; Step S33: Analyze the frame deformation based on the gate leaf vibration data to obtain the frame deformation data; Step S34: Integrate the plate surface damage data and plate surface damage data to obtain the door leaf structure fatigue data; Step S35: Predict the fatigue life of the door leaf based on the fatigue data of the door leaf structure.

9. The coupling method for vibration and water flow excitation of the flood discharge and sediment flushing gate according to claim 8, characterized in that, Step S33 is as follows: Step S331: Extract vibration signals based on the gate leaf vibration data; Step S332: Perform a second integration based on the vibration signal to obtain the frame displacement data; Step S333: Calculate the deflection of the border based on the border displacement data; Step S334: Evaluate the degree of border deformation based on border deflection to obtain border deformation data.

10. A coupling system for vibration and water flow excitation of a flood discharge and sediment flushing gate, characterized in that, For implementing the coupling method of vibration and water flow excitation of the flood discharge and sediment flushing gate as described in claim 1, the coupling system of vibration and water flow excitation of the flood discharge and sediment flushing gate includes: The gate opening simulation module is used to acquire drawings of the multi-cylinder hydraulic jacking structure and construct a three-dimensional model of the multi-cylinder hydraulic jacking; the gate opening is simulated using the three-dimensional model of the multi-cylinder hydraulic jacking to obtain the gate opening parameters. Gate leaf vibration detection is used to calculate instantaneous discharge flow using gate opening parameters and analyze the degree of vortex-induced vibration; based on the degree of vortex-induced vibration, gate leaf vibration detection is performed to obtain gate leaf vibration data; The gate leaf fatigue life prediction is used to perform structural fatigue analysis based on gate leaf vibration data to obtain gate leaf structural fatigue data; and to predict the gate leaf fatigue life based on the gate leaf structural fatigue data. The instantaneous water pressure calculation module is used to calculate the instantaneous water pressure based on the gate opening parameters; calculate the stress on the gate leaf plate using the instantaneous water pressure; and correct the fatigue life of the gate leaf based on the stress on the gate leaf plate to obtain the gate life data.

Images