Physics · Boek 3 · Bachelor Year 1

Universitaire natuurkunde — jaar 1

Universitaire natuurkunde — jaar 1 · Bachelor Year 1

7Overgangsverschijnselen: eerste en tweede orde

Duw de hoek van een geparkeerde auto omlaag en laat los. Een auto in goede staat komt één keer omhoog en staat stil; een auto met versleten schokdempers wipt twee of drie keer voordat hij tot rust komt. In een elektronicalab verschijnen dezelfde twee gedragingen op een oscilloscoop wanneer een condensator zich over een spoel ontlaadt: een vloeiende terugkeer, of een uitdovende naklank. Beide zijn het overgangsverschijnsel tussen de ene evenwichtstoestand en de andere, en beide voldoen aan dezelfde kleine familie lineaire differentiaalvergelijkingen — eerste orde voor één energiereservoir, tweede orde voor twee die uitwisselen. Dit hoofdstuk lost ze eens en voor altijd op, benoemt hun parameters (tijdconstante, eigenfrequentie, kwaliteitsfactor) en laat zien hoe dezelfde getallen een schakeling én een vering beschrijven.

Een spoel, een condensator en een weerstand op een experimenteerbord, en op de oscilloscoop de uitdovende trilling die op een sprong volgt: het pseudoperiodieke regime van dit hoofdstuk.
Een spoel, een condensator en een weerstand op een experimenteerbord, en op de oscilloscoop de uitdovende trilling die op een sprong volgt: het pseudoperiodieke regime van dit hoofdstuk.

7.1 Systemen van de eerste orde

Definitie 7.1 (Lineair systeem van de eerste orde)

Een grootheid x(t)x(t) voldoet aan een lineaire eerste-ordevergelijking met constante coëfficiënten wanneer

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

met τ>0\tau > 0 de tijdconstante en xx_\infty een constante (de waarde die de bron oplegt). Haar eindtoestand is x=xx = x_\infty; het overgangsverschijnsel is de weg daarheen vanaf de beginwaarde x(0)=x0x(0) = x_0.

Stelling 7.2 (Oplossing)

De enige oplossing met x(0)=x0x(0) = x_0 is

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

Het gat tot de eindtoestand krimpt met de factor e10.37\eu^{-1} \approx 0.37 per τ\tau: 63%63\% van de weg na τ\tau, 95%95\% na 3τ3\tau, 99%99\% na 4.6τ4.6\tau. De raaklijn in de oorsprong bereikt xx_\infty op t=τt = \tau.

Bewijs. y=xxy = x - x_\infty voldoet aan τy+y=0\tau y' + y = 0, waarvan de oplossingen y=Ket/τy = K\eu^{-t/\tau} zijn (de wiskunde over lineaire differentiaalvergelijkingen bewijst dat er geen andere zijn); K=x0xK = x_0 - x_\infty. De helling in 00 is (x0x)/τ-(x_0 - x_\infty)/\tau, die het gat in de tijd τ\tau zou dichten.

Het eerste-ordeantwoord op een sprong vanuit x_0 = 0: de raaklijn in de oorsprong snijdt de asymptoot op t =; 63\% van de weg na , 95\% na 3. Elk eerste-orde-overgangsverschijnsel is deze kromme, uitgerekt met  en verschoven met x_0 en x_∈fty.
Het eerste-ordeantwoord op een sprong vanuit x0=0x_0 = 0: de raaklijn in de oorsprong snijdt de asymptoot op t=τt = \tau; 63%63\% van de weg na τ\tau, 95%95\% na 3τ3\tau. Elk eerste-orde-overgangsverschijnsel is deze kromme, uitgerekt met τ\tau en verschoven met x0x_0 en xx_\infty.

Propositie 7.3 (RC- en RL-kringen)

Een condensator CC die via een weerstand RR vanuit een ideale bron EE wordt geladen voldoet aan RC ⁣duC/ ⁣dt+uC=ERC\,\dd u_C/\dd t + u_C = E: τ=RC\tau = RC, uC,=Eu_{C,\infty} = E. Een spoel LL in serie met RR over EE voldoet aan (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. Met de bron weggenomen (E0E \to 0) beschrijven dezelfde vergelijkingen de ontlading naar nul.

Bewijs. Maasregel E=Ri+uCE = Ri + u_C met i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t; maasregel E=Ri+L ⁣di/ ⁣dtE = Ri + L\,\dd i/\dd t.

Propositie 7.4 (Continuïteitsvoorwaarden)

De spanning over een condensator en de stroom door een spoel zijn continue functies van de tijd: op een schakelmoment t0t_0 geldt uC(t0+)=uC(t0)u_C(t_0^+) = u_C(t_0^-) en iL(t0+)=iL(t0)i_L(t_0^+) = i_L(t_0^-). Stromen door weerstanden en spanningen over spoelen mogen springen.

Bewijs. Een sprong van uCu_C zou een oneindige i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t vergen, een sprong van iLi_L een oneindige u=L ⁣di/ ⁣dtu = L\,\dd i/\dd t; beide zijn onmogelijk in een kring met eindige spanningen en stromen (evenzo kunnen de opgeslagen energieën 12CuC2\tfrac12 Cu_C^2 en 12Li2\tfrac12 Li^2 niet ogenblikkelijk veranderen zonder oneindig vermogen).

Methode 7.5 (Een eerste-orde-overgang oplossen)

  1. Schrijf na de schakelaar de maas- en knooppuntregels op, met de componentwetten; herleid tot één vergelijking in uCu_C of iLi_L en breng haar in de vorm τx+x=x\tau x' + x = x_\infty; lees τ\tau en xx_\infty af (het laatste is ook de eindtoestand die men vindt door CC door een open keten en LL door een draad te vervangen).
  2. Bepaal x0x_0 uit de continuïteit van uCu_C of iLi_L over het schakelmoment heen, met de kring ervóór.
  3. Schrijf x=x+(x0x)et/τx = x_\infty + (x_0 - x_\infty)\eu^{-t/\tau}; leid de overige grootheden af door differentiëren en met de maasregel.

Ziet de condensator een netwerk van weerstanden en bronnen, vervang dat netwerk dan eerst door zijn Thévenin-equivalent (ETh,RTh)(E_{\mathrm{Th}}, R_{\mathrm{Th}}) (Hoofdstuk 6): τ=RThC\tau = R_{\mathrm{Th}}C, u=EThu_\infty = E_{\mathrm{Th}}.

Propositie 7.6 (Energiebalans van de RC-lading)

Bij het laden van een condensator van 00 tot EE door een willekeurige weerstand RR levert de bron CE2CE^2, slaat de condensator 12CE2\tfrac12 CE^2 op en dissipeert de weerstand de andere 12CE2\tfrac12 CE^2 — wat RR ook is.

Bewijs. i=(E/R)et/τi = (E/R)\eu^{-t/\tau}; de bron levert 0Ei ⁣dt=E(E/R)τ=CE2\int_0^\infty Ei\,\dd t = E \cdot (E/R)\tau = CE^2; de weerstand neemt 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; het verschil is de opgeslagen 12CE2\tfrac12 CE^2. Een kleinere RR dissipeert dezelfde energie sneller.

Voorbeeld 7.7 (Een relaisspoel)

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}; de contacten trekken aan wanneer ii 1.5A1.5\,\mathrm{A} bereikt, dus op t=τln(10.75)=5.5mst = -\tau\ln(1 - 0.75) = 5.5\,\mathrm{ms}. Bij het openen van de kring moet ii van 2A2\,\mathrm{A} dalen; een vrijloopdiode over de spoel geeft de stroom een weg door RR alleen, τ=4ms\tau = 4\,\mathrm{ms}, in plaats van een vonk.

7.2 Systemen van de tweede orde

Definitie 7.8 (Canonieke tweede-ordevorm)

Een grootheid x(t)x(t) vormt een lineair systeem van de tweede orde wanneer

 ⁣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 ,

met ω0>0\omega_0 > 0 de eigenhoekfrequentie, Q>0Q > 0 de kwaliteitsfactor (dimensieloos) en xx_\infty de eindtoestand. Men schrijft ook ω0/Q=2ξω0\omega_0/Q = 2\xi\omega_0 met ξ=1/2Q\xi = 1/2Q de dempingsverhouding.

Propositie 7.9 (De RLC-kring in serie)

Een condensator CC, een spoel LL en een weerstand RR in serie over een bron EE: de spanning over de condensator voldoet aan de canonieke vorm met

ω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 .

Bewijs. Maasregel: E=L ⁣di/ ⁣dt+Ri+uCE = L\,\dd i/\dd t + Ri + u_C met i=C ⁣duC/ ⁣dti = C\,\dd u_C/\dd t: LCuC+RCuC+uC=ELC\,u_C'' + RC\,u_C' + u_C = E; deel door LCLC en identificeer ω02=1/LC\omega_0^2 = 1/LC, ω0/Q=R/L\omega_0/Q = R/L.

