PracticeVideosBlogFlashcards

13.2 Basic Mechatronics, Sensors & Closed-Loop Controls

Key Takeaways

  • Resistive strain gages quantify mechanical strain via ΔRR=GF⋅ϵ\frac{\Delta R}{R} = GF \cdot \epsilon; configured in a full Wheatstone bridge with four active gages, output voltage quadruples (Vout=Vs⋅GF⋅ϵV_{out} = V_s \cdot GF \cdot \epsilon), maximizing sensitivity while inherently canceling thermal expansion and bending/axial crosstalk.

  • Incremental optical encoders generate quadrature pulse trains (channels A and B shifted by 90∘90^\circ electrical) enabling 4×4\times edge-count resolution multiplication and direction sensing, whereas absolute encoders utilize Gray code tracks to prevent multi-bit transition errors during dynamic reading.

  • Brushed and brushless DC motors exhibit a linear inverse torque-speed relationship ω(T)=ωnl(1−T/Tstall)\omega(T) = \omega_{nl}(1 - T/T_{stall}), producing peak mechanical output power at exactly half the stall torque and half the no-load speed (Pmax⁡=14TstallωnlP_{\max} = \frac{1}{4} T_{stall} \omega_{nl}).

  • In PID control architectures, proportional gain KpK_p reduces rise time and steady-state error, integral gain KiK_i eliminates steady-state offset (ess→0e_{ss} \to 0) by increasing system type at the expense of phase lag and overshoot, and derivative gain KdK_d introduces phase lead to dampen transient oscillations.

  • Closed-loop system stability requires all roots of the characteristic equation 1+C(s)G(s)=01 + C(s)G(s) = 0 to lie strictly within the open left-half complex plane (Re(s)<0\text{Re}(s) < 0), verifiable algebraically via the Routh-Hurwitz stability criterion without factoring high-order polynomials.

Last updated: October 2026

13.2 Basic Mechatronics, Sensors & Closed-Loop Controls

Modern mechanical systems are rarely purely mechanical; they are integrated mechatronic systems that combine precision mechanisms, sensors, electronic signal conditioning, embedded digital processors, and electromechanical actuators. For machine designers, understanding how physical quantities (force, strain, displacement, temperature, and speed) are measured, converted into electrical signals, and regulated by feedback controllers is essential for designing automated machinery, robotic arms, precision positioning stages, and dynamic test equipment.


Mechatronic System Architecture & Signal Conditioning

A complete mechatronic control system functions as a closed information loop:

                      Mechatronic Control Loop
                      
         Reference  +    Error     +------------+  Control  +----------+
  r(t) ------------>(+)----------> | Controller | --------> | Actuator |
         (Setpoint)  ^ -  e(t)     | (DSP/PLC)  |   u(t)    | Driver   |
                     |             +------------+           +----------+
                     |                                            |
                     |           +-----------------------+        v
                     |           |   Mechanical Plant    | <--------------+
                     |           | (Mechanism / Load)    |   Driven       |
                     |           +-----------------------+   Motion       |
                     |                       |                            |
                     |                       v                            |
                     |             +-------------------+                  |
                     +------------ | Sensor / Signal   | <----------------+
                      Measurement  | Conditioning      |   Process Feedback
                          y(t)     +-------------------+

Signal Conditioning & Interfacing Fundamentals

  1. Amplification & Common-Mode Rejection: Sensors frequently output low-level differential voltages (1 mV−50 mV1\,\text{mV} - 50\,\text{mV}). Differential instrumentation amplifiers provide high input impedance and a high Common-Mode Rejection Ratio (CMRR) to amplify microvolt sensor signals while rejecting 60 Hz electrical mains noise and electromagnetic interference (EMI).

  2. Anti-Aliasing Filtering & Sampling (Nyquist Criterion): Analog signals must be converted to digital numbers via an Analog-to-Digital Converter (ADC). To prevent aliasing (where high-frequency noise folds back into the measurement band as false low-frequency signals), an analog low-pass anti-aliasing filter must precede the ADC. By the Nyquist-Shannon sampling theorem, the sampling frequency fsf_s must exceed at least twice the maximum frequency component present in the analog signal:

    fs≥2fmax⁡f_s \ge 2 f_{\max}

    In practical machine control systems, fsf_s is commonly set to 10×10 \times to 20×20 \times the closed-loop system bandwidth.

  3. ADC Resolution & Quantization Error: An NN-bit ADC with an input voltage range VrefV_{\text{ref}} divides the signal into 2N2^N discrete digital levels. The voltage resolution (quantization step size QQ) is:

    Q=Vref2NQ = \frac{V_{\text{ref}}}{2^N}

    The maximum quantization error is ±Q2\pm \frac{Q}{2}. For example, a 16-bit ADC operating over a 0−10 V0-10\,\text{V} range has Q=1065, ⁣536≈0.153 mVQ = \frac{10}{65,\!536} \approx 0.153\,\text{mV}.


