下游常水位水力自动控制渠道运行动态过程与数学模型的研究.pdfVIP

下游常水位水力自动控制渠道运行动态过程及数学模型的研究 王长德管光华 张礼卫 (武汉大学) (深圳市水务局) 摘要：本文通过分析灌溉渠道中AVlO和AVIS下游常水位闸门的工作原理和运行过程，动用自动控 制理论求出两类闸门门后水位偏差和过闸流量之问的传递函数，并将之作离散化处理，再结合渠道 的明渠非恒定流方程联合求解，最后可得出设置AVIO、AVIS下游常水位闸门渠道的运行动态过程。 数学模型由计算机来求解，可为水力自动控制系统的渠道和闸门设计提供依据，并可预测系统的运 行过程与稳定性。 关键词：下游常水位 自动控制渠道非恒定流运行动态过程 自动化灌溉基本上可分为水力自动化和电气自动化两大类，随着计算机技术的发展， 电气自动化由于具有速度快、功能齐全、参数易于修改和优化等优点，得到了广泛发展。 但水力自动化亦凭自身优点，如设备简单、技术易普及和投资小、不要能源等，仍在世 界各国，特别在发展中国家得到了广泛应用。同时，水力自动化可作为电气自动化控制 的子系统，用于各下级渠道，特别是无电能供应的偏僻地区，作为对电气自动化系统的 补充和完善，也极有应用价值。水力自动控制系统，根据其所控制的是闸前水位还是闸 后水位，可分为上游常水位控制系统和下游常水位控制系统，相应所用的典型闸门为： AMIL门上游常水位闸门和AVIO、AVIS下游常水位闸门。 ． 由法国尼尔皮克公司研制的下游常水位水力自动闸门，简称AVI门，按其使用条件 口的渠道内，AVIS闸门用于门前水位变幅较小的渠道内。 1AVIO闸门的运行数学模型 如图1，这种闸门由弧形门叶、浮箱、浮箱套、均衡箱、臂杆、轴承等组成。它的 工作原理是：全关时，浮箱底与下游水面齐平，过闸门总重心的垂线经过轴心，在任一 开度，下游水位在控制水位时，重力矩与浮力矩平衡，闸门保持平衡。此时，重力矩为 开门力矩，浮力矩为关门力矩。 图1 AVIO门工作原理示意图 设在调节过程中，门后渠段需水量的改变为△Q0，门下流量的改变值为△Qi，由此 导致门后水位偏离控制水位，设其偏离值为Ah，则浮力变化值为： AM=j5il蟓△『l(R2一r2)，2， (1) 式中，马为浮箱，缘为水的比重，R为浮箱后半径，r为浮箱前半径。 由肼而产生的闸门运动角加速度： ￡=埘／‘，， (2) 式中，J为闸门的转动惯量。 闸门运动的角速度：to=jeedt。 (3) 闸门转角的改变值为： A≯={eAt (4) 在小变幅的情况下，闸孔高度的改变值6可近似的以下式表示： 6=RoA≯cOS谚o， 式中R为门叶半径，钆为闸门在初始位置的转角，此时门下流量的变化值为： △Qi=GuB2聃‰△≯瓜+瓦QoAh 式中，o-为淹没系数：p为流量系数；Qo为闸门在九开度时的门下流量：B：为相应于 这一开度的闸孔宽度；H。为门前水深。 令△Q=△Qi+△Q2，则： △Ql=咖2Rcos钆△驴瓜， (5) △Q=警 ㈦ 当浮箱在水中运动时，受到的水体阻力(以形状阻力为主)所产生的阻尼力矩为MF： MF=C。p∞2Bl(R4一，．4)／8 式中，e为形状阻力系数，无量纲，其值与箱体形状、姿态、表面上的光滑与粗糙度、 雷诺数的大小有关；p为水的密度；蜀为浮箱宽度。 因MF是国的非线性函数，而∞～般为较小值，在O～咀瞰内可将国2线性化，则上 式变为： MF=eCjp峨(尺4一r4)／8(7) c2为在O～q呱范围内国2线性化的系数，量纲为弧度／秒。如图