IJESR/Sept 2012/ Volume-2/Issue-9/Article No-1/ ISSN International Journal of Engineering & Science Research

Similar documents
EEE 241: Linear Systems

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Parameter Estimation for Dynamic System using Unscented Kalman filter

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Numerical Heat and Mass Transfer

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

DUE: WEDS FEB 21ST 2018

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

A new Approach for Solving Linear Ordinary Differential Equations

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL

Global Sensitivity. Tuesday 20 th February, 2018

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Time-Varying Systems and Computations Lecture 6

1 Derivation of Point-to-Plane Minimization

Note 10. Modeling and Simulation of Dynamic Systems

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Neuro-Adaptive Design - I:

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Scroll Generation with Inductorless Chua s Circuit and Wien Bridge Oscillator

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Identification of Linear Partial Difference Equations with Constant Coefficients

Automatic PID Controller Tuning for Robots with Nonlinear Friction at the Joints

NUMERICAL DIFFERENTIATION

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Chapter Newton s Method

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

Distributed Exponential Formation Control of Multiple Wheeled Mobile Robots

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Kernel Methods and SVMs Extension

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

PREDICTIVE CONTROL BY DISTRIBUTED PARAMETER SYSTEMS BLOCKSET FOR MATLAB & SIMULINK

RELIABILITY ASSESSMENT

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Linear Approximation with Regularization and Moving Least Squares

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Chapter 13: Multiple Regression

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

Tracking with Kalman Filter

829. An adaptive method for inertia force identification in cantilever under moving mass

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

Odd/Even Scroll Generation with Inductorless Chua s and Wien Bridge Oscillator Circuits

Simulation Study of Thermal Sensation Based Control for Single Family Type Distributed Heating Systems Peizhang Chen 1, Fulin Wang 1,*, and Zheliang C

Uncertainty in measurements of power and energy on power networks

Appendix B: Resampling Algorithms

A Fast Computer Aided Design Method for Filters

A Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Finite Element Modelling of truss/cable structures

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

( ) = ( ) + ( 0) ) ( )

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

Comparison of Wiener Filter solution by SVD with decompositions QR and QLP

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Numerical Transient Heat Conduction Experiment

Research Article Green s Theorem for Sign Data

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

1 Convex Optimization

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Air Age Equation Parameterized by Ventilation Grouped Time WU Wen-zhong

Assignment 4. Adsorption Isotherms

Solution Thermodynamics

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Numerical Solution of Ordinary Differential Equations

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Calculation of time complexity (3%)

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Clock-Gating and Its Application to Low Power Design of Sequential Circuits

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

The Chaotic Robot Prediction by Neuro Fuzzy Algorithm (2) = θ (3) = ω. Asin. A v. Mana Tarjoman, Shaghayegh Zarei

Constitutive Modelling of Superplastic AA-5083

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields

Operating conditions of a mine fan under conditions of variable resistance

Color Rendering Uncertainty

Neuro-Adaptive Design II:

Adaptive Consensus Control of Multi-Agent Systems with Large Uncertainty and Time Delays *

Supporting Information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Lecture 12: Discrete Laplacian

Pulse Coded Modulation

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica

Module 3: Element Properties Lecture 1: Natural Coordinates

Canonical transformations

Lecture 21: Numerical methods for pricing American type derivatives

