CN110175403A - It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system - Google Patents
It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system Download PDFInfo
- Publication number
- CN110175403A CN110175403A CN201910447969.6A CN201910447969A CN110175403A CN 110175403 A CN110175403 A CN 110175403A CN 201910447969 A CN201910447969 A CN 201910447969A CN 110175403 A CN110175403 A CN 110175403A
- Authority
- CN
- China
- Prior art keywords
- unit
- simulation
- refrigeration
- condenser
- module
- 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.)
- Pending
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F25—REFRIGERATION OR COOLING; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS; MANUFACTURE OR STORAGE OF ICE; LIQUEFACTION SOLIDIFICATION OF GASES
- F25B—REFRIGERATION MACHINES, PLANTS OR SYSTEMS; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS
- F25B49/00—Arrangement or mounting of control or safety devices
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Mechanical Engineering (AREA)
- Tourism & Hospitality (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Thermal Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Air Conditioning Control Device (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provide it is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system, comprising: user's input module, for input freeze or air-conditioning equipment configuration parameter, structural parameters and external environment correlated inputs parameter;Dynamic simulation computing module, the simulation algorithm needed when for selecting and solve simulation model are simultaneously arranged convergence precision, simulation step length and simulation time and use simulation model, carry out simulation calculation according to the parameter inputted in user's input module;Parameter display module, for showing the simulation result obtained by the dynamic simulation computing module and the corresponding dynamic changing curve of partial simulation calculated result;System overall control module is successively run for control system according to the sequence of user's input module, dynamic simulation computing module and parameter display module.Technical solution of the present invention solves the problems, such as that the computer system of current Refrigeration & Air-Conditioning equipment can not achieve the dynamic property prediction of system.
Description
Technical field
The present invention relates to the Computer-aided Design Technology fields of Refrigeration & Air-Conditioning equipment, specifically, more particularly to one
Kind is applied to the parameterisable dynamic simulation computing system of refrigeration or air-conditioning equipment.
Background technique
Using computer aided technique carry out refrigeration equipment emulation and simulation have become refrigeration product research and development, design
Important means.Computer technology constantly develops in recent years, and the memory space and computing capability of computer have all obtained significantly
It improves, analogue simulation has also obtained the concern of more and more enterprises, scientific research institution, and the design based on computer aided technique is ground
Hair work is also gradually generally got up.
As the improvement of people's living standards, the production and application of refrigeration, air-conditioning equipment are developed rapidly, China is
The use big country of refrigeration, air-conditioning products, annual output accounts for 2/3 or so of whole world total output, while the power consumption that summer uses exists
Many cities are more than the 40% of total electricity consumption.Therefore, the improvement of refrigeration, the design efficiency of air-conditioning equipment and properties of product is China's system
One important topic of cold air-conditioning circle.In the design of actual refrigeration product, generally requires to make a large amount of model machine and be tested,
The deficiency of design is made up by physical simulation, but needs take a considerable amount of time and expense.In addition, multiple refrigeration, sky
Adjusting equipment is that structure is complicated, system is huge, inner parameter coupling is strong, the complicated and diversified complicated refrigeration of external boundary condition for one kind
System really cannot comprehensively reflect system performance in the method that typical condition experiment is designed research to multiple system,
And the required test period is long, and cost is also very high.In this regard, the product design of current refrigeration, air-conditioning equipment is sent out towards digitlization
The trend of exhibition is inevitable, exploitation refrigeration, the performance simulation software of air-conditioning equipment and the product optimization design based on performance simulation
Platform is with a wide range of applications.
Current refrigeration, air-conditioning device computer system can not achieve the dynamic of system mostly based on static Simulation model
State performance prediction and a feasible emulation platform is provided for the efficiency that optimizes the system operation, while also explicitly will be not multi-joint
Formula system design work is taken into account.
Summary of the invention
The dynamic property that can not achieve system according to the computer system of current Refrigeration & Air-Conditioning equipment set forth above is pre-
The technical issues of survey, and provide it is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system, both had
Multiple system dynamic simulation function can realize the refrigeration of parametrization, air-conditioning device simulation calculation again.
The technological means that the present invention uses is as follows:
It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system, comprising:
User's input module, for inputting configuration parameter, structural parameters and the external environment phase of refrigeration or air-conditioning equipment
Close input parameter;
Dynamic simulation computing module, when for selecting to solve simulation model the simulation algorithm that needs and be arranged convergence precision,
Simulation step length and simulation time simultaneously use simulation model, carry out simulation calculation according to the parameter inputted in user's input module;
Parameter display module, for showing the simulation result obtained by the dynamic simulation computing module and part
The corresponding dynamic changing curve of simulation result;
System overall control module, for control system according to user's input module, dynamic simulation computing module and parameter
The sequence of display module is successively run.
It further, further include system initialization module, before system is run each time, system overall control module control
System initialization module operation processed, user's input module, dynamic simulation computing module and parameter display module are initialized as writing from memory
Recognize value;It further include data and curves print module, after parameter display module shows simulation result and dynamic changing curve,
The content that system overall control module control data and curves print module operation print parameters display module is shown.
Further, refrigeration or air-conditioning equipment include evaporator, condenser, refrigeration compressor and expansion valve, configuration parameter
Refer to the quantity of refrigeration compressor, condenser, evaporator and expansion valve, and meet following relationship:
(1) quantity of refrigeration compressor is equal to the quantity of condenser;
(2) quantity of expansion valve is equal to the quantity of evaporator.
Further, the structural parameters of refrigeration compressor include cylinder diameter, number of cylinders and length of stroke;The structure of condenser
Parameter includes bore, pipe outside diameter, length of tube, pipe wall material specific heat at constant pressure and pipe wall material density;The structural parameters of evaporator
Including bore, pipe outside diameter, spacing of fin, fin thickness, effective length, pipe wall material specific heat at constant pressure and pipe wall material density;
The structural parameters of expansion valve include unit refrigerant type, booting electronic valve pretravel value, valve body minimum aperture setting, valve body most
The flow area of big aperture setting and valve body standard-sized sheet is arranged.
Further, external environment correlated inputs parameter refers to the setup parameter and condenser of refrigeration compressor and expansion valve
It is cooling with the state parameter of evaporator periphery medium, including clammy compressor input speed, condenser cooling water flow, condenser
Water goes out temperature, evaporator air inlet mass flow, evaporator inlet air temperature and the practical aperture of electronic valve into temperature, condenser cooling water.
Further, the simulation model of dynamic simulation computing module includes:
(1) Simulation Calculation of evaporator and condenser is indicated by following formula:
Refrigerant quality keeps weighing apparatus equation:
Refrigerant energy keeps weighing apparatus equation:
Tube wall energy conservation equation:
ai: for condenser or the evaporator tube inner wall unit area coefficient of heat transfer, unit is kW/ (m2·K);
ao: for condenser or the evaporator tube outer wall unit area coefficient of heat transfer, unit is kW/ (m2·K);
Di: for condenser or evaporator tube internal diameter, unit m;
Do: for condenser or evaporator tube outer diameter, unit m;
Tb: for the temperature of cooling medium, unit K;
Tr: for the temperature of refrigerant, unit K;
Tt: for condenser or the tube wall temperature of evaporator, unit K;
ρ: for the density of refrigerant, unit kg/m3;
U: being refrigerant along condenser or the flow velocity of evaporator tube length direction, unit m/s;
P: for the refrigerant pressure in condenser or evaporator tube, unit Pa;
H: for the refrigerant specific enthalpy in condenser or evaporator tube, unit kJ/kg;
Cp: for condenser or the specific heat at constant pressure of evaporator tube wall material, unit kJ/kgK;
ρt: for the density of condenser or evaporator tube wall material, unit kg/m3;
At: for the circular ring section of condenser or evaporator tube wall material product, unit m2;
(2) Simulation Calculation of refrigeration compressor is indicated by following formula:
Vt=π D2Szn/240
Td=Ts(Pc/Pe)(n-1)/n
Vt: for the theoretical displacement flow of refrigeration compressor, unit m3/s;
D: for the cylinder bore of refrigeration compressor, unit m;
S: for the piston stroke of refrigeration compressor, unit m;
Z: for the speed of crankshaft of refrigeration compressor, unit r/min;
N: for the number of cylinders of refrigeration compressor;
mc: for the mass flow of refrigeration compressor, unit kg/s;
λ: for the volumetric efficiency of refrigeration compressor;
vs: for the aspiration ratio volume of refrigerant, unit m3/kg;
Td: for the delivery temperature of refrigeration compressor, unit K;
Ts: for the suction temperature of refrigeration compressor, unit K;
N: for the polytropic exponent of refrigerant compression process;
(3) Simulation Calculation of expansion valve is indicated by following formula:
hiv=hov
mv: for the refrigerant flow for flowing through expansion valve, unit kg/s;
Cv: for the discharge coefficient of expansion valve;
Av: for the circulation area of expansion valve, unit m2;
ρ: for the density of expansion valve inlet refrigerant liquid, unit kg/m3;
p1: for the pressure of expansion valve inlet refrigerant, unit Pa;
p2: for the pressure of expansion valve outlet refrigerant, unit Pa;
hiv: for expansion valve import enthalpy, unit kJ/kg;
hov: for expansion valve outlet enthalpy, unit kJ/kg.
Further, dynamic simulation computing module further includes that multiple refrigeration or air-conditioning equipment are established using fluid network method
Pipe net leakage rate, can be stated with following formula:
Discharge relation between node and branch
Energy relationship between node and branch
AI: it indicates to flow into matrix, wherein element aijIndicate that node i is the inflow node of branch j, value is 0 or 1;
A0: indicate outflow matrix, wherein element aijIndicate that node i is the outflow node of branch j, value is 0 or 1.
Further, the simulation result that dynamic simulation module obtains includes: the refrigerant mass flow of refrigeration compressor
Amount and delivery temperature;Overheat segment length, two-phase section length, the supercooling segment length, condensing pressure, condensator outlet enthalpy of condenser
Value, superheat section tube wall temperature, two-phase section tube wall temperature and super cooled sect tube wall temperature;The two-phase segment length of evaporator, overheat segment length
Degree, evaporating pressure, evaporator outlet enthalpy, two-phase section tube wall temperature and superheat section tube wall temperature;And the refrigerant of expansion valve
Mass flow.
Further, in dynamic simulation computing module, selective simulation algorithm includes Euler method, improved Euler method and dragon
Ge Kutafa;The value of convergence precision is 0.0001;The value of simulation step length is 0.02 second;The value of simulation time is 500s.
Compared with the prior art, the invention has the following advantages that
The parameterisable dynamic simulation computing system provided by the invention for being applied to refrigeration or air-conditioning equipment, can instruct skill
Art personnel reasonably select the component of related model, specification to reach pre- in design, the large-scale refrigeration of exploitation or air-conditioning system
The refrigeration effect of phase, predicts the working performance of system under different operating conditions, and for energy saving of system, improve gross data be provided, separately
Outside, the exploitation for large-scale refrigeration, air-conditioning simulation training system provides the template that can refer to.
To sum up, technical solution of the present invention had not only had multiple system dynamic simulation function but also can realize the system of parametrization
Cold, air-conditioning device simulation calculation.Therefore, technical solution of the present invention solves the department of computer science of current Refrigeration & Air-Conditioning equipment
System can not achieve the problem of dynamic property prediction of system.
The present invention can push away extensively in fields such as the Computer-aided Design Technologies of Refrigeration & Air-Conditioning equipment based on the above reasons
Extensively.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to do simply to introduce, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without any creative labor, can be with
It obtains other drawings based on these drawings.
Fig. 1 is the parameterisable dynamic simulation computing system block diagram that the present invention is applied to refrigeration or air-conditioning system.
Fig. 2 is the component configuring condition schematic diagram of refrigeration or air-conditioning system described in embodiment 1.
Fig. 3 is the condensing pressure simulation curve figure that simulation calculation obtains.
Fig. 4 is the condensator outlet enthalpy simulation curve figure that simulation calculation obtains.
Fig. 5 is each phase section length simulation curve figure of condenser that simulation calculation obtains.
Fig. 6 is the evaporating pressure simulation curve figure that simulation calculation obtains.
Fig. 7 is the evaporator outlet enthalpy simulation curve figure that simulation calculation obtains.
Fig. 8 is each phase section length simulation curve figure of evaporator that simulation calculation obtains.
Specific embodiment
In order to enable those skilled in the art to better understand the solution of the present invention, below in conjunction in the embodiment of the present invention
Attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only
The embodiment of a part of the invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people
The model that the present invention protects all should belong in member's every other embodiment obtained without making creative work
It encloses.
It should be noted that description and claims of this specification and term " first " in above-mentioned attached drawing, "
Two " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way
Data be interchangeable under appropriate circumstances, so as to the embodiment of the present invention described herein can in addition to illustrating herein or
Sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that cover
Cover it is non-exclusive include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to
Step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, product
Or other step or units that equipment is intrinsic.
Embodiment 1
As shown in Figure 1, the present invention provides a kind of applied to the calculating of the parameterisable dynamic simulation of refrigeration or air-conditioning equipment
System, comprising:
User's input module, for inputting configuration parameter, structural parameters and the external environment phase of refrigeration or air-conditioning equipment
Close input parameter;
Dynamic simulation computing module, when for selecting to solve simulation model the simulation algorithm that needs and be arranged convergence precision,
Simulation step length and simulation time simultaneously use simulation model, carry out simulation calculation according to the parameter inputted in user's input module;
Parameter display module, for showing the simulation result obtained by the dynamic simulation computing module and part
The corresponding dynamic changing curve of simulation result;
System overall control module, for control system according to user's input module, dynamic simulation computing module and parameter
The sequence of display module is successively run.
It further, further include system initialization module, before system is run each time, system overall control module control
System initialization module operation processed, user's input module, dynamic simulation computing module and parameter display module are initialized as writing from memory
Recognize value;System initialization module is initial by each module by the control of system overall control module before system is run each time
Default value is turned to, guarantees that system each run is not influenced by last operational process;
It further include data and curves print module, when parameter display module shows simulation result and dynamic changing curve
Afterwards, the content that system overall control module control data and curves print module operation print parameters display module is shown.
Further, refrigeration or air-conditioning equipment include evaporator, condenser, refrigeration compressor and expansion valve, configuration parameter
Refer to the quantity of refrigeration compressor, condenser, evaporator and expansion valve, and meet following relationship:
(1) quantity of refrigeration compressor is equal to the quantity of condenser;
(2) quantity of expansion valve is equal to the quantity of evaporator;
As shown in Fig. 2, in the present embodiment, the quantity of refrigeration compressor and condenser is set as 2, expansion valve and evaporator
Quantity be set as 5.
Further, the structural parameters of refrigeration compressor include cylinder diameter, number of cylinders and length of stroke;The structure of condenser
Parameter includes bore, pipe outside diameter, length of tube, pipe wall material specific heat at constant pressure and pipe wall material density;The structural parameters of evaporator
Including bore, pipe outside diameter, spacing of fin, fin thickness, effective length, pipe wall material specific heat at constant pressure and pipe wall material density;
The structural parameters of expansion valve include unit refrigerant type, booting electronic valve pretravel value, valve body minimum aperture setting, valve body most
The flow area of big aperture setting and valve body standard-sized sheet is arranged.
Further, external environment correlated inputs parameter refers to the setup parameter and condenser of refrigeration compressor and expansion valve
It is cooling with the state parameter of evaporator periphery medium, including clammy compressor input speed, condenser cooling water flow, condenser
Water goes out temperature, evaporator air inlet mass flow, evaporator inlet air temperature and the practical aperture of electronic valve into temperature, condenser cooling water.
Further, the simulation model of dynamic simulation computing module includes:
(1) Simulation Calculation of evaporator and condenser is indicated by following formula:
Refrigerant quality keeps weighing apparatus equation:
Refrigerant energy keeps weighing apparatus equation:
Tube wall energy conservation equation:
ai: for condenser or the evaporator tube inner wall unit area coefficient of heat transfer, unit is kW/ (m2·K);
ao: for condenser or the evaporator tube outer wall unit area coefficient of heat transfer, unit is kW/ (m2·K);
Di: for condenser or evaporator tube internal diameter, unit m;
Do: for condenser or evaporator tube outer diameter, unit m;
Tb: for the temperature of cooling medium, unit K;
Tr: for the temperature of refrigerant, unit K;
Tt: for condenser or the tube wall temperature of evaporator, unit K;
ρ: for the density of refrigerant, unit kg/m3;
U: being refrigerant along condenser or the flow velocity of evaporator tube length direction, unit m/s;
P: for the refrigerant pressure in condenser or evaporator tube, unit Pa;
H: for the refrigerant specific enthalpy in condenser or evaporator tube, unit kJ/kg;
Cp: for condenser or the specific heat at constant pressure of evaporator tube wall material, unit kJ/kgK;
ρt: for the density of condenser or evaporator tube wall material, unit kg/m3;
At: for the circular ring section of condenser or evaporator tube wall material product, unit m2;
(2) Simulation Calculation of refrigeration compressor is indicated by following formula:
Vt=π D2Szn/240
Td=Ts(Pc/Pe)(n-1)/n
Vt: for the theoretical displacement flow of refrigeration compressor, unit m3/s;
D: for the cylinder bore of refrigeration compressor, unit m;
S: for the piston stroke of refrigeration compressor, unit m;
Z: for the speed of crankshaft of refrigeration compressor, unit r/min;
N: for the number of cylinders of refrigeration compressor;
mc: for the mass flow of refrigeration compressor, unit kg/s;
λ: for the volumetric efficiency of refrigeration compressor;
vs: for the aspiration ratio volume of refrigerant, unit m3/kg;
Td: for the delivery temperature of refrigeration compressor, unit K;
Ts: for the suction temperature of refrigeration compressor, unit K;
N: for the polytropic exponent of refrigerant compression process;
(3) Simulation Calculation of expansion valve is indicated by following formula:
hiv=hov
mv: for the refrigerant flow for flowing through expansion valve, unit kg/s;
Cv: for the discharge coefficient of expansion valve;
Av: for the circulation area of expansion valve, unit m2;
ρ: for the density of expansion valve inlet refrigerant liquid, unit kg/m3;
p1: for the pressure of expansion valve inlet refrigerant, unit Pa;
p2: for the pressure of expansion valve outlet refrigerant, unit Pa;
hiv: for expansion valve import enthalpy, unit kJ/kg;
hov: for expansion valve outlet enthalpy, unit kJ/kg.
Further, dynamic simulation computing module further includes that multiple refrigeration or air-conditioning equipment are established using fluid network method
Pipe net leakage rate, pipe net leakage rate is for constructing Simulation Calculation, to obtain parametric results;Table can be carried out with following formula
It states:
Discharge relation between node and branch
Energy relationship between node and branch
AI: it indicates to flow into matrix, wherein element aijIndicate that node i is the inflow node of branch j, value is 0 or 1;
A0: indicate outflow matrix, wherein element aijIndicate that node i is the outflow node of branch j, value is 0 or 1.
Further, as shown in table 3, the simulation result that dynamic simulation module obtains includes: the system of refrigeration compressor
Cryogen mass flow and delivery temperature;Overheat segment length, two-phase section length, the supercooling segment length, condensing pressure, condensation of condenser
Device exports enthalpy, superheat section tube wall temperature, two-phase section tube wall temperature and super cooled sect tube wall temperature;The two-phase segment length of evaporator,
Overheat segment length, evaporating pressure, evaporator outlet enthalpy, two-phase section tube wall temperature and superheat section tube wall temperature;And expansion valve
Refrigerant mass fluxes;Wherein, condenser simulation result related to evaporator can export as dynamic changing curve,
Such as condensing pressure, condensator outlet enthalpy, each phase section length of condenser, evaporating pressure, evaporator outlet enthalpy, each phase of evaporator
Section length;Wherein, condenser simulation result related to evaporator can export as dynamic changing curve.
Further, in dynamic simulation computing module, selective simulation algorithm includes Euler method, improved Euler method and dragon
Ge Kutafa;Iterative calculation method selected by this example is Fourth order Runge-Kutta;
Simulation step length: user can carry out the setting of simulation step length according to simulated effect, and the present embodiment is set as 0.02 second;
Simulation time: the convergent that user can calculate according to simulation model carries out the setting of simulation time, once reach
Simulation time, simulation calculation process terminate, and the present embodiment is set as 500s;
Convergence precision: the deviation of twice adjacent calculation result in iterative calculation, when deviation is less than or equal to convergence essence
When spending, even if simulation time does not arrive, calculating will also terminate, and the present embodiment is set as 0.0001.
In the present embodiment, the structural parameters of user's input module input and external environment correlated inputs parameter be respectively such as
Shown in Tables 1 and 2:
1 related elements structural parameters of table
2 external environment correlated inputs parameter of table
Title | Unit | Numerical value |
Refrigeration compressor input speed | rpm | 1400 |
Condenser cooling water flow | m3/h | 8.2 |
Condenser cooling water is into temperature | K | 298.15 |
Condenser cooling water goes out temperature | K | 299.65 |
Evaporator enters the wind mass flow | kg/s | 2.138 |
Evaporator inlet air temperature | K | 255.15 |
The practical aperture of electronic valve | 40% |
The results are shown in Table 3 for final simulation calculation:
3 simulation result parameter of table
As shown in figures 3-8, the corresponding condensing pressure of the present embodiment simulation result, condensator outlet enthalpy are respectively indicated
The dynamic changing curve of value, each phase section length of condenser, evaporating pressure, evaporator outlet enthalpy, each phase section length of evaporator.
Further, the relationship under each state between the temperature, specific enthalpy, pressure, density of refrigerant uses Matlab software
Data fitting is carried out to associated refrigeration agent parameter and is obtained.
Further, when the problems such as data scatter, condenser or evaporator partitioned organization mutation occurs in simulation calculation process
When, parameter display module will provide dummy error prompt.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify to technical solution documented by previous embodiment, or some or all of the technical features are equal
Replacement;And these are modified or replaceed, the model for technical solution of the embodiment of the present invention that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (9)
1. it is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system characterized by comprising
User's input module, configuration parameter, structural parameters and the external environment correlation for inputting refrigeration or air-conditioning equipment are defeated
Enter parameter;
Simultaneously convergence precision, emulation is arranged in dynamic simulation computing module, the simulation algorithm needed when for selecting and solve simulation model
Step-length and simulation time simultaneously use simulation model, carry out simulation calculation according to the parameter inputted in user's input module;
Parameter display module, for showing the simulation result and partial simulation that obtain by the dynamic simulation computing module
The corresponding dynamic changing curve of calculated result;
System overall control module is shown for control system according to user's input module, dynamic simulation computing module and parameter
The sequence of module is successively run.
2. the parameterisable dynamic simulation computing system according to claim 1 for being applied to refrigeration or air-conditioning equipment, special
Sign is, further includes system initialization module, and before system is run each time, system overall control module control system is initial
Change module operation, user's input module, dynamic simulation computing module and parameter display module are initialized as default value;Further include
Data and curves print module, after parameter display module shows simulation result and dynamic changing curve, system is totally controlled
The content that molding block control data and curves print module operation print parameters display module is shown.
3. the parameterisable dynamic simulation computing system according to claim 1 for being applied to refrigeration or air-conditioning equipment, special
Sign is that refrigeration or air-conditioning equipment include evaporator, condenser, refrigeration compressor and expansion valve, and configuration parameter refers to refrigeration pressure
The quantity of contracting machine, condenser, evaporator and expansion valve, and meet following relationship:
(1) quantity of refrigeration compressor is equal to the quantity of condenser;
(2) quantity of expansion valve is equal to the quantity of evaporator.
4. the parameterisable dynamic simulation computing system according to claim 3 for being applied to refrigeration or air-conditioning equipment, special
Sign is that the structural parameters of refrigeration compressor include cylinder diameter, number of cylinders and length of stroke;The structural parameters of condenser include pipe
Internal diameter, pipe outside diameter, length of tube, pipe wall material specific heat at constant pressure and pipe wall material density;The structural parameters of evaporator include bore,
Pipe outside diameter, spacing of fin, fin thickness, effective length, pipe wall material specific heat at constant pressure and pipe wall material density;The structure of expansion valve
Parameter include unit refrigerant type, booting electronic valve pretravel value, valve body minimum aperture setting, valve body maximum opening setting and
The flow area of valve body standard-sized sheet is arranged.
5. the parameterisable dynamic simulation computing system according to claim 4 for being applied to refrigeration or air-conditioning equipment, special
Sign is that external environment correlated inputs parameter refers to outside the setup parameter and condenser and evaporator of refrigeration compressor and expansion valve
Enclose the state parameter of medium, including clammy compressor input speed, condenser cooling water flow, condenser cooling water are into temperature, cold
Condenser cooling water goes out temperature, evaporator air inlet mass flow, evaporator inlet air temperature and the practical aperture of electronic valve.
6. the parameterisable dynamic simulation computing system according to claim 3 for being applied to refrigeration or air-conditioning equipment, special
Sign is that the simulation model of dynamic simulation computing module includes:
(1) Simulation Calculation of evaporator and condenser is indicated by following formula:
Refrigerant quality keeps weighing apparatus equation:
Refrigerant energy keeps weighing apparatus equation:
Tube wall energy conservation equation:
ai: for condenser or the evaporator tube inner wall unit area coefficient of heat transfer, unit is kW/ (m2·K);
ao: for condenser or the evaporator tube outer wall unit area coefficient of heat transfer, unit is kW/ (m2·K);
Di: for condenser or evaporator tube internal diameter, unit m;
Do: for condenser or evaporator tube outer diameter, unit m;
Tb: for the temperature of cooling medium, unit K;
Tr: for the temperature of refrigerant, unit K;
Tt: for condenser or the tube wall temperature of evaporator, unit K;
ρ: for the density of refrigerant, unit kg/m3;
U: being refrigerant along condenser or the flow velocity of evaporator tube length direction, unit m/s;
P: for the refrigerant pressure in condenser or evaporator tube, unit Pa;
H: for the refrigerant specific enthalpy in condenser or evaporator tube, unit kJ/kg;
Cp: for condenser or the specific heat at constant pressure of evaporator tube wall material, unit kJ/kgK;
ρt: for the density of condenser or evaporator tube wall material, unit kg/m3;
At: for the circular ring section of condenser or evaporator tube wall material product, unit m2;
(2) Simulation Calculation of refrigeration compressor is indicated by following formula:
Vt=π D2Szn/240
Td=Ts(Pc/Pe)(n-1)/n
Vt: for the theoretical displacement flow of refrigeration compressor, unit m3/s;
D: for the cylinder bore of refrigeration compressor, unit m;
S: for the piston stroke of refrigeration compressor, unit m;
Z: for the speed of crankshaft of refrigeration compressor, unit r/min;
N: for the number of cylinders of refrigeration compressor;
mc: for the mass flow of refrigeration compressor, unit kg/s;
λ: for the volumetric efficiency of refrigeration compressor;
vs: for the aspiration ratio volume of refrigerant, unit m3/kg;
Td: for the delivery temperature of refrigeration compressor, unit K;
Ts: for the suction temperature of refrigeration compressor, unit K;
N: for the polytropic exponent of refrigerant compression process;
(3) Simulation Calculation of expansion valve is indicated by following formula:
hiv=hov
mv: for the refrigerant flow for flowing through expansion valve, unit kg/s;
Cv: for the discharge coefficient of expansion valve;
Av: for the circulation area of expansion valve, unit m2;
ρ: for the density of expansion valve inlet refrigerant liquid, unit kg/m3;
p1: for the pressure of expansion valve inlet refrigerant, unit Pa;
p2: for the pressure of expansion valve outlet refrigerant, unit Pa;
hiv: for expansion valve import enthalpy, unit kJ/kg;
hov: for expansion valve outlet enthalpy, unit kJ/kg.
7. the parameterisable dynamic simulation computing system according to claim 6 for being applied to refrigeration or air-conditioning equipment, special
Sign is that dynamic simulation computing module further includes the pipe network mould that multiple refrigeration or air-conditioning equipment are established using fluid network method
Type can be stated with following formula:
Discharge relation between node and branch
Energy relationship between node and branch
AI: it indicates to flow into matrix, wherein element aijIndicate that node i is the inflow node of branch j, value is 0 or 1;
A0: indicate outflow matrix, wherein element aijIndicate that node i is the outflow node of branch j, value is 0 or 1.
8. the parameterisable dynamic simulation computing system according to claim 7 for being applied to refrigeration or air-conditioning equipment, special
Sign is that the simulation result that dynamic simulation module obtains includes:
The refrigerant mass fluxes and delivery temperature of refrigeration compressor;
Overheat segment length, two-phase section length, the supercooling segment length, condensing pressure, condensator outlet enthalpy, superheat section pipe of condenser
Wall temperature, two-phase section tube wall temperature and super cooled sect tube wall temperature;
Two-phase segment length, overheat segment length, evaporating pressure, evaporator outlet enthalpy, two-phase section tube wall temperature and the mistake of evaporator
Hot arc tube wall temperature;
And the refrigerant mass fluxes of expansion valve.
9. the parameterisable dynamic simulation computing system according to claim 8 for being applied to refrigeration or air-conditioning equipment, special
Sign is, in dynamic simulation computing module, selective simulation algorithm includes Euler method, improved Euler method and runge kutta method;
The value of convergence precision is 0.0001;The value of simulation step length is 0.02 second;The value of simulation time is 500s.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910447969.6A CN110175403A (en) | 2019-05-27 | 2019-05-27 | It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910447969.6A CN110175403A (en) | 2019-05-27 | 2019-05-27 | It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110175403A true CN110175403A (en) | 2019-08-27 |
Family
ID=67696321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910447969.6A Pending CN110175403A (en) | 2019-05-27 | 2019-05-27 | It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110175403A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110641250A (en) * | 2019-11-05 | 2020-01-03 | 重庆大学 | Intelligent control method of air conditioning system of electric automobile based on human body thermal comfort theory and fuzzy PID control |
CN111125938A (en) * | 2020-01-15 | 2020-05-08 | 华南理工大学 | Suboptimal algorithm-based optimization design method for large central air-conditioning chilled water pipe network |
CN113759749A (en) * | 2020-06-02 | 2021-12-07 | 青岛海信日立空调系统有限公司 | Water chilling unit simulation system |
CN114282389A (en) * | 2021-12-30 | 2022-04-05 | 西安交通大学 | MATLAB-based centrifugal fan system dynamic simulation method and system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005180267A (en) * | 2003-12-18 | 2005-07-07 | Mitsubishi Heavy Ind Ltd | Turbo refrigerator, compressor thereof, and control method thereof |
CN106529021A (en) * | 2016-11-09 | 2017-03-22 | 东南大学 | Air conditioning system simulation method based on feature recognition |
-
2019
- 2019-05-27 CN CN201910447969.6A patent/CN110175403A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005180267A (en) * | 2003-12-18 | 2005-07-07 | Mitsubishi Heavy Ind Ltd | Turbo refrigerator, compressor thereof, and control method thereof |
CN106529021A (en) * | 2016-11-09 | 2017-03-22 | 东南大学 | Air conditioning system simulation method based on feature recognition |
Non-Patent Citations (1)
Title |
---|
杨杰等: "船舶多联式制冷系统的动态建模与仿真", 《船舶工程》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110641250A (en) * | 2019-11-05 | 2020-01-03 | 重庆大学 | Intelligent control method of air conditioning system of electric automobile based on human body thermal comfort theory and fuzzy PID control |
CN111125938A (en) * | 2020-01-15 | 2020-05-08 | 华南理工大学 | Suboptimal algorithm-based optimization design method for large central air-conditioning chilled water pipe network |
CN111125938B (en) * | 2020-01-15 | 2021-07-16 | 华南理工大学 | Suboptimal algorithm-based optimization design method for large central air-conditioning chilled water pipe network |
CN113759749A (en) * | 2020-06-02 | 2021-12-07 | 青岛海信日立空调系统有限公司 | Water chilling unit simulation system |
CN113759749B (en) * | 2020-06-02 | 2024-05-03 | 青岛海信日立空调系统有限公司 | Water chilling unit simulation system |
CN114282389A (en) * | 2021-12-30 | 2022-04-05 | 西安交通大学 | MATLAB-based centrifugal fan system dynamic simulation method and system |
CN114282389B (en) * | 2021-12-30 | 2024-01-16 | 西安交通大学 | Centrifugal fan system dynamic simulation method and system based on MATLAB |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110175403A (en) | It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system | |
Wang et al. | Optimization of refrigeration system with gas-injected scroll compressor | |
Bodys et al. | Full-scale multi-ejector module for a carbon dioxide supermarket refrigeration system: Numerical study of performance evaluation | |
CN102779217B (en) | Computer simulation performance computation method of refrigeration system under frosting working condition | |
Wang et al. | Numerical analysis on the effects of refrigerant injection on the scroll compressor | |
Zhao et al. | Study of the performance of an urban original source heat pump system | |
Wang et al. | Numerical research on the scroll compressor with refrigeration injection | |
Rigola et al. | Numerical simulation and experimental validation of internal heat exchanger influence on CO2 trans-critical cycle performance | |
Yan et al. | Optimization on ejector key geometries of a two-stage ejector-based multi-evaporator refrigeration system | |
Wang et al. | Performance improvement of air source heat pump using gas-injected rotary compressor through port on blade | |
Fu et al. | Steady-state simulation of screw liquid chillers | |
Leducq et al. | Low order dynamic model of a vapor compression cycle for process control design | |
Fu et al. | Dynamic simulation of air-to-water dual-mode heat pump with screw compressor | |
Castaing-Lasvignottes et al. | Dynamic simulation of reciprocating refrigeration compressors and experimental validation | |
Liu et al. | Model-based dynamic optimal control of a CO2 heat pump coupled with hot and cold thermal storages | |
Jin et al. | Intermediate pressure of two-stage compression system under different conditions based on compressor coupling model | |
Rasmussen et al. | Review of compressor models and performance characterizing variables | |
Wei et al. | Frosting performance variations of variable-frequency air source heat pump in different climatic regions | |
Kong et al. | Mass flow rate prediction of direct-expansion solar-assisted heat pump using R290 based on ANN model | |
Zhang et al. | The model predictive control strategy of the transcritical CO2 air conditioning system used in railway vehicles | |
Maddah et al. | Determination of the optimal discharge pressure of the transcritical CO 2 heat pump cycles for heating and cooling performances based on new correlation | |
Wang et al. | Comparative research on air conditioner with gas-injected rotary compressor through injection port on blade | |
Ohkura et al. | Numerical analysis on performance enhancement of a CO2 heat pump water heating system by extracting tepid water | |
Zhang et al. | Performance investigation of a novel EEV-based ejector for refrigerator-freezers | |
Li | Optimal analysis of gas cooler and intercooler for two-stage CO2 trans-critical refrigeration system |
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 |