基于GUI程序的循环水泵特性曲线拟合.doc

基于GUI程序的循环水泵特性曲线拟合 暖研11 马长明 1108521311044 程序编制要求 循环水泵特性曲线拟合：多台同型号水泵并联，通过计算机编程计算，求水泵工作点。 考虑单台或多台循环水泵运行时，比较其流量、扬程、效率、轴功率如何变化，画出水泵和管网特性曲线。 二、已知条件 水泵的几个工况点流量（Q）、扬程（H）、功率（P）和管网的阻力系数（s） 循环水泵曲线拟合原理 从水泵性能曲线图可看出其性能曲线Q一H、一N、一η, 故用三次回归曲线(相比二次回归曲线更精确)对测试数据进行水泵性能回归曲线方程如下： H=H0+A1Q+A2Q2+A3Q3 N=N0+B1Q+B2Q2+B3Q3 (1) η=η0+C1Q+C2Q2+C3Q3 式中： Q——水泵流量，m3/h； H——水泵扬程，mH2O； η——水泵效率，%； N——水泵轴功率，kW H0，A1，A2，A3， N0，B1，B2，B3，η0，C1，C2，C3——回归系数。 对水泵的G一H,回归系数H0，A1，A2，A3可按正规方程组求得,其中n≥4。 同理,也可求出系数N0，B1，B2，B3，η0，C1，C2，C3。 水泵的G-η曲线也可按下式计算： 由最小二乘原理可得三次回归曲线正规方程组，求解该方程组，即可确定上三个方程的系数，从而确定G—H、G—N、G—η—H、—N、—η回归曲线方程； ④ 绘制水泵性能曲线图，横坐标为流量,纵坐标分别为扬程、效率和功率。 多台同型号循环水泵并联运行时Q—H 曲线拟合 n台同型号循环水泵并联运行时Q—H 的拟合曲线如下： H=H0+Q+Q2+Q3 求水泵工作状态点 根据管网中的总阻力和流量Q确定管网的阻力系数S： 式中： S——管网阻力系数，mH2O/(m3h)2； ——管网总阻力，mH2O； 其中流量Q 根据供热系统的热负荷和供回水温度确定。 确定了管网阻力系数以后，就可以绘制出管网特性曲线。根据和水泵曲线求出水泵的运行工况点。function pushbutton1_Callback(hObject, eventdata, handles) A=str2num(get(handles.edit1,string)); % 供热面积 AA=str2num(get(handles.edit2,string)); % 面积热指标 tg=str2num(get(handles.edit3,string)); % 供水温度 th=str2num(get(handles.edit4,string)); % 回水温度 p=str2num(get(handles.edit5,string)); % 管网总阻力 Q=A*AA; G=Q/4187/1000/(tg-th)*3600 % 管网流量 S=p/G^2 % 管网总阻抗 q=0:0.5:500; h=(q.^2).*S; plot(q,h,r); % 绘制管网特性曲线 hold on; set(handles.edit6,string,strcat(num2str(G))); set(handles.edit7,string,strcat(num2str(S))); （2）“拟合并求取状态点”按钮 function pushbutton2_Callback(hObject, eventdata, handles) % 水泵选型 x1=get(findobj(tag,radiobutton1),Value); % 3BL-6 x2=get(findobj(tag,radiobutton2),Value); % 3BL-5A x3=get(findobj(tag,radiobutton3),Value); % wilo-N100/315 x4=get(findobj(tag,radiobutton4),Value); % wilo-N100/400 N=get(findobj(tag,popupmenu1),Value); % 水泵并联台数 if x1==1 q=[30,45,60,70]; % 流量 h=[62,57,50,44.5]; % 扬程 n=[9.3,11,12.3,13.3]; % 功率 e=[55.56,64.77,67.75,65.06]; % 效率 end if x2==1 q=[77,99,119,156]; h=[39.5,39,38.5,36.5]; n=[1