Lesson: Engineering Mathematics I & II
Course: AUTO_1 — Foundation Automotive Technician Program (Beginners in Resource-Constrained African Contexts)
Learning objectives
After studying this topic, the learner will be able to:
- Recognize ordinary differential equations (ODEs) commonly used in automotive problems (first- and second-order linear ODEs).
- Solve simple ODEs by direct methods (separation of variables, integrating factor) and interpret solutions physically.
- Apply Laplace transforms to solve initial-value problems and obtain system transfer functions.
- Use ODE models and Laplace-domain methods to analyze vehicle dynamics (suspension), thermal systems (engine cooling), and simple control responses (actuator/throttle response).
- Perform low-cost experiments to determine model parameters (time constant, damping ratio) using locally available materials and smartphones or simple sensors.
Prerequisites
- Basic algebra, differentiation and integration.
- Basic complex-number concept (real and imaginary parts).
- Familiarity with physical concepts: mass, force, stiffness, damping, heat capacity, convection.
1. Overview — Why differential equations and Laplace transforms?
Many physical systems in vehicles are dynamic: variables change over time (speed, suspension displacement, temperature). Differential equations express the balance of rates (time derivatives) and forces (or thermal flows). Laplace transforms convert time-domain differential equations into algebraic equations in the complex frequency domain (s-domain), which simplifies solving linear initial-value problems and deriving transfer functions used in control and system analysis.
2. Types of ordinary differential equations (ODEs) encountered in automotive contexts
- First-order linear ODE (typical for simple thermal systems, electrical RC circuits, simple actuator dynamics):
dy/dt + a y = f(t) - Second-order linear ODE (typical for mass–spring–damper models such as vehicle suspension, some rotational dynamics):
m d^2y/dt^2 + c dy/dt + k y = F(t) - Linear constant-coefficient ODEs are the primary focus here. Nonlinear models exist (friction, saturation), but are often linearised around an operating point for basic analysis.
3. First-order ODEs — methods and physical example
General form: dy/dt + p(t) y = q(t). Two simple solution methods:
- Separation of variables: if dy/dt = g(y) h(t) you can separate and integrate.
- Integrating factor: for linear first-order ODEs with p(t) known, multiply by μ(t)=exp(∫p(t) dt).
Physical example — Newton’s law of cooling (engine or component cooling):
- Model: dT/dt = −k (T − Ta)
where T(t) is component temperature, Ta is ambient temperature, k > 0 is a constant (s−1). - Rearranged: dT/dt + k T = k Ta
- Solution (direct): T(t) = Ta + (T(0) − Ta) e^(−k t)
- Interpretation: time constant τ = 1/k. After time ≈ 3τ the temperature is close to Ta (≈95%).
Practical note (resource-constrained): measure cooling with a thermometer or thermocouple; record T(t) after engine is turned off and fit an exponential to estimate k (or τ). A smartphone-connected thermocouple or simple digital thermometer works.
4. Second-order ODEs — mass–spring–damper model
Standard form: m y” + c y’ + k y = F(t)
- y(t): displacement
- m: mass
- c: damping coefficient
- k: stiffness
- F(t): external force (road input, impulse, step)
Key parameters:
- Natural (undamped) frequency: ωn = sqrt(k/m) (rad/s)
- Damping ratio: ζ = c / (2 sqrt(m k))
- Characteristic equation (homogeneous): m s^2 + c s + k = 0 → roots s1, s2 (poles)
Time-domain behaviour:
- ζ > 1: overdamped (no oscillation, slow)
- ζ = 1: critically damped (fastest without oscillation)
- 0 < ζ < 1: underdamped (oscillatory decay)
- ζ = 0: undamped (sustained oscillation if no friction)
Applied example — quarter-car suspension (simplified):
- Typical approximate values: m = 250 kg (quarter of vehicle mass over one wheel), k ≈ 15,000 N/m, c chosen to give ζ ≈ 0.3–0.5
- ωn = sqrt(15000 / 250) = sqrt(60) ≈ 7.746 rad/s → natural frequency ≈ 1.233 Hz
- With c chosen for ζ = 0.3 → c = 2 ζ sqrt(m k) ≈ 2 × 0.3 × sqrt(250 × 15000) ≈ 2 × 0.3 × 1936.5 ≈ 1161.9 N·s/m
If the vehicle wheel hits a bump approximated as a step/pulse in force or displacement, the mass–spring–damper response may be oscillatory (ride comfort) depending on ζ.
Practical low-cost experiment: assemble a mass (weights) on a spring (locally made or repurposed spring—e.g., suspension spring from a motorcycle or used car), add dashpot damping using oil viscous damping (a piston in oil) or bicycle tire filled with oil-ish medium. Use a smartphone accelerometer app to record step response and estimate ωn and ζ from oscillation frequency and envelope decay.
5. Laplace transform — basic definitions and properties
Definition (for causal functions f(t), t ≥ 0):
L{f(t)} = F(s) = ∫_0^∞ e^(−s t) f(t) dt, where s is a complex variable (s = σ + jω).
Common transforms:
- L{1} = 1/s
- L{e^(−a t)} = 1/(s + a)
- L{sin(ω t)} = ω/(s^2 + ω^2)
- L{cos(ω t)} = s/(s^2 + ω^2)
- L{t^n} = n! / s^(n+1)
- L{u(t)} (unit step) = 1/s
- L{δ(t)} (Dirac impulse) = 1
Useful properties:
- Linearity: L{a f(t) + b g(t)} = a F(s) + b G(s)
- Differentiation in time: L{f'(t)} = s F(s) − f(0)
L{f”(t)} = s^2 F(s) − s f(0) − f'(0) - Time shift: L{f(t − t0) u(t − t0)} = e^(−s t0) F(s)
- Convolution: L{f * g} = F(s) G(s) (useful for input–output relations)
Interpretation: Laplace transform converts ODEs with initial conditions into algebraic equations in s that include initial-condition terms. Solving for F(s) and taking inverse Laplace gives f(t).
6. Transfer functions and system response
- Transfer function G(s) of a linear time-invariant (LTI) system (zero initial conditions) is ratio of output to input in s-domain:
G(s) = Y(s) / U(s) - For mass–spring–damper with input force F(t) and output displacement y(t) and zero initial conditions:
G(s) = Y(s) / F(s) = 1 / (m s^2 + c s + k) - Poles: roots of denominator. Zeros: roots of numerator.
- Stability: system stable if all poles have negative real parts (Re(s) < 0).
Step response: input U(s) = (A / s) for a step of amplitude A. Multiply by G(s) and invert to obtain y(t).
7. Worked examples
Example 1 — First-order thermal system (Newton’s cooling), direct and Laplace solution
- Model: dT/dt + k T = k Ta, initial T(0) = T0
Direct solution:
T(t) = Ta + (T0 − Ta) e^(−k t)
Laplace method:
Take Laplace:
s T(s) − T0 + k T(s) = k Ta / s
→ T(s) (s + k) = T0 + k Ta / s
→ T(s) = T0 / (s + k) + (k Ta) / (s (s + k))
Inverse-Laplace yields same T(t): Ta + (T0 − Ta) e^(−k t).
Practical numeric example:
- Suppose an engine component cools from 90 °C to ambient Ta = 25 °C; measured time constant τ = 120 s (k = 1/120 ≈ 0.00833 s−1). Then:
T(t) = 25 + 65 e^(−t/120)
Experimental application: use thermometer data to estimate τ: plot ln(T(t) − Ta) vs t, slope = −1/τ. Smartphones with temperature probes or cheap thermocouples yield adequate data.
Example 2 — Quarter-car (second-order) step force response using Laplace
- Model: m y” + c y’ + k y = F0 u(t) (step force F0)
- Take Laplace (zero initial conditions):
(m s^2 + c s + k) Y(s) = F0 / s
Y(s) = (F0 / s) * 1 / (m s^2 + c s + k)
Transfer function: G(s) = 1 / (m s^2 + c s + k)
To invert, factor denominator or express in standard form:
Divide numerator and denominator by m:
G(s) = (1/m) / (s^2 + 2 ζ ωn s + ωn^2)
Step response y(t) with F0 = constant:
- For underdamped (0 < ζ < 1):
y(t) = (F0 / k) [1 − (1 / sqrt(1 − ζ^2)) e^(−ζ ωn t) sin(ωd t + φ)]
where ωd = ωn sqrt(1 − ζ^2) and φ = arctan( sqrt(1 − ζ^2) / ζ )
Numeric illustration:
- m = 250 kg, k = 15000 N/m → ωn ≈ 7.746 rad/s
- ζ = 0.3 → c ≈ 1162 N·s/m
- Step force F0 = 2000 N (example bump impulse approximated as a sustained force)
- Static displacement = F0 / k = 2000 / 15000 ≈ 0.1333 m
The dynamic response will oscillate and settle to 0.1333 m with exponential decay determined by ζ and ωn.
Practical measurement: use smartphone accelerometer or low-cost displacement sensor (string potentiometer made with a wiper and spring) to record response to a known mass or bump; fit data to standard second-order model to estimate m, c, k.
Example 3 — First-order actuator/throttle response using Laplace
- Simple actuator model (electrical + mechanical inertia simplified) approximated as first-order:
τ dy/dt + y = K u(t)
where y is actuator position, u is input command, τ is time constant, K steady-state gain. - Transfer function: G(s) = Y(s)/U(s) = K / (τ s + 1)
- Step input u(t) = U0 u(t) → Y(s) = (K U0) / (s (τ s + 1)). Inverse gives:
y(t) = K U0 [1 − e^(−t/τ)] - Measurement approach: command step input with microcontroller or manual input; record output and estimate τ as time to reach 63% of final value.
Low-cost setup: a DC motor with gear and potentiometer can form an actuator; drive it with a simple PWM signal from a small microcontroller (Arduino/clone) and measure position with the potentiometer. Use smartphone to log data via serial if necessary.
8. Practical considerations for resource-constrained contexts
- Linear approximations: many vehicle components are nonlinear; however, small deviations around operating points are often approximated by linear ODEs for useful first estimates.
- Low-cost sensors: thermocouples, digital thermometers, smartphone accelerometers, and potentiometers are accessible and useful for parameter identification.
- Data acquisition: use mobile phones, low-cost microcontrollers (Arduino, ESP32), or simple logging circuits.
- Use locally available parts for experiments: bicycle springs, reclaimed car springs, weights for masses, viscous damping using household liquids for qualitative experiments.
- Emphasize safety when working with vehicle components, hot engines, and electrical systems.
9. Simple laboratory exercises (competency-based)
-
Cooling experiment (first-order)
- Heat a metal block (or warm engine oil in container) and record cooling temperature vs time using thermometer.
- Fit exponential curve and extract time constant τ.
- Show calculations and physical meaning (heat capacity and heat transfer coefficient relation).
-
Spring–mass–damper experiment (second-order)
- Build spring–mass system using local spring and known mass.
- Provide a step input (displace and release, or apply a steady force) and record displacement or acceleration using smartphone.
- Estimate ωn and ζ from oscillation frequency and decay envelope.
-
Actuator step response (first-order)
- Use a small DC motor and potentiometer to measure position response to step input voltage.
- Estimate time constant τ and steady-state gain K.
Each exercise should include:
- Data collection and plotting.
- Parameter estimation and comparison with model predictions.
- Discussion of discrepancies and sources of error.
10. Summary — key takeaways
- Differential equations model dynamic behaviour: first-order for many thermal and simple actuator responses; second-order for mass–spring–damper systems.
- Laplace transforms are a practical tool to solve linear ODEs with initial conditions and to derive transfer functions used in control and system analysis.
- Important system parameters: time constant (τ), natural frequency (ωn), damping ratio (ζ), and static gain (K).
- Low-cost experiments and smartphone-based measurements can validate models and help learners develop practical skills in resource-constrained settings.
References and further reading (recommended for learners)
- Elementary differential equations and their applications — standard undergraduate texts.
- Introductory control systems texts (for transfer functions and step responses).
- Practical guides on using smartphone sensors for experimental mechanics and thermometry.
End of topic.