Find the closed-loop transfer function

Assignment Help MATLAB Programming
Reference no: EM131220818

Assignment 1 -

Question 1: Design a pumping system that will transfer raw water from Reservoir A to Reservoirs B and C. A control valve shall be located in Line D-C next to Reservoir C to control the share of flow between Reservoir B and C. See Figure 1.

The pump station shall contain a main suction isolation valve and at least three pumps. It shall be located some 200 m downstream of Reservoir A. The pumps shall be arranged and valves be provided and duty/standby pumps be selected as per "good design practice".

System operating modes:

MODE (I) (Design condition)

- Normal flow rate from reservoir A to be 250 litres/second

- Discharge into Reservoir B to be 50 litres/s and into Reservoir C to be 200 litres/s.

MODE (II) (All flow to Reservoir C)

- The same number of pumps running as in MODE (I)

- Control valve in Line D-C set to cause all flow running into Reservoir C, i.e. no water to be discharged into or drained from Reservoir B.

Tasks:

1. Draw a detailed system schematic for MODE (I), allowing for sufficient valves (no non-return valves in Line D-B) and at least three 90 degree bends per pipe section.

2. Show the selected pipeline diameters and their associated head losses (based on good engineering practice) in a table and justify the reasons for your component selection.

3. Show the expected pump operating points for MODE (I) and demonstrate these by showing pump performance and system resistance curves.

4. Determine the amount of the main suction isolation valve to be throttled in MODE (I) in order to match an NPSH safety factor of 1.15 for the pumps.

5. Show pump operating points, system resistance curves and control valve settings for operation in MODE (II).

QUESTION 2: From the information given below, calculate the total life cycle cost for each pipeline diameter and determine the optimum pipeline diameter.

Pipeline system parameters:

System flow rate, Q = 1000 1/s

Pipeline static lift, Hs = 500 m

Pipeline length, Lp = 80 km

Pipe friction factor, f = 0.02

System design life 30 years, NPV factor = 7.50

Cost of pumping power = $0.05/kWh

Hours of operation = 24 hours/day, 365 days/year

Pump set efficiency = 85%

Pipeline maintenance cost = $20,000 pa

Cost per pump station = $ 5,000,000

Number of pump stations = 2

Pipe diameter and costs

Pipe internal diameter Pipe supply and install cost

700 mm $600/m

800 mm $650/m

900 mm $700/m

1000 mm $750/m

In the calculation, you may wish to ignore losses in valves, fittings, inlet and outlet.

QUESTION 3: Write a factsheet or a short essay/report of literature review on the current research and development in pumping technology. You may choose to focus on one of the following topics (other topics should be pre-approved by the course examiner):

  • Strategies for energy efficiency of pumping systems
  • Application and importance of pumping systems

You may need to cover:

  • Introduction - importance of the issue
  • Technology currently being investigated (the working principle, pros and cons)
  • Future research
  • Summary

You may base your essay on the materials assigned for this module and also any other information available on the internet and in the library. Concise and clear to the point presentation style are strongly encouraged. The report should be well structured with appropriate uses of headings (as suggested above) and sub-headings. The key points should be clearly emphasized and the extensive uses of bullet points (dot points) or tables are encouraged.

Marks will be allocated according to the following three criteria:

  • General presentation and layout, grammar, spelling, referencing, logic of the presentation.
  • Accuracy/appropriateness of factual materials.
  • Discussion/interpretation of materials.

Length: A standard essay format of 2-4 A4 pages (including 3-5 references) would be appropriate.

267_Figure1.png

Assignment 2 - Time Response and Controller Design

Please provide your solutions to the below questions. Show all work and present your work logically, justifying all assumptions. The questions marked with an asterisk (*) include a MATLAB component.

Exercise 1: Reduction of subsystems

Find the transfer function T (s) = C(S)/C(R) of the block diagram in Figure given that:

G1 = 2/(s+2), G2 = 1/s, G3 = 5

