河道平面二维水沙数学模型的有限元方法2行政.docVIP

河道平面二维水沙数学模型的有限元方法2行政论文范文大全 河道平面二维水沙数学模型的有限元方法 3　数值解法 　　为书写方便，采用通用变量p来代替方程(13)～(17)中的变量zs、u、v、s和n，用fp表示方程中等号右边项，并用对角矩阵 来代替aij，则方程可化成统一的形式： dpi/dt=fpi。 对该方程的求解，模型采用了“预报-校正-迭代［5］”的计算方法。其方法为：用二阶显式adams公式作为预测公式，梯形公式(隐式)作校正公式，可构造adams二阶pc公式。离散后的方程不需联接，计算过程稳定性好，时间步长可取得较长，同时还有效避免了数 值振荡。具体计算过程如下 预报　　　　　　　　　　　　　　　　(18) 校正　　　　　　　　　　　　　(19) 　迭代给定误差洌杂谒械1≤i≤np，若　　　，则令： ，否则令：= ；转到校正，继续迭代，直到满足精度要求为止。当求出hi、ui、vi、si和ni后，代入河床变形方程便可求得冲淤变形后的河床高程zn+1bi，这样便完成了一个时段的计算。以上角标n表时段，*表预报值，**表校正值，无*标记的表示该时段的计算结果。 在水位求解过程中还采用了kawahara“选择系数集中［4］”的经验处理方法，引入了选择性集中系数： α=0.8～1.0， 其作用相当于对计算结果进行“光滑”处理。 由于系数矩阵为np×np阶矩阵，若直接进行存储会占用大量的内存，从而导致计算无法进行。笔者发现系数矩阵中绝大部分为零元素，为此提出了压缩存储的方法，即只对系数矩阵中的非零元素进行存储和计算，同时用辅助数组指明相应非零元素所在的位置。 　　对于长系列水沙条件的计算，提出了“非恒定-定非-定流”的算法，即将整个非恒定流量过程概化为梯级恒定流，对每个梯级恒定流采用非恒定流方程，以时间步长为迭代参数进行迭代，直到得到恒定的流场。 3　数模计算有关问题处理 4.1　初始条件及边界条件 初始条件：可由初始刻实测资料给出。计算所需的初始条件在实际计算中一般难以全部获得，只能通过估算加以补足，当然估算得越接近实际越好。 边界条件：进口给定流量、含沙量及其级配，出口给定水位；对于不滑动岸边界，取u=0,v= 0；对于滑动岸边界，取边界法线方向流速分量为0。 对于活动边界，采用了动边界模拟技术：对每个膖时段，采用计算的水位及水深值判别和区分水域和陆域计算节点。对陆域计算节点，使其保持一较小富余水深(hmin=0.001m)，并取其糙率为一个接近于无穷大(如1010)的正数。 4.2　三角形网格划分 基于网格生成的基本思想，提出了四边形法。对于单一河道，首先根据河势将河道剖分为若干大的四边形，再分别将这些大的四边形剖分为若干小的四边形，最后将四边形的较短对角线相连，即形成三角形网格。对于分汊河道，将每条汊道当作单一河道处理后再进行拼接，便可形成整个计算域的网格。根据该方法，我们编制了通用的计算程序，只需输入少量的信息，计算机便能自动生成网格，并给出节点坐标和单元关联信息，最后配以屏幕显示及图形自动绘制。 4.3　紊动粘性系数 根据零方程紊流模型，紊动粘性系数由vtαu*h公式确定。其中u*为摩阻流速；α为常数，经调试确定。对于泥沙扩散系数εs ,可近似取ε=vt。 4.3　糙率 通过实测资料反求，并根据河道中不同的部位分块调试糙率。 4.4　床沙级配计算方法 在泥沙冲淤频繁的河段，由于水流与泥沙的相互作用，使某组泥沙发生冲刷时，另一组泥沙可能发生淤积。因此，床沙级配的计算应能反映泥沙冲淤交替过程。本模型将计算河段内河床泥沙概化成三层，即表层(交换层)、次表层(扰动层)和深层。当表层某组泥沙发生冲刷时，处在次表层的该组泥沙将受到扰动，根据扰动强度的大小，确定该层泥沙进入表层参加交换的量。 5　计算实例 河段概况：监利河弯段位于下荆江，上起塔市驿，下止于沙夹边，长约20km。该河段平滩河宽约1400m，最宽处约3200m(乌龟洲)，最窄处约1000m，平滩水位下的平均水深约11m。河段内有一高程为30m(黄海)左右的乌龟洲(洲长约7km，宽约2km)将河道分成左右两汊。本河段内有姚圻脑水文站，该站多年平均流量为11.370m3／s，多年平均含沙量为1.121kg/m3。河床主要由粉质粘土、砂粘土和细砂组成，河床表层床沙中值粒径一般在0.00～0.22mm之间。 图1 计算河段河势及网格布置图 fig.1 the calculated river reach and the grids 计算区域的选取及网格划分：选取塔市驿至沙夹边长约20km的河段作为计算河段。采用三角形网格，将整个计算河段划分为100个断面，6732个三角形网格单元，共有3500个节点，见图1。 地形条件：以1993年10月实测河道地