《Autobank计算原理和结果》由会员分享，可在线阅读，更多相关《Autobank计算原理和结果（7页珍藏版）》请在金锄头文库上搜索。
1、渗流计算原理对于稳定渗流，符合达西定律的非均各向异性二维渗流场，水头势函数满足微分方程斗 k 空+ ) k 纠 + Q 二 0(1-1)dx I x dx 丿 dy ( y dy 丿式中：甲二cp (x,y)为待求水头势函数；x，y 为平面坐标；K ，K 为 x，y 轴方向的渗透系数。xy水头p还必须满足一定的边界条件，经常出现以下几种边界条件：1)在上游边界上水头已知p =p nn（ 1-2 ）2)在逸出边界水头和位置高程相等p =z（ 1-3 ）3)在某边界上渗流量 q 已知k 竺 l + k 竺 l =-q(1-4)x dx x y dy y其中l，l为边界表面向外法线在x, y方向的余
2、弦。xy将渗流场用有限元离散，假定单元渗流场的水头函数势p为多项式，由微分方程及边界条件确定问题的变分形式，可导得出线性方程组：HE = F(1-5)式中H渗透矩阵；p 渗流场水头；F节点渗流量。求解以上方程组可以得到节点水头，据此求得单元的水力坡降，流速等物 理量。求解渗流场的关键是确定浸润线位置， Autobank 采用节点流量平衡法通 过迭代计算自动确定浸润线位置和渗流量。Autobank渗流计算报告1 概要分析类型:二维稳定渗流 工况数量=4工况0,最高水位=1950.4m 工况1,最高水位=1959m 工况2,最高水位=1960.71m 工况3,最高水位=1961.32m2 材料参数
3、材料参数表材料名Kx(m/s)Ky(m/s)给水度坝土1.33e-061.33e-060.01坝基2.301e-062.301e-060.01棱体9.259e-069.259e-060.01灌浆5e-075e-070.013 渗流量流量表工况0H=1950.4(m)流量=4.38 52 4e-0 6(m3/s)0.378885(m3/d)工况1H=1959(m)流量=1.20363e-05(m3/s)1.03994(m3/d)工况2H=1960.71(m)流量=1.40105e-05(m3/s)1.2105(m3/d)工况3H=1961.32(m)流量=1.408 01e-05(m3/s)1.
4、21652(m3/d)4 附图工况3,H=1961.32(m)，水头(m)等值线计算模型有限元网格工况0工况0,H=1950.4(m),流速矢量图工况0,H=1950.4(m),水力坡降等值线工况1工况l,H=1959(m),流速矢量图17.814-9-工况l,H=1959(m),水头(m)等值线水位,990.17=死水位，1950.40m20.823.8 26.7-29. 7-32.7工况1,H=1959(m),水压(m)等值线O01工况1,H=1959(m),水力坡降等值线工况2十设核洪水位位，19961.3ac正常水 J位，1959m死水位，1950.40m-/ :匸 - /； P.-
5、- - _ - 1-:7 ? : : : Z Z J : : . _工况2,H=1960.71(m),流速矢量图工况2,H=1960.71(m),水头(m)等值线-死水位,1950.40m0 00031321.92528.3134.4工况2,H=1960.71(m)，水压(m)等值线0 362。 谦0；=死水位，1950.40m正s &0.362009060工况2,H=1960.71(m)，水力坡降等值线工况3工况3,H=1961.32(m)，流速矢量图工况3,H=1961.32(m)，水头(m)等值线hm19.10. 00031工况3,H=1961. 32(m)，水压(m)等值线22.3 -25.5-28.631.8 - I -.0工况3,H=1961.32(m)，水力坡降等值线
