De dynamica van lineaire systemen die worden aangespoord door fractale Gaussische ruis is een onderwerp dat dieper inzicht biedt in de reacties van structuren en systemen in omstandigheden die ver verwijderd zijn van klassieke ruismodellen, zoals witte ruis. Fractale Gaussische processen, ofwel Fractional Gaussian Noise (fGn), vertonen langdurige afhankelijkheid, wat betekent dat de afhankelijkheid van een bepaalde waarde van de ruis over de tijd niet vervalt zoals bij gewone witte ruis. Dit kan belangrijke implicaties hebben voor de manier waarop we het gedrag van systemen onder dergelijke omstandigheden modelleren en analyseren.

Wanneer we het gedrag van een lineair systeem onder de invloed van fractale Gaussische ruis bestuderen, moeten we begrijpen dat de ruis zelf geen gewone markovische eigenschappen heeft. Het betekent dat het systeem geen geheugenloos gedrag vertoont, wat typisch is voor veel klassieke stochastische processen. In plaats daarvan kunnen we een langdurige afhankelijkheid of ‘long-range dependence’ waarnemen, wat het systeem in staat stelt om de geschiedenis van de ruis over lange periodes heen te ‘onthouden’.

Laten we een systeem overwegen met meerdere vrijheidsgraden (DOF), beschreven door de vergelijking

MX¨+CX˙+KX=RWH(t),MẌ + CẊ + KX = RWH(t),

waar XX de verplaatsingsvector is, bestaande uit X1,X2,...,XnX_1, X_2, ..., X_n, en M,C,KM, C, K respectievelijk de massa-, demping- en stijfheidsmatrices zijn. De vector WH(t)WH(t) vertegenwoordigt de invoer van onafhankelijke fractale Gaussische ruisprocessen met verschillende Hurst-indexen H1,H2,...,HmH_1, H_2, ..., H_m, die aangeven hoe lang de afhankelijkheid in de tijd zich uitstrekt.

In een dergelijk systeem, wanneer we de spectrale analysemethode toepassen, krijgen we een matrix van vermogensdichtheden voor de verplaatsingen, gegeven door

SX(ω)=H(ω)SW(ω)HT(ω),SX(\omega) = H^*(\omega)SW(\omega)H^T(\omega),

waar SW(ω)SW(\omega) de vermogensspectrale dichtheid is van de excitatievector WH(t)WH(t), en H(ω)H(\omega) de frequentieresponsiematrix is van het systeem. Deze formules stellen ons in staat om de respons van het systeem in de frequentiedomein te analyseren, wat essentieel is voor het begrijpen van de reactie van systemen onder fractale Gaussische ruis.

Een belangrijk punt is dat de correlatie tussen paren van incrementele ruiscomponenten, zoals dBHidBHi en dBHjdBHj, geen invloed heeft op de dynamica van het systeem, omdat de ruisprocessen onafhankelijk zijn. Dit maakt de analyse eenvoudiger, omdat we de effectiviteit van de excitatie kunnen isoleren zonder ons zorgen te maken over interacties tussen verschillende componenten van de ruis.

Wanneer we het tijdsdomein beschouwen, kunnen we de oplossing van het stochastische systeem vinden door een geschikte vectorfunctie Y(t)=eAtX(t)Y(t) = e^{ -At}X(t) in te voeren en gebruik te maken van de fractale stochastische differentiaalregel. De oplossing is dan een integraal die de transientterm, die naar nul gaat in een dissipatief systeem na lange tijd, en de stationaire oplossing bevat, die het systeem in een steady-state beschrijft.

Het bovenstaande biedt inzicht in de wiskundige modellen die de basis vormen van de analyses van systemen die worden aangespoord door fractale Gaussische ruis. Zo kan men, bijvoorbeeld, de gemiddelde kwadratische waarden van de verplaatsing en snelheid berekenen, wat van cruciaal belang is voor het ontwerp en de beoordeling van structuren die aan dergelijke ruis zijn blootgesteld.

