CN110135056B - Rapid automatic analysis method for effective water-blocking ratio of pier in river channel - Google Patents

Rapid automatic analysis method for effective water-blocking ratio of pier in river channel Download PDF

Info

Publication number
CN110135056B
CN110135056B CN201910396508.0A CN201910396508A CN110135056B CN 110135056 B CN110135056 B CN 110135056B CN 201910396508 A CN201910396508 A CN 201910396508A CN 110135056 B CN110135056 B CN 110135056B
Authority
CN
China
Prior art keywords
water
pier
section
blocking
effective
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910396508.0A
Other languages
Chinese (zh)
Other versions
CN110135056A (en
Inventor
杨莉玲
杨芳
余顺超
何用
陈文龙
徐峰俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pearl River Hydraulic Research Institute of PRWRC
Original Assignee
Pearl River Hydraulic Research Institute of PRWRC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Pearl River Hydraulic Research Institute of PRWRC filed Critical Pearl River Hydraulic Research Institute of PRWRC
Priority to CN201910396508.0A priority Critical patent/CN110135056B/en
Publication of CN110135056A publication Critical patent/CN110135056A/en
Application granted granted Critical
Publication of CN110135056B publication Critical patent/CN110135056B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Bridges Or Land Bridges (AREA)

Abstract

The invention discloses a rapid automatic analysis method for effective water blocking ratio of piers in a river channel, and belongs to the technical field of demonstration of river dynamics and hydrodynamic influence of wading buildings. The method mainly comprises the following steps: firstly, preprocessing riverway shoreline terrain, determining an analysis water level and a reference overflow section, and calculating the total area of the overflow section; secondly, performing coordinate transformation by taking the water passing section as an axis, and projecting the top of the pier to the water passing section; and then calculating the projection area of the bridge pier on the water passing section, calculating the effective water blocking ratio of the bridge pier by combining the total area of the water passing section, and drawing a water blocking section diagram. The method has the advantages that the defects that the steps of a method for calculating the effective water-blocking ratio of the pier in the river channel are complex, the efficiency is low, the accuracy is greatly influenced by the skill proficiency of a calculator and the like are overcome, the rapid automatic analysis method for the effective water-blocking ratio of the pier in the river channel is provided, the automatic calculation of the water-blocking ratio of the pier in the river channel is realized, and the efficiency and the precision are improved.

Description

Rapid automatic analysis method for effective water-blocking ratio of pier in river channel
Technical Field
The invention relates to a rapid automatic analysis method for effective water blocking ratio of a pier in a river channel, and belongs to the technical field of demonstration of river dynamics and hydrodynamic influence of wading buildings.
Background
With the rapid development of the transportation industry, bridge construction projects are increasing day by day, the piers in the river channel occupy the water passing section of the river channel and block the water flow, and the ratio of the effective area of the water passing section occupied by the piers to the original water passing section area is the effective water blocking ratio of the piers. The effective water-blocking ratio quantitatively reflects the flow-blocking condition of the bridge pier on water flow, and can also qualitatively reflect the influence of the bridge pier on the flow power of the river channel, and the method has very important significance in the optimization and feasibility demonstration of the bridge construction scheme.
The effective water blocking ratio of the bridge pier has a plurality of influence factors, which are closely related to the form, span, river width, water flow conditions and the like of the bridge pier, and are obviously influenced by the intersection angle of the bridge axis and the water flow, and under the influence of the factors, the water blocking condition of each bridge pier is different and needs to be analyzed one by one. The traditional effective water-blocking ratio of the bridge pier is calculated by adopting a manual method: selecting a reference overflow section according to the main flow direction, and calculating the overflow area under the water level by combining the topography of the river channel and the analysis water level; then according to the form of the bridge pier, making vertical lines of the overflowing section by crossing the edge of the bridge pier one by one, calculating the projection area of the bridge pier by combining the analysis water level, and finally calculating the effective water blocking ratio under the analysis water level. The traditional method is time-consuming and labor-consuming, not only is the calculation efficiency low, but also the calculation efficiency and the accuracy are greatly influenced by the skill proficiency of the calculator, so that the reliability of the calculation result is extremely unstable. In order to meet the requirement of pier water-blocking analysis, a rapid automatic analysis method for the effective water-blocking ratio of piers in a river channel is developed, and it is necessary to realize rapid and accurate calculation of the effective water-blocking ratio of piers.
Disclosure of Invention
The invention aims to overcome the defects that the method for calculating the effective water-blocking ratio of the pier in the river channel is complex in steps, low in efficiency, high in accuracy and greatly influenced by the skill proficiency of a calculator and the like, and provides a rapid automatic analysis method for the effective water-blocking ratio of the pier in the river channel, so that the automatic calculation of the water-blocking ratio of the pier in the river channel is realized, and the efficiency and the precision of the method are improved.
The purpose of the invention is realized by the following technical scheme:
a rapid automatic analysis method for effective water-blocking ratio of piers in a river channel comprises the following steps:
step S1, preprocessing a shoreline of a river reach where bridge piers of the bridge comprising N piers are located and underwater topography in a certain range (comprising all piers and possible water blocking analysis sections), generating a topography triangular net, and setting a water blocking ratio analysis water level Z;
step S2, determining the main stream of the river reach and a water-blocking ratio analysis reference overflow section, and calculating the area S of the overflow sectionGeneral assembly。
The specific process is as follows:
and step S21, obtaining the main flow line of the river reach from three aspects, namely historical actually-measured water flow velocity data, model test or calculation results, and rough judgment from the main trough directly according to the beach trough pattern.
Step S22, making the intersection point of the bridge crossing axis and the main flow line as the perpendicular line of the main flow lineLine segments between the banks (denoted L)AB) Cutting into a river bed to obtain a reference water passing section LABNamely the top line segment of the water passing cross section.
Step S23, calculating the cross-sectional area SGeneral assembly: note LABThe coordinates of the two end points are respectively A (a)1,b1)、B(a2,b2) For the top line segment L of the water sectionABNode encryption is performed using distance DL_jmControlling the magnitude of the encryption, typically taking DL_jmThe minimum width of the bridge pier is less than or equal to.
Step S24, the nodes obtained by encryption are interpolated by the terrain triangulation network to obtain the elevation h of each nodeL_jm(m) wherein the coordinates of each node are represented as (X)L(m),YL(m)), m is the node number, and the total number of nodes is recorded as ML. Combining the analysis water level Z to calculate the water depth H of each nodeL_jm(m), further calculating the cross-sectional area:
HL_jm(m)＝Z-hL_jm(m)，
Figure GDA0002326515650000021
step S3, drawing the plane contour of the pier and the top line segment L of the water passing sectionABAnd performing coordinate conversion for the whole.
The specific process is as follows:
step S31, taking the end point A as an axis, and connecting the bridge pier with the top line segment L of the water passing sectionABRotate for whole, rotate contained angle:
Figure GDA0002326515650000022
step S32, the converted plane takes the water flow line as the vertical y axis and the water passing section LABIs a horizontal x-axis, passing through the water section top line LABThe equation is Y ═ b1(b1A known constant).
And S4, projecting the water blocking end point of each pier to the water passing section in the converted plane system to obtain a projection section set of the effective water blocking area of the pier.
The specific process is as follows:
step S41, in the post-conversion plane system, for each bridge pier: extracting X-axis coordinate X of bridge pier contour line nodeLK(N, i), N, i are bridge pier and bridge pier contour line node numbers, respectively, representing the ith node of the nth bridge pier, N ∈ [1, N]I is any integer greater than 3, and the maximum value is 1000 for simplicity; coordinate value X of X axis of contour line nodeLK(n, i) sorting according to size, searching two nodes with the minimum and maximum X-axis coordinate values, namely the water-blocking end points of the bridge pier, wherein the X-axis coordinate values are respectively recorded as XnminAnd XnmaxProjecting the two end points to the cross section of water, the foot of which can be expressed as Sn(Xnmin,b1)、En(Xnmax,b1)。
Step S42, pressing each pair of vertical feet to SnSorting the values, and then sorting the values according to the foot S of the water blocking end point of each piernAnd EnDetermining a water blocking projection line segment corresponding to each bridge pier, namely SnAnd EnThe N bridge piers are summed to obtain N projected line segments LS1E1～LSNEN；
Step S43, searching and comparing the end points of the N sequenced projection line segments (i.e., N pairs of pendants), and removing the overlapped parts after projection to obtain a set of top line segments (referred to as effective water blocking segments for short) of the effective water blocking regions of the piers on the water cross section, wherein the specific process is as follows (see fig. 2):
firstly, setting an initial value of the total number NR of the effective water blocking sections as N;
next, the N pairs of drooping feet are distinguished one by one, and the Kth pair of drooping feet is used ((S)K,EK)K ∈1, N) are examples:
for each pair of the subsequent (J ═ K +1, N) drooping feet in sequence, the following were judged: sJWhether or not it is between SKAnd EJ-1Is there?
If not, [ S ]K,EJ-1]Projecting the effective water blocking area of the pier on the water passing section to obtain an effective water blocking section, and entering a next pair of drop feet to judge that K is J;
if yes, the total number of the effective water blocking sections is NR-1, the data pointer J is equal to J +1, the judgment is repeated for the next pair of the drop feet, and the operation is stopped when the J is equal to N.
Step S433, obtaining NR effective water-blocking sections [ X ]S(r),XE(r)]The coordinates of the starting point and the terminal point are (X) respectively for NR opposite end pointsS(r),b1)、(XE(r),b1) Where r ∈ 1, NR.
Step S5, calculating the projection area S of the bridge pier on the water passing sectionResistance device。
The specific implementation process is as follows:
step S51, inverse transforming the start-stop end point coordinates of the NR effective water blocking sections to the original coordinate system;
and (3) coordinate inverse transformation of each effective water blocking section:
the starting endpoint transformation formula is:
Figure GDA0002326515650000031
the termination endpoint transformation formula is:
Figure GDA0002326515650000041
s52, node encryption is carried out on each projection line segment under the original coordinate system, and the total number of nodes of each segment is recorded as NVrBy distance DR_jm(r) controlling the encryption amplitude by taking DR_jm(r) is less than or equal to the width of the corresponding bridge pier/2.0;
s53, interpolating each node obtained by encrypting NR projection line segments under the original coordinates by a terrain triangulation network, and marking the coordinates of each node as (X)R(r,f),YR(r,f)，r＝1,…,NR，f＝1,…,NVr) To calculate the elevation h of each nodeR_jm(r, f), calculating the water depth H of each node by combining the analysis water level ZR_jm(r, f), and further calculating to obtain the area of the effective water blocking area, wherein the formula is as follows:
HR_jm(r,f)＝Z-hR_jm(r,f)
Figure GDA0002326515650000042
Figure GDA0002326515650000043
and step S6, calculating and analyzing the effective water blocking ratio of the pier at the Z water level, and drawing a water blocking sectional view.
The specific process is as follows:
step S61, according to the effective water blocking ratio calculation formula of the bridge pier: rIs effective＝SResistance device/SGeneral assembly× 100%, calculating the effective water-blocking ratio under the analysis water level Z.
Step S62, according to the top line segment L of the water passing sectionABNode coordinate (X)L(m),YL(m), m ∈ 1, ML) calculates the distance dis (m) of each node to the starting end point a:
the distance DIS (1) from the starting point is equal to 0;
distance between nodes 2 to NL:
Figure GDA0002326515650000044
step S63, according to the top line segment L of the water passing sectionABDistance-elevation ((DIS (m), H)L_jm(m)), m ∈ 1, ML) information draws a cross-sectional view;
step S64, determining each endpoint ((X) of NR effective water-blocking segmentsS(r),XE(r)), distance-elevation information of r ∈ 1, NR), wherein the distance of each end point is the x coordinate value of each end point under the conversion coordinate system, and the elevation h of each end point of the water blocking section isResistance deviceThe design of the top elevation of the pier is realized, and any constant higher than the analysis water level can be selected.
And step S65, drawing a water blocking section data series according to the distance-elevation information, adding the water blocking section data series to the water passing section profile, and generating the water blocking section profile.
Compared with the prior art, the invention has the following advantages and effects:
the invention overcomes the defects that the method for calculating the effective water-blocking ratio of the pier in the river channel is complex in steps, low in efficiency, obvious in accuracy influenced by the skill proficiency of a calculator and the like, provides the rapid automatic analysis method for the effective water-blocking ratio of the pier in the river channel, realizes the automatic calculation of the water-blocking ratio of the pier in the river channel, and improves the efficiency and the precision of the automatic calculation method.
Drawings
FIG. 1 is a step diagram of a rapid automatic analysis method for effective water blocking ratio of a pier in a river channel, which is disclosed by the invention;
FIG. 2 is a flow chart of the method for analyzing the overlapping property of each pair of pier projection line segments to obtain an effective water-blocking segment;
FIG. 3 is a plan view of a pier in accordance with an embodiment of the present invention;
FIG. 4 is an example engineered water underwater topography and scatter plot;
FIG. 5 is an underwater geogrid produced by an engineered body of water in an example;
FIG. 6 is an example flow cross section LABA position map;
FIG. 7 shows the example with the end point A of the water-passing section as the vertex and L as the vertexABA coordinate transformation diagram for the axis;
FIG. 8 is a diagram showing the correspondence between pier contour line nodes No. 6 and projected foot in the example;
FIG. 9 is a projection line segment and detail diagram of all piers in the example;
FIG. 10 is a cross-sectional pre-engineering flow chart of an example;
fig. 11 is a cross-sectional view of pier water blocking in the example.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer and clearer, the present invention is further described in detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Example (b):
effective water-blocking ratio analysis is carried out aiming at a certain bridge design scheme of the Banana south branch water channel, in the embodiment, 2 groups of main piers are distributed on two sides of a main groove close to a west bank, 3 groups of side piers are arranged in a shoal area close to an east bank, all piers are respectively arranged in a separated type up-down amplitude division mode, the plane of a single main pier is rectangular, the size is 9.5 × 14.5.5 m, the plane of the single side pier is a rounded rectangle, the size is 6 × 10m, the designed top elevation of the piers is 13.5m, the plane arrangement of the piers is shown in figure 3, and the set analysis water level is 0 m.
As shown in fig. 1, the method for rapidly and automatically analyzing the effective water blocking ratio of the pier in the river channel specifically comprises the following steps:
step S1, pre-processing the shoreline of the nana-gate river reach where the bridge pier is planned to be built and the underwater topography within a certain range (including all bridge piers and possible water blocking analysis sections, in this example, the range from 238m upstream to 302m downstream of the bridge) to obtain the underwater topography and scattering points of the engineering water area as shown in fig. 4, and generating a topography triangulation network as shown in fig. 5;
step S2, determining a main flow line of a river reach where the project is located and a water-blocking ratio analysis reference overflow section, and calculating the area S of the overflow sectionGeneral assembly。
Step S21, determining the main flow line of the river reach according to the historical flow rate data (mainly referring to the main trough water area), wherein the main flow line of the engineering river reach about 14.68 ° south east.
Step S22, the intersection point of the bridge crossing axis and the main flow line is taken as the vertical line of the main flow line, and the intersection points of the vertical line and the embankment at the two sides of the river bank are respectively marked as A (a)1,b1)、B(a2,b2) Over-line section LABCutting into a river bed to obtain a reference water passing section, and taking its top line L in a plan viewABAs shown in fig. 6. Line segment LABThe coordinates of each end point are known quantities, in this example, a1＝38456410.8002；b1＝2512955.3689；a2＝38457253.0768；b2＝2513175.9547。
Step S23, calculating the cross-sectional area SGeneral assembly: to the top line segment L of the water passing sectionABNode encryption is performed using distance DL_jmControlling the magnitude of the encryption, typically taking DL_jmNot more than the minimum width of the bridge pier, in this example the minimum width of the bridge pier is 6m, and the line segment LABThe uniform encryption is 201(NL) nodes (coordinate is recorded as (X)L(i),YL(i) I 1,201), the encryption amplitude DL_jmAbout 4.3534 m.
Step S24, each node obtained by encryption is interpolated by the terrain triangulation network, wherein the coordinate of each node is recorded as (X)L(m),YL(m)), m is the node number, node totalRecording the number as ML, and calculating the elevation h of each nodeL_jm(m), the water depth H below 0m can be obtainedL_jm(m), further calculating the cross-sectional area:
HL_jm(m)＝0-hL_jm(m)，
Figure GDA0002326515650000061
calculating the cross-sectional area S of the engineering position in the embodimentGeneral assemblyIs 4171.71m2。
Step S3, drawing the plane contour of the pier and the top line segment L of the water passing sectionABThe coordinate conversion is performed for the whole as shown in fig. 7.
Step S31, taking the end point A as an axis, and connecting the bridge pier with the top line segment L of the water passing sectionABRotate for whole, rotate contained angle:
Figure GDA0002326515650000071
this angle is equal to the angle of the main flow line from the forward south direction, 14.68.
Step S32, the converted plane takes the water flow as the vertical y axis and the water passing section as LABHorizontal x-axis, cross-sectional water line LABThe equation is Y ═ b1＝2512955.3689。
And S4, projecting the water blocking end point of each pier to the water passing section in the converted plane system to obtain a projection section set of the effective water blocking area of the pier.
Step S41, in the post-conversion plane system, for each bridge pier: extracting X-axis coordinate (X) of bridge pier contour line nodeLK(n, i), wherein n and i respectively represent the ith node of the nth pier, the number of piers in the river channel is 10 in total in the embodiment, the maximum value of n is 10, i is any integer larger than 3, and the maximum value of i is 100 in the embodiment. Coordinate value X of X axis of contour line nodeLK(n, i) sorting according to size, searching two nodes with the minimum and maximum X-axis coordinate values, namely the water-blocking end points of the bridge pier, wherein the X-axis coordinate values are respectively recorded as XnminAnd XnmaxProjecting the two end points to the cross section of water, the foot of which can be expressed as Sn(Xnmin,b1)、En(Xnmax,b1)。
Taking number 6 pier with more complex pier outline as an example, the coordinate X of each nodeLK(n, i) is shown in Table 1.
Table 16 pier outline node X coordinate
Figure GDA0002326515650000072
Through sequencing analysis, the 35 th node X with the minimum X value in each node of the No. 6 pier outline isLK(6,35), the 13 th node X with the maximum X valueLK(6,13) then there is X6min＝XLK(6,35)＝38457092.24，X6max＝XLKWhen the number 6 pier projection foot is S, the number 13 is 38457098.536(38457092.24,2512955.3689)、E6(38457098.53,2512955.3689) as shown in FIG. 8.
In step S42, if 10 piers are located in the river (N is 10), there are 10 pairs of vertical legs ((S)n,En) N is 1,10), pressing each pair of feet to SnSorting the values in order to obtain 10 projected line segments LS1E1、LS2E2、LS3E3、LS4E4、LS5E5、LS6E6、LS7E7、LS8E8、LS9E9、LS10E10(subscripts S1, E1 are line segment endpoints, i.e., drop legs, and the like), but these line segments may overlap each other and not all form an effective water-blocking section, e.g., LS1E1And LS2E2There is an overlapping part line segment (S)2,E1) As shown in fig. 9.
Step S43, searching and comparing start and stop end points (i.e. all feet) of 10 projection line segments, and removing overlapped parts after projection to obtain a set of top line segments (called effective water blocking segments for short) of the actual water blocking area of the pier on the water cross section, wherein the specific process is as follows, as shown in fig. 2:
firstly, setting the total number NR of the effective water blocking sections as an initial value NR to 10;
next, 10 pairs of drop feet are distinguished one by one:
for (S)1,E1) Judgment, S2Between S1And E1Is there? Is, then LS1E2The effective water blocking sections are effective water blocking sections, the total number NR of the effective water blocking sections is 9, and the data pointer J is 3;
for (S)3,E3) Judgment, S4Between S3And E3Is there? Is, then LS3E4The effective water blocking sections are effective water blocking sections, the total number NR of the effective water blocking sections is 8, and the data pointer J is 5;
for (S)5,E5) Judgment, S6Between S5And E5Is there? Is, then LS5E6The effective water blocking sections are effective water blocking sections, the total number NR of the effective water blocking sections is 7, and the data pointer J is 7;
for (S)7,E7) Judgment, S8Between S7And E7Is there? Is, then LS7E8The effective water blocking sections are effective water blocking sections, the total number NR of the effective water blocking sections is 6, and the data pointer J is 9;
for (S)9,E9) Judgment, S10Between S9And E9Is there? Is, then LS9E10For the effective water-blocking segment, the total number NR of the effective water-blocking segments is 5, the data pointer J is 11, and J>And N, ending the judgment program.
Finally, 5 effective water-blocking sections, L, are obtainedS1E2、LS3E4、LS5E6、LS7E8、LS9E10，NR＝5。
Step S5, calculating the projection area S of the bridge pier on the water passing sectionResistance device；
Step S51, inverse transforming the coordinates of the start and stop end points of the 5 effective water blocking sections to the original coordinate system;
and (3) coordinate inverse transformation of each effective water blocking section:
the starting endpoint transformation formula is:
Figure GDA0002326515650000081
the termination endpoint transformation formula is:
Figure GDA0002326515650000082
in this example, NR is 5, a1＝38456410.8002；b1＝2512955.3689；a2＝38457253.0768；b2＝2513175.9547，θ＝14.68°。
Step S52, node encryption is performed on r (═ 5) projection line segments in the original coordinate system (the total number of nodes per segment is denoted by NV)r R 1,6)), the distance D is usedR_jm(r) controlling the encryption amplitude, DR_jmThe line segment (r) is less than or equal to the width of the corresponding bridge pier/2.0, each projection segment is divided into 5 segments in the example, and the NV is unifiedr＝6，DR_jmThe value of (r) is 3.362m, 3.362m, 1.51m, 1.512m and 1.512m respectively.
Step S53, for 5 projected line segments under the original coordinates, each node (coordinates are expressed as (X) by the terrain triangulation networkR(r,f),YR(r, f)), r ∈ 1,5 and f ∈ 1,6) to obtain the water depth H under the 0m analysis water level as the elevation of each nodeR_jm(r, f), further calculating to obtain the area of the effective water blocking area:
HR_jm(r,f)＝Z-hR_jm(r,f)
Figure GDA0002326515650000091
Figure GDA0002326515650000092
Sresistance device＝111.37+60.62+14.71+17.42+10.4＝214.52(m2)
And step S6, calculating the effective water-blocking ratio of the pier, and drawing a water-blocking sectional view.
Step S61, according to the effective water blocking ratio calculation formula of the bridge pier: rIs effective＝SResistance device/SGeneral assembly× 100, calculating the effective water-blocking ratio of the pier at the water level of 0m as follows:
Ris effective＝214.52/4171.71×100％＝5.14％。
Step S62, according to the top line segment L of the water passing sectionABNode coordinates ((X)L(m),YL(m)), m ∈ 1,201 1,201) calculating each node up toDistance dis (m) (m ∈ 1,201 1,201) of start point a (west bank):
the distance DIS (1) from the starting point is equal to 0;
distance between nodes 2 to 201 (NL):
Figure GDA0002326515650000093
step S63, according to the top line segment L of the water passing sectionABDistance-elevation ((DIS (m), H)L_jm(m)), m ∈ 1, NL) information draws a cross-sectional view through water, as shown in fig. 10;
step S64, determining 5 effective water-blocking sections LS1E2、LS3E4、LS5E6、LS7E8、LS9E10The distance between each end point is the X value of the end point under the transformation coordinate system, such as line segment LS5E6Endpoint S5、E6The distance is the foot S taking the cross section of the water as the X axis5、E6X value of (a), elevation information (h) of all water-blocking sectionsResistance deviceDesigning a top elevation for the pier, which is 13.5m in the example, so that a constant higher than the analysis water level can be simplified and taken, and 5m is taken in the example;
step S65, drawing a water blocking section data series according to the distance-elevation information, adding the water blocking section data series to the water passing section profile to generate a water blocking section profile, as shown in FIG. 8.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.

