An Ag Engineering problem
Cam solved an engineering problem last week. He did a great job. I’ll offer another. This is a diagram of a diversion, which is a very simple and common earth structure for diverting water. The question is to find the flow area (the area below the blue line). If the depth is less than 0.5 foot on the cross slope, the velocity is slow and not counted, thus the vertical blue line. The diversion height is 3 foot (-1.5 to +1.5), but there’s a 0.5 freeboard so the water depth is 2.5 foot. The cross slope (black line) is in percent, 100 foot horizontal to 1 foot vertical. The red slopes are 6:1 (6 foot horizontal to 1 foot vertical). For this section we want the cut/fill ratio to be 1.5. The amount of cut (blow black) is 1.5 times the fill area (above black line). The top width of the diversion is 4 feet and the bottom width of the diversion channel is 12 feet.
Note: The +1.5 and -1.5 was read from the graph, and are approximate. As you can see in the drawing the top is slightly above 1.5 and the bottom of the channel is slight above the -1.5
The exact difference is 3.000 however. Here's a hint. Lines pass through the origin to simplify the process. Moving the whole construction up or down is how the exact cut/fill ratio is obtained.
Question: What is the flow area? This is simple geometry, if you want to go further, then what is the amount of flow in cubic feet per second using Manning’s Equation, given other inputs?
Diversion
Inputs
Diversion Solution
Note: The +1.5 and -1.5 was read from the graph, and are approximate. As you can see in the drawing the top is slightly above 1.5 and the bottom of the channel is slight above the -1.5
The exact difference is 3.000 however. Here's a hint. Lines pass through the origin to simplify the process. Moving the whole construction up or down is how the exact cut/fill ratio is obtained.
Question: What is the flow area? This is simple geometry, if you want to go further, then what is the amount of flow in cubic feet per second using Manning’s Equation, given other inputs?
Diversion
Inputs
Diversion Solution