Naast de technische berekeningen zijn er enkele fundamentele inzichten die van belang zijn voor de lezer. Ten eerste is het belangrijk te begrijpen dat de lange-afhankelijke eigenschap van fractale Gaussische ruis niet alleen invloed heeft op de respons van het systeem, maar ook op de energieverdeling binnen het systeem. Dit kan het gedrag van de kinetische energie ten opzichte van de totale energie beïnvloeden, wat belangrijke gevolgen heeft voor het ontwerp van systemen die veerkrachtig moeten reageren op langdurige externe invloeden.

Ten tweede is het essentieel om het effect van de Hurst-indexen HH te overwegen. Deze index geeft niet alleen de mate van langdurige afhankelijkheid aan, maar beïnvloedt ook de energiedistributie en de structurele respons van het systeem. Het veranderen van HH heeft dus niet alleen gevolgen voor de statistische eigenschappen van de ruis, maar ook voor de mechanische prestaties van het systeem in de tijd.

In de praktijk kan men ook zien hoe de krachtigere aard van de ruis, zoals wanneer HH neigt naar 1, kan resulteren in een afname van de kinetische energie in het systeem. Dit duidt op het feit dat de excitaties steeds meer als statische belastingen gaan functioneren, hetgeen kan leiden tot minder dynamische reactie van het systeem, wat in veel technische toepassingen belangrijk is om te begrijpen.

Hoe werkt stochastisch gemiddeld gedrag in systemen met niet-lineaire demping en ruis?

Het stochastisch gemiddelde van systemen met niet-lineaire demping, zoals beschreven door de bijbehorende stochastische differentiaalvergelijkingen, biedt krachtige methoden voor het analyseren van complexe dynamische systemen. Deze technieken zijn van groot belang voor het begrijpen van de respons van systemen onder externe excitatie, zoals ruis of andere onvoorspelbare invloeden.

Het stochastisch gemiddelde wordt vaak toegepast op systemen die meerdere soorten excitatie ervaren, bijvoorbeeld een systeem dat zowel harmonische als willekeurige excitatie ondergaat. De voornaamste uitdaging bij het gebruik van deze benadering is het integreren van de effecten van zowel de deterministische als de stochastische componenten. Hierbij speelt de potentiaal van de oplossing een cruciale rol. Wanneer we werken met een systeem van twee gekoppelde variabelen xcx_c en xsx_s, die respectievelijk de respons van de systematische en stochastische delen van het systeem representeren, moeten we rekening houden met de compatibiliteitsvoorwaarden tussen de stochastische termen en de deterministische termen. Dit vereist het oplossen van een stelsel van algebraïsche vergelijkingen die voortkomen uit de stochastische differentiaalvergelijkingen.

Een interessant aspect van deze theorie is de manier waarop de variabelen zoals xcx_c en xsx_s kunnen worden gemanipuleerd om een oplossing te verkrijgen die de werkelijkheid van het systeem benadert. Wanneer we een exact afstemmingsscenario (γ=0\gamma = 0 en D=0D = 0) aannemen, komen we uit op een oplossing die de optimale toestand van het systeem beschrijft. Dit betekent dat de systeemfrequentie exact gelijk is aan de resonantiefrequentie ω0\omega_0, wat resulteert in een stabiel gedrag zonder significante stochastische fluctuaties. In dit geval kunnen we de oplossing voor de potentiële functie ϕ(xc,xs)\phi(x_c, x_s) verkrijgen, die de dynamica van het systeem beschrijft in termen van de coördinaten xcx_c en xsx_s, en een geschikte stochastische distributie.

Als de afstemming echter niet exact is, zoals vaak het geval is in praktische toepassingen, moet er een benaderende oplossing worden gevonden. Dit wordt bereikt door de functies hch_c en hsh_s te vervangen door aangepaste functies HcH_c en HsH_s, die de niet-lineaire termen in het systeem nauwkeuriger beschrijven. Het doel is om een systeem te creëren dat exact oplosbaar is, wat leidt tot een benaderende oplossing voor het oorspronkelijke systeem.

