7 章問題解答 7- 予習. 長方形断面であるため, 断面積 と潤辺 S は, 水深, 水路幅 B を用い以下で表される B, S B + 径深 R の算定式に代入すると以下のようになる B R S B + ( B) 分母の /B は河幅が水深に対して十分に広ければ, 非常に小さな値となるため, 上式は R ( B) となり, 径深 R は水深 で近似できる. マニングの式の水深 を等流水深 0 と置き換えると, 以下の式となる Q B0 0 I e n 上式を 0 について整理すると, 以下の等流水深の式を得る nq 0 BI e 演習問題 7-- 式 7-を用いて収束計算により解を導く すなわち, 右辺に代入する水深を仮定値 とし, と求められた左辺 との差 ( ) が収束許容誤差 εに収まればよい ここで, 収束許容誤差を 0.00m とし, 初期の仮定値を.m として計算を行う 回目 ( a b) nq bi e + a b 0.0 00 0 00. 0 (. 0).66m 仮定値との差は,.66. 0.m となり, 収束条件を満たさない 求められた左辺 を新たな仮定値とすれば収束が進み, この作業を繰り返すことで以下のように解が得られる 回目 回目 回目 0.0 00 0 00 0.0 00 0 00 0.0 00 0 00 (.66 0).00m.66 0 (.00 0).00m.00 0 (.00 0).006m.00 0
回目で得られた水深 と仮定値との差は, 0.000m < 0.00m となり, 収束条件を満足し, 水深 は 回目の結果である.006m となる 実用的には, 非常に小さな桁は外すため.0m でよい 7-- 流水断面積 ( 式 7-7) および径深 R( 式 7-) を求めた後, マニングの式により流量 を求める D R ( θ inθ ) (.66 in 0 ).69m D inθ in 0 0.m.66 θ Q R I n.69 0. 0.0 00.6 m 演習問題 B 7--B 低水路と高水敷に分けてそれぞれの流量をマニングの式により算定する. まず, 低水路の流水断面積, 潤辺 S, 径深 R を求め, その後流量を求める 低水路の水面幅は m であることから, 形状要素は以下のようになる 0 7 7 09m S 00 + 7 + 7 0m 09 R S 0 低水路を流れる流量は, マニングの式を用いて, 以下のように求められる Q R I 09 9.09 m n 0.0 000 9.09m 次に高水敷は両岸が同じ形状であることから片側の形状要素は以下のようになる 高水敷を流れる流量は, 両岸を考慮した以下の式により求められる Q ( 6 + 0) 9 m S 0 + 6 + 6.7m 9 R S 6.7.0m R I 9.0 n 0.0 000 対象とする複断面河川を流れる全流量は, 以下のようになる Q Q + Q + 666m m 7--B はじめに, 中心角 αを, 式 7-9 を用いて算定する D θ co θ co D 0.7 co D θ co 7.96rad D D
次に, 流水断面積 および径深 R を求める D D inθ in 7 R 0.9m.96 θ ( θ inθ ) (.96 in 7 ).m マニングの式 ( 式 7-) を変形し, 水路勾配 I を求める nq I R 0.0 0.0. 0.9 7--B 台形断面における水理学的に有利な断面は, 正六角形の下半分の形状となることから, 断面積 および径深 R について式 7-, 式 7-6 の関係を用いて整理すると, 以下のようになる, R これを, マニングの式に当てはめると, 以下の式が得られる. 09 Q R I I n n 上式に既知量をあてはめ水深 を求める n I nq.09i 0.0 00.09 ( 000) 水面幅 B は以下のように算定される.m B b + a + 9.m 以上で回答は得られるが, 確認のため, 算出された水深をもとに流量を算定してみる Q R I I n n 99.m 約 00m / であることが分かり, 得られた水深は正しいといえる 7- 予習. dz : 河床勾配, dz d : 水面勾配, dz d d U + : エネルギー勾配 dx dx dx + + dx dx dx g. 等流水深 ( 式 7-) 0 I e nq
限界水深 ( 式 6-) c q g 上述の両式が釣り合うとき, 等流水深に含まれるエネルギー勾配が限界勾配 Ic となる nq I 上式を Ic について整理すると以下のようになる n g I c q 0 9 9 q g c 演習問題 7-- () 限界水深は以下のように計算される 限界水深 () 等流水深フルード数 c q g. 9. 0.6m qn 0 I. 0.0 /00 qn 0 I. 0.0 /0. 9.. F r.6 9. 0. F r 0...m 0.m () 水面形の概略図および曲線名は以下のようになる 0 c M ゲート M S S ゲート S 水路 /00 水路 /0 c 0 7-- () 等流水深 0, 0 および限界水深 c は, 以下のように算定される Q q 00.00m B 00 0 0 n q 0 I b 0.0 0 0 0 n q 0 I b 0.0 00 0.9 m.9 m
c q g 9. 0.97 m () 水路 Ⅰ の等流水深 0 を用いて共役水深 を求める q g 0.9 9. 0.9 F r.9 0.9 ( F ) (.9 ).m 0 r 上式の結果から > 0 であるため, 跳水は水路 Ⅱ 内で生じる 演習問題 B 7--B 式 7-の下付き に既知点である下流端断面の値を代入し, 断面 Bの初期仮定水深は下流端水深と同じ ( B.m) として計算を行ってみる. 収束確認.0. 0.9667 0. 000 i B > 収束条件に収まらないため, 仮定値を B i とし, 繰り返し計算を進める 何回かの繰り返し計算の後, 断面 Bの水深は.0m と算定される 同様に, 断面 C, 断面 Dの計算を行うと, 断面 Cは.0m, 断面 D は.7m となる 水深は徐々に減少していくが, 水位で考えれば, 下流端で H.m に対し, 断面 Bは H.00m, 断面 Cは H.0m, 断面 Dは H.97m となり, 徐々に上昇することが分かる 計算された値を用いて水位, 河床高の縦断図を描けば, より理解が深まる 7--B 式 7- の下付き に既知点である下流端断面の値を代入し, 断面 Bの初期仮定水深は下流端水深と同じ ( B.m) として計算を行ってみる i ( + Z Z ) B + Q Q g g + + B K K B. + -0.06 + 0.07.6 m X 収束確認.6. 0.6 0. 000 i B > 収束条件に収まらないため, 仮定値を B i とし, 繰り返し計算を進める 何回かの繰り返し計算の後, 断面 Bの水深は.m と算定される 同様に, 断面 C, 断面 Dの計算を行うと, 断面 Cは.7m, 断面 Dは.099m となる 水位は, 下流端でH.m に対し, 断面 BはH.m, 断面 CはH.97m, 断面 DはH.99m となり, 断面 Cまでは下流端からほとんど水位が変化せず, 断面 Dで大きく上昇することが分かる
7- 予習. 上流断面から流入する水の体積と, 下流断面から流出する水の体積の差から,dx 区間内に貯留 される水の体積を計算すると以下のようになる Q Q Qdt Qdt + dxdt dxdt また,dt 時間の体積変化量を, 断面積 を用いて表現すると以下のようになる dx + dtdx dx dtdx 上の 式は, 質量保存則により等しいため, 以下の関係が得られる Q dxdt dtdx 上式を整理すると, 以下に示す連続式が得られる Q + 0. 式の右辺を変換し, 連続式を代入すると以下のようになる QU Q + + U + Q + U Q + Q + U + + Q 以上より, 式が得られる n U + Q gi g g R b S 演習問題 7-- 連続式 ( 式 7-9) を差分形式に変換する B Q + 0 Q Q t B x ここに,Δt: 差分時間であり本問では 分間の事象を知りたいため 60 秒となる. 上式にそれぞれの値を代入すると, 水位変化量 Δ が求められる t Q Q B x 60 00 0 0 000 0. m したがって, 分間に 0.m 上昇したことが分かる 7-- 平均流速を求める Q 000 U 0.9 m 0. 6
シェジーの式による c/u の関係から, 伝播速度を求めると以下のようになる c U 0.9. m マニングの式による c/u の関係から, 伝播速度を求めると以下のようになる c U 0.9. m 演習問題 B 7--B 連続式 ( 式 7-9) を差分形式に変換した式を利用する B Q Qd Q + 0 t B x u 上式で Q u は上流からの流入量,Q d は下流への流出量,Δt は 600 秒,BΔx は貯水面積 0000m である 0 時 ~ 時までの水位変化量 Δ を算定する この 時間の流入 流出量の平均値は, 0 + 00 0 + 0 Q u m, Qd m となる 水位変化量は, 以下のように算定される t Qd Qu 600 B x 0000 時の水位は, 以下のようになる ( ) 0. m t t t H H + 60 + 0. 60. m 他の時刻の水位の計算結果は, 以下の表のようになる 時刻 (r) 0 6 7 9 0 流入量 (m /) 0 00 700 00 00 0 900 700 0 0 0 0 00 流出量 (m /) 0 0 60 90 00 70 000 0 60 0 00 0 0 水位変化量 (m) 0. 0..6 6.6.6 -. -. -. -. -.6 -.9 -.9 水位 (m) 60.00 60. 6.0 6.6 7. 7.6 7.0 70.6 6. 66.00 6. 6.9 60.00 計算結果をグラフで表すと, 下図 ( 右 ) のようになる 参考に, 流入, 流出の時間変化も示す m / m 600 0 00 流入量 7 00 水位 000 流出量 70 00 6 600 60 00 00 0 0 0 6 0 時刻 (r) 0 6 時刻 (r) 0 7--B 三角形断面の形状要素を現す式は, 台形断面の形状要素式 ( 式 7-0~ 式 7-) の底面幅を b0 とすればよい 7
水面幅 B b + a a 潤辺 S b + a a 流水断面積 ( b + a) a 径深 R S a a 任意断面に対する伝播速度 c と断面平均流速 U の関係を表す式 7-0に, 上記の各断面を表す式を代入すると,c/U は以下のように算定される ds S c B d U R S a d a ( a ) a a d a a ( a ) a a a a : ( )