CN108595726B - Riverbed adjusting method based on underwater repose angle of sediment - Google Patents
Riverbed adjusting method based on underwater repose angle of sediment Download PDFInfo
- Publication number
- CN108595726B CN108595726B CN201711497942.5A CN201711497942A CN108595726B CN 108595726 B CN108595726 B CN 108595726B CN 201711497942 A CN201711497942 A CN 201711497942A CN 108595726 B CN108595726 B CN 108595726B
- Authority
- CN
- China
- Prior art keywords
- node
- formula
- adjusting
- stage
- riverbed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000013049 sediment Substances 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 63
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 9
- 239000002689 soil Substances 0.000 abstract description 5
- 238000005293 physical law Methods 0.000 abstract description 2
- 238000006467 substitution reaction Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a riverbed adjusting method based on a sediment underwater repose angle, which has the characteristic of adjusting a riverbed in stages, namely the riverbed adjusting process at a certain point is divided into different stages according to the category, different calculation formulas are adopted to adjust the riverbed at each stage, and one stage is only aimed at one to two adjacent nodes, the riverbed adjusting of each stage meets the soil mass conservation, namely the riverbed adjusting method has the conservation, and in a single stage, each adjusting has the simultaneity, and the conservation and the simultaneity make the whole riverbed adjusting process more accord with the physical reality, comprehensively considers the factors such as the mutual influence between the adjacent nodes and the riverbed soil mass conservation, and the like, so that the adjusting process has good physical simulation, the adjusting result more accords with the physical law, and the riverbed evolution simulation efficiency is improved. Although the invention only provides the adjustment mode of a single node, namely the adjacent point, the invention is easy to be applied to all space grid nodes in an expansion mode and can be better combined with the hydraulic calculation program of the structured grid.
Description
Technical Field
The invention relates to a method for solving erosion and deposition deformation of a riverbed surface, in particular to a riverbed adjusting method.
Background
In current scientific research, a method for solving erosion and deposition deformation of a riverbed surface generally comprises the following steps: the mass of the near-bed sediment suspension and the mass of the bed load are calculated, and the solution is carried out on the basis of a sediment mass conservation equation. Numerical simulation of suspended load and bed load is mostly based on empirical formulas, and most of the empirical formulas establish the corresponding relation between silt concentration and bottom flow rate. When solving for high Reynolds number streams using direct numerical modeling (DNS), the flow rates are often random. The substitution of flow rates into empirical formulas results in large variability in both the amount of suspended and mobile mass in space and time. Therefore, the bottom change of the riverbed is calculated in a DNS (domain name system) by using a traditional sediment numerical simulation mode, the riverbed is often rugged, and the rugged bed surface forms generally do not meet the stability requirement of sediment particles in water. The angle of repose of the silt in the water, as a physical quantity determining the maximum slope of the bottom bed pit and the sand ridge, should be taken into full account in the numerical simulation.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide a method for adjusting a stable bed surface of a riverbed based on an underwater repose angle of sediment, aiming at the problems in the prior art.
The technical scheme is as follows: the invention provides a riverbed adjusting method based on an underwater repose angle of sediment, which comprises the following steps:
(1) calculating the gradient of a certain node and an adjacent node, comparing the gradient with the tangent value of the underwater repose angle of the sediment and assigning a slope ratio coefficient;
(2) calculating a classification coefficient according to the slope ratio coefficient, and classifying the nodes by using the classification coefficient;
(3) calculating ideal adjustment heights of adjacent nodes, and sequencing the ideal adjustment heights;
(4) adjusting the height of the riverbed by stages according to the node type;
(5) and repeating the steps for a plurality of times to obtain the riverbed shape meeting the requirement of the underwater repose angle of the sediment.
Further, step (1) first calculates gradient k of a certain i node and adjacent j nodes under the structural gridj:
kj=(yj-yi)/dj (1)
In the formula, yiIs the riverbed elevation of the i node, yjIs the riverbed elevation of the j node, djIs the horizontal distance between the j node and the i node, j is 1,2,3, 4;
the slope coefficient is obtained by the following formula:
ajis the slope ratio coefficient, and alpha is the underwater repose angle of the sediment.
Further, the classification coefficient C of step (2) is calculated by using the following formula:
in the formula, the essential meanings of the calculated cp and cn are the numbers of positive slope ratio coefficient and negative slope ratio coefficient respectively.
Further, the classification of the step (2) node comprises: s-type nodes, D-type nodes and T-type nodes; wherein,
the S type node is that at least one of the classification coefficients is 0, including S0, S1, S2, S3 and S4, the number represents another classification coefficient except 0, and the physical characteristic of the type node is that the unstable gradient type of the adjacent node is only a positive slope or a negative slope;
the D-class node is that at least one of the classification coefficients is 1 and the other classification coefficient is not 0, including D1, D2 and D3, and the number represents another classification coefficient except 1;
the T-type node has classification coefficients with both terms of 2.
Further, the ideal adjustment height in step (3) is calculated by the following formula:
in the formula,the required adjustment height of the j point is shown, the ideal adjustment height is shown by a superscript 0, and the adjustment heights are sorted from small to large.
Further, step (4) represents a full array of adjacent node labels 1,2,3,4 with the corner mark A, B, C, D, which satisfies the ideal adjustment height dyA 0,dyB 0,dyC 0,dyD 0In order from small to large, the node adjusting method comprises the following steps: :
a. if the node is an S0 type node: dyA 0,dyB 0,dyC 0,dyD 0Are all 0, and do not need to be adjusted;
b. if the node is an S-type node: dyA 0,dyB 0,dyC 0,dyD 0The number of nodes other than 0 is related to its subtype, and the S-type node can be calculated in four stages regardless of the subtype, and the adjustment formula of each node is as follows:
in the first stage, the node A is taken as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, yA,yB,yC,yDRespectively representing the height of the riverbed of four nodes adjacent to the node i, wherein the superscript 0 represents the initial height, namely the height before the adjustment method is applied, and the superscript 1 represents the updated riverbed height at the end of the first stage;
and in the second stage, the node B is taken as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, the superscript 2 represents the updated riverbed height at the end of the second stage, and the superscript 1 represents the initial value of the stage;
and in the third stage, the C node is taken as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, the superscript 3 represents the height of the riverbed after the third section is updated, and the superscript 2 represents the initial value of the stage;
and in the fourth stage, the node D is used as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, the superscript 4 represents the riverbed height after the fourth section is updated, and the superscript 3 represents the initial value of the stage;
c. if D1 type node, dyC 0,dyD 0Is not 0, dyA 0,dyB 0At this time, the two nodes C and D need to be adjusted, and there is only one stage for this type of node adjustment:
in the stage one, the C and D nodes are simultaneously taken as basic adjusting nodes, and an adjusting calculation formula is as follows:
d. if D2 type node, dyB 0,dyC 0,dyD 0And if not 0, the number of the adjusting stages is related to the characteristics of the adjacent nodes:
classification is according to the following calculation:
l=(aB 2)aCaD
(d1) if the value of l is 1, the node B is taken as a basic adjusting node in the stage one, and the calculation formula of the adjustment is as follows:
the corresponding slope ratio coefficient a is then changedBWhen the value is 0, S-type node calculation formulas (7) and (8) are adopted for calculation, the value at the end of the calculation stage of the formula (7) is used as the initial value of the formula (7) when the formula (7) is adopted for calculation, and the value at the end of the calculation stage of the formula (7) is used as the initial value of the formula (8) when the formula (8) is adopted for calculation;
(d2) if the value of l is-1, the node B is taken as a basic adjusting node in the stage one, and the calculation formula of the adjustment is as follows:
the corresponding slope ratio coefficient a is then changedBFor 0, the calculation formula of D1 type node is used to calculate the second stage, when the calculation formula (9) is used, the place whose superscript is 0 needs to be changed to 1, and 1 needs to be changed to 2, so as to indicate that the calculation result substituted into the first stage is further adjusted.
e. If D3 type node, dyA 0,dyB 0,dyC 0,dyD 0All the nodes are not 0, and 4 adjacent nodes are required to be adjusted;
(e1) coefficient of slope ratio aADifferent from the other three symbols
In the first stage, the node A is taken as a basic adjusting node, and the adjusting calculation formula is as follows:
the corresponding slope ratio coefficient a is then changedAIs 0, and is calculated by adopting a step calculation method of S-type nodes in sequence, namely the formulas (6), (7) and (8); using the value at the end of the phase when the formula (6) is used for calculation as the initial value of the formula (6), using the value at the end of the phase of the formula (6) when the formula (7) is used for calculation as the initial value of the formula (7), and using the value at the end of the phase of the formula (7) when the formula (8) is used for calculation as the initial value of the formula (8);
(e2) coefficient of slope ratio aAExist as if they were the other three symbols
In the first stage, the node A is taken as a basic adjusting node, and the adjusting calculation formula is as follows:
the corresponding slope ratio coefficient a is then changedAIs 0, and the phase calculation method of D2 type node is adopted to calculate the remaining phase;
f. if the node is a T-type node, the node A is taken as a basic adjusting node in the stage one, and the adjusting calculation formula is as follows:
the corresponding slope ratio coefficient a is subsequently modifiedAWhen the number is 0 and the remaining stage is calculated by using the calculation method of the D2 type node, the place marked with the 0 is changed to 1, and 1 is changed to 2 when the formula (11) or the formula (12) is used for calculation, so as to indicate that the calculation result substituted into the stage one is further adjusted.
Further, the steps (1) to (5) are repeated until a certain number of cycles is reached or finally the inode becomes a node of type S0.
Has the advantages that: the method is mainly applied to rectangular grids, the slope of each space point and adjacent points of the space point are calculated, the underwater repose angle of the sediment is utilized to judge whether the elevation relation between the space point and each point around the space point is stable, the slope of the space point and the adjacent points of the space point and the tangent value of the underwater repose angle of the sediment are compared aiming at the points which do not meet the stability, the space point and the adjacent points of the space point are classified according to the size relation, and different possible riverbed forms are classified; meanwhile, the invention has the characteristic of adjusting the riverbed by stages, namely the riverbed adjusting process at a certain point is divided into different stages according to the category, different calculation formulas are adopted to adjust the riverbed at each stage, one stage only aims at one to two adjacent nodes, the riverbed adjusting at each stage meets the soil mass conservation, namely the riverbed adjusting process has the conservation, and in a single stage, each adjusting has the simultaneity, the conservation and the simultaneity ensure that the adjusting process of the whole riverbed is more consistent with the physical reality, the factors such as the mutual influence between the adjacent nodes, the soil mass conservation of the riverbed and the like are comprehensively considered, so that the adjusting process has good physical simulation, the adjusting result is more consistent with the physical law, and the improvement of the riverbed evolution simulation efficiency is facilitated. Although the invention only provides the adjustment mode of a single node, namely the adjacent point, the invention is easy to be applied to all space grid nodes in an expansion mode and can be better combined with the hydraulic calculation program of the structured grid.
Drawings
FIG. 1 is a diagram of the physical concept of the present invention;
FIG. 2 is a graph of the sum of the maximum ideal heights of the embodiment as a function of the number of cycles;
fig. 3 is a comparison of riverbed elevation cloud pictures, in which (a) is an original picture and (b) is a cloud picture adjusted by applying the present invention.
Detailed Description
The technical solution of the present invention is described in detail below, but the scope of the present invention is not limited to the embodiments.
A riverbed adjusting method based on a silt underwater repose angle defines that the horizontal distances between an i node and four adjacent nodes are unit length 1, the silt repose angle is 45 degrees, and the heights of the i node and each adjacent node are as follows: assuming that the height of the i-node is 5, the height of the i-node and the height of the surrounding nodes are set to be 12.8, 11.2, 5.5 and 8.9, respectively.
The riverbed adjusting process is as follows:
(1) calculating a slope ratio coefficient (a)1,a2,a3,a4,)=(1,1,0,1)。
(2) The classification coefficient C is calculated as (3, 0), and the classification basis is known as the node of type S3.
(3) The ideal adjustment height is calculated and ordered from small to large as follows, and the physical concept of equation (4) is shown in fig. 1: (dy)0 3,dy0 4,dy0 2,dy0 1) When the number of the nodes is equal to (0, 1.45, 2.6, 3.4), the three adjustment stages sequentially adjust the four nodes 3,4, 2, 1.
(4) The node stage adjustment formula of type S3 is given to adjust in stages:
in the first stage, the node 3 and the node 4 are taken as basic adjusting nodes and calculated by adopting the formula (5), and all the obtained items are 0.
And in the second stage, the node No. 4 is taken as a basic adjusting node and the calculation is carried out by adopting the formula (6):
after the substitution, the height of the node in the second stage 4 which needs to be adjusted is calculated to be 0.725, and the height of the node in the second stage after adjustment is (y)i 2,y1 2,y2 2,y3 2,y4 2)=(7.175,12.075,10.475,5.5,8.175)。
And the third stage adopts the node No. 2 as a basic adjusting node and adopts the formula (7) to calculate:
after the result at the end of the second stage is substituted into the formula, the height of the node No. 2 at the third stage, which needs to be adjusted, can be calculated to be 0.767, and the height of each point at the end of the third stage, which is adjusted to be (y)i 3,y1 3,y2 3,y3 3,y4 3)=(8.709,11.308,9.708,5.5,8.175)。
And the fourth stage takes the node No. 1 as a basic adjusting node and adopts the formula (8) to calculate:
after the result at the third stage is substituted into the formula, the height of the node No. 1 at the fourth stage, which needs to be adjusted, can be calculated to be 0.8, and the height of each point at the fourth stage after adjustment is (y)i 4,y1 4,y2 4,y3 4,y4 4)=(9.509,10.508,9.708,5.5,8.175)。
(6) Summing the elevation of each node after adjustment to 43.4, which is equal to the sum of the height of the original riverbed, shows that the invention has conservation of mass.
(7) Repeating the above steps again with the updated node value, first calculating the slope ratio coefficient as (a)1,a2,a3, a4-) (0, 0, -1, -1); then, the classification coefficient C is calculated as (0, 2) node S2 type. Calculating the ideal adjustment height (dy)0 1,dy0 2dy0 4, dy 0 3) Since the four stages adjust the 1,2, 4, and 3 nodes sequentially (0, 0, 0.167, and 1.505).
(8) In the stage I and the stage II, the nodes 1 and 2 are respectively used as basic nodes for adjustment, the elevation of each point is obtained without change after substitution, in the stage III, the node 4 is used as a basic adjustment node and the calculation is carried out by adopting the formula (7), the height of the node 4 in the stage III, which needs to be adjusted, can be calculated to be 0.112, and the height of each point in the stage III, which is adjusted at the end is (y)i 3,y1 3,y2 3,y3 3,y4 3)=(9.285, 10.508,9.708,5.612,8.287)。
(9) Stage four takes the node No. 3 as basic adjustmentThe node is calculated by adopting the formula (8), the height required to be adjusted of the fourth node and the third node in the stage No. 3 can be calculated to be 1.337, and the height after adjustment of each point in the second stage is (y)i 4,y1 4,y2 4,y3 4,y4 4) = (7.948, 10.508, 9.708, 6.949, 8.287). Similarly, the summation of the nodes of each bed surface after adjustment is 43.4, and the conservation of mass is still kept.
By analyzing the calculation results, the mass of the riverbed soil body is always kept conservative in the riverbed adjusting process, and the whole riverbed tends to be gentle along with the continuation of the adjusting process. In addition, the sum of the ideal adjustment heights can be used to determine whether the river bed is smooth. As shown in fig. 2, which is the sum of the ideal adjustment heights at the initial time of each cycle of the above calculation example, it can be seen from the figure that after 6 cycles, the sum of the ideal adjustment heights changes from 7.45 at the beginning to 0.058 and gradually approaches 0, which indicates that the river bed gradually approaches to be flat after adjustment. In addition, the calculation formulas given in each stage of the method ensure that each adjustment stage does not have the condition of over adjustment, namely the phenomenon that the original positive slope is changed into a negative slope due to adjustment does not occur. Finally, the practical application of the method in the calculation of hydraulics is shown by using the graph 3, the graph 3 is a riverbed height cloud chart before and after the method is applied, in the graph 3, a is an original result of the riverbed deformation calculated by using a numerical simulation result, the riverbed height point is abrupt and discontinuous, and in the graph 3, b is a simulation result of the riverbed adjustment performed by applying the method in the numerical calculation, so that the originally uneven riverbed can be found to be more flat after the method is applied, and the rationality of the application result of the method in the practical calculation of hydraulics is illustrated.
Claims (4)
1. A riverbed adjusting method based on a silt underwater repose angle is characterized by comprising the following steps: the method comprises the following steps:
(1) calculating the gradient of a certain node and an adjacent node, comparing the gradient with the tangent value of the underwater repose angle of the sediment and assigning a slope ratio coefficient;
(2) calculating a classification coefficient according to the slope ratio coefficient, and classifying the nodes by using the classification coefficient;
(3) calculating ideal adjustment heights of adjacent nodes, and sequencing the ideal adjustment heights;
(4) adjusting the height of the riverbed by stages according to the node type;
(5) repeating the steps for a plurality of times to obtain a riverbed shape meeting the requirement of the underwater repose angle of the sediment;
firstly, calculating the gradient k of a certain i node and an adjacent j node under a structural gridj:
kj=(yj-yi)/dj (1)
In the formula, yiIs the riverbed elevation of the i node, yjIs the riverbed elevation of the j node, djIs the horizontal distance between the j node and the i node, j is 1,2,3, 4;
the slope coefficient is obtained by the following formula:
ajis the slope ratio coefficient, and alpha is the underwater repose angle of the sediment;
the ideal adjustment height in the step (3) is calculated by adopting the following formula:
in the formula, dyj 0Representing the required adjustment height of the j point, using a superscript 0 to represent an ideal adjustment height, and sequencing the ideal adjustment height from small to large;
said step (4) is represented by the corner mark A, B, C, D as a full array of adjacent node labels 1,2,3,4, which satisfies the ideal adjustment height dyA 0,dyB 0,dyC 0,dyD 0In order from small to large, the node adjusting method comprises the following steps:
a. if the node is an S0 type node: dyA 0,dyB 0,dyC 0,dyD 0Are all 0, and do not need to be adjusted;
b. if the node is an S-type node: dyA 0,dyB 0,dyC 0,dyD 0The number of nodes other than 0 is related to its subtype, and the S-type node can be calculated in four stages regardless of the subtype, and the adjustment formula of each node is as follows:
in the first stage, the node A is taken as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, yA,yB,yC,yDRespectively representing the height of the riverbed of four nodes adjacent to the node i, wherein the superscript 0 represents the initial height, namely the height before the adjustment method is applied, and the superscript 1 represents the height of the riverbed after the first stage is updated;
and in the second stage, the node B is taken as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, the superscript 2 represents the updated riverbed height at the end of the second stage, and the superscript 1 represents the initial value of the stage;
and in the third stage, the C node is taken as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, the superscript 3 represents the height of the riverbed after the third section is updated, and the superscript 2 represents the initial value of the stage;
and in the fourth stage, the node D is used as a basic adjusting node, and an adjusting calculation formula is as follows:
wherein, the superscript 4 represents the riverbed height after the fourth section is updated, and the superscript 3 represents the initial value of the stage;
c. if D1 type node, dyC 0,dyD 0Is not 0, dyA 0,dyB 0At this time, the two nodes C and D need to be adjusted, and there is only one stage for this type of node adjustment:
in the stage one, the C and D nodes are simultaneously taken as basic adjusting nodes, and an adjusting calculation formula is as follows:
d. if D2 type node, dyB 0,dyC 0,dyD 0And if not 0, the number of the adjusting stages is related to the characteristics of the adjacent nodes:
classification is according to the following calculation:
l=(aB 2)aCaD
d1. if the value of l is 1, the node B is taken as a basic adjusting node in the stage one, and the calculation formula of the adjustment is as follows:
the corresponding slope ratio coefficient a is then changedBWhen the value is 0, S-type node calculation formulas (7) and (8) are adopted for calculation, the value at the end of the stage is used as the initial value of the formula (7) when the formula (7) is adopted for calculation, and the value at the end of the stage of the formula (7) is used as the initial value of the formula (8) when the formula (8) is adopted for calculation;
d2. if the value of l is-1, the node B is taken as a basic adjusting node in the stage one, and the calculation formula of the adjustment is as follows:
the corresponding slope ratio coefficient a is then changedBWhen the calculation is performed by adopting the formula (9), the place with the superscript corresponding to 0 needs to be changed into 1, and 1 needs to be changed into 2, so as to indicate that the calculation result substituted into the first stage is further adjusted;
e. if D3 type node, dyA 0,dyB 0,dyC 0,dyD 0All the nodes are not 0, and 4 adjacent nodes are required to be adjusted;
e1. coefficient of slope ratio aAUnlike the other three symbols:
in the first stage, the node A is taken as a basic adjusting node, and the adjusting calculation formula is as follows:
the corresponding slope ratio coefficient a is then changedAIs 0, and is calculated by adopting a step calculation method of S-type nodes in sequence, namely the formulas (6), (7) and (8); using the value at the end of the phase when the formula (6) is used for calculation as the initial value of the formula (6), using the value at the end of the phase of the formula (6) when the formula (7) is used for calculation as the initial value of the formula (7), and using the value at the end of the phase of the formula (7) when the formula (8) is used for calculation as the initial value of the formula (8);
e2. coefficient of slope ratio aAExist as if they were the other three symbols
In the first stage, the node A is taken as a basic adjusting node, and the adjusting calculation formula is as follows:
the corresponding slope ratio coefficient a is then changedAIs 0, and the phase calculation method of D2 type node is adopted to calculate the remaining phase;
f. if the node is a T-type node, the node A is taken as a basic adjusting node in the stage one, and the adjusting calculation formula is as follows:
the corresponding slope ratio coefficient a is subsequently modifiedAWhen the number is 0 and the remaining stage is calculated by using the calculation method of the D2 type node, the place marked with the 0 is changed to 1, and 1 is changed to 2 when the formula (11) or the formula (12) is used for calculation, so as to indicate that the calculation result substituted into the stage one is further adjusted.
2. The method for adjusting a riverbed based on the underwater angle of repose of sediment according to claim 1, wherein the method comprises the following steps: the classification coefficient C of step (2) is calculated as (cp, cn) using the following formula:
in the formula, the essential meanings of the calculated cp and cn are the numbers of positive slope ratio coefficient and negative slope ratio coefficient respectively.
3. The method for adjusting a riverbed based on the underwater angle of repose of sediment according to claim 1, wherein the method comprises the following steps: the classification of the nodes in the step (2) comprises the following steps: s-type nodes, D-type nodes and T-type nodes; wherein,
the S-type node is a node with at least one of classification coefficients of 0, including S0, S1, S2, S3 and S4, the number represents another classification coefficient except 0, and the physical characteristic of the node is that the unstable gradient type of the adjacent node is only a positive slope or a negative slope;
the D-type node is a node with at least one of the classification coefficients being 1 and the other being different from 0, and comprises D1, D2 and D3, and the number represents another classification coefficient except 1;
the T-type node has classification coefficients with both terms of 2.
4. The method for adjusting a riverbed based on the underwater angle of repose of sediment according to claim 1, wherein the method comprises the following steps: and (5) repeating the steps (1) to (5) until the node i finally becomes the node of the type S0.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711497942.5A CN108595726B (en) | 2017-12-29 | 2017-12-29 | Riverbed adjusting method based on underwater repose angle of sediment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711497942.5A CN108595726B (en) | 2017-12-29 | 2017-12-29 | Riverbed adjusting method based on underwater repose angle of sediment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108595726A CN108595726A (en) | 2018-09-28 |
CN108595726B true CN108595726B (en) | 2021-06-11 |
Family
ID=63633112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711497942.5A Active CN108595726B (en) | 2017-12-29 | 2017-12-29 | Riverbed adjusting method based on underwater repose angle of sediment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108595726B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112179614A (en) * | 2020-08-26 | 2021-01-05 | 河海大学 | Bed surface lifting device for real-time silt flushing test of silt |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103198367A (en) * | 2013-04-11 | 2013-07-10 | 中国水电顾问集团成都勘测设计研究院 | Reservoir bank slump predicted value detection method |
CN103308001A (en) * | 2013-05-17 | 2013-09-18 | 河海大学 | Device and method for measuring underwater AoR (Angle of Repose) of sediment based on optical imaging |
CN103913111A (en) * | 2014-03-13 | 2014-07-09 | 河海大学 | Fine-particle sediment underwater repose angle measuring device and measuring method thereof |
CN104315962A (en) * | 2014-11-17 | 2015-01-28 | 重庆交通大学 | Silt underwater repose angle measurement method |
KR101497993B1 (en) * | 2014-08-29 | 2015-03-05 | 연세대학교 산학협력단 | Method and apparatus for analyzing river sedimentation and flushing using quasi-2-dimensional quasi-steady model |
-
2017
- 2017-12-29 CN CN201711497942.5A patent/CN108595726B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103198367A (en) * | 2013-04-11 | 2013-07-10 | 中国水电顾问集团成都勘测设计研究院 | Reservoir bank slump predicted value detection method |
CN103308001A (en) * | 2013-05-17 | 2013-09-18 | 河海大学 | Device and method for measuring underwater AoR (Angle of Repose) of sediment based on optical imaging |
CN103913111A (en) * | 2014-03-13 | 2014-07-09 | 河海大学 | Fine-particle sediment underwater repose angle measuring device and measuring method thereof |
KR101497993B1 (en) * | 2014-08-29 | 2015-03-05 | 연세대학교 산학협력단 | Method and apparatus for analyzing river sedimentation and flushing using quasi-2-dimensional quasi-steady model |
CN104315962A (en) * | 2014-11-17 | 2015-01-28 | 重庆交通大学 | Silt underwater repose angle measurement method |
Non-Patent Citations (3)
Title |
---|
一种塑料模型沙运动特性试验;唐立模 等;《水利水电科技进展》;20121220;第32卷(第6期);全文 * |
泥沙颗粒水下休止角的研究;吕亭豫 等;《价值工程》;20120718;第31卷(第20期);全文 * |
砂岩颗粒料崩解及水下休止角特性试验研究;姬雪竹 等;《科学技术与工程》;20170918;第17卷(第26期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108595726A (en) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110110413B (en) | Structural topology optimization method based on material field reduction progression expansion | |
CN104850682A (en) | Multiple-point geostatistics modeling method based on position | |
CN105303612B (en) | A kind of extract digital network method based on Triangulated irregular network model | |
CN105808311B (en) | A kind of reservoir modeling Fast Fitting method based on Dimension-reduced Strategy | |
CN110929443B (en) | Two-dimensional flood simulation method based on high-precision terrain generalization | |
CN108595726B (en) | Riverbed adjusting method based on underwater repose angle of sediment | |
CN111507026A (en) | Dual-grid multi-scale finite element method for simulating node Darcy permeation flow rate | |
CN106530397B (en) | A kind of geology face three-dimensional rebuilding method based on sparse section geology contour line | |
CN110924935B (en) | Method, device and equipment for determining bottom hole flowing pressure regulation and control scheme of tight oil reservoir | |
CN110990926B (en) | Urban surface building hydrodynamic simulation method based on area correction rate | |
CN103996221A (en) | Plant organ mesh simplification method targeted for visualization calculation | |
CN102385757B (en) | Semantic restriction texture synthesis method based on geometric space | |
CN110245410B (en) | Multi-parametric variable-based topological optimization design method for thermoelastic structure of multiphase material | |
CN108898679B (en) | Automatic labeling method for serial numbers of parts | |
CN104195979A (en) | Riverway intersection water flow stagnant area wedge cone and building method and application thereof | |
CN102054294A (en) | Surface modeling method based on surface theory and optimal control theory | |
CN108509532A (en) | A kind of accumulation method and apparatus applied to map | |
CN102567594A (en) | Method for simulation modeling of offshore island reef type artificial fish reef cluster flow field | |
CN101833597B (en) | Method for acquiring electric field distribution of large-sized complex aircraft by using fast boundary element method | |
CN116956657A (en) | Layered rock slope stability influence factor sensitivity analysis method | |
CN107784141B (en) | Acceleration method for two-dimensional finite control volume calculation | |
CN113688755B (en) | Multi-flow-direction drainage basin feature extraction method based on hexagonal grid | |
CN110135033A (en) | A kind of clear row's effective sand discharge phase length calculation method in muddy water library of storage | |
CN111931263B (en) | Constant flow simulation method based on optimized solution of planar two-dimensional shallow water equation | |
KR101335497B1 (en) | System and method for fluid simulation using multigrid |
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 |