1. Trang chủ
  2. » Kỹ Thuật - Công Nghệ

Ch 03 solution methods for dynamic models

22 159 0

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

THÔNG TIN TÀI LIỆU

Thông tin cơ bản

Định dạng
Số trang 22
Dung lượng 1,66 MB

Các công cụ chuyển đổi và chỉnh sửa cho tài liệu này

Nội dung

of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien §2.Response Types and Stability • Steady-state: the part of the response that remains with time • Transient: the part of

Trang 1

3 Solution Methods for

Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

System Dynamics 3.01 Solution Methods for Dynamic Models

§1.Differential Equations

- Anordinary differential equation(ODE): an equation containing

ordinary, but not partial, derivatives of the dependent variable

- The subject of system dynamics is time-dependent behavior,the independent variable in our ODEs will be time𝑡

System Dynamics 3.02 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

2.Classification of Differential Equations

- Linear differential equation

System Dynamics 3.03 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations

- The order of differential equation: the order of thehighestderivativeof the dependent variable in the equation

3 𝑥+ 7 𝑥 + 2𝑥 = 5: second order differential equation

- A model can consist of more than one equation

3 𝑥1+ 5𝑥1− 7𝑥2= 5 𝑥

0

𝑡

𝑔 𝑡 𝑑𝑡

System Dynamics 3.04 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations

- Example 3.1.1 Separation of Variables for a Linear Equation

Use separation of variables to solve the following problem for

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations4.Trial Solution Method

Consider 𝑥 + 𝑎𝑥 = 𝑏, 𝑡 ≥ 0

⟹ Find the solution in the form 𝑥 𝑡 = 𝐶 + 𝐷𝑒𝑠𝑡

Subtituiting the solution into the equation

𝑥 + 𝑎𝑥 = 𝑠𝐷𝑒𝑠𝑡+ 𝑎 𝐶 + 𝐷𝑒𝑠𝑡 ⟹ 𝑠 + 𝑎 𝐷𝑒𝑠𝑡+ 𝑎𝐶 = 𝑏Find𝑠, 𝐶

𝑠 + 𝑎 = 0

𝑎𝐶 = 𝑏⟹

𝑠 = −𝑎

𝐶 = 𝑏/𝑎Find𝐷 from the initial condition

𝑥 0 = 𝐶 + 𝐷𝑒0= 𝐶 + 𝐷 ⟹ 𝐷 = 𝑥 0 − 𝐶Therefore

𝑥 𝑡 =𝑏

𝑎+ 𝑥 0 −

𝑏

𝑎 𝑒−𝑎𝑡

System Dynamics 3.06 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

Trang 2

§1.Differential Equations

- Example 3.1.2 Two Distinct, Real Roots

Use trial solution method to solve the following problem for

There are two solutions⟹ the trial form is not suitable

Need an additional term with an arbitrary constant⟹ The

appropriate trial-solution form

𝑥 𝑡 = 𝐶 + 𝐷1𝑒𝑠 1 𝑡+𝐷2𝑒𝑠 2 𝑡

System Dynamics 3.07 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

𝑠1= −2

𝑠2= −5

𝐶 = 2The solution

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations

- Example 3.1.3 Two Repeated, Real Roots

Use trial solution method to solve the following problem for

𝑡 ≥ 0: 5 𝑥 + 20 𝑥 + 20𝑥 = 28, 𝑥 0 = 5, 𝑥 0 = 8

Solution

Trial solution form 𝑥 𝑡 = 𝐶 + 𝐷𝑒𝑠𝑡

Substituting this form into the ODE

5𝑠2+ 20𝑠 + 20 𝐷𝑒𝑠𝑡+ 20𝐶 = 28

⟹ 5𝑠2+ 20𝑠 + 20 = 0

20𝐶 = 28

⟹ 𝑠 = −2, 𝑠 = −2, 𝐶 = 1.4

Need an additional term with an arbitrary constant⟹ The

appropriate trial-solution form

𝑥 𝑡 = 𝐶 + (𝐷1+ 𝐷2𝑡)𝑒𝑠 1 𝑡

System Dynamics 3.09 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

⟹ 𝑠1= −2, 𝐶 = 1.4The solution

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations

- Example 3.1.4 Two Imaginary Roots

Use trial solution method to solve the following problem for

𝑡 ≥ 0: 𝑥 + 16𝑥 = 144, 𝑥 0 = 5, 𝑥 0 = 12

Solution

Trial solution form 𝑥 𝑡 = 𝐶 + 𝐷𝑒𝑠𝑡

Substituting this form into the ODE

System Dynamics 3.11 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

𝑥 𝑡 = 𝐶 + 𝐷1+ 𝐷2𝑐𝑜𝑠4𝑡 + 𝑗 𝐷1− 𝐷2𝑠𝑖𝑛4𝑡

or 𝑥 𝑡 = 𝐶 + 𝐵1𝑐𝑜𝑠4𝑡 + 𝐵2𝑠𝑖𝑛4𝑡Evaluating𝑥(𝑡) and 𝑥(𝑡) at 𝑡 = 0

𝑥 0 = 𝐶 + 𝐵1= 5

𝑥 0 = 4𝐵2= 12 ⟹

𝐵1= −4

𝐵2= +3The solution

𝑥 𝑡 = 9 + 3𝑠𝑖𝑛4𝑡 − 4𝑐𝑜𝑠4𝑡

System Dynamics 3.12 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

Trang 3

§1.Differential Equations

- Example 3.1.5 Motion of a Robot-Arm Link

The equation of motion0.23 + 0.5𝑚𝐿2 𝜃

= 𝑇𝑚− 4.9𝑚𝐿𝑠𝑖𝑛𝜃Solve the equation for the case

𝑇𝑚= 0.5𝑁𝑚,𝑚 = 10𝑘𝑔,𝐿 = 0.3𝑚 Assume the system starts fromrest at 𝜃 = 0 and the angle 𝜃remains small

System Dynamics 3.13 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations

- Example 3.1.6 Two Complex Roots

Use trial solution method to solve the following problem for

𝑡 ≥ 0: 𝑥 + 6 𝑥 + 34𝑥 = 68, 𝑥 0 = 5, 𝑥 0 = 7Solution

Subtituting𝑥 𝑡 = 𝐶 + 𝐷𝑒𝑠𝑡into the ODE

𝑠2+ 6𝑠 + 34 𝐷𝑒𝑠𝑡+ 34𝐶 = 68

⟹ 𝑠2+ 6𝑠 + 34 = 0

34𝐶 = 68⟹ 𝑠 = −3 ± 5𝑗, 𝐶 = 2Therefore

𝑥 𝑡 = 𝐶 + 𝐷1𝑒(−3+5𝑗)𝑡+ 𝐷2𝑒(−3−5𝑗)𝑡

⟹ 𝑥 𝑡 = 𝐶 + 𝑒−3𝑡(𝐷1𝑒5𝑗𝑡+ 𝐷2𝑒−5𝑗𝑡)

⟹ 𝑥 𝑡 = 𝐶 + 𝑒−3𝑡(𝐵1𝑐𝑜𝑠5𝑡 + 𝐵2𝑠𝑖𝑛5𝑡)With the initial conditions𝑥 𝑡 = 2 + 𝑒−3𝑡 3𝑐𝑜𝑠5𝑡 +16

5𝑠𝑖𝑛5𝑡

System Dynamics 3.14 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations

5.Sumary of the trial solution method

Ordinary Differential Equation Solution form

System Dynamics 3.15 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§1.Differential Equations6.Assessment of Solution Behavior

- The characteristic equation can be quickly identified from ODE

by replacing 𝑥 with 𝑠, 𝑥 with 𝑠2, and so forth

- Example

3 𝑥 + 30 𝑥 + 222𝑥 = 148The characteristic equation

𝑠2+ 10𝑠 + 74 = 0

⟹ 𝑠 − −5 − 7𝑗 𝑠 − −5 + 7𝑗 = 0The solution in the form

𝑥 𝑡 = 𝑒−5𝑡 𝐶1𝑠𝑖𝑛7𝑡 + 𝐶2𝑐𝑜𝑠7𝑡 +2

