mike flood一二维耦合计算溃堤.pdfVIP

MIKE FLOOD 一二维耦合计算溃堤 笔记本 ： MI KE 创建 时间 ： 2015/9/23 10:05 更新 时 间 ： 2015/9/23 11:10 作者 ： wushan qin lian @126.com 一二维耦合计算 ，分别用一二维建模 ，然后连接。 一维模型可先按照mike11建模 ，先不考虑溃口对一维模型的影响 ，将模型调试通过 ，率定等。 二维模型尽量不设置开边界 ，溃 口流量通过两个模型连接传递。二维模型优先选用矩形网格。 **************************************************************************************** 一维模型 一维模型建成后 ，在network文件中设置溃口位置和溃口基本信息。 由图可知 ，溃坝采用的是溃坝模型 ，然后通过侧连接 （side structure ）与洪泛区连接。Dam Geometry 设置溃堤堤防的顶高程和溃口宽度 ，breach calc.Method 设置溃堤的计算方法 ，有能量方程法和NWS 溃堤算法。Failure Moment and Mode栏 ，主要设置溃堤控制条件——某时刻溃堤或者某水位溃堤 ， Failure Mode设置采用时间序列控制溃堤过程还是根据堤防冲刷能力控制 ，一般选用时间序列控制。 一维模型里 ，还需要在Bondary文件里设置溃堤接口 设置建筑物边界 ，溃坝模式 ，Branch name的位置 ，需要填SS_河道名称_溃 口累加距 （SS应该是side structure的缩写 ，累加距需要保留一位小数 ，并加 M ），Chai a e一般给50 ，GateID要与network里 的一致。 Dam break level、width、slope分别代表溃堤的高程 ，宽度和溃 口两侧的坡比。可以把这三个要素做到 一个dfs0文件里 ，格式如下 溃堤水位 ，一般比最高水位低一些 （5cm或者其他 ）。增加了溃口以后 ，即便溃口不溃 ，也会影响局部 水面线——有点类似在两个断面中又加了一个断面。所以 ，加了溃 口以后 ，可以在network设置一个较高 水位保证不溃 ，对水位进行以下率定。然后提取溃口位置最高水位 ，修改network里溃堤水位。 选择合适的时间步长 （要协调二维计算 ） 至此 ，一维模型部分就完成了。 **************************************************************************************** 二维模型 （矩形网格 ） 二维模型建模与mike21类似 ，地形网格和糙率网格建好以后 ，新建m21文件 ，导入地形和糙率。 时间步长设置 ，步长应与一维的时间步长一致。 Boundary一般不设置边界。 干湿边界可以调小 ，保证模型运算稳定。 涡粘系数项选用Sma orinsky Formula ，取值0.1——何哥经验总结。 其他设置 ，与mike21无差别。 至此二维模型部分也设置完成 *************************************************************************************** 连接一二维模型 新建一个couple文件 在Linka e Files中载入做好的一二维模型文件 在Link Definitions中设置连接位置 River Chaina e就是溃口在一维模型中的累加距。M21 Coord1中是二维模型的单元格 ，J和K代表行列 数。 Side Structure Link Options中设置Depth Adj ust 为Yes ，在Ex ponential smoothin factor处设置一 个接近1的数字。Ex ponential smoothin factor表示建筑物上下游水位的平滑 ，越趋近于1越不平滑 ， 越接近实际计算 ，但对模型稳定没有好处。 计算点击菜单栏的run ，选择Start simulation就可以了。 *************************************************************************************** 调试模型 计算过程中可能会各种报错 ，也可能计算到一半就跳出了。 出现这种情况 ，可能有以下原因 ： 1、一维模型本身稳定性差