【正文】 UAL39。年 39。年均流量 (萬噸 )39。年 39。深泓 (m)39。距壩里程 (m)39。年 39。高程 (m)39。水位 (m)39。INTO FZV39。庫容 (億立方米 )39。年 39。起點(diǎn)距 (m)39。水面線 (m)39。年 39。距壩里程 (m)39?？傆俜e量 (億立方米 )39。懸移質(zhì)輸沙量 (萬噸 )39。INTO FILE39。MASS NOT EQUAL39。INTO DEPO39。INTO SGB39。IN AREA39。AND39。IN AREA THE WATER LEVEL IS TOO HIGH OR THE HEIGHT OF THE SECTION AT EDGE IS TOO LOW39。INTO AREA39。nr,ip,zmin,fmin39。INTO LEVEL39。讀入天數(shù) :39。逐年水位庫容關(guān)系變化 .txt39。逐年壩前斷面淤積變化 .txt39。逐年水面線變化 .txt39。逐年深泓變化 .txt39。逐年累計淤積量 .txt39。流量沙量 .txt39。河床地形 .txt39。 上述計算結(jié)果 與現(xiàn)實不太 相符 ，其原因是在計算該水庫庫容的泥沙損失時，計算沖往下游的泥沙 量不合理。 B、壩前斷面變化曲線分析： 由壩前斷面變化曲線圖可以看出，隨著時間的增加，壩前斷面各點(diǎn)的總體高程也成增高趨勢。（其中推移質(zhì)輸沙量是按照Y Y Y Dy(i)=0,重新計算 S(j) Y 一維泥沙淤積計算 ~ 5 ~ 懸疑質(zhì)輸沙量的 %計算的）。0 1*10 )( ??? tSSy ???? 公式補(bǔ)充 mgRukS ????????? ?3* K取 ， m取 ， 干密度 39。 二、 基本原理 基本方程 水流連續(xù)方程： 0?????? xQtA ① 水流運(yùn)動方程 ? ?fiigAxhgAAQxtQ ?????????????????? 02 ② 或 034222 ?????????????????? RA nQgxzgAAQxtQ ③ 泥沙連續(xù)方程 ? ? ? ? )( *SSQSxSAt ???????? ?? ④ 河床變形方程 )( *00 SSxGty b ???????? ??? ⑤ 推移質(zhì)平衡輸沙方程 G=G* ⑥ 水流挾沙力公式采用張瑞瑾公式，推移質(zhì)輸沙率公式采用 Mayer_Peter 公式，MAYERPETER 公式中的能坡 J 按均勻流曼寧公式近似計算（每個斷面不同 ）。一維泥沙淤積計算 水庫一維泥沙 淤積計算課程設(shè)計 一維泥沙淤積計算 目錄 一、目的與要求 ..................................................................................................... 1 二、基本原理 ......................................................................................................... 1 基本方程 ...................................................................................................... 1 方程離散 ...................................................................................................... 1 公式補(bǔ)充 ...................................................................................................... 2 三、計算步驟 ......................................................................................................... 2 四、計算框圖 ......................................................................................................... 3 五、計算結(jié)果 ......................................................................................................... 4 歷年輸沙量特征值 ........................................................................................ 4 各年淤積總量 ............................................................................................... 5 各年水位庫容關(guān)系 ........................................................................................ 6 水面線的變化 ............................................................................................... 7 深泓變化 ...................................................................................................... 8 壩前斷面變化 ............................................................................................... 9 六、結(jié)果分析 ....................................................................................................... 12 剖面形態(tài)分析 ............................................................................................. 12 庫容損失合理性分析 ................................................................................... 12 七、計算程序 .......................................................................................................... 13 一維泥沙淤積計算 ~ 1 ~ 一、 目的與要求 通過課程設(shè)計，初步掌握一維數(shù)學(xué)模型建立數(shù)學(xué)模型的 基本過程和計算方法，具備一定的解決實際問題的能力。 泥沙條件：包括懸移質(zhì)，推移質(zhì)的均勻沙模型，推移質(zhì)計算模式為飽和輸沙，懸移質(zhì)計算模式為不飽和輸沙，水流泥沙方程采用非耦合解。0 ????????? ? ?????????? ??t yByBx QSQSx BGBG jjjjbb ? 對于進(jìn)口斷面，推移質(zhì)不考慮，懸移質(zhì)采用單點(diǎn)離散 方程（ 5）可離散為： 39。0 1*10 )( ??? tSSy ???? 修改水各斷面水下河床高程 YHH ??? 0 重新進(jìn)入（ 2）進(jìn)行下一循環(huán) 計算 10年河床變形，計算時段為一天，單位為秒（ s） 1淤積總量年輸出一次，其余每兩年輸出一次計算結(jié)果 四、 計算框圖 開始 一維泥沙淤積計算 ~ 4 ~ 讀入地形資料和糙率 輸出初始庫容、深泓等 i=i+1 讀入一個 Q Q=0 計算水流子程序 計算 S*和 Gb j=j+1 計算 S(j)和 dy(j) 判斷庫尾沖刷與否 Jnpnxt 修改河床地形 判斷年份輸出結(jié)果 I3653 結(jié)束 五、 計算結(jié)果 歷年輸沙量特征值 計算各年懸疑質(zhì)輸沙量、推移質(zhì)輸沙量、年均流量見表 1。該水庫運(yùn)行 10年后，在 距壩約 10000米處，深泓點(diǎn)高程增加 40多米， 計算結(jié)果 表明 ， 該 水庫泥沙 以淤積為主，致使水庫庫容逐年損失，且情況較為嚴(yán)重。 結(jié)果表明，該水庫的調(diào)節(jié)能力較差。 一維泥沙淤積計算 ~ 13 ~ 七、 計算程序 program main parameter(nn=16,mm=60,ndisp=0,npxt=16,ND=3653,NY=10,T=15