3

𝑎 ≠ 0, 𝑎 2 < 4𝑏 : 𝑠 = 𝜎 ± 𝑗𝜔, 𝜎 = −𝑎/2, 𝜔 = 4𝑏 − 𝑎 2 /2 𝑥 𝑡 = 𝑒 𝜎𝑡 𝐶1𝑠𝑖𝑛𝜔𝑡 + 𝐶2𝑐𝑜𝑠𝜔𝑡 +𝑐

System Dynamics 3.16 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

• Steady-state: the part of the response that remains with time

• Transient: the part of the response that disappears with time

• Free: the part of the response that depends on the initial conditions

• Forced: the part of the response due to the forcing function

𝑥 + 𝑎𝑥 = 𝐵, 𝑎 ≠ 0 𝑥 𝑡 =𝑏

𝑎 + 𝐶𝑒 −𝑎𝑡 System Dynamics 3.17 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

1.The time constantThe 1storder model

𝑥 + 𝑎𝑥 = 𝑏

⟹ 𝑥 +1

𝜏𝑥 = 𝑏, 𝜏 ≡ 1/𝑎give the free response

𝑥 𝑡 = 𝑥 0 𝑒−𝑎𝑡

⟹ 𝑥 𝑡 = 𝑥(0)𝑒−𝜏𝑡

𝜏:time constant

• measure of the exponential decay curve

• estimate how long it will take for the transient response todisappear

𝑥 + 𝑎𝑥 = 𝑏, 𝑎 ≠ 0 𝑥 𝑡 =𝑏𝑎+ 𝐶𝑒 −𝑎𝑡 System Dynamics 3.18 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

Trang 4

§2.Response Types and Stability

From figure

• after 𝑡 = 𝜏, 𝑥(𝑡) has decayed to 37% of its initial value

• after𝑡 = 4𝜏, 𝑥(𝑡) has decayed to02%of its initial value

System Dynamics 3.19 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

The time constant is useful also for analyzing the responsewhen the forcing function is a constant

The total response in terms of𝜏 by substituting 𝑎 = 1/𝜏

𝑥 𝑡 = 𝑏𝜏

steady state

+ 𝑥 0 − 𝑏𝜏 𝑒−𝑡𝜏 transient state

= 𝑥𝑠𝑠+ [𝑥 0 − 𝑥𝑠𝑠]𝑒−𝑡𝜏

The response approaches the constant value𝑏𝜏 as 𝑡 → ∞

𝑥𝑠𝑠= lim

𝑡→∞𝑥(𝑡) = 𝑏𝜏

System Dynamics 3.20 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

The forced response (for which𝑥(0) = 0) is plotted in the figure

• at 𝑡 = 𝜏, the response is 63% of the steady-state value

• at𝑡 = 4𝜏, the response is98%of the steady-state value

• at 𝑡 = 5𝜏, the response is 99% of the steady-state value

For most engineering purposes:𝑥(𝑡) reaches steady-state at 𝑡 = 4𝜏

System Dynamics 3.21 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

- Example 3.2.1 Responses for Second-Order, Distinct Roots

Identify the responses of 𝑥 + 7 𝑥 + 10𝑥 = 𝑐 (roots: −2, −5)Solution

System Dynamics 3.22 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

2.The Dominant-Root Approximation

- The time constant concept is not limited to first-order models

It can also be used to estimate the response time of

higher-order models

- Example 3.2.2Responses for Second-Order, Complex Roots

Identify the responses of the equation

System Dynamics 3.23 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

The coefficients𝐵1and𝐵2can be found in terms of arbitraryinitial conditions in the usual way

𝐵1= 𝑥 0 − 𝑐

34, 𝐵2=

34 𝑥 0 + 102𝑥 0 − 3𝑐170

Model’s time constant: τ = 1/3 ⟹ the response is essentially

at steady state for 𝑡 > 4τ = 4/3

System Dynamics 3.24 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

Trang 5

§2.Response Types and Stability

3.Time Constants and Complex Roots

- Example 3.2.3 Responses for Second-Order, Imaginary Roots