13 Comments:
Does anyone have any thoughts about where to begin?
The area to be found is a complex geometric shape. Let's dissect it into simpler shapes, specifically, two trapezoids and two triangles.
dissected image
Triangle A (against the left bank) is a right triangle with a vertical side of 1.0 ft and a horizontal side to be determined. The slope of the bank reaching from the origin to 1.0 is given as 6:1, so the horizontal distance is 6 ft. Therefore the area of triangle A is
1/2 * 1 ft * 6 foot = 3 sq ft
Trapezoid B (traced by the bottom red line, ending at y=0 on the right) has a height of 1.5 ft. The bottom base is 12 ft; the top base is exactly 30 ft. (Not only does it look like 30 from the graph, it can be computed as 30: the red lines are 6:1, so they go 9 ft horiz in 1.5 ft vert; two of those, plus the 12 ft, = 9+9+12 = 30.) The area of trapezoid B is therefore
1/2 * 1.5 ft * (12 ft + 30 ft) = 31.5 sq ft.
The hardest measurement is the altitude of triangle C. I decided to figure it as the intersection of two lines.
The cross-slope line is easy; it starts at the origin and has a slope of 100:1 :
y = x/100
The red line starts at (20,-1.5) and has a slope of 6:1 :
(y+1.5) = (x-21)/6
y = (x-30)/6
These intersect at (1500/47, 15/47). So the altitude of triangle C is .3191 feet (15/47 exact).
Triangle C (from the zero line up to the black cross slope) has a base of 30 ft and an altitude that looks like 1/4 ft. Its area is therefore
1/2 * 0.3191 ft * 30 ft = 4.7865 sq ft
Trapezoid D (from the black cross slope up to the blue horizontal line) has a width of 50 ft (the black line rises half a foot and has a 100:1 slope), a left altitude of 1 ft and a right altitude of 0.5 ft. its area is therefore
1/2 * 50 * (1 + 0.5) = 37.5 sq ft.
The total flow area is therefore
3 + 31.5 + 4.7865 + 37.5 = 76.7865 sq ft.
I don't know Manning's Equation, so I shall forbear to go any farther.
This post has been removed by the author.
Ross,
You are off several square feet.
The diversion is exactly 3 foot high, but I see a minor problem. As you can see the top is not exactly 1.5 and the bottom is not exactly -1.5.
The adjustment up or down from the origin is how the exact cut/fill ratio is obtained. I did not realize I had say that correctly. That is not part of the given, but I used the +/- 1.5 to indicate where they were in the drawing.
I'll make a correction to the post, but without the adjustment you will not have the 1.5 Cut to Fill ratio.
Very good Ross, It was my error in the 1.5 that mislead you. I feel bad about that, so payback for my error,
I will tell you that adjustment up or down of the whole construction is the only way to get the cut/fill ratio right.
Ahhh ... the meaning of the term "cut/fill" ratio did not register, and so I treated it as a straight dissection problem, with straight edged shapes I could recognize.
Recasting this to come up with the cut/fill is going to take me a while... I'm not sure how to measure the cut and fill areas separately. I'll have to redo everything with coordinates.
There is probably an easier way and someone else knows it.
Ross,
I use x,y coordinates. Here's some history. I developed this solution back in the 1980's. In the '90s I programmed the solution into a PC and it used used all across Texas today.
I tried again, but my pencil failed me. So I set up a spreadsheet.
I did some searching and I found that given the coordinates of a convex quadrilateral there is a very simple formula for its area.
I set up my spreadsheet using one cell to represent the amount shifted up and down; set up cells that calculate the fill area and cut area, and their ratio; did a goal seek for the shift amount and found it to be 0.14172510 feet
The fill area is 42.1715 square feet
The cut area is 63.2572 square feet
The cut/fill ratio is 1.5000
The water area above the cut is 105.1364 square feet
Add this to the cut area to get the total flow area, 168.3900 square feet
I am very unsure of this solution, partly because the numbers are so odd and partly because it is so wildly different from my other solution.
My spreadsheet
My labelled diagram
Grrr. My diagram isn't available without logging into google, though I'm not sure why.
Here's a publicly available link:
labelled diagram
H1 is the height of diversion above coordinate (0,0)
SS1, SS2, SS3 are the 3 slopes (left to right)
cs is the cross slope in percent, divided by 100 to be ftft
(x1,y1) is left toe of diversion (far left)
I begin H1 = 1/2 the diversion height (dh), 1.5 in this example.
Then run thru the routine FindDivShape to find the area of the cut (AC)
and the area of the diversion (AD), this ratio gives the cut/fill ratio.
This initial run gives Ac= 36.2872, Ad= 22.0745, C/F= 1.6439
I then add a small number (sm) to H1 run the same equations again and compute
another CF. I use sm = 1 / 10000000000#
I use this numerical derivative to coumpute a new H1
SL = (H1 - Previous_H1) / (CF1 - CF2)
H1 = H1 + (CF - CF1) * SL
After about 3 or 4 adjustments to H1, the correct CF will be found
For this example H1 = 1.54460427161254
The rounded AC and AD are Ac= 34.7917 Ad= 23.1945
The area of Flow (AF) is shown in the Findq routine, Q is the quanity of
Flow in cubic feet per second.
The area of flow may include the AXTRA area (extra area for the out of bank flow)
This example has AXTRA = 12.63 sqft and area in the channel = 67.50 sqft for a total of 80.13 sqft
Using the wetted perimeter and the hyraulic radis gives an Average Velocity of 1.33 ft/sec
The total flow would be 106.35
FindDivShape:
If cs = 0 Then
y1 = 0: x1 = H1 * SS1 + CR + H1 * SS2
Else
y1 = (H1 * (SS1 + SS2) + CR) / (1 / cs - SS1)
x1 = y1 / cs
End If
FindAreaDiv:
AD = (x1 * H1 + (SS2 * H1 + CR) * y1 + CR * H1) / 2
h2 = DH - H1
If cs = 0 Then
y1 = 0: x2 = h2 * SS3 + BW + h2 * SS2
Else
y2 = (h2 * (SS2 + SS3) + BW) / (1 / cs - SS3)
x2 = y2 / cs
End If
FindAreaCut:
AC = (x2 * h2 + (SS2 * h2 + BW) * y2 + BW * h2) / 2
CF1 = AC / AD
Return
Findq:
dist = h2 * SS2 + BW / 2
DCL = h2 + dist * cs
'Compute Wetted Perimeter on Channel
wp = d2 * sf2 + d2 * sf3 + BW
TW = d2 * SS2 + d2 * SS3 + BW
bstor = AF
AXTRA = 0
wpxtra = 0
w7 = d2 - h2 - y2
If w7 > 0 Then
bstor = AF + (w7 / cs) * (w7 / 2) - ((w7 * SS3 * w7) / 2)
TW = (H1 - FB) * SS2 + (H1 - FB) / cs
End If
If w7 >= (0.5) Then
'Compute extra flow out of bank
x3 = (w7 - 0.5) / cs
AXTRA = ((w7 - 0.5) / cs) * ((w7 + 0.5) / 2) - ((w7 * SS3 * w7) / 2)
If AXTRA >= 0 Then
AF = AF + AXTRA
wpxtra = (w7 - 0.5) / cs - (w7 * sf3)
wp = wp + wpxtra
Else
AXTRA = 0
End If
End If
HR = (AF) / (wp)
QF = 1.486 * HR ^ (2 / 3) * AF
VEL = 1.486 * HR ^ (2 / 3) * Sqr(csl) / MN
Q = AF * VEL
Return
Plot Coordinates
x y
-23.97 -0.24 1 Back toe
-13.27 1.54 2 back crown
-9.27 1.54 3 front crown
0.00 0.00 4 cut to fill intersection
8.73 -1.46 5 bottom width begins
20.73 -1.46 6 bottom width ends
31.35 0.31 7slope intersects grd
54.46 0.54 8 flow bottom edge
54.46 1.04 9 flow top edge
-6.27 1.04 10 flow top intersects diversion
-11.27 1.54 11 top center diversion
-11.27 -0.11 12 bottom center diversion
14.73 0.15 13 top center of channel
14.73 -1.46 14 bottom center bottom width
-23.97 -0.24 15 back existing ground
104.46 1.04 16 existing ground @ water line
Ross, you did a lot of hard work, many of your coordinates are not far off, but the area are off quite a bit.
You went out beyond 60 feet on the right, so something there got off a bit. Maybe it was the height (above 0,0) was only 1.54, and I think you had 1.6+.
One of these days I am going to figure out the algebraic derivative, rather than the numerical derivative, but I don't expect it will change much.
Post a Comment
Links to this post:
Create a Link
<< Home