Criteri di Routh-hurwitz per l'analisi delle vibrazioni meccaniche e delle oscillazioni

• Vibrazioni libere[]: Occur quando un sistema è disturbato e poi permesso di oscillare da solo, soggetto solo alle forze interne (ad esempio, un sistema di massa-spring rilasciato da uno spostamento iniziale).
• Vibrazioni forzate[]: Causate da una forza di variazione esterna (ad esempio, una massa rotante sbilanciata in un motore).
• Vibrazioni smorzate[[]: Involve la dissipazione dell'energia attraverso l'attrito, la resistenza dell'aria, o l'umidità del materiale, causando ampiezza di decomposizione nel tempo.
• Vibrazioni non amplificate[[]: Sistemi idealizzati dove non si perde energia – amplificazioni potenziali possono portare alla risonanza in applicazioni reali.

Le equazioni di governo della maggior parte dei problemi di vibrazione meccanica lineare sono equazioni differenziali ordinarie di secondo ordine (ODEs).Per un sistema monogrado-di-freedom (SDOF) la forma standard è:

m d2x/dt2 + c dx/dt + kx = F(t)]

c] è un coefficiente di smorzamento k è rigidità e F(t)]] è la forza esterna.

Criteri di Routh-Hurwitz: Fondazione matematica

Il criterio Routh-Hurwitz si applica all'equazione caratteristica di un sistema LTI, che può essere scritto come:

a0 sn + a1 sn−1 + a2 sn−2 + ... + an = 0]

s]] è la variabile Laplace (o l'igenvalore della matrice dello spazio-stato). Per il sistema essere stabile, tutte le radici devono risiedere nella metà sinistra del piano complesso (parti reali negative). Il criterio utilizza i coefficienti a0, a1, ..., un per determinare questo senza calcolare esplicitamente le radici.

Condizioni necessarie per la stabilità

Prima di applicare l'array completo di Routh, diverse condizioni necessarie (ma non sufficienti) devono contenere:

• Tutti i coefficienti a0...an[] devono essere presenti e positivi (senza coefficienti mancanti).
• Se un coefficiente è negativo o zero, il sistema è instabile o ha radici sull'asse immaginario (stabilità marginale).
• Un coefficiente zero di una potenza non-più alta indica un paio di radici equidistanti dall'origine, spesso portando a oscillazioni di ampiezza costante.

Costruendo il Routh Array

La serie Routh è costruita in modo ricorsivo dai coefficienti polinomiali:

Gli elementi della terza riga ([b1, b2, b3, ...) sono calcolati utilizzando determinanti di 2×2 blocchi dalle due righe precedenti:

b1 = (a1×a2 – a0×a3)/a1]

b2 = (a1×a4 – a0×a5)/a1]

Analogamente, le righe successive sono formate dalle due righe sopra, il processo continua fino a quando non rimane solo un singolo elemento nell'ultima riga.

Condizione di stabilità

Il criterio Routh-Hurwitz afferma:

• Tutte le radici hanno parti reali negative (sistema è stabile) se e solo se tutti gli elementi nella prima colonna del array Routh sono positivi (nonzero).
• Il numero di modifiche di segno nella prima colonna è uguale al numero di radici con parti reali positive (poli non tabella).
• Se la prima colonna contiene zero, si presentano casi speciali (vedi sotto).

Questo semplice controllo dei segni fornisce un potente e non-iterative test per la stabilità. Nell'analisi delle vibrazioni meccaniche, viene applicato di routine per rilevare modalità di oscillazione instabili, quelle che crescerebbero senza limiti con il tempo.

Applicare i criteri ai sistemi di vibrazione meccanica

Sistemi di libertà (SDOF) a singolo livello

Per un tipico oscillatore SDOF umido, l'equazione caratteristica è:

m s2 + c s + k = 0]

Qui, a0 = m[], a1 = c], a2 = k. La serie di routine dell'ordine 2 è banale:

• S2 di fila: m, k
• Voce principale: c, 0
• S0: (c×k – m×0)/c = k

Prima colonna: m, c, k. Poiché tutti i coefficienti sono positivi per un sistema passivo, la condizione per la stabilità è che tutti e tre sono > 0—che è automaticamente vero. Questo corrisponde all'intuizione fisica: massa positiva, smorzamento e rigidità producono un'oscillazione stabile e decaduta. Se l'umidità fosse negativa (]]c < 0), il sistema di controllo potrebbe essere instabile, un comportamento passivo.

