Course 4 Modeling of Dynamic Systems in FD
Course 4 Modeling of Dynamic Systems in FD
1
Rotational Mechanical System TF
Units:
T(t): N-m (newton-meters); 𝜃 𝑡 : rad (radians); 𝜔 𝑡 : rad/s (radians/second); K: N-m/rad
(newton-meters/radian), J: kg-𝑚2 (kilograms-meter𝑠 2 - newton-meters-second𝑠 2 /radian).
2
3
𝜃2 𝑠
Ex.: Find the transfer function 𝑇 𝑠
. The rod is undergoing torsion. A torque is
applied at the left, and the displacement is measured at the right.
The system is considered as a lumped parameter system. So that the torsion acts
like a spring concentrated at one particular point in the rod, with an inertia 𝐽1 to the
left and an inertia 𝐽2 to the right . The damping inside the flexible shaft is negligible.
There are two degrees of freedom, since each inertia can be rotated while the other
is held still. Hence it will take two simultaneous equations to solve the system.
4
Free-body diagram of 𝐽1 :
5
Free-body diagram of 𝐽2 :
6
Summing torques:
𝐽1 𝑠 2 + 𝐷1 𝑠 + 𝐾 𝜃1 𝑠 − 𝐾𝜃2 𝑠 = 𝑇(𝑠)
−𝐾𝜃1 𝑠 + 𝐽2 𝑠 2 + 𝐷2 𝑠 + 𝐾 𝜃2 𝑠 = 0
TF:
𝜃2 (𝑠) 𝐾
=
𝑇(𝑠) ∆
𝐽1 𝑠 2 + 𝐷1 𝑠 + 𝐾 −𝐾
∆=
−𝐾 𝐽2 𝑠 2 + 𝐷2 𝑠 + 𝐾
7
TF for systems with gears
9
Mechanical impedances
driven by gears a. Rotational system driven by gears. Gears
driving a rotational inertia, spring, and viscous
damper.
2
𝑁2
𝐽𝑠 + 𝐷𝑠 + 𝐾 𝜃2 𝑠 = 𝑇1 (𝑠)
𝑁1
2
𝑁1 𝑁2
𝐽𝑠 + 𝐷𝑠 + 𝐾 𝜃 𝑠 = 𝑇1 (𝑠)
𝑁2 1 𝑁1
Then
2 2 2
𝑁1 𝑁1 𝑁1 𝑁1
𝐽 𝑠2 + 𝐷 𝑠+𝐾 𝜃 𝑠 = 𝑇1 (𝑠)
𝑁2 𝑁2 𝑁2 𝑁2 1
2
𝑁2
𝐽𝑒 𝑠 + 𝐷𝑒 𝑠 + 𝐾𝑒 𝜃2 𝑠 = 𝑇1 (𝑠)
𝑁1
𝑁2 2 𝑁2 2
𝐽𝑒 = 𝐽1 + 𝐽2 𝐷𝑒 = 𝐷1 + 𝐷2 𝐾𝑒 = 𝐾2
𝑁1 𝑁1
TF:
𝜽𝟐 (𝒔) 𝑵𝟐 /𝑵𝟏
𝑮 𝒔 = =
𝑻𝟏 (𝒔) 𝑱𝒆 𝒔𝟐 + 𝑫𝒆 𝒔 + 𝑲𝒆
In orderto eliminate gears with large radii, a gear train (see the figure) is used
to implement large gear ratios by cascading smaller gear ratios. Next to each
rotation, the angular displacement relative to 𝜃1 has been calculated. The
equivqlent gear ratio is the product of the individual gear ratios. So:
𝑁1 𝑁3 𝑁5
𝜃4 = 𝜃
𝑁2 𝑁4 𝑁6 1
Gears with loss
Ex.: Find the TF 𝜃1 (𝑠)/𝑇1 (𝑠).
We reflect all the impedance to the input shaft, 𝜃1 . The equation of motion:
𝐽𝑒 𝑠 2 + 𝐷𝑒 𝑠 𝜃1 𝑠 = 𝑇1 (𝑠)
𝑁1 2 𝑁1 𝑁3 2 𝑁1 2
𝐽𝑒 = 𝐽1 + 𝐽2 + 𝐽3 + 𝐽4 + 𝐽5 ; 𝐷𝑒 = 𝐷1 + 𝐷2
𝑁2 𝑁2 𝑁4 𝑁2
𝜽𝟏 (𝒔) 𝟏
TF: 𝑮 𝒔 = =
𝑻𝟏 (𝒔) 𝑱𝒆 𝒔𝟐 +𝑫𝒆 𝒔
𝑅𝑎 +𝐿𝑎 𝑠 𝑇𝑚 (𝑠)
+ 𝐾𝑏 𝑠𝜃𝑚 𝑠 = 𝐸𝑎 (𝑠) (d)
𝐾𝑡
We must now find 𝑇𝑚 𝑠 in term of 𝜃𝑚 (𝑠). The figure shows a typical equivalent
mechanical loading on a motor. 𝐽𝑚 is the equivalent inertia at the armature and
includes both the armature inertia and the load inertia reflected to the armature.
𝐷𝑚 is the equivalent viscous damping at the armature and includes both the
armature viscous damping and the load viscous damping reflected to the
armature. Then
𝑇𝑚 𝑠 = 𝐽𝑚 𝑠 2 + 𝐷𝑚 𝑠 𝜃𝑚 (𝑠) (e)
(e) → (d):
𝑅𝑎 + 𝐿𝑎 𝑠 𝐽𝑚 𝑠 2 + 𝐷𝑚 𝑠 𝜃𝑚 (𝑠)
+ 𝐾𝑏 𝑠𝜃𝑚 𝑠 = 𝐸𝑎 (𝑠) 17
𝐾𝑡
𝐿𝑎 , the armature inductance, is assumed to be small compared to the armature
resistance, 𝑅𝑎 , which is usual for a dc motor, so
𝑅𝑎
𝐽 𝑠 + 𝐷𝑚 + 𝐾𝑏 𝑠𝜃𝑚 𝑠 = 𝐸𝑎 (𝑠)
𝐾𝑡 𝑚
TF:
𝑲𝒕
𝜽𝒎 (𝒔) ൗ(𝑹 𝑱 )
𝒂 𝒎
𝑮 𝒔 = =
𝑬𝒂 (𝒔) 𝒔 𝒔 + 𝟏 𝑫 + 𝑲𝒕 𝑲𝒃
𝑱𝒎 𝒎 𝑹𝒂
𝜃𝑚 (𝑠) 𝐾
=
𝐸𝑎 (𝑠) 𝑠(𝑠 + 𝛼)
𝐾𝑡 1 𝐾𝑡 𝐾𝑏
with 𝐾 = ൗ(𝑅𝑎𝐽𝑚 ) 𝛼= 𝐷𝑚 +
𝐽𝑚 𝑅𝑎
Electrical constants.
Substituting (a) and (c) to (b) with 𝐿𝑎 = 0 gives
𝑅𝑎
𝑇 𝑠 + 𝐾𝑏 𝑠𝜃𝑚 𝑠 = 𝐸𝑎 (𝑠)
𝐾𝑡 𝑚
By the inverse Laplace transform, we get
19
𝑅𝑎 𝑑𝜃𝑚 (𝑡)
𝑇𝑚 𝑡 + 𝐾𝑏 = 𝐸𝑎 (𝑡)
𝐾𝑡 𝑑𝑡
𝑅𝑎
𝑇 𝑡 + 𝐾𝑏 𝜔𝑚 (𝑡) = 𝑒𝑎 (𝑡)
𝐾𝑡 𝑚
If a dc voltage 𝑒𝑎 is applied, the motor will turn at a constant angular velocity
𝜔𝑚 with a constant torque 𝑇𝑚 . So when a motor is operating at a steady
state with a dc input,
𝑅𝑎
𝑇𝑚 + 𝐾𝑏 𝜔𝑚 = 𝑒𝑎
𝐾𝑡 Torque – speed curve
or
𝐾𝑏 𝐾𝑡 𝐾𝑡
𝑇𝑚 = − 𝜔 + 𝑒
𝑅𝑎 𝑚 𝑅𝑎 𝑎
𝐾𝑡 𝑲𝒕 𝑻𝒔𝒕𝒂𝒍𝒍
𝑇𝑠𝑡𝑎𝑙𝑙 = 𝑒 → =
𝑅𝑎 𝑎 𝑹𝒂 𝑹𝒂
𝑒𝑎 𝒆𝒂
𝜔𝑛𝑜−𝑙𝑜𝑎𝑑 = → 𝑲𝒃 =
𝐾𝑏 𝝎𝒏𝒐−𝒍𝒐𝒂𝒅
𝐾𝑡
Electrical constants , 𝐾𝑏
𝑅𝑎
𝐾𝑡 𝑇𝑠𝑡𝑎𝑙𝑙 500 𝑒𝑎 100
From the curve: = = = 5 and 𝐾𝑏 = = =2
𝑅𝑎 𝑒𝑎 100 𝜔𝑛𝑜−𝑙𝑜𝑎𝑑 50
21
Then
𝐾𝑡 5ൗ
𝜃𝑚 (𝑠) ൗ(𝑅 𝐽 ) 0.417
= 𝑎 𝑚
= 12 =
𝐸𝑎 (𝑠) 𝑠 𝑠 + 1 𝐷 + 𝐾𝑡 𝐾𝑏 1 𝑠(𝑠 + 1.667)
𝑠 𝑠+ 10 + (5)(2)
𝐽𝑚 𝑚 𝑅𝑎 12
𝑁1 1
We know that = , so the TF:
𝑁2 10
𝜃𝐿 (𝑠) 0.0417
=
𝐸𝑎 (𝑠) 𝑠(𝑠 + 1.667)
22
HW
• Problems no. 33 and 44 Chapter 2 [Nise]
23