Skip to main content

Chapter 2 · 6 hours

Component Modeling

IOE past exam questions

Past questions and answers

41 questions set from this chapter, 3 of them more than once. Most asked first.

  • Asked 3 times
  • 2079 Bhadra · 8+2 marks
  • 2076 Asoj · 10 marks
  • 2070 Chaitra · 6 marks

Find transfer function for the following mechanical system considering displacement of mass M2 as output of the system. Also develop force current analogous circuit. [Figure: wall – spring K1 – mass M1 (force f(t) applied); M1 connected to mass M2 (displacement x) through spring K2 in parallel with a series combination of damper B and spring K3]

Answer

Assumptions: the masses have no friction with the ground; x1x_1 is the displacement of M1M_1, xx the displacement of M2M_2, and x3x_3 the displacement of the junction between damper BB and spring K3K_3 (B on the M1M_1 side).

Equations of motion (Laplace, zero initial conditions)

Mass M1M_1:

F(s)=(M1s2+K1)X1+K2(X1−X)+Bs(X1−X3)F(s) = (M_1s^2 + K_1)X_1 + K_2(X_1 - X) + Bs(X_1 - X_3)

Junction (x3x_3, massless):   Bs(X3−X1)+K3(X3−X)=0\;Bs(X_3 - X_1) + K_3(X_3 - X) = 0

Mass M2M_2:   M2s2X+K2(X−X1)+K3(X−X3)=0\;M_2s^2X + K_2(X - X_1) + K_3(X - X_3) = 0

Solving

The series BB–K3K_3 branch acts as an equivalent spring BsK3Bs+K3\dfrac{BsK_3}{Bs+K_3}. Together with K2K_2 in parallel, the coupling between the masses is

Kc(s)=K2+BK3sBs+K3=N(s)Bs+K3,N(s)=B(K2+K3)s+K2K3K_c(s) = K_2 + \frac{BK_3s}{Bs+K_3} = \frac{N(s)}{Bs+K_3}, \quad N(s) = B(K_2+K_3)s + K_2K_3

The two mass equations become (M1s2+K1+Kc)X1−KcX=F(M_1s^2+K_1+K_c)X_1 - K_cX = F and −KcX1+(M2s2+Kc)X=0-K_cX_1 + (M_2s^2+K_c)X = 0. So

X(s)F(s)=Kc(M1s2+K1+Kc)(M2s2+Kc)−Kc2=Kc(M1s2+K1)(M2s2+Kc)+KcM2s2\frac{X(s)}{F(s)} = \frac{K_c}{(M_1s^2+K_1+K_c)(M_2s^2+K_c) - K_c^2} = \frac{K_c}{(M_1s^2+K_1)(M_2s^2+K_c) + K_cM_2s^2}

Multiplying top and bottom by (Bs+K3)(Bs+K_3) and expanding:

X(s)F(s)=B(K2+K3)s+K2K3Δ(s)\frac{X(s)}{F(s)} = \frac{B(K_2+K_3)s + K_2K_3}{\Delta(s)} Δ(s)=BM1M2s5+K3M1M2s4+B[K1M2+(K2+K3)(M1+M2)]s3+K3[K1M2+K2(M1+M2)]s2+BK1(K2+K3)s+K1K2K3\begin{aligned} \Delta(s) ={}& BM_1M_2s^5 + K_3M_1M_2s^4 + B\left[K_1M_2 + (K_2+K_3)(M_1+M_2)\right]s^3 \\ &+ K_3\left[K_1M_2 + K_2(M_1+M_2)\right]s^2 + BK_1(K_2+K_3)s + K_1K_2K_3 \end{aligned}

Force–voltage (F-V) analogy

Force → voltage, mass → inductance L=ML = M, damper → resistance R=BR = B, spring → capacitance C=1/KC = 1/K, velocity → loop current. Mechanical parallel → electrical series and vice versa, so the coupling becomes C2C_2 in series with (R∥C3R \parallel C_3).

  +--[L1]--[C1]--+--[L2]--+
  |              |        |
 (f) i1->       C2   <-i2 |
  |              |        |
  |           +--+--+     |
  |           R     C3    |
  |           +--+--+     |
  |              |        |
  +--------------+--------+
 L1=M1, C1=1/K1, C2=1/K2, R=B, C3=1/K3, L2=M2

Force–current (F-I) analogy

Force → current source, mass → capacitance C=MC = M, damper → conductance (R=1/BR = 1/B), spring → inductance L=1/KL = 1/K, velocity → node voltage. Each displacement point becomes a node.

     v1                v3             v2
  +---o---+---[R]------o-----[L3]---+---o
  |   |   |                         |   |
  |   |   +----------[L2]-----------+   |
  |   |   |                             |
 (f)  C1  L1                            C2
  |   |   |                             |
  +---+---+-----------------------------+
            ground (reference)
 C1=M1, L1=1/K1, L2=1/K2, R=1/B, L3=1/K3, C2=M2

Node equation at v1v_1, for example: f=M1v˙1+K1 ⁣∫v1dt+K2 ⁣∫(v1−v2)dt+B(v1−v3)f = M_1\dot v_1 + K_1\!\int v_1dt + K_2\!\int(v_1-v_2)dt + B(v_1-v_3), which matches the force equation of M1M_1 with v=x˙v = \dot x.

  • Asked 2 times
  • 2075 Asoj · 4 marks
  • 2068 Chaitra · 6 marks

Find the transfer function of given circuit. [Figure: op-amp circuit; input Ei is applied through resistor R1 to the inverting (−) input and through capacitor C to the non-inverting (+) input; R2 connects the non-inverting input to ground; feedback resistor R1 from output to the inverting input; output Eo]

Answer

Assumption: ideal op-amp (infinite gain and input impedance), so V−=V+V_- = V_+ and no current enters the inputs.

        R1 (feedback)
     +---/\/\/---+
     |           |
Ei --+-/\/\/--(-)|\
        R1       | >--+-- Eo
Ei ---||--+---(+)|/
       C  |
         R2
          |
         GND

Non-inverting input

CC and R2R_2 form a voltage divider:

V+=EiR2R2+1sC=EisCR21+sCR2V_+ = E_i\frac{R_2}{R_2 + \dfrac{1}{sC}} = E_i\frac{sCR_2}{1 + sCR_2}

Inverting input (KCL)

Current through input R1R_1 = current through feedback R1R_1, with V−=V+V_- = V_+:

Ei−V+R1=V+−EoR1  ⇒  Eo=2V+−Ei\frac{E_i - V_+}{R_1} = \frac{V_+ - E_o}{R_1} \;\Rightarrow\; E_o = 2V_+ - E_i

Transfer function

Eo=Ei(2sCR21+sCR2−1)=Ei sCR2−1sCR2+1Eo(s)Ei(s)=sR2C−1sR2C+1=−1−sR2C1+sR2C\begin{aligned} E_o &= E_i\left(\frac{2sCR_2}{1+sCR_2} - 1\right) = E_i\,\frac{sCR_2 - 1}{sCR_2 + 1} \\ \frac{E_o(s)}{E_i(s)} &= \frac{sR_2C - 1}{sR_2C + 1} = -\frac{1 - sR_2C}{1 + sR_2C} \end{aligned}

This is an all-pass (phase-shift) network: ∣Eo/Ei∣=1|E_o/E_i| = 1 at all frequencies, and phase ϕ=180∘−2tan⁡−1(ωR2C)\phi = 180^\circ - 2\tan^{-1}(\omega R_2C), going from 180∘180^\circ at ω=0\omega = 0 to 0∘0^\circ at high frequency (90° at ω=1/R2C\omega = 1/R_2C). It has a zero at s=+1/R2Cs = +1/R_2C (non-minimum phase) and a pole at s=−1/R2Cs = -1/R_2C.

  • Asked 2 times
  • 2079 Bhadra · 8 marks
  • 2076 Chaitra · 8 marks

Find the transfer function θ2(s)/T(s) for the mechanical rotational system of figure below. Also draw the T-V and T-I analogy circuit of the system. [Figure: fixed wall – viscous damper D1 – inertia J1 (torque T(t) applied, angle θ1); J1 – damper D – torsional spring K – inertia J2 (angle θ2); J2 – damper D2 – fixed wall]

Answer

Assumptions: zero initial conditions; θ3\theta_3 is the angle of the junction between damper DD and spring KK (D on the J1J_1 side). The series DD–KK coupling is massless.

Equations (Laplace form)

Inertia J1J_1:

T(s)=(J1s2+D1s)θ1+Ds(θ1−θ3)T(s) = (J_1s^2 + D_1s)\theta_1 + Ds(\theta_1 - \theta_3)

Junction:   Ds(θ3−θ1)+K(θ3−θ2)=0\;Ds(\theta_3 - \theta_1) + K(\theta_3 - \theta_2) = 0

Inertia J2J_2:   (J2s2+D2s)θ2+K(θ2−θ3)=0\;(J_2s^2 + D_2s)\theta_2 + K(\theta_2 - \theta_3) = 0

Solving

Damper and spring in series transmit the same torque, so they act as one element with "stiffness"

Ke(s)=DsKDs+KK_e(s) = \frac{DsK}{Ds + K}

Then (J1s2+D1s+Ke)θ1−Keθ2=T(J_1s^2+D_1s+K_e)\theta_1 - K_e\theta_2 = T and −Keθ1+(J2s2+D2s+Ke)θ2=0-K_e\theta_1 + (J_2s^2+D_2s+K_e)\theta_2 = 0, giving

θ2T=Ke(J1s2+D1s)(J2s2+D2s)+Ke[(J1+J2)s2+(D1+D2)s]\frac{\theta_2}{T} = \frac{K_e}{(J_1s^2+D_1s)(J_2s^2+D_2s) + K_e\left[(J_1+J_2)s^2 + (D_1+D_2)s\right]}

Multiply top and bottom by (Ds+K)(Ds+K) and cancel the common ss:

θ2(s)T(s)=DKs[(Ds+K)(J1s+D1)(J2s+D2)+DK{(J1+J2)s+D1+D2}]\frac{\theta_2(s)}{T(s)} = \frac{DK}{s\left[(Ds+K)(J_1s+D_1)(J_2s+D_2) + DK\{(J_1+J_2)s + D_1 + D_2\}\right]}

Expanded:

θ2(s)T(s)=DKs[a3s3+a2s2+a1s+a0]\frac{\theta_2(s)}{T(s)} = \frac{DK}{s\left[a_3s^3 + a_2s^2 + a_1s + a_0\right]} a3=DJ1J2,a2=D(D1J2+D2J1)+KJ1J2a1=DD1D2+K[D(J1+J2)+D1J2+D2J1]a0=K[D(D1+D2)+D1D2]\begin{aligned} a_3 &= DJ_1J_2, \quad a_2 = D(D_1J_2 + D_2J_1) + KJ_1J_2 \\ a_1 &= DD_1D_2 + K\left[D(J_1+J_2) + D_1J_2 + D_2J_1\right] \\ a_0 &= K\left[D(D_1+D_2) + D_1D_2\right] \end{aligned}

Torque–voltage (T-V) analogy

Torque → voltage, JJ → LL, DD → RR, KK → C=1/KC = 1/K, angular velocity → loop current. The series DD–KK pair becomes R∥CR \parallel C shared by the two loops.

  +--[L1]--[R1]--+-----+--[L2]--[R2]--+
  |              |     |              |
 (T)  i1 ->      R     C     <- i2    |
  |              |     |              |
  +--------------+-----+--------------+
 L1=J1, R1=D1, R=D, C=1/K, L2=J2, R2=D2

Torque–current (T-I) analogy

Torque → current source, JJ → CC, DD → conductance (R=1/DR = 1/D), KK → L=1/KL = 1/K, angular velocity → node voltage.

      ω1          ω3           ω2
  +----o----+--[R]--o--[L]--+----o
  |    |    |               |    |
 (T)   C1   R1              C2   R2
  |    |    |               |    |
  +----+----+---------------+----+
              ground
 C1=J1, R1=1/D1, R=1/D, L=1/K, C2=J2, R2=1/D2
  • 2082 Baisakh · 8+2 marks

Obtain transfer function for the system as below considering displacement of mass M2 as output. Also develop Force-Voltage analogous circuit. [Figure: wall – spring K1 – mass m1 (force F applied, displacement x1, friction fc1 with ground) – spring K2 in series with damper B1 – mass m2 (displacement x2, friction fc2 with ground) – damper B2 – wall]

Answer

Assumptions: fc1f_{c1}, fc2f_{c2} are viscous friction coefficients of m1m_1, m2m_2 with the ground; x3x_3 is the displacement of the massless junction between K2K_2 and B1B_1 (K2K_2 on the m1m_1 side); zero initial conditions.

Equations of motion (Laplace form)

Mass m1m_1:

F(s)=(m1s2+fc1s+K1)X1+K2(X1−X3)F(s) = (m_1s^2 + f_{c1}s + K_1)X_1 + K_2(X_1 - X_3)

Junction:   K2(X3−X1)+B1s(X3−X2)=0\;K_2(X_3 - X_1) + B_1s(X_3 - X_2) = 0

Mass m2m_2:

0=(m2s2+fc2s+B2s)X2+B1s(X2−X3)0 = (m_2s^2 + f_{c2}s + B_2s)X_2 + B_1s(X_2 - X_3)

Solving

K2K_2 and B1B_1 in series carry the same force, so they act as one coupling element

Zc(s)=K2B1sK2+B1sZ_c(s) = \frac{K_2B_1s}{K_2 + B_1s}

Let Z1=m1s2+fc1s+K1Z_1 = m_1s^2 + f_{c1}s + K_1 and Z2=m2s2+(fc2+B2)sZ_2 = m_2s^2 + (f_{c2}+B_2)s. Then

[Z1+Zc−Zc−ZcZ2+Zc][X1X2]=[F0]\begin{bmatrix} Z_1 + Z_c & -Z_c \\ -Z_c & Z_2 + Z_c \end{bmatrix} \begin{bmatrix} X_1 \\ X_2 \end{bmatrix} = \begin{bmatrix} F \\ 0 \end{bmatrix} X2F=Zc(Z1+Zc)(Z2+Zc)−Zc2=ZcZ1Z2+Zc(Z1+Z2)\frac{X_2}{F} = \frac{Z_c}{(Z_1+Z_c)(Z_2+Z_c) - Z_c^2} = \frac{Z_c}{Z_1Z_2 + Z_c(Z_1 + Z_2)}

Substituting ZcZ_c, multiplying by (K2+B1s)(K_2 + B_1s) and cancelling the common factor ss:

X2(s)F(s)=K2B1(K2+B1s)(m1s2+fc1s+K1)(m2s+fc2+B2)+K2B1[(m1+m2)s2+(fc1+fc2+B2)s+K1]\frac{X_2(s)}{F(s)} = \frac{K_2B_1}{(K_2 + B_1s)(m_1s^2 + f_{c1}s + K_1)(m_2s + f_{c2} + B_2) + K_2B_1\left[(m_1+m_2)s^2 + (f_{c1}+f_{c2}+B_2)s + K_1\right]}

The denominator is a fourth-order polynomial; its leading term is B1m1m2s4B_1m_1m_2s^4 and its constant term is K1K2(B1+B2+fc2)K_1K_2(B_1 + B_2 + f_{c2}).

Force–voltage (F-V) analogy

Force → voltage, mass → LL, friction/damper → RR, spring → C=1/KC = 1/K, velocity → loop current. Mechanical series (K2K_2, B1B_1) → electrical parallel (C2∥RB1C_2 \parallel R_{B1}), shared by both loops.

  +--[L1]--[Rf1]--[C1]--+--[L2]--[Rf2]--[RB2]--+
  |                     |                      |
 (F)  i1 ->          +--+--+          <- i2    |
  |                  C2   RB1                  |
  |                  +--+--+                   |
  |                     |                      |
  +---------------------+----------------------+
 L1=m1, Rf1=fc1, C1=1/K1, C2=1/K2, RB1=B1
 L2=m2, Rf2=fc2, RB2=B2

Loop equations, with i=x˙i = \dot x:

F=L1di1dt+Rf1i1+1C1∫i1dt+vp0=L2di2dt+(Rf2+RB2)i2−vp\begin{aligned} F &= L_1\frac{di_1}{dt} + R_{f1}i_1 + \frac{1}{C_1}\int i_1dt + v_p \\ 0 &= L_2\frac{di_2}{dt} + (R_{f2} + R_{B2})i_2 - v_p \end{aligned}

where vpv_p is the voltage across C2∥RB1C_2 \parallel R_{B1} carrying (i1−i2)(i_1 - i_2) in total; these match the mechanical equations.

  • 2081 Bhadra · 6+2 marks

Obtain transfer function for the following system considering angular position of object second (θ2) as output. Also develop Torque-Current Analogous Circuit. [Figure: torque T(t) applied to inertia J1 (angle θ1(t)) which has viscous friction D1 to the fixed frame; J1 is connected through a shaft of torsional stiffness K to inertia J2 (angle θ2(t)), which has viscous friction D2 to the fixed frame]

Answer

Assumptions: zero initial conditions; D1D_1, D2D_2 are viscous friction coefficients to the fixed frame.

Equations (Laplace form)

Inertia J1J_1:

T(s)=(J1s2+D1s+K)θ1−Kθ2T(s) = (J_1s^2 + D_1s + K)\theta_1 - K\theta_2

Inertia J2J_2:

0=−Kθ1+(J2s2+D2s+K)θ20 = -K\theta_1 + (J_2s^2 + D_2s + K)\theta_2

