Solution:
Magnitude of the mechanical components chosen from Table 1:
Mass,
m1 = 5 kg
m2 = 4 kg
Spring Compliance (Reciprocal of Spring Constant K),
k1 = 2000 N/m
k2 = 4000 N/m
k3 = 4000 N/m
Dashpot Damping Co-efficient,
c1 = 30 kg/s
c2 = 10 kg/s
c3 = 10 kg/s
Let the displacement and velocity for m1 be x1 & u1 and for m2 be x2 & u2 respectively.
Mathematical model of the mechanical system:
0 ) ( ) ( ) 0 ( ) 0 ( ) 0 ( ) ( ) (
2 3 1 2 2
0
2 2 3
0
1 2 1 2 2 2 2
= + +
(
+
(
+ +
} }
u c u u c x u k x x u u k u
dt
d
m
t t
Or, 0 ) ( ) ( ) ( ) ( ) (
2 3 1 2 2 2 3 1 2 2 2 2
= + + + x c x x c x k x x k x m
Or, 0 ) ( 10 ) ( 10 ) ( 4000 ) ( 4000 ) ( 4
2 1 2 2 1 2 2
= + + + x x x x x x x
0 ) ( ) ( ) 0 ( ) 0 ( ) 0 ( ) ( ) (
1 1 2 1 2
0
1 1 1
0
2 1 2 1 2 1 1
= + +
(
}
+ +
(
}
+ + u c u u c x u k x x u u k u
dt
d
m
t t
Or, 0 ) ( ) ( ) ( ) ( ) (
1 1 2 1 2 1 1 2 1 2 1 1
= + + + + x c x x c x k x x k x m
Or, 0 ) ( 30 ) ( 10 ) ( 2000 ) ( 4000 ) ( 5
1 2 1 1 2 1 1
= + + + + x x x x x x x
Please note that no force is shown on the mechanical system.
Considering ForceVoltage analogy between the Mechanical and Electrical systems,
Force (f) = Voltage (v)
Velocity (u) = Current (i)
Displacement (x) = Charge (q)
Inductor = Mass
L1 = m1
L2 = m2
Capacitor = Spring Constant (reciprocal of spring Compliance)
C1 = 1/k1
C2 = 1/k2
C3 = 1/k3
Resistor = Dashpot Damping Co-efficient
R1 = c1
R2 = c2
R3 = c3
Mathematical model of the equivalent electrical system:
0 ) ( ) ( ) 0 (
1
) 0 ( ) 0 ( ) (
1
) (
2 3 1 2 2
0
2 2
3 0
1 2 1 2
2
2 2
= + +
(
+ +
(
+ +
} }
i R i i R q i
C
q q i i
C
i
dt
d
L
t t
Or, 0 ) ( ) (
) ( ) (
) (
2 3 1 2 2
3
2
2
1 2
2 2
= + + +
+ q R q q R
C
q
C
q q
q L
0 ) ( ) ( ) 0 (
1
) 0 ( ) 0 ( ) (
1
) (
1 1 2 1 2
0
1 1
1
0
2 1 2 1
2
1 1
= + +
(
}
+ +
(
}
+ + i R i i R q i
C
q q i i
C
i
dt
d
L
t t
Or, 0 ) ( ) (
) ( ) (
) (
1 1 2 1 2
1
1
2
2 1
1 1
= + + +
+ q R q q R
C
q
C
q q
q L
Magnitude of the electrical components:
L1 = 5000 mH = 5 H
L2 = 4000 mH = 4 H
C1 = 500 F = 0.0005 F
C2 = 250 F = 0.00025 F
C3 = 250 F = 0.00025 F
R1 = 30
R2 = 10
R3 = 10
Equivalent electrical circuit,
Please note that as there is no force on the mechanical system, there is no voltage source for the
electrical circuit also.
MATLAB Simulink Model:
MATLAB Simulink Output:
Y axis: Displacement (m) X axis: Time (s)