Een belangrijk gevolg van deze benaderingen is de behoefte om de residuele fout tussen de oorspronkelijke en de benaderde systemen zo klein mogelijk te maken. Dit wordt gedaan door de zogenaamde gewogen residuenmethode toe te passen, waarbij een geschikte gewogen functie wordt gekozen om de fouten in de oplossing te minimaliseren. De resulterende set lineaire algebraïsche vergelijkingen wordt numeriek opgelost, wat resulteert in een nauwkeurige benadering van de oplossing voor het systeem.

Daarnaast moeten we ook rekening houden met de rol van externe excitatie, zoals witte ruis of Poisson-ruis, die de dynamica van het systeem beïnvloedt. In dergelijke gevallen kan het systeem worden gemodelleerd door stochastische differentiaalvergelijkingen die de invloed van de ruis expliciet opnemen. De toepassing van Poisson-ruis biedt een ander niveau van complexiteit in de dynamische analyse, omdat de ruis niet meer als een continue proces wordt beschouwd, maar als een discrete, toevallige excitatie die de systeemrespons beïnvloedt.

In systemen die zowel deterministische als stochastische excitatie ervaren, zoals bij de niet-lineaire demping of in systemen die reageren op Poisson-ruis, kan de Fokker-Planck-vergelijking (FPK) worden toegepast om de stationaire kansdichtheidsfunctie (PDF) van de responsamplitude te verkrijgen. Deze kansdichtheidsfunctie beschrijft de verdeling van de systeemrespons over tijd, wat essentieel is voor het begrijpen van het langetermijngedrag van het systeem.

Wat belangrijk is om te begrijpen, is dat stochastisch gemiddelde benaderingen in staat zijn om de dynamica van complexe systemen op een pragmatische manier te modelleren, waarbij zowel de stochastische als de deterministische aspecten van het systeem effectief worden geïntegreerd. De uiteindelijke resultaten zijn niet altijd exact, maar ze bieden een bruikbare benadering voor het voorspellen van de systeemrespons onder verschillende condities van excitatie.

Hoe de Stationaire Kansenverdeling en Momentwaarden van Stochastische Systemen te Bepalen

In de studie van stochastische systemen, vooral die welke worden aangedreven door fractionele Gaussiaanse ruis, wordt de dynamica van de systemen vaak beschreven door middel van complete elliptische integralen van de tweede en eerste soort, E2(x)E_2(x) en K(x)K(x), respectievelijk. Deze functies spelen een cruciale rol bij het berekenen van de stationaire kansverdelingen en de momentwaarden van systemen, zoals geïllustreerd in de numerieke simulaties van de oorspronkelijke en geaverageerde systemen (vergelijkingen 4.450, 4.447 en 4.454).

De stationaire kansverdeling p(λ)p(\lambda) en de momentwaarden, zoals het gemiddelde E[Λ]E[\Lambda] en het gemiddelde kwadraat E[Λ2]E[\Lambda^2], kunnen worden berekend door numerieke methoden en simulaties die de stochastische karakteristieken van het systeem vastleggen. De resultaten van deze simulaties tonen aan dat er een sterke overeenkomst is tussen de theoretische en numerieke benaderingen, zoals weergegeven in de figuren 4.32 en 4.33. Het is van belang dat men de invloed van de parameter kk begrijpt, die een belangrijke rol speelt bij de vorming van de kansverdelingen p(x)p(x) en p(x˙)p(\dot{x}) voor de verplaatsing XX en de snelheid X˙\dot{X} van het systeem.

De stationaire kansverdeling p(x,x˙)p(x, \dot{x}), die de gezamenlijke kansverdeling van de verplaatsing en de snelheid beschrijft, kan worden berekend via de methode van stochastische averaging, zoals gedetailleerd beschreven in de vergelijkingen (4.447) en (4.455). Dit is van cruciaal belang voor het begrijpen van het gedrag van stochastische systemen, vooral bij niet-lineaire en visco-elastische systemen die onder invloed staan van beide additieve en multiplicatieve stochastische ruis.

