Physics · Libro 3 · Bachelor Year 1

Física universitaria — primer año

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.
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)x(t) obedece una ecuación lineal de primer orden con coeficientes constantes cuando

τ ⁣dx ⁣dt+x=x,\tau\,\frac{\dd x}{\dd t} + x = x_\infty ,

con τ>0\tau > 0 la constante de tiempo y xx_\infty una constante (el valor impuesto por la fuente). Su régimen permanente es x=xx = x_\infty; el régimen transitorio es la aproximación a él desde el valor inicial x(0)=x0x(0) = x_0.

Teorema 7.2 (Solución)

La única solución con x(0)=x0x(0) = x_0 es

x(t)=x+(x0x)et/τ.x(t) = x_\infty + (x_0 - x_\infty)\,\eu^{-t/\tau} .

La distancia al régimen permanente se reduce en el factor e10.37\eu^{-1} \approx 0.37 cada τ\tau: se ha recorrido el 63%63\% del camino en τ\tau, el 95%95\% en 3τ3\tau y el 99%99\% en 4.6τ4.6\tau. La tangente en el origen alcanza xx_\infty en t=τt = \tau.

Demostración. y=xxy = x - x_\infty obedece τy+y=0\tau y' + y = 0, cuyas soluciones son y=Ket/τy = K\eu^{-t/\tau} (el curso de matemáticas sobre ecuaciones diferenciales lineales demuestra que no hay otras); K=x0xK = x_0 - x_\infty. La pendiente en 00 vale (x0x)/τ-(x_0 - x_\infty)/\tau, que cerraría la distancia en un tiempo τ\tau.

La respuesta de primer orden a un escalón desde x_0 = 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 x_0 y x_∈fty.
La respuesta de primer orden a un escalón desde x0=0x_0 = 0: la tangente inicial corta la asíntota en t=τt = \tau; el 63%63\% del camino en τ\tau y el 95%95\% en 3τ3\tau. Todo transitorio de primer orden es esta curva, dilatada por τ\tau y desplazada por x0x_0 y xx_\infty.

Proposición 7.3 (Circuitos RC y RL)

Un condensador CC que se carga a través de un resistor RR desde una fuente ideal EE obedece RC ⁣duC/ ⁣dt+uC=ERC\,\dd u_C/\dd t + u_C = E: τ=RC\tau = RC, uC,=Eu_{C,\infty} = E. Una bobina LL en serie con RR sometida a EE obedece (L/R) ⁣di/ ⁣dt+i=E/R(L/R)\,\dd i/\dd t + i = E/R: τ=L/R\tau = L/R, i=E/Ri_\infty = E/R. Retirada la fuente (E0E \to 0), esas mismas ecuaciones describen la descarga hacia cero.

Demostración. Ley de las mallas E=Ri+uCE = Ri + u_C con i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t; ley de las mallas E=Ri+L ⁣di/ ⁣dtE = Ri + L\,\dd i/\dd t.

Proposición 7.4 (Condiciones de continuidad)

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 t0t_0, uC(t0+)=uC(t0)u_C(t_0^+) = u_C(t_0^-) e iL(t0+)=iL(t0)i_L(t_0^+) = i_L(t_0^-). Las corrientes en los resistores y las tensiones en las bobinas sí pueden saltar.

Demostración. Un salto de uCu_C exigiría i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t infinita, y un salto de iLi_L exigiría u=L ⁣di/ ⁣dtu = L\,\dd i/\dd t infinita; ambas cosas son imposibles en un circuito de tensiones y corrientes finitas (dicho de otro modo, las energías almacenadas 12CuC2\tfrac12 Cu_C^2 y 12Li2\tfrac12 Li^2 no pueden cambiar instantáneamente sin una potencia infinita).

Método 7.5 (Resolver un transitorio de primer orden)

  1. 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 uCu_C o en iLi_L y pásese a la forma τx+x=x\tau x' + x = x_\infty; léanse τ\tau y xx_\infty (este último es también el régimen permanente que se obtiene sustituyendo CC por un circuito abierto y LL por un hilo).
  2. Hállese x0x_0 por la continuidad de uCu_C o de iLi_L en el instante de conmutación, usando el circuito anterior a ella.
  3. Escríbase x=x+(x0x)et/τx = x_\infty + (x_0 - x_\infty)\eu^{-t/\tau}; 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)(E_{\mathrm{Th}}, R_{\mathrm{Th}}) (Capítulo 6): τ=RThC\tau = R_{\mathrm{Th}}C, u=EThu_\infty = E_{\mathrm{Th}}.

Proposición 7.6 (Balance energético de la carga RC)

Al cargar un condensador de 00 a EE a través de un resistor RR cualquiera, la fuente aporta CE2CE^2, el condensador almacena 12CE2\tfrac12 CE^2 y el resistor disipa el otro 12CE2\tfrac12 CE^2 — sea cual sea RR.

Demostración. i=(E/R)et/τi = (E/R)\eu^{-t/\tau}; la fuente aporta 0Ei ⁣dt=E(E/R)τ=CE2\int_0^\infty Ei\,\dd t = E \cdot (E/R)\tau = CE^2; el resistor toma 0Ri2 ⁣dt=(E2/R)0e2t/τ ⁣dt=(E2/R)(τ/2)=12CE2\int_0^\infty Ri^2\,\dd t = (E^2/R)\int_0^\infty \eu^{-2t/\tau}\dd t = (E^2/R)(\tau/2) = \tfrac12 CE^2; la diferencia es la energía almacenada 12CE2\tfrac12 CE^2. Una RR menor disipa la misma energía más deprisa.

