For å forstå hvordan stivhetsmatrisene for forskjellige ramme- og trestangselementer kan utledes, er det viktig å gå gjennom en detaljert fremgangsmåte, slik som den som brukes for planar rammeelementer. I denne sammenhengen vil vi fokusere på hvordan de grunnleggende prinsippene for virtualt arbeid og de relevante integrasjonsprosedyrene fører til en forståelse av de matematiske modellene som brukes i strukturanalyse.
La (x, y) være et sett med ortogonale kartesiske koordinater, hvor x-aksen representerer sentrumaksen til det planare rammeelementet. I et slikt element er det kun den aksiale strekkingen som er betydelig, mens alle andre komponenter av strekkingen kan settes til null. Dette forenkler ligningene, og vi kan dermed bruke den grunnleggende uttrykket for virtualt arbeid:
hvor er Youngs modul for materialet, og integrasjonen tas over hele volumet til elementet. For et vilkårlig punkt N i elementet, kan aksial strekk relateres til aksial forskyvning ved å bruke den partielle deriverte:
Ved å anta Bernoulli–Euler-hypotesen, som sier at et tverrsnitt som er plan og normalt til den langsgående aksen av balken forblir plan og normalt etter deformasjon, kan vi relatere forskyvningene og til forskyvningene og til sentrum C i balkens tverrsnitt. Dette gir uttrykkene:
hvor primen () betegner den deriverte med hensyn til koordinaten x. Ved å erstatte i uttrykket for aksial strekk , får vi en ligning som uttrykker strekk som en funksjon av tverrdeformasjonene:
Neste trinn innebærer å erstatte denne uttrykket for i den opprinnelige ligningen for virtualt arbeid, noe som gir oss en generalisert form for arbeidsekvasjonen:
hvor er lengden på elementet, er tverrsnittsarealet, og er treghetsmomentet for tverrsnittet rundt z-aksen. Denne ligningen reflekterer det faktum at både aksial og bøyningsdeformasjoner spiller en rolle i elementets totale stivhet.
Et viktig aspekt ved utledningen av virtualt arbeid i denne sammenhengen er betingelsen for ortogonalitet i de sentrale koordinatene, som er:
Videre må vi beregne det eksterne virtuelle arbeidet , som kan deles opp i to komponenter: og , som representerer de eksterne kreftene som virker på elementene ved endene A og B. Ved å anta at overfladetrykkene og er konsentrert ved endene, kan vi beregne de eksterne kreftene ved å integrere over tverrsnittet til elementene på sluttene A og B:
Ved å bruke ligningene for forskyvningene og kan vi uttrykke det eksterne arbeidet i form av nodelige krefter. Den totale virtuelle arbeidsmengden for elementet er summen av de to komponentene:
Til slutt får vi den endelige uttrykksformen for virtualt arbeid i et planlagt rammeelement:
hvor er forskyvningsvektoren, og er kraftvektoren som er knyttet til nodene i elementet. Denne ligningen er grunnlaget for å etablere stivhetsmatrisene for rammeelementene.
Det er også viktig å merke seg at konvensjonen som benyttes for å bestemme tegnene på kreftene i elementet kan variere. I denne sammenhengen anses et direkte trykk som positivt når det virker mot høyre eller oppover, mens et bøyningsmoment er positivt når det virker mot klokken. Dette er en forskjell fra den vanlige mekanikk-konvensjonen, hvor aksialkraften er positiv når den strekker ut elementet.
Det bør også understrekes at en finitt elementformulering basert på virtualt arbeid er like rasjonell som en basert på de styrende differensialligningene for problemet som vurderes. Selv om denne formaliseringen kan utføres uten nødvendigvis å ha kunnskap om de tilknyttede differensialligningene, gir en forhåndsvisning av de relevante differensialligningene nyttige ledetråder for valg av interpolasjonsfunksjoner og nodal grader av frihet i finitte elementmodellen.
Hvordan beregne krefter i romlige truss-strukturer med ikke-lineær analyse
I en truss-struktur er krefter og deformasjoner som følge av belastningene på elementene nært knyttet til både de geometriske egenskapene og materialets elastiske egenskaper. For å forstå hvordan kreftene på de enkelte elementene i en romlig truss-struktur kan beregnes, er det nødvendig å benytte en kombinasjon av stivhetsmatriser som tar hensyn til både strekking og stivhet i den geometriske formen til strukturen.
For et trusselement, representert i et tredimensjonalt rom, består kreftene på punktene C1 og C2 av tre direkte krefter på hvert punkt. For C1, som illustrert i ligningene (4.61) og (4.62), er kreftene i x, y og z-retningene på de to endene av elementet uttrykt som og på henholdsvis de to endene av elementet.
Det er viktig å merke seg at for et trusselement i likevekt ved C1, vil alle de tverrgående skjærkreftene forsvinne, det vil si at . De aksiale kreftene på de to endene av stangen er like i størrelse, men motsatt i tegn, altså .
Videre, for å beskrive de nødvendige stivhetsmatrisene som brukes til å analysere slike elementer, er det viktig å forstå deres struktur. I romlige trusselementer er stivhetsmatrisene som [ke], [kg], [s1], [s2], og [s3] avgjørende for korrekt modellering av både elastisk og geometrisk stivhet. Disse matrisene har dimensjonene 6x6 for trusselementet, og deres utledning er en avgjørende del av den numeriske analysen.
For eksempel kan den elastiske stivhetsmatrisen uttrykkes som en matrise som inkluderer materialets elastisitetsmodul og elementenes geometri. Den geometriske stivhetsmatrisen beskriver elementenes stivhet som følge av deres geometriske konfigurasjon, og matrisene [s1], [s2], og [s3] tar hensyn til de høyere ordens stivhetene som kan være relevante i mer komplekse belastningssituasjoner.
I det nåværende tredimensjonale tilfellet kan man derivere de nødvendige ligningene for de elementære kreftene som følger av de forskjellige stivhetsmatrisene. For eksempel kan , , og de høyere ordens matrisene uttrykkes på en måte som bevarer både elastisk stivhet og geometriske effekter under belastning. Dette kan resultere i en endelig kraftuttrykk som er viktig for å forstå hvordan strukturen oppfører seg under belastning.
En viktig betraktning i denne prosessen er hvordan deformasjoner og strekking påvirker de resulterende kreftene. Når deformasjonene er små, kan den constitutive loven for materialet (som representerer aksial deformasjon) brukes til å forenkle beregningene. Dette kan redusere kompleksiteten ved å utelate visse termer som påvirker de høyere ordens deformasjonene.
Derimot, når de geometriske og materialmessige deformasjonene blir store, eller når elementene undergår betydelige rotasjoner, kan de forenklede uttrykkene som ikke tar hensyn til rotasjonsbevegelsene (som for eksempel ligning (4.80)) være utilstrekkelige. Det er derfor nødvendig å bruke en mer kompleks tilnærming for å ivareta de nøyaktige effektene av rotasjoner, spesielt når strukturen er utsatt for store deformasjoner eller rotasjoner i post-buklingstilstanden.
Når man ser på trusselementers oppførsel under store belastninger, bør en viktig forståelse være at elastisiteten alene kanskje ikke er tilstrekkelig for å modellere oppførselen til et trusselement under ikke-lineære forhold. Effekten av strekning og stivhet i strukturen, sammen med de geometriske endringene som skjer når strukturen deformeres, spiller en avgjørende rolle. For eksempel kan stivhetsmatrisene som er nevnt ovenfor, inkludert både elastisk stivhet og geometrisk stivhet, kombineres for å gi en mer nøyaktig beskrivelse av systemets respons.
Når det gjelder selve prosedyren for krefterberegning, er det også viktig å merke seg at analysen for romlige trusselementer kan gjennomføres ved hjelp av inkrementelle iterasjonsteknikker. For å sikre at den ikke-lineære analysen gir nøyaktige resultater, er det den korrigerende fasen av beregningen som er avgjørende. Denne fasen tar sikte på å finne de eksakte kreftene ved hjelp av matriseuttrykkene som er avledet gjennom de inkrementelle metodene, og gir en løsning som er i samsvar med de fysiske prinsippene for kontinuerlig mekanikk.
Det er også viktig å forstå at små og store rotasjoner i trusselementer kan føre til signifikante forskjeller i hvordan kreftene oppstår og balanseres. Dette er særlig relevant når man ser på trusselementer som er utsatt for store belastninger eller når rotasjoner blir tilstrekkelig store til å påvirke strukturell stabilitet. Dette gjør det nødvendig å tilpasse analysen ved å inkludere effektene av de forskjellige stivhetsmatrisene, som kan endres avhengig av rotasjonene i systemet.
Hvordan forstå prinsippet om virtuelle forskyvninger og dets rolle i analyse av ikke-lineære rammestrukturer?
Prinsippet om virtuelle forskyvninger er en grunnleggende mekanisk metode som muliggjør analyse av strukturer under komplekse belastninger, særlig i ikke-lineære rammestrukturer hvor deformasjoner og krefter påvirker hverandre gjensidig. Dette prinsippet bygger på konseptet at for en struktur i likevekt, vil det virtuelle arbeidet utført av indre krefter under en vilkårlig, infinitesimal, tenkt forskyvning være lik det virtuelle arbeidet utført av ytre krefter under samme forskyvning. Det innebærer en balanse mellom indre og ytre krefter uttrykt gjennom energimetoder, noe som åpner for en systematisk formulering av problemstillingen som kan løses numerisk.
I analysen av ikke-lineære rammestrukturer brukes prinsippet som utgangspunkt for å formulere differensiallikninger og endelig element-metode, der strukturen deles inn i mindre elementer med tilhørende stivhetsmatriser. Den inkrementelle Lagrange-formuleringen, både total og oppdatert variant, er sentrale i dette arbeidet for å håndtere geometri- og materialnonlineariteter. Total Lagrange-formulering refererer til analyse basert på den opprinnelige uforandrede konfigurasjonen av strukturen, mens oppdatert Lagrange-formulering benytter den siste kjente deformerede tilstanden som referanse for neste beregningssteg, noe som gir bedre tilpasning til store deformasjoner.
Lineær analyse utgjør et uunnværlig fundament for å forstå og kunne gjennomføre ikke-lineære beregninger. Gjennom derivering av elementstivhetsmatriser for forskjellige ramme- og buelementer etableres et grunnlag som sikrer både fysisk realisme og numerisk stabilitet. Det er essensielt å teste elementenes kvalitet gjennom kriterier som konvergens, patch-test og egenverdianalyse for å sikre at de grunnleggende prinsippene om bevegelsesfrihet og stivhetsforhold ivaretas korrekt. Videre sikrer rigid kroppsbevegelsestest at elementer oppfører seg fysisk riktig ved ren rigid bevegelse, uten kunstige deformasjoner eller feil i kraftbalansen.
Kinematiske hypoteser som legges til grunn ved formuleringen av elementer må være nøye vurdert for å ivareta balanse mellom elastiske og geometriske bidrag til stivheten. Dette innebærer også en forståelse av hvordan elementkrefter kan tolkes som følge av både strekking og rigid kroppsrotasjon, der begge påvirker strukturens respons. Stivhetsmatrisene utformes ofte i symmetrisk form, noe som sikrer konsistens i beregningene og er spesielt viktig ved implementering av iterative løsningsmetoder som Newton-Raphson og displacementskontrollmetoder.
Den iterative inkrementelle analysen er en nødvendig prosedyre for ikke-lineære problemer hvor løsningen oppdateres i små steg. For hver inkrementell fase gjennomføres en korrektorfase som sikrer at både geometri og kreftene oppdateres korrekt i samsvar med den nye tilstanden. Metoder som Newton-Raphson, buelengde-metoden og arbeidkontrollmetoden gir ulike tilnærminger til å sikre konvergens og stabilitet i analysen, spesielt når strukturen nærmer seg kritiske tilstander som buckling eller store rotasjoner.
Forståelsen av disse prinsippene krever også en innsikt i hvordan ikke-lineære strukturer oppfører seg under belastning, spesielt i stadiene før og etter buckling. Her spiller både statikk og kinematikk en rolle for å kunne beskrive hvordan nodale momenter og krefter genereres og påvirker stabilitet. Momentene som oppstår som følge av både stressresultanter og ytre påkjenninger må inkluderes i formuleringen for å oppnå en realistisk beskrivelse av strukturen.
Det er viktig å være klar over at numerisk analyse ikke bare er en teknisk øvelse, men en prosess hvor fysisk forståelse må kombineres med matematiske metoder for å sikre at resultatene er meningsfulle og anvendbare i praksis. Valg av passende elementer, vurdering av stivhetsmatriser, testing av rigid kroppsadferd og korrekt implementering av inkrementelle metoder danner grunnlaget for pålitelighet i beregningene. Det gir ingeniøren et verktøy til å predikere komplekse strukturelle responsmønstre, unngå uønskede kollapser og optimalisere design.
Endelig bør man forstå at selv om avanserte metoder som total og oppdatert Lagrange-formulering, samt inkrementell-iterative metoder, gir muligheter for detaljert og nøyaktig analyse, så fordrer dette et sterkt grunnlag i klassisk mekanikk, numerisk metode og strukturteori. Det å mestre samspillet mellom kinematikk, statikk og materialrespons i numeriske modeller er nøkkelen til å løse utfordringer i moderne rammestrukturanalyse, særlig når ikke-linearitet dominerer.
Hvordan fungerer inkrementell-iterativ ikke-lineær analyse i strukturer, og hva er forskjellene mellom Newton–Raphson og forskyvningskontrollmetoder?
Inkrementell-iterativ analyse er en sentral metode innenfor numerisk løsning av ikke-lineære strukturelle problemer. Denne tilnærmingen kombinerer små belastningsøkninger (inkrementer) med flere iterasjoner per trinn for å oppnå en stabil og nøyaktig løsning på strukturets respons. Notasjonen for iterasjonene er slik at venstre indeks betegner inkrementnummeret, mens høyre indeks indikerer iterasjonssteget innenfor det inkrementet. For tydelighet brukes symbolet Δ for å representere inkrementelle endringer i variabler.
Den generelle formen for de ikke-lineære likningene som styrer strukturen i den j-te iterasjonen av det i-te inkrementet kan skrives som
hvor er vektorene med forflytningsinkrementer, er tangentstivhetsmatrisen, og og er henholdsvis påførte krefter og indre krefter fra forrige iterasjon. Tangentstivheten består av både elastisk og geometrisk stivhet, samt momentstivhet for romrammer, noe som gjør modellen svært fleksibel for forskjellige strukturelle systemer.
En vesentlig forbedring i denne metoden, sammenlignet med klassiske Newton–Raphson, er at lasten kan variere innenfor iterasjonene i samme inkrement, i stedet for å holdes konstant. Dette er avgjørende for å kunne løse problemer som oppstår nær grensepunkter (limit points), hvor Newton–Raphson-metoden ofte divergerer.
Lastevektoren kan dekomponeres i et referanselastvektor multiplisert med en lastfaktor , samt en lastøkning. Dette gjør at kan behandles som en ekstra ukjent i systemet, sammen med forflytningsinkrementene, og danner grunnlaget for den såkalte N+1 teorien. Problemet deles ofte opp i to delsystemer, der det ene løser for responsen på referanselasten og det andre for balanse mellom påførte og indre krefter, noe som gir både numerisk effektivitet og fleksibilitet i løsningen.
Den klassiske Newton–Raphson-metoden øker lasten med et konstant trinn i første iterasjon av hvert inkrement, men holder den konstant i påfølgende iterasjoner. Dette kan oppsummeres som:

Deutsch
Francais
Nederlands
Svenska
Norsk
Dansk
Suomi
Espanol
Italiano
Portugues
Magyar
Polski
Cestina
Русский