PHYS 705: Classical Mechanics. Calculus of Variations II

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Transcription:

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 Internatonal Journal of Engneerng & Scence Research ADAPIVE SELF-UNING DECOUPLED CONROL OF EMPERAURE AND RELAIVE HUMIDIY FOR A HVAC SYSEM Ahmad Parvaresh*, Seyed Mohammad Al Mohammad 2 Department of electrcal engneerng of Shahd Bahonar Unversty of Kerman, Kerman. 2 Department of electrcal engneerng of Shahd Bahonar Unversty of Kerman, Kerman. ABSRAC hs paper presents an adaptve self-tunng decoupled controller for control of temperature and relatve humdty of a heatng, ventlatng, and ar condtonng (HVAC) system n order to acheve comfort condtons for occupants. It s mportant to be able to control temperature and relatve humdty ndependently. Because of nteracton of these varables n some ndustral processes, t s necessary to use decouplng technques to acheve accurate control. In ths regards, we have used recursve least square (RLS) algorthm for dentfcaton of the HVAC system to self-tunng of proportonal-ntegral-dervatve (PID) controller for temperature and relatve humdty. Smulaton results demonstrate that the adaptve selftunng decoupled controller s applcable for ths HVAC system. Keywords: Adaptve control, RLS algorthm, HVAC system, PID controller, self-tunng decoupled control.. INRODUCION Heatng, ventlatng, and ar-condtonng (HVAC) systems requre control of envronmental varables such as temperature, relatve humdty, pressure etc. the most of the controllers that have used n HVAC systems are of Proportonal-Integral-Dervatve (PID) type. HVAC systems used n lvng and ndustral buldngs should fulfll thermal comfort needs and ndoor ar qualty. he comfort of the people n ther lvng ndoor envronment s partally dependent on some factors, as qualty, humdty, pressure, and temperature of ar nsde the buldng. he purpose of the HVAC system of a buldng s to provde complete thermal comfort for ts occupants. Hence, t s necessary to understand the thermal and humdty aspects of the human body n order to desgn an effectve HVAC system. HVAC system s a non-lnear and tme varant system. It s dffcult to acheve desred trackng control performance. Over the past few years, several methods for determnng of parameters of appled controllers on these systems have been developed. As ndcated n [], more than 95% of the control loops are of PID type n process control. Over the years, there are many formulas derved to tune the PID controllers for stable processes, such as Zegler-Nchols, Cohen-Coon, nternal model control, ntegral absolute error optmum (ISE, IAE, and IAE), and recently proposed tunng methods[2,3]. Some researchers have been focused on decoupled control of temperature and relatve humdty of HVAC systems. For nstance, [4-6]. Some studes as [7,8] are about applcaton of PID controller n the HVAC systems. We know HVAC systems are non-lnear and tme *Correspondng Author www.jesr.org 900

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 varant, so, drect controls of these systems are dffcult. Identfcaton approaches are very useful methods for acheve approxmate models. In ths regards, some of researchers use dfferent dentfcaton methods for apply on ths systems as [9]. RLS dentfcaton s one of these methods. In ths paper we developed and smulated an adaptve self-tunng algorthm for tunng PID parameters for a decoupled HVAC system n order to control of temperature and relatve humdty of a thermal zone, where thermal comfort of the system s evaluated. hs paper organzed as follow: In the frst part of ths paper the HVAC system and thermal space model s descrbed. he second part descrbed the non-nteractng control n a HVAC system. After descrbe of PID controller n a separate secton adaptve self-tunng algorthm are presented. In the last part smulaton result are presented. 2. HE HVAC SYSEM AND HERMAL SPACE MODEL We consder the sngle-zone HVAC system shown n Fgure. Components of ths system nclude thermal space, heatng/coolng col, humdfer/dehumdfer, mxng box, ar flter, supply and return fans, flters, dampers, and ductwork. In ths system, ntally, fresh ar enters and mxes wth 75% of the return ar, and remanng ar s exhausted. hen, mxed ar passes through the heat exchanger components and fnally by supply fan enter to thermal zone. Supply ar satsfes the thermal space. By changng of thermal load, the system controller smultaneously vares volumetrc flow rate of ar and water, so that the desred setponts n temperature and relatve humdty are mantaned. he dfferental equatons descrbng the dynamc behavor of the HVAC system n Fgure, can be conclude from energy-mass equatons as follow. dw3 W2 W3 M z = + () dt V ρ V z a z d2 ( ) f 0.25( ) f = + dt V V 3 2 a o 3 a he h (0.25W + 0.75 W W ) f ρ h f - C V ρc V w o 3 2 a w w he pa he pa he (2) 3 ( 2 ( 3 ) hfg W2 W3 ) fa ( Qz hfgm z ) a d f = + (3) dt V C V ρv C z pa z z pa Copyrght 202 Publshed by IJESR. All rghts reserved 90

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 Fgure. HVAC system model 3. NON-INERACIVE CONROL In order to control of mult-nput mult-output (MIMO) systems, we can use many dfferent methods. One of these methods s non-nteractve control method. In ths method the feedback s used transform the MIMO system, from the nput-output pont of vew, to an aggregate of ndependent sngle nput SISO channels. Assume a system n state-spacee form wth m nput u and n output y j that = =,2,,m and j=,2,,n. for ths system we can wrte x& = f ( x) + g ( x) u = y= h ( x) n m (4) he problem s to fnd a feedbackk control law as follow m u ψ ( x) δ v = + j j j= (5) Such that for the closed-loop system x& = f ( x ) + g ( x) ψ ( x) + ( g ( x) δ ( x) v ) y = h ( x) m m m j j = j= = (6) he resultng closed-loop systemm s a decoupled system n whch the th outputt s controlled only by th nput. he non-nteractng control law s gven by u A = ( x) b( x) + A ( x) v (7) Copyrght 202 Publshed by IJESR. All rghts reserved 902

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 Where A(x) s decouplng matrx, for defnton and calculaton see [4], v s the new nput for control of control varables ndependently, and b(x) s defned as m h ( f ( x )) = x f ( x) x m h 2 ( f ( x )) = x b( x) = f ( x) x M m h n ( f ( x )) = x f ( x) x For presented HVAC system we want to control temperature and relatve humdty of the zone. In ths regards, we need two control loops and two control sgnals. After calculaton of A(x) and b(x) we can obtan the non-nteractng control sgnals as followng expresson (8) u u = 2 ( αx4β4g( x) +.388 αx4β4g3 ( x)) α x g ( x).388α x β v B ( x) 4 4 4 g3( x) g( x) v2 B2 ( x) (9) x = W, x =, x =, x = fa And Where 3 2 2 3 3 4 g = (5000α W + 07α α x.388 α x +.388α W 2 2 2 2 2.776 0 α x 3.85 0 α x +29.7 0 α ) 4 4 3 2 2 3 2 g = (.5 0 β x β x + ( β β 2. 0 β ) x 4 4 2 3 2 2 3 3 + β 0.25β W + β W + 0.06 β ) 2 o 3 o 3 2 3 g ( x) = ( α ( W 0.024)) + 2 0 4 α x + α x + (2.776 3 2 2 2 2 0 4 α ) x 2 3 α = 60 / V α = 60 h / C V α = / ρ C V α = / ρ V z 2 fg pa z 3 a pa z 4 a z β = 60 / V β = 5 / V β =60 h / C V β = 6000 / ρ C V he 2 he 3 w pa he 4 a pa he B ( x) = ( α x 2.776 0 α x )( f ( x) + g ( x) x ) 4 4 2 4 4 +(.388 α x ) g ( x) x + ( 3.85 0 α x ) 4 4 2 4 2 4 ( f ( x) + g ( x) x ) + ( x g ( x)) 3 3 4 4 B ( x) = (2 0 α x )( f ( x) + g ( x) x ) + ( α x )( g ( x) x ) 4 2 2 4 4 4 2 4 + + + 4 (2.776 0 α2 α)( f3( x) g3( x) x4 ) x4 ( g3( x) x4) (0) () (2) (3) (4) (5) Copyrght 202 Publshed by IJESR. All rghts reserved 903

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 In Eq. 4 and Eq. 5 f ( x) and f ( x ) are 3 f ( x) = 5000α M.388 α ( Q h M ) 4 z 3 z fg z f ( x) = α ( Q h M ) 3 3 z fg z (6) 4. PID CONROLLER A PID controller s a knd of lnear controller, as t composes the control error accordng to the settng value and process output and then makes the controller output value based on the lnear resultant of error s proporton, ntegral and dervatve. For a PID controller, the control sgnal at tme t s determned from Where u( t) s the controller output, parameter, t ( ) ( ) ( ) de( t) (7) u t = K pe t + K e t dt+ K 0 d dt K p s the proporton parameter, K s the ntegral K d s dervatve parameter, e( t ) s the error sgnal at tme t defned as e( t) = ysp ( t) y( t), where y SP s the set pont for process output at tme t and y( t ) s process output at tme t. 5. ADAPIVE SELF-UNING CONROL ALGORIHM hs s a tradtonal adaptve flter algorthm that s useful for dentfcaton of control systems. It estmates the current parameter vector ˆ( θ k) based on the prevous estmated vector ˆ( θ k ), as follows [5,6]: ˆ θ ( k) = f ( ˆ θ ( k ), D( k), k) (8) Where, D(k) denotes data avalable at tme (k), and f(.,.,.) denotes an algebrac functon, the form of whch determnes the specfc algorthm. In the case of dynamc system, data D(k) normally consder the form of present and past observaton of the system outputs and nputs. For mult-parameter system, ths form can be represented as follows: Where, y( k) = ψ ( k) θ (9) ψ ( k) = [ y( k ),..., y( k m), u( k),..., u( k m)] (20) = [ a,..., a, b,..., b ] (2) θ m m+ he estmaton of the parameters vectorθ s performed n a way such that the estmated ˆr θ mnmzes the cost ndex J (r) where r denotes the number of sets of measurement, r ψ ˆ θ 2 (22) k= J ( r) = ( y( k) ( k) ( k)) Equaton (5) can be wrtten n a recursve form as: Copyrght 202 Publshed by IJESR. All rghts reserved 904

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 ˆ θ ( k) = ˆ θ ( k ) + P( k) ψ ( k)( y( k) ψ ( k) ˆ θ ( k )) P( k ) u( k) ψ ( k) P( k ) P( k ) = [ P( k ) ] + ψ ( k) P( k ) ψ ( k) (23) (24) We use from ths algorthm for dentfcaton of HVAC system ntroducedd n prevous secton. Fgure 2 shows the dagram sketch of RLS-PID control of decoupled HVAC system. Fg. 2. Dagram sketch of RLS-PID control of decoupled HVAC system 6. SIMULAION RESULS For nstance, we consder a system by nonlnear equaton as follow 2 y ( t ) y( t) = + u( t) 2 + y ( t ) (25) Where u( t) = 0.4sn(2 πt / 50) for t=0,,2,, 00 and u( t) = 0.3sn(2 πt / 50) + 0.sn(2 πt / 5) for t= 0, 02,, 200. Fgure (3) shows the proposed dentfcaton result of the system that ntroduced n Eq. 25. he normalzed error of the dentfcaton s 0.29 that s a performance measure defned as E ( t) = ( e e ) / ( y y ), where s a transpose of the vector. n After dentfcaton model s completed, the trackng operaton takes command of the RLS process control to track the desred set-ponts r ( t) = 25 u( t) + 50 u( t 200), that t= 0,,, 200 for emperature and r. ( t ) = 0 u( t) + 20 u( t 00), that t=0,,, 2000 for Relatve R H Humdty. he results of trackng performance are llustrated n fgure (4) and fgure (5) wth the mean square error (MSE) 0.343 for the plant and for RLS self tunng algorthm response of 0.023. Copyrght 202 Publshed by IJESR. All rghts reserved 905

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 0.8 Output Plant Identfcaton 0.6 Ampltude 0.4 0.2 0-0.2-0.4-0.6 0 50 00 50 200 250 me(s) Fg. 3. RLS algorthm dentfcaton of a nonlnear system 350 300 Set pont emperature RLS algorthm 250 emperature(c) 200 50 00 50 0-50 0 20 40 60 80 00 20 40 60 80 200 me(s) Fg. 4. Output results of trackng operaton for emperature Copyrght 202 Publshed by IJESR. All rghts reserved 906

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 40 20 Set Pont Relatve Humdty RLS algorthm 00 Relatve Humdty 80 60 40 20 0-20 0 20 40 60 80 00 20 40 60 80 200 me(s) Fg. 5. Output results of trackng operaton for Relatve Humdty 7. CONCLUSION AND RECOMMENDAIONS In ths paper, the performance of a HVAC system that uses a nonlnear decouplng algorthm for control of temperature and relatve humdty s studed. We show that decouplng control law works satsfactorly. We used RLS algorthm for dentfcaton of the system and results of the dentfcaton are used for tunng of PID controller parameters. Smulaton results were used to demonstrate of the proposed adaptve self-tunng decoupled control. REFERENCES [] Astrom KJ, Hagglund. PID Controllers: heory, Desgn, and unng, Instrument Socety of Amerca. Research rangle Park, NC, 995. [2] Wang YG, Shao HH. PID autotuner based on gan and phase margn specfcatons. Ind.Eng.Chem.Res. 999; 38: 3007-302. [3] Wang YG, Shao HH. Optmal tunng for PI controller. Automatca 2000; 36: 47-52. [4] Rentel-Gomez C, Velez-Reyes M. Decoupled control of temperature and relatve humdty usng a varable ar volume HVAC system and non-nteractng control. In: Proceedng of the IEEE nternatonal conference on control applcatons, Mexco Cty, Mexco; 200.p.47 5. [5] Rentel CH. Nonlnear decouplng of temperature and relatve humdty for an ar condtonng system. MS. hess, Unversty of Puerto Rco, Mayaguez Campus; December 996. [6] Cu J, Watanabe, Ryu Y, Akash Y, Nshyama N. Numercal smulaton on smultaneous control process of ndoor ar temperature and humdty. In: Sxth nternatonal IBPSA conference, proceedngs, 999; 2; 005 2. Copyrght 202 Publshed by IJESR. All rghts reserved 907

IJESR/Sept 202/ Volume-2/Issue-9/Artcle No-/900-908 ISSN 2277-2685 [7] B Q, Wen-Jan C. Advanced controller auto-tunng and ts applcaton n HVAC systems. Control Engneerng Practce 2000; 8: 633-644. [8] Guang Q, Zaheeruddn M. Real-tme tunng of PI controllers n HVAC systems. Int J Energy Res 2004; 28: 33 27. [9] B Q, Ca WJ, Lee EL, Wang QG, Hang CC, Zhang Y. Robust dentfcaton of frstorder plus dead-tme model from step response. Control Engneerng Practce, 999; 7: 7-77. Copyrght 202 Publshed by IJESR. All rghts reserved 908