Ejemplo 7.7 (La bobina de un relé)

L=20mHL = 20\,\mathrm{mH}, R=5.0ΩR = 5.0\,\Omega, E=10VE = 10\,\mathrm{V}: τ=4.0ms\tau = 4.0\,\mathrm{ms}, i=2.0Ai_\infty = 2.0\,\mathrm{A}; los contactos se cierran cuando ii llega a 1.5A1.5\,\mathrm{A}, es decir en t=τln(10.75)=5.5mst = -\tau\ln(1 - 0.75) = 5.5\,\mathrm{ms}. Al abrir el circuito, ii debe caer desde 2A2\,\mathrm{A}; un diodo de rueda libre en paralelo con la bobina le ofrece un camino a través de RR solo, con τ=4ms\tau = 4\,\mathrm{ms}, 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)x(t) es un sistema lineal de segundo orden cuando

 ⁣d2x ⁣dt2+ω0Q ⁣dx ⁣dt+ω02x=ω02x,\frac{\dd^2 x}{\dd t^2} + \frac{\omega_0}{Q}\,\frac{\dd x}{\dd t} + \omega_0^2\,x = \omega_0^2\,x_\infty ,

con ω0>0\omega_0 > 0 la frecuencia angular propia, Q>0Q > 0 el factor de calidad (adimensional) y xx_\infty el régimen permanente. También se escribe ω0/Q=2ξω0\omega_0/Q = 2\xi\omega_0, con ξ=1/2Q\xi = 1/2Q el factor de amortiguamiento.

Proposición 7.9 (El circuito RLC serie)

Un condensador CC, una bobina LL y un resistor RR en serie sometidos a una fuente EE: la tensión del condensador obedece la forma canónica con

ω0=1LC,Q=1RLC,uC,=E.\omega_0 = \frac{1}{\sqrt{LC}}, \qquad Q = \frac{1}{R}\sqrt{\frac{L}{C}}, \qquad u_{C,\infty} = E .

Demostración. Ley de las mallas: E=L ⁣di/ ⁣dt+Ri+uCE = L\,\dd i/\dd t + Ri + u_C con i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t: LCuC+RCuC+uC=ELC\,u_C'' + RC\,u_C' + u_C = E; divídase por LCLC e identifíquense ω02=1/LC\omega_0^2 = 1/LC y ω0/Q=R/L\omega_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.
Las dos caras del sistema de segundo orden: el circuito RLC serie (ω0=1/LC\omega_0 = 1/\sqrt{LC}, Q=L/C/RQ = \sqrt{L/C}/R) y el conjunto masa–resorte–amortiguador (ω0=k/m\omega_0 = \sqrt{k/m}, Q=km/αQ = \sqrt{km}/\alpha). La misma ecuación y los mismos tres regímenes.

Proposición 7.10 (El oscilador mecánico)

Una masa mm sobre un resorte de rigidez kk con un amortiguador viscoso de coeficiente α\alpha (fuerza αx˙-\alpha\dot x), separada xx de la posición de equilibrio, obedece mx¨+αx˙+kx=0m\ddot x + \alpha\dot x + kx = 0: la forma canónica con

ω0=km,Q=kmα.\omega_0 = \sqrt{\frac{k}{m}}, \qquad Q = \frac{\sqrt{km}}{\alpha} .

La correspondencia xqx \leftrightarrow q, x˙i\dot x \leftrightarrow i, mLm \leftrightarrow L, αR\alpha \leftrightarrow R, k1/Ck \leftrightarrow 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-kx y la fuerza de amortiguamiento; divídase por mm.

Teorema 7.11 (Los tres regímenes)

Las soluciones de la ecuación homogénea x+(ω0/Q)x+ω02x=0x'' + (\omega_0/Q)x' + \omega_0^2 x = 0 las gobiernan las raíces de r2+(ω0/Q)r+ω02=0r^2 + (\omega_0/Q)r + \omega_0^2 = 0, de discriminante Δ=ω02(1/Q24)\Delta = \omega_0^2(1/Q^2 - 4):

  • Q<12Q < \tfrac12, aperiódico: dos raíces reales negativas r±=ω02Q(114Q2)r_\pm = -\frac{\omega_0}{2Q}\big(1 \mp \sqrt{1 - 4Q^2}\big), x=Aer+t+Bertx = A\eu^{r_+t} + B\eu^{r_-t}, un retorno sin oscilación;
  • Q=12Q = \tfrac12, crítico: una raíz doble r=ω0r = -\omega_0, x=(A+Bt)eω0tx = (A + Bt)\eu^{-\omega_0 t}, el retorno más rápido sin sobreoscilación;
  • Q>12Q > \tfrac12, seudoperiódico: raíces complejas r=1/τ±iωr = -1/\tau \pm \iu\omega con

    τ=2Qω0,ω=ω0114Q2,x=et/τ(Acosωt+Bsinωt),\tau = \frac{2Q}{\omega_0}, \qquad \omega = \omega_0\sqrt{1 - \frac{1}{4Q^2}}, \qquad x = \eu^{-t/\tau}\big(A\cos\omega t + B\sin\omega t\big),

    una oscilación amortiguada de seudoperiodo T=2π/ωT = 2\pi/\omega dentro de la envolvente ±A2+B2et/τ\pm\sqrt{A^2 + B^2}\,\eu^{-t/\tau}.