Sistemi multi-digitale (MDOF)

I sistemi meccanici reali, come sospensioni dei veicoli, rotori delle turbine o strutture edili, hanno molti gradi di libertà. Il loro caratteristico polinomio può essere di ordine 4, 6 o superiore. Il Routh-Hurwitz Criterion diventa particolarmente prezioso qui: calcolare gli autovalori di una grande matrice è computazionalmente intensiva, ma la costruzione della serie Routh è semplice e può rivelare margini di instabilità rapidamente.

Considera un sistema ammortizzatore a vibrazioni a quattro gradi di libertà con il caratteristico polinomio:

a0 s4 + a1 s3 + a2 s2 + a3 s + a4 = 0]

La gamma Routh per questa quartica prevede l'elaborazione b1, b2, c1, ecc I segni di primo livello indicano stabilità. In pratica, gli ingegneri spesso automatizzano la costruzione di array nel software, ma la comprensione del processo manuale costruisce l'intuizione.

Casi speciali nell'analisi di Routh-Hurwitz

Due casi speciali spesso si presentano in problemi di vibrazione meccanica:

1. Zero nella prima colonna

Se un elemento nella prima colonna è zero mentre gli elementi rimanenti della riga non sono zero, il sistema è marginalmente stabile o instabile. Uno zero nella prima colonna impedisce di completare direttamente l'array. Il rimedio standard è quello di sostituire lo zero con un piccolo numero positivo ε (epsilon) e quindi esaminare il limite come ε→0. Il segno cambia nella colonna risultante indica il numero di radici instabili.

Ad esempio, un polinomio come s3 + 2s2 + s + 2 = 0 produce uno zero nella prima colonna della riga s1. Utilizzando il metodo ε si rivelano due cambiamenti di segno, che significa due radici nella metà piano destra—non installabile.

2. Intera fila di zero

Quando un'intera fila nella matrice consiste di zero, il polinomio contiene radici simmetriche sull'origine (ad esempio, coppie di radici puramente immaginarie o radici reali con segni opposti). Ciò indica stabilità marginale (oscillazioni dismesse a ampiezza costante) o instabilità. La procedura è quella di formare un polinomio ausiliario dalla fila sopra la riga zero, differenziarla e sostituire la riga zero con i coefficienti del derivato.

Questo scenario è comune in sistemi meccanici non smorzati, dove il caratteristico polinomio comprende coppie di radici immaginarie che rappresentano frequenze naturali senza smorzamento. Ad esempio, un sistema a due maglie non smorzato potrebbe produrre un polinomio s4 + 5s2 + 4 = 0]]. La fila zero a s2 richiede la costruzione delle radici immaginario di dmar.

Esempio di lavoro: stabilità di un braccio robotico flessibile

Applicare il criterio Routh-Hurwitz ad un modello semplificato di un braccio robotico flessibile, noto per esporre vibrazioni instabili in determinate condizioni. L’equazione caratteristica derivata dalla dinamica flessibile è:

s4 + 3s3 + 5s2 + 9s + 4 = 0]

Step 1: Condizioni necessarie. Tutti i coefficienti sono positivi e presenti → condizioni necessarie soddisfatte (ma non sufficienti).

Step 2: Construct the Routh array.

• 1, 5, 4
• S3: 3, 9, 0
• S2 riga: Computo b1 = (3×5 – 1×9)/3 = (15 – 9)/3 = 2; b2 = (3×4 – 1×0)/3 = 12/3 = 4
• S1: c1 = (2×9 – 3×4)/2 = (18 – 12)/2 = 3; c2 = (2×0 – 3×0)/2 = 0
• S0: d1 = (3×4 – 2×0)/3 = 12/3 = 4

I valori della prima colonna sono: 1, 3, 2, 3, 4 — tutti positivi. Nessun cambiamento di segno. Pertanto, il sistema è stabile, il che significa che il braccio flessibile non esporrà oscillazioni crescenti. Se il coefficiente di smorzamento (relativo al termine s3) è stato ridotto, potrebbero apparire cambiamenti di segno, indicando l'inizio dell'instabilità del flusso.