De twee gezichten van het tweede-ordesysteem: de RLC-kring in serie (_0 = 1/√LC, Q = √L/C/R) en het stelsel massa–veer–demper (_0 = √k/m, Q = √km/). Dezelfde vergelijking, dezelfde drie regimes.
De twee gezichten van het tweede-ordesysteem: de RLC-kring in serie (ω0=1/LC\omega_0 = 1/\sqrt{LC}, Q=L/C/RQ = \sqrt{L/C}/R) en het stelsel massa–veer–demper (ω0=k/m\omega_0 = \sqrt{k/m}, Q=km/αQ = \sqrt{km}/\alpha). Dezelfde vergelijking, dezelfde drie regimes.

Propositie 7.10 (De mechanische oscillator)

Een massa mm aan een veer met veerconstante kk met een viskeuze demper met coëfficiënt α\alpha (kracht αx˙-\alpha\dot x), over xx uit het evenwicht verplaatst, voldoet aan mx¨+αx˙+kx=0m\ddot x + \alpha\dot x + kx = 0: de canonieke vorm met

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

De overeenkomst xqx \leftrightarrow q, x˙i\dot x \leftrightarrow i, mLm \leftrightarrow L, αR\alpha \leftrightarrow R, k1/Ck \leftrightarrow 1/C beeldt elk resultaat van het ene systeem op het andere af.

Bewijs. De tweede wet van Newton (Hoofdstuk 12) met de veerkracht kx-kx en de dempingskracht; deel door mm.

Stelling 7.11 (De drie regimes)

De oplossingen van de homogene vergelijking x+(ω0/Q)x+ω02x=0x'' + (\omega_0/Q)x' + \omega_0^2 x = 0 worden beheerst door de wortels van r2+(ω0/Q)r+ω02=0r^2 + (\omega_0/Q)r + \omega_0^2 = 0, met discriminant Δ=ω02(1/Q24)\Delta = \omega_0^2(1/Q^2 - 4):

  • Q<12Q < \tfrac12, aperiodisch: twee reële negatieve wortels 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}, een terugkeer zonder trilling;
  • Q=12Q = \tfrac12, kritisch: een dubbele wortel r=ω0r = -\omega_0, x=(A+Bt)eω0tx = (A + Bt)\eu^{-\omega_0 t}, de snelste terugkeer zonder doorschot;
  • Q>12Q > \tfrac12, pseudoperiodiek: complexe wortels r=1/τ±iωr = -1/\tau \pm \iu\omega met

    τ=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),

    een gedempte trilling met pseudoperiode T=2π/ωT = 2\pi/\omega binnen de omhullende ±A2+B2et/τ\pm\sqrt{A^2 + B^2}\,\eu^{-t/\tau}.

De volledige oplossing bij een constant rechterlid is xx_\infty plus de homogene oplossing; AA en BB volgen uit x(0)x(0) en x(0)x'(0).

Bewijs. De karakteristieke vergelijking van een lineaire vergelijking met vaste coëfficiënten (wiskunde, lineaire vergelijkingen van de tweede orde): de wortels zijn r=ω0/2Q±12Δr = -\omega_0/2Q \pm \tfrac12\sqrt\Delta. Voor Q>12Q > \tfrac12 is Δ=iω041/Q2\sqrt\Delta = \iu\omega_0\sqrt{4 - 1/Q^2}, wat het reële deel ω0/2Q=1/τ-\omega_0/2Q = -1/\tau en het imaginaire deel ±ω011/4Q2\pm\omega_0\sqrt{1 - 1/4Q^2} geeft; de reële oplossingen zijn de combinaties van et/τcosωt\eu^{-t/\tau}\cos\omega t en et/τsinωt\eu^{-t/\tau}\sin\omega t.

Antwoord op een sprong vanuit rust van een tweede-ordesysteem, voor drie kwaliteitsfactoren. Kleine Q: een traag kruipen; Q = 1/2: de snelste terugkeer zonder doorschot; grote Q: naklank die ongeveer Q trillingen duurt.
Antwoord op een sprong vanuit rust van een tweede-ordesysteem, voor drie kwaliteitsfactoren. Kleine QQ: een traag kruipen; Q=12Q = \tfrac12: de snelste terugkeer zonder doorschot; grote QQ: naklank die ongeveer QQ trillingen duurt.

Propositie 7.12 (QQ aflezen op een pseudoperiodiek beeld)