La solución completa con segundo miembro constante es xx_\infty más la solución homogénea; AA y BB se obtienen de x(0)x(0) y x(0)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±12Δr = -\omega_0/2Q \pm \tfrac12\sqrt\Delta. Para Q>12Q > \tfrac12, Δ=iω041/Q2\sqrt\Delta = \iu\omega_0\sqrt{4 - 1/Q^2}, lo que da la parte real ω0/2Q=1/τ-\omega_0/2Q = -1/\tau y la parte imaginaria ±ω011/4Q2\pm\omega_0\sqrt{1 - 1/4Q^2}; las soluciones reales son las combinaciones de et/τcosωt\eu^{-t/\tau}\cos\omega t y et/τsinωt\eu^{-t/\tau}\sin\omega 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 = 1/2: el retorno más rápido sin sobreoscilación; Q grande: una oscilación que dura alrededor de Q periodos.
Respuesta a un escalón, partiendo del reposo, de un sistema de segundo orden, para tres factores de calidad. QQ pequeño: una subida lenta; Q=12Q = \tfrac12: el retorno más rápido sin sobreoscilación; QQ grande: una oscilación que dura alrededor de QQ periodos.

Proposición 7.12 (Leer QQ en un registro seudoperiódico)

En régimen seudoperiódico, dos máximos sucesivos están en la razón xn+1/xn=eT/τx_{n+1}/x_n = \eu^{-T/\tau}; el decremento logarítmico

δ=lnxnxn+1=Tτ=πQ11/4Q2πQ(Q1)\delta = \ln\frac{x_n}{x_{n+1}} = \frac{T}{\tau} = \frac{\pi}{Q\sqrt{1 - 1/4Q^2}} \approx \frac{\pi}{Q} \quad (Q \gg 1)

da QQ; el número de oscilaciones visibles antes de que la amplitud baje del 5%5\% es del orden de QQ. Para Q1Q \gg 1, ωω0\omega \approx \omega_0 y T2π/ω0T \approx 2\pi/\omega_0.

Demostración. La envolvente et/τ\eu^{-t/\tau} se multiplica por eT/τ\eu^{-T/\tau} en cada periodo; T/τ=(2π/ω)(ω0/2Q)T/\tau = (2\pi/\omega)(\omega_0/2Q). La envolvente llega a e35%\eu^{-3} \approx 5\% en 3τ=6Q/ω0QT3\tau = 6Q/\omega_0 \approx Q\,T, es decir al cabo de unos QQ periodos.

Decaimiento libre seudoperiódico, x = x_0 -t/ t: los máximos sucesivos se reducen en el factor constante -T/, y (x_n/x_n+1) = T/ π/Q mide el factor de calidad a partir del registro.
Decaimiento libre seudoperiódico, x=x0et/τcosωtx = x_0\eu^{-t/\tau}\cos\omega t: los máximos sucesivos se reducen en el factor constante eT/τ\eu^{-T/\tau}, y ln(xn/xn+1)=T/τπ/Q\ln(x_n/x_{n+1}) = T/\tau \approx \pi/Q mide el factor de calidad a partir del registro.

Ejemplo 7.13 (Una oscilación RLC)

L=10mHL = 10\,\mathrm{mH}, C=100nFC = 100\,\mathrm{nF}, R=100ΩR = 100\,\Omega: ω0=1/109=3.16×104rad/s\omega_0 = 1/\sqrt{10^{-9}} = 3.16 \times 10^{4}\,\mathrm{rad}/\mathrm{s} (f0=5.03kHzf_0 = 5.03\,\mathrm{kHz}), Q=L/C/R=316/100=3.16Q = \sqrt{L/C}/R = 316/100 = 3.16: seudoperiódico, τ=2Q/ω0=0.20ms\tau = 2Q/\omega_0 = 0.20\,\mathrm{ms}, ω=0.987ω0\omega = 0.987\,\omega_0, T=0.20msT = 0.20\,\mathrm{ms}, tres oscilaciones visibles. La resistencia crítica vale Rc=2L/C=632ΩR_c = 2\sqrt{L/C} = 632\,\Omega; por encima de ella la descarga es aperiódica.

Método 7.14 (Resolver un transitorio de segundo orden)

  1. Escríbase la ecuación; pásese a la forma canónica; léanse ω0\omega_0, QQ y xx_\infty.
  2. Decídase el régimen a partir de QQ y escríbase la solución general como xx_\infty más la solución homogénea.
  3. Hállense x(0)x(0) y x(0)x'(0) por la continuidad de uCu_C y de iLi_L (para el RLC: uC(0+)=uC(0)u_C(0^+) = u_C(0^-) y uC(0+)=i(0)/Cu_C'(0^+) = i(0^-)/C); despéjense las dos constantes.

Ejemplo 7.15 (Condiciones iniciales)