Identify the responses of 𝑥 + 16𝑥 = 𝑐

Solution

The characteristic roots𝑠 = ±4𝑗, the solution form

𝑥 𝑡 = 𝑐

16+ 𝐵1𝑐𝑜𝑠4𝑡 + 𝐵2𝑠𝑖𝑛4𝑡

no terms that disappear as𝑡 → ∞ ⟹ no transient response

The free and forced responses

System Dynamics 3.25 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability4.Stability

- Unstable: the free response approaches∞ as 𝑡 → ∞

- Stable: the free response approaches0

- Neutral stability: the borderline between stable and unstable

The free response does not approach both∞ and 0

System Dynamics 3.26 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

- The stability properties of a linear model are determined from

its characteristics roots

- The first-order model 𝑥 + 𝑎𝑥 = 𝑓(𝑡)

The characteristic equation 𝑠 + 𝑎 = 0

The free response 𝑥 𝑡 = 𝑥(0)𝑒−𝑎𝑡

𝑠 = −𝑎 < 0: 𝑡 ⟶ ∞, 𝑥(𝑡) ⟶ 0: the model is stable

𝑠 = −𝑎 > 0: 𝑡 ⟶ ∞, 𝑥(𝑡) ⟶ ∞ : the model is unstable

𝑠 = 𝑎 = 0: 𝑡 ⟶ ∞, 𝑥(𝑡) ↛ 0, ∞: the model is neural stability

System Dynamics 3.27 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

- The second-order models with the same initial conditions

System Dynamics 3.28 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

- Stability Test for Linear Constant-Coefficient Models

• A constant-coefficient linear model isstableif and only if all

of its characteristic roots have negative real parts

• The model isneutrally stableif one or more roots have a

zero real part, and the remaining roots have negative real

parts

• The model isunstableif any root has a positive real part

System Dynamics 3.29 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability5.A Physical Example

- Pendulum motion equation (𝜃 ≈ 0)

𝑚𝐿2 𝜃 + 𝑐 𝜃 + 𝑚𝑔𝐿𝜃 = 0

- Equilibrium point: 𝜃 = 0

𝑠 =−𝑐 ± 𝑐

2− 4𝑚2𝐿3𝑔2𝑚𝐿2

• 2𝑚𝐿 𝐿𝑔 > 𝑐 > 0: damped oscillating

• 𝑐 > 2𝑚𝐿 𝐿𝑔: no oscillating

• 𝑐 = 0: oscillating about𝜃 = 0

System Dynamics 3.30 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

Roots

Trang 6

§2.Response Types and Stability

- Pendulum motion equation (𝜃 ≈ 0)

𝑚𝐿2 𝜃 + 𝑐 𝜃 + 𝑚𝑔𝐿𝜃 = 0

- Equilibrium point: 𝜃 = 𝜋

𝑠 =−𝑐 ± 𝑐

2+ 4𝑚2𝐿3𝑔2𝑚𝐿2

- Pendulum will not oscillate about𝜃 = 𝜋 but will

continue to fall away if disturbed⟹ the system

is unstable

- The model is based on the assumption that𝜃 ≈

0 ⟹ cannot draw any conclusions from the

model regarding the behavior when 𝜃 is not

near𝜋

System Dynamics 3.31 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

System Dynamics 3.32 Solution Methods for Dynamic Models

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

System Dynamics 3.33 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

7.Stability and Equilibrium

- Anequilibrium: a state of no change

The pendulum in the figure is in equilibrium at

𝜃 = 0 and when perfectly balanced at 𝜃 = 𝜋

• the equilibrium at 𝜃 = 0: stable

• the equilibrium at 𝜃 = 𝜋: unstable

⟹ the same physical system can have different

stability characteristics at different equilibria

- Stability is not a property of the system alone, but is a

property of a specific equilibrium of the system

- When we speak of the stability properties of a model, we are

actually speaking of the stability properties of the specific

equilibrium on which the model is based

System Dynamics 3.34 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§2.Response Types and Stability

- The figure shows a ball on a surface that has a valley and a hill

- The bottom of the valley is an equilibrium, and if the ball isdisplaced slightly from this position, it will

