0% found this document useful (0 votes)
22 views4 pages

Laplace Transform & DC Motor Modelling

This tutorial document covers modeling physical systems using Laplace transforms and transfer functions in MATLAB. It contains 4 problems: 1) Build Simulink models of a DC motor and compare responses with and without inductance. 2) Take the Laplace transforms of 5 functions. 3) Derive and simulate the transfer function of a mass-spring-damper system. 4) Derive and simulate the transfer function of a torsional vibration system.

Uploaded by

Sam Stide
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
22 views4 pages

Laplace Transform & DC Motor Modelling

This tutorial document covers modeling physical systems using Laplace transforms and transfer functions in MATLAB. It contains 4 problems: 1) Build Simulink models of a DC motor and compare responses with and without inductance. 2) Take the Laplace transforms of 5 functions. 3) Derive and simulate the transfer function of a mass-spring-damper system. 4) Derive and simulate the transfer function of a torsional vibration system.

Uploaded by

Sam Stide
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Control ENGD3113 2022/2023

Tutorial 2 – Modelling of physical systems and


introduction to Laplace transform
Dr. Yunda Yan

Objectives
The objective of this tutorial is to:

• learn how to build a mathematical model of physical systems in Simulink

• learn how to compute the Laplace transform and use its properties;

• understand the concept of transfer functions as a model to physical systems and


its implementation in Matlab.

• implement and analyse numerically transfer functions and time reponse using
Matlab.

Problem 1
1. Explain the working principle of a DC motor depicted in Figure 1.

Figure 1: Physical layout of a D.C. motor.

2. Explain the components and variables shown of the schematic diagram in Figure
2 and map them to the physical D.C. motor shown in Figure 1.

Page 1
Control ENGD3113 2022/2023

Figure 2: Equivalent schematic of the D.C. motor.

3. The mathematical model of the D.C. motor is represented by equations (1) and
(2).
Electrical Part:
dia Ra Kb 1
= − ia − ωm + ea (1)
dt La La La
Mechanical Part:
dωm Km b
= ia − ωm (2)
dt Im Im
where Kb is the back emf constant, Km is the motor constant and b is the friction
coefficient.

• Build a Simulink model of the D.C. motor, and simulate its behaviour for
the constant input voltage and the following values of the motor parameters:
Ra = 1, La = 1 mH, Im = 1 N.m.s2 /rad, b = 20 N.m.s/rad and Km = Kb =
5 N.m/A, ea = 10V .
• Build a Simulink model of the D.C. motor ignoring the inductance, simu-
late its behaviour and compare the results with the ones obtained in 3).

Problem 2
Find the Laplace transforms of the following functions:

• delayed unit-step function f1 (t) = u(t − T ), t ≥ 0.

• f2 (t) = δ (t) − 43 e−t u(t) + 31 e2t u(t), where δ (t) is the Dirac-delta function and
u(t) is the unit-step function.

• f3 (t) = e−at cos(ω t), t ≥ 0.

Page 2
Control ENGD3113 2022/2023

• f4 (t) = e−at sin(ω t), t ≥ 0.

• f5 (t) = e−2t u(t) + e−t cos(3t)u(t).

Problem 3
• Calculate the transfer function from ( f → x), assuming the input is the force
and the output is the mass’s displacement, of the rectilinear vibration system
described by the differential equation given below:

mẍ + cẋ + kx = f (t).

with x(0) = ẋ(0) = 0.

• Compare the calculated transfer function to the standard second order oscillatory
system transfer function,

Kωn2
G(s) = ,
s + 2ζ ωn s + ωn2
2

and evaluate parameters K, ωn , ζ as functions of the parameters m, c, k of the


physical system.

• By assuming numerical values for m = 0.5 kg, c = 5N.s/m, k = 300N/m calcu-


late K, ωn , ζ and define the transfer function of the system G(s) in Matlab.

• Simulate the step response of the system resulting from the sudden change of the
force f (t) using the command “step” in Matlab. Investigate how the response
changes for different values of c.

• Compare the final value of the response from the simulations to the final value
of the response obtained from the final value theorem.

Problem 4
• Calculate the transfer function from (τ → θ ), assuming the input is the torque
and the output is the angular position, of the torsional vibration system described
by the differential equation given below:

J θ̈ + cθ̇ = τ.

Page 3
Control ENGD3113 2022/2023

with θ (0) = θ̇ (0) = 0.

• Compare the calculated transfer function to the standard form,

K
G(s) = ,
s(sT + 1)

and evaluate parameters K, T as functions of the parameters J, c of the physical


system.

• By assuming numerical values for J = 0.0214, c = 0.086 (these values are com-
patible with the torsional laboratory apparatus) calculate K, T and define the
transfer function of the system in Matlab.