Sensors & Transducers

1. Resistive Strain Gages & Wheatstone Bridges

A resistive foil strain gage consists of a thin metallic grid bonded to a flexible backing. When subjected to mechanical strain ϵ\epsilon, its electrical resistance changes proportionally due to physical elongation and piezoresistive effects:

ΔRR=GF⋅ϵ\frac{\Delta R}{R} = GF \cdot \epsilon

where:

  • RR = nominal gage resistance (commonly 120 Ω,350 Ω,or 1, ⁣000 Ω120\,\Omega, 350\,\Omega, \text{or } 1,\!000\,\Omega).
  • ϵ=ΔLL\epsilon = \frac{\Delta L}{L} = engineering strain (dimensionless, commonly expressed in microstrain με=10−6\mu\varepsilon = 10^{-6}).
  • GFGF = Gage Factor (dimensionless sensitivity factor; typically GF≈2.0−2.1GF \approx 2.0 - 2.1 for metallic constantan foil gages, and GF≈100−150GF \approx 100 - 150 for semiconductor silicon gages).
                     Wheatstone Bridge Circuit
                                (+) V_s (Excitation)
                                    *
                                   / \
                       R_1 (Arm 1)/   \ R_2 (Arm 2)
                                 /     \
                                *-------* 
                             (-)  V_out  (+)
                                 \     /
                       R_4 (Arm 4)\   / R_3 (Arm 3)
                                   \ /
                                    *
                                   GND

Wheatstone Bridge Circuit Formulations

Because ΔR/R\Delta R / R is tiny (typically <0.1%< 0.1\%), a Wheatstone bridge converts small fractional resistance changes into measurable differential voltages VoutV_{\text{out}}. The general balanced bridge relation is R1R3=R2R4R_1 R_3 = R_2 R_4. For small variations, the bridge output is:

VoutVs≈14(ΔR1R1−ΔR2R2+ΔR3R3−ΔR4R4)=GF4(ϵ1−ϵ2+ϵ3−ϵ4)\frac{V_{\text{out}}}{V_s} \approx \frac{1}{4} \left( \frac{\Delta R_1}{R_1} - \frac{\Delta R_2}{R_2} + \frac{\Delta R_3}{R_3} - \frac{\Delta R_4}{R_4} \right) = \frac{GF}{4} (\epsilon_1 - \epsilon_2 + \epsilon_3 - \epsilon_4)
  • Quarter-Bridge (1 active gage in arm 1, 3 fixed resistors):

    Vout=VsGF⋅ϵ4V_{\text{out}} = V_s \frac{GF \cdot \epsilon}{4}

    Susceptible to thermal drift because temperature changes cause thermal resistance changes that mimic mechanical strain. A dummy unstrained gage in arm 2 exposed to the same ambient temperature cancels temperature drift.

  • Half-Bridge (2 active gages):

    • Bending Application: Gage 1 mounted on top surface in tension (+ϵ+\epsilon), Gage 2 on bottom surface in compression (−ϵ-\epsilon): Vout=VsGF4[(ϵ)−(−ϵ)]=VsGF⋅ϵ2V_{\text{out}} = V_s \frac{GF}{4} [(\epsilon) - (-\epsilon)] = V_s \frac{GF \cdot \epsilon}{2} Doubles sensitivity, and inherently cancels both axial tension/compression and thermal expansion.
  • Full-Bridge (4 active gages):

    • Bending Beam: Gages 1 and 3 in tension (+ϵ+\epsilon), Gages 2 and 4 in compression (−ϵ-\epsilon): Vout=VsGF4[(ϵ)−(−ϵ)+(ϵ)−(−ϵ)]=Vs⋅GF⋅ϵV_{\text{out}} = V_s \frac{GF}{4} [(\epsilon) - (-\epsilon) + (\epsilon) - (-\epsilon)] = V_s \cdot GF \cdot \epsilon Quadruples output sensitivity, completely rejects temperature drift, cancels axial loads, and provides maximum signal-to-noise ratio.

2. Linear Variable Differential Transformers (LVDTs)

An LVDT is an electromechanical transducer that produces an electrical output proportional to linear displacement:

  • Construction: One primary coil centered between two identical secondary coils wound in series-opposition around a cylindrical hollow bobbin. A moveable high-permeability ferromagnetic core attaches to the moving machine member.
  • Operation: The primary coil is excited with an AC carrier voltage (1−10 kHz1-10\,\text{kHz}, 3−15 Vrms3-15\,\text{V}_{\text{rms}}). When the core is centered at the null position, equal magnetic flux links both secondary coils; because they are wired in series-opposition, the net output voltage is zero (Vout=Vs1−Vs2=0V_{\text{out}} = V_{s1} - V_{s2} = 0).
  • Displacement: Moving the core from null increases flux in one secondary while decreasing flux in the other. The amplitude of VoutV_{\text{out}} scales linearly with core displacement, while its phase (0∘0^\circ vs 180∘180^\circ relative to primary excitation) indicates direction.
  • Features: Frictionless non-contact operation, infinite theoretical resolution, high linearity (<0.25%< 0.25\%), and rugged resistance to shock, vibration, and contaminated environments.

