CN109657327B - The evaluation method of pump installation outlet passage comprehensive performance - Google Patents
The evaluation method of pump installation outlet passage comprehensive performance Download PDFInfo
- Publication number
- CN109657327B CN109657327B CN201811529400.6A CN201811529400A CN109657327B CN 109657327 B CN109657327 B CN 109657327B CN 201811529400 A CN201811529400 A CN 201811529400A CN 109657327 B CN109657327 B CN 109657327B
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Flow Control (AREA)
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
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, In
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.
(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 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 One Dimensional Projection value z (i);R is the related coefficient of z (i) and y (i);W is a2(j) information
The calculating formula of entropy, S, 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 (8) 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*(i, j) is after nondimensionalization is handled and best projection side
To a*(j) the projection value z that i-th of standard scores cutpoint can be acquired in calculating formula (8) is substituted into together*(i), then by z*(i) it substitutes into and calculates
In formula (14), the calculated value y of the degree of membership of outlet passage comprehensive performance to be evaluated is obtained*(i), further according to y*(i) 4 sons are fallen in
Concrete condition in section obtains the opinion rating of the outlet passage comprehensive 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>0.7
When, 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 One Dimensional Projection value z (i) in visitor
The variation information for possessing { x (i, j) } in sight as much as possible reaches the standard deviation of z (i) and the absolute value of the related coefficient of y (i)
It is 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,
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 One Dimensional Projection value z (i);R is the related coefficient of z (i) and y (i);W is a2(j) information
The calculating formula of entropy, S, 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 (8) 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*
(j) the projection value z that i-th of standard scores cutpoint can be acquired in calculating formula (8) is substituted into together*(i), then by z*(i) calculating formula is substituted into
(14) in, the calculated value y of the degree of membership of outlet passage comprehensive performance to be evaluated is obtained*(i), further according to y*(i) 4 sub-districts are fallen in
Interior concrete condition, obtains the opinion rating of the outlet passage comprehensive 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 max Q (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 (5)
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;The comprehensive performance evaluation index CPEI membership function of outlet passage divides mark
It is quasi- 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
It is good;When 0.8 >=y > 0.7, during the overall performance rating of outlet passage is;When 0.7 >=y, the overall performance rating of outlet passage is
Difference;
(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, α (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 One Dimensional 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 (8) 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*(i, j) is after nondimensionalization is handled and a of best projection direction*
(j) the projection value z that i-th of standard scores cutpoint can be acquired in calculating formula (8) is substituted into together*(i), then by z*(i) calculating formula is substituted into
(14) in, the calculated value y of the degree of membership of outlet passage comprehensive performance to be evaluated is obtained*(i), further according to y*(i) 4 sub-districts are fallen in
Interior concrete condition, obtains the opinion rating of the outlet passage comprehensive performance.
2. 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.
3. 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.
4. 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.
5. evaluation method described in any one of -4 according to claim 1, it is characterised in that: in step (9), by 9 dimension datas
When synthesis is at One Dimensional Projection value, makes One Dimensional Projection value z (i) in the variation information for objectively possessing { x (i, j) } as much as possible, make z
(i) absolute value of the related coefficient of standard deviation and y (i) reaches as big as possible.
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 CN109657327A (en) | 2019-04-19 |
CN109657327B true 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) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN112270164B (en) * | 2020-11-10 | 2021-06-08 | 湖南凯利特泵业有限公司 | Pump station operation information visualization method for energy efficiency assessment |
CN113033115B (en) * | 2021-03-04 | 2023-04-21 | 惠达卫浴股份有限公司 | Toilet flushing efficiency evaluation method |
CN113987948B (en) * | 2021-11-03 | 2022-11-04 | 江苏四联水务科技有限公司 | Intelligent measuring and calculating method and system for outlet water flow of pump station |
Citations (7)
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 |
-
2018
- 2018-12-13 CN CN201811529400.6A patent/CN109657327B/en active Active
Patent Citations (7)
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)
Title |
---|
Theoretical, numerical and experimental prediction of pump as turbine performance;Sun-Sheng Yang 等;《Renewable Energy》;20120604;第48卷;第507-513页 * |
不同转速时贯流泵装置出水流道水力性能分析;杨帆 等;《水动力学研究与进展》;20171130;第32卷(第6期);第770-785页 * |
低扬程泵装置水动力特性及多目标优化关键技术研究;杨帆;《万方数据库.学位论文库》;20131231;全文 * |
双向竖井贯流泵进出水流道优化研究;谢荣盛 等;《农业机械学报》;20151031;第46卷(第10期);第68-74页 * |
轴流泵装置出水流道内流脉动与流动噪声试验分析;杨帆 等;《农业机械学报》;20180331;第49卷(第3期);第144-149、165页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109657327A (en) | 2019-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109657327B (en) | The evaluation method of pump installation outlet passage comprehensive performance | |
CN106931603B (en) | Central air conditioning cooling water system energy efficiency monitoring system based on internet of things | |
CN113236581B (en) | Intelligent parallel pump system and optimal adjustment method | |
CN103383121B (en) | Regional air conditioner dispersion secondary pumping system | |
CN102486261B (en) | Comprehensive evaluation method of gas pipeline scheduling scheme | |
CN1186572C (en) | Variable pressure difference and variable flow rate control method and system for air conditioning water system | |
CN112215461B (en) | Long-distance water diversion project water delivery scheduling state analysis and early warning system | |
CN111784536B (en) | Method for estimating water level overrun time in open channel according to actually measured water level change condition | |
CN108537449A (en) | Meter and river are passed the flood period the reservoir coordinated scheduling strategy acquisition methods of demand | |
CN109284873A (en) | A kind of wisdom drainage pumping plant carrying out pump group Optimized Operation and early warning disposition based on SWMM model | |
CN107917089B (en) | Water circulating pump pumps leaf aperture optimal regulation method | |
CN105718701A (en) | Pump type optimized design method based on water pump characteristic three-dimensional spatial patterns of pump station | |
CN105864016A (en) | Variable-water-volume operation regulating method for open multi-water-pump transmission and distribution system | |
CN112664464A (en) | Intelligent efficient optimization control method for water pump | |
CN103984999B (en) | A kind of industrial circulating cooling water energy integrated optimization method | |
CN105608287B (en) | Design method for water outlet flow passage three-dimensional-form flow surface of high-performance and large-flow pump station | |
CN115169803A (en) | Irrigation district water transmission and distribution gate group combined scheduling method and system | |
CN103093056A (en) | Optimize system and method of hydropower station water turbine design | |
CN1757933A (en) | Blade type optimized designing method of turbine compression fluid machine | |
CN103870620B (en) | Parametric modeling method for diffuser and exit guide blade Combined design | |
CN108446452A (en) | A kind of mixed-flow pump impeller Robust Optimal Design | |
CN107038295A (en) | A kind of hydraulic ram inner flow passage is evaluated and optimization method | |
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 |
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 |