《SMS81地表水模拟软件使用说明与开发》由会员分享，可在线阅读，更多相关《SMS81地表水模拟软件使用说明与开发（8页珍藏版）》请在金锄头文库上搜索。
1、SMS8.1 地表水模拟软件使用说明与开发SMS 是英文 Surface Water Modeling System( 地表水模拟系统)首个字母的组合，是美国 陆军工程兵水利工程实验室(United States Army Corps of Engineers Hydraulics Laboratory)和 扬伯明翰大学(Brigham Young University)等合作开发的商业软件，可用来模拟水体的流场 和浓度场。它由FESWMS-2DH、RMA2、RMA4、SED2D等软件包组成。本章对将使用的 RMA2和RMA4软件包及其强大的后处理功能作较详细地说明。1 水动力模型的建立水动力
2、模型建立的步骤如图 1所示：图 1 RMA2 二维水动力模型的工作流程1.1 输入底图采用SMS进行地表水模拟时，首先要输入较详细的底图。SMS8.1软件输入底图的途径 主要有两种：一、打开 tiff 格式的图形文件，将地图中随机的三点的坐标输入定位，使底图与实际地 形吻合；二、此外SMS8.1软件还支持AutoCAD转换成的R12的.dxf格式的电子地图，可以直接 打开.dxf格式的地图。1.2 定义节点具体操作如下：1、 点击图标增加节点，利用图标飞选择节点，以便对节点进行修改；2、选中两节点，在 SMS8.1 版本中在窗口下方偏右出现 distance 项为两节点距离，如 图2所示。已知
3、两点距离，可在Node/lnterpolation Option项中调整插入的节点数, 根据实际情况插入节点。3、 节点坐标包括x、y、z值可在坐标框中输入，确定节点的确切位置及高程，如图3 所示。Node info: 2 selected; First id = 8495; Second id = 8968; Distance = 3903.132791 m;图2两节点间距离显示图X|29750. 1001Y|41325.20635z|30.7图3 节点坐标显示图1.3 构建网格1.3.1 网格的构建应满足以下条件：1、据水力特征流速大小和过水能力的大小使网格疏密有致；2、构建模型者所估计的
4、流线平行；3、三角形与四边形注意过渡。边界、流场复杂区域用三角形，流速均匀、航道、湖区 等地使用四边形。1.3.2 构建方法：选中四个节点，点1、手动添加网格单元：选中三个节点，点击，形成形如图标的六节点三角形网格;形成八节点的四边形网格；点分割四边形网格或将两相邻的三角形网格合成为四边形网格；点将内含两个三角形网格的四边形网格调转方向，使得网格满足适应流线，避免出现网格质量的原则。节点如图 4 所示，2、通过节点线构造三角形、四边形网格：点连接节点线，三角形网格三边界，四边形网格四边界，从起点开始，终点双击结束，形如图 5 所示；最终，6 条边界的节点线如图6 所示。图5图6选择节点线，三边
5、界须构造三角形网格单元，选中三条边界，在 Element/TriangularPatch构造三角形单兀，点Preview预览，如图7；Cancel图7四边形网格构造与三角形网格构造大体相同，选其四条边界，在 Element/RectangularPatch中构造，如图8；另外可根据图中的四条边界End1、End2、Sidel、Side2中Distribute命令可在中间插入若干个节点，此方法适合构建大规模、结构简单的四边形网格。图83、通过 map 模型，构造大规模的网格按畔切换至地图模型，点匸在底图边界上绘制节点线，起始点单击，在终点双击结 束此条节点线的绘制,如图 9。绘制另一条相邻的节点
6、线时，软件将自动识别节点之间的耦合性如图 10。不管地形规则或不规则，都要形成一个图块，图块至少需要四条节点线，形如图 11 所示。构建模型时，为了不同的材质、地形可设置互相连接的若干个块；连接完所 有的块，点feature objects/build polygons，这时，模型自动建立块，但是从界面上看不出来。 点因再选择你所要构造的块，出现如图12所示的效果。接下来点feature objectiveattribute，出现如图13所示，您可以在Mesh Type中选择网格的类型，分paving、 patch、adaptive tessellation 分别代表三角形网格、四边形网格、密
7、集型三角形网格。点 preview 预览。假如构造的块含有4条以上的节点线，点斥选中某个你认为不必要的点（此点是节点线的起点或终点）在node option中选择merge就可以把此点转换成中间点，可进行下一 步骤。最后点 feature objects/map-2D mesh 把地图生成二维网格。图9图 10图 11图 12网格构造完后，要对模型定义底高程。确定模型所用的单位（米或英尺及其它），在 editcurrent coordination 中改变 unit 选项调整单位。图 13 网格的类型1.4 模型质量检查整个计算区域的网格划分完毕后，为使网格收敛，点击图标选中Mesh qual
8、ity选项, 对网格进行检查，检查内容包括：1、三角形内角满足10Wa 150，四边形内角满足30Wa 150；2、两单元之间的地形最大坡度10%；3、相邻有限元的面积比是否控制在 0.52 之间；4、三角形、四边形任意两条边长度的最大比值不超过 10；5、节点地面高程大于水面高程时，是否打开“干/湿”边界条件；6、为保证能量守恒，尽量使系统的边界光滑；7、三角形、四边形有限元是否光滑、平整；8、为节省计算时间，需要对有限元网格重新编号。对于不符合要求的有限元网格如图 14 所示，可根据检查结果进行调整，具体步骤如下1、首先将选项 Nodes/Locked 前面的对号去掉，以便可以任意移动节点
9、；点击图标*豪移动节点，直到该节点处不出现质量问题为止；2、对于一些节点需改变其底高程，可在图 3 节点坐标显示图中输入适当的值，直到不 出现网格质量；3、点击图标一按逆时针方向，点击需光滑的节点线起始点，按住Ctrl键点击终点，然后点击图标二J选中需光滑的节点线，点击Nodestrings/Smooth命令即可；4、此外，对于四边形单元应尽量使其长边与短边的比例增大，并使长边方向与水流方 向一致。在敏感区域（如水流变化较大处）需适当增加单元数，以使网格收敛。修改后直至不再出现网格质量问题，方可进行下一步操作。图 14 网格质量显示图1.5 定义材质不同的材质有不同的糙率，定义材质步骤如下：1
10、、添加新材质，选择 Edit/Material Date 点击 New 创建新材质，点击对话框中的 IDName 和 Pattern 可以定义不同的特性。2、选择Element/Option,点击Set default material可以设置默认的材质。若需更改单兀的材质，可以点击工具选择需要修改的单元。然后选择Element/Assign Element Type重新选择材质。3、最后在RMA2/Material Properties中设置材质的不同的参数（曼宁系数和涡粘性系数等）。1.6 设置边界条件1、利用1/%工具，按下Shift键点击起始点，最后双击终点即可创建节点串。然后利用工具
11、点击节点串，查看边界线方向。入流边界方向应朝向计算域内，出流边界方向应朝向计算域外。若不符合要求，可以选择 Nodestrings/Reverse Direction 更改方向。 置入流和出流边界条件。选项Boundary Condition Type中可以设置边界的类型Specified flow 可以定义流量边界；Water surface elevation可以定义水位高程边界；Reflecting boundary可2、入流与出流边界方向确定后，利用工具点击节点串选择RMA2/Assign BC分别设以定义反射边界； Rating curve 可以定义非恒定边界。一般进口处可以设置入流
12、量和流速 出口处设置水位。1.7 模型控制参数设置模型控制参数的设置是用来定义水动力模型的初始条件的。选择 RMA2/Model Control 打开模型控制对话框，参数设置具体步骤如下：1、General 选项：计算机类型 Machine Type 选择第一项； Water Properties 选择中， Temperature 设置水温，Density 设置水的密度，Specify initial water surface for 定义模型出口 初始水位，Specify ID node initial condition设置初始水深和流速估计值；其它选项取默认值, 如图 15 所示。图1
13、5 RMA2模型控制General选项对话框2、Timing选项：模拟类型Simulation Type中恒定计算选择Steady state，动态模拟则选 择 Dynamic；设置迭代次数 Iteration For Flow Calculations/Initial；在 Depth Convergence Parameters/Steady state depth选项中输入模拟计算的收敛标准（一般取0.0001）；非恒定计算 选择Dynamic，迭代次数和收敛标准设置同上，并在Computation Time中输入各项值。3、File 选项：该选项中设计到文件的存储功能。若是冷启动，则保持
14、选项中的默认值 即可；若为热启动，则需在进行初次计算时打开 Write hotstart 按钮，计算完毕后，系统会自 动生成“*.hot”文件，若需要利用同一地形继续进行计算，则可以打开Hotstart input按钮， 找到“ *.hot ”文件，便可进行热启动计算。4、Materials 选项：该选项主要控制干/湿地的打开和关闭。5、Weather 选项：该选项可以控制风效应、降雨、蒸发等因素对计算的影响。1.8 RMA2 模型检查上述操作完成后，便可进行RMA2模型检查了。选择RMA2/Model Check,然后点击Run Check即可运行模型检查。若出现错误，可根据提示进行修改；若
15、无错误提示，剩下的可能 提示就是对网格进行重新编号。此时可以利用个工具点击出口处节点串，选择 Nodestrings/Renumber 命令即可对网格重新编号，文件存档。1.9 模型运行对网格进行重新编号后，便可进行模型计算。选择RMA2/Run RMA2命令进行模型计 算，若模型计算结果收敛，则计算框内会出现“RMA2 finished”字样；若模型计算发散， 则会出现“Stop depth convergence exceed 25.0”字样；此时必须重新对网格进行检查，可以 通过修改网格、修改参数设置等措施使网格最终收敛。具体操作如下：1、打开网格质量显示 Mesh quality ，查看网格是否还有未发现的质量问题；2、网格总体方向应与水流流向大体保持一致，特别是过渡地区，三角形与四边形应合 理选取，网格疏密程度应与水流特性相符，水流较复杂、流速变化较大的区域应加密网格， 反之则适当减小网格密度；3、检查边界条件的设置：出口水位是否太低，入口流量是否太大，在模型控制中，初 始水位的设置可以进行调整，方便模型收敛，入口与出口边界应设置过渡段防止局部区域流 速过高；4