• oscillate forever about the bottom if there is no friction:

neutrally stable

• return to the bottom if friction is present:stable

- If displace the ball so much to the left that it lies outside thevalley, it will never return:locally stablebutglobally unstable

- If the system returns to its equilibrium for any initialdisplacement:globally stable

- The equilibrium on the hilltop isglobally unstable

System Dynamics 3.35 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- The Laplace transforms the problem in time-domain to

problem in s-domain, then applying the solution in s-domain,

and finally using inverse transform to converse the solution

back to the time-domain

- The Laplace transformℒ{𝑥 𝑡 } of a function 𝑥(𝑡) is defined as

ℒ 𝑥 𝑡 = lim

𝑇→∞ 0

𝑇

𝑥 𝑡 𝑒−𝑠𝑡𝑑𝑡but is usually expressed more compactly as

ℒ 𝑥 𝑡 =

0

𝑥 𝑡 𝑒−𝑠𝑡𝑑𝑡

System Dynamics 3.36 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Notation for the Laplace and inverse Laplace transforms

𝑋 𝑠 = ℒ 𝑥 𝑡

𝑥 𝑡 = ℒ−1{𝑋(𝑠)}

-Table of Laplace transform pairs

Trang 7

System Dynamics 3.37 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

System Dynamics 3.38 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method1.Transforms of Common Functions

- Example 3.3.1 Transform of a Constant

Suppose 𝑥(𝑡) = 𝑐 , a constant, for 𝑡 ≥ 0 Determine itsLaplace transform

SolutionFrom the transform definition

System Dynamics 3.39 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- The step function

Unit-step function

𝑢𝑠𝑡 = 0 𝑡 < 0

1 𝑡 > 0Step function

𝑥(𝑡) = 𝑀𝑢𝑠(𝑡)

𝑋 𝑠 = ℒ 𝑀𝑢𝑠𝑡 = 𝑀/𝑠

System Dynamics 3.40 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.2 The Exponential Function

Derive the Laplace transform of the exponential function𝑥(𝑡) = 𝑒−𝑎𝑡, where𝑎 is a constant

SolutionFrom the transform definition

System Dynamics 341 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

2.Properties of the Laplace Transform

System Dynamics 3.42 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.3 The Sine and Cosine Functions

Derive the Laplace transforms of𝑒−𝑎𝑡𝑠𝑖𝑛𝜔𝑡 and 𝑒−𝑎𝑡𝑐𝑜𝑠𝜔𝑡,where𝑎 and 𝜔 are constants

SolutionRecall: - the Euler identity

𝑒𝑗𝜃= 𝑐𝑜𝑠𝜃 + 𝑗𝑠𝑖𝑛𝜃, with 𝜃 = 𝜔𝑡

- the relation1

𝑥 − 𝑗𝑦=

𝑥 + 𝑗𝑦(𝑥 − 𝑗𝑦)(𝑥 + 𝑗𝑦)=

Trang 8

System Dynamics 3.43 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

Comparing the above two equation, the Laplace transforms

of the exponentially decaying sine and cosine functions are

ℒ 𝑒−𝑎𝑡𝑐𝑜𝑠𝜔𝑡 = 𝑠 + 𝑎

(𝑠 + 𝑎)2+𝜔2

ℒ 𝑒−𝑎𝑡𝑠𝑖𝑛𝜔𝑡 = 𝜔

(𝑠 + 𝑎)2+𝜔2

System Dynamics 3.44 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

-Shifting along the 𝑠-axisor multiplication by an exponential

- Example 3.3.4 The Function𝑡𝑒−𝑎𝑡

Derive the Laplace transform of the function𝑡𝑒−𝑎𝑡

= 1(𝑠 + 𝑎)2

System Dynamics 3.45 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.5 The Function𝑡𝑐𝑜𝑠𝜔𝑡

Derive the Laplace transform of the function𝑡𝑐𝑜𝑠𝜔𝑡

System Dynamics 3.46 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.6 The Shifted Step Function

If the discontinuity in the unit-step functionoccurs at𝑡 = 𝐷

System Dynamics 3.47 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

b From the time-shifting property

The pulse= A unit-step + A shifted, negative unit-step

⟹ℒ 𝑃 𝑡 =1

𝑠 1 − 𝑒

−𝑠𝐷

System Dynamics 3.48 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method3.The Derivative Property

Applying integration by parts to the definition of the transform

Trang 9

System Dynamics 3.49 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

4.The Initial Value Theorem

Use to find the value of the function 𝑥(𝑡) at 𝑡 = 0+(a time

infinitesimally greater than0) with given the transform 𝑋(𝑠)

System Dynamics 3.50 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method5.The Final Value Theorem

Use to find the limit of the function𝑥(𝑡) as 𝑡 → ∞

𝑓 ∞ = lim

𝑡→∞𝑓(𝑡) = lim

𝑠→0𝑠𝐹(𝑠)Example

𝑋 𝑠 = 7(𝑠 + 4)2+49

⟹ 𝑥 ∞ = lim

𝑠→0𝑠𝑋(𝑠) = lim

𝑠→0

7𝑠(𝑠 + 4)2+49= 0This is confirmed by evaluating the inverse transform

𝑥 𝑡 = ℒ−1𝑋 𝑠 = 𝑒−4𝑡sin7𝑡

⟹ 𝑥 ∞ = 𝑒−(4×0)sin(7 × 0) = 0

The final value theorem does not apply to a periodic function

System Dynamics 3.51 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

6.Solving Equations with the Laplace Transform

Consider the linear first-order equation

𝑥 + 𝑎𝑥 = 𝑓(𝑡)

𝑓(𝑡): the input 𝑎: a constant

Multiply both sides of the equation by𝑒−𝑠𝑡and then integrate

System Dynamics 3.52 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.8 Step Response of a First-Order Equation

Suppose that the input𝑓(𝑡) of the equation 𝑥 + 𝑎𝑥 = 𝑓(𝑡) is astep function of magnitude𝑀 whose transform is 𝐹(𝑠) = 𝑀/𝑠

Obtain the expression for the complete responseSolution

The forced response is obtained from

System Dynamics 3.53 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.9 Ramp Response of a First-Order Equation

Determine the complete response of the following model,

which has a ramp input 𝑥 + 3𝑥 = 5𝑡, 𝑥 0 = 10

System Dynamics 3.54 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

Comparing the numerators

𝐶2+ 𝐶3= 0

𝐶1+ 3𝐶2= 03𝐶1= 5

𝐶1= 5/3

𝐶2= −5/9

𝐶3= 5/9The forced response

by𝑥(𝑡) = 5𝑡/3 − 5/9

Trang 10

System Dynamics 3.55 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.10 Transform Inversion for Complex Factors

Invert the following transform

⟹ 𝑥 𝑡 = 8𝑒−2𝑡𝑐𝑜𝑠7𝑡 −3

7𝑒

−2𝑡𝑠𝑖𝑛7𝑡

System Dynamics 3.56 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

- Example 3.3.11 Step Response of a Second-Order Equation

Obtain the complete response of the following model

𝑥 + 4 𝑥 + 53𝑥 = 15𝑢𝑠𝑡 𝑥 0 = 8, 𝑥 0 = −19Solution

Transforming the equation gives

𝑠2𝑋 𝑠 − 𝑠𝑥 0 − 𝑥 0 + 4 𝑠𝑋 𝑠 − 𝑥 0 + 53𝑋 𝑠 =15

𝑠Solve for 𝑋(𝑠) using the given initial conditions

𝑋 𝑠 =𝑥 0 𝑠 + 𝑥 0 + 4𝑥(0)

𝑠2+ 4𝑠 + 53 +

15𝑠(𝑠2+ 4𝑠 + 53)

= 8𝑠 + 13

𝑠2+ 4𝑠 + 53+

15𝑠(𝑠2+ 4𝑠 + 53)The first term on the right of eq.(1) corresponds to the freeresponse8𝑒−2𝑡𝑐𝑜𝑠7𝑡 − (3/7)𝑒−2𝑡𝑠𝑖𝑛7𝑡 (Ex 3.3.10)

(1)

(2)

System Dynamics 3.57 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

The second term on the right of eq.(1) corresponds to the

forced response It can be expressed as follows

=(𝐶1+ 𝐶2) 𝑠 + 2

2+ (4𝐶1+2𝐶2+ 7𝐶3)𝑠 +53𝐶1𝑠[ 𝑠 + 22+ 72]

Comparing numerators on the left and right sides

System Dynamics 3.58 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

The complete response is the sum of the free and forcedresponses given by equations

(3)

System Dynamics 3.59 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform Method

Alternatively, combine the terms on the right side of eq.(1)

𝑋 𝑠 =𝑥 0 𝑠

2+ 𝑥 0 + 4𝑥 0 𝑠 + 15𝑠(𝑠2+ 4𝑠 + 53) =

8𝑠2+ 13𝑠 + 15𝑠[ 𝑠 + 22+ 72]

𝐶1+ 𝐶2= 84𝐶1+ 2𝐶2+ 7𝐶3= 13

System Dynamics 3.60 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§3.The Laplace Transform MethodStep response of a second-order equation with complex roots

Using a transform of the following form

𝐴𝑠2+ 𝐵𝑠 + 𝐶𝑠[ 𝑠 + 𝑎2+ 𝑏2]= 𝐶1

𝑥 𝑡 = 𝐶1+ 𝐶2𝑒−𝑎𝑡𝑐𝑜𝑠𝑏𝑡 + 𝐶3𝑒−𝑎𝑡𝑠𝑖𝑛𝑏𝑡

Trang 11

System Dynamics 3.61 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§4.Transfer Function

- Response of a linear system= Free response + Forced response

• To focus the analysis on the effects of the input only by

taking the initial conditions to be zero temporarily

• Then, the free response due to any nonzero initial conditions

can be added to the result

- Consider the model 𝑥 + 𝑎𝑥 = 𝑓(𝑡)

Transforming both sides of the equation

𝑠𝑋 𝑠 + 𝑎𝑋 𝑠 = 𝐹 𝑠obtain thetransfer function

𝑇 𝑠 =𝑋(𝑠)𝐹(𝑠)=

1

𝑠 + 𝑎

- A Transfer Function is the ratio of the output of a system to the

input of a system, in the Laplace domain considering its initial

conditions and equilibrium point to be zero

System Dynamics 3.62 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§4.Transfer Function

The TF concept is extremely useful for several reasons

-Transfer Functions and Software

• Matlab accepts a description based on the TF

• Simulink uses TF as the basis graphical system descriptioncalled the block diagram

-ODE Equivalence

The transfer function⟺ The ODE

-The Transfer Function and Characteristic Roots

• The characteristic polynomial is the denominator of the TF

• The roots of characteristics equation give the informationabout the stability of the system

System Dynamics 3.63 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§4.Transfer Function

- Example 3.4.1 Two Inputs and One Output

Obtain the transfer functions𝑋(𝑠)/𝐹(𝑠) and 𝑋(𝑠)/𝐺(𝑠) for the

temporarily setting the other inputs equal to zero

𝑋(𝑠)

𝐹(𝑠)=

65𝑠2+ 30𝑠 + 40,

𝑋(𝑠)𝐺(𝑠)=

205𝑠2+ 30𝑠 + 40

System Dynamics 3.64 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

§4.Transfer Function

- Example 3.4.2 A System of Equations

a.Obtain the transfer functions𝑋(𝑠)/𝑉(𝑠) and 𝑌(𝑠)/𝑉(𝑠) of thefollowing system of equations

𝑥 = −3𝑥 + 2𝑦

𝑦 = −9𝑦 − 4𝑥 + 3𝑣(𝑡)b.Obtain the forced response for𝑥(𝑡) and 𝑦(𝑡) if the input is𝑣(𝑡) = 5𝑢𝑠(𝑡)

System Dynamics 3.65 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

𝑠 + 32

System Dynamics 3.66 Fluid and Thermal Systems

HCM City Univ of Technology, Faculty of Mechanical Engineering Nguyen Tan Tien

Ngày đăng: 07/11/2017, 14:54

TỪ KHÓA LIÊN QUAN