Geometriske stivhetsmatriser for både stive bjelker og stive trekantede plateelementer (TPE) har en iboende asymmetri som stammer fra underliggende delmatriser knyttet til nodale momenter i tre dimensjoner. Denne asymmetrien er et resultat av hvordan nodale momenter roterer i rommet og kan sammenlignes med den inducerte momentmatrisen kjent fra 3D bjelkeelementer. Denne egenskapen ved matriser er imidlertid begrenset til elementnivå, og for strukturen som helhet oppstår et symmetrisk stivhetssystem når betingelsene for likevekt i de roterte konfigurasjonene blir tilfredsstilt.

Det sentrale i denne sammenhengen er at de anti-symmetriske komponentene i stivhetsmatrisene, som kan uttrykkes gjennom permutasjonssymbolet og nodale momenter, kansellerer ut i summen når elementene møtes i et felles knutepunkt og likevektsbetingelsene oppfylles. Dette betyr at på strukturelt nivå vil totalstivhetsmatrisen være symmetrisk, selv om de individuelle elementmatrisene er asymmetriske. Denne symmetrien er essensiell for numerisk stabilitet og effektivitet i beregningsmetoder, spesielt innenfor ikke-lineær analyse.

For å transformere de nodale momentene og de anti-symmetriske matrisene fra lokale til globale koordinater brukes transformasjonsmatriser, som opprettholder de tensorielle relasjonene og sikrer korrekt sammenstilling i den globale strukturmatrisen. Dette innebærer at momentvektorene i lokale systemer først roteres til det globale koordinatsystemet, før deres bidrag summeres i hver knutepunkt. Likevektsbetingelsen uttrykkes ved at summen av momentene rundt alle globale akser i hvert knutepunkt er null.

Når det gjelder praktisk analyse av ikke-lineære problemer, benyttes den stive plateelementets geometriske stivhetsmatrise, som er kvalifisert for stive legemer, for å representere rotasjonseffekter nøyaktig. Denne kan kombineres med en elastisk stivhetsmatrise som representerer membran- og bøyevirkninger, typisk sammensatt av eksisterende elementformuleringer fra litteraturen. Slik kombinert gir den numeriske modellen mulighet til å simulere ikke-lineære og post-bucklingsfenomener i plater og skall med høy presisjon.

Videre anvendes en inkrementell-iterativ oppdateringsmetode av Lagrange-typen, hvor rotasjonseffekter håndteres fullt ut ved hvert trinn, mens de naturlige deformasjonene behandles lineært. Analysen deles opp i en prediktorfase, der strukturelle forskyvninger beregnes som svar på lastøkninger, og en korrektorfase som forbedrer løsningen. Denne metoden sikrer at rotasjonsvirkninger integreres sømløst i analyseprosessen, hvilket er avgjørende for å fange opp kompleks ikke-lineær oppførsel i strukturer med betydelige rotasjoner.

Det er avgjørende å forstå at selv om elementmatrisene isolert sett kan virke kompliserte og asymmetriske, oppnår strukturen som helhet en stabil og symmetrisk respons gjennom de fysiske likevektsbetingelsene ved knutepunktene. Dette fundamentet er viktig for å kunne konstruere pålitelige og effektive numeriske metoder som kan håndtere avanserte mekaniske fenomener i tredimensjonale plater og skall.

Det er videre vesentlig for leseren å innse betydningen av koordinattransformasjoners rolle i å sikre korrekt samspill mellom lokale elementbidrag og den globale strukturelle responsen. Det understreker nødvendigheten av en dyp forståelse av tensorielle relasjoner og symmetriegenskaper i stivhetsmatriser, noe som ofte er avgjørende for suksessfull implementering av ikke-lineære analysemetoder i avansert konstruksjonsteknikk.

Hvordan Beregne Stivhetsmatriser og Oppdatere Geometriske Egenskaper i Romrammeelementer

Stivhetsmatriser utgjør en fundamental del av beregningene som benyttes i strukturanalyser av romrammesystemer. Disse matrisene er nødvendige for å representere de elastiske egenskapene til materialer og hvordan de reagerer på ytre belastninger. I denne konteksten er det viktig å forstå hvordan integralmatriser benyttes for å beregne stivhetsmatrisene for forskjellige typer elementer i rammesystemer. En slik matrise, som ofte brukes i analysen, har følgende form:

1Kstvg=iv{nsg}{nth}Thdi\int 1 K_{stv} g = iv \{ ns g \} \{ nt h \} T h di