In het pseudoperiodieke regime staan twee opeenvolgende maxima in de verhouding xn+1/xn=eT/τx_{n+1}/x_n = \eu^{-T/\tau}; het logaritmisch decrement

δ=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)

geeft QQ; het aantal trillingen dat zichtbaar is voordat de amplitude onder 5%5\% zakt is ongeveer QQ. Voor Q1Q \gg 1 is ωω0\omega \approx \omega_0 en T2π/ω0T \approx 2\pi/\omega_0.

Bewijs. De omhullende et/τ\eu^{-t/\tau} wordt per periode met eT/τ\eu^{-T/\tau} vermenigvuldigd; T/τ=(2π/ω)(ω0/2Q)T/\tau = (2\pi/\omega)(\omega_0/2Q). De omhullende bereikt e35%\eu^{-3} \approx 5\% op 3τ=6Q/ω0QT3\tau = 6Q/\omega_0 \approx Q\,T, dus na ongeveer QQ perioden.

Vrije pseudoperiodieke uitdemping, x = x_0 -t/ t: opeenvolgende maxima krimpen met de constante factor -T/, en (x_n/x_n+1) = T/ π/Q meet de kwaliteitsfactor van het beeld af.
Vrije pseudoperiodieke uitdemping, x=x0et/τcosωtx = x_0\eu^{-t/\tau}\cos\omega t: opeenvolgende maxima krimpen met de constante factor eT/τ\eu^{-T/\tau}, en ln(xn/xn+1)=T/τπ/Q\ln(x_n/x_{n+1}) = T/\tau \approx \pi/Q meet de kwaliteitsfactor van het beeld af.

Voorbeeld 7.13 (Een RLC-naklank)

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: pseudoperiodiek, τ=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}, drie zichtbare trillingen. De kritische weerstand is Rc=2L/C=632ΩR_c = 2\sqrt{L/C} = 632\,\Omega; daarboven is de ontlading aperiodisch.

Methode 7.14 (Een tweede-orde-overgang oplossen)

  1. Schrijf de vergelijking op; breng haar in canonieke vorm; lees ω0\omega_0, QQ en xx_\infty af.
  2. Bepaal het regime uit QQ en schrijf de algemene oplossing als xx_\infty plus de homogene oplossing.
  3. Bepaal x(0)x(0) en x(0)x'(0) uit de continuïteit van uCu_C en iLi_L (voor de RLC: uC(0+)=uC(0)u_C(0^+) = u_C(0^-) en uC(0+)=i(0)/Cu_C'(0^+) = i(0^-)/C); los de twee constanten op.

Voorbeeld 7.15 (Beginvoorwaarden)

De condensator uit het vorige voorbeeld is tot U0U_0 geladen en wordt op t=0t = 0 op LL en RR gesloten (geen bron): x=0x_\infty = 0, uC(0)=U0u_C(0) = U_0, uC(0)=i(0)/C=0u_C'(0) = i(0)/C = 0, want de spoelstroom was nul. Dus A=U0A = U_0 en 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}] — voor Q1Q \gg 1 eenvoudig U0et/τcosω0tU_0\eu^{-t/\tau}\cos\omega_0 t, en de stroom i=C ⁣duC/ ⁣dti = -C\,\dd u_C/\dd t piekt rond U0C/L=U0/(ω0L)U_0\sqrt{C/L} = U_0/(\omega_0 L).

Opmerking 7.16 (Energie)

De opgeslagen energie 12CuC2+12Li2\tfrac12 Cu_C^2 + \tfrac12 Li^2 (of 12kx2+12mx˙2\tfrac12 kx^2 + \tfrac12 m\dot x^2) neemt af met de snelheid Ri2Ri^2 (αx˙2\alpha\dot x^2): differentieer en gebruik de vergelijking. Per pseudoperiode krimpt de energie met de factor e2T/τ\eu^{-2T/\tau}; bij Q=3Q = 3 gaat ongeveer 88%88\% ervan verloren per trilling. Een hoge QQ betekent een traag lek ten opzichte van de trilling — een goede klok, een scherpe resonantie (Hoofdstuk 8) — en een lage QQ een snelle terugkeer: de juiste keuze voor een vering of een meternaald.

7.3 Oefeningen

Oefening 7.1

Een RC-kring: R=10kΩR = 10\,\mathrm{k}\Omega, C=47nFC = 47\,\mathrm{nF}, geladen vanuit EE. Bereken τ\tau, uC(τ)/Eu_C(\tau)/E en de tijd om 99%99\% van EE te bereiken.