Claims (6)

1. A rapid automatic analysis method for effective water-blocking ratio of piers in a river channel is characterized by comprising the following steps:
s1, preprocessing a shoreline of a river reach where a bridge comprising N piers is located and an underwater terrain within a certain range to generate a terrain triangulation network, and setting a water blocking ratio analysis water level Z, wherein the certain range refers to an area comprising all piers and possible water blocking analysis sections;
s2, determining a main flow line of a river reach where the pier is located and a water-blocking ratio analysis reference overflow section, and calculating the area S of the overflow sectionGeneral assembly；
S3, drawing the plane contour of the pier and the top line segment L of the water passing sectionABPerforming coordinate conversion as a whole;
s4, in the converted plane system, projecting the water blocking end point of each pier to the water passing section to obtain a projection section set of the effective water blocking area of the pier, wherein the process is as follows:
s41, extracting X-axis coordinates X of bridge pier contour line nodes for each bridge pier in the converted plane systemLK(N, i), N, i are bridge pier and bridge pier contour line node numbers, respectively, representing the ith node of the nth bridge pier, N ∈ [1, N]I is any integer greater than 3, and the maximum value is 1000 for simplicity; coordinate value X of X axis of contour line nodeLK(n, i) sorting according to size, searching two nodes with the minimum and maximum X-axis coordinate values, namely the water-blocking end points of the bridge pier, wherein the X-axis coordinate values are respectively recorded as XnminAnd XnmaxProjecting the two end points to the cross section of water, the foot of which can be expressed as Sn(Xnmin,b1)、En(Xnmax,b1)；
S42, pressing each pair of vertical feet to SnSorting the values, and then sorting the values according to the foot S of the water blocking end point of each piernAnd EnDetermining a water blocking projection line segment corresponding to each bridge pier, namely SnAnd EnThe N bridge piers are summed to obtain N projected line segments LS1E1～LSNEN；
S43, searching and comparing the end points of the sequenced N projection line segments, namely N pairs of pendulous feet, and removing the overlapped parts after projection to obtain a top line segment set of an effective water blocking area of the pier on the water passing section;
s5, calculating the projection area S of the bridge pier on the water passing sectionResistance device；
S6, calculating and analyzing effective water-blocking ratio R of pier under Z water levelIs effectiveDrawing a water-blocking profile, whereinThe effective water-blocking ratio calculation formula of the pier is as follows: rIs effective＝SResistance device/SGeneral assembly×100％。
2. The method according to claim 1, wherein the step S2 comprises determining a main stream line of a river course where the pier is located and a water blocking ratio analysis reference flow cross section, and calculating a flow cross section area SGeneral assemblyThe process of (2) is as follows:
s21, determining a main flow line of the river reach, wherein the method for determining the main flow line of the river reach comprises the following steps: (1) measuring water flow velocity data through history actual measurement, (2) passing through model test or calculation results, (3) judging according to the beach tank pattern and then by the main tank trend;
s22, making the intersection point of the bridge-crossing axis and the main flow line as the vertical line of the main flow line, making a section from the line segment of the vertical line between the river banks to the river bed, and obtaining the section as the reference water passing section, wherein the line segment of the vertical line between the river banks is marked as LABThe top line segment of the water passing section is obtained;
s23, calculating the cross-sectional area SGeneral assembly: note LABThe coordinates of the two end points are respectively A (a)1,b1)、B(a2,b2) For the top line segment L of the water sectionABNode encryption is performed using distance DL_jmControlling the magnitude of the encryption, taking DL_jmThe minimum width of the bridge pier is less than or equal to;
s24, interpolating each node obtained by encryption by the terrain triangulation network to obtain the elevation h of each nodeL_jm(m) wherein the coordinates of each node are represented as (X)L(m),YL(m)), m is the node number, and the total number of the nodes is recorded as ML; combining the analysis water level Z to calculate the water depth H of each nodeL_jm(m), calculating the area of the cross section:
HL_jm(m)＝Z-hL_jm(m)，
Figure FDA0002423260670000021
3. according to the rightThe method for rapidly and automatically analyzing the effective water blocking ratio of the pier in the river channel according to claim 2, wherein in the step S3, a plane contour diagram of the pier and a top line segment L of a water passing section are usedABThe process of performing the coordinate conversion as a whole is as follows:
s31, taking the end point A as an axis, and connecting the pier with the top line segment L of the water passing sectionABRotate for whole, rotate contained angle:
Figure FDA0002423260670000022
s32, taking the water flow as the vertical y axis and the water passing section L as the converted plane coordinate systemABIs a horizontal x-axis, passing through the water section top line LABEquation expressed as Y ═ b1，b1Is a known constant.
4. The method for rapidly and automatically analyzing the effective water blocking ratio of piers in a river channel according to claim 1, wherein the process of step S43 is as follows:
s431, setting the total number NR of the effective water blocking sections as an initial value N;
s432, the N pairs of the drooping feet are judged one by one, and the Kth pair of the drooping feet is used (S)K,EK) K ∈ 1, N is an example, the pair of the subsequent (J-K +1, N) pair of the foot hanging in turn is judged, and S is judgedJWhether or not it is between SKAnd EJ-1In the above-mentioned manner,
if not, [ S ]K,EJ-1]Projecting the effective water blocking area of the pier on the water passing section to obtain an effective water blocking section, and entering a next pair of drop feet to judge that K is J;
if yes, the total number of the effective water blocking sections is NR-1, a data pointer J is equal to J +1, the judgment is repeated for the next pair of feet, and the judgment is ended when J is equal to N;
s433, obtaining NR effective water-blocking sections [ X ]S(r),XE(r)]The coordinates of the starting point and the terminal point are (X) respectively for NR opposite end pointsS(r),b1)、(XE(r),b1) Where r is 1, …, NR.
5. According to claim 1The method for quickly and automatically analyzing the effective water blocking ratio of the pier in the river channel is characterized in that in the step S5, the projection area S of the pier on the water passing section is calculatedResistance deviceThe process of (2) is as follows:
s51, converting the start-stop end point coordinates of the NR effective water blocking sections into an original coordinate system in an inverse manner;
and (3) coordinate inverse transformation of each effective water blocking section:
the starting endpoint transformation formula is:
Figure FDA0002423260670000031
r∈1,NR；
the termination endpoint transformation formula is:
Figure FDA0002423260670000032
r∈1,NR；
s52, node encryption is carried out on each projection line segment under the original coordinate system, and the total number of nodes of each segment is recorded as NVrBy distance DR_jm(r) controlling the encryption amplitude by taking DR_jm(r) is less than or equal to the width of the corresponding bridge pier/2.0;
s53, interpolating each node obtained by encrypting NR projection line segments under the original coordinates by a terrain triangulation network, and marking the coordinates of each node as (X)R(r,f),YR(r,f)，r＝1,…,NR，f＝1,…,NVr) To calculate the elevation h of each nodeR_jm(r, f), calculating the water depth H of each node by combining the analysis water level ZR_jm(r, f), and further calculating to obtain the area of the effective water blocking area, wherein the formula is as follows:
HR_jm(r,f)＝Z-hR_jm(r,f)
Figure FDA0002423260670000033
Figure FDA0002423260670000034
6. the method for rapidly and automatically analyzing the effective water-blocking ratio of piers in a river channel according to claim 1, wherein the process of calculating the effective water-blocking ratio of piers and drawing a water-blocking profile in step S6 is as follows:
s61, according to the effective water blocking ratio calculation formula of the pier: rIs effective＝SResistance device/SGeneral assembly× 100%, calculating the effective water-blocking ratio under the analysis water level Z;
s62, according to the top line segment L of the water passing sectionABNode coordinate (X)L(m),YL(m), m ∈ 1, ML) calculates the distance dis (m) of each node to the starting end point a:
the distance DIS (1) from the starting point is equal to 0;
distance between nodes 2 to NL:
Figure FDA0002423260670000041
s63, according to the top line segment L of the water passing sectionABDistance-elevation (DIS (m), H)L_jm(m) information drawing a cross-sectional view of the water;
s64, determining each endpoint of NR effective water-blocking sections ((X)S(r),XE(r)), distance-elevation information of r ∈ 1, NR), wherein the distance of each end point is the x coordinate value of each end point under the conversion coordinate system, and the elevation h of each end point of the water blocking section isResistance deviceDesigning a top elevation for the bridge pier, and selecting any constant higher than the analysis water level;
and S65, drawing a water blocking section data series according to the distance-elevation information, adding the water blocking section data series to the water passing section profile, and generating the water blocking section profile.
CN201910396508.0A 2019-05-14 2019-05-14 Rapid automatic analysis method for effective water-blocking ratio of pier in river channel Active CN110135056B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910396508.0A CN110135056B (en) 2019-05-14 2019-05-14 Rapid automatic analysis method for effective water-blocking ratio of pier in river channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910396508.0A CN110135056B (en) 2019-05-14 2019-05-14 Rapid automatic analysis method for effective water-blocking ratio of pier in river channel