3. Optical Encoders (Incremental vs. Absolute)

Optical encoders convert rotational displacement into digital pulses using a patterned glass or polymer disc and photodetectors.

  • Incremental Encoders:

    • Utilize two optical tracks: Channel A and Channel B, offset by 90∘90^\circ electrical phase (quadrature).
    • Direction Sensing: If Channel A leads Channel B, shaft rotation is clockwise; if B leads A, rotation is counterclockwise.
    • Resolution Multiplication (4×4\times Decoding): By counting both rising and falling edges of channels A and B, the effective resolution is quadrupled (4×PPR4 \times \text{PPR} pulses per revolution).
    • Index Channel (Z): Produces one pulse per complete revolution, serving as a physical reference marker for homing.
    • Limitation: Loses absolute shaft position upon loss of power; requires a homing routine on power-up.
  • Absolute Encoders:

    • Disc features multiple concentric tracks, where each angular sector corresponds to a unique digital binary word.
    • Retains absolute position even during power outages; no homing sequence required.
    • Gray Code vs. Binary: In natural binary code, transitioning between consecutive numbers often requires multiple bits to flip simultaneously (e.g., 0112→1002011_2 \to 100_2 flips all 3 bits). Slight mechanical disc misalignment causes intermediate read glitches (e.g., momentarily reading 1112111_2). Gray code is a reflected binary code designed so that only one bit changes between any two consecutive values, eliminating transition glitches.

4. Temperature Sensors

                               Temperature Sensor Comparison
  +-----------------------------------------------------------------------------------------+
  | Sensor Type  | Transduction Principle          | Operating Range     | Key Characteristic |
  |--------------+---------------------------------+---------------------+--------------------|
  | Thermocouple | Seebeck thermoelectric effect   | -200C to +1800C     | Rugged, self-power |
  | RTD (Pt100)  | Positive temp-coefficient (PTC) | -200C to +650C      | High accuracy/line |
  | Thermistor   | Negative temp-coefficient (NTC) | -50C to +150C       | High sensitivity   |
  +-----------------------------------------------------------------------------------------+
  • Thermocouples:

    • Two dissimilar metal wires joined at one end form a measuring (hot) junction. The Seebeck effect generates an open-circuit voltage proportional to the temperature difference between hot junction ThotT_{\text{hot}} and reference (cold) junction TrefT_{\text{ref}}: V=SAB(Thot−Tref)V = S_{AB}(T_{\text{hot}} - T_{\text{ref}}).
    • Requires Cold-Junction Compensation (CJC) using a local IC temperature sensor on the isothermal terminal block.
    • Common types: Type K (Chromel-Alumel, wide industrial range −200∘C-200^\circ\text{C} to +1250∘C+1250^\circ\text{C}, ≈41 μV/∘C\approx 41\,\mu\text{V}/^\circ\text{C}), Type J (Iron-Constantan), Type T (Copper-Constantan).
  • Resistance Temperature Detectors (RTDs):

    • Operating principle: Pure metals exhibit a linear increase in electrical resistance with temperature: R(T)=R0[1+α(T−T0)]R(T) = R_0 [1 + \alpha(T - T_0)]
    • Industrial standard: Pt100 (R0=100.0 ΩR_0 = 100.0\,\Omega at 0∘C0^\circ\text{C}, α=0.00385 Ω/Ω/∘C\alpha = 0.00385\,\Omega/\Omega/^\circ\text{C}) and Pt1000 (1, ⁣000 Ω1,\!000\,\Omega at 0∘C0^\circ\text{C}).
    • Highly linear, stable, and accurate (±0.1∘C\pm 0.1^\circ\text{C}). 3-wire and 4-wire lead configurations cancel lead-wire resistance in remote installations.
  • Thermistors:

    • Semiconductor metal oxides with a Negative Temperature Coefficient (NTC): R(T)=R0exp⁡[β(1T−1T0)]R(T) = R_0 \exp\left[ \beta \left(\frac{1}{T} - \frac{1}{T_0}\right) \right]
    • Extremely sensitive to small temperature changes over narrow ranges, but highly non-linear.

Basic Circuit Relations for Mechatronics

