Bezig met laden van [MathJax]/jax/output/CommonHTML/jax.js
Cookz
Artikelen: 0
Berichten: 31
Lid geworden op: wo 17 feb 2021, 09:54

dubbele resonantie frequentie

Stel je hebt volgend systeempje, een locomotief, een wagon en in de wagon een lading. Alles is met veren aan elkaar verbonden. Er is geen demping. (de veer tussen de wagon en de trein heeft een waarde k).

Afbeelding

Als ik de vergelijkingen opstel krijg ik:
m¨x1+kx1=F+kx2
m¨x2+2kx2=kx1+k2x3+k2x4
m2¨x3+kx3=k2x2+k2x4
m2¨x4+kx4=k2x2+k2x3
Of in matrix vorm:
M__¨x_+K__x_=F_
[m0000m0000m20000m2]¨x_+[kk00k2kk2k20k2kk20k2k2k]x_=[F000]
Ok en dan nu de resonantie frequenties berekenen:,
Door de kracht weg te laten en het te transformeren naar het frequentie domein:
ω2M__X_+K__X_=0
M__1K__X_=ω2X_
En dus door de eigenwaarde van de matrix links uit te rekenen krijg ik:
ω2=[0km3km3km]
Dat betekent dat er een "dubbele" resonantie frequentie is op 3km.
Mijn probleem is nu dat bij die eigenwaarde ook eigenvectoren horen.
Die eigenvectoren zijn de modevormen (hoe het systeem beweegt in resonantie),
maar welke modevorm gaat er dan voor komen?
Het lijkt me niet zo intuïtief hoe het systeem nu gaat bewegen op die frequentie.
Of heb ik ergens een fout gemaakt waardoor mijn oplossing niet klopt?
CoenCo
Technicus
Artikelen: 0
Berichten: 1.211
Lid geworden op: di 18 okt 2011, 00:17

Re: dubbele resonantie frequentie

Ziet er volgens mij netjes uit.
Laten we beginnen met het bepalen van die eigenvectoren, dan gaan we daarna eens kijken of ze optreden.
Lukt het je om de eigenvectoren te bepalen?
Weet je op intuitie/ervaring al wat de eigenvector is met hoekfrequentie 0?
Gebruikersavatar
Xilvo
Moderator
Artikelen: 0
Berichten: 11.197
Lid geworden op: vr 30 mar 2018, 16:51

Re: dubbele resonantie frequentie

Hier de eigenwaardes en eigenvectoren (in kolommen) voor k=1, m=1

[ 3.00000000e+00 1.00000000e+00 -2.48646608e-16 3.00000000e+00]

[[ -5.00000000e-01 -1.00000000e+00 -1.00000000e+00 -3.22580645e-02]
[ 1.00000000e+00 -3.75323885e-16 -1.00000000e+00 6.45161290e-02]
[ -5.00000000e-01 1.00000000e+00 -1.00000000e+00 -1.00000000e+00]
[ -5.00000000e-01 1.00000000e+00 -1.00000000e+00 9.35483871e-01]]

De eerste mode met ω2=3k/m is blijkbaar een trilling waarbij de wagon tegengesteld aan, en met twee maal zo grote amplitude, beweegt als locomotief en massa's in de wagon.

De mode met ω2=k/m is een trilling waarbij de wagon stilstaat en de massa's erin met gelijke amplitude maar in tegenfase met de locomotief trillen.
Cookz
Artikelen: 0
Berichten: 31
Lid geworden op: wo 17 feb 2021, 09:54

Re: dubbele resonantie frequentie

Het uitrekenen van de eigenvectoren was nu niet meteen mijn probleem,
maar dus zoals je bij Xilvo kan zien:
(laat me het even duidelijk opschrijven, de eigenvectoren zouden onafhankelijk moeten zijn van m en k)

Eigenwaarden:
ω2=[0km3km3km]

