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.
7.1 Systemen van de eerste orde
Definitie 7.1(Lineair systeem van de eerste orde)
Een grootheid x(t) voldoet aan een lineaire eerste-ordevergelijking met constante coëfficiënten wanneer
τdtdx+x=x∞,
met τ>0 de tijdconstante en x∞ een constante (de waarde die de bron oplegt). Haar eindtoestand is x=x∞; het overgangsverschijnsel is de weg daarheen vanaf de beginwaarde x(0)=x0.
Stelling 7.2(Oplossing)
De enige oplossing met x(0)=x0 is
x(t)=x∞+(x0−x∞)e−t/τ.
Het gat tot de eindtoestand krimpt met de factor e−1≈0.37 per τ: 63% van de weg na τ, 95% na 3τ, 99% na 4.6τ. De raaklijn in de oorsprong bereikt x∞ op t=τ.
Bewijs.y=x−x∞ voldoet aan τy′+y=0, waarvan de oplossingen y=Ke−t/τ zijn (de wiskunde over lineaire differentiaalvergelijkingen bewijst dat er geen andere zijn); K=x0−x∞. De helling in 0 is −(x0−x∞)/τ, die het gat in de tijd τ zou dichten. ∎
Het eerste-ordeantwoord op een sprong vanuit x0=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 x0 en x∞.
Propositie 7.3(RC- en RL-kringen)
Een condensator C die via een weerstand R vanuit een ideale bron E wordt geladen voldoet aan RCduC/dt+uC=E: τ=RC, uC,∞=E. Een spoel L in serie met R over E voldoet aan (L/R)di/dt+i=E/R: τ=L/R, i∞=E/R. Met de bron weggenomen (E→0) beschrijven dezelfde vergelijkingen de ontlading naar nul.
Bewijs. Maasregel E=Ri+uC met i=CduC/dt; maasregel E=Ri+Ldi/dt. ∎
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 t0 geldt uC(t0+)=uC(t0−) en iL(t0+)=iL(t0−). Stromen door weerstanden en spanningen over spoelen mogen springen.
Bewijs. Een sprong van uC zou een oneindige i=CduC/dt vergen, een sprong van iL een oneindige u=Ldi/dt; beide zijn onmogelijk in een kring met eindige spanningen en stromen (evenzo kunnen de opgeslagen energieën 21CuC2 en 21Li2 niet ogenblikkelijk veranderen zonder oneindig vermogen). ∎
Methode 7.5(Een eerste-orde-overgang oplossen)
Schrijf na de schakelaar de maas- en knooppuntregels op, met de componentwetten; herleid tot één vergelijking in uC of iL en breng haar in de vorm τx′+x=x∞; lees τ en x∞ af (het laatste is ook de eindtoestand die men vindt door C door een open keten en L door een draad te vervangen).
Bepaal x0 uit de continuïteit van uC of iL over het schakelmoment heen, met de kring ervóór.
Schrijf x=x∞+(x0−x∞)e−t/τ; 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) (Hoofdstuk 6): τ=RThC, u∞=ETh.
Propositie 7.6(Energiebalans van de RC-lading)
Bij het laden van een condensator van 0 tot E door een willekeurige weerstand R levert de bron CE2, slaat de condensator 21CE2 op en dissipeert de weerstand de andere 21CE2 — wat R ook is.
Bewijs.i=(E/R)e−t/τ; de bron levert ∫0∞Eidt=E⋅(E/R)τ=CE2; de weerstand neemt ∫0∞Ri2dt=(E2/R)∫0∞e−2t/τdt=(E2/R)(τ/2)=21CE2; het verschil is de opgeslagen 21CE2. Een kleinere R dissipeert dezelfde energie sneller. ∎
Voorbeeld 7.7(Een relaisspoel)
L=20mH, R=5.0Ω, E=10V: τ=4.0ms, i∞=2.0A; de contacten trekken aan wanneer i1.5A bereikt, dus op t=−τln(1−0.75)=5.5ms. Bij het openen van de kring moet i van 2A dalen; een vrijloopdiode over de spoel geeft de stroom een weg door R alleen, τ=4ms, in plaats van een vonk.
7.2 Systemen van de tweede orde
Definitie 7.8(Canonieke tweede-ordevorm)
Een grootheid x(t) vormt een lineair systeem van de tweede orde wanneer
dt2d2x+Qω0dtdx+ω02x=ω02x∞,
met ω0>0 de eigenhoekfrequentie, Q>0 de kwaliteitsfactor (dimensieloos) en x∞ de eindtoestand. Men schrijft ook ω0/Q=2ξω0 met ξ=1/2Q de dempingsverhouding.
Propositie 7.9(De RLC-kring in serie)
Een condensator C, een spoel L en een weerstand R in serie over een bron E: de spanning over de condensator voldoet aan de canonieke vorm met
ω0=LC1,Q=R1CL,uC,∞=E.
Bewijs. Maasregel: E=Ldi/dt+Ri+uC met i=CduC/dt: LCuC′′+RCuC′+uC=E; deel door LC en identificeer ω02=1/LC, ω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.
Propositie 7.10(De mechanische oscillator)
Een massa m aan een veer met veerconstante k met een viskeuze demper met coëfficiënt α (kracht −αx˙), over x uit het evenwicht verplaatst, voldoet aan mx¨+αx˙+kx=0: de canonieke vorm met
ω0=mk,Q=αkm.
De overeenkomst x↔q, x˙↔i, m↔L, α↔R, k↔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 en de dempingskracht; deel door m. ∎
Stelling 7.11(De drie regimes)
De oplossingen van de homogene vergelijking x′′+(ω0/Q)x′+ω02x=0 worden beheerst door de wortels van r2+(ω0/Q)r+ω02=0, met discriminant Δ=ω02(1/Q2−4):
Q<21, aperiodisch: twee reële negatieve wortels r±=−2Qω0(1∓1−4Q2), x=Aer+t+Ber−t, een terugkeer zonder trilling;
Q=21, kritisch: een dubbele wortel r=−ω0, x=(A+Bt)e−ω0t, de snelste terugkeer zonder doorschot;
Q>21, pseudoperiodiek: complexe wortels r=−1/τ±iω met
τ=ω02Q,ω=ω01−4Q21,x=e−t/τ(Acosωt+Bsinωt),
een gedempte trilling met pseudoperiodeT=2π/ω binnen de omhullende ±A2+B2e−t/τ.
De volledige oplossing bij een constant rechterlid is x∞ plus de homogene oplossing; A en B volgen uit x(0) en 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±21Δ. Voor Q>21 is Δ=iω04−1/Q2, wat het reële deel −ω0/2Q=−1/τ en het imaginaire deel ±ω01−1/4Q2 geeft; de reële oplossingen zijn de combinaties van e−t/τcosωt en e−t/τsinωt. ∎
Antwoord op een sprong vanuit rust van een tweede-ordesysteem, voor drie kwaliteitsfactoren. Kleine Q: een traag kruipen; Q=21: de snelste terugkeer zonder doorschot; grote Q: naklank die ongeveer Q trillingen duurt.
Propositie 7.12(Q aflezen op een pseudoperiodiek beeld)
geeft Q; het aantal trillingen dat zichtbaar is voordat de amplitude onder 5% zakt is ongeveer Q. Voor Q≫1 is ω≈ω0 en T≈2π/ω0.
Bewijs. De omhullende e−t/τ wordt per periode met e−T/τ vermenigvuldigd; T/τ=(2π/ω)(ω0/2Q). De omhullende bereikt e−3≈5% op 3τ=6Q/ω0≈QT, dus na ongeveer Q perioden. ∎
Vrije pseudoperiodieke uitdemping, x=x0e−t/τcosωt: opeenvolgende maxima krimpen met de constante factor e−T/τ, en ln(xn/xn+1)=T/τ≈π/Q meet de kwaliteitsfactor van het beeld af.
Voorbeeld 7.13(Een RLC-naklank)
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: pseudoperiodiek, τ=2Q/ω0=0.20ms, ω=0.987ω0, T=0.20ms, drie zichtbare trillingen. De kritische weerstand is Rc=2L/C=632Ω; daarboven is de ontlading aperiodisch.
Methode 7.14(Een tweede-orde-overgang oplossen)
Schrijf de vergelijking op; breng haar in canonieke vorm; lees ω0, Q en x∞ af.
Bepaal het regime uit Q en schrijf de algemene oplossing als x∞ plus de homogene oplossing.
Bepaal x(0) en x′(0) uit de continuïteit van uC en iL (voor de RLC: uC(0+)=uC(0−) en uC′(0+)=i(0−)/C); los de twee constanten op.
Voorbeeld 7.15(Beginvoorwaarden)
De condensator uit het vorige voorbeeld is tot U0 geladen en wordt op t=0 op L en R gesloten (geen bron): x∞=0, uC(0)=U0, uC′(0)=i(0)/C=0, want de spoelstroom was nul. Dus A=U0 en −A/τ+Bω=0, B=U0/(ωτ)=U0/4Q2−1: uC=U0e−t/τ[cosωt+sin(ωt)/4Q2−1] — voor Q≫1 eenvoudig U0e−t/τcosω0t, en de stroom i=−CduC/dt piekt rond U0C/L=U0/(ω0L).
Opmerking 7.16(Energie)
De opgeslagen energie 21CuC2+21Li2 (of 21kx2+21mx˙2) neemt af met de snelheid Ri2 (αx˙2): differentieer en gebruik de vergelijking. Per pseudoperiode krimpt de energie met de factor e−2T/τ; bij Q=3 gaat ongeveer 88% ervan verloren per trilling. Een hoge Q betekent een traag lek ten opzichte van de trilling — een goede klok, een scherpe resonantie (Hoofdstuk 8) — en een lage Q een snelle terugkeer: de juiste keuze voor een vering of een meternaald.
7.3 Oefeningen
Oefening 7.1★
Een RC-kring: R=10kΩ, C=47nF, geladen vanuit E. Bereken τ, uC(τ)/E en de tijd om 99% van E te bereiken.
Welke weerstand geeft bij dezelfde L en C het kritische regime? Wat is Q dan? Wat gebeurt er voor R=2kΩ?
Oplossing
Oplossing van Oefening 7.4.
Rc=2L/C=632Ω, Q=21. Bij 2kΩ is Q=0.16: aperiodisch, een traag kruipen dat door de kleine wortel ∣r+∣≈ω0Q=5×103s−1 wordt beheerst (tijdschaal 0.2ms, maar geen trilling).
Oefening 7.5★★
Een geladen condensator C=1.0µF ontlaadt zich over een onbekende R; de spanning halveert elke 3.5ms. Toon aan dat de halveringstijd τln2 is en bepaal R.
Toon door directe integratie aan dat het laden van een condensator van 0 tot E via R precies 21CE2 in R dissipeert, onafhankelijk van R. Wat wordt daarmee het rendement van de overdracht, en hoe zou het beter kunnen?
Oplossing
Oplossing van Oefening 7.6.
i=(E/R)e−t/τ: ∫0∞Ri2dt=(E2/R)(τ/2)=21CE2, met geen R in zicht. De bron levert CE2: rendement 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/n verliest maar 21C(E/n)2.
Oefening 7.7★★
Op een oscilloscoop toont een vrije RLC-uitdemping opeenvolgende maxima van 2.0V, 1.5V, 1.12V met tussenpozen van 0.40ms. Leid Q en ω0 af, en met C=1.0µF de waarden van L en R.
Oplossing
Oplossing van Oefening 7.7.
δ=ln(2.0/1.5)=0.29 (en ln(1.5/1.12)=0.29: consistent); Q≈π/δ=11. T=0.40ms, dus ω0≈2π/T=1.57×104rad/s (de correctie 1/8Q2 is 0.1%). L=1/(ω02C)=4.1mH; R=ω0L/Q=5.8Ω.
Oefening 7.8★★
Een condensator geladen tot U0=100V (C=10µF) wordt op een spoel L=1.0mH met kleine weerstand geschakeld (Q≫1). Schrijf uC(t) en i(t) op en schat de piekstroom. Waar zit de energie op het moment van de piek?
Oplossing
Oplossing van Oefening 7.8.
uC≈U0e−t/τcosω0t; i=−CduC/dt≈CU0ω0sinω0t (voor t≪τ), piek U0C/L=100×0.1=10A, bereikt een kwart periode na het sluiten, wanneer uC=0: alle energie (21CU02=50mJ) is dan magnetisch, 21Li2=50mJ.
Oefening 7.9★★
Een kwart auto: m=400kg op een veer k=2.0×104N/m. Bereken ω0, f0 en de kritische dempingscoëfficiënt. Geef met een demper α=2000Ns/m de waarde van Q, het regime, de pseudoperiode en de uitdooftijd τ.
Een spoel L=0.50H voert I0=2.0A; de bron wordt losgekoppeld en de spoel blijft op een ontlaadweerstand R=100Ω staan. Geef i(t), τ, de piekspanning over de spoel en de gedissipeerde energie. Wat zou er zonder die weerstand gebeuren?
Oplossing
Oplossing van Oefening 7.10.
i=I0e−t/τ, τ=L/R=5.0ms; piekspanning RI0=200V op t=0+; energie 21LI02=1.0J gedissipeerd in R. Zonder die weerstand heeft de stroom geen weg: di/dt is enorm, de spoelspanning loopt op tot kilovolts, en er slaat een vonk over de schakelaar.
Oefening 7.11★★★
Een condensator C=1.0µF wordt vanuit E=12V geladen via R1=10kΩ, terwijl R2=20kΩ er overheen staat. Bepaal met een Thévenin-equivalent de tijdconstante en de eindspanning; schrijf uC(t) op vanaf uC(0)=0.
Oplossing
Oplossing van Oefening 7.11.
Gezien vanuit C: ETh=ER2/(R1+R2)=8.0V, RTh=R1∥R2=6.67kΩ; τ=RThC=6.7ms; uC=8.0(1−e−t/τ) V.
Oefening 7.12★★★
Kritisch regime, sprong vanuit rust (x(0)=x′(0)=0): toon aan dat x=x∞[1−(1+ω0t)e−ω0t], en bepaal numeriek de tijd om 95% van x∞ te bereiken (in eenheden van 1/ω0). Vergelijk met de insteltijd van de omhullende voor Q=5 en met die van de trage wortel voor Q=0.1; concludeer waarom “kritisch” het doel van de ingenieur is.
Oplossing
Oplossing van Oefening 7.12.
Dubbele wortel −ω0: x=x∞+(A+Bt)e−ω0t; x(0)=0 geeft A=−x∞, x′(0)=0 geeft B=ω0A: x=x∞[1−(1+ω0t)e−ω0t]. 95%: (1+s)e−s=0.05, s≈4.7: t95=4.7/ω0. Voor Q=5 heeft de omhullende e−t/τ met τ=10/ω03τ=30/ω0 nodig; voor Q=0.1 is de trage wortel ∣r+∣≈ω0Q en is 3/∣r+∣=30/ω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=10mH, C=100nF; de delen II–IV een kwartautomodel: een massa m=350kg rustend op een veer k=22kN/m met een schokdemper die −αx˙ uitoefent. Het doorschot van een sprongantwoord in het pseudoperiodieke regime (gegeven): xmax/x∞−1=exp(−πξ/1−ξ2), ξ=1/2Q.
Deel I — De schakeling.
Schrijf de maasregel op voor de RLC in serie over een bron E, en de differentiaalvergelijking voor uC.
Breng haar in canonieke vorm; druk ω0 en Q uit.
Bereken ω0, f0 en de kritische weerstand Rc.
Schrijf de karakteristieke vergelijking en haar discriminant op; benoem de drie regimes in termen van Q.
Geef in het pseudoperiodieke geval τ en ω en de vorm van uC(t).
Bereken voor R=100Ω de waarden van Q, τ, de pseudoperiode en het aantal zichtbare trillingen.
De condensator, tot U0 geladen, wordt op L en R geschakeld op t=0. Geef uC(0+) en duC/dt(0+), met de natuurkundige reden voor elk.
Deel II — De kwartauto.
Meet x vanaf de evenwichtsstand, schrijf de tweede wet van Newton voor de carrosserie op en toon aan dat het gewicht wegvalt.
Breng de vergelijking in canonieke vorm; geef ω0 en Q in termen van m, k, α, en noem de elektromechanische overeenkomsten.
Bereken ω0, f0 en de kritische coëfficiënt αc.
De ontwerper monteert α=4000Ns/m: bereken Q en ξ, en benoem het regime.
Bereken τ en de pseudoperiode.
Bereken het doorschot na een sprong (een stoeprand): is de auto comfortabel?
Welke weerstand R zou de schakeling van deel I dezelfde Q geven als deze vering?
Deel III — De versleten demper. Er is olie weggelekt: α=1200Ns/m.
Bereken Q, de pseudoperiode en τ.
Met welke factor krimpen opeenvolgende maxima? Hoeveel keer wippen is zichtbaar voordat de amplitude onder 5% zakt?
De veertest: men duwt de hoek x0=5.0cm omlaag en laat los. Geef de beginvoorwaarden, de energie die in de veer is opgeslagen, en zeg waar die energie blijft.
Schat de pieksnelheid van de carrosserie tijdens de eerste terugkeer (≈x0ω0 voor Q≫1) en de piekkracht in de demper.
Waarom is precies kritische demping (Q=21) niet de keuze van de ontwerper, hoewel zij geen doorschot geeft? (Denk aan de laatste millimeters van de terugkeer.)
Bereken het doorschot voor de versleten demper en vergelijk met deel II.
Deel IV — Q uit een beeld aflezen. Een versnellingsopnemer op de carrosserie registreert de vrije uitdemping.
Voor de versleten auto zijn de opeenvolgende maxima 1.00, 0.25, 0.062 (willekeurige eenheden). Bereken het logaritmisch decrement, leid ξ exact af uit δ=2πξ/1−ξ2 en daarna Q; vergelijk met deel III.
Voor de gezonde auto is geen tweede maximum zichtbaar; alleen het doorschot, 3.8%, is af te lezen. Keer de doorschotformule om om ξ en Q te vinden.
Welk deel van de mechanische energie dissipeert de versleten demper per pseudoperiode?
Een monteur zonder auto bij de hand bouwt de RLC van deel I met een R zo gekozen dat Q=2.3: welke R, en waarom heeft het oscilloscoopbeeld, uitgezet tegen ω0t, dezelfde vorm als dat van de auto?
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”.
19. Bij Q=21 verloopt de terugkeer als (1+ω0t)e−ω0t, die asymptotisch kruipt: de laatste millimeters duren lang. Een iets te weinig gedempt systeem (ξ≈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%, tegen 3.8%: de versleten auto wipt halverwege terug voorbij het evenwicht.
21.δ=ln4=1.39 (twee keer, consistent); ξ=δ/4π2+δ2=1.39/6.43=0.215, Q=1/2ξ=2.3: de versleten waarde.
22.−ln0.038=3.27=πξ/1−ξ2, dus ξ=3.27/π2+3.272=0.72, Q=0.69.
23. Energie ∝amplitude2: overblijvende fractie per periode 0.252=6%, dus 94% gedissipeerd per pseudoperiode.
24.R=316/2.3=137Ω. In de variabele ω0t bevat de canonieke vergelijking alleen Q: twee systemen met dezelfde Q tekenen dezelfde kromme.
25. De kwaliteitsfactorQ (of ξ) alleen bepaalt de vorm; een ontwerper streeft naar Q≈0.7 (ξ≈0.7); meer dan één zichtbare wip betekent Q ruim boven 1: de demper is versleten.