Sensor, actuator and controller hardware rests on a few circuit relations:

  • Ohm's law and power: V=IRV = IR; P=VI=I2R=V2/RP = VI = I^2 R = V^2/R.
  • Series and parallel resistors: Req=R1+R2+…R_{eq} = R_1 + R_2 + \dots in series; 1/Req=1/R1+1/R2+…1/R_{eq} = 1/R_1 + 1/R_2 + \dots in parallel.
  • Voltage divider: Vout=VinR2/(R1+R2)V_{\text{out}} = V_{\text{in}} R_2/(R_1 + R_2), used to scale signals and in potentiometer position sensors.
  • RC circuits: time constant τ=RC\tau = RC; a capacitor charges to 63.2%63.2\% of its final voltage in one τ\tau and about 99%99\% in 5τ5\tau. A first-order RC low-pass filter has cutoff frequency fc=1/(2πRC)f_c = 1/(2\pi R C) and is the simplest anti-aliasing filter.
  • RL circuits: τ=L/R\tau = L/R. Current in solenoids and motor windings rises with this time constant, which is why stepper torque falls at high step rates.
  • Inductive loads (relays, solenoids, motor windings) produce a large voltage spike (V=L di/dtV = L\,di/dt) when switched off, so a flyback (freewheeling) diode is placed across the coil to protect the switching transistor.
  • Switching and drive circuits: transistors and MOSFETs let low-power logic switch loads; an H-bridge of four switches reverses a DC motor; pulse-width modulation (PWM) sets the average voltage Vavg=D VsupplyV_{\text{avg}} = D\,V_{\text{supply}} for duty cycle DD. Relays, contactors and solid-state relays interface controller outputs with motors and heaters.

Mini-example: A 10 kΩ10\,\text{k}\Omega resistor and a 0.1 μF0.1\,\mu\text{F} capacitor give τ=RC=1 ms\tau = RC = 1\,\text{ms} and fc=1/(2π×10,000×0.1×10−6)=159 Hzf_c = 1/(2\pi \times 10{,}000 \times 0.1 \times 10^{-6}) = 159\,\text{Hz}. A 24 V24\,\text{V} supply switched at 40%40\% duty cycle delivers an average of 9.6 V9.6\,\text{V} to a motor.


Electromechanical Actuators

1. DC Motors (Brushed vs. Brushless)

In a direct current (DC) motor, magnetic fields interact to produce torque:

V=IaRa+Eb=IaRa+KeωV = I_a R_a + E_b = I_a R_a + K_e \omega T=KtIaT = K_t I_a

where:

  • VV = applied armature terminal voltage (V)
  • IaI_a = armature current (A)
  • RaR_a = armature winding resistance (Ω\Omega)
  • Eb=KeωE_b = K_e \omega = counter-electromotive force / back-EMF (V)
  • ω\omega = rotor angular velocity (rad/s)
  • KeK_e = back-EMF constant (V⋅s/rad\text{V}\cdot\text{s/rad})
  • KtK_t = torque sensitivity constant (N⋅m/A\text{N}\cdot\text{m/A})

(In consistent SI units, Ke=KtK_e = K_t numerically).

Torque-Speed Characteristic

Combining the electrical and torque equations yields the linear torque-speed curve for a permanent magnet DC motor:

ω(T)=VKe−RaKtKeT=ωnl(1−TTstall)\omega(T) = \frac{V}{K_e} - \frac{R_a}{K_t K_e} T = \omega_{nl} \left(1 - \frac{T}{T_{\text{stall}}}\right)
  • No-Load Speed (ωnl\omega_{nl} at T=0T = 0): ωnl=VKe\omega_{nl} = \frac{V}{K_e}
  • Stall Torque (TstallT_{\text{stall}} at ω=0\omega = 0): Tstall=KtVRaT_{\text{stall}} = \frac{K_t V}{R_a}
                  DC Motor Torque-Speed Curve
       Speed w ^
          w_nl +\ 
               | \ 
               |  \   Operating Characteristic Curve
               |   \ 
        w_nl/2 +----+----\  <-- Maximum Power Point: P_max
               |    |     \ 
               +----+------+---------> Torque T
               0  T_stall/2 T_stall

Mechanical Output Power

The mechanical output power is P(T)=TωP(T) = T \omega:

P(T)=Tωnl(1−TTstall)P(T) = T \omega_{nl} \left(1 - \frac{T}{T_{\text{stall}}}\right)

Differentiating with respect to TT and setting to zero demonstrates that maximum mechanical power output occurs at exactly half the stall torque and half the no-load speed:

TPmax⁡=12Tstall,ωPmax⁡=12ωnlT_{P\max} = \frac{1}{2} T_{\text{stall}}, \qquad \omega_{P\max} = \frac{1}{2} \omega_{nl} Pmax⁡=TPmax⁡ωPmax⁡=14Tstallωnl=V24RaP_{\max} = T_{P\max} \omega_{P\max} = \frac{1}{4} T_{\text{stall}} \omega_{nl} = \frac{V^2}{4 R_a}

Brushed vs. Brushless DC Motors (BLDC)

  • Brushed DC: Mechanical commutator rings and graphite brushes continuously reverse armature coil currents. Simple two-lead drive, but subject to brush arcing, wear, acoustic/RF noise, and limited operational life.
  • Brushless DC (BLDC): Permanent magnets on rotor, multiphase electromagnetic windings on stator. Electronic inverter commutates currents based on rotor position (Hall effect sensors or sensorless back-EMF sensing). Offers higher efficiency, superior thermal dissipation, higher continuous torque, and virtually maintenance-free operation.