• Simulate the step response of the system resulting from a sudden change of
the torque using the command “step” in Matlab. Investigate how the response
changes for different J.

• Compare the final value of the response from the simulations to the final value
of the response obtained from the final value theorem.

Page 4

Common questions

Powered by AI

Simulating a rectilinear vibration system's step response in Matlab involves defining the system's transfer function parameters from the differential equation m\(\ddot{x}\)+c\(\dot{x}\)+kx=f(t). Utilizing Matlab, the 'step' function computes the system's time-domain response to an instantaneous input. Variations in the damping coefficient (c) critically affect dynamic behavior: increasing c reduces oscillations and quickens settling, leading to an overdamped response. Conversely, a low c promotes underdamping with significant oscillations. This simulation aids in identifying optimal damping for stability and performance, essential for control strategy refinement .

The Laplace transform is pivotal in control engineering as it facilitates the transformation of differential equations governing physical systems into algebraic equations, simplifying analysis and design. This transformation enables the handling of complex time-domain systems in the s-domain, allowing for the use of transfer functions and stability analysis. It streamlines the process of system response characterization, helping engineers design controllers and predict system behavior efficiently, especially in a Matlab environment .

The step response of a torsional vibration system indicates how quickly and accurately the system reaches a new equilibrium in response to a sudden torque change, influenced significantly by system inertia (J). When J increases, it leads to slower response and increased overshoot due to higher stored energy. Observing these responses informs control strategy design by dictating necessary damping or gain adjustments to achieve desired settling times and minimized overshoot, ensuring system stability and minimizing energy loss, thus optimizing performance for various operational scenarios .

Ignoring inductance in the Simulink model of a DC motor simplifies the mathematical model but may lead to inaccurate transient response predictions, particularly under high-frequency operations or rapid changes in input voltage. The omission can affect current surges and voltage stability assumptions. Analytically, this simplification should be justified by assessing the time constant formed by the motor’s inductance and comparing it to the system's operational time scales. If inductance is negligible relative to resistive and mechanical effects, the simplification might be valid, highlighting the need for critical evaluation in model validation .

The back emf constant (Kb) signifies the electromagnetic force counteracting the applied voltage in a DC motor, directly affecting speed and torque. As Kb increases, the voltage drop due to back emf increases with rotational speed, reducing current and thus torque, stabilizing speed against load variations. It ensures operational efficiency and responds dynamically to changes in load, thus affecting the motor's ability to maintain consistent performance under varying mechanical conditions .

In a rectilinear vibration system, the damping ratio (ζ) determines the rate at which oscillations decrease, impacting system stability and transient response. A higher ζ implies an overdamped system with no oscillations, while a ζ less than 1 results in an underdamped system with pronounced oscillations. The natural frequency (ωn) indicates the frequency at which the system would oscillate without damping, determining the system’s responsiveness to dynamic inputs. Together, ζ and ωn provide critical insights into system design and performance tuning, affecting settling time, overshoot, and resonance tendencies .

The friction coefficient (b) affects the opposing torque experienced in the mechanical movement of the DC motor, as represented by the mechanical part of the motor's differential equation. A higher b value increases the damping effect, reducing oscillations and settling time, but may also reduce velocity under constant input conditions. In Simulink, b impacts the system's time response characteristics, such as overshoot and stability, thereby influencing performance engineering decisions such as speed regulation and torque control .

The working principle of a DC motor involves the conversion of electrical energy into mechanical energy through electromagnetic interaction within the motor's components. This principle is applied in Simulink by building a mathematical model represented by electrical and mechanical differential equations. The electrical equation models the armature circuit's dynamics, while the mechanical equation models the rotational motion. The system’s behavior is simulated under specific parameter values to analyze the motor's response to varying inputs .

The final value theorem allows calculating the steady-state value of a system's output directly from its transfer function in the Laplace domain, provided the system is stable. In control engineering, it helps predict the system's long-term behavior and validate simulation results by ensuring that the expected end-state matches the predicted steady-state value. By comparing simulation outcomes to analytic calculations, engineers can validate model accuracy and detect discrepancies arising from assumptions or numerical methods used in the simulation process .

To derive the transfer function for a torsional vibration system from the differential equation J \(\ddot{\theta}\) + c \(\dot{\theta}\) = \(\tau\), you apply the Laplace transform, assuming initial conditions are zero (θ(0) = ˙θ(0) = 0). In the Laplace domain, this becomes Js²Θ(s) + csΘ(s) = Τ(s), leading to the transfer function G(s) = Θ(s)/Τ(s) = 1/(Js^2 + cs). This transfer function is crucial for analyzing system dynamics, predicting response to inputs, and designing controllers to meet desired performance specifications .

You might also like