I mekanikkens verden møter man ofte materialer som ikke bare har elastiske eller viskøse egenskaper, men som utviser begge deler samtidig. Slike materialer kalles visko-elastiske materialer. For å forstå hvordan disse materialene reagerer på påkjenninger, er det avgjørende å ha en dypere forståelse av de krefter som oppstår under deres deformasjon, og hvordan de påvirker systemets respons.

En av de viktigste aspektene ved visko-elastiske materialer er deres evne til å gjennomgå både deformasjoner og flytfenomener når de utsettes for eksterne belastninger. Dette fører til to karakteristiske fenomener: avslapning og krype. Avslapning skjer når stressen i materialet synker over tid ved konstant deformasjon, mens krype innebærer at deformasjonen øker når stressen holdes konstant. Disse fenomenene er essensielle for å beskrive det mekaniske responsen til materialet over tid.

For å beskrive hvordan visko-elastiske materialer reagerer på belastning, benyttes en constitutive lov, som typisk uttrykker forholdet mellom stress og strain. I tilfelle av lineære elastiske materialer er denne loven simpel: σ = Eε, hvor σ er stressen og ε er strainen. For Newtonske væsker, som har viskøse egenskaper, kan loven skrives som σ = ηε̇, hvor η er viskositeten og ε̇ er strainens tidsderiverte.

Imidlertid, for visko-elastiske materialer, er constitutive loven langt mer kompleks. I stedet for en enkel lineær sammenheng, uttrykkes lovene for visko-elastisitet gjennom både elastiske og viskøse komponenter som er sammenkoblet på en spesifikk måte. Disse lovene kan uttrykkes som en kombinasjon av stress og strain i differensialform. Et vanlig eksempel på en slik modell er Kelvin-Voigt-modellen, hvor elastiske og viskøse komponenter er kombinert i parallel, og som beskrives av ligningen ησ̇ + Eσ = Eηε̇.

En annen modell, Maxwell-modellen, kombinerer elastiske og viskøse komponenter i serie og har en annen constitutive lov: ηε̇ + Eε = σ. Det finnes også mer komplekse modeller som Burgers-modellen, som benytter flere elastiske og viskøse elementer, og beskriver materialers oppførsel på en mer detaljert måte.

En annen viktig parameter i studiet av visko-elastiske materialer er det såkalte avslapningmodul, G(t), og krypekompatibiliteten, J(t). Disse er funksjoner som beskriver materialets respons på påkjenningene over tid. G(t) representerer materialets evne til å redusere stress under konstant strain, mens J(t) beskriver hvordan strain utvikler seg over tid når stressen holdes konstant. Begge disse funksjonene er viktige for å karakterisere visko-elastiske materialer, da de gir innsikt i hvordan materialet vil oppføre seg under varierende belastninger.

I tillegg til de grunnleggende ligningene for visko-elastiske materialer, kan man bruke Laplace-transformasjon for å forenkle beregningene. Ved å transformere den differensielle constitutive loven til Laplace-domenet, får man en enklere måte å beregne både avslapningmodul og krypekompatibilitet på. Dette er spesielt nyttig i systemer der man trenger å analysere materialets respons på et dynamisk nivå.

Et viktig poeng å merke seg er at komponentmodellene som er nevnt, som Kelvin-Voigt, Maxwell og Burgers, ikke nødvendigvis representerer de fysiske strukturene til visko-elastiske materialer. De er matematiske representasjoner som forenkler materialets komplekse oppførsel, og de kan være mer eller mindre nøyaktige avhengig av materialet og dets spesifikke egenskaper. Det er viktig å bruke de riktige modellene for forskjellige materialer og eksperimentelt verifisere modellens prediksjoner.

Videre er det også avgjørende å forstå at hysterese, som kan oppstå i forbindelse med materialers respons på dynamiske belastninger, kan føre til at materialer viser uventede krefter, spesielt i systemer med genetiske effekter som de som beskrives i relaterte hysteresismodeller. I disse tilfellene vil materialets respons ikke bare være avhengig av den påførte belastningen, men også av tidligere tilstander, noe som kan føre til langvarige effekter som påvirker systemets ytelse.

Materialer som har både elastiske og viskøse egenskaper finner vi i mange praktiske anvendelser, som plast, gummi, betong i byggteknikk, og til og med i bioteknologi som blod og muslinger. Derfor er det viktig å ha en grundig forståelse av både de teoretiske og praktiske aspektene ved visko-elastisitet for å kunne forutsi og optimere materialers ytelse i forskjellige miljøer og under ulike forhold.

Hvordan Beregne Stasjonære PDF-er og Middelverdi for Quasi-Hamiltonianske Systemer