G4 = ((s+5)/(s+3)(s+10)), G5 = 2, G6 = s, G7 = 5/s

2380_Fig1.png

Also, write a MATLAB script using feedback, series and parallel commands and confirm that your handwritten answer and MATLAB agree.

Exercise 2: What's the damping?

Consider the motor/load feedback control system shown in Figure 2 below. The load damping DL is unknown, but when R(s) is a step, the response of θL has an overshoot of 20%. Do the following:

a) Find the transfer function θL(s)/R(s) (you'll have a DL in your equation since you don't know it yet)

b) Calculate the damping ratio that corresponds to 20% overshoot

c) Find DL using the information from a and b).

2290_Fig2.png

Exercise 3: A proportional controller

For the velocity control system in Figure 3:

a) Find the closed-loop transfer function T(s) = c(s)/R(s).

b) Find a value of Kp that will yield less than 10% overshoot for the closed-loop system. (Note: ignore the zero dynamics to calculate Kp initially).

c) Is the second order component of the closed-loop system dominant?

d) What is the magnitude of step input that will give a steady-state velocity of c(∞) = 50 m/s?

e) Using your Kp from part b), write a MATLAB script that calculates the closed-loop transfer function, T(s) = c(s)/R(s).

f) Simulate the step response of T(s). Is the overshoot 10% as you designed? Discuss.

126_Fig3.png

Exercise 4: Linearization of a vehicle moving through air

Consider the simple model of a vehicle moving through air

mv· = F - kv2

where v is the velocity of the vehicle, F is the forward traction force, and k is a drag parameter that depends on the properties of the air and the geometry of the vehicle. The parameters are m = 1000kg and k = 10N·s2/m2. Do the following:

a) Linearize the model about the steady-state F- = 1000N.

b) Find the transfer function v'(s)/F'(s) for the linearized system in a)

c) Consider a vehicle that is at the steady-state you found in a). A step change in the traction force is introduced, F'(s) = 200/s. Use your transfer function model to predict time response of the velocity v(t).

Exercise 5: Stability

Consider a unity feedback (H(s) = 1) configuration Figure 4. For the following systems, determine the range of K that will maintain stability.

(a) Gc(s) = K, Ga(s) = 1/s(s2+1), and Gp(s) = 1/(s2+s+1)

(b) Gc(s) = Ks, Ga(s) = 1/(s+1), and Gp(s) =  1/((s+2)(s+5))

2350_Fig4.png

Exercise 6: Disc drive problem

Recall the disc drive problem introduced in tutorial. Last time we developed the transfer functions seen in Figure 5. Do the following:

a) Find the transfer function Gd(s) = Po(s)/Td(s) when Gc(s) = Kc.

b) For Gc(s) = Kp, find the range of Kp for which the closed-loop system is stable.

c) Write a MATLAB script using feedback, series and parallel commands to obtain Po(s)/Pt(s) and Po(s)/Td(s) when Kc = 12.14.

957_Fig5.png

Exercise 7: Disc drive controller

Recall the disc drive problem from Tutorials, where we demonstrated that the open-loop system can be written as

1160_Fig6.png

a) Consider the controller that we designed in tutorial:

Gc(s) = Kc = 12.14

Find the steady-state error to a ramp input with this controller. If we wish to reduce this error to 0.005, can we do it with a different Kc? (Hint: consider your stability limits!)

b) Well try a more complex controller of the form

Gc(s) = Kc(s + a)

which is sometimes called a proportional-derivative controller. Find the closed-loop transfer function, T(s) = Po(s)/Pt(s).

c) Find the conditions (inequalities) on Kc and a such that the closed-loop system is stable. Use MATLAB to plot the stability boundaries (again, inequalities) on a Kc vs a plot.

d) Find values of Kc and a such that the steady-state error to a ramp is less than 0.005.

 Exercise 8: Using Root Locus

Consider the feedback control system in Figure 7. In this exercise, we'll walk through designing Gc(s) with different levels of complexity.