Bij het berekenen van de marginaal kansverdelingen p(x)p(x) en p(x˙)p(\dot{x}), evenals de momentwaarden E[X2]E[X^2] en E[X˙2]E[\dot{X}^2], moeten we verder kijken dan alleen de numerieke oplossingen. De resultaten van de simulaties, gepresenteerd in figuren 4.34 en 4.35, tonen aan dat er enkele afwijkingen zijn in de momentwaarden voor de snelheid X˙\dot{X}, vooral bij een Hurst-index van 0.5. Dit suggereert dat de kansverdeling voor xx niet altijd een Gaussische verdeling volgt, wat het geval is in het lineaire systeem met k=0k = 0. Dit verschil kan worden toegeschreven aan de complexe dynamiek die ontstaat door de interactie van de stochastische ruis en de niet-lineaire krachten die op het systeem inwerken.

Bovendien zijn er variaties in de stationaire kansverdeling p(x)p(x) voor verschillende waarden van de parameter kk, zoals zichtbaar is in figuur 4.36. Dit bevestigt dat de aard van de kansverdeling sterk afhankelijk is van de systeemparameters en de dynamische eigenschappen van het systeem. Dergelijke resultaten zijn van essentieel belang voor het ontwerpen en analyseren van systemen die worden beïnvloed door stochastische ruis, vooral wanneer deze systemen een niet-lineaire respons vertonen.

Hoewel de stochastische benaderingen in deze context zeer nuttig zijn, is het belangrijk voor de lezer te begrijpen dat de gedetailleerde berekeningen van kansverdelingen en momentwaarden slechts een deel van het verhaal vormen. De numerieke simulaties en analytische benaderingen kunnen alleen zo accuraat zijn als de aannames die aan de basis liggen van het systeemmodel. De keuze van het stochastische model, de parameters en de methode van stochastische averaging hebben allemaal invloed op de uiteindelijke resultaten. In veel gevallen kunnen er belangrijke afwijkingen optreden wanneer de systeemparameters buiten het bereik vallen waarin de stochastische benaderingen valide zijn.

Voor een vollediger begrip is het ook van belang te beseffen dat de simulaties en de verkregen kansverdelingen in de praktijk kunnen worden beïnvloed door diverse externe factoren, zoals ruisbronnen, variaties in de systeemparameters, of zelfs numerieke onnauwkeurigheden in de simulatie zelf. Het is van cruciaal belang om deze factoren in overweging te nemen bij de interpretatie van de resultaten.

Hoe Invloed van Impuls en Elastische Wanden de Statische Respons van een Twee-DoF Systeem Beïnvloedt

In complexe trillingssystemen waar massa's in contact komen met elastische muren, speelt de interactie tussen de massa’s en de wanden een cruciale rol in het bepalen van de statistische distributies van de verplaatsingen en energieniveaus van het systeem. Bij het onderzoeken van dergelijke systemen met twee graden van vrijheid (2-DOF) en stochastische processen is de invloed van stijfheid van de wanden, de afstand tussen de massa en de wanden, evenals de mate van demping van groot belang voor de uiteindelijke dynamische respons van het systeem.

De onderzochte systemen bestaan uit twee massa’s (m1 en m2) en twee wanden (links en rechts), waarbij de parameters van de elasticiteit van de wanden en de afstanden tussen de massa’s en de wanden variëren. Wanneer de impact van massa m2 tegen de elastische muren toeneemt, blijkt uit de numerieke simulaties dat de kansdichtheidsfunctie (PDF) van de verplaatsing Q2 steeds verder afwijkt van een normale (Gaussian) verdeling. Dit gebeurt onder de volgende omstandigheden: (i) wanneer de stijfheid van de elastische wanden (Bl en Br) toeneemt, (ii) wanneer de afstand tussen de massa m2 en de wanden (δl en δr) afneemt, en (iii) wanneer de verhouding tussen de excitatie-intensiteit en de demping toeneemt. In dergelijke gevallen wordt de PDF van de verplaatsing Q2 breder en vertoont deze meer niet-Gaussiaanse kenmerken.

