CN109657327A - The evaluation method of pump installation outlet passage comprehensive performance - Google Patents

The evaluation method of pump installation outlet passage comprehensive performance Download PDF

Info

Publication number
CN109657327A
CN109657327A CN201811529400.6A CN201811529400A CN109657327A CN 109657327 A CN109657327 A CN 109657327A CN 201811529400 A CN201811529400 A CN 201811529400A CN 109657327 A CN109657327 A CN 109657327A
Authority
CN
China
Prior art keywords
outlet passage
formula
value
pump installation
outlet
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811529400.6A
Other languages
Chinese (zh)
Other versions
CN109657327B (en
Inventor
杨帆
胡文竹
高慧
刘超
施伟
钱军
黄璐
陈晔
王铁力
张海军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangzhou University
Original Assignee
Yangzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yangzhou University filed Critical Yangzhou University
Priority to CN201811529400.6A priority Critical patent/CN109657327B/en
Publication of CN109657327A publication Critical patent/CN109657327A/en
Application granted granted Critical
Publication of CN109657327B publication Critical patent/CN109657327B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Abstract

The invention discloses a kind of evaluation methods of pump installation outlet passage comprehensive performance, CFD Three Dimensional Steady numerical result of this method based on pump installation, determine the subitem evaluation index of 5 hydraulic performances of outlet passage, and scheme is set according to the first of outlet passage, determine the subitem evaluation index of 4 structures based on pump installation outlet passage three dimensional design;It converts to 9 subitem evaluation indexes of outlet passage determined above, the objective function of outlet passage comprehensive performance is obtained, to obtain the opinion rating of outlet passage comprehensive performance.The present invention overcomes in the prior art only to pumping plant outlet passage hydraulic performance evaluation method the drawbacks of, be more advantageous to pumping plant just set stage preferred outlet passage, the energy consumption of pumping plant can be further decreased, achieved the purpose that energy-saving.

Description

The evaluation method of pump installation outlet passage comprehensive performance
Technical field
The present invention relates to Design of Pump Station field, especially a kind of evaluation method of pump installation outlet passage comprehensive performance.
Background technique
Outlet passage is the important water flowing out structure of low-lift-head pumping station, for the rank of connection water pump guide vane body outlet and discharge bay Road is connected, main function is to be diffused the higher speed water flow of diffuser outlet, recycles rotation function;The waterpower of outlet passage Performance directly affects the height of low-lift pump device energy efficiency, and the lift of pump installation is lower, and outlet passage is to pump installation energy The influence of performance is bigger.
Currently, in pumping plant Preliminary design, outlet passage it is preferred usually using the size of hydraulic loss as judging quota, Designer or scientific research personnel often independently seek its hydraulic loss size for outlet passage as computing object, this calculation method is neglected It has omited outlet passage and has interfered with each other influence with what is flowed in water pump, selected using this calculated result as the multi-scheme ratio for judging outlet passage There are apparent technological deficiencies.The pith that outlet passage is designed as pump station structure, if the pumping plant concept phase only It is aobvious to ignore influence of the outlet passage structure size to pumping plant construction investment using the hydraulic loss of outlet passage as judging basis So also have obvious insufficient.
Summary of the invention
Goal of the invention: the purpose of the present invention is to provide a kind of evaluation methods of pump installation outlet passage comprehensive performance, should Evaluation method is easy to be reliable, is conducive to pumping plant and just sets stage preferred outlet passage, further decreases the energy consumption of pumping plant.
Technical solution: in order to achieve the above object, the evaluation method packet of pump installation outlet passage comprehensive performance of the present invention Include following steps:
(1) according to the first three-dimensional mathematical model for setting scheme constructs pump installation of pump installation, according to the design conditions of pump installation Condition carries out the calculating of Three Dimensional Steady numerical value to pump installation using CFD (Computational Fluid Dynamics) method and asks Solution, is below 1.0 × 10 with the residual error convergence precision of each physical quantity-5, and by the variation of monitoring lift, when the variation of lift becomes When definite value, i.e., numerical value calculating meets convergent requirement under the conditions of the design conditions;
(2) it according to the Three Dimensional Steady numerical result of pump installation, calculates the opposite hydraulic loss of outlet passage, go out water flow The opposite average speed of the pressure-recovery factor in road, the axial velocity profile uniformity of outlet passage outlet and outlet passage outlet 4 subitem evaluation indexes;
(3) according to the Three Dimensional Steady numerical result of pump installation, judged by outlet passage interior three-dimensional streamline, if stream Line direction can determine whether that the region is flowed back on the contrary, there is convoluted water flow i.e. with the mainstream direction of motion, determine therefrom that out water flow The quantity in road internal reflux region;
(4) scheme is set according to the first of outlet passage, the relative equivalent angle of flare θ of outlet passage and outlet passage exit face Radius R is drawn in changeout
(5) outlet passage bending section number determines the complexity of outlet passage structure;According to the cutout side of outlet passage The evaluation index for being provided with cutout mode and inspection hole without determining outlet passage of formula and inspection hole;
(6) determine each cross-sectional area of outlet passage along the situation of change of journey;If the two is linear relationship variation, this evaluation Index is 100;If the two is 6 order polynomial curvilinear motion relationships, this evaluation index is 50;If both sides relation curve is between line Between property variation relation and 6 order polynomial curvilinear motion relationships, then according to the evenness of both sides relation curvilinear motion and two The complexity of person's variation relation fitting formula determines evaluation index numerical value;
(7) with the opposite hydraulic loss of above-mentioned outlet passage, the pressure-recovery factor of outlet passage, outlet passage outlet The axial velocity profile uniformity, outlet passage outlet opposite average speed, the quantity in outlet passage internal reflux region, out The equivalent angle of flare of water flow passage, the complexity of outlet passage structure, the cutout mode and inspection hole of outlet passage, outlet passage Each cross-sectional area changes evaluation index of this 9 subitem evaluation indexes as building outlet passage overall performance rating along journey CPEI(Comprehensive Performance Evaluation Indicators)。
Wherein, in step (7), the comprehensive performance evaluation index CPEI membership function of outlet passage is divided into 4 grades, Segmentation point value y (i) on section [0,1] takes 1.0,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2 and 0.1 respectively;It is corresponding Each subitem evaluation index value be defined as x*(i, j), wherein i=1~10, j=1~9.
Further, the calculating formula of each subitem evaluation index is as follows in step (2):
The opposite hydraulic loss Δ h of outlet passageR:
ΔhR=Δ h/H (1)
In formula: Δ h is the hydraulic loss of outlet passage, unit m;H is the net lift of pump installation, unit m;
The pressure-recovery factor ξ of outlet passage:
ξ=Ecout/Ecin (2)
In formula: EcoutFor the stagnation pressure of outlet passage exit face, unit Pa;EcinIt is single for the stagnation pressure of outlet passage inlet face Position is Pa;
The axial velocity profile uniformity V of outlet passage outletu+Are as follows:
In formula:For the axial velocity arithmetic mean of instantaneous value of outlet passage exit face, unit m/s;vaiFor outlet passage The axial velocity of each grid node on exit face, unit m/s;N is the sum of grid node on outlet passage exit face;B is The number of grid node on outlet passage exit face, the integer that b is 1 to N;
The opposite average speed v of outlet passage outletaRAre as follows:
VaR=Va/vmax (4)
In formula: vaFor the average speed of outlet passage exit face, unit m/s;vmaxFor the limitation speed of outlet passage outlet Degree, when outlet passage cutout facility is flap valve, vmaxTake 2.0m/s;When outlet passage cutout facility is non-flap valve, vmaxIt takes 1.5m/s。
Further, in step (4), the calculating formula of the relative equivalent angle of flare θ of outlet passage is
In formula: RinFor the radius of outlet passage inlet face, unit m;RoutRadius is drawn in change for outlet passage exit face, Unit is m;L is the length of outlet passage, unit m;
Radius R is drawn in the change of outlet passage exit faceoutCalculating formula is
In formula: AoutFor the area of outlet passage exit face, unit m2;π is constant, takes 3.1415926.
Further, in step (6), determine each cross-sectional area of outlet passage along the specific steps of the situation of change of journey are as follows: According to each cross-sectional area evolution with distance relation curve of outlet passage, by EXCEL software by each cross-sectional area numerical value of outlet passage List is carried out with the relationship along journey, then Trendline is added to curve, and show Trendline formula and discriminant coefficient, is thereby determined that out Situation of change of each cross-sectional area of water flow passage along journey.
It (8) is the variation of the dimension and unified each index of each subitem evaluation index data of solution outlet passage comprehensive performance Range, by evaluation index x*(i, j) carries out dimensionless processing by following calculating formula and has obtained standardized value x (i, j):
(9) 9 dimension datas { x (i, j) | j=1~9 } are comprehensive at One Dimensional Projection value z*(i):
In formula: a is unit length projecting direction vector,
Projecting direction is according to probability a2(j) it is distributed, takes { a2(j) } the maximum probability distribution of information entropy;Project target function Q (a) are as follows:
Q (a)=S | R | W (9)
In formula: S is the standard deviation of projection value z (i);R is the related coefficient of z (i) and y (i);W is a2(j) comentropy, S, the calculating formula of R and W is respectively as follows:
WhereinFor the average value of y (i),For the average value of z (i);
(10) best projection direction is estimated using solution projection target function maximization problems:
MaxQ (a)=S | R | W (13)
Constraint condition:
(11) the best projection direction a that optimizing index function is obtained*(j) calculating formula (2) are substituted into and obtains i-th of standard The projection value z of cut-point*(i), according to z*(i)~y (i) scatter plot constructs mathematical model, which is outlet passage The evaluation model of comprehensive performance:
y*(i)=f [z*(i)] (14)
(12) to the indices value x of pump installation outlet passage*(j) after nondimensionalization is handled and a of best projection direction* Projection value z can be acquired by substituting into calculating formula (8) together*, then by z*It substitutes into calculating formula (11), it is comprehensive to obtain outlet passage to be evaluated Close the calculated value y of the degree of membership of performance*, further according to y*The concrete condition in 4 subintervals is fallen in, it is comprehensive to obtain the outlet passage The opinion rating of performance.
Further, the comprehensive performance evaluation index CPEI membership function criteria for classifying of outlet passage are as follows: when y > 0.9, The overall performance rating of outlet passage is excellent;When 0.9 >=y > 0.8, the overall performance rating of outlet passage is good;0.8 >=y > When 0.7, during the overall performance rating of outlet passage is;When 0.7 >=y, the overall performance rating of outlet passage is poor.
Further, in step (9), when by 9 dimension datas synthesis at One Dimensional Projection value, make projection value z*(i) objective On possess the variation information of { x (i, j) } as much as possible, make z*(i) absolute value of the related coefficient of standard deviation z (i) and y (i) Reach as big as possible.
The utility model has the advantages that compared with prior art, the present invention has following remarkable advantage:
(1) evaluation method of outlet passage comprehensive performance provided by the invention is conducive to set the stage at the beginning of pumping plant outlet passage Scientifically and rationally select outlet passage, avoid it is current only using outlet passage hydraulic loss as judge outlet passage superiority and inferiority only One standard avoids and only ignores water flow out using the hydraulic loss of outlet passage as judging basis in the pumping plant concept phase Influence of the road structure size to pumping plant construction investment and the obvious deficiency having;
(2) the method overcome previous outlet passage calculate have ignored flowed in outlet passage and water pump interfere with each other influence The drawbacks of, the pith that outlet passage is designed as pump station structure overcomes and comments at present pumping plant outlet passage hydraulic performance The technical deficiency of valence method breaches the calculation method constraint of traditional outlet passage water-power performance parameter;The application of this method is more Be conducive to pumping plant and just set stage preferred outlet passage, further decreases the energy consumption of pumping plant, achieve the purpose that energy-saving.
Detailed description of the invention
Fig. 1 is the flow chart of evaluation method of the present invention;
Fig. 2 is the pump installation three-dimensional model diagram of pumping plant.
Specific embodiment
Technical solution of the present invention is described further with reference to the accompanying drawing.
It is as shown in Figure 1 the flow chart of the outlet passage comprehensive performance evaluation method based on pump installation.The evaluation method packet Include following steps:
(1) according to the first three-dimensional mathematical model for setting scheme constructs pump installation of pump installation, according to the design conditions of pump installation Condition carries out the calculating of Three Dimensional Steady numerical value to pump installation using CFD (Computational Fluid Dynamics) method and asks Solution, is below 1.0 × 10 with the residual error convergence precision of each physical quantity-5, and by the variation of monitoring lift, when the variation of lift becomes When definite value, i.e., numerical value calculating meets convergent requirement under the conditions of the design conditions;
(2) it according to the Three Dimensional Steady numerical result of pump installation, calculates the opposite hydraulic loss of outlet passage, go out water flow The opposite average speed of the pressure-recovery factor in road, the axial velocity profile uniformity of outlet passage outlet and outlet passage outlet The calculating formula of 4 subitem evaluation indexes, each subitem evaluation index is as follows:
The opposite hydraulic loss Δ h of outlet passageR:
ΔhR=Δ h/H (1)
In formula: Δ h is the hydraulic loss of outlet passage, unit m;H is the net lift of pump installation, unit m.
The pressure-recovery factor ξ of outlet passage:
ξ=Ecout/Ecin (2)
In formula: EcoutFor the stagnation pressure of outlet passage exit face, unit Pa;EcinIt is single for the stagnation pressure of outlet passage inlet face Position is Pa.
The axial velocity profile uniformity V of outlet passage outletu+Are as follows:
In formula:For the axial velocity arithmetic mean of instantaneous value of outlet passage exit face, unit m/s;VaiFor outlet passage The axial velocity of each grid node on exit face, unit m/s;N is the sum of grid node on outlet passage exit face;B is The number of grid node on outlet passage exit face, the integer that b is 1 to N.
The opposite average speed V of outlet passage outletaRAre as follows:
VaR=Va/vmax (4)
In formula: VaFor the average speed of outlet passage exit face, unit m/s;vmaxFor the limitation speed of outlet passage outlet Degree, when outlet passage cutout facility is flap valve, vmaxTake 2.0m/s;When outlet passage cutout facility is non-flap valve, vmaxIt takes 1.5m/s。
(3) according to the Three Dimensional Steady numerical result of pump installation, judged by outlet passage interior three-dimensional streamline, if stream Line direction can determine whether that the region is flowed back on the contrary, there is convoluted water flow i.e. with the mainstream direction of motion, determine therefrom that out water flow The quantity in road internal reflux region.
(4) scheme is set according to the first of outlet passage, the calculating formula of the relative equivalent angle of flare θ of outlet passage is
In formula: RinFor the radius of outlet passage inlet face, unit m;RoutRadius is drawn in change for outlet passage exit face, Unit is m;L is the length of outlet passage, unit m.
Radius R is drawn in the change of outlet passage exit faceoutCalculating formula is
In formula: AoutFor the area of outlet passage exit face, unit m2;π is constant, takes 3.1415926.
(5) complexity that outlet passage structure is determined according to outlet passage bending section number, if bending section number is more, The construction degree of flow passage structure is more complicated;According to being provided with without determining outlet passage for the cutout mode of outlet passage and inspection hole Cutout mode and inspection hole evaluation index;
Wherein, commonly the mode of stopping has 3 kinds to outlet passage: stop gate, flap valve, air admission valve, part pumping plant are adopted It is combined use with wherein 2 kinds, such as stop gate and air admission valve, air admission valve and flap valve.Inspection hole is to be discharged For checking that outlet passage cutout facility and structural strength use, easy access personnel are entered in outlet passage for setting on runner The maintenance that inspection hole is convenient for outlet passage is arranged in portion, outlet passage.
The determination that evaluation index is carried out according to the combination of the cutout mode of outlet passage and inspection hole, according to flap valve and nothing Inspection hole, then the safety coefficient of the outlet passage is minimum, evaluation index 1;According to stop gate, air admission valve and inspection Hole, then the safety coefficient of the outlet passage is up to 10, remaining boundary is between the two.
(6) determine each cross-sectional area of outlet passage along the situation of change of journey;If the two is linear relationship variation, this evaluation Index is 100;If the two is 6 order polynomial curvilinear motion relationships, this evaluation index is 50;If both sides relation curve is between line Between property variation relation and 6 order polynomial curvilinear motion relationships, then according to the evenness of both sides relation curvilinear motion and two The complexity of person's variation relation fitting formula determines evaluation index numerical value.
According to each cross-sectional area evolution with distance relation curve of outlet passage, by EXCEL software by each section of outlet passage Area figures and list is carried out along the relationship of journey, then Trendline is added to curve, and shows Trendline formula and discriminant coefficient, example Such as, if each cross-sectional area numerical value and the relation curve along journey are linear relationship y=kx, and discriminant coefficient is greater than 0.95, then evaluates Index value is 100;If power exponent relationship, and discriminant coefficient is greater than 0.95, then evaluation number is 95;If fit correlation is 6 times Multinomial, discriminant coefficient are greater than 0.95, then evaluation number is 60;If not first 9 kinds, then evaluation number is 50;Determine therefrom that water outlet Situation of change of each cross-sectional area of runner along journey.The corresponding relationship of specific evaluation number and distinguishing rule is as shown in the table.
Each cross-sectional area evolution with distance evaluation number of 1 outlet passage of table table corresponding with judgment basis
(7) with the opposite hydraulic loss of above-mentioned outlet passage, the pressure-recovery factor of outlet passage, outlet passage outlet The axial velocity profile uniformity, outlet passage outlet opposite average speed, the quantity in outlet passage internal reflux region, out The equivalent angle of flare of water flow passage, the complexity of outlet passage structure, the cutout mode and inspection hole of outlet passage, outlet passage Each cross-sectional area changes evaluation index of this 9 subitem evaluation indexes as building outlet passage overall performance rating along journey CPEI (Comprehensive Performance Evaluation Indicators), the comprehensive performance of outlet passage is commented Valence index CPEI membership function is divided into 4 grades, and classification standard is as shown in table 2 below, the segmentation point value y (i) on section [0,1] 1.0,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2 and 0.1 are taken respectively.The evaluation index value of corresponding each subitem is defined as x*(i, j), wherein i=1~10, j=1~9.
The cut-point of the evaluation degree of membership of each subitem index of outlet passage comprehensive performance is as shown in table 3.
The overall performance rating of 2 outlet passage of table divides table
The degree of membership cut-point of each subitem index of 3 outlet passage comprehensive performance of table
The opinion rating of the cut-point and outlet passage comprehensive performance of 9 subitem evaluation indexes of above-mentioned outlet passage is drawn Divide is with reference to " Code for design of pumping cstations " (GB50265-2010), the related academic paper research achievement of outlet passage hydraulic performance etc. It is determined.
It (8) is the variation of the dimension and unified each index of each subitem evaluation index data of solution outlet passage comprehensive performance Range, by evaluation index x*(i, j) carries out dimensionless processing by following calculating formula and has obtained standardized value x (i, j):
(9) 9 dimension datas { x (i, j) | j=1~9 } are comprehensive at One Dimensional Projection value z*(i):
In formula: a is unit length projecting direction vector,
In comprehensive projection value, projection value z should be made*(i) the variation letter for possessing { x (i, j) } as much as possible is objectively being needed Breath, it is ensured that z*(i) absolute value of the related coefficient of standard deviation z (i) and y (i) reaches as big as possible.Meanwhile projecting direction according to Probability a2(j) it is distributed, according to Maximum Entropy Principle Method, entropy means that more greatly the constraint of human factor addition and hypothesis are fewer, Take { a2(j) } the maximum probability distribution of information entropy.It projects target function Q (a) are as follows:
Q (a)=S | R | W (9)
In formula: S is the standard deviation of projection value z (i);R is the related coefficient of z (i) and y (i);W is a2(j) comentropy, S, the calculating formula of R and W is respectively as follows:
WhereinFor the average value of y (i),For the average value of z (i);
(10) best projection direction is estimated using solution projection target function maximization problems:
MaxQ (a)=S | R | W (13)
Constraint condition:
(11) the best projection direction a that optimizing index function is obtained*(j) substituting into calculating formula (2) can be obtained i-th of mark The projection value z of quasi- cut-point*(i), according to z*(i)~y (i) scatter plot constructs mathematical model, which is water flow The evaluation model of trace integration performance:
y*(i)=f [z*(i)] (14)
(12) to the indices value x of pump installation outlet passage*(j) after nondimensionalization is handled and a of best projection direction* Projection value z can be acquired by substituting into calculating formula (8) together*, then by z*It substitutes into calculating formula (11), outlet passage to be evaluated can be obtained The calculated value y of the degree of membership of comprehensive performance*, further according to y*The concrete condition in 4 subintervals is fallen in, it is comprehensive to judge the outlet passage Close the opinion rating of performance.
It is a specific embodiment using the evaluation method below.The pump installation threedimensional model of certain pumping plant as shown in Fig. 2, The pump installation carries out the calculating of Three Dimensional Steady numerical value when to design conditions, when the residual error convergence precision of each physical quantity is below 1.0 × 10-5, and when the variation of pump installation lift tends to definite value, it is believed that pump installation Three Dimensional Steady numerical value calculates full under the conditions of the design conditions Sufficient convergent requirement.Based on numerical result, by formula (1), formula (2), formula (3) and formula (4) to the opposite waterpower of outlet passage What loss, the axial velocity profile uniformity of the pressure-recovery factor of outlet passage, outlet passage outlet and outlet passage exported Opposite average speed is solved;The quantity in outlet passage internal reflux region is determined when to design conditions;It is filled according to pump The Preliminary design set calculates the relative equivalent angle of flare of outlet passage, determines the complexity of outlet passage structure, outlet passage Cutout mode and inspection hole setting, determine each cross-sectional area of outlet passage along the situation of change of journey.Using formula (7)~formula (13) each subitem evaluation parameter of the outlet passage is handled.It is computed, 9 son evaluations of the outlet passage of the pump installation refer to It is as shown in table 4 to mark parameter.
9 subsystem assessment indicator parameters of 4 outlet passage of table
Standardized value can be obtained after subitem evaluation index value each in table 3 is carried out no quantization processing by calculating formula (7), then pass through Formula (8)~calculating formula (13) can be calculated maxQ (a)=2.10, use genetic algorithm that best projection direction can be obtained as a*(j) =(0.3327,0.3369,0.3281,0.3336,0.336,0.3343,0.3395,0.3336,0.3274), by above-mentioned calculating As a result the projection value z that 10 standard scores cutpoints can be obtained in calculating formula (8) is substituted into*(i)=(3.0021,2.7553,2.5051, 2.2583,2.0082,1.7614,1.4325,1.0766,0.6475,0.1501).Using least square method to z* (i) and y (i) Carry out fitting of a polynomial, the mathematical regression equation after fitting are as follows: y*(i)=- 0.0118 [z* (i)]4+0.0688[z*(i)]3- 0.0725[z*(i)]2+0.2283z*(i)+0.067。
The coefficient R of the regression equation2=0.9999.To every evaluation index (such as table of the pump installation outlet passage 4) no quantization processing result is carried out are as follows: x (j)=(0.9806,0.7889,0.9472,0.8944,0.8311,1.0000, 0.9412,1.0000,1.0000), by each subitem evaluation index nondimensionalization processing result of the outlet passage and best projection side To calculating formula (8) are substituted into together, the projection value that each subitem evaluation index current value can be obtained is z*=2.7953, which is substituted into The degree of membership y for the comprehensive performance that can obtain the pump installation outlet passage is computed into the mathematical regression equation of above-mentioned acquisition*= 0.921.According to the overall performance rating demarcation interval of 2 outlet passage of table, the outlet passage overall performance rating of the pump installation is It is excellent.

Claims (7)

1. a kind of evaluation method of pump installation outlet passage comprehensive performance, it is characterised in that: the following steps are included:
(1) according to the first three-dimensional mathematical model for setting scheme constructs pump installation of pump installation, according to the design conditions condition of pump installation, Three Dimensional Steady numerical value is carried out to pump installation using computational fluid dynamics and calculates solution;
(2) according to the Three Dimensional Steady numerical result of pump installation, the opposite hydraulic loss of outlet passage, outlet passage are calculated Opposite average speed 4 of the axial velocity profile uniformity and outlet passage outlet that pressure-recovery factor, outlet passage export Subitem evaluation index;
(3) according to the Three Dimensional Steady numerical result of pump installation, outlet passage is determined by outlet passage interior three-dimensional streamline The quantity in internal reflux region;
(4) scheme is set according to the first of outlet passage, calculates the relative equivalent angle of flare θ and outlet passage exit face of outlet passage Radius R is drawn in changeout
(5) foundation just sets the complexity that outlet passage bending section number in scheme determines outlet passage structure;According to water flow out The cutout mode in road and the evaluation index for being provided with cutout mode and inspection hole without determining outlet passage of inspection hole;
(6) determine each cross-sectional area of outlet passage along the situation of change of journey;If the two is linear relationship variation, this evaluation index It is 100;If the two is 6 order polynomial curvilinear motion relationships, this evaluation index is 50;If both sides relation curve becomes between linear Between change relationship and 6 order polynomial curvilinear motion relationships, then become according to the evenness of both sides relation curvilinear motion and the two The complexity of change relationship fitting formula determines evaluation index numerical value;
(7) axis exported with the opposite hydraulic loss of above-mentioned outlet passage, the pressure-recovery factor of outlet passage, outlet passage The opposite average speed exported to the VELOCITY DISTRIBUTION uniformity, outlet passage, the quantity in outlet passage internal reflux region, out water flow The equivalent angle of flare in road, the complexity of outlet passage structure, the cutout mode and inspection hole of outlet passage, outlet passage are each disconnected Face area changes evaluation index CPEI of this 9 subitem evaluation indexes as building outlet passage overall performance rating along journey; The comprehensive performance evaluation index CPEI membership function of outlet passage is divided into 4 grades, the segmentation point value y on section [0,1] (i) 1.0,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2 and 0.1 are taken respectively;Corresponding each subitem evaluation index value definition For x*(i, j), wherein i=1~10, j=1~9.
2. evaluation method according to claim 1, it is characterised in that: comprehensive performance evaluation index CPEI of outlet passage etc. The grade degree of membership criteria for classifying are as follows: when y > 0.9, the overall performance rating of outlet passage is excellent;When 0.9 >=y > 0.8, outlet passage Overall performance rating be it is good;When 0.8 >=y > 0.7, during the overall performance rating of outlet passage is;When 0.7 >=y, outlet passage Overall performance rating be it is poor.
3. evaluation method according to claim 1, it is characterised in that: the calculating formula of each subitem evaluation index in step (2) It is as follows:
The opposite hydraulic loss Δ h of outlet passageR:
ΔhR=Δ h/H (1)
In formula: Δ h is the hydraulic loss of outlet passage, unit m;H is the net lift of pump installation, unit m;
The pressure-recovery factor ξ of outlet passage:
ξ=Ecout/Ecin (2)
In formula: EcoutFor the stagnation pressure of outlet passage exit face, unit Pa;EcinFor the stagnation pressure of outlet passage inlet face, unit is Pa;
The axial velocity profile uniformity V of outlet passage outletu+Are as follows:
In formula:For the axial velocity arithmetic mean of instantaneous value of outlet passage exit face, unit m/s;vaiFor outlet passage exit face The axial velocity of upper each grid node, unit m/s;N is the sum of grid node on outlet passage exit face;B is water flow The number of grid node on road exit face, the integer that b is 1 to N;
The opposite average speed v of outlet passage outletaRAre as follows:
vaR=va/vmax (4)
In formula: vaFor the average speed of outlet passage exit face, unit m/s;vmaxFor the limitation speed of outlet passage outlet.
4. evaluation method according to claim 1, it is characterised in that: in step (4), the relative equivalent of outlet passage is spread The calculating formula of angle θ is
In formula: RinFor the radius of outlet passage inlet face, unit m;RoutDraw radius, unit for the change of outlet passage exit face For m;L is the length of outlet passage, unit m;
Radius R is drawn in the change of outlet passage exit faceoutCalculating formula is
In formula: AoutFor the area of outlet passage exit face, unit m2;π is constant.
5. evaluation method according to claim 1, it is characterised in that: in step (6), determine each cross-sectional area of outlet passage Along the specific steps of the situation of change of journey are as follows: soft by EXCEL according to each cross-sectional area evolution with distance relation curve of outlet passage Part carries out list by each cross-sectional area numerical value of outlet passage and along the relationship of journey, then adds Trendline to curve, and show trend Line formula and discriminant coefficient thereby determine that each cross-sectional area of outlet passage along the situation of change of journey.
6. evaluation method according to any one of claims 1-5, it is characterised in that: further comprising the steps of:
(8) by evaluation index x*(i, j) carries out dimensionless processing by following calculating formula and has obtained standardized value x (i, j):
(9) 9 dimension datas { x (i, j) | j=1~9 } are comprehensive at One Dimensional Projection value z*(i):
In formula: a be unit length projecting direction vector, a (j) > 0,
Projecting direction is according to probability a2(j) it is distributed, takes { a2(j) } the maximum probability distribution of information entropy;It projects target function Q (a) Are as follows:
Q (a)=S | R | W (9)
In formula: S is the standard deviation of projection value z (i);R is the related coefficient of z (i) and y (i);W is a2(j) comentropy, S, R and W Calculating formula be respectively as follows:
WhereinFor the average value of y (i),For the average value of z (i);
(10) best projection direction is estimated using solution projection target function maximization problems:
MaxQ (a)=S | R | W (13)
Constraint condition:
(11) the best projection direction a that optimizing index function is obtained*(j) calculating formula (2) are substituted into and obtains i-th of Standard Segmentation The projection value z of point*(i), according to z*(i)~y (i) scatter plot constructs mathematical model, which is that outlet passage is comprehensive The evaluation model of performance:
y*(i)=f [z*(i)] (14)
(12) to the indices value x of pump installation outlet passage*(j) after nondimensionalization is handled and a of best projection direction*Together Projection value z can be acquired by substituting into calculating formula (8)*, then by z*It substitutes into calculating formula (11), it is comprehensive to obtain outlet passage to be evaluated The calculated value y of the degree of membership of energy*, further according to y*The concrete condition in 4 subintervals is fallen in, the outlet passage comprehensive performance is obtained Opinion rating.
7. evaluation method according to claim 6, it is characterised in that: in step (9), by the comprehensive synthesizing one-dimensional of 9 dimension datas When projection value, make projection value z*(i) in the variation information for objectively possessing { x (i, j) } as much as possible, make z*(i) standard deviation z (i) and the absolute value of the related coefficient of y (i) reaches as big as possible.
CN201811529400.6A 2018-12-13 2018-12-13 The evaluation method of pump installation outlet passage comprehensive performance Active CN109657327B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811529400.6A CN109657327B (en) 2018-12-13 2018-12-13 The evaluation method of pump installation outlet passage comprehensive performance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811529400.6A CN109657327B (en) 2018-12-13 2018-12-13 The evaluation method of pump installation outlet passage comprehensive performance

Publications (2)

Publication Number Publication Date
CN109657327A true CN109657327A (en) 2019-04-19
CN109657327B CN109657327B (en) 2019-10-29

Family

ID=66114260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811529400.6A Active CN109657327B (en) 2018-12-13 2018-12-13 The evaluation method of pump installation outlet passage comprehensive performance

Country Status (1)

Country Link
CN (1) CN109657327B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110593187A (en) * 2019-09-11 2019-12-20 扬州大学 Series water outlet flow passage of front-located shaft type tubular pump device and application method thereof
CN112270164A (en) * 2020-11-10 2021-01-26 湖南凯利特泵业有限公司 Pump station operation information visualization method for energy efficiency assessment
CN113033115A (en) * 2021-03-04 2021-06-25 惠达卫浴股份有限公司 Toilet flushing efficiency evaluation method
CN113987948A (en) * 2021-11-03 2022-01-28 江苏四联自动化科技有限公司 Intelligent measuring and calculating method and system for outlet water flow of pump station

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009155903A (en) * 2007-12-26 2009-07-16 Caterpillar Japan Ltd Front control method of working machine
CN104265652A (en) * 2014-07-25 2015-01-07 扬州大学 Collaborative optimization design method for large-flow low-lift pump station
CN104615831A (en) * 2015-02-10 2015-05-13 扬州大学 Shaping algorithm for quartic curvilinear equation of three-dimensional shape of low hump type water draining flow way and application of shaping algorithm
CN104636542A (en) * 2014-12-31 2015-05-20 扬州大学 Method for predicting energy performance of pump through adjustable guide vane on basis of CFD
CN105179303A (en) * 2015-10-24 2015-12-23 扬州大学 Axial flow pump impeller all-operating-condition design method
CN105608287A (en) * 2016-01-12 2016-05-25 扬州大学 Design method for water outlet flow passage three-dimensional-form flow surface of high-performance and large-flow pump station
CN105740501A (en) * 2016-01-20 2016-07-06 扬州大学 Steady numerical calculation-based hydraulic stability analysis method for pump device

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009155903A (en) * 2007-12-26 2009-07-16 Caterpillar Japan Ltd Front control method of working machine
CN104265652A (en) * 2014-07-25 2015-01-07 扬州大学 Collaborative optimization design method for large-flow low-lift pump station
CN104636542A (en) * 2014-12-31 2015-05-20 扬州大学 Method for predicting energy performance of pump through adjustable guide vane on basis of CFD
CN104615831A (en) * 2015-02-10 2015-05-13 扬州大学 Shaping algorithm for quartic curvilinear equation of three-dimensional shape of low hump type water draining flow way and application of shaping algorithm
CN105179303A (en) * 2015-10-24 2015-12-23 扬州大学 Axial flow pump impeller all-operating-condition design method
CN105608287A (en) * 2016-01-12 2016-05-25 扬州大学 Design method for water outlet flow passage three-dimensional-form flow surface of high-performance and large-flow pump station
CN105740501A (en) * 2016-01-20 2016-07-06 扬州大学 Steady numerical calculation-based hydraulic stability analysis method for pump device

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
SUN-SHENG YANG 等: "Theoretical, numerical and experimental prediction of pump as turbine performance", 《RENEWABLE ENERGY》 *
杨帆 等: "不同转速时贯流泵装置出水流道水力性能分析", 《水动力学研究与进展》 *
杨帆 等: "轴流泵装置出水流道内流脉动与流动噪声试验分析", 《农业机械学报》 *
杨帆: "低扬程泵装置水动力特性及多目标优化关键技术研究", 《万方数据库.学位论文库》 *
谢荣盛 等: "双向竖井贯流泵进出水流道优化研究", 《农业机械学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110593187A (en) * 2019-09-11 2019-12-20 扬州大学 Series water outlet flow passage of front-located shaft type tubular pump device and application method thereof
CN110593187B (en) * 2019-09-11 2020-10-20 扬州大学 Series water outlet flow passage of front-located shaft type tubular pump device and application method thereof
CN112270164A (en) * 2020-11-10 2021-01-26 湖南凯利特泵业有限公司 Pump station operation information visualization method for energy efficiency assessment
CN112270164B (en) * 2020-11-10 2021-06-08 湖南凯利特泵业有限公司 Pump station operation information visualization method for energy efficiency assessment
CN113033115A (en) * 2021-03-04 2021-06-25 惠达卫浴股份有限公司 Toilet flushing efficiency evaluation method
CN113033115B (en) * 2021-03-04 2023-04-21 惠达卫浴股份有限公司 Toilet flushing efficiency evaluation method
CN113987948A (en) * 2021-11-03 2022-01-28 江苏四联自动化科技有限公司 Intelligent measuring and calculating method and system for outlet water flow of pump station

Also Published As

Publication number Publication date
CN109657327B (en) 2019-10-29

Similar Documents

Publication Publication Date Title
CN109657327B (en) The evaluation method of pump installation outlet passage comprehensive performance
CN113236581B (en) Intelligent parallel pump system and optimal adjustment method
CN103383121B (en) Regional air conditioner dispersion secondary pumping system
CN105805043B (en) A kind of design method of the non-adjustable axial-flow pump impeller with deviated splitter vane feature
CN108268742A (en) Industrial circulating water impeller of pump double optimization energy-conserving reconstruction method
CN108537449A (en) Meter and river are passed the flood period the reservoir coordinated scheduling strategy acquisition methods of demand
CN1186572C (en) Variable pressure difference and variable flow rate control method and system for air conditioning water system
CN105574288B (en) Method for designing water inlet conduit three-dimensional body flow surface of high-performance large-flow pump station
CN107917089B (en) Water circulating pump pumps leaf aperture optimal regulation method
CN111784536B (en) Method for estimating water level overrun time in open channel according to actually measured water level change condition
CN105864016A (en) Variable-water-volume operation regulating method for open multi-water-pump transmission and distribution system
CN105608287B (en) Design method for water outlet flow passage three-dimensional-form flow surface of high-performance and large-flow pump station
CN105718701A (en) Pump type optimized design method based on water pump characteristic three-dimensional spatial patterns of pump station
CN115169803A (en) Irrigation district water transmission and distribution gate group combined scheduling method and system
CN107038295B (en) A kind of hydraulic ram inner flow passage evaluation and optimization method
CN109882446A (en) A kind of design method of Low Specific Speed Centrifugal Pump Impellers splitterr vanes
CN103870620B (en) Parametric modeling method for diffuser and exit guide blade Combined design
CN1757933A (en) Blade type optimized designing method of turbine compression fluid machine
Zhou et al. Influence of different lateral bending angles on the flow pattern of pumping station lateral inflow
CN109344535A (en) A kind of Diffusion Vanes of Axial Pumps body automatic optimizing design method based on self-adaption flow-field
CN110516369A (en) One kind being based on Joukowski aerofoil profile axial-flow pump impeller the Automation Design method
CN113653649B (en) Interstage flow channel structure for improving performance of secondary impeller of multi-stage pump
CN106934104A (en) A kind of method for designing of the siphon outlet conduit of pumping plant
CN104005983A (en) Designing method of three operating points of impeller of high-specific-speed axial flow pump
Wang et al. Influences of flow rate and baffle spacing on hydraulic characteristics of a novel baffle dropshaft

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant