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 PDF

Info

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
Application number
CN201910447969.6A
Other languages
Chinese (zh)
Inventor
杨杰
高洪涛
邢鹏飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian Maritime University
Original Assignee
Dalian Maritime University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian Maritime University filed Critical Dalian Maritime University
Priority to CN201910447969.6A priority Critical patent/CN110175403A/en
Publication of CN110175403A publication Critical patent/CN110175403A/en
Pending legal-status Critical Current

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F25REFRIGERATION OR COOLING; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS; MANUFACTURE OR STORAGE OF ICE; LIQUEFACTION SOLIDIFICATION OF GASES
    • F25BREFRIGERATION MACHINES, PLANTS OR SYSTEMS; COMBINED HEATING AND REFRIGERATION SYSTEMS; HEAT PUMP SYSTEMS
    • F25B49/00Arrangement or mounting of control or safety devices
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting 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

It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system
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.
CN201910447969.6A 2019-05-27 2019-05-27 It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system Pending CN110175403A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
杨杰等: "船舶多联式制冷系统的动态建模与仿真", 《船舶工程》 *

Cited By (7)

* Cited by examiner, † Cited by third party
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