El condensador del ejemplo anterior se carga a U0U_0 y, en t=0t = 0, se cierra sobre LL y RR (sin fuente): x=0x_\infty = 0, uC(0)=U0u_C(0) = U_0, uC(0)=i(0)/C=0u_C'(0) = i(0)/C = 0, pues la corriente de la bobina era nula. Por tanto A=U0A = U_0 y A/τ+Bω=0-A/\tau + B\omega = 0, B=U0/(ωτ)=U0/4Q21B = U_0/(\omega\tau) = U_0/\sqrt{4Q^2 - 1}: uC=U0et/τ[cosωt+sin(ωt)/4Q21]u_C = U_0\eu^{-t/\tau}[\cos\omega t + \sin(\omega t)/\sqrt{4Q^2 - 1}] — para Q1Q \gg 1, simplemente U0et/τcosω0tU_0\eu^{-t/\tau}\cos\omega_0 t, y la corriente i=C ⁣duC/ ⁣dti = -C\,\dd u_C/\dd t tiene un máximo cercano a U0C/L=U0/(ω0L)U_0\sqrt{C/L} = U_0/(\omega_0 L).

Observación 7.16 (Energía)

La energía almacenada 12CuC2+12Li2\tfrac12 Cu_C^2 + \tfrac12 Li^2 (o 12kx2+12mx˙2\tfrac12 kx^2 + \tfrac12 m\dot x^2) decrece a razón de Ri2Ri^2 (αx˙2\alpha\dot x^2): derívese y úsese la ecuación. En cada seudoperiodo la energía se reduce en el factor e2T/τ\eu^{-2T/\tau}; con Q=3Q = 3 se pierde alrededor del 88%88\% por oscilación. Un QQ alto significa una fuga lenta comparada con la oscilación — un buen reloj, una resonancia aguda (Capítulo 8) — y un QQ 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ΩR = 10\,\mathrm{k}\Omega, C=47nFC = 47\,\mathrm{nF}, cargado desde EE. Calcúlense τ\tau, uC(τ)/Eu_C(\tau)/E y el tiempo para llegar al 99%99\% de EE.

Solución

Solución de Ejercicio 7.1.

τ=RC=104×47×109=0.47ms\tau = RC = 10^4 \times 47\times10^{-9} = 0.47\,\mathrm{ms}; uC(τ)/E=1e1=0.63u_C(\tau)/E = 1 - \eu^{-1} = 0.63; 99%99\%: et/τ=0.01\eu^{-t/\tau} = 0.01, t=τln100=4.6τ=2.2mst = \tau\ln 100 = 4.6\tau = 2.2\,\mathrm{ms}.

Ejercicio 7.2

Un circuito RL: L=20mHL = 20\,\mathrm{mH}, R=5.0ΩR = 5.0\,\Omega, E=10VE = 10\,\mathrm{V}. Calcúlense τ\tau, la corriente final, la corriente en t=τt = \tau 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\tau = L/R = 4.0\,\mathrm{ms}; i=E/R=2.0Ai_\infty = E/R = 2.0\,\mathrm{A}; i(τ)=0.63×2.0=1.26Ai(\tau) = 0.63 \times 2.0 = 1.26\,\mathrm{A}; en 0+0^+, i=0i = 0 (continuidad), luego uL=ERi=10Vu_L = E - Ri = 10\,\mathrm{V}.

Ejercicio 7.3

Un RLC serie: L=10mHL = 10\,\mathrm{mH}, C=100nFC = 100\,\mathrm{nF}, R=100ΩR = 100\,\Omega. Calcúlense ω0\omega_0, f0f_0 y QQ; nómbrese el régimen; dense τ\tau y el seudoperiodo.

Solución

Solución de Ejercicio 7.3.

ω0=1/LC=3.16×104rad/s\omega_0 = 1/\sqrt{LC} = 3.16 \times 10^{4}\,\mathrm{rad}/\mathrm{s}, f0=5.03kHzf_0 = 5.03\,\mathrm{kHz}; Q=L/C/R=316/100=3.2>12Q = \sqrt{L/C}/R = 316/100 = 3.2 > \tfrac12: seudoperiódico; τ=2Q/ω0=0.20ms\tau = 2Q/\omega_0 = 0.20\,\mathrm{ms}; ω=ω011/4Q2=0.987ω0\omega = \omega_0\sqrt{1 - 1/4Q^2} = 0.987\omega_0, T=0.20msT = 0.20\,\mathrm{ms}.

Ejercicio 7.4

Con los mismos LL y CC, ¿qué resistencia da el régimen crítico? ¿Cuánto vale entonces QQ? ¿Qué ocurre para R=2kΩR = 2\,\mathrm{k}\Omega?

Solución

Solución de Ejercicio 7.4.

Rc=2L/C=632ΩR_c = 2\sqrt{L/C} = 632\,\Omega, Q=12Q = \tfrac12. A 2kΩ2\,\mathrm{k}\Omega, Q=0.16Q = 0.16: aperiódico, una subida lenta gobernada por la raíz pequeña r+ω0Q=5×103s1\abs{r_+} \approx \omega_0 Q = 5 \times 10^{3}\,\mathrm{s}^{-1} (escala de tiempo 0.2ms0.2\,\mathrm{ms}, pero sin oscilación).