Oplossing

Oplossing van Oefening 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}.

Oefening 7.2

Een RL-kring: L=20mHL = 20\,\mathrm{mH}, R=5.0ΩR = 5.0\,\Omega, E=10VE = 10\,\mathrm{V}. Bereken τ\tau, de eindstroom, de stroom op t=τt = \tau en de spanning over de spoel vlak na het sluiten van de schakelaar.

Oplossing

Oplossing van Oefening 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}; op 0+0^+ is i=0i = 0 (continuïteit), dus uL=ERi=10Vu_L = E - Ri = 10\,\mathrm{V}.

Oefening 7.3

Een RLC in serie: L=10mHL = 10\,\mathrm{mH}, C=100nFC = 100\,\mathrm{nF}, R=100ΩR = 100\,\Omega. Bereken ω0\omega_0, f0f_0, QQ; benoem het regime; geef τ\tau en de pseudoperiode.

Oplossing

Oplossing van Oefening 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: pseudoperiodiek; τ=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}.

Oefening 7.4

Welke weerstand geeft bij dezelfde LL en CC het kritische regime? Wat is QQ dan? Wat gebeurt er voor R=2kΩR = 2\,\mathrm{k}\Omega?

Oplossing

Oplossing van Oefening 7.4.

Rc=2L/C=632ΩR_c = 2\sqrt{L/C} = 632\,\Omega, Q=12Q = \tfrac12. Bij 2kΩ2\,\mathrm{k}\Omega is Q=0.16Q = 0.16: aperiodisch, een traag kruipen dat door de kleine wortel r+ω0Q=5×103s1\abs{r_+} \approx \omega_0 Q = 5 \times 10^{3}\,\mathrm{s}^{-1} wordt beheerst (tijdschaal 0.2ms0.2\,\mathrm{ms}, maar geen trilling).

Oefening 7.5 ★★

Een geladen condensator C=1.0µFC = 1.0\,\text{µ}\mathrm{F} ontlaadt zich over een onbekende RR; de spanning halveert elke 3.5ms3.5\,\mathrm{ms}. Toon aan dat de halveringstijd τln2\tau\ln 2 is en bepaal RR.

Oplossing

Oplossing van Oefening 7.5.

u=U0et/τu = U_0\eu^{-t/\tau} halveert wanneer 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.

Oefening 7.6 ★★

Toon door directe integratie aan dat het laden van een condensator van 00 tot EE via RR precies 12CE2\tfrac12 CE^2 in RR dissipeert, onafhankelijk van RR. Wat wordt daarmee het rendement van de overdracht, en hoe zou het beter kunnen?

Oplossing

Oplossing van Oefening 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, met geen RR in zicht. De bron levert CE2CE^2: rendement 50%50\% voor elke weerstand. Beter: laden via een spoel (de stroom is dan niet evenredig met het spanningsverschil) of in meerdere spanningsstappen — elke stap van E/nE/n verliest maar 12C(E/n)2\tfrac12 C(E/n)^2.

Oefening 7.7 ★★

Op een oscilloscoop toont een vrije RLC-uitdemping opeenvolgende maxima van 2.0V2.0\,\mathrm{V}, 1.5V1.5\,\mathrm{V}, 1.12V1.12\,\mathrm{V} met tussenpozen van 0.40ms0.40\,\mathrm{ms}. Leid QQ en ω0\omega_0 af, en met C=1.0µFC = 1.0\,\text{µ}\mathrm{F} de waarden van LL en RR.

Oplossing

Oplossing van Oefening 7.7.