Publications (2)

Publication Number Publication Date
CN110135056A CN110135056A (en) 2019-08-16
CN110135056B true CN110135056B (en) 2020-06-30

Family

ID=67573772

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910396508.0A Active CN110135056B (en) 2019-05-14 2019-05-14 Rapid automatic analysis method for effective water-blocking ratio of pier in river channel

Country Status (1)

Country Link
CN (1) CN110135056B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111597643B (en) * 2020-05-28 2023-08-04 中国水利水电夹江水工机械有限公司 Modeling analysis method suitable for stress and deformation of crane trolley
CN113804916B (en) * 2021-09-17 2022-09-30 河海大学 Frequency domain space-time image velocity measurement method based on maximum flow velocity prior information
CN114397476B (en) * 2021-11-15 2022-10-14 河海大学 Flow velocity effectiveness identification and correction method for frequency domain space-time image velocity measurement
CN115293072B (en) * 2022-09-29 2023-02-10 珠江水利委员会珠江水利科学研究院 Water-blocking analysis method, device, equipment and medium based on estuary open water area bridge
CN116341237B (en) * 2023-03-17 2023-12-19 长江水利委员会水文局长江中游水文水资源勘测局 Wading building water blocking analysis method and system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6057516A (en) * 1983-09-09 1985-04-03 Hitachi Ltd Magnetic head
JP5088650B2 (en) * 2008-05-19 2012-12-05 東洋建設株式会社 Reinforcement method and structure of existing pile vent pier
CN106978789B (en) * 2017-04-10 2018-12-14 李彬 A kind of bridge based on river energy blocks water on river and influences method of discrimination
CN107356285B (en) * 2017-06-16 2019-12-24 珠江水利委员会珠江水利科学研究院 River mouth bridge engineering water-blocking ratio measuring method and device
CN109060056B (en) * 2018-08-20 2019-08-27 长江水利委员会长江科学院 A kind of river cross-section method of calculating flux of contactless radar flow measurement