Ejercicio 7.5 ★★

Un condensador cargado C=1.0µFC = 1.0\,\text{µ}\mathrm{F} se descarga a través de una RR desconocida; la tensión se reduce a la mitad cada 3.5ms3.5\,\mathrm{ms}. Muéstrese que el tiempo de semirreducción vale τln2\tau\ln 2 y hállese RR.

Solución

Solución de Ejercicio 7.5.

u=U0et/τu = U_0\eu^{-t/\tau} se reduce a la mitad cuando et/τ=12\eu^{-t/\tau} = \tfrac12: t1/2=τln2t_{1/2} = \tau\ln 2. τ=3.5/0.693=5.05ms\tau = 3.5/0.693 = 5.05\,\mathrm{ms}, R=τ/C=5.1kΩR = \tau/C = 5.1\,\mathrm{k}\Omega.

Ejercicio 7.6 ★★

Muéstrese por integración directa que cargar un condensador de 00 a EE a través de RR disipa exactamente 12CE2\tfrac12 CE^2 en RR, con independencia de RR. ¿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)et/τi = (E/R)\eu^{-t/\tau}: 0Ri2 ⁣dt=(E2/R)(τ/2)=12CE2\int_0^\infty Ri^2\,\dd t = (E^2/R)(\tau/2) = \tfrac12 CE^2, sin rastro de RR. La fuente aporta CE2CE^2: rendimiento del 50%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/nE/n solo pierde 12C(E/n)2\tfrac12 C(E/n)^2.

Ejercicio 7.7 ★★

En un osciloscopio, una descarga libre de un RLC muestra máximos sucesivos de 2.0V2.0\,\mathrm{V}, 1.5V1.5\,\mathrm{V} y 1.12V1.12\,\mathrm{V} separados 0.40ms0.40\,\mathrm{ms}. Dedúzcanse QQ, ω0\omega_0 y, con C=1.0µFC = 1.0\,\text{µ}\mathrm{F}, los valores de LL y RR.

Solución

Solución de Ejercicio 7.7.

δ=ln(2.0/1.5)=0.29\delta = \ln(2.0/1.5) = 0.29 (y ln(1.5/1.12)=0.29\ln(1.5/1.12) = 0.29: coherente); Qπ/δ=11Q \approx \pi/\delta = 11. T=0.40msT = 0.40\,\mathrm{ms}, luego ω02π/T=1.57×104rad/s\omega_0 \approx 2\pi/T = 1.57 \times 10^{4}\,\mathrm{rad}/\mathrm{s} (la corrección 1/8Q21/8Q^2 es del 0.1%0.1\%). L=1/(ω02C)=4.1mHL = 1/(\omega_0^2 C) = 4.1\,\mathrm{mH}; R=ω0L/Q=5.8ΩR = \omega_0 L/Q = 5.8\,\Omega.

Ejercicio 7.8 ★★

Un condensador cargado a U0=100VU_0 = 100\,\mathrm{V} (C=10µFC = 10\,\text{µ}\mathrm{F}) se conmuta sobre una bobina L=1.0mHL = 1.0\,\mathrm{mH} de resistencia pequeña (Q1Q \gg 1). Escríbanse uC(t)u_C(t) e i(t)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.

uCU0et/τcosω0tu_C \approx U_0\eu^{-t/\tau}\cos\omega_0 t; i=C ⁣duC/ ⁣dtCU0ω0sinω0ti = -C\,\dd u_C/\dd t \approx CU_0\omega_0\sin\omega_0 t (para tτt \ll \tau), con máximo U0C/L=100×0.1=10AU_0\sqrt{C/L} = 100 \times 0.1 = 10\,\mathrm{A}, alcanzado un cuarto de periodo después de cerrar, cuando uC=0u_C = 0: toda la energía (12CU02=50mJ\tfrac12 CU_0^2 = 50\,\mathrm{mJ}) es entonces magnética, 12Li2=50mJ\tfrac12 Li^2 = 50\,\mathrm{mJ}.

Ejercicio 7.9 ★★

Un cuarto de coche: m=400kgm = 400\,\mathrm{kg} sobre un resorte k=2.0×104N/mk = 2.0 \times 10^{4}\,\mathrm{N}/\mathrm{m}. Calcúlense ω0\omega_0, f0f_0 y el coeficiente de amortiguamiento crítico. Con un amortiguador α=2000Ns/m\alpha = 2000\,\mathrm{N}\,\mathrm{s}/\mathrm{m}, dense QQ, el régimen, el seudoperiodo y el tiempo de decaimiento τ\tau.

Solución

Solución de Ejercicio 7.9.

ω0=k/m=7.1rad/s\omega_0 = \sqrt{k/m} = 7.1\,\mathrm{rad}/\mathrm{s}, f0=1.1Hzf_0 = 1.1\,\mathrm{Hz}; αc=2km=28×106=5.7×103Ns/m\alpha_c = 2\sqrt{km} = 2\sqrt{8\times10^6} = 5.7 \times 10^{3}\,\mathrm{N}\,\mathrm{s}/\mathrm{m}. α=2000\alpha = 2000: Q=km/α=2828/2000=1.4Q = \sqrt{km}/\alpha = 2828/2000 = 1.4, seudoperiódico; T=2π/(ω011/4Q2)=0.89/0.935=0.95sT = 2\pi/(\omega_0\sqrt{1 - 1/4Q^2}) = 0.89/0.935 = 0.95\,\mathrm{s}; τ=2Q/ω0=0.40s\tau = 2Q/\omega_0 = 0.40\,\mathrm{s}.