Her representerer indekseringene gg og hh gradene av interpolasjonsfunksjoner, mens superskriftene ss og tt betegner ordren av derivasjoner. Eksponenten vv representerer multipliseringsfaktoren. Disse integralmatriserne blir brukt i alle beregninger relatert til stivhetsmatriser gjennom boka, og ulike matriser kan benyttes avhengig av den spesifikke strukturen og analysekravene.

Stivhetsmatrisene som presenteres her er beregnet ved hjelp av forskjellige integrerte matriser. Eksempler på de mest brukte integralmatrisene er som følger:

[K02133]=[363301013363327][K021_{33}] = \begin{bmatrix} -3 & -6 & 3 & 3 \\ 0 & -1 & 0 & 1 \\ 33 & 6 & -33 & 27 \end{bmatrix}

Disse matrisene er nødvendige for å beskrive stivheten til et romrammeelement under belastning, og de brukes til å representere elastiske egenskaper i både lineære og ikke-lineære analyser. Matrisene som listes opp, er kun noen av de mest benyttede; imidlertid finnes det mange andre matriser som kan konstrueres ved å bruke en grunnleggende formel:

[Ktsvhg]T[Ktsv hg]^{T}

Matematisk sett gir dette oss et uttrykk for hvordan stivhetsmatrisene kan beregnes på en systematisk måte, og hvordan disse matriser kan tilpasses ulike typer strukturelle analyser.

Det er viktig å merke seg at de integralmatriser som benyttes, kan variere sterkt avhengig av typen rammeelement og de fysiske egenskapene til materialet som analyseres. Dette innebærer at de valgte interpolasjonsfunksjonene må være tilpasset for å reflektere de spesifikke kravene til hver enkelt struktur.

Når man ser på romrammestrukturer, blir analysen betydelig mer kompleks enn for andre, mer enkle strukturer, da hver ende av et rammeelement har tre roterende frihetsgrader. Tradisjonelt har det blitt antatt at rotasjonene ved endene av et rammeelement er små, noe som forenkler beregningene betydelig. Dette gjør at rotasjonsvinklene kan addere i små trinn for hvert analysesteg. Imidlertid gjelder ikke denne forenklingen når rotasjonene er store, da loven om kommutativitet ikke holder.

Når rotasjonene er store, må vi bruke teorien om endelige rotasjoner for å oppdatere enderotasjonene til hver rammeelement. Dette krever en betydelig mer komplisert beregningsmetode som tar hensyn til den ikke-kommutative naturen av store rotasjoner.

I slike tilfeller brukes de såkalte "konvokert koordinater", som beskriver hvordan referanseaksene ved hvert node av rammeelementet oppdateres når det er gjort en stor rotasjon. Dette innebærer at referanseaksene for node A i den opprinnelige konfigurasjonen C0C_0 kan representeres ved de ortogonale enhetsvektorene ξa0,ηa0,ζa0\xi_a^0, \eta_a^0, \zeta_a^0. Etterhvert som rotasjoner skjer, blir disse aksene oppdatert til de nye verdiene ved konfigureringen C2C_2, og rotasjonene beregnes ved hjelp av Euler’s rotasjonsformel:

r=cos(ϕ)r+sin(ϕ)(n×r)+(1cos(ϕ))(nr)nr' = \cos(\phi)r + \sin(\phi) (\mathbf{n} \times r) + (1 - \cos(\phi))(\mathbf{n} \cdot r) \mathbf{n}

Dette beskriver hvordan rotasjonene påvirker posisjonen til noden og orienteringen til aksene i rommet.

Denne teknikken for å oppdatere geometri og aksene til romrammeelementer er essensiell når man jobber med store deformasjonsanalyser, som er nødvendige for å forstå og forutsi atferden til strukturer under realistiske belastninger.

Når elementene i romrammestrukturen deformeres, er det viktig å oppdatere både elementaksene og nodalaksene for å sikre nøyaktige beregninger av de naturlige deformasjonene og kreftene i elementene. For å gjøre dette, introduseres flere sett med ortogonale akser for å følge deformasjonene til hvert element i strukturen.

Avslutningsvis, for å forstå den nøyaktige oppførselen til rammeelementene under belastning, er det viktig å beherske metodene for oppdatering av geometri og aksene til de enkelte elementene. Dette vil tillate mer presise og realistiske analyser, spesielt når det gjelder store deformasjonsanalyser der små rotasjoner ikke kan brukes som en forenkling.