Eigenvectoren:
[V1_V2_V3_V4_]=[0.50.57740.37800.00950.500.75590.01890.50.57740.37800.69740.50.57740.37800.7164]
De eerste 2 vectoren interesseren me nu niet.
Die 2 laatste zijn veel interessanter. Dus dit is wat de computer terug geeft als je de eigenvectoren vraagt.
Maar omdat je nu zit met een dubbele eigenwaarde op 3km zijn er niet gewoon 2 eigenvectoren.
Maar elke combinatie αV3_+βV4_ is een eigenvector van de frequentie 3km. En dus de vraag: hoe weet je welke beweging er nu gaat voorkomen op die frequentie?
Gebruikersavatar
Xilvo
Moderator
Artikelen: 0
Berichten: 11.197
Lid geworden op: vr 30 mar 2018, 16:51

Re: dubbele resonantie frequentie

Klopt V4 wel?

Die houdt het zwaartepunt niet op z'n plaats dus is het op z'n minst een lineaire combinatie van een trilling en de translatie, met ω=0.
Cookz
Artikelen: 0
Berichten: 31
Lid geworden op: wo 17 feb 2021, 09:54

Re: dubbele resonantie frequentie

oei ja sorry. -0.7164 moet 0.7164 zijn
Gebruikersavatar
Xilvo
Moderator
Artikelen: 0
Berichten: 11.197
Lid geworden op: vr 30 mar 2018, 16:51

Re: dubbele resonantie frequentie

Dan klopt het wel.
Mijn V1 komt overeen met jouw V3, mijn V4 niet met jouw V4.
Jouw V4 is dus een lineaire combinatie van mijn V1 en V4, en omgekeerd, uiteraard.

Iedere lineaire combinatie van de eigenvectoren met ω2=3k/m is weer een eigenvector voor die frequentie.

Hoe de zaak precies gaat trillen hangt af van hoe het systeem aangestoten wordt.
Gebruikersavatar
Xilvo
Moderator
Artikelen: 0
Berichten: 11.197
Lid geworden op: vr 30 mar 2018, 16:51

Re: dubbele resonantie frequentie

Tenslotte, een eigenvector voor de ontaarde frequentie loodrecht op [-0,5 1 -0,5 -0,5] is [0 0 -1 1]
Er is een trillingswijze waarbij de massa's binnen de wagon gelijk bewegen en gelijk met de locomotief. De wagon beweegt daarmee in tegenfase, met een dubbele amplitude.
En er is een trillingswijze waarbij locomotief en wagon stil staan en de massa's binnen de wagon tegen elkaar in bewegen.
Cookz
Artikelen: 0
Berichten: 31
Lid geworden op: wo 17 feb 2021, 09:54

Re: dubbele resonantie frequentie

Xilvo schreef: do 24 jun 2021, 13:46 En er is een trillingswijze waarbij locomotief en wagon stil staan en de massa's binnen de wagon tegen elkaar in bewegen.
Dus eigenlijk kan je besluiten dat deze modevorm niet gaat optreden voor een exciterende kracht op de locomotief?
Want er moet minstens een beweging zijn tussen de locomotief en de wagon opdat de inwerkende kracht wordt doorgegeven. Of maak ik hier een fout?
Gebruikersavatar
Xilvo
Moderator
Artikelen: 0
Berichten: 11.197
Lid geworden op: vr 30 mar 2018, 16:51

Re: dubbele resonantie frequentie

Als, vanuit een situatie waarin alles in rust is, de wagon gaat bewegen (doordat de locomotief er aan trek of om welke andere reden ook), dan ondervinden de massa's in de wagon identieke krachten (de linker via de veer links, de rechter via de veer rechts). Ze krijgen een gelijke versnelling en zullen t.o.v. elkaar niet bewegen. Die mode waarin ze tegen elkaar in bewegen zal niet aangeslagen worden.

Terug naar “Klassieke mechanica”