Ejercicio 7.10 ★★★

Una bobina L=0.50HL = 0.50\,\mathrm{H} lleva I0=2.0AI_0 = 2.0\,\mathrm{A}; se desconecta la fuente y la bobina queda sobre un resistor de descarga R=100ΩR = 100\,\Omega. Dense i(t)i(t), τ\tau, 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=I0et/τi = I_0\eu^{-t/\tau}, τ=L/R=5.0ms\tau = L/R = 5.0\,\mathrm{ms}; tensión máxima RI0=200VRI_0 = 200\,\mathrm{V} en t=0+t = 0^+; energía 12LI02=1.0J\tfrac12 LI_0^2 = 1.0\,\mathrm{J} disipada en RR. Sin el resistor la corriente no tiene camino:  ⁣di/ ⁣dt\dd i/\dd t 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µFC = 1.0\,\text{µ}\mathrm{F} se carga desde E=12VE = 12\,\mathrm{V} a través de R1=10kΩR_1 = 10\,\mathrm{k}\Omega mientras R2=20kΩR_2 = 20\,\mathrm{k}\Omega 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)u_C(t) a partir de uC(0)=0u_C(0) = 0.

Solución

Solución de Ejercicio 7.11.

Visto desde CC: ETh=ER2/(R1+R2)=8.0VE_{\mathrm{Th}} = ER_2/(R_1 + R_2) = 8.0\,\mathrm{V}, RTh=R1R2=6.67kΩR_{\mathrm{Th}} = R_1 \parallel R_2 = 6.67\,\mathrm{k}\Omega; τ=RThC=6.7ms\tau = R_{\mathrm{Th}}C = 6.7\,\mathrm{ms}; uC=8.0(1et/τ)u_C = 8.0(1 - \eu^{-t/\tau}) V.

Ejercicio 7.12 ★★★

