Física universitaria — primer año · Bachelor Year 1
7Regímenes transitorios: primer y segundo orden
Basta empujar hacia abajo la esquina de un coche aparcado y soltar. Un coche en buen estado sube una vez y se detiene; un coche con los amortiguadores gastados cabecea dos o tres veces antes de calmarse. En un laboratorio de electrónica aparecen esos mismos dos comportamientos en un osciloscopio cuando un condensador se descarga a través de una bobina: un retorno suave, o una oscilación que se extingue. Ambos son el transitorio entre un régimen permanente y otro, y ambos obedecen la misma pequeña familia de ecuaciones diferenciales lineales — de primer orden para un único almacén de energía, de segundo orden para dos que intercambian. Este capítulo las resuelve de una vez por todas, nombra sus parámetros (constante de tiempo, frecuencia propia, factor de calidad) y muestra cómo los mismos números describen un circuito y una suspensión.
Una bobina, un condensador y un resistor sobre una placa de pruebas, y en el osciloscopio la oscilación amortiguada que sigue a un escalón: el régimen seudoperiódico de este capítulo.
7.1 Sistemas de primer orden
Definición 7.1(Sistema lineal de primer orden)
Una magnitud x(t) obedece una ecuación lineal de primer orden con coeficientes constantes cuando
τdtdx+x=x∞,
con τ>0 la constante de tiempo y x∞ una constante (el valor impuesto por la fuente). Su régimen permanente es x=x∞; el régimen transitorio es la aproximación a él desde el valor inicial x(0)=x0.
Teorema 7.2(Solución)
La única solución con x(0)=x0 es
x(t)=x∞+(x0−x∞)e−t/τ.
La distancia al régimen permanente se reduce en el factor e−1≈0.37 cada τ: se ha recorrido el 63% del camino en τ, el 95% en 3τ y el 99% en 4.6τ. La tangente en el origen alcanza x∞ en t=τ.
Demostración.y=x−x∞ obedece τy′+y=0, cuyas soluciones son y=Ke−t/τ (el curso de matemáticas sobre ecuaciones diferenciales lineales demuestra que no hay otras); K=x0−x∞. La pendiente en 0 vale −(x0−x∞)/τ, que cerraría la distancia en un tiempo τ. ∎
La respuesta de primer orden a un escalón desde x0=0: la tangente inicial corta la asíntota en t=τ; el 63% del camino en τ y el 95% en 3τ. Todo transitorio de primer orden es esta curva, dilatada por τ y desplazada por x0 y x∞.
Proposición 7.3(Circuitos RC y RL)
Un condensador C que se carga a través de un resistorR desde una fuente ideal E obedece RCduC/dt+uC=E: τ=RC, uC,∞=E. Una bobina L en serie con R sometida a E obedece (L/R)di/dt+i=E/R: τ=L/R, i∞=E/R. Retirada la fuente (E→0), esas mismas ecuaciones describen la descarga hacia cero.
La tensión en bornes de un condensador y la corriente por una bobina son funciones continuas del tiempo: en un instante de conmutación t0, uC(t0+)=uC(t0−) e iL(t0+)=iL(t0−). Las corrientes en los resistores y las tensiones en las bobinas sí pueden saltar.
Demostración. Un salto de uC exigiría i=CduC/dt infinita, y un salto de iL exigiría u=Ldi/dt infinita; ambas cosas son imposibles en un circuito de tensiones y corrientes finitas (dicho de otro modo, las energías almacenadas 21CuC2 y 21Li2 no pueden cambiar instantáneamente sin una potencia infinita). ∎
Método 7.5(Resolver un transitorio de primer orden)
Escríbanse las leyes de mallas y de nudos después de la conmutación, junto con las leyes de los componentes; redúzcase a una única ecuación en uC o en iL y pásese a la forma τx′+x=x∞; léanse τ y x∞ (este último es también el régimen permanente que se obtiene sustituyendo C por un circuito abierto y L por un hilo).
Hállese x0 por la continuidad de uC o de iL en el instante de conmutación, usando el circuito anterior a ella.
Escríbase x=x∞+(x0−x∞)e−t/τ; dedúzcanse las demás magnitudes por derivación y por la ley de las mallas.
Cuando el condensador ve una red de resistores y fuentes, sustitúyase antes esa red por su equivalente de Thévenin (ETh,RTh) (Capítulo 6): τ=RThC, u∞=ETh.
Proposición 7.6(Balance energético de la carga RC)
Al cargar un condensador de 0 a E a través de un resistorR cualquiera, la fuente aporta CE2, el condensador almacena 21CE2 y el resistor disipa el otro 21CE2 — sea cual sea R.
Demostración.i=(E/R)e−t/τ; la fuente aporta ∫0∞Eidt=E⋅(E/R)τ=CE2; el resistor toma ∫0∞Ri2dt=(E2/R)∫0∞e−2t/τdt=(E2/R)(τ/2)=21CE2; la diferencia es la energía almacenada 21CE2. Una R menor disipa la misma energía más deprisa. ∎
Ejemplo 7.7(La bobina de un relé)
L=20mH, R=5.0Ω, E=10V: τ=4.0ms, i∞=2.0A; los contactos se cierran cuando i llega a 1.5A, es decir en t=−τln(1−0.75)=5.5ms. Al abrir el circuito, i debe caer desde 2A; un diodo de rueda libre en paralelo con la bobina le ofrece un camino a través de R solo, con τ=4ms, en vez de una chispa.
7.2 Sistemas de segundo orden
Definición 7.8(Forma canónica de segundo orden)
Una magnitud x(t) es un sistema lineal de segundo orden cuando
Un condensador C, una bobina L y un resistorR en serie sometidos a una fuente E: la tensión del condensador obedece la forma canónica con
ω0=LC1,Q=R1CL,uC,∞=E.
Demostración.Ley de las mallas: E=Ldi/dt+Ri+uC con i=CduC/dt: LCuC′′+RCuC′+uC=E; divídase por LC e identifíquense ω02=1/LC y ω0/Q=R/L. ∎
Las dos caras del sistema de segundo orden: el circuito RLC serie (ω0=1/LC, Q=L/C/R) y el conjunto masa–resorte–amortiguador (ω0=k/m, Q=km/α). La misma ecuación y los mismos tres regímenes.
Proposición 7.10(El oscilador mecánico)
Una masa m sobre un resorte de rigidez k con un amortiguador viscoso de coeficiente α (fuerza −αx˙), separada x de la posición de equilibrio, obedece mx¨+αx˙+kx=0: la forma canónica con
ω0=mk,Q=αkm.
La correspondencia x↔q, x˙↔i, m↔L, α↔R, k↔1/C lleva todo resultado de un sistema al otro.
Demostración. La segunda ley de Newton (Capítulo 12) con la fuerza del resorte −kx y la fuerza de amortiguamiento; divídase por m. ∎
Teorema 7.11(Los tres regímenes)
Las soluciones de la ecuación homogénea x′′+(ω0/Q)x′+ω02x=0 las gobiernan las raíces de r2+(ω0/Q)r+ω02=0, de discriminante Δ=ω02(1/Q2−4):
Q<21, aperiódico: dos raíces reales negativas r±=−2Qω0(1∓1−4Q2), x=Aer+t+Ber−t, un retorno sin oscilación;
Q=21, crítico: una raíz doble r=−ω0, x=(A+Bt)e−ω0t, el retorno más rápido sin sobreoscilación;
Q>21, seudoperiódico: raíces complejas r=−1/τ±iω con
τ=ω02Q,ω=ω01−4Q21,x=e−t/τ(Acosωt+Bsinωt),
una oscilación amortiguada de seudoperiodoT=2π/ω dentro de la envolvente ±A2+B2e−t/τ.
La solución completa con segundo miembro constante es x∞ más la solución homogénea; A y B se obtienen de x(0) y x′(0).
Demostración. Ecuación característica de una ecuación lineal con coeficientes constantes (curso de matemáticas, ecuaciones lineales de segundo orden): las raíces son r=−ω0/2Q±21Δ. Para Q>21, Δ=iω04−1/Q2, lo que da la parte real −ω0/2Q=−1/τ y la parte imaginaria ±ω01−1/4Q2; las soluciones reales son las combinaciones de e−t/τcosωt y e−t/τsinωt. ∎
Respuesta a un escalón, partiendo del reposo, de un sistema de segundo orden, para tres factores de calidad. Q pequeño: una subida lenta; Q=21: el retorno más rápido sin sobreoscilación; Q grande: una oscilación que dura alrededor de Q periodos.
Proposición 7.12(Leer Q en un registro seudoperiódico)
da Q; el número de oscilaciones visibles antes de que la amplitud baje del 5% es del orden de Q. Para Q≫1, ω≈ω0 y T≈2π/ω0.
Demostración. La envolvente e−t/τ se multiplica por e−T/τ en cada periodo; T/τ=(2π/ω)(ω0/2Q). La envolvente llega a e−3≈5% en 3τ=6Q/ω0≈QT, es decir al cabo de unos Q periodos. ∎
Decaimiento libre seudoperiódico, x=x0e−t/τcosωt: los máximos sucesivos se reducen en el factor constante e−T/τ, y ln(xn/xn+1)=T/τ≈π/Q mide el factor de calidad a partir del registro.
Ejemplo 7.13(Una oscilación RLC)
L=10mH, C=100nF, R=100Ω: ω0=1/10−9=3.16×104rad/s (f0=5.03kHz), Q=L/C/R=316/100=3.16: seudoperiódico, τ=2Q/ω0=0.20ms, ω=0.987ω0, T=0.20ms, tres oscilaciones visibles. La resistencia crítica vale Rc=2L/C=632Ω; por encima de ella la descarga es aperiódica.
Método 7.14(Resolver un transitorio de segundo orden)
Escríbase la ecuación; pásese a la forma canónica; léanse ω0, Q y x∞.
Decídase el régimen a partir de Q y escríbase la solución general como x∞ más la solución homogénea.
Hállense x(0) y x′(0) por la continuidad de uC y de iL (para el RLC: uC(0+)=uC(0−) y uC′(0+)=i(0−)/C); despéjense las dos constantes.
Ejemplo 7.15(Condiciones iniciales)
El condensador del ejemplo anterior se carga a U0 y, en t=0, se cierra sobre L y R (sin fuente): x∞=0, uC(0)=U0, uC′(0)=i(0)/C=0, pues la corriente de la bobina era nula. Por tanto A=U0 y −A/τ+Bω=0, B=U0/(ωτ)=U0/4Q2−1: uC=U0e−t/τ[cosωt+sin(ωt)/4Q2−1] — para Q≫1, simplemente U0e−t/τcosω0t, y la corriente i=−CduC/dt tiene un máximo cercano a U0C/L=U0/(ω0L).
Observación 7.16(Energía)
La energía almacenada 21CuC2+21Li2 (o 21kx2+21mx˙2) decrece a razón de Ri2 (αx˙2): derívese y úsese la ecuación. En cada seudoperiodo la energía se reduce en el factor e−2T/τ; con Q=3 se pierde alrededor del 88% por oscilación. Un Q alto significa una fuga lenta comparada con la oscilación — un buen reloj, una resonancia aguda (Capítulo 8) — y un Q bajo, un retorno rápido: lo que conviene a una suspensión o a la aguja de un aparato de medida.
7.3 Ejercicios
Ejercicio 7.1★
Un circuito RC: R=10kΩ, C=47nF, cargado desde E. Calcúlense τ, uC(τ)/E y el tiempo para llegar al 99% de E.
Un circuito RL: L=20mH, R=5.0Ω, E=10V. Calcúlense τ, la corriente final, la corriente en t=τ y la tensión en bornes de la bobina justo después de cerrar el interruptor.
Solución
Solución de Ejercicio 7.2.
τ=L/R=4.0ms; i∞=E/R=2.0A; i(τ)=0.63×2.0=1.26A; en 0+, i=0 (continuidad), luego uL=E−Ri=10V.
Ejercicio 7.3★
Un RLC serie: L=10mH, C=100nF, R=100Ω. Calcúlense ω0, f0 y Q; nómbrese el régimen; dense τ y el seudoperiodo.
Con los mismos L y C, ¿qué resistencia da el régimen crítico? ¿Cuánto vale entonces Q? ¿Qué ocurre para R=2kΩ?
Solución
Solución de Ejercicio 7.4.
Rc=2L/C=632Ω, Q=21. A 2kΩ, Q=0.16: aperiódico, una subida lenta gobernada por la raíz pequeña ∣r+∣≈ω0Q=5×103s−1 (escala de tiempo 0.2ms, pero sin oscilación).
Ejercicio 7.5★★
Un condensador cargado C=1.0µF se descarga a través de una R desconocida; la tensión se reduce a la mitad cada 3.5ms. Muéstrese que el tiempo de semirreducción vale τln2 y hállese R.
Solución
Solución de Ejercicio 7.5.
u=U0e−t/τ se reduce a la mitad cuando e−t/τ=21: t1/2=τln2. τ=3.5/0.693=5.05ms, R=τ/C=5.1kΩ.
Ejercicio 7.6★★
Muéstrese por integración directa que cargar un condensador de 0 a E a través de R disipa exactamente 21CE2 en R, con independencia de R. ¿En qué queda el rendimiento de la transferencia y cómo podría mejorarse?
Solución
Solución de Ejercicio 7.6.
i=(E/R)e−t/τ: ∫0∞Ri2dt=(E2/R)(τ/2)=21CE2, sin rastro de R. La fuente aporta CE2: rendimiento del 50% con cualquier resistor. Mejor: cargar a través de una bobina (la corriente ya no es proporcional a la diferencia de tensión) o en varios escalones de tensión — cada escalón de E/n solo pierde 21C(E/n)2.
Ejercicio 7.7★★
En un osciloscopio, una descarga libre de un RLC muestra máximos sucesivos de 2.0V, 1.5V y 1.12V separados 0.40ms. Dedúzcanse Q, ω0 y, con C=1.0µF, los valores de L y R.
Solución
Solución de Ejercicio 7.7.
δ=ln(2.0/1.5)=0.29 (y ln(1.5/1.12)=0.29: coherente); Q≈π/δ=11. T=0.40ms, luego ω0≈2π/T=1.57×104rad/s (la corrección 1/8Q2 es del 0.1%). L=1/(ω02C)=4.1mH; R=ω0L/Q=5.8Ω.
Ejercicio 7.8★★
Un condensador cargado a U0=100V (C=10µF) se conmuta sobre una bobina L=1.0mH de resistencia pequeña (Q≫1). Escríbanse uC(t) e i(t), y estímese la corriente máxima. ¿Dónde está la energía en ese máximo?
Solución
Solución de Ejercicio 7.8.
uC≈U0e−t/τcosω0t; i=−CduC/dt≈CU0ω0sinω0t (para t≪τ), con máximo U0C/L=100×0.1=10A, alcanzado un cuarto de periodo después de cerrar, cuando uC=0: toda la energía (21CU02=50mJ) es entonces magnética, 21Li2=50mJ.
Ejercicio 7.9★★
Un cuarto de coche: m=400kg sobre un resorte k=2.0×104N/m. Calcúlense ω0, f0 y el coeficiente de amortiguamiento crítico. Con un amortiguador α=2000Ns/m, dense Q, el régimen, el seudoperiodo y el tiempo de decaimiento τ.
Una bobina L=0.50H lleva I0=2.0A; se desconecta la fuente y la bobina queda sobre un resistor de descarga R=100Ω. Dense i(t), τ, la tensión máxima en bornes de la bobina y la energía disipada. ¿Qué ocurriría sin el resistor?
Solución
Solución de Ejercicio 7.10.
i=I0e−t/τ, τ=L/R=5.0ms; tensión máxima RI0=200V en t=0+; energía 21LI02=1.0J disipada en R. Sin el resistor la corriente no tiene camino: di/dt es enorme, la tensión de la bobina llega a los kilovoltios y salta un arco en el interruptor.
Ejercicio 7.11★★★
Un condensador C=1.0µF se carga desde E=12V a través de R1=10kΩ mientras R2=20kΩ está conectado en paralelo con él. Usando un equivalente de Thévenin, hállense la constante de tiempo y la tensión final; escríbase uC(t) a partir de uC(0)=0.
Solución
Solución de Ejercicio 7.11.
Visto desde C: ETh=ER2/(R1+R2)=8.0V, RTh=R1∥R2=6.67kΩ; τ=RThC=6.7ms; uC=8.0(1−e−t/τ) V.
Ejercicio 7.12★★★
Régimen crítico, escalón desde el reposo (x(0)=x′(0)=0): muéstrese que x=x∞[1−(1+ω0t)e−ω0t] y hállese numéricamente el tiempo necesario para alcanzar el 95% de x∞ (en unidades de 1/ω0). Compárese con el tiempo de establecimiento de la envolvente para Q=5 y con el de la raíz lenta cuando Q=0.1; conclúyase por qué el régimen crítico es el objetivo del ingeniero.
Solución
Solución de Ejercicio 7.12.
Raíz doble −ω0: x=x∞+(A+Bt)e−ω0t; x(0)=0 da A=−x∞ y x′(0)=0 da B=ω0A: x=x∞[1−(1+ω0t)e−ω0t]. Para el 95%: (1+s)e−s=0.05, s≈4.7: t95=4.7/ω0. Para Q=5 la envolvente e−t/τ con τ=10/ω0 necesita 3τ=30/ω0; para Q=0.1 la raíz lenta es ∣r+∣≈ω0Q y 3/∣r+∣=30/ω0. Ambos son seis veces más lentos: el amortiguamiento crítico es el retorno sin sobreoscilación más rápido, que es lo que quiere un aparato de medida, un cierrapuertas o una suspensión.
7.4 Problema: diseñar una suspensión
Problema 7.1
Problema de fin de semana — la misma ecuación diferencial en un banco de laboratorio y bajo un coche: cuántos rebotes permite un amortiguador gastado, qué elige en su lugar el diseñador y cómo un solo número distingue los dos casos
La Parte I estudia un circuito RLC serie, L=10mH, C=100nF; las Partes II–IV, un modelo de cuarto de coche: una masa m=350kg apoyada en un resorte k=22kN/m y un amortiguador que ejerce −αx˙. Sobreoscilación de una respuesta a un escalón en régimen seudoperiódico (dato): xmax/x∞−1=exp(−πξ/1−ξ2), ξ=1/2Q.
Parte I — El circuito.
Escríbase la ley de las mallas para el RLC serie sometido a una fuente E, y la ecuación diferencial de uC.
Pásese a la forma canónica; exprésense ω0 y Q.
Calcúlense ω0, f0 y la resistencia crítica Rc.
Escríbanse la ecuación característica y su discriminante; enúnciense los tres regímenes en función de Q.
En el caso seudoperiódico, dense τ y ω y la forma de uC(t).
Para R=100Ω calcúlense Q, τ, el seudoperiodo y el número de oscilaciones visibles.
El condensador, cargado a U0, se conmuta sobre L y R en t=0. Dense uC(0+) y duC/dt(0+), con la razón física de cada uno.
Parte II — El cuarto de coche.
Midiendo x desde la posición de equilibrio, escríbase la segunda ley de Newton para la carrocería y muéstrese que el peso desaparece.
Pásese la ecuación a la forma canónica; dense ω0 y Q en función de m, k y α, y enumérense las correspondencias electromecánicas.
Calcúlense ω0, f0 y el coeficiente crítico αc.
El diseñador monta α=4000Ns/m: calcúlense Q y ξ, y nómbrese el régimen.
Calcúlense τ y el seudoperiodo.
Calcúlese la sobreoscilación tras un escalón (un bordillo): ¿es cómodo el coche?
¿Qué resistencia R daría al circuito de la Parte I el mismo Q que esta suspensión?
Parte III — El amortiguador gastado. Se ha perdido aceite: α=1200Ns/m.
Calcúlense Q, el seudoperiodo y τ.
¿En qué factor se reducen los máximos sucesivos? ¿Cuántos rebotes se ven antes de que la amplitud baje del 5%?
La prueba del rebote: se empuja la esquina hacia abajo x0=5.0cm y se suelta. Dense las condiciones iniciales, la energía almacenada en el resorte y dígase adónde va esa energía.
Estímense la velocidad máxima de la carrocería durante el primer retorno (≈x0ω0 para Q≫1) y la fuerza máxima del amortiguador.
¿Por qué el amortiguamiento exactamente crítico (Q=21) no es la elección del diseñador, pese a no dar sobreoscilación? (Piénsese en los últimos milímetros del retorno.)
Calcúlese la sobreoscilación del amortiguador gastado y compárese con la Parte II.
Parte IV — Diagnosticar Q a partir de un registro. Un acelerómetro sobre la carrocería registra el decaimiento libre.
Para el coche gastado, los máximos sucesivos valen 1.00, 0.25 y 0.062 (unidades arbitrarias). Calcúlese el decremento logarítmico, dedúzcase ξ exactamente de δ=2πξ/1−ξ2 y después Q; compárese con la Parte III.
En el coche sano no se ve un segundo máximo; solo puede leerse la sobreoscilación, del 3.8%. Invirtiendo la fórmula de la sobreoscilación, hállense ξ y Q.
¿Qué fracción de la energía mecánica disipa por seudoperiodo el amortiguador gastado?
Un técnico que no tiene un coche a mano monta el RLC de la Parte I con una R elegida para que Q=2.3: ¿qué R, y por qué el registro del osciloscopio, representado frente a ω0t, tiene exactamente la misma forma que el del coche?
Resúmase: qué número adimensional decide la forma de cualquier transitorio de segundo orden, qué valor busca quien diseña una suspensión y qué debe concluir un conductor de “más de un rebote”.
16.e−T/τ=e−1.39=0.25: cada rebote es la cuarta parte del anterior; 0.253=1.6%: tres rebotes visibles.
17.x(0)=−x0, x˙(0)=0; energía 21kx02=0.5×22000×0.0025=28J, disipada en forma de calor en el aceite del amortiguador.
18.v≈x0ω0=0.05×7.93=0.40m/s; fuerza del amortiguador αv≈1200×0.40=480N.
19. Con Q=21 el retorno es (1+ω0t)e−ω0t, que se aproxima asintóticamente: los últimos milímetros tardan mucho. Un sistema ligeramente subamortiguado (ξ≈0.7) sobrepasa unos pocos por ciento, pero entra antes en una banda de tolerancia y se queda.
20.exp(−π×0.22/0.976)=exp(−0.70)=50%, frente al 3.8%: el coche gastado rebota media amplitud por encima del equilibrio.
21.δ=ln4=1.39 (dos veces, coherente); ξ=δ/4π2+δ2=1.39/6.43=0.215, Q=1/2ξ=2.3: el valor del amortiguador gastado.
22.−ln0.038=3.27=πξ/1−ξ2, luego ξ=3.27/π2+3.272=0.72, Q=0.69.
23. La energía es ∝amplitud2: fracción restante por periodo 0.252=6%, luego se disipa el 94% por seudoperiodo.
24.R=316/2.3=137Ω. En la variable ω0t la ecuación canónica solo contiene Q: dos sistemas con el mismo Q trazan la misma curva.
25. El factor de calidadQ (o ξ) fija por sí solo la forma; un diseñador apunta a Q≈0.7 (ξ≈0.7); más de un rebote visible significa un Q bastante mayor que 1: el amortiguador está gastado.