2. Stepper Motors

Stepper motors convert discrete digital pulses into precise incremental shaft rotations without closed-loop feedback.

  • Step Angle (θs\theta_s): θs=360∘2 m Z\theta_s = \frac{360^\circ}{2\,m\,Z} for a hybrid stepper driven in full steps, where mm = number of phases (typically 2) and ZZ = number of rotor teeth (commonly 50). Each tooth pitch is crossed in 2m2m full steps, giving the standard hybrid stepper angle: θs=360∘2×2×50=1.8∘(  ⟹  200 steps per revolution)\theta_s = \frac{360^\circ}{2 \times 2 \times 50} = 1.8^\circ \quad (\implies 200\,\text{steps per revolution})
  • Microstepping: Modulating phase currents using Pulse-Width Modulation (PWM) into sinusoidal proportions divides each fundamental full step into fractions (1/4,1/16,1/64,1/2561/4, 1/16, 1/64, 1/256), smoothing shaft torque and suppressing mid-frequency mechanical resonance.
  • Torque Behavior: Stepper motors generate high holding torque at zero speed, but usable torque drops rapidly at higher stepping frequencies due to winding inductance (L/RL/R time constant) limiting current rise.

Closed-Loop Feedback Control Systems

Open-Loop vs. Closed-Loop Systems

  • Open-Loop: Control input is applied without measuring system output. Simple, low cost, unconditionally stable, but unable to compensate for external disturbances (load changes, friction variations) or plant parameter drift.
  • Closed-Loop (Feedback): Sensor monitors the output y(t)y(t), subtracts it from reference setpoint r(t)r(t) to generate error signal e(t)=r(t)−y(t)e(t) = r(t) - y(t). Provides high disturbance rejection, reduced sensitivity to plant parameter changes, and stabilization of unstable systems.

Block Diagram Algebra & Transfer Functions

Consider a negative feedback system with forward controller C(s)C(s), plant G(s)G(s), and unity feedback (H(s)=1H(s) = 1):

T(s)=Y(s)R(s)=C(s)G(s)1+C(s)G(s)T(s) = \frac{Y(s)}{R(s)} = \frac{C(s) G(s)}{1 + C(s) G(s)}

The open-loop loop transfer function is L(s)=C(s)G(s)L(s) = C(s) G(s). The characteristic equation governing system stability is:

1+L(s)=1+C(s)G(s)=01 + L(s) = 1 + C(s) G(s) = 0

PID Control Fundamentals

The Proportional-Integral-Derivative (PID) controller is the most widely deployed feedback algorithm in machine design:

u(t)=Kpe(t)+Ki∫0te(τ)dτ+Kdde(t)dtu(t) = K_p e(t) + K_i \int_0^t e(\tau) d\tau + K_d \frac{de(t)}{dt}

In the Laplace domain:

C(s)=Kp+Kis+Kds=Kds2+Kps+KisC(s) = K_p + \frac{K_i}{s} + K_d s = \frac{K_d s^2 + K_p s + K_i}{s}
                                PID Tuning Influences
  +----------------------------------------------------------------------------------------+
  | Parameter | Rise Time     | Overshoot     | Settling Time | Steady-State Error | Stability |
  |-----------+---------------+---------------+---------------+--------------------+-----------|
  | Increase K_p | Decreases  | Increases     | Minor Effect  | Decreases (offset) | Degrades  |
  | Increase K_i | Decreases  | Increases     | Increases     | Eliminates to ZERO | Degrades  |
  | Increase K_d | Minor Effect| Decreases     | Decreases     | Minor Effect       | Improves  |
  +----------------------------------------------------------------------------------------+
  • Proportional Gain (KpK_p): Provides immediate restoring action proportional to current error. Accelerates response and reduces steady-state error, but high KpK_p induces excessive overshoot and oscillation.
  • Integral Gain (KiK_i): Integrates past error over time. Adds a free integrator (pole at s=0s = 0), raising system type by 1 and forcing steady-state error to zero (ess=0e_{ss} = 0) under step commands and constant disturbance loads. However, introduces −90∘-90^\circ phase lag, which reduces damping and can cause integral windup during actuator saturation.
  • Derivative Gain (KdK_d): Anticipates future error based on rate of change (de/dtde/dt). Adds phase lead, increasing effective damping, suppressing overshoot, and accelerating settling time. However, amplifies high-frequency measurement noise; practical implementations require a first-order low-pass filter on the derivative term: Cd(s)=Kdsτfs+1C_d(s) = \frac{K_d s}{\tau_f s + 1}.

Stability Analysis & Routh-Hurwitz Stability Criterion

A closed-loop linear system is Bounded-Input Bounded-Output (BIBO) stable if and only if all roots of the characteristic equation lie strictly in the open left-half complex plane (Re(s)<0\text{Re}(s) < 0).

For a general characteristic polynomial:

A(s)=ansn+an−1sn−1+an−2sn−2+⋯+a1s+a0=0A(s) = a_n s^n + a_{n-1} s^{n-1} + a_{n-2} s^{n-2} + \dots + a_1 s + a_0 = 0
  1. Necessary Condition (Stodola): All coefficients aia_i must be strictly positive (ai>0a_i > 0) and non-zero.
  2. Routh-Hurwitz Array: Form the array from polynomial coefficients: snanan−2an−4…sn−1an−1an−3an−5…sn−2b1b2b3…sn−3c1c2c3…\begin{array}{c|cccc} s^n & a_n & a_{n-2} & a_{n-4} & \dots \\ s^{n-1} & a_{n-1} & a_{n-3} & a_{n-5} & \dots \\ s^{n-2} & b_1 & b_2 & b_3 & \dots \\ s^{n-3} & c_1 & c_2 & c_3 & \dots \end{array} where: b1=an−1an−2−anan−3an−1,b2=an−1an−4−anan−5an−1b_1 = \frac{a_{n-1} a_{n-2} - a_n a_{n-3}}{a_{n-1}}, \qquad b_2 = \frac{a_{n-1} a_{n-4} - a_n a_{n-5}}{a_{n-1}} c1=b1an−3−an−1b2b1c_1 = \frac{b_1 a_{n-3} - a_{n-1} b_2}{b_1} Criterion: The number of roots with positive real parts (unstable poles) equals the number of sign changes in the first column of the Routh array. For absolute stability, there must be zero sign changes in the first column.

For a Third-Order System (s3+a2s2+a1s+a0=0s^3 + a_2 s^2 + a_1 s + a_0 = 0):

s31a1s2a2a0s1a2a1−a0a20s0a00\begin{array}{c|cc} s^3 & 1 & a_1 \\ s^2 & a_2 & a_0 \\ s^1 & \frac{a_2 a_1 - a_0}{a_2} & 0 \\ s^0 & a_0 & 0 \end{array}

The system is strictly stable if and only if:

a2>0,a1>0,a0>0,anda2a1>a0a_2 > 0, \quad a_1 > 0, \quad a_0 > 0, \quad \text{and} \quad a_2 a_1 > a_0

Step-by-Step Worked Sensor & Control Scenario

Problem Statement

An automated machine spindle drive requires precision torque measurement and closed-loop speed control.

Part A: Wheatstone Bridge Strain Gage Transducer

A solid steel shaft of diameter d=30.0 mmd = 30.0\,\text{mm} (0.030 m0.030\,\text{m}) transmits a steady working torque T=250.0 N⋅mT = 250.0\,\text{N}\cdot\text{m}. The shear modulus of elasticity of the steel is G=79.3 GPaG = 79.3\,\text{GPa}. Four identical metallic foil strain gages with nominal resistance R=350.0 ΩR = 350.0\,\Omega and Gage Factor GF=2.10GF = 2.10 are bonded to the shaft surface at ±45∘\pm 45^\circ helix angles to measure the principal tensile and compressive strains induced by torsion in a full Wheatstone bridge configuration. The bridge excitation voltage is Vs=10.0 VV_s = 10.0\,\text{V}.

Calculate:

  1. The maximum torsional surface shear stress τmax⁡\tau_{\max} and the resulting surface principal strain ϵ\epsilon.
  2. The differential output voltage VoutV_{\text{out}} produced by the full Wheatstone bridge in millivolts.

Part B: Closed-Loop Speed Control & Stability

The DC spindle drive and load have the plant transfer function:

G(s)=25s(s+5)G(s) = \frac{25}{s(s + 5)}

A proportional-derivative (PD) controller C(s)=Kp+KdsC(s) = K_p + K_d s is placed in the forward path with unity negative feedback (H(s)=1H(s) = 1).

Calculate:

  1. The closed-loop characteristic equation in terms of KpK_p and KdK_d.
  2. The required gains KpK_p and KdK_d to achieve a closed-loop natural frequency ωn=10.0 rad/s\omega_n = 10.0\,\text{rad/s} and a damping ratio ζ=0.707\zeta = 0.707.
  3. If an integral gain KiK_i is added to form a full PID controller (C(s)=Kds2+Kps+KisC(s) = \frac{K_d s^2 + K_p s + K_i}{s}), apply the Routh-Hurwitz criterion to determine the maximum upper bound on KiK_i for which the closed-loop system remains strictly stable.

Step-by-Step Solution

Step 1: Torsional Stress & Principal Strain

Compute the polar moment of inertia JJ of the solid circular shaft:

J=πd432=π(0.030 m)432=π×8.1×10−732=7.9522×10−8 m4J = \frac{\pi d^4}{32} = \frac{\pi (0.030\,\text{m})^4}{32} = \frac{\pi \times 8.1 \times 10^{-7}}{32} = 7.9522 \times 10^{-8}\,\text{m}^4

Compute maximum surface shear stress τmax⁡\tau_{\max}:

τmax⁡=T(d/2)J=16Tπd3=16×250.0π(0.030)3=4, ⁣0008.4823×10−5=4.7157×107 Pa=47.16 MPa\tau_{\max} = \frac{T (d/2)}{J} = \frac{16 T}{\pi d^3} = \frac{16 \times 250.0}{\pi (0.030)^3} = \frac{4,\!000}{8.4823 \times 10^{-5}} = 4.7157 \times 10^7\,\text{Pa} = 47.16\,\text{MPa}

Compute the surface engineering shear strain γ\gamma:

γ=τmax⁡G=47.157×106 Pa79.3×109 Pa=5.9467×10−4 rad\gamma = \frac{\tau_{\max}}{G} = \frac{47.157 \times 10^6\,\text{Pa}}{79.3 \times 10^9\,\text{Pa}} = 5.9467 \times 10^{-4}\,\text{rad}

In pure torsion, the principal normal stresses occur at ±45∘\pm 45^\circ to the shaft axis with magnitudes σ1=+τmax⁡\sigma_1 = +\tau_{\max} and σ2=−τmax⁡\sigma_2 = -\tau_{\max}. The corresponding principal normal strains are:

ϵ1=+γ2=+2.9734×10−4=+297.3 με\epsilon_1 = +\frac{\gamma}{2} = +2.9734 \times 10^{-4} = +297.3\,\mu\varepsilon ϵ2=−γ2=−2.9734×10−4=−297.3 με\epsilon_2 = -\frac{\gamma}{2} = -2.9734 \times 10^{-4} = -297.3\,\mu\varepsilon

Step 2: Full-Bridge Output Voltage

In a full-bridge torsion setup, gages 1 and 3 experience tensile principal strain ϵ1=ϵ3=+ϵ\epsilon_1 = \epsilon_3 = +\epsilon, and gages 2 and 4 experience compressive principal strain ϵ2=ϵ4=−ϵ\epsilon_2 = \epsilon_4 = -\epsilon:

VoutVs=GF4(ϵ1−ϵ2+ϵ3−ϵ4)=GF4[ϵ−(−ϵ)+ϵ−(−ϵ)]=GF4(4ϵ)=GF⋅ϵ\frac{V_{\text{out}}}{V_s} = \frac{GF}{4} (\epsilon_1 - \epsilon_2 + \epsilon_3 - \epsilon_4) = \frac{GF}{4} [\epsilon - (-\epsilon) + \epsilon - (-\epsilon)] = \frac{GF}{4} (4\epsilon) = GF \cdot \epsilon Vout=Vs⋅GF⋅ϵ=(10.0 V)×(2.10)×(2.9734×10−4)=6.244×10−3 V=6.24 mVV_{\text{out}} = V_s \cdot GF \cdot \epsilon = (10.0\,\text{V}) \times (2.10) \times (2.9734 \times 10^{-4}) = 6.244 \times 10^{-3}\,\text{V} = 6.24\,\text{mV}

The full-bridge circuit outputs 6.24 mV6.24\,\text{mV}, which can be directly fed into an instrumentation amplifier.

Step 3: PD Control Parameter Tuning

The closed-loop transfer function is:

T(s)=C(s)G(s)1+C(s)G(s)=(Kp+Kds)25s(s+5)1+(Kp+Kds)25s(s+5)=25Kds+25Kps(s+5)+25(Kds+Kp)T(s) = \frac{C(s) G(s)}{1 + C(s) G(s)} = \frac{(K_p + K_d s) \frac{25}{s(s + 5)}}{1 + (K_p + K_d s) \frac{25}{s(s + 5)}} = \frac{25 K_d s + 25 K_p}{s(s + 5) + 25(K_d s + K_p)}

The characteristic equation is:

s2+(5+25Kd)s+25Kp=0s^2 + (5 + 25 K_d) s + 25 K_p = 0

Compare to the standard second-order canonical form s2+2ζωns+ωn2=0s^2 + 2 \zeta \omega_n s + \omega_n^2 = 0:

  1. Match the constant term: ωn2=25Kp  ⟹  (10.0)2=25Kp  ⟹  100=25Kp  ⟹  Kp=4.0\omega_n^2 = 25 K_p \implies (10.0)^2 = 25 K_p \implies 100 = 25 K_p \implies K_p = 4.0
  2. Match the linear damping coefficient: 2ζωn=5+25Kd  ⟹  2(0.707)(10.0)=5+25Kd  ⟹  14.14=5+25Kd2 \zeta \omega_n = 5 + 25 K_d \implies 2 (0.707)(10.0) = 5 + 25 K_d \implies 14.14 = 5 + 25 K_d 25Kd=14.14−5=9.14  ⟹  Kd=9.1425=0.3656 s25 K_d = 14.14 - 5 = 9.14 \implies K_d = \frac{9.14}{25} = 0.3656\,\text{s}