δ=ln(2.0/1.5)=0.29\delta = \ln(2.0/1.5) = 0.29 (en ln(1.5/1.12)=0.29\ln(1.5/1.12) = 0.29: consistent); Qπ/δ=11Q \approx \pi/\delta = 11. T=0.40msT = 0.40\,\mathrm{ms}, dus ω02π/T=1.57×104rad/s\omega_0 \approx 2\pi/T = 1.57 \times 10^{4}\,\mathrm{rad}/\mathrm{s} (de correctie 1/8Q21/8Q^2 is 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.

Oefening 7.8 ★★

Een condensator geladen tot U0=100VU_0 = 100\,\mathrm{V} (C=10µFC = 10\,\text{µ}\mathrm{F}) wordt op een spoel L=1.0mHL = 1.0\,\mathrm{mH} met kleine weerstand geschakeld (Q1Q \gg 1). Schrijf uC(t)u_C(t) en i(t)i(t) op en schat de piekstroom. Waar zit de energie op het moment van de piek?

Oplossing

Oplossing van Oefening 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 (voor tτt \ll \tau), piek U0C/L=100×0.1=10AU_0\sqrt{C/L} = 100 \times 0.1 = 10\,\mathrm{A}, bereikt een kwart periode na het sluiten, wanneer uC=0u_C = 0: alle energie (12CU02=50mJ\tfrac12 CU_0^2 = 50\,\mathrm{mJ}) is dan magnetisch, 12Li2=50mJ\tfrac12 Li^2 = 50\,\mathrm{mJ}.

Oefening 7.9 ★★

Een kwart auto: m=400kgm = 400\,\mathrm{kg} op een veer k=2.0×104N/mk = 2.0 \times 10^{4}\,\mathrm{N}/\mathrm{m}. Bereken ω0\omega_0, f0f_0 en de kritische dempingscoëfficiënt. Geef met een demper α=2000Ns/m\alpha = 2000\,\mathrm{N}\,\mathrm{s}/\mathrm{m} de waarde van QQ, het regime, de pseudoperiode en de uitdooftijd τ\tau.

Oplossing

Oplossing van Oefening 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, pseudoperiodiek; 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}.

Oefening 7.10 ★★★

Een spoel L=0.50HL = 0.50\,\mathrm{H} voert I0=2.0AI_0 = 2.0\,\mathrm{A}; de bron wordt losgekoppeld en de spoel blijft op een ontlaadweerstand R=100ΩR = 100\,\Omega staan. Geef i(t)i(t), τ\tau, de piekspanning over de spoel en de gedissipeerde energie. Wat zou er zonder die weerstand gebeuren?

Oplossing

Oplossing van Oefening 7.10.

i=I0et/τi = I_0\eu^{-t/\tau}, τ=L/R=5.0ms\tau = L/R = 5.0\,\mathrm{ms}; piekspanning RI0=200VRI_0 = 200\,\mathrm{V} op t=0+t = 0^+; energie 12LI02=1.0J\tfrac12 LI_0^2 = 1.0\,\mathrm{J} gedissipeerd in RR. Zonder die weerstand heeft de stroom geen weg:  ⁣di/ ⁣dt\dd i/\dd t is enorm, de spoelspanning loopt op tot kilovolts, en er slaat een vonk over de schakelaar.

Oefening 7.11 ★★★

Een condensator C=1.0µFC = 1.0\,\text{µ}\mathrm{F} wordt vanuit E=12VE = 12\,\mathrm{V} geladen via R1=10kΩR_1 = 10\,\mathrm{k}\Omega, terwijl R2=20kΩR_2 = 20\,\mathrm{k}\Omega er overheen staat. Bepaal met een Thévenin-equivalent de tijdconstante en de eindspanning; schrijf uC(t)u_C(t) op vanaf uC(0)=0u_C(0) = 0.

Oplossing

Oplossing van Oefening 7.11.

Gezien vanuit 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.

Oefening 7.12 ★★★