Solving (Cramer's rule)

Δ=(J1s2+D1s+K)(J2s2+D2s+K)−K2\Delta = (J_1s^2 + D_1s + K)(J_2s^2 + D_2s + K) - K^2 θ2=K T(s)Δ\theta_2 = \frac{K\,T(s)}{\Delta}

Expanding Δ\Delta:

Δ=J1J2s4+(J1D2+J2D1)s3+[K(J1+J2)+D1D2]s2+K(D1+D2)s\Delta = J_1J_2s^4 + (J_1D_2 + J_2D_1)s^3 + \left[K(J_1+J_2) + D_1D_2\right]s^2 + K(D_1+D_2)s θ2(s)T(s)=Ks[J1J2s3+(J1D2+J2D1)s2+{K(J1+J2)+D1D2}s+K(D1+D2)]\frac{\theta_2(s)}{T(s)} = \frac{K}{s\left[J_1J_2s^3 + (J_1D_2 + J_2D_1)s^2 + \{K(J_1+J_2) + D_1D_2\}s + K(D_1+D_2)\right]}

The factor ss in the denominator shows that a constant torque makes the whole system rotate continuously (angle grows with time).

Torque–current (T-I) analogy

Torque → current source, inertia JJ → capacitance C=JC = J, viscous friction DD → conductance (R=1/DR = 1/D), stiffness KK → inductance L=1/KL = 1/K, angular velocity ω\omega → node voltage. Each inertia is a node.

      ω1                 ω2
  +----o----+---[L]---+----o
  |    |    |         |    |
 (T)   C1   R1        C2   R2
  |    |    |         |    |
  +----+----+---------+----+
          ground
 C1=J1, R1=1/D1, L=1/K, C2=J2, R2=1/D2

Node equations:

T=C1dω1dt+ω1R1+1L∫(ω1−ω2)dt0=C2dω2dt+ω2R2+1L∫(ω2−ω1)dt\begin{aligned} T &= C_1\frac{d\omega_1}{dt} + \frac{\omega_1}{R_1} + \frac{1}{L}\int(\omega_1 - \omega_2)dt \\ 0 &= C_2\frac{d\omega_2}{dt} + \frac{\omega_2}{R_2} + \frac{1}{L}\int(\omega_2 - \omega_1)dt \end{aligned}

which match the mechanical equations with ω=θ˙\omega = \dot\theta.

  • 2081 Baisakh · 2+8 marks

The given mechanical system has force F(t) as input and X1 and X2 as displacement output. Draw the equivalent F-V analogous electrical circuit and determine the transfer function with X1 as output. [Figure: wall – spring K1 = 1 N/m – mass M1 = 1 kg (displacement X1(t), friction with ground B1 = 2 Ns/m); M1 connected to mass M2 = 1 kg (displacement X2(t)) through spring K2 = 1 N/m in parallel with damper B2 = 1 Ns/m in series with spring K3 = 2 N/m; force F(t) applied to M2; M2 has friction with ground B3 = 1 Ns/m]

Answer

Assumptions: damper B2B_2 and spring K3K_3 are in series and this pair is in parallel with K2K_2 between M1M_1 and M2M_2; X3X_3 is the junction between B2B_2 and K3K_3; zero initial conditions. Values: M1=M2=1M_1 = M_2 = 1 kg, K1=K2=1K_1 = K_2 = 1 N/m, K3=2K_3 = 2 N/m, B1=2B_1 = 2, B2=1B_2 = 1, B3=1B_3 = 1 N·s/m.

F-V analogous circuit

Force → voltage, mass → LL, damper → RR, spring → C=1/KC = 1/K, velocity → loop current. The coupling (K2K_2 parallel to B2B_2–K3K_3 series) becomes C2C_2 in series with (R2∥C3R_2 \parallel C_3).

  +--[L1]--[R1]--[C1]--+--[L2]--[R3]--+
  |                    |              |
  |     <- i1         C2      i2 ->  (F)
  |                    |              |
  |                 +--+--+           |
  |                 R2    C3          |
  |                 +--+--+           |
  |                    |              |
  +--------------------+--------------+
 L1=M1=1 H, R1=B1=2 Ω, C1=1/K1=1 F
 C2=1/K2=1 F, R2=B2=1 Ω, C3=1/K3=0.5 F
 L2=M2=1 H, R3=B3=1 Ω, source = F(t)

Loop currents i1=x˙1i_1 = \dot x_1, i2=x˙2i_2 = \dot x_2; the shared branch carries (i2−i1)(i_2 - i_1).

Equations of motion

Equivalent of series B2B_2–K3K_3: B2sK3B2s+K3=2ss+2\dfrac{B_2sK_3}{B_2s+K_3} = \dfrac{2s}{s+2}. Total coupling:

Kc(s)=1+2ss+2=3s+2s+2K_c(s) = 1 + \frac{2s}{s+2} = \frac{3s+2}{s+2} M1:(s2+2s+1)X1+Kc(X1−X2)=0M2:(s2+s)X2+Kc(X2−X1)=F(s)\begin{aligned} M_1: &\quad (s^2 + 2s + 1)X_1 + K_c(X_1 - X_2) = 0 \\ M_2: &\quad (s^2 + s)X_2 + K_c(X_2 - X_1) = F(s) \end{aligned}

Transfer function X1/FX_1/F

From the first equation, X2=s2+2s+1+KcKcX1X_2 = \dfrac{s^2+2s+1+K_c}{K_c}X_1. Substituting in the second:

X1F=Kc(s2+2s+1)(s2+s)+Kc(2s2+3s+1)\frac{X_1}{F} = \frac{K_c}{(s^2+2s+1)(s^2+s) + K_c(2s^2+3s+1)}

Using (s2+2s+1)(s2+s)=s(s+1)3(s^2+2s+1)(s^2+s) = s(s+1)^3 and 2s2+3s+1=(2s+1)(s+1)2s^2+3s+1 = (2s+1)(s+1), then multiplying top and bottom by (s+2)(s+2):

X1F=3s+2s(s+2)(s+1)3+(3s+2)(2s+1)(s+1)=3s+2(s+1)[s(s+2)(s+1)2+(3s+2)(2s+1)]=3s+2(s+1)(s4+4s3+11s2+9s+2)\begin{aligned} \frac{X_1}{F} &= \frac{3s+2}{s(s+2)(s+1)^3 + (3s+2)(2s+1)(s+1)} \\ &= \frac{3s+2}{(s+1)\left[s(s+2)(s+1)^2 + (3s+2)(2s+1)\right]} \\ &= \frac{3s+2}{(s+1)(s^4 + 4s^3 + 11s^2 + 9s + 2)} \end{aligned}

Answer:

X1(s)F(s)=3s+2s5+5s4+15s3+20s2+11s+2\frac{X_1(s)}{F(s)} = \frac{3s+2}{s^5 + 5s^4 + 15s^3 + 20s^2 + 11s + 2}
  • 2080 Bhadra · 6 marks

Find the transfer function of armature controlled dc motor and also draw the block diagram of same. [Figure: armature circuit with supply va, armature current ia, resistance Ra and inductance La feeding the motor; field winding Lf with constant field voltage vf (field current if); motor shaft (angle θ) drives a load with inertia J and friction B]

Answer

In an armature-controlled DC motor the field current ifi_f is kept constant (constant vfv_f) and the speed/position is controlled by changing the armature voltage vav_a.

Assumptions and symbols

  • Ra,LaR_a, L_a: armature resistance and inductance; iai_a: armature current
  • ebe_b: back emf; KbK_b: back-emf constant
  • TmT_m: motor torque; KtK_t: torque constant
  • J,BJ, B: inertia and viscous friction of motor plus load; θ\theta: shaft angle
  • Flux ϕ∝if\phi \propto i_f is constant, so torque and back emf are linear.

Governing equations

  1. Armature circuit (KVL):
va=Raia+Ladiadt+ebv_a = R_a i_a + L_a \frac{di_a}{dt} + e_b
  1. Back emf: eb=Kbdθdte_b = K_b \dfrac{d\theta}{dt}
  2. Torque: Tm=Kt iaT_m = K_t\, i_a
  3. Load (Newton's law): Tm=Jd2θdt2+BdθdtT_m = J\dfrac{d^2\theta}{dt^2} + B\dfrac{d\theta}{dt}

Laplace transform (zero initial conditions)

Ia(s)=Va(s)−Eb(s)Ra+sLaTm(s)=KtIa(s)θ(s)=Tm(s)s(Js+B)Eb(s)=Kb s θ(s)\begin{aligned} I_a(s) &= \frac{V_a(s) - E_b(s)}{R_a + sL_a} \\ T_m(s) &= K_t I_a(s) \\ \theta(s) &= \frac{T_m(s)}{s(Js + B)} \\ E_b(s) &= K_b\, s\,\theta(s) \end{aligned}

Block diagram

 Va +    ┌──────┐ Ia ┌──┐ Tm ┌─────┐ ω  ┌───┐ θ
 ──►(Σ)─►│  1   ├───►│Kt├───►│  1  ├─┬─►│1/s├─►
    ▲ -  │Ra+sLa│    └──┘    │Js+B │ │  └───┘
    │    └──────┘            └─────┘ │
    │ Eb        ┌────┐               │
    └───────────┤ Kb │◄──────────────┘
                └────┘

The back emf forms an internal negative feedback loop on speed ω=sθ\omega = s\theta.

Transfer function

Reducing the inner loop (forward path G=Kt(Ra+sLa)(Js+B)G = \dfrac{K_t}{(R_a+sL_a)(Js+B)}, feedback KbK_b):

ω(s)Va(s)=Kt(Ra+sLa)(Js+B)+KtKb\frac{\omega(s)}{V_a(s)} = \frac{K_t}{(R_a + sL_a)(Js + B) + K_t K_b}

Multiplying by the integrator 1/s1/s:

θ(s)Va(s)=Kts[(Ra+sLa)(Js+B)+KtKb]\frac{\theta(s)}{V_a(s)} = \frac{K_t}{s\left[(R_a + sL_a)(Js + B) + K_t K_b\right]}

If LaL_a is negligible, this becomes θ(s)Va(s)=Kms(1+sτm)\dfrac{\theta(s)}{V_a(s)} = \dfrac{K_m}{s(1 + s\tau_m)} with Km=KtRaB+KtKbK_m = \dfrac{K_t}{R_aB + K_tK_b} and τm=RaJRaB+KtKb\tau_m = \dfrac{R_aJ}{R_aB + K_tK_b}, a type-1 second-order system.

  • 2080 Bhadra · 4 marks

Develop F-V and F-I analogy circuit of the mechanical system shown below. [Figure: force f applied to mass M (displacement X1, on rollers on ground); M connected through spring K1 to a point of displacement X2; from that point, damper B2 in parallel with spring K2 connect to the fixed wall]

Answer

Take x1x_1 (mass) and x2x_2 (junction of K1K_1, K2K_2, B2B_2) as the two displacements. The junction has no mass.

Equations of motion

Mx¨1+K1(x1−x2)=fK1(x2−x1)+B2x˙2+K2x2=0\begin{aligned} M\ddot{x}_1 + K_1(x_1 - x_2) &= f \\ K_1(x_2 - x_1) + B_2\dot{x}_2 + K_2 x_2 &= 0 \end{aligned}

F-V analogy (mesh/loop)

f→ef \to e, M→LM \to L, B→RB \to R, K→1/CK \to 1/C, x˙→i\dot x \to i. Each displacement gives one mesh; K1K_1 is common to both.

Ldi1dt+1C1∫(i1−i2) dt=e1C1∫(i2−i1) dt+R2i2+1C2∫i2 dt=0\begin{aligned} L\frac{di_1}{dt} + \frac{1}{C_1}\int (i_1 - i_2)\,dt &= e \\ \frac{1}{C_1}\int (i_2 - i_1)\,dt + R_2 i_2 + \frac{1}{C_2}\int i_2\,dt &= 0 \end{aligned}

with L=ML = M, C1=1/K1C_1 = 1/K_1, R2=B2R_2 = B_2, C2=1/K2C_2 = 1/K_2.

      L = M               R2 = B2
  +---LLLL----+---------/\/\/----+
  |           |                  |
 (~) e      C1=1/K1           C2=1/K2
  |    i1     |       i2         |
  +-----------+------------------+

F-I analogy (node)

f→if \to i, M→CM \to C, B→1/RB \to 1/R, K→1/LK \to 1/L, x˙→v\dot x \to v. Each displacement gives one node; K1K_1 joins the nodes.

Cdv1dt+1L1∫(v1−v2) dt=i1L1∫(v2−v1) dt+v2R2+1L2∫v2 dt=0\begin{aligned} C\frac{dv_1}{dt} + \frac{1}{L_1}\int (v_1 - v_2)\,dt &= i \\ \frac{1}{L_1}\int (v_2 - v_1)\,dt + \frac{v_2}{R_2} + \frac{1}{L_2}\int v_2\,dt &= 0 \end{aligned}

with C=MC = M, L1=1/K1L_1 = 1/K_1, R2=1/B2R_2 = 1/B_2, L2=1/K2L_2 = 1/K_2.

    v1          L1=1/K1         v2
  +---+--------LLLL-------+--------+
  |   |                   |        |
 (i) C=M              R2=1/B2  L2=1/K2
  |   |                   |        |
  +---+-------------------+--------+  ground
  • 2080 Baisakh · 8 marks

Find the transfer X2(s)/F(s) for the given mechanical system. Also develop F-I analogy circuit. [Figure: wall – spring K1 – mass M (force f applied, displacement X1, friction fc1 with ground) – spring K2 in series with damper B1 – mass M (displacement X2, friction fc2 with ground) – damper B2 – wall]

Answer

Let x1x_1, x2x_2 be the displacements of the two masses (M1M_1, M2M_2) and x3x_3 the displacement of the massless junction between K2K_2 and B1B_1 (spring on the M1M_1 side).

Differential equations

M1x¨1+fc1x˙1+K1x1+K2(x1−x3)=fK2(x3−x1)+B1(x˙3−x˙2)=0M2x¨2+(fc2+B2)x˙2+B1(x˙2−x˙3)=0\begin{aligned} M_1\ddot x_1 + f_{c1}\dot x_1 + K_1 x_1 + K_2(x_1 - x_3) &= f \\ K_2(x_3 - x_1) + B_1(\dot x_3 - \dot x_2) &= 0 \\ M_2\ddot x_2 + (f_{c2} + B_2)\dot x_2 + B_1(\dot x_2 - \dot x_3) &= 0 \end{aligned}

Eliminating the junction

From the second equation, X3=K2X1+B1sX2K2+B1sX_3 = \dfrac{K_2X_1 + B_1sX_2}{K_2 + B_1 s}. The force through the series spring–damper is then

K2(X1−X3)=Keq(X1−X2),Keq(s)=K2B1sK2+B1sK_2(X_1 - X_3) = K_{eq}(X_1 - X_2), \qquad K_{eq}(s) = \frac{K_2 B_1 s}{K_2 + B_1 s}

So, writing a=M1s2+fc1s+K1a = M_1s^2 + f_{c1}s + K_1 and b=M2s2+(fc2+B2)sb = M_2s^2 + (f_{c2}+B_2)s:

(a+Keq)X1−KeqX2=F−KeqX1+(b+Keq)X2=0\begin{aligned} (a + K_{eq})X_1 - K_{eq}X_2 &= F \\ -K_{eq}X_1 + (b + K_{eq})X_2 &= 0 \end{aligned}

Transfer function

X2(s)F(s)=Keq(a+Keq)(b+Keq)−Keq2=Keqab+Keq(a+b)\frac{X_2(s)}{F(s)} = \frac{K_{eq}}{(a+K_{eq})(b+K_{eq}) - K_{eq}^2} = \frac{K_{eq}}{ab + K_{eq}(a+b)}

Multiplying through by (K2+B1s)(K_2 + B_1 s) and cancelling the common factor ss:

X2(s)F(s)=B1K2(K2+B1s)(M1s2+fc1s+K1)(M2s+fc2+B2)+B1K2[(M1+M2)s2+(fc1+fc2+B2)s+K1]\frac{X_2(s)}{F(s)} = \frac{B_1K_2}{(K_2 + B_1s)(M_1s^2 + f_{c1}s + K_1)(M_2s + f_{c2} + B_2) + B_1K_2\left[(M_1+M_2)s^2 + (f_{c1}+f_{c2}+B_2)s + K_1\right]}

The denominator is a 4th-order polynomial in ss.

F-I analogy

f→if \to i, M→CM \to C, friction/damper B→B \to conductance 1/R1/R, K→1/LK \to 1/L, velocity →\to node voltage. Nodes v1,v3,v2v_1, v_3, v_2 correspond to x˙1,x˙3,x˙2\dot x_1, \dot x_3, \dot x_2.

MechanicalElectrical (F-I)
M1M_1, M2M_2C1=M1C_1 = M_1, C2=M2C_2 = M_2 (node to ground)
fc1f_{c1}, fc2f_{c2}, B2B_2R=1/fc1R = 1/f_{c1}, 1/fc21/f_{c2}, 1/B21/B_2 (node to ground)
K1K_1L1=1/K1L_1 = 1/K_1 (node 1 to ground)
K2K_2L2=1/K2L_2 = 1/K_2 (node 1 to node 3)
B1B_1RB1=1/B1R_{B1} = 1/B_1 (node 3 to node 2)
  v1               v3               v2
 +--+---+----+-LLLL-+--/\/\/--+--+----+----+
 |  |   |    | 1/K2   1/B1    |  |    |    |
(i) C1 R1   L1               C2 R2   R3    |
 |  =M1 =1/ =1/              =M2 =1/ =1/B2 |
 |  |  fc1  K1                |  fc2  |    |
 +--+---+----+----------------+--+----+----+ gnd

Node equation at node 1: C1dv1dt+v1R1+1L1∫v1dt+1L2∫(v1−v3)dt=iC_1\dfrac{dv_1}{dt} + \dfrac{v_1}{R_1} + \dfrac{1}{L_1}\int v_1dt + \dfrac{1}{L_2}\int (v_1 - v_3)dt = i, and similarly for nodes 3 and 2 (matching the three mechanical equations).

  • 2080 Baisakh · 4 marks

Find transfer function. [Figure: op-amp circuit; input ein applied through R1 to one op-amp input, with capacitor C connected from that input to the op-amp output; the other op-amp input is connected through R2 to the common (ground) line; output eout taken at the op-amp output]

Answer

The circuit is an inverting integrator (the input, through R1R_1, and the capacitor both go to the inverting input). R2R_2 only balances the input bias current; for an ideal op-amp it carries no current, so the non-inverting input is at 0 V.

Assumptions (ideal op-amp)

  • Infinite input impedance: no current into the op-amp terminals.
  • Virtual ground: V−=V+=0V_- = V_+ = 0 (no current in R2R_2, so no drop across it).

Derivation

Current through R1R_1:

I1(s)=Ein(s)−0R1I_1(s) = \frac{E_{in}(s) - 0}{R_1}

This whole current flows through CC (impedance 1/sC1/sC):

I1(s)=0−Eout(s)1/(sC)=−sC Eout(s)I_1(s) = \frac{0 - E_{out}(s)}{1/(sC)} = -sC\,E_{out}(s)

Equating:

Ein(s)R1=−sC Eout(s)\frac{E_{in}(s)}{R_1} = -sC\,E_{out}(s)

Transfer function

Eout(s)Ein(s)=−1sR1C\frac{E_{out}(s)}{E_{in}(s)} = -\frac{1}{sR_1C}

In time domain, eout(t)=−1R1C∫ein dte_{out}(t) = -\dfrac{1}{R_1C}\displaystyle\int e_{in}\,dt. So the circuit is an integrator with time constant R1CR_1C and a sign inversion; it gives a pole at s=0s = 0 and is used to realise the integral (I) action of a controller. Usually R2R_2 is chosen equal to R1R_1 to reduce output offset.

  • 2078 Bhadra · 8 marks

Find the transfer function θL(s)/T(s) of the mechanical rotational system shown below. [Figure: torque T(t) applied to gear N1 = 11, which drives gear N2 = 33; N2 is on a shaft with inertia 1 kg-m² and a viscous damper 2 N-m-s/rad, followed by a torsional spring 3 N-m/rad connected to gear N3 = 50; N3 drives gear N4 = 10, whose shaft (angle θL(t)) has viscous damping 0.04 N-m-s/rad to the frame]

Answer

Reflect everything to one shaft. Mechanical impedances move across a gear pair by (Ndest/Nsource)2(N_{dest}/N_{source})^2, and torques by Ndest/NsourceN_{dest}/N_{source}.

Step 1: Reflect the input torque to the N2N_2 shaft

T2=T N2N1=T⋅3311=3TT_2 = T\,\frac{N_2}{N_1} = T\cdot\frac{33}{11} = 3T

Step 2: Reflect the load damping to the N3N_3 shaft

Deq=0.04(N3N4)2=0.04(5010)2=1 N-m-s/radD_{eq} = 0.04\left(\frac{N_3}{N_4}\right)^2 = 0.04\left(\frac{50}{10}\right)^2 = 1\ \text{N-m-s/rad}

Step 3: Equivalent system

Let θ2\theta_2 = angle of the N2N_2 shaft (with J=1J = 1, D=2D = 2) and θ3\theta_3 = angle of gear N3N_3 (end of the spring, K=3K = 3, carrying the reflected damping 1):

 3T ─► [J=1, D=2] ──spring K=3── (θ3) ── D=1 ─► frame
        θ2

Equations (Laplace):

(s2+2s+3) θ2−3 θ3=3T−3 θ2+(s+3) θ3=0\begin{aligned} (s^2 + 2s + 3)\,\theta_2 - 3\,\theta_3 &= 3T \\ -3\,\theta_2 + (s + 3)\,\theta_3 &= 0 \end{aligned}

Step 4: Solve

From the second equation θ2=(s+3)3θ3\theta_2 = \dfrac{(s+3)}{3}\theta_3. Substituting:

[(s2+2s+3)(s+3)3−3]θ3=3T[(s3+5s2+9s+9)−9]θ3=9Tθ3(s)T(s)=9s(s2+5s+9)\begin{aligned} \left[\frac{(s^2+2s+3)(s+3)}{3} - 3\right]\theta_3 &= 3T \\ \left[(s^3 + 5s^2 + 9s + 9) - 9\right]\theta_3 &= 9T \\ \frac{\theta_3(s)}{T(s)} &= \frac{9}{s(s^2 + 5s + 9)} \end{aligned}

Step 5: Output shaft

θL=θ3 N3N4=5 θ3\theta_L = \theta_3\,\frac{N_3}{N_4} = 5\,\theta_3

Answer:

θL(s)T(s)=45s(s2+5s+9)\frac{\theta_L(s)}{T(s)} = \frac{45}{s(s^2 + 5s + 9)}

The pole at s=0s=0 appears because no spring ties the system to the frame, so a constant torque gives a steadily increasing angle.

  • 2078 Kartik · 6+2 marks

Consider displacement of mass M2 as output and find transfer function of the mechanical system as in figure. Also find F-I analogous circuit. [Figure: mass M1 is a hollow frame resting on the ground, with force f applied; mass M2 sits inside M1 and is connected to M1 through spring k1; M1 is connected to the right wall through spring k2]

Answer

Let x1x_1 = displacement of the outer frame M1M_1 and x2x_2 = displacement of the inner mass M2M_2. Assumption: no friction (none is shown) between M1M_1 and ground or between M2M_2 and M1M_1.

Free body equations

  • M1M_1: force ff acts; spring k2k_2 to the wall resists by k2x1k_2x_1; spring k1k_1 pulls by k1(x1−x2)k_1(x_1 - x_2).
  • M2M_2: only spring k1k_1 acts, k1(x2−x1)k_1(x_2 - x_1).
M1x¨1+k1(x1−x2)+k2x1=fM2x¨2+k1(x2−x1)=0\begin{aligned} M_1\ddot x_1 + k_1(x_1 - x_2) + k_2x_1 &= f \\ M_2\ddot x_2 + k_1(x_2 - x_1) &= 0 \end{aligned}

Laplace form

(M1s2+k1+k2)X1−k1X2=F−k1X1+(M2s2+k1)X2=0\begin{aligned} (M_1s^2 + k_1 + k_2)X_1 - k_1X_2 &= F \\ -k_1X_1 + (M_2s^2 + k_1)X_2 &= 0 \end{aligned}

From the second equation, X1=M2s2+k1k1X2X_1 = \dfrac{M_2s^2 + k_1}{k_1}X_2. Substituting:

[(M1s2+k1+k2)(M2s2+k1)−k12]X2=k1F\left[(M_1s^2 + k_1 + k_2)(M_2s^2 + k_1) - k_1^2\right]X_2 = k_1F

Transfer function

X2(s)F(s)=k1M1M2s4+[M1k1+M2(k1+k2)]s2+k1k2\frac{X_2(s)}{F(s)} = \frac{k_1}{M_1M_2s^4 + \left[M_1k_1 + M_2(k_1 + k_2)\right]s^2 + k_1k_2}

With no damping the poles lie on the jωj\omega axis (undamped oscillation at two natural frequencies).

F-I analogous circuit

f→if\to i, M→CM\to C, k→1/Lk\to 1/L, velocity →\to node voltage. Note that M2M_2 is referred to the ground reference (absolute velocity), so its capacitor goes to ground.

   v1          L1 = 1/k1        v2
  +----+---+-----LLLL-------+------+
  |    |   |                |      |
 (i)  C1  L2               C2      |
  |  =M1 =1/k2             =M2     |
  |    |   |                |      |
  +----+---+----------------+------+ gnd
C1dv1dt+1L2∫v1dt+1L1∫(v1−v2)dt=iC2dv2dt+1L1∫(v2−v1)dt=0\begin{aligned} C_1\frac{dv_1}{dt} + \frac{1}{L_2}\int v_1dt + \frac{1}{L_1}\int(v_1 - v_2)dt &= i \\ C_2\frac{dv_2}{dt} + \frac{1}{L_1}\int(v_2 - v_1)dt &= 0 \end{aligned}
  • 2078 Kartik · 4 marks

Draw the block diagram of circuit shown below. [Figure: input Vi; series R1 to a node; capacitor C1 from that node to ground; series R2 to the output node; capacitor C2 from output node to ground; Vo across C2; mesh currents i1(t), i2(t)]

Answer

Write one equation per element so each gives one block, then join the blocks.

Equations (Laplace domain)

Let V1V_1 = voltage across C1C_1.

I1(s)=Vi(s)−V1(s)R1V1(s)=I1(s)−I2(s)sC1I2(s)=V1(s)−Vo(s)R2Vo(s)=I2(s)sC2\begin{aligned} I_1(s) &= \frac{V_i(s) - V_1(s)}{R_1} \\ V_1(s) &= \frac{I_1(s) - I_2(s)}{sC_1} \\ I_2(s) &= \frac{V_1(s) - V_o(s)}{R_2} \\ V_o(s) &= \frac{I_2(s)}{sC_2} \end{aligned}

Block diagram

 Vi─►(Σ)►[1/R1]─►(Σ)►[1/sC1]─┬─►(Σ)►[1/R2]─┬─►[1/sC2]─┬─►Vo
      ▲-          ▲-         │   ▲-        │          │
      │           └── I2 ────┼───┼─────────┘          │
      └── V1 ────────────────┘   └── Vo ──────────────┘
  • V1V_1 is fed back (negative) to the first summer, I2I_2 to the second and VoV_o to the third.

Resulting transfer function

Reducing the three interacting loops (or by Mason's rule) gives

Vo(s)Vi(s)=1R1R2C1C2s2+(R1C1+R1C2+R2C2)s+1\frac{V_o(s)}{V_i(s)} = \frac{1}{R_1R_2C_1C_2s^2 + (R_1C_1 + R_1C_2 + R_2C_2)s + 1}

The term R1C2R_1C_2 shows the loading of the second stage on the first; the result is not simply the product 1(1+sR1C1)(1+sR2C2)\dfrac{1}{(1+sR_1C_1)(1+sR_2C_2)}.

  • 2076 Chaitra · 10 marks

The given mechanical system has force f(t) as input and x1 and x2 as displacement outputs. Draw equivalent f-v and f-I analogous electrical circuit and determine the transfer function with x2 as output. [Figure: fixed wall – spring k – mass M2 (displacement x2) – damper B1 – mass M1 (displacement x1, force F(t) applied); M1 is also connected to the fixed wall through damper B2]

Answer

Displacements: x1x_1 for M1M_1 (where FF acts) and x2x_2 for M2M_2. Assume no ground friction other than the dampers shown.

Differential equations

  • M1M_1: inertia, damper B1B_1 (relative to M2M_2), damper B2B_2 to the wall.
  • M2M_2: inertia, damper B1B_1, spring kk to the wall.
M1x¨1+B1(x˙1−x˙2)+B2x˙1=F(t)M2x¨2+B1(x˙2−x˙1)+kx2=0\begin{aligned} M_1\ddot x_1 + B_1(\dot x_1 - \dot x_2) + B_2\dot x_1 &= F(t) \\ M_2\ddot x_2 + B_1(\dot x_2 - \dot x_1) + kx_2 &= 0 \end{aligned}

F-V analogous circuit (mesh)

F→eF\to e, M→LM\to L, B→RB\to R, k→1/Ck\to 1/C, x˙→i\dot x\to i, x→qx \to q.

L1di1dt+R1(i1−i2)+R2i1=eL2di2dt+R1(i2−i1)+1C∫i2 dt=0\begin{aligned} L_1\frac{di_1}{dt} + R_1(i_1 - i_2) + R_2 i_1 &= e \\ L_2\frac{di_2}{dt} + R_1(i_2 - i_1) + \frac{1}{C}\int i_2\,dt &= 0 \end{aligned}

(L1=M1L_1 = M_1, L2=M2L_2 = M_2, R1=B1R_1 = B_1, R2=B2R_2 = B_2, C=1/kC = 1/k)

     L1=M1   R2=B2         L2=M2
  +--LLLL--/\/\/--+--------LLLL--+
  |               |              |
 (~) e         R1 = B1        C = 1/k
  |      i1       |      i2      |
  +---------------+--------------+

F-I analogous circuit (node)

F→iF\to i, M→CM\to C, B→1/RB\to 1/R, k→1/Lk\to 1/L, x˙→v\dot x\to v.

C1dv1dt+v1−v2R1+v1R2=iC2dv2dt+v2−v1R1+1L∫v2 dt=0\begin{aligned} C_1\frac{dv_1}{dt} + \frac{v_1 - v_2}{R_1} + \frac{v_1}{R_2} &= i \\ C_2\frac{dv_2}{dt} + \frac{v_2 - v_1}{R_1} + \frac{1}{L}\int v_2\,dt &= 0 \end{aligned}

(C1=M1C_1 = M_1, C2=M2C_2 = M_2, R1=1/B1R_1 = 1/B_1, R2=1/B2R_2 = 1/B_2, L=1/kL = 1/k)

   v1          R1 = 1/B1        v2
  +---+----+----/\/\/----+----+----+
  |   |    |             |    |    |
 (i) C1   R2            C2    L    |
  |  =M1 =1/B2          =M2  =1/k  |
  +---+----+-------------+----+----+ gnd

Transfer function X2(s)/F(s)X_2(s)/F(s)

Laplace form:

[M1s2+(B1+B2)s]X1−B1sX2=F−B1sX1+(M2s2+B1s+k)X2=0\begin{aligned} \left[M_1s^2 + (B_1 + B_2)s\right]X_1 - B_1sX_2 &= F \\ -B_1sX_1 + (M_2s^2 + B_1s + k)X_2 &= 0 \end{aligned}

By Cramer's rule:

X2F=B1s[M1s2+(B1+B2)s](M2s2+B1s+k)−B12s2\frac{X_2}{F} = \frac{B_1s}{\left[M_1s^2 + (B_1+B_2)s\right](M_2s^2 + B_1s + k) - B_1^2s^2}

Cancelling the common factor ss:

X2(s)F(s)=B1M1M2s3+(B1M1+B1M2+B2M2)s2+(B1B2+kM1)s+(B1+B2)k\frac{X_2(s)}{F(s)} = \frac{B_1}{M_1M_2s^3 + (B_1M_1 + B_1M_2 + B_2M_2)s^2 + (B_1B_2 + kM_1)s + (B_1 + B_2)k}
  • 2075 Chaitra · 8 marks

Find the transfer function θ(s)/T(s) for the mechanical rotational system shown below. Also develop T-I analogous circuit. [Figure: torque T(t) applied to inertia J1; J1 coupled through torsional spring k1 in parallel with damper b1 to a shaft point, which has viscous damper b2 to the fixed frame; that point is coupled through torsional spring k2 to inertia J2, whose angle is θ(t)]

Answer

Let θ1\theta_1 = angle of J1J_1, θ3\theta_3 = angle of the massless shaft point (where b2b_2 is attached) and θ\theta = angle of J2J_2. Bearings of J1J_1, J2J_2 are taken as frictionless (only the shown elements act).

Differential equations (Laplace form)

J1s2θ1+(b1s+k1)(θ1−θ3)=T(b1s+k1)(θ3−θ1)+b2s θ3+k2(θ3−θ)=0J2s2θ+k2(θ−θ3)=0\begin{aligned} J_1s^2\theta_1 + (b_1s + k_1)(\theta_1 - \theta_3) &= T \\ (b_1s + k_1)(\theta_3 - \theta_1) + b_2s\,\theta_3 + k_2(\theta_3 - \theta) &= 0 \\ J_2s^2\theta + k_2(\theta - \theta_3) &= 0 \end{aligned}

In matrix form:

[J1s2+b1s+k1−(b1s+k1)0−(b1s+k1)(b1+b2)s+k1+k2−k20−k2J2s2+k2][θ1θ3θ]=[T00]\begin{bmatrix} J_1s^2 + b_1s + k_1 & -(b_1s+k_1) & 0 \\ -(b_1s+k_1) & (b_1+b_2)s + k_1 + k_2 & -k_2 \\ 0 & -k_2 & J_2s^2 + k_2 \end{bmatrix} \begin{bmatrix}\theta_1\\ \theta_3\\ \theta\end{bmatrix} = \begin{bmatrix}T\\0\\0\end{bmatrix}

Transfer function

By Cramer's rule, the numerator for θ\theta is the product of the off-diagonal couplings, (b1s+k1)k2(b_1s+k_1)k_2:

θ(s)T(s)=k2(b1s+k1)Δ\frac{\theta(s)}{T(s)} = \frac{k_2(b_1s + k_1)}{\Delta}

where

Δ=(J1s2+b1s+k1){[(b1+b2)s+k1+k2](J2s2+k2)−k22}−(b1s+k1)2(J2s2+k2)\Delta = (J_1s^2 + b_1s + k_1)\left\{\left[(b_1+b_2)s + k_1 + k_2\right](J_2s^2 + k_2) - k_2^2\right\} - (b_1s + k_1)^2(J_2s^2 + k_2)

Expanded:

Δ=  J1J2(b1+b2)s5+[J1J2(k1+k2)+J2b1b2]s4+[(J1+J2)b1k2+J1b2k2+J2b2k1]s3+[(J1+J2)k1k2+b1b2k2]s2+b2k1k2 s\begin{aligned} \Delta =\; & J_1J_2(b_1 + b_2)s^5 + \left[J_1J_2(k_1 + k_2) + J_2b_1b_2\right]s^4 \\ & + \left[(J_1 + J_2)b_1k_2 + J_1b_2k_2 + J_2b_2k_1\right]s^3 \\ & + \left[(J_1 + J_2)k_1k_2 + b_1b_2k_2\right]s^2 + b_2k_1k_2\,s \end{aligned}

The factor ss in Δ\Delta shows a free integration: nothing ties the system to the frame by a spring.

T-I analogy

Torque →\to current, J→CJ \to C, b→1/Rb \to 1/R, k→1/Lk \to 1/L, angular velocity ω→\omega \to node voltage.

MechanicalElectrical
TTcurrent source ii
J1J_1, J2J_2C1=J1C_1 = J_1, C2=J2C_2 = J_2
b1b_1, b2b_2R1=1/b1R_1 = 1/b_1, R2=1/b2R_2 = 1/b_2
k1k_1, k2k_2L1=1/k1L_1 = 1/k_1, L2=1/k2L_2 = 1/k_2
ω1,ω3,ω\omega_1, \omega_3, \omegav1,v3,v2v_1, v_3, v_2
  v1       L1 ∥ R1       v3     L2=1/k2   v2
 +--+----+--LLLL----+--+-----+---LLLL---+---+
 |  |    |          |  |     |          |   |
(i) C1   +--/\/\/---+  |    R2         C2   |
 | =J1                 |   =1/b2      =J2   |
 +--+------------------+-----+----------+---+ gnd

(L1L_1 and R1R_1 are in parallel between nodes 1 and 3.)

  • 2075 Asoj · 6+2 marks

For the mechanical system shown below find the transfer function X2(s)/F(s). Draw the force voltage analogy. [Figure: wall – spring k1 – mass M2 (displacement x2, friction B2 with ground); M2 connected to mass M1 (displacement x1) through damper B1 in parallel with spring k2; force F applied to M1; M1 has friction B3 with ground]

Answer

Displacements: x1x_1 for M1M_1 (force FF applied, ground friction B3B_3) and x2x_2 for M2M_2 (ground friction B2B_2, spring k1k_1 to wall). B1∥k2B_1 \parallel k_2 couples the masses.

Equations of motion

M1x¨1+B3x˙1+B1(x˙1−x˙2)+k2(x1−x2)=FM2x¨2+B2x˙2+k1x2+B1(x˙2−x˙1)+k2(x2−x1)=0\begin{aligned} M_1\ddot x_1 + B_3\dot x_1 + B_1(\dot x_1 - \dot x_2) + k_2(x_1 - x_2) &= F \\ M_2\ddot x_2 + B_2\dot x_2 + k_1x_2 + B_1(\dot x_2 - \dot x_1) + k_2(x_2 - x_1) &= 0 \end{aligned}

Laplace form:

[M1s2+(B1+B3)s+k2]X1−(B1s+k2)X2=F−(B1s+k2)X1+[M2s2+(B1+B2)s+k1+k2]X2=0\begin{aligned} \left[M_1s^2 + (B_1 + B_3)s + k_2\right]X_1 - (B_1s + k_2)X_2 &= F \\ -(B_1s + k_2)X_1 + \left[M_2s^2 + (B_1 + B_2)s + k_1 + k_2\right]X_2 &= 0 \end{aligned}

Transfer function

X2(s)F(s)=B1s+k2Δ\frac{X_2(s)}{F(s)} = \frac{B_1s + k_2}{\Delta} Δ=[M1s2+(B1+B3)s+k2][M2s2+(B1+B2)s+k1+k2]−(B1s+k2)2\Delta = \left[M_1s^2 + (B_1+B_3)s + k_2\right]\left[M_2s^2 + (B_1+B_2)s + k_1 + k_2\right] - (B_1s + k_2)^2

Expanded:

Δ=  M1M2s4+[B1(M1+M2)+B2M1+B3M2]s3+[B1B2+B1B3+B2B3+M1(k1+k2)+M2k2]s2+[B1k1+B2k2+B3(k1+k2)]s+k1k2\begin{aligned} \Delta =\; & M_1M_2s^4 + \left[B_1(M_1 + M_2) + B_2M_1 + B_3M_2\right]s^3 \\ & + \left[B_1B_2 + B_1B_3 + B_2B_3 + M_1(k_1 + k_2) + M_2k_2\right]s^2 \\ & + \left[B_1k_1 + B_2k_2 + B_3(k_1 + k_2)\right]s + k_1k_2 \end{aligned}

F-V analogy

F→eF\to e, M→LM\to L, B→RB\to R, k→1/Ck\to 1/C, x˙→i\dot x \to i. Mesh 1 (i1i_1) for M1M_1, mesh 2 (i2i_2) for M2M_2; the coupling B1∥k2B_1 \parallel k_2 becomes the common branch R1R_1 in series with C2C_2 (parallel mechanical elements carry the same velocity difference, so they become series elements in the shared branch).

    L1=M1  R3=B3            L2=M2   R2=B2
 +--LLLL--/\/\/--+----------LLLL---/\/\/--+
 |               |                        |
(~) e         R1 = B1                  C1 = 1/k1
 |               |                        |
 |    i1      C2 = 1/k2        i2         |
 +---------------+------------------------+
L1di1dt+R3i1+R1(i1−i2)+1C2∫(i1−i2)dt=eL2di2dt+R2i2+1C1∫i2dt+R1(i2−i1)+1C2∫(i2−i1)dt=0\begin{aligned} L_1\frac{di_1}{dt} + R_3i_1 + R_1(i_1 - i_2) + \frac{1}{C_2}\int(i_1 - i_2)dt &= e \\ L_2\frac{di_2}{dt} + R_2i_2 + \frac{1}{C_1}\int i_2dt + R_1(i_2 - i_1) + \frac{1}{C_2}\int(i_2 - i_1)dt &= 0 \end{aligned}
  • 2075 Asoj · 4 marks

Develop block diagram model for the circuit shown in figure below. [Figure: input Vi; series R1 and L1 to a node; capacitor C from that node to ground; series L2 to the output node; R2 from output node to ground; Vo across R2]

Answer

Let V1V_1 = voltage across CC, I1I_1 = current in R1R_1–L1L_1, I2I_2 = current in L2L_2–R2R_2.

Element equations (Laplace)

I1(s)=Vi(s)−V1(s)R1+sL1V1(s)=1sC[I1(s)−I2(s)]I2(s)=V1(s)R2+sL2Vo(s)=R2 I2(s)\begin{aligned} I_1(s) &= \frac{V_i(s) - V_1(s)}{R_1 + sL_1} \\ V_1(s) &= \frac{1}{sC}\left[I_1(s) - I_2(s)\right] \\ I_2(s) &= \frac{V_1(s)}{R_2 + sL_2} \\ V_o(s) &= R_2\,I_2(s) \end{aligned}

Block diagram

Here Z1=R1+sL1Z_1 = R_1 + sL_1 and Z2=R2+sL2Z_2 = R_2 + sL_2.

 Vi─►(Σ)►[1/Z1]─►(Σ)►[1/sC]─┬─►[1/Z2]─┬─►[R2]─►Vo
      ▲-          ▲-        │          │
      │           └── I2 ───┼──────────┘
      └── V1 ───────────────┘

Reduction

Inner loop (between I1I_1 and V1V_1, feedback via 1R2+sL2\frac{1}{R_2+sL_2}):

V1I1=1sC1+1sC(R2+sL2)=R2+sL2sC(R2+sL2)+1\frac{V_1}{I_1} = \frac{\frac{1}{sC}}{1 + \frac{1}{sC(R_2+sL_2)}} = \frac{R_2 + sL_2}{sC(R_2+sL_2) + 1}

Outer loop with forward 1R1+sL1\frac{1}{R_1+sL_1} and unity feedback of V1V_1, then multiply by R2R2+sL2\frac{R_2}{R_2+sL_2}:

Vo(s)Vi(s)=R2(R1+sL1)+(R2+sL2)+sC(R1+sL1)(R2+sL2)\frac{V_o(s)}{V_i(s)} = \frac{R_2}{(R_1 + sL_1) + (R_2 + sL_2) + sC(R_1 + sL_1)(R_2 + sL_2)}

Expanded:

Vo(s)Vi(s)=R2CL1L2s3+C(L1R2+L2R1)s2+(CR1R2+L1+L2)s+R1+R2\frac{V_o(s)}{V_i(s)} = \frac{R_2}{CL_1L_2s^3 + C(L_1R_2 + L_2R_1)s^2 + (CR_1R_2 + L_1 + L_2)s + R_1 + R_2}
  • 2074 Chaitra · 8 marks

Find the transfer function θm(s)/Va(s) of the system below by constructing the block diagram. [Figure: armature-controlled DC motor; armature voltage Va drives armature current Ia through resistance ra and inductance La; field current If = constant; back emf Eb; motor torque constant Km; motor shaft (angle θm, torque τ) drives a load of inertia J with viscous friction f to the fixed frame]

Answer

With IfI_f constant the motor is armature controlled; flux is constant, so torque ∝Ia\propto I_a and back emf ∝\propto speed.

Equations

Va(s)=(ra+sLa)Ia(s)+Eb(s)(armature KVL)τ(s)=KmIa(s)(torque)τ(s)=(Js2+fs) θm(s)(load)Eb(s)=Kb s θm(s)(back emf)\begin{aligned} V_a(s) &= (r_a + sL_a)I_a(s) + E_b(s) && \text{(armature KVL)} \\ \tau(s) &= K_m I_a(s) && \text{(torque)} \\ \tau(s) &= (Js^2 + fs)\,\theta_m(s) && \text{(load)} \\ E_b(s) &= K_b\,s\,\theta_m(s) && \text{(back emf)} \end{aligned}

(KbK_b = back-emf constant; in SI units Kb=KmK_b = K_m.)

Constructing the block diagram

Each equation is written as output = block × input:

  • Ia=1ra+sLa(Va−Eb)I_a = \dfrac{1}{r_a + sL_a}(V_a - E_b)
  • τ=KmIa\tau = K_m I_a
  • ωm=1Js+fτ\omega_m = \dfrac{1}{Js + f}\tau, and θm=1sωm\theta_m = \dfrac{1}{s}\omega_m
  • Eb=Kb ωmE_b = K_b\,\omega_m (feedback)
 Va  +     ┌────────┐ Ia ┌────┐ τ ┌────────┐ ωm  ┌───┐ θm
 ───►(Σ)──►│   1    ├───►│ Km ├──►│   1    ├──┬─►│1/s├──►
      ▲ -  │ ra+sLa │    └────┘   │ Js + f │  │  └───┘
      │    └────────┘             └────────┘  │
      │ Eb           ┌────┐                   │
      └──────────────┤ Kb │◄──────────────────┘
                     └────┘

Reduction

  1. Combine the three forward blocks in cascade:
G(s)=Km(ra+sLa)(Js+f)G(s) = \frac{K_m}{(r_a + sL_a)(Js + f)}
  1. Eliminate the back-emf feedback loop (H=KbH = K_b):
ωm(s)Va(s)=G1+GKb=Km(ra+sLa)(Js+f)+KmKb\frac{\omega_m(s)}{V_a(s)} = \frac{G}{1 + GK_b} = \frac{K_m}{(r_a + sL_a)(Js + f) + K_mK_b}
  1. Multiply by the integrator 1/s1/s:
θm(s)Va(s)=Kms[(ra+sLa)(Js+f)+KmKb]\frac{\theta_m(s)}{V_a(s)} = \frac{K_m}{s\left[(r_a + sL_a)(Js + f) + K_mK_b\right]}

Expanded:

θm(s)Va(s)=Kms[LaJs2+(raJ+Laf)s+raf+KmKb]\frac{\theta_m(s)}{V_a(s)} = \frac{K_m}{s\left[L_aJs^2 + (r_aJ + L_af)s + r_af + K_mK_b\right]}

Simplified form

Neglecting LaL_a (electrical time constant much smaller than mechanical):

θm(s)Va(s)=Ks(τms+1),K=Kmraf+KmKb,τm=raJraf+KmKb\frac{\theta_m(s)}{V_a(s)} = \frac{K}{s(\tau_m s + 1)}, \quad K = \frac{K_m}{r_af + K_mK_b}, \quad \tau_m = \frac{r_aJ}{r_af + K_mK_b}

The back emf acts like extra viscous friction KmKb/raK_mK_b/r_a, which is why the effective time constant is reduced.

  • 2074 Asoj · 6 marks

Find the transfer function X2(s)/F(s) for the mechanical system of figure below. Also draw the F-V and F-I analogy circuit of the system. [Figure: wall – spring K1 – mass M1 (displacement X1); M1 connected to mass M2 (displacement X2) through damper D1 in parallel with spring K2; M2 connected to mass M3 (displacement X3) through spring K3; force F(t) applied to M3; ground friction Fv1 = Fv2 = Fv3 = 0]

Answer

Displacements X1,X2,X3X_1, X_2, X_3 of M1,M2,M3M_1, M_2, M_3; no ground friction (fv1=fv2=fv3=0f_{v1} = f_{v2} = f_{v3} = 0).

Equations of motion (Laplace)

(M1s2+D1s+K1+K2)X1−(D1s+K2)X2=0−(D1s+K2)X1+(M2s2+D1s+K2+K3)X2−K3X3=0−K3X2+(M3s2+K3)X3=F\begin{aligned} (M_1s^2 + D_1s + K_1 + K_2)X_1 - (D_1s + K_2)X_2 &= 0 \\ -(D_1s + K_2)X_1 + (M_2s^2 + D_1s + K_2 + K_3)X_2 - K_3X_3 &= 0 \\ -K_3X_2 + (M_3s^2 + K_3)X_3 &= F \end{aligned}

Transfer function by Cramer's rule

X2(s)F(s)=K3(M1s2+D1s+K1+K2)Δ\frac{X_2(s)}{F(s)} = \frac{K_3(M_1s^2 + D_1s + K_1 + K_2)}{\Delta} Δ=  (M1s2+D1s+K1+K2)[(M2s2+D1s+K2+K3)(M3s2+K3)−K32]−(D1s+K2)2(M3s2+K3)\begin{aligned} \Delta =\; & (M_1s^2 + D_1s + K_1 + K_2)\left[(M_2s^2 + D_1s + K_2 + K_3)(M_3s^2 + K_3) - K_3^2\right] \\ & - (D_1s + K_2)^2(M_3s^2 + K_3) \end{aligned}

(Δ\Delta is a 6th-order polynomial; its constant term is K1K2K3K_1K_2K_3.)

F-V analogy (mesh)

F→eF\to e, M→LM\to L, D→RD\to R, K→1/CK\to 1/C, velocity →\to mesh current. Elements between two masses are shared by their meshes.

  L1=M1       L2=M2          L3=M3
 +-LLLL-+-----LLLL----+------LLLL----+
 |      |             |              |
C1     R1=D1        C3=1/K3         (~) e
=1/K1   |             |              |
 |     C2=1/K2        |              |
 | i1   |      i2     |      i3      |
 +------+-------------+--------------+

Mesh 1: L1,C1L_1, C_1 plus shared (R1+C2R_1 + C_2); mesh 2: L2L_2 plus shared branches; mesh 3: L3L_3, source ee, shared C3C_3.

F-I analogy (node)

F→iF\to i, M→CM\to C, D→1/RD\to 1/R, K→1/LK\to 1/L, velocity →\to node voltage.

  v1    L2∥R1      v2    L3=1/K3   v3
 +--+--+-LLLL-+--+--+----LLLL---+--+--+
 |  |  |      |  |  |           |  |  |
C1  L1 +-/\/\-+  | C2          C3  | (i)
=M1 =1/K1  R1    | =M2         =M3 |  |
 |  |      =1/D1 |  |           |  |  |
 +--+------------+--+-----------+--+--+ gnd

L2=1/K2L_2 = 1/K_2 and R1=1/D1R_1 = 1/D_1 are in parallel between nodes 1 and 2.

  • 2074 Asoj · 4 marks

Find transfer function of an op-amp model as below. [Figure: inverting amplifier; Vin through R1 to the inverting (−) input; feedback from output to the inverting input through Rs in parallel with capacitor C; non-inverting (+) input connected to ground through R2; output Vout]

Answer

This is an inverting amplifier with a feedback impedance Zf=Rs∥CZ_f = R_s \parallel C; R2R_2 on the non-inverting input only compensates bias current.

Assumptions (ideal op-amp)

  • No current into the input terminals, so no drop across R2R_2 and V+=0V_+ = 0.
  • Virtual ground: V−=V+=0V_- = V_+ = 0.

Impedances

Z1=R1,Zf=Rs∥1sC=Rs⋅1sCRs+1sC=Rs1+sRsCZ_1 = R_1, \qquad Z_f = R_s \parallel \frac{1}{sC} = \frac{R_s\cdot\frac{1}{sC}}{R_s + \frac{1}{sC}} = \frac{R_s}{1 + sR_sC}

Derivation

Current through R1R_1 equals current through ZfZ_f:

Vin(s)−0R1=0−Vout(s)Zf\frac{V_{in}(s) - 0}{R_1} = \frac{0 - V_{out}(s)}{Z_f} Vout(s)Vin(s)=−ZfZ1\frac{V_{out}(s)}{V_{in}(s)} = -\frac{Z_f}{Z_1}

Transfer function

Vout(s)Vin(s)=−RsR1⋅11+sRsC\frac{V_{out}(s)}{V_{in}(s)} = -\frac{R_s}{R_1}\cdot\frac{1}{1 + sR_sC}

This is a first-order lag (inverting low-pass) circuit:

  • DC gain =−Rs/R1= -R_s/R_1
  • pole at s=−1RsCs = -\dfrac{1}{R_sC}, time constant τ=RsC\tau = R_sC

For example, with R1=10 kΩR_1 = 10\ \text{k}\Omega, Rs=100 kΩR_s = 100\ \text{k}\Omega, C=1 μFC = 1\ \mu\text{F}, the gain is −10/(1+0.1s)-10/(1 + 0.1s). If Rs→∞R_s \to \infty the circuit becomes a pure integrator −1/(sR1C)-1/(sR_1C).

  • 2073 Shrawan · 8 marks

Find transfer function (consider displacement of mass M2 as output) for the given mechanical system. Also develop force-current analogous circuit. [Figure: wall – spring k1 – mass M1 (force F(t) applied, displacement x1, friction B1 with ground) – damper B2 in series with spring k2 – mass M2 (displacement x2, friction B3 with ground) – spring k3 – wall]

Answer

Displacements: x1x_1 (M1M_1), x2x_2 (M2M_2) and x3x_3 = junction between damper B2B_2 (on the M1M_1 side) and spring k2k_2 (on the M2M_2 side).

Differential equations

M1x¨1+B1x˙1+k1x1+B2(x˙1−x˙3)=F(t)B2(x˙3−x˙1)+k2(x3−x2)=0M2x¨2+B3x˙2+k3x2+k2(x2−x3)=0\begin{aligned} M_1\ddot x_1 + B_1\dot x_1 + k_1x_1 + B_2(\dot x_1 - \dot x_3) &= F(t) \\ B_2(\dot x_3 - \dot x_1) + k_2(x_3 - x_2) &= 0 \\ M_2\ddot x_2 + B_3\dot x_2 + k_3x_2 + k_2(x_2 - x_3) &= 0 \end{aligned}

Eliminating the junction

The series B2B_2–k2k_2 combination acts like one element:

Keq(s)=B2s⋅k2B2s+k2K_{eq}(s) = \frac{B_2s\cdot k_2}{B_2s + k_2}

With a=M1s2+B1s+k1a = M_1s^2 + B_1s + k_1 and b=M2s2+B3s+k3b = M_2s^2 + B_3s + k_3:

(a+Keq)X1−KeqX2=F−KeqX1+(b+Keq)X2=0\begin{aligned} (a + K_{eq})X_1 - K_{eq}X_2 &= F \\ -K_{eq}X_1 + (b + K_{eq})X_2 &= 0 \end{aligned}

Transfer function

X2F=Keqab+Keq(a+b)\frac{X_2}{F} = \frac{K_{eq}}{ab + K_{eq}(a + b)}

Multiplying numerator and denominator by (B2s+k2)(B_2s + k_2):

X2(s)F(s)=B2k2 s(B2s+k2)(M1s2+B1s+k1)(M2s2+B3s+k3)+B2k2s[(M1+M2)s2+(B1+B3)s+k1+k3]\frac{X_2(s)}{F(s)} = \frac{B_2k_2\,s}{(B_2s + k_2)(M_1s^2 + B_1s + k_1)(M_2s^2 + B_3s + k_3) + B_2k_2s\left[(M_1+M_2)s^2 + (B_1+B_3)s + k_1 + k_3\right]}

The denominator is 5th order (leading term B2M1M2s5B_2M_1M_2s^5, constant term k1k2k3k_1k_2k_3).

Force-current analogy

F→iF\to i, M→CM\to C, B→1/RB\to 1/R, k→1/Lk\to 1/L, velocity →\to node voltage.

ElementAnalogConnected between
M1M_1, B1B_1, k1k_1C1C_1, 1/B11/B_1, 1/k11/k_1node 1 and ground
B2B_2R=1/B2R = 1/B_2node 1 and node 3
k2k_2L=1/k2L = 1/k_2node 3 and node 2
M2M_2, B3B_3, k3k_3C2C_2, 1/B31/B_3, 1/k31/k_3node 2 and ground
  v1                v3                v2
 +--+--+--+--/\/\/--+--LLLL--+--+--+--+
 |  |  |  |  1/B2      1/k2  |  |  |
(i) C1 R  L                  C2 R  L
 | =M1 1/ 1/                =M2 1/ 1/
 |  |  B1 k1                 |  B3 k3
 +--+--+--+------------------+--+--+ gnd

Node equation at node 1: C1dv1dt+B1v1+k1∫v1dt+B2(v1−v3)=iC_1\dfrac{dv_1}{dt} + B_1v_1 + k_1\displaystyle\int v_1dt + B_2(v_1 - v_3) = i (conductances written directly), and similarly for nodes 3 and 2.

  • 2072 Chaitra · 8 marks

The given mechanical system has force f(t) as input and x1 and x2 as displacement outputs. Draw equivalent F-V analogous circuit and determine the transfer functions X1(s)/F(s) and X2(s)/F(s). [Figure: force f(t) applied to mass M1 (displacement x1, friction D1 with ground); M1 connected to mass M2 (displacement x2, friction D2 with ground) through spring k1 in parallel with damper b1; M2 connected to the fixed wall through damper b2 in parallel with spring k2]

Answer

Displacements x1x_1 (M1M_1, ground friction D1D_1, force ff) and x2x_2 (M2M_2, ground friction D2D_2). Coupling: k1∥b1k_1 \parallel b_1; M2M_2 to wall: k2∥b2k_2 \parallel b_2.

Differential equations

M1x¨1+D1x˙1+b1(x˙1−x˙2)+k1(x1−x2)=f(t)M2x¨2+D2x˙2+b2x˙2+k2x2+b1(x˙2−x˙1)+k1(x2−x1)=0\begin{aligned} M_1\ddot x_1 + D_1\dot x_1 + b_1(\dot x_1 - \dot x_2) + k_1(x_1 - x_2) &= f(t) \\ M_2\ddot x_2 + D_2\dot x_2 + b_2\dot x_2 + k_2x_2 + b_1(\dot x_2 - \dot x_1) + k_1(x_2 - x_1) &= 0 \end{aligned}

Laplace form:

[M1s2+(D1+b1)s+k1]X1−(b1s+k1)X2=F−(b1s+k1)X1+[M2s2+(D2+b1+b2)s+k1+k2]X2=0\begin{aligned} \left[M_1s^2 + (D_1 + b_1)s + k_1\right]X_1 - (b_1s + k_1)X_2 &= F \\ -(b_1s + k_1)X_1 + \left[M_2s^2 + (D_2 + b_1 + b_2)s + k_1 + k_2\right]X_2 &= 0 \end{aligned}

Transfer functions (Cramer's rule)

Δ=[M1s2+(D1+b1)s+k1][M2s2+(D2+b1+b2)s+k1+k2]−(b1s+k1)2\Delta = \left[M_1s^2 + (D_1+b_1)s + k_1\right]\left[M_2s^2 + (D_2+b_1+b_2)s + k_1 + k_2\right] - (b_1s + k_1)^2 X1(s)F(s)=M2s2+(D2+b1+b2)s+k1+k2Δ,X2(s)F(s)=b1s+k1Δ\frac{X_1(s)}{F(s)} = \frac{M_2s^2 + (D_2 + b_1 + b_2)s + k_1 + k_2}{\Delta}, \qquad \frac{X_2(s)}{F(s)} = \frac{b_1s + k_1}{\Delta}

Expanded:

Δ=  M1M2s4+[M1(D2+b1+b2)+M2(D1+b1)]s3+[D1D2+(D1+D2)b1+D1b2+b1b2+M1(k1+k2)+M2k1]s2+[(D1+D2)k1+D1k2+b1k2+b2k1]s+k1k2\begin{aligned} \Delta =\; & M_1M_2s^4 + \left[M_1(D_2 + b_1 + b_2) + M_2(D_1 + b_1)\right]s^3 \\ & + \left[D_1D_2 + (D_1 + D_2)b_1 + D_1b_2 + b_1b_2 + M_1(k_1 + k_2) + M_2k_1\right]s^2 \\ & + \left[(D_1 + D_2)k_1 + D_1k_2 + b_1k_2 + b_2k_1\right]s + k_1k_2 \end{aligned}

F-V analogous circuit

f→ef\to e, M→LM\to L, D,b→RD,b\to R, k→1/Ck\to 1/C, x˙→i\dot x\to i.

MechanicalElectrical
M1M_1, M2M_2L1L_1, L2L_2
D1D_1, D2D_2RD1R_{D1}, RD2R_{D2}
b1b_1 (coupling)Rb1R_{b1} in common branch
k1k_1 (coupling)C1=1/k1C_1 = 1/k_1 in common branch
b2b_2, k2k_2Rb2R_{b2}, C2=1/k2C_2 = 1/k_2 in mesh 2
   L1   R_D1             L2   R_D2   R_b2
 +-LLLL-/\/\-+---------+-LLLL-/\/\--/\/\-+
 |           |                           |
(~) e      R_b1                       C2=1/k2
 |           |                           |
 |   i1    C1=1/k1          i2           |
 +-----------+---------------------------+
L1di1dt+RD1i1+Rb1(i1−i2)+1C1∫(i1−i2)dt=eL2di2dt+(RD2+Rb2)i2+1C2∫i2dt+Rb1(i2−i1)+1C1∫(i2−i1)dt=0\begin{aligned} L_1\frac{di_1}{dt} + R_{D1}i_1 + R_{b1}(i_1 - i_2) + \frac{1}{C_1}\int(i_1 - i_2)dt &= e \\ L_2\frac{di_2}{dt} + (R_{D2} + R_{b2})i_2 + \frac{1}{C_2}\int i_2dt + R_{b1}(i_2 - i_1) + \frac{1}{C_1}\int(i_2 - i_1)dt &= 0 \end{aligned}
  • 2071 Chaitra · 6 marks

Develop block diagram model for the system below. [Figure: input Vin; series R1 to node 1; L1 from node 1 to ground; series R2 from node 1 to node 2; L2 from node 2 to ground; output VL across L2]

Answer

Let V1V_1 = voltage of node 1 (across L1L_1), I1I_1 = current in R1R_1, I2I_2 = current in R2R_2 and L2L_2, and VLV_L = output across L2L_2.

Element equations (Laplace)

I1(s)=Vin(s)−V1(s)R1V1(s)=sL1[I1(s)−I2(s)]I2(s)=V1(s)−VL(s)R2VL(s)=sL2 I2(s)\begin{aligned} I_1(s) &= \frac{V_{in}(s) - V_1(s)}{R_1} \\ V_1(s) &= sL_1\left[I_1(s) - I_2(s)\right] \\ I_2(s) &= \frac{V_1(s) - V_L(s)}{R_2} \\ V_L(s) &= sL_2\,I_2(s) \end{aligned}

Block diagram

 Vin─►(Σ)►[1/R1]─►(Σ)►[sL1]─┬─►(Σ)►[1/R2]─┬─►[sL2]─┬─►VL
       ▲-          ▲-       │   ▲-        │        │
       │           └── I2 ──┼───┼─────────┘        │
       └── V1 ──────────────┘   └── VL ────────────┘

Reduction

Each block pair forms a simple loop; let

  • loop La=−sL1R1L_a = -\dfrac{sL_1}{R_1} (via V1V_1),
  • loop Lb=−sL1R2L_b = -\dfrac{sL_1}{R_2} (via I2I_2),
  • loop Lc=−sL2R2L_c = -\dfrac{sL_2}{R_2} (via VLV_L).

Only LaL_a and LcL_c do not touch. Forward path P=1R1⋅sL1⋅1R2⋅sL2P = \dfrac{1}{R_1}\cdot sL_1\cdot\dfrac{1}{R_2}\cdot sL_2, with Δ1=1\Delta_1 = 1. Using Mason's rule (same result as stepwise reduction):

VLVin=s2L1L2R1R21+sL1R1+sL1R2+sL2R2+s2L1L2R1R2\frac{V_L}{V_{in}} = \frac{\frac{s^2L_1L_2}{R_1R_2}}{1 + \frac{sL_1}{R_1} + \frac{sL_1}{R_2} + \frac{sL_2}{R_2} + \frac{s^2L_1L_2}{R_1R_2}}

Multiplying by R1R2R_1R_2:

VL(s)Vin(s)=L1L2s2L1L2s2+(L1R1+L1R2+L2R1)s+R1R2\frac{V_L(s)}{V_{in}(s)} = \frac{L_1L_2s^2}{L_1L_2s^2 + (L_1R_1 + L_1R_2 + L_2R_1)s + R_1R_2}

This is a second-order high-pass network: the output is zero at DC (inductors short) and approaches VinV_{in} at high frequency.

  • 2071 Shrawan · 8 marks

Find the transfer function X2(s)/F(s) for the mechanical system of figure below. Also draw the F-V and F-I analogy circuit of the system. [Figure: wall – spring K1 – mass M1 (force F(t) applied, displacement x1(t), friction D1 with ground); M1 connected to mass M2 (displacement x2(t), friction D2 with ground) through spring K2 in parallel with damper D3; M2 – spring K3 – wall]

Answer

Displacements x1x_1 (M1M_1, force FF, ground friction D1D_1) and x2x_2 (M2M_2, ground friction D2D_2). K1K_1 ties M1M_1 to the left wall, K3K_3 ties M2M_2 to the right wall, and K2∥D3K_2 \parallel D_3 couples them.

Equations of motion (Laplace)

[M1s2+(D1+D3)s+K1+K2]X1−(D3s+K2)X2=F−(D3s+K2)X1+[M2s2+(D2+D3)s+K2+K3]X2=0\begin{aligned} \left[M_1s^2 + (D_1 + D_3)s + K_1 + K_2\right]X_1 - (D_3s + K_2)X_2 &= F \\ -(D_3s + K_2)X_1 + \left[M_2s^2 + (D_2 + D_3)s + K_2 + K_3\right]X_2 &= 0 \end{aligned}

(Sum of impedances at each mass on the diagonal; coupling impedance off the diagonal.)

Transfer function

X2(s)F(s)=D3s+K2Δ\frac{X_2(s)}{F(s)} = \frac{D_3s + K_2}{\Delta} Δ=[M1s2+(D1+D3)s+K1+K2][M2s2+(D2+D3)s+K2+K3]−(D3s+K2)2\Delta = \left[M_1s^2 + (D_1+D_3)s + K_1 + K_2\right]\left[M_2s^2 + (D_2+D_3)s + K_2 + K_3\right] - (D_3s + K_2)^2

Expanded:

Δ=  M1M2s4+[D1M2+D2M1+D3(M1+M2)]s3+[D1D2+D1D3+D2D3+M1(K2+K3)+M2(K1+K2)]s2+[D1(K2+K3)+D2(K1+K2)+D3(K1+K3)]s+K1K2+K1K3+K2K3\begin{aligned} \Delta =\; & M_1M_2s^4 + \left[D_1M_2 + D_2M_1 + D_3(M_1 + M_2)\right]s^3 \\ & + \left[D_1D_2 + D_1D_3 + D_2D_3 + M_1(K_2 + K_3) + M_2(K_1 + K_2)\right]s^2 \\ & + \left[D_1(K_2 + K_3) + D_2(K_1 + K_2) + D_3(K_1 + K_3)\right]s \\ & + K_1K_2 + K_1K_3 + K_2K_3 \end{aligned}

F-V analogy (mesh)

F→eF\to e, M→LM\to L, D→RD\to R, K→1/CK\to 1/C, velocity →\to current.

  L1=M1 R=D1  C=1/K1        L2=M2 R=D2  C=1/K3
 +-LLLL-/\/\--||---+-------+-LLLL-/\/\--||---+
 |                 |                         |
(~) e           R3 = D3                      |
 |                 |                         |
 |    i1        C2 = 1/K2         i2         |
 +-----------------+-------------------------+

F-I analogy (node)

F→iF\to i, M→CM\to C, D→1/RD\to 1/R, K→1/LK\to 1/L, velocity →\to voltage.

  v1        L=1/K2 ∥ R=1/D3        v2
 +--+--+--+---+--LLLL--+---+--+--+--+
 |  |  |  |   +-/\/\/--+   |  |  |  |
(i) C1 R  L               C2  R  L  |
 | =M1 1/ 1/              =M2 1/ 1/ |
 |  |  D1 K1               |  D2 K3 |
 +--+--+--+----------------+--+--+--+ gnd

The mesh (F-V) and node (F-I) equations have exactly the same form as the two mechanical equations above.

  • 2070 Chaitra (old course) · 8 marks

For armature controlled separately excited DC motor, identify the necessary differential equations governing its behaviour and hence derive the dynamic model of such motor.

Answer

In a separately excited DC motor under armature control, the field current ifi_f is held constant, so the air-gap flux ϕ\phi is constant, and the armature voltage vav_a is the control input.

Symbols

RaR_a, LaL_a: armature resistance, inductance; iai_a: armature current; ebe_b: back emf; TmT_m: developed torque; TLT_L: load torque; JJ: total inertia; BB: viscous friction; ω=θ˙\omega = \dot\theta: speed.

Governing differential equations

  1. Armature circuit (KVL):
va=Raia+Ladiadt+ebv_a = R_ai_a + L_a\frac{di_a}{dt} + e_b
  1. Back emf (proportional to flux and speed; flux constant):
eb=Kϕ ω=Kb ωe_b = K\phi\,\omega = K_b\,\omega
  1. Developed torque:
Tm=Kϕ ia=Kt iaT_m = K\phi\,i_a = K_t\,i_a
  1. Mechanical (Newton's law for rotation):
Tm=Jdωdt+Bω+TLT_m = J\frac{d\omega}{dt} + B\omega + T_L
  1. Position: dθdt=ω\dfrac{d\theta}{dt} = \omega

In SI units Kt=Kb=KK_t = K_b = K.

Dynamic model in state-space form

State variables x1=iax_1 = i_a, x2=ωx_2 = \omega, x3=θx_3 = \theta; inputs vav_a and TLT_L:

[i˙aω˙θ˙]=[−RaLa−KbLa0KtJ−BJ0010][iaωθ]+[1La00−1J00][vaTL]\begin{bmatrix}\dot i_a\\ \dot\omega\\ \dot\theta\end{bmatrix} = \begin{bmatrix} -\frac{R_a}{L_a} & -\frac{K_b}{L_a} & 0 \\ \frac{K_t}{J} & -\frac{B}{J} & 0 \\ 0 & 1 & 0 \end{bmatrix} \begin{bmatrix} i_a\\ \omega\\ \theta\end{bmatrix} + \begin{bmatrix} \frac{1}{L_a} & 0 \\ 0 & -\frac{1}{J} \\ 0 & 0\end{bmatrix} \begin{bmatrix} v_a \\ T_L\end{bmatrix}

Output: y=ωy = \omega (speed) or y=θy = \theta (position).

Transfer-function model (with TL=0T_L = 0)

Taking Laplace transforms:

Ia(s)=Va(s)−KbΩ(s)Ra+sLa,Ω(s)=KtIa(s)Js+B\begin{aligned} I_a(s) &= \frac{V_a(s) - K_b\Omega(s)}{R_a + sL_a}, \qquad \Omega(s) = \frac{K_tI_a(s)}{Js + B} \end{aligned} Ω(s)Va(s)=Kt(Ra+sLa)(Js+B)+KtKb,θ(s)Va(s)=Kts[(Ra+sLa)(Js+B)+KtKb]\frac{\Omega(s)}{V_a(s)} = \frac{K_t}{(R_a + sL_a)(Js + B) + K_tK_b}, \qquad \frac{\theta(s)}{V_a(s)} = \frac{K_t}{s\left[(R_a + sL_a)(Js + B) + K_tK_b\right]}

Block diagram of the dynamic model

 Va +    ┌───────┐ Ia ┌────┐ Tm +   ┌──────┐ ω  ┌───┐ θ
 ──►(Σ)─►│  1    ├───►│ Kt ├───►(Σ)►│  1   ├─┬─►│1/s├─►
    ▲ -  │Ra+sLa │    └────┘    ▲-  │Js+B  │ │  └───┘
    │    └───────┘              │   └──────┘ │
    │ Eb         ┌────┐        TL            │
    └────────────┤ Kb ├◄─────────────────────┘
                 └────┘

Remarks

  • Electrical time constant τa=La/Ra\tau_a = L_a/R_a; mechanical time constant τm=J/B\tau_m = J/B. Usually τa≪τm\tau_a \ll \tau_m, so LaL_a may be neglected, giving ΩVa=Km1+sτ\dfrac{\Omega}{V_a} = \dfrac{K_m}{1 + s\tau} with Km=KtRaB+KtKbK_m = \dfrac{K_t}{R_aB + K_tK_b}, τ=RaJRaB+KtKb\tau = \dfrac{R_aJ}{R_aB + K_tK_b}.
  • The back emf provides inherent speed feedback, which makes the motor more stable than a field-controlled motor.
  • 2070 Chaitra (old course) · 10 marks

Draw the free body diagram, write the differential equations and find the mentioned transfer functions X2(s)/F(s) and X1(s)/F(s) of the system below. [Figure: wall – spring k1 – mass M1 (displacement x1, friction B1 with ground) – spring k2 – mass M2 (displacement x2, friction B2 with ground); force F(t) applied to M2]

Answer

Displacements x1x_1 (M1M_1) and x2x_2 (M2M_2), both positive to the right; F(t)F(t) acts on M2M_2.

Free body diagrams

            M1                          M2
 k1x1    ┌──────┐  k2(x1-x2)  k2(x2-x1) ┌──────┐
 ◄───────┤      ├──────►      ◄─────────┤      ├──► F(t)
 M1x1''  │      │                M2x2'' │      │
 ◄───────┤  M1  │                ◄──────┤  M2  │
 B1x1'   │      │                B2x2'  │      │
 ◄───────┤      │                ◄──────┤      │
         └──────┘                       └──────┘
  • On M1M_1: inertia M1x¨1M_1\ddot x_1, friction B1x˙1B_1\dot x_1, spring k1x1k_1x_1 (all opposing motion), spring k2k_2 force k2(x1−x2)k_2(x_1 - x_2) opposing.
  • On M2M_2: applied FF; opposing M2x¨2M_2\ddot x_2, B2x˙2B_2\dot x_2, k2(x2−x1)k_2(x_2 - x_1).

Differential equations

M1x¨1+B1x˙1+k1x1+k2(x1−x2)=0M2x¨2+B2x˙2+k2(x2−x1)=F(t)\begin{aligned} M_1\ddot x_1 + B_1\dot x_1 + k_1x_1 + k_2(x_1 - x_2) &= 0 \\ M_2\ddot x_2 + B_2\dot x_2 + k_2(x_2 - x_1) &= F(t) \end{aligned}

Laplace form

(M1s2+B1s+k1+k2)X1−k2X2=0−k2X1+(M2s2+B2s+k2)X2=F\begin{aligned} (M_1s^2 + B_1s + k_1 + k_2)X_1 - k_2X_2 &= 0 \\ -k_2X_1 + (M_2s^2 + B_2s + k_2)X_2 &= F \end{aligned}

Determinant:

Δ=(M1s2+B1s+k1+k2)(M2s2+B2s+k2)−k22\Delta = (M_1s^2 + B_1s + k_1 + k_2)(M_2s^2 + B_2s + k_2) - k_2^2

Expanded:

Δ=  M1M2s4+(B1M2+B2M1)s3+[B1B2+M1k2+M2(k1+k2)]s2+[B1k2+B2(k1+k2)]s+k1k2\begin{aligned} \Delta =\; & M_1M_2s^4 + (B_1M_2 + B_2M_1)s^3 + \left[B_1B_2 + M_1k_2 + M_2(k_1 + k_2)\right]s^2 \\ & + \left[B_1k_2 + B_2(k_1 + k_2)\right]s + k_1k_2 \end{aligned}

Transfer functions (Cramer's rule)

X2(s)F(s)=M1s2+B1s+k1+k2Δ\frac{X_2(s)}{F(s)} = \frac{M_1s^2 + B_1s + k_1 + k_2}{\Delta} X1(s)F(s)=k2Δ\frac{X_1(s)}{F(s)} = \frac{k_2}{\Delta}

Check: at steady state (s→0s\to 0) with a constant force, X1/F=1/k1X_1/F = 1/k_1 and X2/F=(k1+k2)/(k1k2)=1/k1+1/k2X_2/F = (k_1+k_2)/(k_1k_2) = 1/k_1 + 1/k_2, which is correct for two springs in series.

  • 2070 Chaitra (old course) · 6 marks

Draw the block diagram and reduce it to calculate Vo(s)/Vi(s) for the following network. [Figure: input Vi; series L1 to node 1; R1 from node 1 to ground; series L2 from node 1 to node 2; R2 from node 2 to ground; Vo across R2]

Answer

Let I1I_1 = current in L1L_1, V1V_1 = voltage across R1R_1 (node 1), I2I_2 = current in L2L_2 and R2R_2, VoV_o = voltage across R2R_2.

Element equations (Laplace)

I1=Vi−V1sL1V1=R1(I1−I2)I2=V1−VosL2Vo=R2I2\begin{aligned} I_1 &= \frac{V_i - V_1}{sL_1} \\ V_1 &= R_1(I_1 - I_2) \\ I_2 &= \frac{V_1 - V_o}{sL_2} \\ V_o &= R_2I_2 \end{aligned}

Block diagram

 Vi +   ┌─────┐ I1 +   ┌────┐ V1 +   ┌─────┐ I2  ┌────┐ Vo
 ──►(Σ)►│1/sL1├──►(Σ)─►│ R1 ├─┬─►(Σ)►│1/sL2├──┬─►│ R2 ├─┬─►
    ▲-  └─────┘    ▲-  └────┘ │   ▲-  └─────┘  │  └────┘ │
    │              │          │   │            │         │
    │              └──────────┼───┼────────────┘ I2      │
    └────── V1 ──────────────┘   └──────── Vo ──────────┘

Reduction step by step

  1. Last loop (I2→VoI_2 \to V_o fed back): forward R2sL2\dfrac{R_2}{sL_2}, unity feedback:
VoV1=R2/sL21+R2/sL2=R2sL2+R2,I2V1=1sL2+R2\frac{V_o}{V_1} = \frac{R_2/sL_2}{1 + R_2/sL_2} = \frac{R_2}{sL_2 + R_2}, \qquad \frac{I_2}{V_1} = \frac{1}{sL_2 + R_2}
  1. Middle loop: forward R1R_1 from (I1−I2)(I_1 - I_2) to V1V_1, feedback 1sL2+R2\dfrac{1}{sL_2 + R_2}:
V1I1=R11+R1sL2+R2=R1(sL2+R2)sL2+R1+R2\frac{V_1}{I_1} = \frac{R_1}{1 + \frac{R_1}{sL_2 + R_2}} = \frac{R_1(sL_2 + R_2)}{sL_2 + R_1 + R_2}
  1. First loop: forward 1sL1⋅V1I1\dfrac{1}{sL_1}\cdot\dfrac{V_1}{I_1}, unity feedback:
V1Vi=R1(sL2+R2)sL1(sL2+R1+R2)+R1(sL2+R2)\frac{V_1}{V_i} = \frac{R_1(sL_2 + R_2)}{sL_1(sL_2 + R_1 + R_2) + R_1(sL_2 + R_2)}
  1. Multiply by VoV1=R2sL2+R2\dfrac{V_o}{V_1} = \dfrac{R_2}{sL_2 + R_2}:
Vo(s)Vi(s)=R1R2L1L2s2+(L1R1+L1R2+L2R1)s+R1R2\frac{V_o(s)}{V_i(s)} = \frac{R_1R_2}{L_1L_2s^2 + (L_1R_1 + L_1R_2 + L_2R_1)s + R_1R_2}

At DC the gain is 1 (inductors act as shorts), so this is a second-order low-pass network.

  • 2069 Chaitra · 8+4 marks

Write differential equations governing the mechanical system shown in figure below and find X2(s)/F(s). Also tabulating the necessary analogies draw the Force-Current and Force-Voltage electrical analogous circuit. [Figure: force f(t) applied to mass M1 (displacement x1, friction B1 with ground); M1 connected to mass M2 (displacement x2, on rollers) through damper B12 in parallel with spring K12; M2 connected to the fixed wall through damper B2 in parallel with spring K2]

Answer

Displacements: x1x_1 for M1M_1 (force ff, ground friction B1B_1) and x2x_2 for M2M_2 (on rollers, so no ground friction). B12∥K12B_{12} \parallel K_{12} couples the masses; B2∥K2B_2 \parallel K_2 ties M2M_2 to the wall.

Differential equations

M1x¨1+B1x˙1+B12(x˙1−x˙2)+K12(x1−x2)=f(t)M2x¨2+B12(x˙2−x˙1)+K12(x2−x1)+B2x˙2+K2x2=0\begin{aligned} M_1\ddot x_1 + B_1\dot x_1 + B_{12}(\dot x_1 - \dot x_2) + K_{12}(x_1 - x_2) &= f(t) \\ M_2\ddot x_2 + B_{12}(\dot x_2 - \dot x_1) + K_{12}(x_2 - x_1) + B_2\dot x_2 + K_2x_2 &= 0 \end{aligned}

Transfer function X2(s)/F(s)X_2(s)/F(s)

Laplace form:

[M1s2+(B1+B12)s+K12]X1−(B12s+K12)X2=F−(B12s+K12)X1+[M2s2+(B12+B2)s+K12+K2]X2=0\begin{aligned} \left[M_1s^2 + (B_1 + B_{12})s + K_{12}\right]X_1 - (B_{12}s + K_{12})X_2 &= F \\ -(B_{12}s + K_{12})X_1 + \left[M_2s^2 + (B_{12} + B_2)s + K_{12} + K_2\right]X_2 &= 0 \end{aligned}

From the second equation, X1=M2s2+(B12+B2)s+K12+K2B12s+K12X2X_1 = \dfrac{M_2s^2 + (B_{12}+B_2)s + K_{12} + K_2}{B_{12}s + K_{12}}X_2. Substituting into the first:

X2(s)F(s)=B12s+K12Δ\frac{X_2(s)}{F(s)} = \frac{B_{12}s + K_{12}}{\Delta} Δ=[M1s2+(B1+B12)s+K12][M2s2+(B12+B2)s+K12+K2]−(B12s+K12)2\Delta = \left[M_1s^2 + (B_1+B_{12})s + K_{12}\right]\left[M_2s^2 + (B_{12}+B_2)s + K_{12} + K_2\right] - (B_{12}s + K_{12})^2

Expanded:

Δ=  M1M2s4+[B1M2+B12(M1+M2)+B2M1]s3+[B1B12+B1B2+B12B2+K12(M1+M2)+K2M1]s2+[B1(K12+K2)+B12K2+B2K12]s+K12K2\begin{aligned} \Delta =\; & M_1M_2s^4 + \left[B_1M_2 + B_{12}(M_1 + M_2) + B_2M_1\right]s^3 \\ & + \left[B_1B_{12} + B_1B_2 + B_{12}B_2 + K_{12}(M_1 + M_2) + K_2M_1\right]s^2 \\ & + \left[B_1(K_{12} + K_2) + B_{12}K_2 + B_2K_{12}\right]s + K_{12}K_2 \end{aligned}

Table of analogies

Mechanical (translational)F-V (mesh)F-I (node)
Force ffVoltage eeCurrent ii
Mass MMInductance LLCapacitance CC
Damper BBResistance RRConductance 1/R1/R
Spring KKElastance 1/C1/CReciprocal inductance 1/L1/L
Displacement xxCharge qqFlux linkage ψ\psi
Velocity x˙\dot xCurrent iiVoltage vv
Elements in parallel (same x˙\dot x)SeriesParallel

Force-current analogous circuit

Node 1 ↔ x˙1\dot x_1, node 2 ↔ x˙2\dot x_2; C1=M1C_1 = M_1, C2=M2C_2 = M_2, R1=1/B1R_1 = 1/B_1, R12=1/B12R_{12} = 1/B_{12}, L12=1/K12L_{12} = 1/K_{12}, R2=1/B2R_2 = 1/B_2, L2=1/K2L_2 = 1/K_2.

  v1        L12 ∥ R12          v2
 +--+---+---+--LLLL--+---+---+---+---+
 |  |   |   +-/\/\/--+   |   |   |   |
(i) C1  R1               C2  R2  L2  |
 |  |   |                |   |   |   |
 +--+---+----------------+---+---+---+ gnd
C1dv1dt+v1R1+v1−v2R12+1L12∫(v1−v2)dt=iC2dv2dt+v2−v1R12+1L12∫(v2−v1)dt+v2R2+1L2∫v2dt=0\begin{aligned} C_1\frac{dv_1}{dt} + \frac{v_1}{R_1} + \frac{v_1 - v_2}{R_{12}} + \frac{1}{L_{12}}\int(v_1 - v_2)dt &= i \\ C_2\frac{dv_2}{dt} + \frac{v_2 - v_1}{R_{12}} + \frac{1}{L_{12}}\int(v_2 - v_1)dt + \frac{v_2}{R_2} + \frac{1}{L_2}\int v_2dt &= 0 \end{aligned}

Force-voltage analogous circuit

Mesh 1 ↔ x˙1\dot x_1, mesh 2 ↔ x˙2\dot x_2; L1=M1L_1 = M_1, L2=M2L_2 = M_2, R1=B1R_1 = B_1, R12=B12R_{12} = B_{12}, C12=1/K12C_{12} = 1/K_{12}, R2=B2R_2 = B_2, C2=1/K2C_2 = 1/K_2.

   L1    R1               L2    R2    C2
 +-LLLL-/\/\-+----------+-LLLL-/\/\---||--+
 |           |                            |
(~) e      R12                            |
 |           |                            |
 |   i1     C12            i2             |
 +-----------+----------------------------+
L1di1dt+R1i1+R12(i1−i2)+1C12∫(i1−i2)dt=eL2di2dt+R12(i2−i1)+1C12∫(i2−i1)dt+R2i2+1C2∫i2dt=0\begin{aligned} L_1\frac{di_1}{dt} + R_1i_1 + R_{12}(i_1 - i_2) + \frac{1}{C_{12}}\int(i_1 - i_2)dt &= e \\ L_2\frac{di_2}{dt} + R_{12}(i_2 - i_1) + \frac{1}{C_{12}}\int(i_2 - i_1)dt + R_2i_2 + \frac{1}{C_2}\int i_2dt &= 0 \end{aligned}

Both sets have the same form as the mechanical equations, confirming the analogies.

  • 2068 Chaitra · 10 marks

For an electromechanical system shown below, derive an expression for VL(s) [?] considering it as armature controlled dc motor. Motor: i) Moment of inertia = Jm (ii) Frictional coefficient = Dm (iii) Torsional [rest not printed]. Load: i) Moment of inertia = Jm [?] (ii) Frictional coefficient = Dm [?] (iii) Torsional [rest not printed]. [Figure: armature circuit with supply ea, armature current Ia, resistance Ra, inductance La and back emf eb; field winding Lf, Rf with field current If; motor shaft (angle θm) drives the load (angle θL) through a gear pair N1:N2]

Answer

The printed question is incomplete. It is read here as the standard problem: find the load-angle transfer function θL(s)/Ea(s)\theta_L(s)/E_a(s) of an armature-controlled DC motor driving a load through a gear pair N1:N2N_1{:}N_2. Motor data: inertia JmJ_m, friction DmD_m; load data: inertia JLJ_L, friction DLD_L; shafts are taken as rigid.

Assumptions

  • Field current IfI_f (through LfL_f, RfR_f) is constant, so flux is constant.
  • Torque Tm=KtIaT_m = K_tI_a; back emf eb=Kbθ˙me_b = K_b\dot\theta_m.
  • Gears are ideal (no loss, no backlash): θL=N1N2θm\theta_L = \dfrac{N_1}{N_2}\theta_m.

Step 1: Electrical equation

ea=Raia+Ladiadt+eb  ⇒  Ia(s)=Ea(s)−Kbsθm(s)Ra+sLae_a = R_ai_a + L_a\frac{di_a}{dt} + e_b \;\Rightarrow\; I_a(s) = \frac{E_a(s) - K_bs\theta_m(s)}{R_a + sL_a}

Step 2: Reflect the load to the motor shaft

Load impedances move to the motor side multiplied by (N1/N2)2(N_1/N_2)^2:

Jeq=Jm+JL(N1N2)2,Deq=Dm+DL(N1N2)2J_{eq} = J_m + J_L\left(\frac{N_1}{N_2}\right)^2, \qquad D_{eq} = D_m + D_L\left(\frac{N_1}{N_2}\right)^2

Step 3: Mechanical equation at the motor shaft

Tm=KtIa=Jeqθ¨m+Deqθ˙m  ⇒  KtIa(s)=s(Jeqs+Deq)θm(s)T_m = K_tI_a = J_{eq}\ddot\theta_m + D_{eq}\dot\theta_m \;\Rightarrow\; K_tI_a(s) = s(J_{eq}s + D_{eq})\theta_m(s)

Step 4: Motor-angle transfer function

Substituting IaI_a:

Kt[Ea−Kbsθm]=s(Ra+sLa)(Jeqs+Deq)θmK_t\left[E_a - K_bs\theta_m\right] = s(R_a + sL_a)(J_{eq}s + D_{eq})\theta_m θm(s)Ea(s)=Kts[(Ra+sLa)(Jeqs+Deq)+KtKb]\frac{\theta_m(s)}{E_a(s)} = \frac{K_t}{s\left[(R_a + sL_a)(J_{eq}s + D_{eq}) + K_tK_b\right]}

Step 5: Load-angle transfer function

θL(s)Ea(s)=N1N2⋅Kts[(Ra+sLa)(Jeqs+Deq)+KtKb]\frac{\theta_L(s)}{E_a(s)} = \frac{N_1}{N_2}\cdot\frac{K_t}{s\left[(R_a + sL_a)(J_{eq}s + D_{eq}) + K_tK_b\right]}

Expanded denominator: s[LaJeqs2+(RaJeq+LaDeq)s+RaDeq+KtKb]s\left[L_aJ_{eq}s^2 + (R_aJ_{eq} + L_aD_{eq})s + R_aD_{eq} + K_tK_b\right].

The load speed is ΩL(s)=sθL(s)\Omega_L(s) = s\theta_L(s), so ΩL(s)Ea(s)=(N1/N2)Kt(Ra+sLa)(Jeqs+Deq)+KtKb\dfrac{\Omega_L(s)}{E_a(s)} = \dfrac{(N_1/N_2)K_t}{(R_a + sL_a)(J_{eq}s + D_{eq}) + K_tK_b}.

Block diagram

 Ea +   ┌──────┐Ia┌──┐Tm┌───────┐ωm┌───┐θm┌─────┐ θL
 ──►(Σ)►│  1   ├─►│Kt├─►│   1   ├┬►│1/s├─►│N1/N2├──►
    ▲-  │Ra+sLa│  └──┘  │Jeq s+ ││ └───┘  └─────┘
    │   └──────┘        │  Deq  ││
    │ Eb      ┌────┐    └───────┘│
    └─────────┤ Kb │◄────────────┘
              └────┘

Simplified form (La≈0L_a \approx 0)

θL(s)Ea(s)=Ks(τs+1),K=(N1/N2)KtRaDeq+KtKb,τ=RaJeqRaDeq+KtKb\frac{\theta_L(s)}{E_a(s)} = \frac{K}{s(\tau s + 1)}, \quad K = \frac{(N_1/N_2)K_t}{R_aD_{eq} + K_tK_b}, \quad \tau = \frac{R_aJ_{eq}}{R_aD_{eq} + K_tK_b}

A large reduction ratio (N2≫N1N_2 \gg N_1) makes the load inertia and friction look much smaller to the motor, which is why gears are used for matching.

  • 2068 Baisakh (old course) · 6+2 marks

Write differential equation and obtain transfer function of the mechanical system as shown below considering θ3 as output. Also draw torque voltage analogy network. [Figure: torque T applied to inertia J1 (angle θ1, friction B1); J1 – torsional spring K1 – inertia J2 (angle θ2, friction B2); J2 – torsional spring K2 – inertia J3 (angle θ3); J3 – torsional spring K3 – fixed wall]

Answer

Angles θ1,θ2,θ3\theta_1, \theta_2, \theta_3 of J1,J2,J3J_1, J_2, J_3; torque TT on J1J_1. J3J_3 has no friction (none shown); K3K_3 ties J3J_3 to the wall.

Differential equations

J1θ¨1+B1θ˙1+K1(θ1−θ2)=TJ2θ¨2+B2θ˙2+K1(θ2−θ1)+K2(θ2−θ3)=0J3θ¨3+K2(θ3−θ2)+K3θ3=0\begin{aligned} J_1\ddot\theta_1 + B_1\dot\theta_1 + K_1(\theta_1 - \theta_2) &= T \\ J_2\ddot\theta_2 + B_2\dot\theta_2 + K_1(\theta_2 - \theta_1) + K_2(\theta_2 - \theta_3) &= 0 \\ J_3\ddot\theta_3 + K_2(\theta_3 - \theta_2) + K_3\theta_3 &= 0 \end{aligned}

Laplace form (matrix)

[J1s2+B1s+K1−K10−K1J2s2+B2s+K1+K2−K20−K2J3s2+K2+K3][θ1θ2θ3]=[T00]\begin{bmatrix} J_1s^2 + B_1s + K_1 & -K_1 & 0 \\ -K_1 & J_2s^2 + B_2s + K_1 + K_2 & -K_2 \\ 0 & -K_2 & J_3s^2 + K_2 + K_3 \end{bmatrix} \begin{bmatrix}\theta_1\\ \theta_2\\ \theta_3\end{bmatrix} = \begin{bmatrix} T\\0\\0\end{bmatrix}

Transfer function

By Cramer's rule, the numerator is the product of the couplings K1K2K_1K_2:

θ3(s)T(s)=K1K2Δ\frac{\theta_3(s)}{T(s)} = \frac{K_1K_2}{\Delta} Δ=  (J1s2+B1s+K1)[(J2s2+B2s+K1+K2)(J3s2+K2+K3)−K22]−K12(J3s2+K2+K3)\begin{aligned} \Delta =\; & (J_1s^2 + B_1s + K_1)\left[(J_2s^2 + B_2s + K_1 + K_2)(J_3s^2 + K_2 + K_3) - K_2^2\right] \\ & - K_1^2(J_3s^2 + K_2 + K_3) \end{aligned}

Δ\Delta is 6th order, with leading term J1J2J3s6J_1J_2J_3s^6 and constant term K1K2K3K_1K_2K_3, so the steady-state gain for a constant torque is 1/K31/K_3 (all of TT is finally held by K3K_3).

Torque-voltage analogy

T→eT\to e, J→LJ\to L, B→RB\to R, K→1/CK\to 1/C, θ˙→i\dot\theta\to i (mesh currents i1,i2,i3i_1, i_2, i_3).

L1di1dt+R1i1+1C1∫(i1−i2)dt=eL2di2dt+R2i2+1C1∫(i2−i1)dt+1C2∫(i2−i3)dt=0L3di3dt+1C2∫(i3−i2)dt+1C3∫i3dt=0\begin{aligned} L_1\frac{di_1}{dt} + R_1i_1 + \frac{1}{C_1}\int(i_1 - i_2)dt &= e \\ L_2\frac{di_2}{dt} + R_2i_2 + \frac{1}{C_1}\int(i_2 - i_1)dt + \frac{1}{C_2}\int(i_2 - i_3)dt &= 0 \\ L_3\frac{di_3}{dt} + \frac{1}{C_2}\int(i_3 - i_2)dt + \frac{1}{C_3}\int i_3dt &= 0 \end{aligned}
  L1  R1         L2  R2         L3      C3
 +LLLL/\/\-+----LLLL/\/\-+-----LLLL----||--+
 |         |             |                 |
(~) e    C1=1/K1       C2=1/K2             |
 |   i1    |     i2      |       i3        |
 +---------+-------------+-----------------+

(Lk=JkL_k = J_k, Rk=BkR_k = B_k, Ck=1/KkC_k = 1/K_k.)

  • 2067 Asar (old course) · 8 marks

Find the transfer function X1(s)/F(s) of the mechanical system shown in figure 1. Also find the force-voltage analogy of the same system. [Figure 1: wall – spring K1 – mass M1 (displacement x1(t)) – spring K2 – mass M2 (displacement x2(t), force F(t) applied); ground friction B1 = B2 = 0]

Answer

Displacements x1x_1 (M1M_1) and x2x_2 (M2M_2, where FF acts). No friction (B1=B2=0B_1 = B_2 = 0).

Equations of motion

M1x¨1+K1x1+K2(x1−x2)=0M2x¨2+K2(x2−x1)=F(t)\begin{aligned} M_1\ddot x_1 + K_1x_1 + K_2(x_1 - x_2) &= 0 \\ M_2\ddot x_2 + K_2(x_2 - x_1) &= F(t) \end{aligned}

Laplace form

(M1s2+K1+K2)X1−K2X2=0−K2X1+(M2s2+K2)X2=F\begin{aligned} (M_1s^2 + K_1 + K_2)X_1 - K_2X_2 &= 0 \\ -K_2X_1 + (M_2s^2 + K_2)X_2 &= F \end{aligned}

From the first equation, X2=M1s2+K1+K2K2X1X_2 = \dfrac{M_1s^2 + K_1 + K_2}{K_2}X_1. Substituting into the second:

[(M1s2+K1+K2)(M2s2+K2)K2−K2]X1=F\left[\frac{(M_1s^2 + K_1 + K_2)(M_2s^2 + K_2)}{K_2} - K_2\right]X_1 = F

Transfer function

X1(s)F(s)=K2(M1s2+K1+K2)(M2s2+K2)−K22\frac{X_1(s)}{F(s)} = \frac{K_2}{(M_1s^2 + K_1 + K_2)(M_2s^2 + K_2) - K_2^2} X1(s)F(s)=K2M1M2s4+[M1K2+M2(K1+K2)]s2+K1K2\frac{X_1(s)}{F(s)} = \frac{K_2}{M_1M_2s^4 + \left[M_1K_2 + M_2(K_1 + K_2)\right]s^2 + K_1K_2}

Only even powers of ss appear because there is no damping; the poles are on the jωj\omega axis (two undamped natural frequencies).

Force-voltage analogy

F→eF\to e, M→LM\to L, K→1/CK\to 1/C, velocity →\to current.

L1di1dt+1C1∫i1dt+1C2∫(i1−i2)dt=0L2di2dt+1C2∫(i2−i1)dt=e\begin{aligned} L_1\frac{di_1}{dt} + \frac{1}{C_1}\int i_1dt + \frac{1}{C_2}\int(i_1 - i_2)dt &= 0 \\ L_2\frac{di_2}{dt} + \frac{1}{C_2}\int(i_2 - i_1)dt &= e \end{aligned}

with L1=M1L_1 = M_1, L2=M2L_2 = M_2, C1=1/K1C_1 = 1/K_1, C2=1/K2C_2 = 1/K_2.

    L1 = M1                 L2 = M2
 +---LLLL----+-------------LLLL----+
 |           |                     |
C1=1/K1    C2=1/K2                (~) e
 |    i1     |         i2          |
 +-----------+---------------------+
  • 2066 Bhadra (old course) · 8 marks

For the mechanical system shown in fig. 1(b), draw free body diagram, write complete differential equations and identify the transfer function X1(s)/F1(s). [Figure: wall – spring K2 – mass M2 (displacement x2(t), friction B2 with ground); M2 connected to mass M1 (displacement x1(t), friction B1 with ground) through spring K12 in parallel with damper B12; force f1(t) applied to M1]

Answer

Displacements x1x_1 (M1M_1, friction B1B_1, force f1f_1) and x2x_2 (M2M_2, friction B2B_2, spring K2K_2 to the wall). K12∥B12K_{12} \parallel B_{12} couples them.

Free body diagrams

 Opposing forces (to the left)        Applied
 M1x1'', B1x1'
 K12(x1-x2), B12(x1'-x2')   ┌────┐
        ◄───────────────────┤ M1 ├────► f1(t)
                            └────┘
 M2x2'', B2x2', K2x2
 K12(x2-x1), B12(x2'-x1')   ┌────┐
        ◄───────────────────┤ M2 │     (none)
                            └────┘

(For each mass, the applied force acts in the positive direction and all reaction forces oppose it.)

Differential equations

M1x¨1+B1x˙1+B12(x˙1−x˙2)+K12(x1−x2)=f1(t)M2x¨2+B2x˙2+K2x2+B12(x˙2−x˙1)+K12(x2−x1)=0\begin{aligned} M_1\ddot x_1 + B_1\dot x_1 + B_{12}(\dot x_1 - \dot x_2) + K_{12}(x_1 - x_2) &= f_1(t) \\ M_2\ddot x_2 + B_2\dot x_2 + K_2x_2 + B_{12}(\dot x_2 - \dot x_1) + K_{12}(x_2 - x_1) &= 0 \end{aligned}

Laplace form

[M1s2+(B1+B12)s+K12]X1−(B12s+K12)X2=F1−(B12s+K12)X1+[M2s2+(B2+B12)s+K12+K2]X2=0\begin{aligned} \left[M_1s^2 + (B_1 + B_{12})s + K_{12}\right]X_1 - (B_{12}s + K_{12})X_2 &= F_1 \\ -(B_{12}s + K_{12})X_1 + \left[M_2s^2 + (B_2 + B_{12})s + K_{12} + K_2\right]X_2 &= 0 \end{aligned}

Transfer function

By Cramer's rule:

X1(s)F1(s)=M2s2+(B2+B12)s+K12+K2Δ\frac{X_1(s)}{F_1(s)} = \frac{M_2s^2 + (B_2 + B_{12})s + K_{12} + K_2}{\Delta} Δ=[M1s2+(B1+B12)s+K12][M2s2+(B2+B12)s+K12+K2]−(B12s+K12)2\Delta = \left[M_1s^2 + (B_1+B_{12})s + K_{12}\right]\left[M_2s^2 + (B_2+B_{12})s + K_{12} + K_2\right] - (B_{12}s + K_{12})^2

Expanded:

Δ=  M1M2s4+[B1M2+B2M1+B12(M1+M2)]s3+[B1B2+B1B12+B2B12+K12(M1+M2)+K2M1]s2+[B1(K12+K2)+B2K12+B12K2]s+K12K2\begin{aligned} \Delta =\; & M_1M_2s^4 + \left[B_1M_2 + B_2M_1 + B_{12}(M_1 + M_2)\right]s^3 \\ & + \left[B_1B_2 + B_1B_{12} + B_2B_{12} + K_{12}(M_1 + M_2) + K_2M_1\right]s^2 \\ & + \left[B_1(K_{12} + K_2) + B_2K_{12} + B_{12}K_2\right]s + K_{12}K_2 \end{aligned}

Check: at s=0s = 0, X1/F1=(K12+K2)/(K12K2)=1/K12+1/K2X_1/F_1 = (K_{12} + K_2)/(K_{12}K_2) = 1/K_{12} + 1/K_2, the compliance of the two springs in series.

  • 2066 Jestha (old course) · 8 marks

Draw free body diagram, write complete differential equations and find the transfer function X1(s)/F(s) for the dynamic system shown below. [Figure: wall – spring K1 – mass M1 (displacement x1, force f(t) applied) – damper B – mass M2 (displacement x2) – spring K2 – wall]

Answer

Displacements x1x_1 (M1M_1, force ff) and x2x_2 (M2M_2); no ground friction. K1K_1 ties M1M_1 to the left wall, BB couples the masses and K2K_2 ties M2M_2 to the right wall.

Free body diagrams

          M1                            M2
 K1x1  ┌──────┐ f(t)    B(x2'-x1') ┌──────┐
 ◄─────┤      ├──►      ◄──────────┤      │
 M1x1''│  M1  │         M2x2''     │  M2  │
 ◄─────┤      │         ◄──────────┤      │
 B(x1'-x2')   │         K2x2       │      │
 ◄─────┤      │         ◄──────────┤      │
       └──────┘                    └──────┘

Differential equations

M1x¨1+K1x1+B(x˙1−x˙2)=f(t)M2x¨2+K2x2+B(x˙2−x˙1)=0\begin{aligned} M_1\ddot x_1 + K_1x_1 + B(\dot x_1 - \dot x_2) &= f(t) \\ M_2\ddot x_2 + K_2x_2 + B(\dot x_2 - \dot x_1) &= 0 \end{aligned}

Laplace form

(M1s2+Bs+K1)X1−BsX2=F−BsX1+(M2s2+Bs+K2)X2=0\begin{aligned} (M_1s^2 + Bs + K_1)X_1 - BsX_2 &= F \\ -BsX_1 + (M_2s^2 + Bs + K_2)X_2 &= 0 \end{aligned}

From the second equation X2=BsM2s2+Bs+K2X1X_2 = \dfrac{Bs}{M_2s^2 + Bs + K_2}X_1. Substituting:

[(M1s2+Bs+K1)−B2s2M2s2+Bs+K2]X1=F\left[(M_1s^2 + Bs + K_1) - \frac{B^2s^2}{M_2s^2 + Bs + K_2}\right]X_1 = F

Transfer function

X1(s)F(s)=M2s2+Bs+K2(M1s2+Bs+K1)(M2s2+Bs+K2)−B2s2\frac{X_1(s)}{F(s)} = \frac{M_2s^2 + Bs + K_2}{(M_1s^2 + Bs + K_1)(M_2s^2 + Bs + K_2) - B^2s^2}

Expanding the denominator:

X1(s)F(s)=M2s2+Bs+K2M1M2s4+B(M1+M2)s3+(M1K2+M2K1)s2+B(K1+K2)s+K1K2\frac{X_1(s)}{F(s)} = \frac{M_2s^2 + Bs + K_2}{M_1M_2s^4 + B(M_1 + M_2)s^3 + (M_1K_2 + M_2K_1)s^2 + B(K_1 + K_2)s + K_1K_2}

Check: for a steady force (s→0s\to 0), X1=F/K1X_1 = F/K_1; the damper carries no static force, so K2K_2 does not share the load.

  • 2065 Shrawan (old course) · 8 marks

Derive the transfer function X1(s)/F(s) for the system shown below. [Figure: vertical system; mass M1 (displacement x1(t)) hangs from a fixed support through spring K1; mass M2 (displacement x2(t)) hangs below M1 through damper B1; damper B2 connects M2 to the fixed support above; force F(t) acts downward on M2]

Answer

Measure x1x_1, x2x_2 downward from the static equilibrium position, so gravity is balanced by the initial spring stretch and drops out of the equations.

Forces

  • M1M_1: spring K1K_1 (to support) gives K1x1K_1x_1; damper B1B_1 (to M2M_2) gives B1(x˙1−x˙2)B_1(\dot x_1 - \dot x_2); inertia M1x¨1M_1\ddot x_1.
  • M2M_2: applied F(t)F(t) downward; damper B2B_2 (to support) gives B2x˙2B_2\dot x_2; B1B_1 gives B1(x˙2−x˙1)B_1(\dot x_2 - \dot x_1); inertia M2x¨2M_2\ddot x_2.

Differential equations

M1x¨1+K1x1+B1(x˙1−x˙2)=0M2x¨2+B2x˙2+B1(x˙2−x˙1)=F(t)\begin{aligned} M_1\ddot x_1 + K_1x_1 + B_1(\dot x_1 - \dot x_2) &= 0 \\ M_2\ddot x_2 + B_2\dot x_2 + B_1(\dot x_2 - \dot x_1) &= F(t) \end{aligned}

Laplace form

(M1s2+B1s+K1)X1−B1sX2=0−B1sX1+[M2s2+(B1+B2)s]X2=F\begin{aligned} (M_1s^2 + B_1s + K_1)X_1 - B_1sX_2 &= 0 \\ -B_1sX_1 + \left[M_2s^2 + (B_1 + B_2)s\right]X_2 &= F \end{aligned}

Solving

From the first equation, X2=M1s2+B1s+K1B1sX1X_2 = \dfrac{M_1s^2 + B_1s + K_1}{B_1s}X_1. Substituting:

[(M1s2+B1s+K1) s(M2s+B1+B2)B1s−B1s]X1=F\left[\frac{(M_1s^2 + B_1s + K_1)\,s(M_2s + B_1 + B_2)}{B_1s} - B_1s\right]X_1 = F X1(s)F(s)=B1ss(M1s2+B1s+K1)(M2s+B1+B2)−B12s2\frac{X_1(s)}{F(s)} = \frac{B_1s}{s(M_1s^2 + B_1s + K_1)(M_2s + B_1 + B_2) - B_1^2s^2}

Cancelling ss:

X1(s)F(s)=B1(M1s2+B1s+K1)(M2s+B1+B2)−B12s\frac{X_1(s)}{F(s)} = \frac{B_1}{(M_1s^2 + B_1s + K_1)(M_2s + B_1 + B_2) - B_1^2s}

Final form

X1(s)F(s)=B1M1M2s3+[B1(M1+M2)+B2M1]s2+(B1B2+K1M2)s+K1(B1+B2)\frac{X_1(s)}{F(s)} = \frac{B_1}{M_1M_2s^3 + \left[B_1(M_1 + M_2) + B_2M_1\right]s^2 + (B_1B_2 + K_1M_2)s + K_1(B_1 + B_2)}

Check: under a constant force, M2M_2 finally moves at a steady velocity v2=F/(B1+B2)v_2 = F/(B_1 + B_2) (it is held only by dampers), and M1M_1 settles where K1x1=B1v2K_1x_1 = B_1v_2, i.e. x1=B1FK1(B1+B2)x_1 = \dfrac{B_1F}{K_1(B_1 + B_2)}, which is the value of the transfer function at s=0s = 0.

  • 2081 Bhadra · 4 marks

Draw block diagram of a general armature-controlled dc motor with feedback taken via a tachometer using P-controller and describe each variable at different points before or after the blocks used.

Answer

A tachometer is a small DC generator on the motor shaft whose output voltage is proportional to speed, vt=KT ωv_t = K_T\,\omega. Feeding it back and comparing with a reference voltage gives a closed-loop speed control system with a proportional (P) controller.

Block diagram

 Vr +   e  ┌──┐Va +   ┌──────┐Ia┌──┐Tm┌────┐ ω
 ──►(Σ)───►│Kp├─►(Σ)─►│  1   ├─►│Kt├─►│ 1  ├─┬──► ω
    ▲ -    └──┘   ▲ - │Ra+sLa│  └──┘  │Js+B│ │
    │             │   └──────┘        └────┘ │
    │             │ Eb      ┌──┐             │
    │             └─────────┤Kb├◄────────────┤
    │                       └──┘             │
    │ Vt           ┌──┐                      │
    └──────────────┤KT├◄─────────────────────┘
                   └──┘

Variables at each point

PointVariableMeaning
Reference inputVrV_rvoltage set for the desired speed
After first summere=Vr−Vte = V_r - V_tspeed error voltage
After KpK_pVa=KpeV_a = K_pearmature voltage from the P-controller/power amplifier
After second summerVa−EbV_a - E_bnet voltage across armature impedance
After 1Ra+sLa\frac{1}{R_a + sL_a}IaI_aarmature current
After KtK_tTm=KtIaT_m = K_tI_amotor torque
After 1Js+B\frac{1}{Js + B}ω\omegashaft speed (output)
Inner feedback KbK_bEb=KbωE_b = K_b\omegaback emf
Outer feedback KTK_TVt=KTωV_t = K_T\omegatachometer voltage

Closed-loop transfer function

With Gm(s)=Kt(Ra+sLa)(Js+B)+KtKbG_m(s) = \dfrac{K_t}{(R_a + sL_a)(Js + B) + K_tK_b} (motor with back emf):

Ω(s)Vr(s)=KpGm(s)1+KpKTGm(s)\frac{\Omega(s)}{V_r(s)} = \frac{K_pG_m(s)}{1 + K_pK_TG_m(s)}

Increasing KpK_p reduces the steady-state speed error and the effect of load-torque disturbances, but a P-controller alone always leaves a small steady-state error for a step speed command.

  • 2081 Bhadra · 6+2 marks

Find the transfer function X2(s)/F(s) of the mechanical system shown below by constructing the free body diagram and writing necessary mathematical equations. Also draw F-V analogy circuit. [Figure: force F applied to mass M1 (displacement x1, friction fc1 with ground); M1 connected to mass M2 (displacement x2, friction fc2 with ground) by two parallel paths: spring K1 in series with damper B1 (junction displacement x3), and spring K2 directly; M2 connected to the wall through damper B2]

Answer

Displacements: x1x_1 (M1M_1), x2x_2 (M2M_2), and x3x_3 = junction of K1K_1 (on the M1M_1 side) and B1B_1 (on the M2M_2 side). K2K_2 connects the masses directly.

Free body equations

  • M1M_1: FF applied; opposed by M1x¨1M_1\ddot x_1, fc1x˙1f_{c1}\dot x_1, K1(x1−x3)K_1(x_1 - x_3), K2(x1−x2)K_2(x_1 - x_2).
  • Junction (massless): K1(x3−x1)+B1(x˙3−x˙2)=0K_1(x_3 - x_1) + B_1(\dot x_3 - \dot x_2) = 0.
  • M2M_2: opposed by M2x¨2M_2\ddot x_2, (fc2+B2)x˙2(f_{c2} + B_2)\dot x_2, B1(x˙2−x˙3)B_1(\dot x_2 - \dot x_3), K2(x2−x1)K_2(x_2 - x_1).
M1x¨1+fc1x˙1+K1(x1−x3)+K2(x1−x2)=FK1(x3−x1)+B1(x˙3−x˙2)=0M2x¨2+(fc2+B2)x˙2+B1(x˙2−x˙3)+K2(x2−x1)=0\begin{aligned} M_1\ddot x_1 + f_{c1}\dot x_1 + K_1(x_1 - x_3) + K_2(x_1 - x_2) &= F \\ K_1(x_3 - x_1) + B_1(\dot x_3 - \dot x_2) &= 0 \\ M_2\ddot x_2 + (f_{c2} + B_2)\dot x_2 + B_1(\dot x_2 - \dot x_3) + K_2(x_2 - x_1) &= 0 \end{aligned}

Eliminating x3x_3

The series K1K_1–B1B_1 branch acts as Keq=K1B1sK1+B1sK_{eq} = \dfrac{K_1B_1s}{K_1 + B_1s}, in parallel with K2K_2. Total coupling:

Zc=K2+K1B1sK1+B1s=K1K2+B1(K1+K2)sK1+B1sZ_c = K_2 + \frac{K_1B_1s}{K_1 + B_1s} = \frac{K_1K_2 + B_1(K_1 + K_2)s}{K_1 + B_1s}

With a=M1s2+fc1sa = M_1s^2 + f_{c1}s and b=M2s2+(fc2+B2)sb = M_2s^2 + (f_{c2} + B_2)s:

(a+Zc)X1−ZcX2=F−ZcX1+(b+Zc)X2=0\begin{aligned} (a + Z_c)X_1 - Z_cX_2 &= F \\ -Z_cX_1 + (b + Z_c)X_2 &= 0 \end{aligned} X2F=Zcab+Zc(a+b)\frac{X_2}{F} = \frac{Z_c}{ab + Z_c(a + b)}

Transfer function

Multiplying by (K1+B1s)(K_1 + B_1s) and writing P=K1K2+B1(K1+K2)sP = K_1K_2 + B_1(K_1 + K_2)s:

X2(s)F(s)=K1K2+B1(K1+K2)ss{(K1+B1s) s (M1s+fc1)(M2s+fc2+B2)+P[(M1+M2)s+fc1+fc2+B2]}\frac{X_2(s)}{F(s)} = \frac{K_1K_2 + B_1(K_1 + K_2)s}{s\left\{(K_1 + B_1s)\,s\,(M_1s + f_{c1})(M_2s + f_{c2} + B_2) + P\left[(M_1 + M_2)s + f_{c1} + f_{c2} + B_2\right]\right\}}

The pole at s=0s = 0 appears because neither mass is tied to the wall by a spring.

F-V analogy

F→eF\to e, M→LM\to L, B,fc→RB, f_c\to R, K→1/CK\to 1/C. Mechanical elements in parallel (same velocity difference) become series elements, and elements in series (same force) become parallel. So the shared branch between mesh 1 (M1M_1) and mesh 2 (M2M_2) is C=1/K2C = 1/K_2 in series with the parallel pair 1/K1∥B11/K_1 \parallel B_1; the small loop inside that pair is mesh 3 (junction x3x_3).

   L1=M1  R=fc1           L2=M2  R=fc2+B2
 +-LLLL---/\/\--+--------LLLL---/\/\--+
 |              |                     |
(~) e        C=1/K2                   |
 |    i1        |           i2        |
 |          +---+---+                 |
 |        C=1/K1  R=B1                |
 |          | i3    |                 |
 |          +---+---+                 |
 |              |                     |
 +--------------+---------------------+

Mesh equations: L1i˙1+Rc1i1+1C2∫(i1−i2)+1C1∫(i1−i3)=eL_1\dot i_1 + R_{c1}i_1 + \frac{1}{C_2}\int(i_1 - i_2) + \frac{1}{C_1}\int(i_1 - i_3) = e; 1C1∫(i3−i1)+RB1(i3−i2)=0\frac{1}{C_1}\int(i_3 - i_1) + R_{B1}(i_3 - i_2) = 0; L2i˙2+(Rc2+RB2)i2+1C2∫(i2−i1)+RB1(i2−i3)=0L_2\dot i_2 + (R_{c2} + R_{B2})i_2 + \frac{1}{C_2}\int(i_2 - i_1) + R_{B1}(i_2 - i_3) = 0.

  • 2081 Baisakh · 10 marks

Find the transfer function θ1(s)/T(s) for the given rotational mechanical system. Also draw analogous electrical networks. [Figure: fixed wall – torsional spring K1 – inertia J1 (torque T(t) applied, angle θ1); J1 connected to inertia J2 (angle θ2) by damper B1 in series with spring K2, in parallel with spring K3; J2 – damper B2 – fixed wall]

Answer

Angles: θ1\theta_1 (J1J_1, torque TT), θ2\theta_2 (J2J_2) and θ3\theta_3 = junction of B1B_1 (on the J1J_1 side) and K2K_2 (on the J2J_2 side). K3K_3 couples J1J_1 and J2J_2 directly; K1K_1 and B2B_2 go to the frame.

Differential equations (Laplace form)

J1s2θ1+K1θ1+B1s(θ1−θ3)+K3(θ1−θ2)=TB1s(θ3−θ1)+K2(θ3−θ2)=0J2s2θ2+B2sθ2+K2(θ2−θ3)+K3(θ2−θ1)=0\begin{aligned} J_1s^2\theta_1 + K_1\theta_1 + B_1s(\theta_1 - \theta_3) + K_3(\theta_1 - \theta_2) &= T \\ B_1s(\theta_3 - \theta_1) + K_2(\theta_3 - \theta_2) &= 0 \\ J_2s^2\theta_2 + B_2s\theta_2 + K_2(\theta_2 - \theta_3) + K_3(\theta_2 - \theta_1) &= 0 \end{aligned}

Reducing the series branch

B1B_1 in series with K2K_2 is equivalent to B1sK2B1s+K2\dfrac{B_1sK_2}{B_1s + K_2}; with K3K_3 in parallel the total coupling is

Zc=K3+B1K2sB1s+K2=QB1s+K2,Q=K3(B1s+K2)+B1K2sZ_c = K_3 + \frac{B_1K_2s}{B_1s + K_2} = \frac{Q}{B_1s + K_2}, \quad Q = K_3(B_1s + K_2) + B_1K_2s

With a=J1s2+K1a = J_1s^2 + K_1 and b=J2s2+B2sb = J_2s^2 + B_2s:

(a+Zc)θ1−Zcθ2=T−Zcθ1+(b+Zc)θ2=0⇒θ1T=b+Zcab+Zc(a+b)\begin{aligned} (a + Z_c)\theta_1 - Z_c\theta_2 &= T \\ -Z_c\theta_1 + (b + Z_c)\theta_2 &= 0 \end{aligned} \qquad\Rightarrow\qquad \frac{\theta_1}{T} = \frac{b + Z_c}{ab + Z_c(a + b)}

Transfer function

Multiplying by (B1s+K2)(B_1s + K_2):

θ1(s)T(s)=(B1s+K2)(J2s2+B2s)+Q(B1s+K2)(J1s2+K1)(J2s2+B2s)+Q[(J1+J2)s2+B2s+K1]\frac{\theta_1(s)}{T(s)} = \frac{(B_1s + K_2)(J_2s^2 + B_2s) + Q}{(B_1s + K_2)(J_1s^2 + K_1)(J_2s^2 + B_2s) + Q\left[(J_1 + J_2)s^2 + B_2s + K_1\right]}

Expanded:

θ1T=B1J2s3+(B1B2+J2K2)s2+(B1K2+B1K3+B2K2)s+K2K3Δ\frac{\theta_1}{T} = \frac{B_1J_2s^3 + (B_1B_2 + J_2K_2)s^2 + (B_1K_2 + B_1K_3 + B_2K_2)s + K_2K_3}{\Delta} Δ=  B1J1J2s5+(B1B2J1+J1J2K2)s4+[B1J1(K2+K3)+B1J2(K1+K2+K3)+B2J1K2]s3+[B1B2(K1+K2+K3)+J1K2K3+J2K2(K1+K3)]s2+[B1K1(K2+K3)+B2K2(K1+K3)]s+K1K2K3\begin{aligned} \Delta =\; & B_1J_1J_2s^5 + (B_1B_2J_1 + J_1J_2K_2)s^4 \\ & + \left[B_1J_1(K_2 + K_3) + B_1J_2(K_1 + K_2 + K_3) + B_2J_1K_2\right]s^3 \\ & + \left[B_1B_2(K_1 + K_2 + K_3) + J_1K_2K_3 + J_2K_2(K_1 + K_3)\right]s^2 \\ & + \left[B_1K_1(K_2 + K_3) + B_2K_2(K_1 + K_3)\right]s + K_1K_2K_3 \end{aligned}

Analogous electrical networks

MechanicalT-V (mesh)T-I (node)
TTeeii
JJLLCC
BBRR1/R1/R
KK1/C1/C1/L1/L
ω\omegaiivv

T-I network (nodes v1,v3,v2v_1, v_3, v_2):

  v1                 v3              v2
 +--+---+---+-/\/\/--+---LLLL---+--+---+
 |  |   |   |  1/B1      1/K2   |  |   |
(i) C1  L1  +------LLLL---------+  C2  R2
 | =J1 =1/K1       1/K3          =J2 =1/B2
 +--+---+------------------------+--+---+ gnd

T-V network (mesh currents i1,i2i_1, i_2 and inner loop i3i_3):

   L1=J1  C1=1/K1          L2=J2  R2=B2
 +-LLLL----||----+---------LLLL---/\/\-+
 |               |                     |
(~) e         C3=1/K3                  |
 |    i1         |          i2         |
 |           +---+---+                 |
 |         R1=B1  C2=1/K2              |
 |           | i3    |                 |
 |           +---+---+                 |
 |               |                     |
 +---------------+---------------------+

K3K_3 (parallel coupling) becomes a series element C3C_3 in the shared branch, and the series pair B1B_1–K2K_2 becomes the parallel pair R1∥C2R_1 \parallel C_2. The mesh equations are exactly the three equations above with θ→q\theta \to q.

  • 2080 Bhadra · 8+2 marks

Find the transfer function X2(s)/F(s). Also draw the F-V analogous circuit. [Figure: force F applied to mass M1 (displacement x1, ground friction Bv1); M1 connected to M1's right side: damper B1 in series with spring K3 to mass M2 (displacement x2, ground friction Bv2); M1 also connected through spring K1 to the ground/fixed frame; M2 connected to the right wall through spring K2 and through damper B2 to the ground/fixed frame]

Answer

Displacements: x1x_1 (M1M_1), x2x_2 (M2M_2) and x3x_3 = junction between damper B1B_1 (on the M1M_1 side) and spring K3K_3 (on the M2M_2 side).

Equations of motion (Laplace)

(M1s2+Bv1s+K1)X1+B1s(X1−X3)=FB1s(X3−X1)+K3(X3−X2)=0[M2s2+(Bv2+B2)s+K2]X2+K3(X2−X3)=0\begin{aligned} (M_1s^2 + B_{v1}s + K_1)X_1 + B_1s(X_1 - X_3) &= F \\ B_1s(X_3 - X_1) + K_3(X_3 - X_2) &= 0 \\ \left[M_2s^2 + (B_{v2} + B_2)s + K_2\right]X_2 + K_3(X_2 - X_3) &= 0 \end{aligned}

Eliminating the junction

B1B_1 in series with K3K_3:

Keq(s)=B1sK3B1s+K3K_{eq}(s) = \frac{B_1sK_3}{B_1s + K_3}

With a=M1s2+Bv1s+K1a = M_1s^2 + B_{v1}s + K_1 and b=M2s2+(Bv2+B2)s+K2b = M_2s^2 + (B_{v2} + B_2)s + K_2:

X2F=Keq(a+Keq)(b+Keq)−Keq2=Keqab+Keq(a+b)\frac{X_2}{F} = \frac{K_{eq}}{(a + K_{eq})(b + K_{eq}) - K_{eq}^2} = \frac{K_{eq}}{ab + K_{eq}(a + b)}

Transfer function

Multiplying by (B1s+K3)(B_1s + K_3):

X2(s)F(s)=B1K3s(B1s+K3)(M1s2+Bv1s+K1)[M2s2+(Bv2+B2)s+K2]+B1K3s[(M1+M2)s2+(Bv1+Bv2+B2)s+K1+K2]\frac{X_2(s)}{F(s)} = \frac{B_1K_3s}{(B_1s + K_3)(M_1s^2 + B_{v1}s + K_1)\left[M_2s^2 + (B_{v2} + B_2)s + K_2\right] + B_1K_3s\left[(M_1 + M_2)s^2 + (B_{v1} + B_{v2} + B_2)s + K_1 + K_2\right]}

The denominator is 5th order (leading term B1M1M2s5B_1M_1M_2s^5, constant term K1K2K3K_1K_2K_3). The zero at s=0s = 0 shows that a steady force cannot move M2M_2, because the damper B1B_1 transmits no static force.

F-V analogous circuit

F→eF\to e, M→LM\to L, B→RB\to R, K→1/CK\to 1/C, velocity →\to mesh current. B1B_1 and K3K_3 are in series (same force), so in the F-V circuit they become the parallel pair R=B1∥C=1/K3R = B_1 \parallel C = 1/K_3, shared by mesh 1 (M1M_1) and mesh 2 (M2M_2); the loop inside the pair is mesh 3 (junction x3x_3).

  L1  R=Bv1 C=1/K1      L2  R=Bv2+B2 C=1/K2
 +LLLL-/\/\--||---+-----LLLL---/\/\----||--+
 |                |                        |
(~) e       +-----+-----+                  |
 |    i1    |           |      i2          |
 |        R=B1   i3   C=1/K3               |
 |          |           |                  |
 |          +-----+-----+                  |
 |                |                        |
 +----------------+------------------------+

Mesh equations:

L1di1dt+Rv1i1+1C1∫i1dt+R1(i1−i3)=eR1(i3−i1)+1C3∫(i3−i2)dt=0L2di2dt+(Rv2+R2)i2+1C2∫i2dt+1C3∫(i2−i3)dt=0\begin{aligned} L_1\frac{di_1}{dt} + R_{v1}i_1 + \frac{1}{C_1}\int i_1dt + R_1(i_1 - i_3) &= e \\ R_1(i_3 - i_1) + \frac{1}{C_3}\int(i_3 - i_2)dt &= 0 \\ L_2\frac{di_2}{dt} + (R_{v2} + R_2)i_2 + \frac{1}{C_2}\int i_2dt + \frac{1}{C_3}\int(i_2 - i_3)dt &= 0 \end{aligned}

(R1=B1R_1 = B_1 is shared by meshes 1 and 3; C3=1/K3C_3 = 1/K_3 by meshes 3 and 2.)

  • 2080 Baisakh · 8 marks

Find the transfer function of mechanical system given below taking output as velocity of mass M2; also draw F-V and F-I analogy electrical networks. [Figure: left wall – damper B1 – mass m1 (displacement x1) – spring k1 – mass m2 (displacement x2, force f(t) applied); m2 also connected to the left wall through spring k2; m2 connected to the right wall through damper B2 and spring k3 in parallel]

Answer

Displacements x1x_1 (m1m_1) and x2x_2 (m2m_2, force ff). Output: velocity of m2m_2, v2=x˙2v_2 = \dot x_2, so V2(s)=sX2(s)V_2(s) = sX_2(s).

Equations of motion

  • m1m_1: damper B1B_1 to the left wall, spring k1k_1 to m2m_2.
  • m2m_2: spring k1k_1 to m1m_1, k2k_2 to the left wall, B2∥k3B_2 \parallel k_3 to the right wall.
m1x¨1+B1x˙1+k1(x1−x2)=0m2x¨2+B2x˙2+(k2+k3)x2+k1(x2−x1)=f(t)\begin{aligned} m_1\ddot x_1 + B_1\dot x_1 + k_1(x_1 - x_2) &= 0 \\ m_2\ddot x_2 + B_2\dot x_2 + (k_2 + k_3)x_2 + k_1(x_2 - x_1) &= f(t) \end{aligned}

Laplace form:

(m1s2+B1s+k1)X1−k1X2=0−k1X1+(m2s2+B2s+k1+k2+k3)X2=F\begin{aligned} (m_1s^2 + B_1s + k_1)X_1 - k_1X_2 &= 0 \\ -k_1X_1 + (m_2s^2 + B_2s + k_1 + k_2 + k_3)X_2 &= F \end{aligned}

Transfer function

X2F=m1s2+B1s+k1Δ,Δ=(m1s2+B1s+k1)(m2s2+B2s+k1+k2+k3)−k12\frac{X_2}{F} = \frac{m_1s^2 + B_1s + k_1}{\Delta}, \qquad \Delta = (m_1s^2 + B_1s + k_1)(m_2s^2 + B_2s + k_1 + k_2 + k_3) - k_1^2

Velocity output:

V2(s)F(s)=s(m1s2+B1s+k1)Δ\frac{V_2(s)}{F(s)} = \frac{s(m_1s^2 + B_1s + k_1)}{\Delta}

with

Δ=  m1m2s4+(B1m2+B2m1)s3+[B1B2+m1(k1+k2+k3)+m2k1]s2+[B1(k1+k2+k3)+B2k1]s+k1(k2+k3)\begin{aligned} \Delta =\; & m_1m_2s^4 + (B_1m_2 + B_2m_1)s^3 + \left[B_1B_2 + m_1(k_1 + k_2 + k_3) + m_2k_1\right]s^2 \\ & + \left[B_1(k_1 + k_2 + k_3) + B_2k_1\right]s + k_1(k_2 + k_3) \end{aligned}

F-V analogy (mesh, velocity →\to current)

   L1=m1  R1=B1             L2=m2   R2=B2   C=1/(k2+k3)
 +-LLLL---/\/\--+----------LLLL----/\/\------||-----+
 |              |                                   |
 |           C1=1/k1                               (~) e
 |     i1       |              i2                   |
 +--------------+-----------------------------------+

(k2k_2 and k3k_3 both act on x2x_2 alone, so they combine into one capacitor 1/(k2+k3)1/(k_2+k_3), i.e. C2=1/k2C_2 = 1/k_2 and C3=1/k3C_3 = 1/k_3 in series.) The required output V2V_2 corresponds to mesh current i2i_2.

F-I analogy (node, velocity →\to voltage)

  v1          L1=1/k1           v2
 +--+---+------LLLL-------+---+---+---+---+
 |  |   |                 |   |   |   |   |
 C1 R1  |                 C2  R2  L2  L3 (i)
=m1 =1/B1                =m2 =1/B2 =1/k2 =1/k3
 +--+---+-----------------+---+---+---+---+ gnd

The output velocity v2v_2 is the voltage of node 2.

  • 2078 Bhadra · 8 marks

Find transfer function and develop F-V and F-I analogous circuits of the following figure. [Figure: mass M2 is a hollow frame resting on the ground (viscous friction B2 between M2 and ground), connected to the left wall through spring K2; mass M1 is a plate sliding inside M2 with viscous friction B1 between M1 and M2; force F(t) pulls M1 through spring K1]

Answer

Displacements: x1x_1 (plate M1M_1), x2x_2 (frame M2M_2). The force FF acts on the free end of spring K1K_1; since that end has no mass, the spring transmits the full force, K1(x0−x1)=FK_1(x_0 - x_1) = F, so FF effectively acts on M1M_1. (Here x0x_0 is the free-end displacement.)

Equations of motion

  • M1M_1: pulled by FF; viscous friction B1B_1 relative to M2M_2.
  • M2M_2: friction B1B_1 from M1M_1, ground friction B2B_2, spring K2K_2 to the wall.
M1x¨1+B1(x˙1−x˙2)=FM2x¨2+B1(x˙2−x˙1)+B2x˙2+K2x2=0\begin{aligned} M_1\ddot x_1 + B_1(\dot x_1 - \dot x_2) &= F \\ M_2\ddot x_2 + B_1(\dot x_2 - \dot x_1) + B_2\dot x_2 + K_2x_2 &= 0 \end{aligned}

Laplace form:

(M1s2+B1s)X1−B1sX2=F−B1sX1+[M2s2+(B1+B2)s+K2]X2=0\begin{aligned} (M_1s^2 + B_1s)X_1 - B_1sX_2 &= F \\ -B_1sX_1 + \left[M_2s^2 + (B_1 + B_2)s + K_2\right]X_2 &= 0 \end{aligned}

Transfer functions

Determinant (after cancelling a common ss):

Δ′=(M1s+B1)[M2s2+(B1+B2)s+K2]−B12s=M1M2s3+[B1(M1+M2)+B2M1]s2+(B1B2+K2M1)s+B1K2\Delta' = (M_1s + B_1)\left[M_2s^2 + (B_1 + B_2)s + K_2\right] - B_1^2s = M_1M_2s^3 + \left[B_1(M_1 + M_2) + B_2M_1\right]s^2 + (B_1B_2 + K_2M_1)s + B_1K_2 X2(s)F(s)=B1Δ′\frac{X_2(s)}{F(s)} = \frac{B_1}{\Delta'} X1(s)F(s)=M2s2+(B1+B2)s+K2s Δ′\frac{X_1(s)}{F(s)} = \frac{M_2s^2 + (B_1 + B_2)s + K_2}{s\,\Delta'}

and for the point where the force is applied, X0(s)=X1(s)+F(s)/K1X_0(s) = X_1(s) + F(s)/K_1.

F-V analogous circuit

F→eF\to e, M→LM\to L, B→RB\to R, K→1/CK\to 1/C. Mesh i0i_0 = free-end velocity, i1i_1 = velocity of M1M_1, i2i_2 = velocity of M2M_2. K1K_1 is shared by meshes 0 and 1; B1B_1 (relative friction) by meshes 1 and 2.

              L1=M1             L2=M2  R2=B2  C2=1/K2
 +------+-----LLLL----+---------LLLL---/\/\----||--+
 |      |             |                            |
(~) e  C=1/K1       R1=B1                          |
 |  i0  |      i1     |            i2              |
 +------+-------------+----------------------------+

F-I analogous circuit

F→iF\to i, M→CM\to C, B→1/RB\to 1/R, K→1/LK\to 1/L. Node 0 = free end, node 1 = M1M_1, node 2 = M2M_2; C1=M1C_1 = M_1, C2=M2C_2 = M_2, R2=1/B2R_2 = 1/B_2, L2=1/K2L_2 = 1/K_2.

  v0    L=1/K1    v1   R1=1/B1    v2
 +---+---LLLL---+---+---/\/\/---+---+---+---+
 |                  |               |   |   |
(i)                C1              C2  R2  L2
 |                  |               |   |   |
 +------------------+---------------+---+---+ gnd
  • 2078 Kartik · 4+2+2 marks

Find transfer function for the system as in figure. Also develop F-V and F-I analogy circuits. [Figure: force f(t) applied at a free end (displacement x1(t)) connected through spring K1 to mass M2 (displacement x2(t)); M2 connected to the fixed wall through spring K2 in parallel with damper B2]

Answer

Displacements: x1x_1 (free end where ff acts) and x2x_2 (mass M2M_2). The free end has no mass.

Equations of motion

At the massless end, the spring force equals the applied force:

K1(x1−x2)=f(t)K_1(x_1 - x_2) = f(t)

For M2M_2:

M2x¨2+B2x˙2+K2x2+K1(x2−x1)=0M_2\ddot x_2 + B_2\dot x_2 + K_2x_2 + K_1(x_2 - x_1) = 0

Transfer functions

Laplace form:

K1X1−K1X2=F−K1X1+(M2s2+B2s+K1+K2)X2=0\begin{aligned} K_1X_1 - K_1X_2 &= F \\ -K_1X_1 + (M_2s^2 + B_2s + K_1 + K_2)X_2 &= 0 \end{aligned}

Adding the two equations: (M2s2+B2s+K2)X2=F(M_2s^2 + B_2s + K_2)X_2 = F, so

X2(s)F(s)=1M2s2+B2s+K2\frac{X_2(s)}{F(s)} = \frac{1}{M_2s^2 + B_2s + K_2}

and since X1=X2+F/K1X_1 = X_2 + F/K_1:

X1(s)F(s)=M2s2+B2s+K1+K2K1(M2s2+B2s+K2)\frac{X_1(s)}{F(s)} = \frac{M_2s^2 + B_2s + K_1 + K_2}{K_1(M_2s^2 + B_2s + K_2)}

The spring K1K_1 only passes the force on; it does not change the dynamics of M2M_2.

F-V analogy

f→ef\to e, M→LM\to L, B→RB\to R, K→1/CK\to 1/C. Mesh 1 (velocity x˙1\dot x_1) and mesh 2 (x˙2\dot x_2) share C1=1/K1C_1 = 1/K_1.

                   L2=M2   R2=B2   C2=1/K2
 +--------+--------LLLL----/\/\-----||---+
 |        |                              |
(~) e   C1=1/K1                          |
 |   i1   |              i2              |
 +--------+------------------------------+
1C1∫(i1−i2)dt=e,L2di2dt+R2i2+1C2∫i2dt+1C1∫(i2−i1)dt=0\frac{1}{C_1}\int(i_1 - i_2)dt = e, \qquad L_2\frac{di_2}{dt} + R_2i_2 + \frac{1}{C_2}\int i_2dt + \frac{1}{C_1}\int(i_2 - i_1)dt = 0

F-I analogy

f→if\to i, M→CM\to C, B→1/RB\to 1/R, K→1/LK\to 1/L.

  v1      L1=1/K1      v2
 +--+------LLLL------+---+---+---+
 |                   |   |   |   |
(i)                 C2  R2  L2   |
 |                 =M2 =1/B2 =1/K2
 +-------------------+---+---+---+ gnd
1L1∫(v1−v2)dt=i,C2dv2dt+v2R2+1L2∫v2dt+1L1∫(v2−v1)dt=0\frac{1}{L_1}\int(v_1 - v_2)dt = i, \qquad C_2\frac{dv_2}{dt} + \frac{v_2}{R_2} + \frac{1}{L_2}\int v_2dt + \frac{1}{L_1}\int(v_2 - v_1)dt = 0

Questions from Old Question Collection (EE 602) (IOE Control System exam papers (EE 602 and older course), 2065 to 2082) and Old Question Collection (BEI, EE 504) (IOE BEI Control System (EE 504) exam papers, 2076 to 2081). Answers are written for this site; check them against your class notes.

Chapter titles and hours from the IOE syllabus ↗