1115_Fig7.png

To this ends, do the following by hand (unless otherwise stated):

a) Sketch (by hand) the root locus and find the dosed loop poles when Gc(s) = Kc = 1. Also: find the steady-state error to a step and ramp inputs, ζ and the settling time.

b) In order to improve the transient response, a PD controller of the form

Gc(s) = Kc(s + a)

is being considered. Determine the values of K and a so that the closed loop system has an overshoot of 1.6 and a settling time of 2s for a step input. (Hint: use a to ensure that the poles are on the root locus). You may ignore the effect of the zero dynamics for this part.

c) What is the steady-state error to a ramp?

d) We now require that steady-state error to a ramp is eliminated. Your boss has told you to "just add an integrator" to the controller to eliminate the steady-state error. Sketch the root locus and use it to tell your boss why this won't work.

e) A smarter way to eliminate the steady-state error is to use a RID controller of the form

Gc(s) =  (Kc(s + a)(s + b)/s)

where a is the same as in b) and b is close to the origin. Select a suitable value for b so that the desired poles from b) are on the root locus. If necessary (it may or may not be!), adjust the value of Kc to ensure that the desired closed-loop poles are achieved. You may use MATIAB to check the location of the closed-loop poles.

f) Simulate the closed-loop step response in MATLAB using the PID controller from e). Twiddle Kc a and b (if necessary) to achieve your targets (i.e. 16% overshoot and a settling time of 2s).

Exercise 9: Lead-Lag design using Root Locus

Recall the disc drive problem from Tutorials, where we demonstrated that the open-loop system can be written as

1007_Fig8.png

We will now try to design a lead-lag compensator with the requirements that

  • Overshoot ≤ 10
  • Ts ≤ 75MS
  • eramp(∞)≤ 0.001

Do the following (you may use MATLAB at your leisure, but be sure to explain your logic for your design choices):

a) Use MATLAB to draw the root locus when Gc = Kc.

b) Use MATLAB to draw the region where the dominant closed-loop poles must be to satisfy the transient requirements. Comment on your ability to achieve these requirements with a gain-only controller.

c) Design a lead compensator(s) to meet the transient requirements (i.e. overshoot and settling time).

d) Design a lag compensator to achieve the steady-state tracking requirement.

e) Use MATLAB to compute the resulting closed-loop poles and discuss second order dominance.

Some Hints:

  • You may need to place more than one lead compensator for part c)
  • When assessing second order dominance of the closed-loop system, be sure to cancel poles and zeros (i.e. use minreal).

Reference no: EM131220818

Write program that representing a month-day and year

How to write a program to determine the day of the week - write a program that inputs three integers representing a month, day, and year. For example, 5 12 2015 represents May

Write a program in matlab using the galerkin method

The assignment consists in writing a program in MATLAB Using the galerkin method for finite elements that can solve for: Displacements in x and y of a solid bar that is cons

Calculate and display the altitude and velocity of rocket

Write a program that will solve the problem specified above( for the rocket). However, lit user of your program is to insert the starting height, and the time the rocket tra

What is the uniform background lifetime

You should find that the calculated lifetime is much more spatially uniform than the noisy intensity data, with a single well-defined region of spatial contrast. What is the

What is the average growth rates of the economy

What is the average growth rates of the economy along the path (kH, kL, kH kL, kH, ...) and how it compares with the average growth rate along the the balanced growth path?

Value for the side of the octagonal prism

Value for the side of the octagonal prism, and a scalar value for the height of the octagonal prism. Generate the x,y,z arrays that define the 9 faces of the pyramid (octago

How much did he spend each day

A hungry college student goes to the cafeteria and buys lunch. The next day he spends twice as much. The third day he spends $1 less than he did the second day. At the end o

Generate a sample of the exponential distribution

MS455/MS555 Project - Simulation in Finance. Use the inverse transforms method to generate a sample of the exponential distribution with parameter θ, which is given by its de

Reviews

Write a Review

 
Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd