Differentialekvationen för en systemrörelse med stokastiska parametriska och externa excitationer kan skrivas som

X¨+2ζω0X˙+ω02[1+ξ1(t)]X=ξ2(t),\ddot{X} + 2 \zeta \omega_0 \dot{X} + \omega_0^2 [1 + \xi_1(t)] X = \xi_2(t),

där ξ1(t)\xi_1(t) och ξ2(t)\xi_2(t) är två stationära bredbandsprocesser med spektrala densiteter Sij(ω)S_{ij}(\omega). Denna formel kan exempelvis beskriva första modläge hos en pelare utsatt för slumpmässiga axiella och transversala laster. Genom transformationen X=AcosθX = A \cos \theta, X˙=Aω0sinθ\dot{X} = -A \omega_0 \sin \theta, där θ=ω0t+ϕ\theta = \omega_0 t + \phi, härleds de långsamt varierande amplitud- och fasfunktionerna, vilka i närvaro av svag dämpning och svaga excitationer kan analyseras med stokastisk medelvärdesbildning.

Antaganden om korta korrelationstider för excitationerna i förhållande till systemets relaxationstid leder till att amplituden A(t)A(t) och fasen Θ(t)\Theta(t) är långsamt varierande processer. Medelvärdesbildning över en kvasi-period ger uttryck för den glättade driv- och diffusionskoefficienterna, som är funktioner av amplituden. Amplituden följer en Itô-stokastisk differentialekvation vars parametrar beror på spektrala värden vid de karakteristiska frekvenserna ω0\omega_0 och 2ω02\omega_0. Av denna analys framgår att systemets respons beror på den parametriska excitationens spektrala densitet vid dubbla systemets egenfrekvens, samt på den externa excitationens spektrala densitet vid systemets egenfrekvens. Korrelationsstrukturen mellan de två excitationerna har däremot ingen direkt påverkan.

Den stationära sannolikhetsfördelningen (PDF) för amplituden visar ett kraftigt beroende av dämpningen och spektrala tätheter. Om den parametriska excitationen saknas, reduceras PDF:en till en Rayleigh-fördelning för amplituden, vilket är ett klassiskt resultat för linjära oscillatorer med bredbandsstörning. Om däremot den externa excitationen saknas och den parametriska excitationen är närvarande utan tillräcklig dämpning, blir systemets respons ostabil och amplituden kan växa obegränsat.

Vidare kan metoden tillämpas på primär-sekundära system, där ett huvudsystem är direkt exciterat av bredbandsstörningar medan ett sekundärt system är kopplat men inte direkt exciterat. Den sekundära systemets dynamik påverkar huvudsystemet genom modifieringar av dess effektiva massa, dämpning och styvhet. Det sekundära systemets rörelse kan lösas i steady state och bidrar därigenom till en justerad ekvation för huvudsystemets rörelse med modifierade dynamiska egenskaper.

Det är avgörande att i sådana samverkande system undvika resonans mellan primär- och sekundärsystemets egenfrekvenser, eftersom detta kan leda till starka och potentiellt instabila vibrationer. Den stokastiska medelvärdesbildningen ger möjlighet att hantera dessa svaga icke-linjäriteter och svaga excitationer och utvinna statistiska egenskaper hos systemets respons.

I denna teori spelar sambandet mellan korrelations- och spektrala funktioner en central roll, eftersom det möjliggör övergången från tidsdomänsbeskrivning till frekvensdomänsanalys där de relevanta spektrala värdena vid ω0\omega_0 och 2ω02\omega_0 kan isoleras och användas för prediktion av systemets dynamiska beteende.

Det är viktigt att förstå att antagandet om korta korrelationstider och svaga excitationer ligger till grund för tillämpligheten av stokastisk medelvärdesbildning. Om dessa villkor inte uppfylls kan resultaten bli missvisande, och mer avancerade metoder eller numeriska simuleringar kan krävas för att korrekt beskriva systemets respons.

Vidare påverkar icke-linjär dämpning i huvudsystemet starkt den stationära statistiken för systemets rörelse, vilket kan ge upphov till mer komplexa responsmönster än i linjära system med endast linjär dämpning. Den modifierade frekvensen och dämpningen i närvaro av sekundärsystemet kräver särskild uppmärksamhet i systemdesign och analys, särskilt när det gäller strukturell hållbarhet och vibrationskontroll.

Det bör också noteras att den stokastiska modellen för parametrisk excitation är mer generell än en enkel harmonisk parametrisk svängning, och att dess effekt på systemet är mer subtil, beroende på den bredbandskarakteristik som beskrivs av spektrala tätheter snarare än enskilda frekvenser.

Slutligen är det av vikt att systemets parametrar och excitationernas spektrala egenskaper noggrant bestäms eller uppskattas, då detta direkt påverkar precisionen i prediktionerna av systemets statistiska respons. Detta är särskilt relevant i ingenjörstillämpningar där säkerhetsmarginaler och livslängdsberäkningar baseras på tillförlitliga stokastiska modeller.

Hur kan stokastisk genomsnittsmetod tillämpas på quasi-icke-integrerbara Hamiltoniansystem?

Inom dynamiken av Hamiltoniansystem, särskilt för system som inte är fullt integrerbara, är det en utmaning att hantera de komplexa fluktuationerna och de kaotiska rörelserna som uppstår när externa störningar, såsom vit Gaussisk brus, introduceras. Den stokastiska genomsnittsmetoden ger en kraftfull teknik för att approximera den stationära sannolikhetsfördelningen (PDF) för ett system som påverkas av dessa störningar.

För ett Hamiltoniansystem med ett stort värde av Hamiltonianen är rörelsen ofta kaotisk, vilket innebär att de klassiska metoderna för att beräkna rörelser blir otillräckliga. Men när man lägger till stokastiska element till ett quasi-icke-integrerbart system, till exempel vitt brus, får man ett mer komplext system som fortfarande kan analyseras med vissa förenklingar, såsom den stokastiska genomsnittsmetoden.

För att beräkna den stationära sannolikhetsfördelningen för ett sådant system, börjar man med att formulera den genomsnittsberäknade Itô-ekvationen, som ofta skrivs på en form där drift- och diffusionskoefficienterna avgörs. De här koefficienterna beror på Hamiltonianens värde och på systemets egenskaper, såsom energi, potentiell energi och de olika parametrarna som styr dynamiken. Metoden innebär att man använder sig av flera integraler som i sin tur gör det möjligt att approximera den slutgiltiga fördelningen av systemets parametrar.

I fall av icke-integrerbara Hamiltoniansystem, som exempelvis det två-dimensionella icke-linjära systemet som ges av ekvationerna (5.19) i texten, blir analysen mer utmanande när det gäller att beräkna de nödvändiga sannolikhetsfördelningarna. För sådana system kan Hamiltonianen ofta inte lösas exakt med konventionella metoder, men genom att använda den stokastiska genomsnittsmetoden kan man ändå erhålla en användbar approximation av systemets dynamik.

En central aspekt i denna metod är att den gör det möjligt att beräkna en approximativ stationär PDF för Hamiltonianen, vilket kan användas för att förutsäga systemets beteende på lång sikt. För de givna exempelparametrarna, som i systemet (5.19), kan resultat som ges av den stokastiska genomsnittsmetoden jämföras med Monte Carlo-simuleringar för att verifiera deras noggrannhet. Dessa simuleringar ger en numerisk lösning som kan användas för att bedöma systemets respons under olika förhållanden.

För mer komplexa system, till exempel med fler fria parametrar eller högre dimensionsräknare, uppstår ofta svårigheter i att fullfölja flera domänintegraler som är nödvändiga för att beräkna genomsnittsdrift och diffusionskoefficienter. För att hantera dessa svårigheter har en metod föreslagits som innebär en transformering av koordinater i elliptiska system, vilket gör det möjligt att konvertera ett (2n-1)-dimensionalt integralområde till ett mer hanterbart n-dimensionellt område. Detta gör beräkningarna betydligt enklare och mer effektiva, vilket i sin tur gör att stokastisk genomsnittsmetod kan tillämpas även på mer komplexa system med många frihetsgrader.

Det är även viktigt att förstå att trots de sofistikerade matematiska verktygen och metoderna som används för att approximera den stationära sannolikhetsfördelningen, är dessa resultat i grunden fortfarande förenklingar av den verkliga dynamiken hos systemen. Därför är det avgörande för läsaren att inse att även om den stokastiska genomsnittsmetoden ger mycket användbara approximationer, kan det finnas avvikelser från verkliga system i vissa extrema eller ovanliga förhållanden. Därför är det klokt att använda metoden som ett komplement till andra mer exakta simuleringar eller analytiska metoder för att få en heltäckande bild av systemets beteende.

Hur Quasi-Partiellt Integrerbara Hamiltonianska System Fungerar

I studier av Hamiltonianska system, särskilt i kvasi-partiellt integrerbara system, är det avgörande att förstå deras dynamiska beteende och hur de kan beskrivas genom stokastiska metoder. Dessa system, även om de inte är fullt integrerbara, visar ändå vissa reglerbundna beteenden på grund av de svaga interna resonanser som kan uppkomma mellan systemets olika frekvenser.

För ett system där det finns β (1 ≤ β ≤ r − 2) svaga interna resonanser, kan relationerna mellan de olika frekvenserna ωη beskrivas genom summan

η=1r1kvωη=Ovε\sum_{η=1}^{r-1} k_v ω_η = O_v ε

Här representerar k_v resonansfaktorerna och O_v är en konstant, medan ε är en liten parameter som styr systemets uppträdande vid olika tidsintervall. De resulterande variablerna, som ∑r−1 kvη ωη, definieras för att fånga de små variationerna av systemets faser, vilka inte är helt integrerbara men ändå uppvisar ett visst grad av regelbundenhet.

För att hantera dessa system introduceras en uppsättning av nya vinkelvariabler, ϕ_v, definierade som:

ϕv=η=1r1kvηϕ_v = ∑_{η=1}^{r-1} k_v η

Denna transformation till nya variabler gör det möjligt att uttrycka systemets rörelser på ett mer hanterbart sätt, särskilt när det gäller att lösa de stochastiska differentialekvationerna (SIDEs) som styr systemets beteende.

De sidesteppade ekvationerna som beskriver förändringen av variablerna ϕ_v kan härledas från de ursprungliga Hamiltonianska ekvationerna. Detta görs genom att använda de relationer som introduceras av ϕ_v, och det leder till en uppsättning av ekvationer som tar hänsyn till både de långsamma och snabba variationerna i systemets dynamik. De långsamma processerna som beskriver Iη och H_r (de totala energi- och impulsvariablerna) förändras gradvis över tid, medan de snabba processerna, som hänför sig till de generella koordinaterna Q_i och momenta P_i, varierar snabbare och måste hanteras med en annan metod, ofta genom en stokastisk medelvärdesoperation.

I det här sammanhanget, genom att använda stokastiska genomsnitt, konvergerar systemet mot en (r+β)-dimensionell vektor-Markov-process när ε → 0. Detta innebär att systemets beteende kan approximeras som en Markov-process där framtida tillstånd endast beror på det nuvarande tillståndet, inte på historiken, vilket förenklar beräkningarna och tillämpningen av ekvationerna.

En annan central aspekt i hanteringen av dessa system är trunceringen av de oändliga termerna i de stokastiska differentialekvationerna. Genom att ignorera högre ordningens termer (ε^u och högre) erhålls approximativa lösningar som kan ge användbara resultat för praktiska tillämpningar. Den resulterande uppsättningen av ekvationer, som beskrivs genom (6.280)-(6.282), innefattar endast de relevanta termerna som påverkar systemets dynamik på de tids- och skalnivåer som är av intresse.

Vidare kan den komplexa interaktionen mellan de långsamt och snabbt varierande processerna beskrivas genom så kallade "störningstermer", där Gη,l och Gr,l representerar störningarna som kommer från de snabba variablerna. Dessa störningar kan införlivas i den slutliga modellen, vilket gör att systemet fortfarande kan bevara sina långsiktiga egenskaper trots de externa påverkningarna.

Det är också viktigt att förstå hur de olika komponenterna i systemet interagerar genom resonanser. Dessa resonanser leder till att vissa delar av systemet blir mer känsliga för förändringar än andra, och att dessa förändringar kan ha långtgående effekter på hela systemets dynamik. Att korrekt modellera och förstå resonansfenomen är därför avgörande för att kunna förutsäga hur systemet kommer att utvecklas under olika betingelser.

Avslutningsvis, när man arbetar med sådana system är det viktigt att förstå att det finns ett behov av både analytiska och numeriska metoder för att kunna hantera de komplexa och ofta icke-linjära relationerna mellan variablerna. Att förenkla dessa relationer genom approximativa metoder, samtidigt som man bevarar de mest avgörande dynamiska egenskaperna, är nyckeln till att kunna förstå och förutsäga beteendet hos quasi-partiellt integrerbara Hamiltonianska system.

Hur stokastiska genomsnitt metoder tillämpas på kvasi-Hamiltoniska system och deras dynamik

För att analysera dynamiken i kvasi-Hamiltoniska system som är utsatta för stokastiska effekter är det nödvändigt att förstå de underliggande matematiska modellerna och tillvägagångssätten. Den stokastiska process som representeras genom Differential Equations med stochastisk påverkan (SIDEs) innebär en komplex interaktion mellan deterministiska och stokastiska komponenter, vilket gör att systemets beteende kan vara svårt att förutsäga utan att använda medelvärdesmetoder och approximationer.

I sådana system är det vanligen nödvändigt att applicera metoder för stokastisk medelvärdesberäkning för att förenkla de komplexa systemets dynamik. Dessa metoder gör det möjligt att förenkla systemets rörelse genom att reducera antalet variabler och eliminera snabbare skalan av systemets fluktuationer. Detta uppnås genom att medelvärdera den stokastiska komponenten i systemet, vilket leder till förenklade uttryck som kan hanteras matematiskt.

En typisk kvasi-Hamiltonisk modell som används inom denna kontext involverar komplexa system med flera frihetsgrader, som de system som beskrivs av dynamiska ekvationer för rörelse med både Gaussisk och Poisson-vitt brus. I dessa fall modelleras systemets rörelse av ett uppsättning differentialekvationer där varje variabel representerar en osjälvständig del av systemets dynamik. För att hantera de stokastiska påverkningarna används metoder för att göra systemet mer hanterbart genom att reducera komplexiteten i de stokastiska termerna.

För att effektivt analysera denna dynamik, är det avgörande att förstå hur de stokastiska medelvärdesmetoderna tillämpas på dessa system. Det handlar inte bara om att hantera de stokastiska fluktuationerna genom att använda Fokker-Planck-ekvationer (FPK), utan också om att hantera de korrelationer som kan uppstå mellan de olika variablerna i systemet. Stokastiska medelvärdesmetoder för kvasi-Hamiltoniska system är alltså inte bara en teknisk förenkling utan en nödvändig metod för att förstå det övergripande beteendet i systemet när det påverkas av externa stokastiska faktorer.

Med hjälp av stokastiska medelvärden kan man approximera lösningarna till ekvationerna för olika systemkonfigurationer, såsom det som beskriver ett 4-frihetsgrads-system utsatt för både Gaussiskt och Poisson-vitt brus. Dessa system kan karakteriseras av en uppsättning av stokastiska differentialekvationer som beskriver hur systemets positioner och rörelser beror på både deterministiska krafter och externa, stokastiska bruseffekter. Genom att tillämpa metodiken för stokastisk medelvärde, kan systemets rörelser och dynamik förenklas och därmed göras hanterbara trots komplexiteten hos de stokastiska påverkningarna.

En annan viktig aspekt är hur de stokastiska systemens korrelationer påverkar de slutgiltiga lösningarna. Eftersom de stokastiska krafter som verkar på varje frihetsgrad kan vara korrelerade, krävs en noggrann hantering av dessa korrelationer för att undvika felaktiga approximationer. Detta gör att den stokastiska medelvärdesmetoden blir en kraftfull teknik för att få en bättre förståelse av dynamiken i system som annars skulle vara svåra att analysera.

Slutligen är det viktigt att notera att de gränsvillkor som appliceras på de stokastiska differentialekvationerna spelar en avgörande roll i lösningen av systemet. Dessa gränsvillkor beror på systemets karaktäristika, vilket innebär att förståelsen för hur dessa gränsvillkor tillämpas kan ge djupare insikter i systemets stabilitet och dynamik. Ett korrekt hanterat gränsvillkor kan avsevärt förenkla lösningen och ge mer exakta förutsägelser om systemets beteende under olika förhållanden.