Omgekeerd, wanneer de invloed van de impact afneemt, nadert de PDF van de verplaatsing steeds meer een Gaussiaanse verdeling, zoals duidelijk zichtbaar is in de grafieken voor kleinere waarden van δ. Dit suggereert dat de systematische dampings- en stijfheidsinstellingen essentieel zijn bij het voorspellen van de dynamische respons van het systeem, vooral wanneer de elastische interactie tussen de massa en de wanden significant is.

Het gebruik van stochastische methoden zoals het gemiddelde van quasi-NI-systemen biedt inzicht in de lange-termijn dynamiek van dergelijke systemen, ondanks dat directe analytische oplossingen moeilijk te verkrijgen zijn. De stochastische gemiddelde methode, toegepast op quasi-niet-integrabele Hamiltoniaanse systemen, maakt het mogelijk om de stationaire oplossingen te benaderen door middel van stochastische differentiaalvergelijkingen. In dit kader kunnen belangrijke statistieken zoals de drift- en diffusiecoëfficiënten worden afgeleid, wat cruciaal is voor het begrijpen van de fluctuaties in de energie en de verplaatsingen van de massa’s.

De numerieke resultaten, verkregen door zowel de stochastische methode als Monte Carlo-simulaties, tonen aan dat wanneer de stijfheid van de wanden groot is, de excitatie-intensiteit sterk is, of de massa dicht bij de wanden is, de invloed van de impact sterk is, en de stochastische methode goede resultaten oplevert. Echter, wanneer de impact verwaarloosbaar is, zoals in systemen met lagere stijfheid of grotere afstanden tussen de massa’s en de wanden, neemt de nauwkeurigheid van de stochastische gemiddelde methode af, wat duidt op de noodzaak om alternatieve benaderingen te overwegen in dergelijke gevallen.

De berekeningen suggereren ook dat bij het verwaarlozen van de impact de dynamica van het systeem kan worden gemodelleerd met een lineaire benadering. Dit leidt tot een vereenvoudigde formulering van het systeem als een vectorvergelijking, waar de massa’s en de demping door matrices worden voorgesteld. Het gebruik van deze lineaire benadering maakt het mogelijk om de stabiliteit en respons van het systeem te bestuderen in gevallen waar de impact minimaal is, en biedt een eenvoudiger pad voor de analytische oplossing van de bewegingsvergelijkingen van het systeem.

Aanvullende overwegingen:

Het is belangrijk te begrijpen dat de invloed van de impact op de dynamische respons van het systeem afhankelijk is van verschillende parameters, zoals de stijfheid van de wanden, de afstand tussen de massa en de wanden, en de mate van demping. Deze factoren kunnen de resultaten van stochastische berekeningen aanzienlijk beïnvloeden, vooral wanneer de wanden zich dichter bij de massa bevinden of de stijfheid hoog is. Het nauwkeurig instellen van deze parameters is cruciaal voor het verkrijgen van realistische simulaties van het systeem.

Daarnaast moet men in gedachten houden dat stochastische benaderingen, hoewel zeer nuttig, bepaalde aannames en beperkingen met zich meebrengen. Bijvoorbeeld, in systemen waar de impact klein is, kan de stochastische methode onnauwkeurig zijn, en het kan noodzakelijk zijn om andere wiskundige technieken te gebruiken om de dynamica van het systeem beter te modelleren.

Hoe Kwasi-Integrabele Hamiltoniaanse Systemen Werken in de Stochastische Benadering

In de studie van dynamische systemen, vooral in de context van quasi-integrabele Hamiltoniaanse systemen, wordt vaak verondersteld dat de onderliggende systemen, beschreven door Hamiltoniaanse vergelijkingen, in staat zijn om onder bepaalde omstandigheden integrabel te zijn. Dit betekent dat er een rijke structuur bestaat die de evolutie van het systeem volledig kan beschrijven via een set van behoudswetten of constante grootheden. In gevallen van resonantie kunnen deze systemen nog verder geanalyseerd worden door gebruik te maken van stochastische benaderingen, die helpen bij het begrijpen van de invloed van externe verstoringen, zoals ruis of stochastische krachten, op het systeem.

