Calcolatore della Risposta all’Impulso del Sistema del Secondo Ordine
Guida Completa alla Risposta all’Impulso dei Sistemi del Secondo Ordine
La risposta all’impulso di un sistema dinamico del secondo ordine è un concetto fondamentale nell’ingegneria dei controlli automatici. Questo articolo fornisce una trattazione approfondita, dalle basi teoriche alle applicazioni pratiche, con particolare attenzione ai parametri che influenzano la risposta del sistema.
1. Fondamenti Teorici
Un sistema del secondo ordine è descritto dall’equazione differenziale:
a₂(d²y/dt²) + a₁(dy/dt) + a₀y = b₀u(t)
Dove u(t) è l’ingresso (nel caso dell’impulso, una delta di Dirac) e y(t) è l’uscita. La funzione di trasferimento standard è:
G(s) = ωₙ² / (s² + 2ζωₙs + ωₙ²)
Dove:
- ζ (zeta): rapporto di smorzamento (adimensionale)
- ωₙ (omega_n): frequenza naturale non smorzata [rad/s]
Classificazione dei Sistemi
In base al valore di ζ, i sistemi del secondo ordine si classificano in:
- Sottosmorzato (0 < ζ < 1): risposta oscillatoria
- Criticamente smorzato (ζ = 1): risposta più rapida senza oscillazioni
- Sovrasmorzato (ζ > 1): risposta lenta senza oscillazioni
- Non smorzato (ζ = 0): oscillazioni costanti
Parametri Chiave
I principali indicatori di prestazione sono:
- Tempo di picco (tₚ): tempo per raggiungere il primo massimo
- Sovraelongazione (%OS): scostamento percentuale dal valore finale
- Tempo di assestamento (tₛ): tempo per entrare nella banda del ±2% o ±5%
- Tempo di salita (tᵣ): tempo per passare dal 10% al 90% del valore finale
2. Risposta all’Impulso Unitario
Per un sistema del secondo ordine sottosmorzato (0 < ζ < 1), la risposta all'impulso unitario è data da:
h(t) = (ωₙ/√(1-ζ²)) * e^(-ζωₙt) * sin(ω_d t)
Dove ω_d = ωₙ√(1-ζ²) è la frequenza naturale smorzata.
Analisi dei Parametri
| Parametro | Formula | Descrizione |
|---|---|---|
| Tempo di picco (tₚ) | tₚ = π/(ωₙ√(1-ζ²)) | Tempo in cui si verifica il primo picco massimo |
| Valore di picco (Mₚ) | Mₚ = e^(-ζπ/√(1-ζ²)) | Valore massimo della risposta normalizzato |
| Sovraelongazione (%OS) | %OS = 100 * e^(-ζπ/√(1-ζ²)) | Scostamento percentuale oltre il valore finale |
| Tempo di assestamento (tₛ) | tₛ ≈ 4/(ζωₙ) (criterio 2%) | Tempo per entrare nella banda del ±2% |
3. Applicazioni Pratiche
I sistemi del secondo ordine modellano numerosi fenomeni fisici:
- Sistemi meccanici: massa-molla-smorzatore (sospensioni automobilistiche, edifici antisismici)
- Sistemi elettrici: circuiti RLC (filtri, oscillatori)
- Sistemi idraulici/pneumatici: attuatori, ammortizzatori
- Controlli automatici: regolatori PID, sistemi di guida
Esempio: Sospensione Automobilistica
In una sospensione, ζ ≈ 0.2-0.4 per comfort e ωₙ ≈ 1-2 Hz. Una sovraelongazione eccessiva (>20%) causa instabilità, mentre un tempo di assestamento troppo lungo (>2s) riduce la maneggevolezza.
Esempio: Filtro Elettronico
In un filtro passa-basso RLC, ζ = 1/2Q (dove Q è il fattore di qualità). Per ζ = 0.707 si ottiene una risposta di Butterworth (massima piattezza in banda passante).
4. Confronto tra Diverse Configurazioni
| Parametro | ζ = 0.1 (Basso Smorzamento) | ζ = 0.707 (Smorzamento Ottimale) | ζ = 1.0 (Critico) | ζ = 2.0 (Alto Smorzamento) |
|---|---|---|---|---|
| Sovraelongazione (%) | 72.9 | 4.3 | 0 | 0 |
| Tempo di picco (normalizzato) | 3.24 | 4.49 | — | — |
| Tempo di assestamento (normalizzato) | 75.7 | 5.6 | 4.0 | 8.0 |
| Tempo di salita (normalizzato) | 1.7 | 2.8 | 3.4 | 5.3 |
| Applicazioni tipiche | Oscillatori, orologi | Controlli generici, filtri | Sistemi critici (aerei) | Sistemi lenti (valvole) |
5. Metodologie di Progetto
La scelta di ζ e ωₙ dipende dagli obiettivi di progetto:
- Minimizzazione del tempo di assestamento:
- Aumentare ζ (fino a 0.7-0.8 per un compromesso)
- Aumentare ωₙ (ma attenzione al rumore e alla saturazione)
- Eliminazione della sovraelongazione:
- ζ ≥ 1 (sistema critico o sovrasmorzato)
- Utilizzare compensatori in retroazione (es. rete derivativa)
- Ottimizzazione della banda passante:
- ωₙ determina la velocità di risposta
- ζ influenza la risonanza (picco nella risposta in frequenza)
6. Strumenti di Analisi
Oltre all’analisi temporale, si utilizzano:
- Diagrammi di Bode: risposta in frequenza (modulo e fase)
- Luogo delle radici: analisi della stabilità al variare dei parametri
- Criterio di Nyquist: stabilità in sistemi retroazionati
- Simulazione numerica: MATLAB, Python (SciPy), o strumenti come il calcolatore sopra
7. Errori Comuni e Soluzioni
Problema: Sovraelongazione Eccessiva
Cause:
- ζ troppo basso (< 0.4)
- ωₙ troppo alta rispetto alla dinamica del sistema
Soluzioni:
- Aumentare lo smorzamento (es. aggiungere un ammortizzatore)
- Ridurre ωₙ (aumentare la massa o l’inerzia)
- Implementare un controllore derivativo
Problema: Tempo di Assestamento Lento
Cause:
- ζ troppo alto (> 0.9)
- ωₙ troppo bassa
Soluzioni:
- Ridurre ζ a 0.6-0.8
- Aumentare ωₙ (es. ridurre la massa)
- Usare un controllore proporzionale-integrale
8. Approfondimenti Matematici
La risposta all’impulso può essere derivata dalla trasformata di Laplace inversa della funzione di trasferimento. Per un sistema sottosmorzato:
L⁻¹{ωₙ² / (s² + 2ζωₙs + ωₙ²)} = (ωₙ/√(1-ζ²)) * e^(-ζωₙt) * sin(ω_d t)
Dove ω_d = ωₙ√(1-ζ²) è la frequenza delle oscillazioni smorzate. Il termine esponenziale e^(-ζωₙt) determina il decadimento dell’ampiezza, mentre il termine sinusoidale descrive l’oscillazione.
Per ζ ≥ 1, la risposta non è oscillatoria e assume forme diverse:
- ζ = 1 (critico): h(t) = ωₙ² t e^(-ωₙt)
- ζ > 1 (sovrasmorzato): h(t) = (ωₙ²/(2√(ζ²-1))) * (e^(-s₁t) – e^(-s₂t)), dove s₁,₂ = ωₙ(ζ ± √(ζ²-1))
9. Validazione Sperimentale
Per validare un modello del secondo ordine:
- Test all’impulso: applicare un ingresso a impulso e misurare l’uscita
- Stima dei parametri:
- Misurare la sovraelongazione per ricavare ζ
- Misurare il periodo delle oscillazioni per ricavare ω_d, quindi ωₙ
- Confronto con la risposta teorica: utilizzare strumenti come il calcolatore sopra
Esempio: in un sistema con sovraelongazione del 15% e periodo di oscillazione di 0.5s:
- ζ ≈ 0.5 (da %OS = 100*e^(-ζπ/√(1-ζ²)))
- ω_d ≈ 2π/0.5 = 12.57 rad/s
- ωₙ ≈ ω_d/√(1-ζ²) ≈ 14.43 rad/s
10. Risorse Esterne
Per approfondire: