logo Use SA10RAM to get 10%* Discount.
Order Now logo

Ask This Question To Be Solved By Our ExpertsGet A+ Grade Solution Guaranteed

expert
Tracy BartramEnglish
(5/5)

986 Answers

Hire Me
expert
Paul BurlinggComputer science
(5/5)

791 Answers

Hire Me
expert
Avinash KumarData mining
(/5)

550 Answers

Hire Me
expert
Dillon RamosResume writing
(5/5)

511 Answers

Hire Me
Others
(5/5)

Write your own Python code to implement both numerical methods to simulate the system.

INSTRUCTIONS TO CANDIDATES
ANSWER ALL QUESTIONS

PART A (Modelling)

 

Question 1 (2 marks)

Provide a real-life example where a general mass-spring-damper system can be found.

 

Question 2 (10 marks)

Using a free-body-diagram, show that the differential equations representing the system in Figure 1 are given by:

𝑑2π‘₯1(𝑑) 𝑑π‘₯1(𝑑) 𝑑π‘₯2(𝑑)

2 + 2 + π‘₯2(𝑑) − 2 = 0

𝑑𝑑2 𝑑𝑑 1 𝑑𝑑 Eq.(2)

𝑑2π‘₯2(𝑑) 𝑑π‘₯2(𝑑) 𝑑π‘₯1(𝑑)

2 + 2 − 2 = 𝑓(𝑑)

𝑑𝑑2 𝑑𝑑 𝑑𝑑 Eq.(3)

 

Question 3 (10 marks)

Linearise the system for the point π‘₯1(0) = 1, π‘₯2(0) = 0, 𝑓(0) = 0, and show that the linearised differential equations are given by:

𝑑2𝛿π‘₯1(𝑑) 𝑑𝛿π‘₯1(𝑑) 𝑑𝛿π‘₯2(𝑑)

2 + 2 + (1 + 2𝛿π‘₯1(𝑑)) − 2 = 0

𝑑𝑑2 𝑑𝑑 𝑑𝑑 Eq.(4)

𝑑2𝛿π‘₯2(𝑑) 𝑑𝛿π‘₯2(𝑑) 𝑑𝛿π‘₯1(𝑑)

2 + 2 − 2 = 𝛿𝑓(𝑑)

𝑑𝑑2 𝑑𝑑 𝑑𝑑 Eq.(5)

 

Hints: 𝛿π‘₯1(𝑑) = π‘₯1(𝑑) − 1, where 1 is the linearisation point.

 

PART B (Analytical and numerical methods)

 

Remark: The remaining questions will be based on the linearised system given in Equation Eq.(6) and Equation Eq.(7).

 

For the remaining questions, replace the linear variable 𝛿π‘₯1(𝑑) with π‘₯1(𝑑), 𝛿π‘₯2(𝑑) in Eq.(4) and Eq.(5) with π‘₯2(𝑑) and 𝛿𝑓(𝑑) with 𝑓(𝑑) (i.e. ignoring the linearisation point). In other words, the remaining questions are based on the result of the linearisation process, which is given by:

𝑑2π‘₯1(𝑑) 𝑑π‘₯1(𝑑) 𝑑π‘₯2(𝑑)

2 + 2 + (1 + 2π‘₯1(𝑑)) − 2 = 0

𝑑𝑑2 𝑑𝑑 𝑑𝑑 Eq.(6)

𝑑2π‘₯2(𝑑) 𝑑π‘₯2(𝑑) 𝑑π‘₯1(𝑑)

2 + 2 − 2 = 𝑓(𝑑)

𝑑𝑑2 𝑑𝑑 𝑑𝑑 Eq.(7)

 

 

Question 4 (6 marks)

Apply Laplace transform on the linearised system in equations Eq.(6)-Eq.(7). Assume zero initial conditions, i.e. π‘₯1(0) = 0, π‘₯Μ‡1(0) = 0, π‘₯2(0) = 0, π‘₯Μ‡2(0) = 0.

 

Question 5 (7 marks)

Using the results from Question (4),

a) Obtain the transfer function for the linearised system (i.e. equations Eq.(6) and

Eq.(7)). Note that the input of the system is 𝑓(𝑑) and the output is π‘₯2(𝑑).

(5 marks)

b) What is the order of the system? Justify your answer.

(2 marks)

 

Question 6 (20 marks)

Using the transfer function from Question (5)

a) Obtain and clearly list all the poles and zeros of the linearised system. (Remark: there should be two zeros and four poles. Two of the poles are at s = −0.12256 ± 0.74486i and one at s = 0. You need to find the other one real pole.)

(Remark 2: the location of the remaining pole is between -2 and 0).

 

 

b) Determine if the system is stable. Justify your answer.

 

c) Draw the “𝑠-plane”.

 

(12 marks)

 

(3 marks)

 

(5 marks)

 

 

Question 7 (5 marks)

Use Final Value Theorem to predict the value of π‘₯2(𝑑) (if available) when the input is a unit impulse. Justify your answer.

 

Question 8 (40 marks)

Using Equation (6) and Equation (7), create a simulation of the system using both:

a) Euler method, and

b) Heun (Runge-Kutta second order) method Assume that the input force 𝑓(𝑑) is a unit impulse.

For each method (20 marks × πŸ = πŸ’πŸŽ marks):

i. Write your own Python code that creates the simulation of the system

(10 marks)

 

ii. Select an appropriate step size β„Ž = Δ𝑑. Explain your reasoning.

(2 marks)

 

iii. Select an appropriate end time for the simulation. Explain your reasoning.

(2 marks)

 

iv. Calculate the error (𝐸𝑑) and percentage error (|πœ–π‘‘|%) and discuss the results.

Remark: Use the final value from Question (7) as the exact numerical solution

(3 marks)

 

v. Plot the output responses (i.e. plot the output position π‘₯2(𝑑) against time in seconds) and discuss the results in terms of general system behaviour (e.g. underdamped or overdamped etc).

(3 marks)

 

Important remark for Question 8:

- Write your own Python code to implement both numerical methods to simulate the system.

- Use a single Python file to implement both methods.

- You can use an existing/example program to check your solution.

- However, do not submit the existing/example program as your submission.

- To avoid compatibility issues when marking the Python code, please use Spyder

4.2.5 (Python 3.8). (You can download Spyder through Anaconda navigator – see the webpage: https://www.anaconda.com/products/individual).

 

(5/5)
Attachments:

Related Questions

. The fundamental operations of create, read, update, and delete (CRUD) in either Python or Java

CS 340 Milestone One Guidelines and Rubric  Overview: For this assignment, you will implement the fundamental operations of create, read, update,

. Develop a program to emulate a purchase transaction at a retail store. ThisΒ  program will have two classes, a LineItem class and a Transaction class

Retail Transaction Programming Project  Project Requirements:  Develop a program to emulate a purchase transaction at a retail store. This

. The following program contains five errors. Identify the errors and fix them

7COM1028   Secure Systems Programming   Referral Coursework: Secure

. Accepts the following from a user: Item Name Item Quantity Item Price Allows the user to create a file to store the sales receipt contents

Create a GUI program that:Accepts the following from a user:Item NameItem QuantityItem PriceAllows the user to create a file to store the sales receip

. The final project will encompass developing a web service using a software stack and implementing an industry-standard interface. Regardless of whether you choose to pursue application development goals as a pure developer or as a software engineer

CS 340 Final Project Guidelines and Rubric  Overview The final project will encompass developing a web service using a software stack and impleme

Get Free Quote!

343 Experts Online