CN108595726B - Riverbed adjusting method based on underwater repose angle of sediment - Google Patents

Riverbed adjusting method based on underwater repose angle of sediment Download PDF

Info

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
Application number
CN201711497942.5A
Other languages
Chinese (zh)
Other versions
CN108595726A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201711497942.5A priority Critical patent/CN108595726B/en
Publication of CN108595726A publication Critical patent/CN108595726A/en
Application granted granted Critical
Publication of CN108595726B publication Critical patent/CN108595726B/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
    • 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)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (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

Riverbed adjusting method based on underwater repose angle of sediment
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:
Figure BDA0001531983520000021
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:
Figure BDA0001531983520000022
Figure BDA0001531983520000023
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 content of the first and second substances,
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:
Figure BDA0001531983520000024
in the formula (I), the compound is shown in the specification,
Figure BDA0001531983520000025
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:
Figure BDA0001531983520000031
Figure BDA0001531983520000032
Figure BDA0001531983520000033
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:
Figure BDA0001531983520000034
Figure BDA0001531983520000035
Figure BDA0001531983520000036
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:
Figure BDA0001531983520000037
Figure BDA0001531983520000038
Figure BDA0001531983520000039
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:
Figure BDA00015319835200000310
Figure BDA00015319835200000311
Figure BDA00015319835200000312
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:
Figure BDA0001531983520000041
Figure BDA0001531983520000042
Figure BDA0001531983520000043
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:
Figure BDA0001531983520000044
Figure BDA0001531983520000045
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:
Figure BDA0001531983520000046
Figure BDA0001531983520000047
Figure BDA0001531983520000048
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:
Figure BDA0001531983520000051
Figure BDA0001531983520000052
Figure BDA0001531983520000053
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:
Figure BDA0001531983520000054
Figure BDA0001531983520000055
Figure BDA0001531983520000056
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:
Figure BDA0001531983520000057
Figure BDA0001531983520000058
Figure BDA0001531983520000059
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):
Figure BDA0001531983520000071
Figure BDA0001531983520000072
Figure BDA0001531983520000073
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:
Figure BDA0001531983520000074
Figure BDA0001531983520000075
Figure BDA0001531983520000076
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:
Figure BDA0001531983520000077
Figure BDA0001531983520000078
Figure BDA0001531983520000079
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:
Figure FDA0003013478100000011
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:
Figure FDA0003013478100000012
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:
Figure FDA0003013478100000021
Figure FDA0003013478100000022
Figure FDA0003013478100000023
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:
Figure FDA0003013478100000024
Figure FDA0003013478100000025
Figure FDA0003013478100000026
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:
Figure FDA0003013478100000027
Figure FDA0003013478100000028
Figure FDA0003013478100000029
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:
Figure FDA00030134781000000210
Figure FDA00030134781000000211
Figure FDA00030134781000000212
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:
Figure FDA0003013478100000031
Figure FDA0003013478100000032
Figure FDA0003013478100000033
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:
Figure FDA0003013478100000034
Figure FDA0003013478100000035
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:
Figure FDA0003013478100000036
Figure FDA0003013478100000037
Figure FDA0003013478100000038
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:
Figure FDA0003013478100000041
Figure FDA0003013478100000042
Figure FDA0003013478100000043
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:
Figure FDA0003013478100000044
Figure FDA0003013478100000045
Figure FDA0003013478100000046
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:
Figure FDA0003013478100000047
Figure FDA0003013478100000048
Figure FDA0003013478100000049
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:
Figure FDA0003013478100000051
Figure FDA0003013478100000052
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 content of the first and second substances,
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.
CN201711497942.5A 2017-12-29 2017-12-29 Riverbed adjusting method based on underwater repose angle of sediment Active CN108595726B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
Sukumar et al. Brittle fracture in polycrystalline microstructures with the extended finite element method
CN105303612B (en) A kind of extract digital network method based on Triangulated irregular network model
CN109829217B (en) Fractured fracture oil reservoir productivity simulation method and device
Semenychev et al. Modelling and forecasting the trends of life cycle curves in the production of non-renewable resources
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
Jiang et al. Parallel contributing area calculation with granularity control on massive grid terrain datasets
CN110990926B (en) Urban surface building hydrodynamic simulation method based on area correction rate
CN103996221A (en) Plant organ mesh simplification method targeted for visualization calculation
CN104195979A (en) Riverway intersection water flow stagnant area wedge cone and building method and application thereof
CN108509532A (en) A kind of accumulation method and apparatus applied to map
CN110543691A (en) planar fracture grouting quasi-three-dimensional simulation method and system based on finite volume method
CN107564100B (en) Method for generating smooth entity by contour line
CN113688755B (en) Multi-flow-direction drainage basin feature extraction method based on hexagonal grid
CN110069843A (en) Based on the considerations of the Structural Topology Optimization Design method of ANSYS row wave effect
KR101335497B1 (en) System and method for fluid simulation using multigrid
CN111047687B (en) Three-dimensional T-spline-based heterogeneous material solid modeling method
CN111125805A (en) Method for identifying river section information of data-free area based on digital elevation model and river flow direction
CN114398746B (en) Construction method of equivalent drainage pipe network for earth surface overflow
Korobov Scaling properties of the area distribution functions and kinetic curves of dense plane discrete Poisson-Voronoi tessellations
CN111931263B (en) Constant flow simulation method based on optimized solution of planar two-dimensional shallow water equation
CN106875480A (en) A kind of method of city three-dimensional data tissue
CN113238295B (en) Method, device and equipment for analyzing original river channel slope of ancient river sediment

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