Voor het geval van interne resonantie, zoals beschreven in de oorspronkelijke set van vergelijkingen, zijn er significante aanwijzingen dat het systeem kan worden vereenvoudigd door middel van het gebruik van geschikte gemiddelden, vooral in de tijdsafhankelijke en ruimtegeïntegreerde vorm. Dit biedt een diep inzicht in hoe de Hamiltoniaanse dynamica van een systeem kan worden geanalyseerd door gebruik te maken van gemiddelde fraccionaire stochastische differentiaalvergelijkingen (SDEs).

Een fundamenteel resultaat in deze context is dat voor een systeem dat als quasi-integrabel kan worden beschouwd, de dynamica van de systeemvariabelen kan worden gereduceerd tot een aantal langzaam veranderende vectorprocessen die de toestand van het systeem volledig beschrijven in termen van de actievariabelen en de resonantiefasen. Het is belangrijk om te begrijpen dat deze benaderingen gebruik maken van het principe van stochastische gemiddelde, waarbij de snel variërende componenten van het systeem worden gemiddeld over de tijd of de fase-ruimte.

Wanneer we kijken naar de specifieke effecten van externe excitatie, zoals de aanwezigheid van stochastische krachten die het systeem beïnvloeden, kunnen we de stochastische dynamica van de systeemvariabelen beschrijven door middel van een set van fraccionaire SDEs. Deze stochastische modellen kunnen gedetailleerd worden geanalyseerd door Monte Carlo-simulaties, waarmee we de stationaire kansdichtheidsfunctie van het systeem kunnen berekenen. Dit biedt waardevolle inzichten in de langetermijngedragingen van het systeem onder de invloed van de externe stochastische krachten.

In het geval van systemen die zich in een kwasi-integrabele en resonante toestand bevinden, kunnen we een nieuwe reeks SDEs afleiden die de effecten van zowel de interne resonantie als de externe excitatie beschrijven. De modelgrootheden kunnen met behulp van tijdgemiddelde benaderingen worden geëvalueerd, die de evolutie van de systeemvariabelen effectief in kaart brengen door de interacties tussen de verschillende resonante componenten van het systeem te integreren.

Deze aanpak is niet alleen theoretisch interessant, maar ook praktisch van belang voor het begrijpen van fysische systemen die onderhevig zijn aan resonantie, zoals in de studie van mechanische systemen, elektronische circuits, en zelfs biologische systemen waar complexe dynamica optreden door resonantie-effecten. Het is essentieel om te begrijpen dat de gereduceerde SDEs die uit de oorspronkelijke Hamiltoniaanse beschrijvingen worden afgeleid, de essentiële eigenschappen van de lange-termijndynamica van het systeem bewaren, zelfs in de aanwezigheid van stochastische verstoringen.

Het uiteindelijke doel is het verkrijgen van een gedetailleerd begrip van de stationaire verdeling van de systeemvariabelen, wat de kansverdeling van de fasen en de actiewaarden beschrijft. Dit helpt niet alleen om het systeem te begrijpen, maar ook om voorspellingen te doen over het gedrag van het systeem onder variërende externe invloeden.

Wanneer we kijken naar simulaties van dergelijke systemen, wordt het duidelijk dat de stochastische benaderingen die hier worden gepresenteerd, een krachtige tool vormen voor het simuleren van realistische fysische systemen. De simulatie van de fraccionaire SDEs, in vergelijking met de oorspronkelijke Hamiltoniaanse dynamica, levert vaak resultaten die goed overeenkomen, wat aantoont dat deze benaderingen robuust zijn voor het modelleren van dynamische systemen in complexe omgevingen.

In de praktijk is het ook belangrijk te realiseren dat de verwaarlozing van niet-resonante en snel variërende componenten, die vaak in traditionele modellen worden gemist, de nauwkeurigheid van voorspellingen kan beïnvloeden. Het gebruik van gemiddelde technieken biedt echter een solide basis voor het verkrijgen van resultaten die dicht bij de werkelijke observaties liggen.