I et stasjonært system kan sannsynlighetsfordelingen (PDF) p(q1,q2,p1,p2)p(q_1, q_2, p_1, p_2) for kvasi-Hamiltonianske systemer beskrives gjennom forskjellige metoder, som f.eks. stokastisk gjennomsnitt. Ved å bruke Monte Carlo-simuleringer og stokastiske differensialligninger (SDE), kan både den stasjonære PDF-en og middelverdiene E[Q12]E[Q_1^2] og E[P12]E[P_1^2] beregnes for et system som oppfyller likningene som beskrives i (7.11). Disse beregningene kan presenteres som integrasjoner over de relevante variablene i systemet, for eksempel:

E[Q12]=q12p(q1,q2)dq1dq2E[Q_1^2] = \int_{ -\infty}^{\infty} \int_{ -\infty}^{\infty} q_1^2 p(q_1, q_2) dq_1 dq_2

og

E[Q22]=q22p(q1,q2)dq1dq2.E[Q_2^2] = \int_{ -\infty}^{\infty} \int_{ -\infty}^{\infty} q_2^2 p(q_1, q_2) dq_1 dq_2.

Slike beregninger gjør det mulig å analysere systemet i dybden, der stasjonære PDF-er kan studeres både analytisk og numerisk.

Simuleringer viser at de stasjonære PDF-ene for systemet kan være veldig like for både den originale dynamiske modellen og den gjennomsnittlige stokastiske modellen (basert på den fraksjonelle SDE-en). For eksempel, i Figurene 7.1 og 7.2, sammenlignes simuleringene fra den originale systemet med de som er oppnådd fra den forenklede stokastiske differensialligningen, hvor de estimerte middelverdiene for energien E[H]E[H] og E[H2]E[H^2] viser tilnærmet likhet.

Som det er vist i Figurene 7.3 og 7.4, kan Monte Carlo-simuleringer brukes til å generere konturplot av de stasjonære PDF-ene for de individuelle systemvariablene p(q1,q2)p(q_1, q_2) og middelverdiene av kvadrerte posisjoner og momenta E[Q12]E[Q_1^2] og E[Q22]E[Q_2^2]. Dette gir en visuell fremstilling av hvordan systemet oppfører seg under ulike parametere. Forenklingen ved å bruke stokastisk gjennomsnitt i stedet for å simulere det originale systemet direkte, fører til betydelig kortere beregningstid uten å miste nøyaktigheten i resultatene.

En av hovedfordelene med stokastisk gjennomsnitt er at den tillater en raskere beregning av systemets stasjonære tilstand. Dette er spesielt viktig for komplekse systemer med mange frie parametere, som de som involverer flere uavhengige integrerte oscillatorer eller forstyrrede systemer.

Når det gjelder de mer avanserte kvasi-integrerbare Hamiltonianske systemene, kan systemet beskrives som integrerbart og ikke-resonant, med n uavhengige første integraler i involusjon. Dette gir et mer komplekst rammeverk for å forstå hvordan systemets energiparametere H1,H2,...,HnH_1, H_2, ..., H_n utvikler seg over tid. I slike tilfeller kan de relevante stokastiske differensialligningene for systemet være av den generelle formen:

nHiPjdPj=ϵmmjkdt+ϵ1/2lgjldBHl(t)\sum_{n} \frac{\partial H_i}{\partial P_j} dP_j = -\epsilon \sum_{m} m_{jk} dt + \epsilon^{1/2} \sum_{l} g_{jl} dB_{Hl}(t)

der gjlg_{jl} representerer støykomponentene som påvirker systemets dynamikk. Slike systemer kan ytterligere forenkles ved å bruke spesifikke moment-beskrivelser eller angulære variabler som IiI_i og ϕi\phi_i for integrerbare systemer, som gjør det mulig å bruke tidsgjennomsnitt i stedet for romgjennomsnitt for å finne løsninger.

Når det gjelder simulering av de kvasi-Hamiltonianske systemene, må man være oppmerksom på at resultatene kan variere avhengig av hvordan de numeriske metodene er implementert. Bruken av simuleringer for å finne den stasjonære PDF-en for p(h)p(h) gjør at man kan finne en god tilnærming til den originale dynamikken for systemet, som kan være for kompleks til å løses eksakt.

Et annet viktig aspekt ved modelleringen er hvordan de stokastiske differensiallikningene for IiI_i og ϕi\phi_i kan brukes til å beskrive dynamikken i integrerbare Hamiltonianske systemer. Når de stasjonære PDF-ene er kjent, kan de brukes til å generere mer realistiske simuleringer for systemets atferd, og gjøre det lettere å sammenligne med eksperimentelle data.

Til slutt, en nøkkel for forståelsen av slike systemer er at det ikke er nok å bare stole på de analoge beregningene av den stasjonære tilstanden. For å få et nøyaktig bilde av systemets dynamikk og forstå hvordan det reagerer på ulike forstyrrelser, må man ha en grundig forståelse av både de teoretiske prinsippene som ligger bak modellene, og hvordan man effektivt kan implementere og validere numeriske metoder for å simulere disse prosessene.