Kritisch regime, sprong vanuit rust (x(0)=x(0)=0x(0) = x'(0) = 0): toon aan dat x=x[1(1+ω0t)eω0t]x = x_\infty[1 - (1 + \omega_0 t)\eu^{-\omega_0 t}], en bepaal numeriek de tijd om 95%95\% van xx_\infty te bereiken (in eenheden van 1/ω01/\omega_0). Vergelijk met de insteltijd van de omhullende voor Q=5Q = 5 en met die van de trage wortel voor Q=0.1Q = 0.1; concludeer waarom “kritisch” het doel van de ingenieur is.

Oplossing

Oplossing van Oefening 7.12.

Dubbele wortel ω0-\omega_0: x=x+(A+Bt)eω0tx = x_\infty + (A + Bt)\eu^{-\omega_0 t}; x(0)=0x(0) = 0 geeft A=xA = -x_\infty, x(0)=0x'(0) = 0 geeft B=ω0AB = \omega_0 A: x=x[1(1+ω0t)eω0t]x = x_\infty[1 - (1 + \omega_0 t)\eu^{-\omega_0 t}]. 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. Voor Q=5Q = 5 heeft de omhullende et/τ\eu^{-t/\tau} met τ=10/ω0\tau = 10/\omega_0 3τ=30/ω03\tau = 30/\omega_0 nodig; voor Q=0.1Q = 0.1 is de trage wortel r+ω0Q\abs{r_+} \approx \omega_0 Q en is 3/r+=30/ω03/\abs{r_+} = 30/\omega_0. Beide zijn zes keer trager: kritische demping is de snelste terugkeer zonder doorschot, en dat is wat een meter, een deurdranger of een vering wil.

7.4 Probleem: Een vering ontwerpen

Probleem 7.1

Weekendopgave — dezelfde differentiaalvergelijking op een labtafel en onder een auto: hoeveel keer een versleten schokdemper laat wippen, wat de ontwerper in plaats daarvan kiest, en hoe één getal de twee uit elkaar houdt

Deel I bestudeert een RLC-kring in serie, L=10mHL = 10\,\mathrm{mH}, C=100nFC = 100\,\mathrm{nF}; de delen II–IV een kwartautomodel: een massa m=350kgm = 350\,\mathrm{kg} rustend op een veer k=22kN/mk = 22\,\mathrm{kN}/\mathrm{m} met een schokdemper die αx˙-\alpha\dot x uitoefent. Het doorschot van een sprongantwoord in het pseudoperiodieke regime (gegeven): xmax/x1=exp(πξ/1ξ2)x_{\max}/x_\infty - 1 = \exp(-\pi\xi/\sqrt{1 - \xi^2}), ξ=1/2Q\xi = 1/2Q.

Deel I — De schakeling.

  1. Schrijf de maasregel op voor de RLC in serie over een bron EE, en de differentiaalvergelijking voor uCu_C.
  2. Breng haar in canonieke vorm; druk ω0\omega_0 en QQ uit.
  3. Bereken ω0\omega_0, f0f_0 en de kritische weerstand RcR_c.
  4. Schrijf de karakteristieke vergelijking en haar discriminant op; benoem de drie regimes in termen van QQ.
  5. Geef in het pseudoperiodieke geval τ\tau en ω\omega en de vorm van uC(t)u_C(t).
  6. Bereken voor R=100ΩR = 100\,\Omega de waarden van QQ, τ\tau, de pseudoperiode en het aantal zichtbare trillingen.
  7. De condensator, tot U0U_0 geladen, wordt op LL en RR geschakeld op t=0t = 0. Geef uC(0+)u_C(0^+) en  ⁣duC/ ⁣dt(0+)\dd u_C/\dd t(0^+), met de natuurkundige reden voor elk.

Deel II — De kwartauto.

  1. Meet xx vanaf de evenwichtsstand, schrijf de tweede wet van Newton voor de carrosserie op en toon aan dat het gewicht wegvalt.
  2. Breng de vergelijking in canonieke vorm; geef ω0\omega_0 en QQ in termen van mm, kk, α\alpha, en noem de elektromechanische overeenkomsten.
  3. Bereken ω0\omega_0, f0f_0 en de kritische coëfficiënt αc\alpha_c.
  4. De ontwerper monteert α=4000Ns/m\alpha = 4000\,\mathrm{N}\,\mathrm{s}/\mathrm{m}: bereken QQ en ξ\xi, en benoem het regime.
  5. Bereken τ\tau en de pseudoperiode.
  6. Bereken het doorschot na een sprong (een stoeprand): is de auto comfortabel?
  7. Welke weerstand RR zou de schakeling van deel I dezelfde QQ geven als deze vering?

Deel III — De versleten demper. Er is olie weggelekt: α=1200Ns/m\alpha = 1200\,\mathrm{N}\,\mathrm{s}/\mathrm{m}.

  1. Bereken QQ, de pseudoperiode en τ\tau.
  2. Met welke factor krimpen opeenvolgende maxima? Hoeveel keer wippen is zichtbaar voordat de amplitude onder 5%5\% zakt?
  3. De veertest: men duwt de hoek x0=5.0cmx_0 = 5.0\,\mathrm{cm} omlaag en laat los. Geef de beginvoorwaarden, de energie die in de veer is opgeslagen, en zeg waar die energie blijft.
  4. Schat de pieksnelheid van de carrosserie tijdens de eerste terugkeer (x0ω0\approx x_0\omega_0 voor Q1Q \gg 1) en de piekkracht in de demper.
  5. Waarom is precies kritische demping (Q=12Q = \tfrac12) niet de keuze van de ontwerper, hoewel zij geen doorschot geeft? (Denk aan de laatste millimeters van de terugkeer.)
  6. Bereken het doorschot voor de versleten demper en vergelijk met deel II.

Deel IV — QQ uit een beeld aflezen. Een versnellingsopnemer op de carrosserie registreert de vrije uitdemping.

  1. Voor de versleten auto zijn de opeenvolgende maxima 1.001.00, 0.250.25, 0.0620.062 (willekeurige eenheden). Bereken het logaritmisch decrement, leid ξ\xi exact af uit δ=2πξ/1ξ2\delta = 2\pi\xi/\sqrt{1 - \xi^2} en daarna QQ; vergelijk met deel III.
  2. Voor de gezonde auto is geen tweede maximum zichtbaar; alleen het doorschot, 3.8%3.8\%, is af te lezen. Keer de doorschotformule om om ξ\xi en QQ te vinden.
  3. Welk deel van de mechanische energie dissipeert de versleten demper per pseudoperiode?
  4. Een monteur zonder auto bij de hand bouwt de RLC van deel I met een RR zo gekozen dat Q=2.3Q = 2.3: welke RR, en waarom heeft het oscilloscoopbeeld, uitgezet tegen ω0t\omega_0 t, dezelfde vorm als dat van de auto?
  5. Vat samen: welk dimensieloos getal de vorm van elk tweede-orde-overgangsverschijnsel bepaalt, welke waarde een ontwerper van een vering nastreeft, en wat een bestuurder moet afleiden uit “meer dan één keer wippen”.
Oplossing

Oplossing van Probleem 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 met ω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 aperiodisch, Q=12Q = \tfrac12 kritisch, Q>12Q > \tfrac12 pseudoperiodiek.

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}, ongeveer drie trillingen.