Régimen crítico, escalón desde el reposo (x(0)=x(0)=0x(0) = x'(0) = 0): muéstrese que x=x[1(1+ω0t)eω0t]x = x_\infty[1 - (1 + \omega_0 t)\eu^{-\omega_0 t}] y hállese numéricamente el tiempo necesario para alcanzar el 95%95\% de xx_\infty (en unidades de 1/ω01/\omega_0). Compárese con el tiempo de establecimiento de la envolvente para Q=5Q = 5 y con el de la raíz lenta cuando Q=0.1Q = 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-\omega_0: x=x+(A+Bt)eω0tx = x_\infty + (A + Bt)\eu^{-\omega_0 t}; x(0)=0x(0) = 0 da A=xA = -x_\infty y x(0)=0x'(0) = 0 da B=ω0AB = \omega_0 A: x=x[1(1+ω0t)eω0t]x = x_\infty[1 - (1 + \omega_0 t)\eu^{-\omega_0 t}]. Para el 95%95\%: (1+s)es=0.05(1 + s)\eu^{-s} = 0.05, s4.7s \approx 4.7: t95=4.7/ω0t_{95} = 4.7/\omega_0. Para Q=5Q = 5 la envolvente et/τ\eu^{-t/\tau} con τ=10/ω0\tau = 10/\omega_0 necesita 3τ=30/ω03\tau = 30/\omega_0; para Q=0.1Q = 0.1 la raíz lenta es r+ω0Q\abs{r_+} \approx \omega_0 Q y 3/r+=30/ω03/\abs{r_+} = 30/\omega_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=10mHL = 10\,\mathrm{mH}, C=100nFC = 100\,\mathrm{nF}; las Partes II–IV, un modelo de cuarto de coche: una masa m=350kgm = 350\,\mathrm{kg} apoyada en un resorte k=22kN/mk = 22\,\mathrm{kN}/\mathrm{m} y un amortiguador que ejerce αx˙-\alpha\dot x. Sobreoscilación de una respuesta a un escalón en régimen seudoperiódico (dato): xmax/x1=exp(πξ/1ξ2)x_{\max}/x_\infty - 1 = \exp(-\pi\xi/\sqrt{1 - \xi^2}), ξ=1/2Q\xi = 1/2Q.

Parte I — El circuito.

  1. Escríbase la ley de las mallas para el RLC serie sometido a una fuente EE, y la ecuación diferencial de uCu_C.
  2. Pásese a la forma canónica; exprésense ω0\omega_0 y QQ.
  3. Calcúlense ω0\omega_0, f0f_0 y la resistencia crítica RcR_c.
  4. Escríbanse la ecuación característica y su discriminante; enúnciense los tres regímenes en función de QQ.
  5. En el caso seudoperiódico, dense τ\tau y ω\omega y la forma de uC(t)u_C(t).
  6. Para R=100ΩR = 100\,\Omega calcúlense QQ, τ\tau, el seudoperiodo y el número de oscilaciones visibles.
  7. El condensador, cargado a U0U_0, se conmuta sobre LL y RR en t=0t = 0. Dense uC(0+)u_C(0^+) y  ⁣duC/ ⁣dt(0+)\dd u_C/\dd t(0^+), con la razón física de cada uno.

Parte II — El cuarto de coche.

  1. Midiendo xx desde la posición de equilibrio, escríbase la segunda ley de Newton para la carrocería y muéstrese que el peso desaparece.
  2. Pásese la ecuación a la forma canónica; dense ω0\omega_0 y QQ en función de mm, kk y α\alpha, y enumérense las correspondencias electromecánicas.
  3. Calcúlense ω0\omega_0, f0f_0 y el coeficiente crítico αc\alpha_c.
  4. El diseñador monta α=4000Ns/m\alpha = 4000\,\mathrm{N}\,\mathrm{s}/\mathrm{m}: calcúlense QQ y ξ\xi, y nómbrese el régimen.
  5. Calcúlense τ\tau y el seudoperiodo.
  6. Calcúlese la sobreoscilación tras un escalón (un bordillo): ¿es cómodo el coche?
  7. ¿Qué resistencia RR daría al circuito de la Parte I el mismo QQ que esta suspensión?

Parte III — El amortiguador gastado. Se ha perdido aceite: α=1200Ns/m\alpha = 1200\,\mathrm{N}\,\mathrm{s}/\mathrm{m}.

  1. Calcúlense QQ, el seudoperiodo y τ\tau.
  2. ¿En qué factor se reducen los máximos sucesivos? ¿Cuántos rebotes se ven antes de que la amplitud baje del 5%5\%?
  3. La prueba del rebote: se empuja la esquina hacia abajo x0=5.0cmx_0 = 5.0\,\mathrm{cm} y se suelta. Dense las condiciones iniciales, la energía almacenada en el resorte y dígase adónde va esa energía.
  4. Estímense la velocidad máxima de la carrocería durante el primer retorno (x0ω0\approx x_0\omega_0 para Q1Q \gg 1) y la fuerza máxima del amortiguador.
  5. ¿Por qué el amortiguamiento exactamente crítico (Q=12Q = \tfrac12) no es la elección del diseñador, pese a no dar sobreoscilación? (Piénsese en los últimos milímetros del retorno.)
  6. Calcúlese la sobreoscilación del amortiguador gastado y compárese con la Parte II.

Parte IV — Diagnosticar QQ a partir de un registro. Un acelerómetro sobre la carrocería registra el decaimiento libre.

  1. Para el coche gastado, los máximos sucesivos valen 1.001.00, 0.250.25 y 0.0620.062 (unidades arbitrarias). Calcúlese el decremento logarítmico, dedúzcase ξ\xi exactamente de δ=2πξ/1ξ2\delta = 2\pi\xi/\sqrt{1 - \xi^2} y después QQ; compárese con la Parte III.
  2. En el coche sano no se ve un segundo máximo; solo puede leerse la sobreoscilación, del 3.8%3.8\%. Invirtiendo la fórmula de la sobreoscilación, hállense ξ\xi y QQ.
  3. ¿Qué fracción de la energía mecánica disipa por seudoperiodo el amortiguador gastado?
  4. Un técnico que no tiene un coche a mano monta el RLC de la Parte I con una RR elegida para que Q=2.3Q = 2.3: ¿qué RR, y por qué el registro del osciloscopio, representado frente a ω0t\omega_0 t, tiene exactamente la misma forma que el del coche?
  5. 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”.
Solución

Solución de Problema 7.1.

1. E=L ⁣di/ ⁣dt+Ri+uCE = L\,\dd i/\dd t + Ri + u_C, i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t: LCuC+RCuC+uC=ELC\,u_C'' + RC\,u_C' + u_C = E.

2. uC+(ω0/Q)uC+ω02uC=ω02Eu_C'' + (\omega_0/Q)u_C' + \omega_0^2u_C = \omega_0^2E con ω0=1/LC\omega_0 = 1/\sqrt{LC}, Q=L/C/RQ = \sqrt{L/C}/R.

3. ω0=3.16×104rad/s\omega_0 = 3.16 \times 10^{4}\,\mathrm{rad}/\mathrm{s}, f0=5.03kHzf_0 = 5.03\,\mathrm{kHz}, Rc=2L/C=632ΩR_c = 2\sqrt{L/C} = 632\,\Omega.

4. r2+(ω0/Q)r+ω02=0r^2 + (\omega_0/Q)r + \omega_0^2 = 0, Δ=ω02(1/Q24)\Delta = \omega_0^2(1/Q^2 - 4): Q<12Q < \tfrac12 aperiódico, Q=12Q = \tfrac12 crítico, Q>12Q > \tfrac12 seudoperiódico.

5. τ=2Q/ω0\tau = 2Q/\omega_0, ω=ω011/4Q2\omega = \omega_0\sqrt{1 - 1/4Q^2}, uC=E+et/τ(Acosωt+Bsinωt)u_C = E + \eu^{-t/\tau}(A\cos\omega t + B\sin\omega t).

6. Q=3.16Q = 3.16, τ=0.20ms\tau = 0.20\,\mathrm{ms}, T=2π/ω=0.20msT = 2\pi/\omega = 0.20\,\mathrm{ms}, unas tres oscilaciones.

7. uC(0+)=U0u_C(0^+) = U_0 (uCu_C continua: corriente finita); uC(0+)=i(0+)/C=0u_C'(0^+) = i(0^+)/C = 0 (iLi_L continua y nula antes).

8. mx¨=k(x+xeq)αx˙+mgm\ddot x = -k(x + x_{\text{eq}}) - \alpha\dot x + mg con kxeq=mgkx_{\text{eq}} = mg: mx¨+αx˙+kx=0m\ddot x + \alpha\dot x + kx = 0.

9. ω0=k/m\omega_0 = \sqrt{k/m}, Q=km/αQ = \sqrt{km}/\alpha; xqx \leftrightarrow q, x˙i\dot x \leftrightarrow i, mLm \leftrightarrow L, αR\alpha \leftrightarrow R, k1/Ck \leftrightarrow 1/C.

10. ω0=22000/350=7.93rad/s\omega_0 = \sqrt{22000/350} = 7.93\,\mathrm{rad}/\mathrm{s}, f0=1.26Hzf_0 = 1.26\,\mathrm{Hz}; αc=2km=5.55×103Ns/m\alpha_c = 2\sqrt{km} = 5.55 \times 10^{3}\,\mathrm{N}\,\mathrm{s}/\mathrm{m}.

11. Q=2775/4000=0.69Q = 2775/4000 = 0.69, ξ=0.72\xi = 0.72: seudoperiódico, justo por debajo del crítico.

12. τ=2Q/ω0=0.18s\tau = 2Q/\omega_0 = 0.18\,\mathrm{s}; ω=ω01ξ2=5.5rad/s\omega = \omega_0\sqrt{1 - \xi^2} = 5.5\,\mathrm{rad}/\mathrm{s}, T=1.1sT = 1.1\,\mathrm{s}.

13. Sobreoscilación exp(π×0.72/0.69)=exp(3.3)=3.8%\exp(-\pi \times 0.72/0.69) = \exp(-3.3) = 3.8\%: una pequeña subida por encima del equilibrio y se acabó — cómodo.

14. R=L/C/Q=316/0.69=455ΩR = \sqrt{L/C}/Q = 316/0.69 = 455\,\Omega.

15. Q=2775/1200=2.3Q = 2775/1200 = 2.3, ξ=0.22\xi = 0.22; ω=0.976ω0\omega = 0.976\omega_0, T=0.81sT = 0.81\,\mathrm{s}; τ=2Q/ω0=0.58s\tau = 2Q/\omega_0 = 0.58\,\mathrm{s}.

16. eT/τ=e1.39=0.25\eu^{-T/\tau} = \eu^{-1.39} = 0.25: cada rebote es la cuarta parte del anterior; 0.253=1.6%0.25^3 = 1.6\%: tres rebotes visibles.

17. x(0)=x0x(0) = -x_0, x˙(0)=0\dot x(0) = 0; energía 12kx02=0.5×22000×0.0025=28J\tfrac12 kx_0^2 = 0.5 \times 22000 \times 0.0025 = 28\,\mathrm{J}, disipada en forma de calor en el aceite del amortiguador.

18. vx0ω0=0.05×7.93=0.40m/sv \approx x_0\omega_0 = 0.05 \times 7.93 = 0.40\,\mathrm{m}/\mathrm{s}; fuerza del amortiguador αv1200×0.40=480N\alpha v \approx 1200 \times 0.40 = 480\,\mathrm{N}.

19. Con Q=12Q = \tfrac12 el retorno es (1+ω0t)eω0t(1 + \omega_0 t)\eu^{-\omega_0 t}, que se aproxima asintóticamente: los últimos milímetros tardan mucho. Un sistema ligeramente subamortiguado (ξ0.7\xi \approx 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%\exp(-\pi \times 0.22/0.976) = \exp(-0.70) = 50\%, frente al 3.8%3.8\%: el coche gastado rebota media amplitud por encima del equilibrio.

21. δ=ln4=1.39\delta = \ln 4 = 1.39 (dos veces, coherente); ξ=δ/4π2+δ2=1.39/6.43=0.215\xi = \delta/\sqrt{4\pi^2 + \delta^2} = 1.39/6.43 = 0.215, Q=1/2ξ=2.3Q = 1/2\xi = 2.3: el valor del amortiguador gastado.

22. ln0.038=3.27=πξ/1ξ2-\ln 0.038 = 3.27 = \pi\xi/\sqrt{1 - \xi^2}, luego ξ=3.27/π2+3.272=0.72\xi = 3.27/\sqrt{\pi^2 + 3.27^2} = 0.72, Q=0.69Q = 0.69.

23. La energía es \propto amplitud2^2: fracción restante por periodo 0.252=6%0.25^2 = 6\%, luego se disipa el 94%94\% por seudoperiodo.

24. R=316/2.3=137ΩR = 316/2.3 = 137\,\Omega. En la variable ω0t\omega_0 t la ecuación canónica solo contiene QQ: dos sistemas con el mismo QQ trazan la misma curva.

25. El factor de calidad QQ (o ξ\xi) fija por sí solo la forma; un diseñador apunta a Q0.7Q \approx 0.7 (ξ0.7\xi \approx 0.7); más de un rebote visible significa un QQ bastante mayor que 11: el amortiguador está gastado.

Términos definidos en este capítulo

Ver los 393 términos del glosario