Questo esempio illustra come i Criteri Routh-Hurwitz possono essere utilizzati per guidare le scelte dei parametri di progettazione per garantire una robusta stabilità dei sistemi meccanici.

Applicazioni pratiche in Ingegneria

Il criterio Routh-Hurwitz è ampiamente impiegato nella progettazione e nell'analisi di sistemi meccanici e aerospaziali in cui il controllo dell'oscillazione è fondamentale.

• Progettazione sospensioni Veicoli[[[]: Ottimizzazione dei tassi di smorzamento e di molla per fornire un giro confortevole evitando l'instabilità (ad esempio, shimmy nei sistemi di sterzo).
• Le dinamiche di rotore[]: Assicurare che i macchinari rotanti (turbine, compressori) funzionino sotto velocità critiche dove si può verificare l'instabilità.
• Strutture aerospaziali[[[]]: Analizzando i margini di fluttuazione delle ali e delle superfici di controllo; il criterio è uno strumento standard nell'aeroelasticità per i controlli di stabilità rapidi.
• Sistemi di isolamento sismico[[]: Progettazione isolatori di base che impediscono la risonanza degli edifici durante i terremoti, assicurando che le radici caratteristiche del sistema rimangano nel mezzo piano sinistro.
• Controllo delle vibrazioni attivo[[]: Nei sistemi di controllo dei feedback, il polinomio caratteristico a ciclo chiuso deve soddisfare le condizioni di Routh-Hurwitz per una soppressione stabile delle vibrazioni.

Il criterio serve anche come precursore di metodi di analisi della stabilità più avanzati come il criterio di Nyquist o la trama del locus della radice. La sua natura algebrica lo rende ideale per studi parametrici: variando un coefficiente (rappresentando un parametro fisico come l'umidità), gli ingegneri possono determinare la soglia di instabilità.

Limitazioni e Alternative

Mentre potente, il criterio Routh-Hurwitz ha limitazioni. Si applica solo a sistemi lineari e time-invariant. Molti moderni sistemi meccanici includono non linearità (ad esempio, backlash attrito, isteresi materiale) o parametri di tempo-varying (ad esempio, rigidità variabile nelle strutture di morfing dettagliate). In tali casi, l'analisi di stabilità non lineare (metodo di Lyapunov) o la simulazione è necessaria.

I metodi alternativi per l'analisi della stabilità includono:

• Root locus tecnica[[]: Mostra come le radici si muovono come un parametro cambia.
• Trame di bode e Nyquist[[]: Fornire la percezione del dominio di frequenza ai margini di stabilità (gain margin, margine di fase).
• Metodo determinante della matrice di Hurwitz[[]: Equivalente al array Routh ma utilizza determinanti delle sottomatrice Hurwitz.
• Computazione numerica dell'autovalore[: Computes direttamente i poli utilizzando software come MATLAB o la libreria di controllo di Python.

Nonostante queste alternative, il Criterion Routh-Hurwitz rimane uno strumento di prima linea per la sua semplicità computazionale e la chiara comprensione che offre nel rapporto tra coefficienti polinomiali e stabilità.

Conclusioni

Comprendere e applicare il criterio Routh-Hurwitz è essenziale per qualsiasi ingegnere impegnato nell'analisi delle vibrazioni meccaniche e delle oscillazioni. Fornisce un percorso robusto e algebrico per determinare la stabilità del sistema senza risolvere equazioni caratteristiche di alto grado.

Dalle semplici oscillazioni SDOF alle complesse strutture multi-grado-di-libero, il criterio aiuta a impostare i parametri di progettazione per smorzamento, rigidità e distribuzione di massa. È ampiamente insegnato in curricula di ingegneria ed è implementato nella maggior parte degli strumenti di ingegneria assistiti dal computer. Tuttavia, la comprensione del metodo manuale approfondisce l'intuizione di come i parametri di sistema influenzano la stabilità - conoscenza che rimane inestimabile nella progettazione di sistemi meccanici sicuri e affidabili.

Altri riferimenti e letture:

• criterio di stabilità di Routh–Hurwitz – Wikipedia
• Criteri di Routh-Hurwitz – Kiel University Control Engineering
• MIT OpenCourseWare: Engineering Dynamics (contenuti di analisi della vibrazione)