7. uC(0+)=U0u_C(0^+) = U_0 (uCu_C continu: eindige stroom); uC(0+)=i(0+)/C=0u_C'(0^+) = i(0^+)/C = 0 (iLi_L continu en ervoor nul).

8. mx¨=k(x+xeq)αx˙+mgm\ddot x = -k(x + x_{\text{eq}}) - \alpha\dot x + mg met 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: pseudoperiodiek, net onder kritisch.

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. Doorschot exp(π×0.72/0.69)=exp(3.3)=3.8%\exp(-\pi \times 0.72/0.69) = \exp(-3.3) = 3.8\%: één kleine uitwijking voorbij het evenwicht, en klaar — comfortabel.

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: elke wip is een kwart van de vorige; 0.253=1.6%0.25^3 = 1.6\%: drie zichtbare wippen.

17. x(0)=x0x(0) = -x_0, x˙(0)=0\dot x(0) = 0; energie 12kx02=0.5×22000×0.0025=28J\tfrac12 kx_0^2 = 0.5 \times 22000 \times 0.0025 = 28\,\mathrm{J}, als warmte gedissipeerd in de olie van de demper.

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}; dempkracht αv1200×0.40=480N\alpha v \approx 1200 \times 0.40 = 480\,\mathrm{N}.

19. Bij Q=12Q = \tfrac12 verloopt de terugkeer als (1+ω0t)eω0t(1 + \omega_0 t)\eu^{-\omega_0 t}, die asymptotisch kruipt: de laatste millimeters duren lang. Een iets te weinig gedempt systeem (ξ0.7\xi \approx 0.7) schiet enkele procenten door maar komt sneller binnen een tolerantieband en blijft daar.

20. exp(π×0.22/0.976)=exp(0.70)=50%\exp(-\pi \times 0.22/0.976) = \exp(-0.70) = 50\%, tegen 3.8%3.8\%: de versleten auto wipt halverwege terug voorbij het evenwicht.

21. δ=ln4=1.39\delta = \ln 4 = 1.39 (twee keer, consistent); ξ=δ/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: de versleten waarde.

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

23. Energie \propto amplitude2^2: overblijvende fractie per periode 0.252=6%0.25^2 = 6\%, dus 94%94\% gedissipeerd per pseudoperiode.

24. R=316/2.3=137ΩR = 316/2.3 = 137\,\Omega. In de variabele ω0t\omega_0 t bevat de canonieke vergelijking alleen QQ: twee systemen met dezelfde QQ tekenen dezelfde kromme.

25. De kwaliteitsfactor QQ (of ξ\xi) alleen bepaalt de vorm; een ontwerper streeft naar Q0.7Q \approx 0.7 (ξ0.7\xi \approx 0.7); meer dan één zichtbare wip betekent QQ ruim boven 11: de demper is versleten.

Begrippen gedefinieerd in dit hoofdstuk

Bekijk alle 393 begrippen in de begrippenlijst