matlab计算系统过渡过程时间,MATLAB在水电站水力过渡过程计算仿真中的应用

一、引言水电站水力过渡过程是水、机、电系统相互影响、相互制约的复杂过渡过程,也是水电站实际工作过程中不可避免的一种特殊运行工况。特别是水轮机甩负荷计算是水电站和水力机组设计不可缺少的一项工作。虽然采用数解法和图解法可以进行计算,但精度已不能满足现代水电站的技术要求。随着计算机应用的推广,在大中型水电站和压力引水系统比较复杂的电站设计中越来越多地采用电子计算机进行水力过渡过程计算仿真。MATLAB是集数值计算、符号运算及图形处理等强大功能于一体的科学计算语言,作为一种大型数学软件,它包含了丰富的函数资源,使人们解决各种学科的数学问题的编程简易可行。结合在系统研究中的应用实践,以下将介绍MATLAB在水电站水力过渡过程计算仿真中的应用。所有方针结果均在MATLAB6.5环境下实现。二、数学模型建立水电站水力过渡过程过渡过程计算仿真涉及到:管道的瞬变流计算;上、下水库端边界计算;水轮机边界计算;系统初始工况计算等。1.单一特性管的瞬变计算在文献[1]中,对管道系统瞬变流的基本方程有详细具体的推导,在此不赘述,只列出其最终的适用公式,即正、负特征方程:QP=CP?CaHP(1)QP=Cn+CaHP(2)式中,Ca=gA/a;Hp、Qp分别为t时刻管段第i节点处的压头和流量;Cp、Cn分别与t-t时刻的压头和流量有关,对t时刻是已知量。其中,g为重力加速度,A为管道截面面积,a为水击波速。2.上、下水库端边界条件水库水位按恒定值考虑,若进、出口损失和流速水头均被忽略,则有水头方程HP=Hres(3)式中,Hres为水库水面在基准面以上的高度。3.串联管联结处数学模型对于串联管联结处的计算数学模型,有的文献虽然考虑了联结处流速水头的差别和局部阻力损失,但是未考虑水流方向的影响[1],或者局部阻力系数水流方向不同的区别[2]。为此,作者推导出了考虑水流方向影响及局部阻力系数区别的总水头方程,得到下列数学模型。(1)正、负特征方程Qpi,n+1=Cpi?CaiHpi,n+1(4)Qpi+1,1=Cni+1+Cai+1Hpi+1,1(5)(2)连续方程,+1+1,1=pinpiQQ(6)(3)总水头方程2121,12,12(1)2,11,11,1+=++iPPpiiPpingAQQHQgAHkinII(7)上式中,第一个下标指管段编号,第二个下标指管段中计算节点编号。如Hpi,n+1表示第i段管第n+1号节点的压头;Qpi,n+1表示第i段管第n+1号节点的流量;A为管道截面面积;k为局部阻力系数,其值按流体力学的有关图表通过样条插值得到。4.水轮机特性处理在水电站过渡过程数值仿真中,需要利用水轮机的全特性曲线,包括飞逸特性曲线和小开度的等值线。但是,目前转轮模型试验工作尚无法满足这一要求,大多数转轮没有这种全特性。对于紫坪铺水电站,文献[8]提供了模型转轮飞逸特性曲线,n=f(a,H),Q=f(a,H)及=f(a,H)等曲线。作者据此计算出了真机Q11=Q1(1a,n11)及Mt11=Mt11(a,n11)曲线,并内插补充了水轮机小开度区的不足部分。其中,n为水轮机转速,n11为单位转速,nr为额定转速;Q为水轮机流量,Q11为单位流量;H为水轮机水头;a为水轮机导叶开度;Mt11为水轮机单位力矩;为水轮机效率。图1混流式水轮机模型综合特性曲线图2混流式水轮机模型综合特性曲线的三维曲面为了便于插值并保证足够精度,在水轮机特性曲线处理中,作者应用Matlab语言中的meshgrid函数对根据各开度下Q11=Q11(a,n11)及Mt11=Mt11(a,n11)曲线