Step 4: PID Stability Boundary via Routh-Hurwitz

With PID controller C(s)=Kds2+Kps+KisC(s) = \frac{K_d s^2 + K_p s + K_i}{s}, the loop transfer function is:

L(s)=C(s)G(s)=Kds2+Kps+Kis⋅25s(s+5)=25(Kds2+Kps+Ki)s2(s+5)L(s) = C(s) G(s) = \frac{K_d s^2 + K_p s + K_i}{s} \cdot \frac{25}{s(s + 5)} = \frac{25(K_d s^2 + K_p s + K_i)}{s^2(s + 5)}

The new third-order closed-loop characteristic equation is:

s2(s+5)+25(Kds2+Kps+Ki)=0  ⟹  s3+(5+25Kd)s2+25Kps+25Ki=0s^2(s + 5) + 25(K_d s^2 + K_p s + K_i) = 0 \implies s^3 + (5 + 25 K_d) s^2 + 25 K_p s + 25 K_i = 0

Substitute Kp=4.0K_p = 4.0 and Kd=0.3656K_d = 0.3656:

  • a2=5+25(0.3656)=5+9.14=14.14a_2 = 5 + 25(0.3656) = 5 + 9.14 = 14.14
  • a1=25Kp=25(4.0)=100.0a_1 = 25 K_p = 25(4.0) = 100.0
  • a0=25Kia_0 = 25 K_i

The characteristic equation is:

s3+14.14s2+100.0s+25Ki=0s^3 + 14.14 s^2 + 100.0 s + 25 K_i = 0

Applying the Routh-Hurwitz stability criterion for a 3rd-order system (a2a1>a0a_2 a_1 > a_0 with all ai>0a_i > 0):

25Ki>0  ⟹  Ki>025 K_i > 0 \implies K_i > 0 a2a1>a0  ⟹  (14.14)×(100.0)>25Ki  ⟹  1, ⁣414>25Kia_2 a_1 > a_0 \implies (14.14) \times (100.0) > 25 K_i \implies 1,\!414 > 25 K_i Ki<1, ⁣41425=56.56 s−1K_i < \frac{1,\!414}{25} = 56.56\,\text{s}^{-1}

Therefore, for the system to remain strictly stable, the integral gain must satisfy:

0<Ki<56.56 s−10 < K_i < 56.56\,\text{s}^{-1}

Exceeding Ki=56.56K_i = 56.56 will drive the closed-loop poles across the imaginary axis into the right-half plane, causing sustained oscillations and unstable runaway.

Test Your Knowledge

A structural test engineer mounts strain gages on a cantilever beam subjected to pure bending. To achieve maximum voltage sensitivity while automatically compensating for thermal expansion and cancelling any stray axial loads, which Wheatstone bridge configuration should be selected?

A

A quarter-bridge with one active gage on the tension surface and an unstrained reference dummy gage.

B

A half-bridge with two active gages mounted side-by-side along the upper tension surface.

C

A half-bridge with one active gage on the tension surface and one passive precision resistor on the compression surface.

D

A full-bridge with two active gages on the upper tension surface and two active gages on the lower compression surface.

Test Your Knowledge

An industrial servo positioning system exhibits unacceptable steady-state position error when holding against a constant gravitational disturbance torque, but its transient response is already lightly damped with significant overshoot. How should the controller parameters be adjusted to eliminate the steady-state error without degrading high-frequency noise rejection?

A

Increase proportional gain K_p substantially and remove the derivative gain K_d entirely.

B

Add integral gain K_i to remove the steady-state offset, and add derivative gain K_d to restore damping.

C

Increase derivative gain K_d toward infinity while removing the proportional gain K_p from the controller completely.

D

Switch to open-loop control and disable all of the feedback error terms to stop the overshoot.

Test Your Knowledge

A closed-loop motor velocity control system has the characteristic equation s^3 + 6 s^2 + 11 s + K = 0, where K is the adjustable loop gain. According to the Routh-Hurwitz stability criterion, what is the allowable range of gain K for which the closed-loop system remains strictly stable?

A

0 < K < 11

B

6 < K < 11

C

0 < K < 66

D

K > 66

Test Your Knowledge

A 10 kΩ10\text{ k}\Omega resistor and a 0.47 μF0.47\,\mu\text{F} capacitor form a first-order RC low-pass filter ahead of an analog-to-digital converter. What is the filter's cutoff frequency?

A

213 Hz213\text{ Hz}

B

3.39 Hz3.39\text{ Hz}

C

338.6 Hz338.6\text{ Hz}

D

33.9 Hz33.9\text{ Hz}

Sections you finish are checked off in the contents.