Nonlinear Modeling of Butterfly Valves and Flow Rate Control Using the Circle Criterion Bode Plot
|
|
- Calvin Preston
- 5 years ago
- Views:
Transcription
1 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, WeC.6 Nonlinear Modeling of Butterfly Valves and Flow Rate Control Using the Circle Criterion Bode Plot J.D. Taylor, Bruno Sinopoli, and William Messner Abstract Butterfly valves are widely used control components in piping systems. This paper derives the nonlinear mathematical relationship between flow rate, pressure difference, and valve angle in a butterfly valve. We use the Circle Criterion Bode plot and loop shaping to design a flow rate controller in a motor-gearbox-butterfly valve system. The design employs a complex proportional-integral-lead (CPIL) compensator, and we show the benefit of using a damping ratio less than unity. Finally, we demonstrate the effectiveness of the design through simulation. I. INTRODUCTION The distribution and control of water has been a significant engineering challenge throughout history. Notable examples include the irrigation techniques developed in ancient Egypt, the aqueducts constructed by the Roman Empire, and the water-wheel and steam powered mechanisms of the industrial revolution. Water is still used extensively in modern engineering systems. For instance, chilled water networks are often used to cool critical electronics in large-scale data centers and aboard naval ships. The robustness of these systems is of particular importance since loss of cooling water can lead to device failure in a matter of minutes. Rapid rerouting of water flow in piping systems subject to damage must consider the coupled hydro-mechanical dynamics of the system components. In this paper we consider the modeling and control of butterfly valves, which are widely used in chilled water systems. This paper presents a new derivation of the model of flow past a butterfly valve as a function of valve angle and pressure difference. It combines the mathematically formidable free streamline theory [] with the empirical relations for viscous pressure losses in expanding flows (Borda-Carnot) and within pipes (Darcy-Weisbach). The paper also presents an application of the Circle Criterion Bode plot (CBode plot) [] for controlling a motor-gearbox-butterfly valve system with flow rate feedback, and it shows the utility of the complexproportional-integral-lead (CPIL) compensator for this type of loop shaping design. II. VALVE MODELING A schematic of the butterfly valve used to derive the nonlinear relationship between volumetric flow rate, pressure difference, and valve plate angle of attack, α, is presented in Fig.. This derivation is carried out in -D; however, it has been verified experimentally that this approach can The authors are with the Departments of Mechanical, Electrical and Computer Engineering, Carnegie Mellon University, Pittsburgh, PA 53 jdtaylor@andrew.cmu.edu,brunos@ece.cmu.edu, bmessner@andrew.cmu.edu D Stagnation Point V Fig.. d dc α Vj θ L Wake Region d Vj dc Schematic of Butterfly Valve accurately predict the fluid resistance of 3-D valves by integrating over cross-sections [3]. The Borda-Carnot equation, Eq., is an empirical relation describing viscous pressure losses in a fluid undergoing an abrupt expansion in flow area. These losses are due primarily to turbulence in the wake regions. V D P = ξ ρ( V ) () where P is the pressure lost due to the flow expansion and ρ is the fluid density. V = V j V is the reduction in average flow velocity between the vena contracta, V j, and valve outlet, V. ξ is an empirical, dimensionless loss coefficient with a value ranging from ξ. For an abrupt and wide expansion (assumed for the butterfly valve) the loss coefficient is equal to. Note that upstream of the vena contracta, the flow accelerates and pressure losses are negligible in this region. Bernoulli s equation can be modified to account for viscous losses in incompressible flows, Eq.. P + ρv = P + ρv + P () It is convenient to define the pressure loss coefficient, ζ. ζ = P ρv (3) In this analysis, the inlet pressure and flow velocity, P and V, are assumed to be uniform and known. The flow is incompressible, and the velocities at all the cross-sections are related by continuity, Eq. 4. D V =d V =(d c + d c )V j = D V (4) where D,D and V,V are the diameters and average velocities of the inlet and outlet respectively, d and V are the width and average velocity of the orifices between the //$6. AACC 967
2 valve plate and wall. d is maximum when the valve is fully open (α = o ) and minimum when the valve is fully closed (α =9 o ). d c and d c are the widths of the upper and lower jets at their minimum points (i.e. vena contracta) and V j is the corresponding jet velocity. Following [4] we take the jet velocities, V j, to be equal. However, the values of d c and d c are distinct due to the different angles at which the fluid impinges on the surface of the valve plate. We assume in this case that the inlet and outlet diameters are equal D = D D implying V = V V. Assuming that the valve plate length is equal to the diameter of the valve, the flow cross-section as a function of the angle of attack, α, is given by d = D ( sinα) (5) We now define the coefficients of contraction, C c and C c, to be the ratios of the minimum width of the jets to the width of the orifices. C c d c d C c d c d (6) Now from Eqs. 4-6 we have ( ) Cc + C c V = ( sinα)v j (7) The pressure loss coefficient due to the flow expansion, ζ e, from the vena contracta to the outlet of the butterfly valve is given by Eq. 8. [( ζ e = (C c + C c )( sinα) ) ] (8) This coefficient characterizes the pressure loss due to flow expansion. However, at small angles when the valve is nearly open and the wake region is very small, C c and C c approach unity, and the predicted pressure loss is zero. This result is unphysical, and an additional term due to pipe friction must be included. This pressure drop is described by the well known Darcy-Weisbach equation, Eq. 9. P = f L D ρv (9) here f is a dimensionless laminar or turbulent friction factor, L and D are the valve length and diameter respectively, and ρ and V are the density and velocity of the fluid. The friction factor may in general be determined from the Moody chart; however, for laminar flows it is given by the simple formula f =64/Re in which the Reynolds number is Re = ρv D/µ. Therefore, the pressure loss coefficient due to laminar pipe friction, ζ f, is readily determined to be ζ f = 64µL ρv D () The fluid resistance, R, is defined as the ratio of the overall static pressure difference, P P across the valve to the mass flow rate through the valve, ṁ. Since the inlet and outlet diameters are equal, there is no convective acceleration, and the static pressure difference is equal to the total pressure lost due to both flow expansion and pipe friction, P P = P. Note that the fluid resistance is a nonlinear function of both the valve attack angle and the fluid velocity. R = P ṁ =(ζ e + ζ f ) V D () In order to evaluate Eq. 9, the coefficients of contraction for the particular system must be specified. These parameters usually are determined empirically since analytic solutions are generally not available. Fortunately, however, the coefficients of contraction for butterfly valves have both been derived analytically and verified experimentally, [], [4], [5], [6]. These analytical derivations use Helmholtz free streamline theory to calculate the velocity field using successive conformal transformations. Though free streamline theory applies only to inviscid incompressible flow, it is well suited to this problem due to the significance of pressure and inertia in establishing the flow field, while shear and gravitational effects are secondary [5]. In this analysis, the upper and lower coefficients of contraction derived by Sarpkaya in [4] will be used. The interpolated coefficients are shown graphically in Fig.. Coefficients of Contraction C c C c Angle of Attack, α [deg] Fig.. Theoretical upper and lower contraction coefficients versus angle of attack for butterfly valves [4]. The fluid resistance is shown versus angle of attack in Fig. 3. The resistance is small but nonzero for small α due to the pipe friction, it grows exponentially in the approximate range 5 α 75, and it diverges at α =9 at which point the flow ceases. III. CIRCLE CRITERION BODE PLOT The well known Circle Criterion provides a sufficient condition for the stability of a nonlinear feedback system of the form given in Fig. 4, where L(s) is the transfer function of a linear time invariant (LTI) system and Ψ( ) is a memoryless and possibly time-varying sector nonlinearity [7]. The nonlinearity is said to belong to the sector [α, β], if for any input u to Ψ( ), αu uψ(u) βu. 968
3 Fluid Resistance, R [Pa (kg/s) ] Angle of Attack, α [deg] Fig. 3. Fluid resistance as a function of the angle of attack at constant velocity for butterfly valve. + - L( s ) Ψ( ) βu Ψ(u) Fig. 4. Block diagram of system with linear part L(s) and memory-less nonlinearity, Ψ( ). An example of a nonlinear function belonging to the sector [α, β]. Circle Criterion: For a SISO transfer function L(s) and a memory-less function Ψ( ) belonging to the sector [α, β], sufficient conditions for determining the stability of the closed-loop system are the following []. ) β>α, a) The Nyquist plot of L(jω) remains outside the disc with center on the real axis and intersecting the real axis at β and α for α>. b) If L(s) has N u unstable poles, the Nyquist plot of L(jω) encircles the disc anti-clockwise N u times. c) If α =, the Nyquist plot of L(jω) remains to the right of the vertical line x = β. ) β> >α, The Nyquist plot of L(jω) remains inside the disc with center on the real axis and intersecting the real axis at β and α. 3) α<β, Replace L(jω) with L(jω) and Ψ( ) with Ψ( ) which belongs to sector [ β, α]. Return to Case. Fig. shows the Nyquist plot and disc of a system satisfying Case of the Circle Criterion. The Circle Criterion, however, is difficult to use for loopshaping controller design because frequency is a hidden variable on the Nyquist plot. To address this difficulty, Messner and Xia developed the Circle Criterion Bode (CBode) plot αu u []. The CBode plot maps the disc from the Nyquist plot onto the Bode magnitude and phase plots by transforming the Circle Criterion into inequality relationships between the open-loop gain and phase that involve α and β. For a first quadrant nonlinearity (Case ), the inequality is ( Re [L(jω)] ( )) + () > ( α + β ( α β )) () That is, the distance of L(jω) from the disc center should always be greater than the disc radius. These relations are used to define forbidden regions on the magnitude plot based on the open-loop phase, and forbidden regions on the phase plot based on open-loop magnitude. In the CBode plots shown in Figs. 8 and 9, there are intersections between the open-loop frequency response and the forbidden regions depicted in gray, and therefore the Circle Criterion is not satisfied. In Fig. there are no intersections, so the Circle Criterion is satisfied, and the closed-loop system is stable. The objective of loop-shaping controller design using the CBode plot is to apply compensators to shape the open-loop response to eliminate intersections between the open-loop response and the forbidden regions. A useful property of the CBode magnitude plots is that the boundary functions depend only on the open-loop phase. Thus, multiplying the open-loop by a constant gain will move open-loop magnitude up or down while leaving the forbidden regions on the magnitude plot unchanged []. IV. CONTROLLER DESIGN In this section we illustrate the use of the CBode plot for controller design for an output nonlinearity. The design specifications for this system based on speed of response are Stability in the presence of the valve nonlinearity Zero steady-state error for a step disturbance Open-loop db crossover of rad/s Phase margin of 6 TABLE I Motor and Gearbox Parameters Motor Torque Constant k τ = N-m/A Motor Coil Resistance R = Ω Motor Coil Inductance L =.3 H Motor-Gearbox Inertia J =.4 kg-m Motor-Gearbox Damping b =.7 N-m/s Gearbox Gear Ratio N = 5 The transfer function for this system for a voltage input and valve angle output is P(s)= Nkt s(n (Js+ b)(ls + R)+kt ) where the parameters are given in Table I. (3) 969
4 A. Choosing the Optimum Output Variable Offset One of the important aspects of Circle Criterion design is choosing the output variable offset, which affects the slopes of the sector boundaries. Proper choice of this offset can reduce the ratio of slopes of the two bounding lines thus reducing the size of the Circle Criterion disc. The closer the ratio β n = β/α is to one, the closer the system is to being linear, and the closer the disc is to a point. Fig. 5 shows the flow rate as a function of the valve opening angle, θ =9 α, and the sector lines for the optimum offset of 6. Since the nonlinearity is part of the plant rather than the sensor, after selecting the offset, it is important to normalize the sector by multiplying the plant gain by α and multiplying the sector slopes by /α. This ensures that the linear portion of the plant accounts for the minimum gain. The slope of the lower boundary will be α n =and the slope of upper boundary is β n = β/α. The point +j is the critical point for stability, such that when the compensated openloop system has positive gain margin and has no intersection with the forbidden regions, the system satisfies the Circle Criterion..5 θ = θ = 6.5 θ = θ =. θ = 6. θ = Fig. 7. Circle criterion discs for output offsets θ =, θ =6, and θ = Flow Rate, Q [m 3 /s] Ψ( ) α 6 β CBode plot for offset θ = Fig. 5. The nonlinear function Ψ( ) belongs to the sector [α θ,β θ ] depending on the choice of operating point (6 shown). Fig. 6 shows the value of /β n = α/β as a function of the offset angle. The optimal offset θ maximizes this value. It has the effect of moving the right endpoint of the disc on the real axis at /β n = α/β away from the imaginary axis, as illustrated in Fig CBode plot for offset θ = α / β Fig. 6. α/β as a function of offset angle θ (c) CBode plot for offset θ =9 Fig. 8. CBode plots for different output variable offsets. θ =. θ =6. (c) θ =9. 97
5 Figs. 8-8(c) show the CBode plots for valve angle offsets of, 6, and 9 respectively with the gain normalized to at frequency rad/s. The gain of the plant has been multiplied by the corresponding α and the boundary slopes have been multiplied by /α. Fig. 7 shows the discs of the corresponding offsets. The offset of θ =6 gives the smallest value of β n = β/α. B. Loop-Shaping Controller Design The optimal output variable offset is 6, and so we begin the design by considering the CBode plot of Fig. 8. The idea is to design a compensator that provides integral action for zero steady-state error while lifting the phase plot over the phase forbidden region. The compensator must also provide sufficient phase to achieve a phase margin of 6 at ω db = rad/s. The recently invented complex proportional-integral-lead compensator (CPIL) is well suited for this task [8]. Eq. 4 shows the structure of the CPIL transfer function, C cpil (s)= s +ζω z s + ω z s(s + p) (4) where ω z = ω m ( ζtan(φ m )+ ) ζ tan (φ m )+, (5) p = ω m +sin (φ m ) sin (φ m ), (6) The phase contribution at ω m is 3φm 45. The phase angle must satisfy <φ m < 9. The zeros of this compensator are complex when ζ<. Lower damping ratios provide higher gains at low frequencies and a steeper phase peak. The highest point of the phase forbidden region in Fig. 8 is about 8 above the open-loop phase at frequency 4 rad/s. Therefore, as a first attempt, we choose a CPIL compensator with ζ =.7 providing φ m =8 at frequency ω m =4rad/s. The resulting CBode plot is shown in Fig. 9. The intersections with the forbidden regions are still present but have been reduced compared with the uncompensated system, Fig. 8. The intersection between the phase forbidden region and the phase response occur at the upper frequencies of the forbidden region. The amount of phase added seems to be enough, but the phase should be added at a slightly higher frequency. As a second attempt we keep the parameters of the CPIL the same, but change ω m to 64 rad/s. Fig. shows that with this compensator, all intersections with the forbidden regions are eliminated : : 44.4 : : 4-3 Fig.. CBode plot of open-loop after compensation with CPIL compensator with ζ =.7 providing φ m =8 at frequency ω m =64rad/s. The compensator transfer function in this case is given by C CPIL (s)= s +7.43s +8. s(s + ) (7) Fig. shows the corresponding Nyquist plot and the Circle Criterion disc, confirming that the Circle Criterion is, in fact, satisfied. Checking the phase at ω =rad/s confirms that the phase margin requirement is also satisfied; therefore, the closed-loop nonlinear feedback system is stable Fig. 9. CBode plot of open-loop after compensation with CPIL compensator with ζ =.7 providing φ m =8 at frequency ω m =4rad/s. Fig.. Nyquist plot of the compensated system along with the Circle Criterion. The system response does not intersect the circle with endpoints - and α/β, guaranteeing the system stability and verifying the CBode approach. is the overall Nyquist plot and is a zoomed in view near the origin. To show the benefit of using smaller damping ratios, Fig. shows the CBode plot for a CPIL compenator with ζ =. providing φ m =8 at frequency ω m =7rad/s, which also satisfies the Circle Criterion and the phase margin requirement; however, the gain at ω =rad/s is about 5.5 db higher than for the ζ =.7 compensator. 97
6 : : 39.5 : : 4-3 Fig.. CBode plot of open-loop after compensation with CPIL compensator with ζ =. providing φ m =8 at frequency ω m =7. V. CONCLUSION This paper derived the nonlinear relationship between flow rate, pressure difference, and valve angle in a butterfly valve from first principles and data in the literature. This derivation is easier to understand and more compact than previous approaches. This paper also applied the Circle Criterion Bode plot to the problem of flow rate control using voltage input to to motor-gear-valve feedback system. The controller design method employed iterative loop-shaping, and a complex proportional-integral-lead compensator was chosen with a damping ratio less than one, which provides larger gain at low frequencies compared to proportional-integral-lead compensators with real zeros. The Nyquist plot confirmed that the Circle Criterion was satisfied, and a simulation indicated an improvement in dynamic response. In order to verify the stability and performance of the compensated system, a simulation of the dynamic system response to a step input was created using MATLAB Simulink. The time-domain response of the uncompensated and CPIL compensated unity-gain feedback systems are presented in Fig. 3. Note that both systems are stable, which is consistent with fact that the Circle Criterion is a conservative indicator of stability (i.e. sufficient but not necessary). However, the compensated system does exhibit a marked improvement in dynamic response. Flow Rate [m 3 /s].5 x Reference Flow Rate REFERENCES [] W. Hassenpelug, Free-streamlines, Computers and Mathematics with Applications, vol. 36, no., pp. 69 9, 998. [] L. Xia and W. Messner, Loop shaping with the circle criterion-bode plot with application to active tape steering, in American Control Conference, 7. ACC 7, 7, pp [3] T. Kimura, T. Tanaka, K. Fujimoto, and K. Ogawa, Hydrodynamic characteristics of a butterfly valve: Prediction of pressure loss characteristics, ISA Transactions, vol. 34, no. 4, pp , 995. [4] T. Sarpkaya, Oblique impact of a bounded stream on a plane lamina, Journal of the Franklin Institute, vol. 67, pp. 9 4, 959. [5], Torque and cavitation characteristics of butterfly valves, ASME Journal of Applied Mechanics, vol. 8, pp. 5 58, 96. [6] J. Park and M. Chung, Study on hydrodynamic torque of a butterfly valve, Journal of Fluids Engineering, vol. 8, p. 9, 6. [7] H. Khalil, Nonlinear systems. Prentice-Hall,. [8] W. Messner, Classical control revisited: variations on a theme, AMC 8. th IEEE International Workshop, pp. 5, March Time [s].5 x 3 Flow Rate [m 3 /s].5 Reference Flow Rate Time [s] Fig. 3. Simulated time-domain responses to a step input of the uncompensated and CPIL compensated unity-gain feedback systems. 97
Lesson 37 Transmission Of Air In Air Conditioning Ducts
Lesson 37 Transmission Of Air In Air Conditioning Ducts Version 1 ME, IIT Kharagpur 1 The specific objectives of this chapter are to: 1. Describe an Air Handling Unit (AHU) and its functions (Section 37.1).
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationHomework 7 - Solutions
Homework 7 - Solutions Note: This homework is worth a total of 48 points. 1. Compensators (9 points) For a unity feedback system given below, with G(s) = K s(s + 5)(s + 11) do the following: (c) Find the
More informationChapter 3 Bernoulli Equation
1 Bernoulli Equation 3.1 Flow Patterns: Streamlines, Pathlines, Streaklines 1) A streamline, is a line that is everywhere tangent to the velocity vector at a given instant. Examples of streamlines around
More informationFACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4)
FACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4) 1 1.0 Objectives The objective of this experiment is to calculate loss coefficient (K
More informationIterative Controller Tuning Using Bode s Integrals
Iterative Controller Tuning Using Bode s Integrals A. Karimi, D. Garcia and R. Longchamp Laboratoire d automatique, École Polytechnique Fédérale de Lausanne (EPFL), 05 Lausanne, Switzerland. email: alireza.karimi@epfl.ch
More informationExperiment (4): Flow measurement
Experiment (4): Flow measurement Introduction: The flow measuring apparatus is used to familiarize the students with typical methods of flow measurement of an incompressible fluid and, at the same time
More informationOutline. Classical Control. Lecture 1
Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction
More informationRichiami di Controlli Automatici
Richiami di Controlli Automatici Gianmaria De Tommasi 1 1 Università degli Studi di Napoli Federico II detommas@unina.it Ottobre 2012 Corsi AnsaldoBreda G. De Tommasi (UNINA) Richiami di Controlli Automatici
More informationIndex. Index. More information. in this web service Cambridge University Press
A-type elements, 4 7, 18, 31, 168, 198, 202, 219, 220, 222, 225 A-type variables. See Across variable ac current, 172, 251 ac induction motor, 251 Acceleration rotational, 30 translational, 16 Accumulator,
More information2 Internal Fluid Flow
Internal Fluid Flow.1 Definitions Fluid Dynamics The study of fluids in motion. Static Pressure The pressure at a given point exerted by the static head of the fluid present directly above that point.
More informationClassify a transfer function to see which order or ramp it can follow and with which expected error.
Dr. J. Tani, Prof. Dr. E. Frazzoli 5-059-00 Control Systems I (Autumn 208) Exercise Set 0 Topic: Specifications for Feedback Systems Discussion: 30.. 208 Learning objectives: The student can grizzi@ethz.ch,
More informationMass of fluid leaving per unit time
5 ENERGY EQUATION OF FLUID MOTION 5.1 Eulerian Approach & Control Volume In order to develop the equations that describe a flow, it is assumed that fluids are subject to certain fundamental laws of physics.
More informationFEEDBACK CONTROL SYSTEMS
FEEDBAC CONTROL SYSTEMS. Control System Design. Open and Closed-Loop Control Systems 3. Why Closed-Loop Control? 4. Case Study --- Speed Control of a DC Motor 5. Steady-State Errors in Unity Feedback Control
More informationChapter 7 The Energy Equation
Chapter 7 The Energy Equation 7.1 Energy, Work, and Power When matter has energy, the matter can be used to do work. A fluid can have several forms of energy. For example a fluid jet has kinetic energy,
More informationFrequency Response Techniques
4th Edition T E N Frequency Response Techniques SOLUTION TO CASE STUDY CHALLENGE Antenna Control: Stability Design and Transient Performance First find the forward transfer function, G(s). Pot: K 1 = 10
More informationQ1 Give answers to all of the following questions (5 marks each):
FLUID MECHANICS First Year Exam Solutions 03 Q Give answers to all of the following questions (5 marks each): (a) A cylinder of m in diameter is made with material of relative density 0.5. It is moored
More informationFLUID MECHANICS PROF. DR. METİN GÜNER COMPILER
FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 5. FLOW IN PIPES 5.1.3. Pressure and Shear Stress
More informationManufacturing Equipment Control
QUESTION 1 An electric drive spindle has the following parameters: J m = 2 1 3 kg m 2, R a = 8 Ω, K t =.5 N m/a, K v =.5 V/(rad/s), K a = 2, J s = 4 1 2 kg m 2, and K s =.3. Ignore electrical dynamics
More informationH(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at )
.7 Quiz Solutions Problem : a H(s) = s a a) Calculate the zero order hold equivalent H eq (z). H eq (z) = z z G(s) Z{ } s G(s) a Z{ } = Z{ s s(s a ) } G(s) A B Z{ } = Z{ + } s s(s + a) s(s a) G(s) a a
More informationME332 FLUID MECHANICS LABORATORY (PART II)
ME332 FLUID MECHANICS LABORATORY (PART II) Mihir Sen Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, IN 46556 Version: April 2, 2002 Contents Unit 5: Momentum transfer
More informationThe basic principle to be used in mechanical systems to derive a mathematical model is Newton s law,
Chapter. DYNAMIC MODELING Understanding the nature of the process to be controlled is a central issue for a control engineer. Thus the engineer must construct a model of the process with whatever information
More informationFrequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability
Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods
More informationAppendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2)
Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2) For all calculations in this book, you can use the MathCad software or any other mathematical software that you are familiar
More informationME 375 Final Examination Thursday, May 7, 2015 SOLUTION
ME 375 Final Examination Thursday, May 7, 2015 SOLUTION POBLEM 1 (25%) negligible mass wheels negligible mass wheels v motor no slip ω r r F D O no slip e in Motor% Cart%with%motor%a,ached% The coupled
More informationMAE 143B - Homework 9
MAE 43B - Homework 9 7.2 2 2 3.8.6.4.2.2 9 8 2 2 3 a) G(s) = (s+)(s+).4.6.8.2.2.4.6.8. Polar plot; red for negative ; no encirclements of, a.s. under unit feedback... 2 2 3. 4 9 2 2 3 h) G(s) = s+ s(s+)..2.4.6.8.2.4
More informationCompensator Design to Improve Transient Performance Using Root Locus
1 Compensator Design to Improve Transient Performance Using Root Locus Prof. Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning
More informationRobust fixed-order H Controller Design for Spectral Models by Convex Optimization
Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Alireza Karimi, Gorka Galdos and Roland Longchamp Abstract A new approach for robust fixed-order H controller design by
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationLecture 6 Classical Control Overview IV. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 6 Classical Control Overview IV Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lead Lag Compensator Design Dr. Radhakant Padhi Asst.
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationReynolds, an engineering professor in early 1880 demonstrated two different types of flow through an experiment:
7 STEADY FLOW IN PIPES 7.1 Reynolds Number Reynolds, an engineering professor in early 1880 demonstrated two different types of flow through an experiment: Laminar flow Turbulent flow Reynolds apparatus
More informationequation 4.1 INTRODUCTION
4 The momentum equation 4.1 INTRODUCTION It is often important to determine the force produced on a solid body by fluid flowing steadily over or through it. For example, there is the force exerted on a
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
More informationImplementation Issues for the Virtual Spring
Implementation Issues for the Virtual Spring J. S. Freudenberg EECS 461 Embedded Control Systems 1 Introduction One of the tasks in Lab 4 is to attach the haptic wheel to a virtual reference position with
More informationEXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER
EXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER 1.1 AIM: To determine the co-efficient of discharge of the orifice meter 1.2 EQUIPMENTS REQUIRED: Orifice meter test rig, Stopwatch 1.3 PREPARATION 1.3.1
More informationME332 FLUID MECHANICS LABORATORY (PART I)
ME332 FLUID MECHANICS LABORATORY (PART I) Mihir Sen Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, IN 46556 Version: January 14, 2002 Contents Unit 1: Hydrostatics
More information(Continued on next page)
(Continued on next page) 18.2 Roots of Stability Nyquist Criterion 87 e(s) 1 S(s) = =, r(s) 1 + P (s)c(s) where P (s) represents the plant transfer function, and C(s) the compensator. The closedloop characteristic
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationChapter 10 Flow in Conduits
Chapter 10 Flow in Conduits 10.1 Classifying Flow Laminar Flow and Turbulent Flow Laminar flow Unpredictable Turbulent flow Near entrance: undeveloped developing flow In developing flow, the wall shear
More informationStabilizing the dual inverted pendulum
Stabilizing the dual inverted pendulum Taylor W. Barton Massachusetts Institute of Technology, Cambridge, MA 02139 USA (e-mail: tbarton@mit.edu) Abstract: A classical control approach to stabilizing a
More informationLOSSES DUE TO PIPE FITTINGS
LOSSES DUE TO PIPE FITTINGS Aim: To determine the losses across the fittings in a pipe network Theory: The resistance to flow in a pipe network causes loss in the pressure head along the flow. The overall
More informationR10 JNTUWORLD B 1 M 1 K 2 M 2. f(t) Figure 1
Code No: R06 R0 SET - II B. Tech II Semester Regular Examinations April/May 03 CONTROL SYSTEMS (Com. to EEE, ECE, EIE, ECC, AE) Time: 3 hours Max. Marks: 75 Answer any FIVE Questions All Questions carry
More information1-Reynold s Experiment
Lect.No.8 2 nd Semester Flow Dynamics in Closed Conduit (Pipe Flow) 1 of 21 The flow in closed conduit ( flow in pipe ) is differ from this occur in open channel where the flow in pipe is at a pressure
More information10.52 Mechanics of Fluids Spring 2006 Problem Set 3
10.52 Mechanics of Fluids Spring 2006 Problem Set 3 Problem 1 Mass transfer studies involving the transport of a solute from a gas to a liquid often involve the use of a laminar jet of liquid. The situation
More informationThe Bernoulli Equation
The Bernoulli Equation The most used and the most abused equation in fluid mechanics. Newton s Second Law: F = ma In general, most real flows are 3-D, unsteady (x, y, z, t; r,θ, z, t; etc) Let consider
More informationLinear Control Systems Lecture #3 - Frequency Domain Analysis. Guillaume Drion Academic year
Linear Control Systems Lecture #3 - Frequency Domain Analysis Guillaume Drion Academic year 2018-2019 1 Goal and Outline Goal: To be able to analyze the stability and robustness of a closed-loop system
More information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
More informationELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 2010/2011 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques
CAIRO UNIVERSITY FACULTY OF ENGINEERING ELECTRONICS & COMMUNICATIONS DEP. 3rd YEAR, 00/0 CONTROL ENGINEERING SHEET 5 Lead-Lag Compensation Techniques [] For the following system, Design a compensator such
More information5 ENERGY EQUATION OF FLUID MOTION
5 ENERGY EQUATION OF FLUID MOTION 5.1 Introduction In order to develop the equations that describe a flow, it is assumed that fluids are subject to certain fundamental laws of physics. The pertinent laws
More informationAnswers to questions in each section should be tied together and handed in separately.
EGT0 ENGINEERING TRIPOS PART IA Wednesday 4 June 014 9 to 1 Paper 1 MECHANICAL ENGINEERING Answer all questions. The approximate number of marks allocated to each part of a question is indicated in the
More informationDigital Control Systems
Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist
More informationThe loop shaping paradigm. Lecture 7. Loop analysis of feedback systems (2) Essential specifications (2)
Lecture 7. Loop analysis of feedback systems (2). Loop shaping 2. Performance limitations The loop shaping paradigm. Estimate performance and robustness of the feedback system from the loop transfer L(jω)
More informationControl Systems I Lecture 10: System Specifications
Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture
More informationChapter 8: Flow in Pipes
8-1 Introduction 8-2 Laminar and Turbulent Flows 8-3 The Entrance Region 8-4 Laminar Flow in Pipes 8-5 Turbulent Flow in Pipes 8-6 Fully Developed Pipe Flow 8-7 Minor Losses 8-8 Piping Networks and Pump
More informationFrequency domain analysis
Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationChapter Four fluid flow mass, energy, Bernoulli and momentum
4-1Conservation of Mass Principle Consider a control volume of arbitrary shape, as shown in Fig (4-1). Figure (4-1): the differential control volume and differential control volume (Total mass entering
More informationExam. 135 minutes + 15 minutes reading time
Exam January 23, 27 Control Systems I (5-59-L) Prof. Emilio Frazzoli Exam Exam Duration: 35 minutes + 5 minutes reading time Number of Problems: 45 Number of Points: 53 Permitted aids: Important: 4 pages
More informationChapter 8: Flow in Pipes
Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate the major and minor losses associated with pipe flow in piping networks
More informationFriction Factors and Drag Coefficients
Levicky 1 Friction Factors and Drag Coefficients Several equations that we have seen have included terms to represent dissipation of energy due to the viscous nature of fluid flow. For example, in the
More informationControl Systems I. Lecture 9: The Nyquist condition
Control Systems I Lecture 9: The Nyquist condition adings: Guzzella, Chapter 9.4 6 Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Emilio Frazzoli Institute
More informationShell/Integral Balances (SIB)
Shell/Integral Balances (SIB) Shell/Integral Balances Shell or integral (macroscopic) balances are often relatively simple to solve, both conceptually and mechanically, as only limited data is necessary.
More informationr + - FINAL June 12, 2012 MAE 143B Linear Control Prof. M. Krstic
MAE 43B Linear Control Prof. M. Krstic FINAL June, One sheet of hand-written notes (two pages). Present your reasoning and calculations clearly. Inconsistent etchings will not be graded. Write answers
More informationFLUID MECHANICS PROF. DR. METİN GÜNER COMPILER
FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 4. ELEMENTARY FLUID DYNAMICS -THE BERNOULLI EQUATION
More informationMajor and Minor Losses
Abstract Major and Minor Losses Caitlyn Collazo, Team 2 (1:00 pm) A Technovate fluid circuit system was used to determine the pressure drop across a pipe section and across an orifice. These pressure drops
More informationThe Design of Gating System 3. Theoretical considerations in gating design
MME 345 Lecture 16 The Design of Gating System 3. Theoretical considerations in gating design Ref: [1] ASM Metal Handbook, Vol. 15: Casting, ASM International [] Taylor, Flemings, and Wulff. Foundry engineering,
More informationTransient response via gain adjustment. Consider a unity feedback system, where G(s) = 2. The closed loop transfer function is. s 2 + 2ζωs + ω 2 n
Design via frequency response Transient response via gain adjustment Consider a unity feedback system, where G(s) = ωn 2. The closed loop transfer function is s(s+2ζω n ) T(s) = ω 2 n s 2 + 2ζωs + ω 2
More informationPositioning Servo Design Example
Positioning Servo Design Example 1 Goal. The goal in this design example is to design a control system that will be used in a pick-and-place robot to move the link of a robot between two positions. Usually
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using
More informationCHAPTER 3 BASIC EQUATIONS IN FLUID MECHANICS NOOR ALIZA AHMAD
CHAPTER 3 BASIC EQUATIONS IN FLUID MECHANICS 1 INTRODUCTION Flow often referred as an ideal fluid. We presume that such a fluid has no viscosity. However, this is an idealized situation that does not exist.
More informationExperiment No.4: Flow through Venturi meter. Background and Theory
Experiment No.4: Flow through Venturi meter Background and Theory Introduction Flow meters are used in the industry to measure the volumetric flow rate of fluids. Differential pressure type flow meters
More informationChapter 2. Classical Control System Design. Dutch Institute of Systems and Control
Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral
More informationAutomatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: Student ID number... Signature...
Automatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: 29..23 Given and family names......................solutions...................... Student ID number..........................
More informationECSE 4962 Control Systems Design. A Brief Tutorial on Control Design
ECSE 4962 Control Systems Design A Brief Tutorial on Control Design Instructor: Professor John T. Wen TA: Ben Potsaid http://www.cat.rpi.edu/~wen/ecse4962s04/ Don t Wait Until The Last Minute! You got
More informationRobust Performance Example #1
Robust Performance Example # The transfer function for a nominal system (plant) is given, along with the transfer function for one extreme system. These two transfer functions define a family of plants
More informationCollocated versus non-collocated control [H04Q7]
Collocated versus non-collocated control [H04Q7] Jan Swevers September 2008 0-0 Contents Some concepts of structural dynamics Collocated versus non-collocated control Summary This lecture is based on parts
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationFE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)
Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
More informationApplied Fluid Mechanics
Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More information100 (s + 10) (s + 100) e 0.5s. s 100 (s + 10) (s + 100). G(s) =
1 AME 3315; Spring 215; Midterm 2 Review (not graded) Problems: 9.3 9.8 9.9 9.12 except parts 5 and 6. 9.13 except parts 4 and 5 9.28 9.34 You are given the transfer function: G(s) = 1) Plot the bode plot
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationObjectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation
Objectives Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Conservation of Mass Conservation of Mass Mass, like energy, is a conserved
More informationViscous Flow in Ducts
Dr. M. Siavashi Iran University of Science and Technology Spring 2014 Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 5. 2. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More informationControl of Electromechanical Systems
Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance
More informationC(s) R(s) 1 C(s) C(s) C(s) = s - T. Ts + 1 = 1 s - 1. s + (1 T) Taking the inverse Laplace transform of Equation (5 2), we obtain
analyses of the step response, ramp response, and impulse response of the second-order systems are presented. Section 5 4 discusses the transient-response analysis of higherorder systems. Section 5 5 gives
More informationMAS107 Control Theory Exam Solutions 2008
MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve
More informationSYSTEMS VS. CONTROL VOLUMES. Control volume CV (open system): Arbitrary geometric space, surrounded by control surfaces (CS)
SYSTEMS VS. CONTROL VOLUMES System (closed system): Predefined mass m, surrounded by a system boundary Control volume CV (open system): Arbitrary geometric space, surrounded by control surfaces (CS) Many
More informationLab 3: Model based Position Control of a Cart
I. Objective Lab 3: Model based Position Control of a Cart The goal of this lab is to help understand the methodology to design a controller using the given plant dynamics. Specifically, we would do position
More informationPart A: 1 pts each, 10 pts total, no partial credit.
Part A: 1 pts each, 10 pts total, no partial credit. 1) (Correct: 1 pt/ Wrong: -3 pts). The sum of static, dynamic, and hydrostatic pressures is constant when flow is steady, irrotational, incompressible,
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using the
More informationApplied Fluid Mechanics
Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and
More informationFluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh
Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create
More informationDimensions represent classes of units we use to describe a physical quantity. Most fluid problems involve four primary dimensions
BEE 5330 Fluids FE Review, Feb 24, 2010 1 A fluid is a substance that can not support a shear stress. Liquids differ from gasses in that liquids that do not completely fill a container will form a free
More informationDESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD
206 Spring Semester ELEC733 Digital Control System LECTURE 7: DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD For a unit ramp input Tz Ez ( ) 2 ( z ) D( z) G( z) Tz e( ) lim( z) z 2 ( z ) D( z)
More informationFE Exam Fluids Review October 23, Important Concepts
FE Exam Fluids Review October 3, 013 mportant Concepts Density, specific volume, specific weight, specific gravity (Water 1000 kg/m^3, Air 1. kg/m^3) Meaning & Symbols? Stress, Pressure, Viscosity; Meaning
More informationTopic # Feedback Control Systems
Topic #19 16.31 Feedback Control Systems Stengel Chapter 6 Question: how well do the large gain and phase margins discussed for LQR map over to DOFB using LQR and LQE (called LQG)? Fall 2010 16.30/31 19
More information