Also Published As

Publication number Publication date
CN110135056A (en) 2019-08-16

Similar Documents

Publication Publication Date Title
CN110135056B (en) Rapid automatic analysis method for effective water-blocking ratio of pier in river channel
JP6736731B2 (en) Hydrological model watershed scale determination method based on energy process similarity
CN112417573B (en) GA-LSSVM and NSGA-II shield tunneling multi-objective optimization method based on existing tunnel construction
CN110570031A (en) Method for predicting river channel split ratio of two-stage branched river network in plain
CN104652347B (en) Mountain area non-static water body water level affects population assessed in relation method with flooding
CN111898303A (en) River basin water level and waterlogging forecasting method based on weather forecasting and hydrodynamic simulation
CN110906992A (en) River flow measuring method based on horizontal ADCP measuring vertical line flow velocity distribution
CN106649987B (en) A kind of anemometer tower sets up the quantitative analysis method of scheme
CN104933268A (en) Flood analyzing method based on one-dimensional unsteady flow numerical model
CN116151013A (en) Method for pushing out design flood of small-river-basin urban river channel
CN117454724B (en) Method for determining geometric parameters of urban river type flood control dike through CFD
CN109063370B (en) Method for predicting water depth average water flow direction of curved riverway of flood plain
CN110263428B (en) Riverbed evolution analysis method based on flow weighted average streamline length index
CN115950679B (en) River sediment sampling method
CN110751726B (en) River engineering quality detection method
CN113327323B (en) Water body environment terrain construction method based on scatter data
CN113609665B (en) Method for estimating flood of power engineering in stagnant flood storage area under data-free condition
CN115293037A (en) Hydrodynamic modeling method of river and lake composite system based on machine learning
CN111611641A (en) Port construction method considering residual flow
Liu et al. Comparison and analysis of calculation of Bridge backwater based on Mike21 hydrodynamic model
He et al. Modeling bedload transport trajectories along a Sine-generated channel
CN116341283B (en) Underwater topography simulation method and system based on deep water silt sampling
CN116467790B (en) Discrimination method for optimizing plane layout of breakwater in half-shelter harbor area
CN112948915B (en) Generalization processing method for vertical wading building in numerical simulation test
CN114169155B (en) Method for calibrating representative line by depending on section flow velocity

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant