Beregningen av deflasjonen til Bernoulli-bjelker under ulike lastekombinasjoner er et sentralt tema innenfor strukturell analyse. For å forstå hvordan bjelkene reagerer på ytre belastninger, er det viktig å begynne med de grunnleggende betingelsene som styrer systemet, og deretter bruke spesifikke funksjoner for å beskrive lastens distribusjon.

En vanlig situasjon for bjelker som er støttet på begge ender, er å finne forholdene mellom de indre reaksjonene og de eksterne kreftene som virker på bjelken. Grunnleggende randbetingelser for slike problemer kan uttrykkes som:

  • Qy(x=0)=R1,Qy(x=L)=+R2Q_y(x=0) = -R_1, Q_y(x=L) = +R_2

  • Mz(x=0)=0,Mz(x=L)=0M_z(x=0) = 0, M_z(x=L) = 0

  • uy(x=0)=0,uy(x=L)=0u_y(x=0) = 0, u_y(x=L) = 0

Disse betingelsene reflekterer det faktum at bjelken er festet i begge ender, og at både vertikal forskyvning og moment er null ved støttepunktene. Etter å ha etablert disse, kan vi gå videre til å finne uttrykkene for deflasjonen og momentene langs bjelken.

For laster som har en trekantet form, kan vi benytte oss av Föppl-funksjoner for å uttrykke hvordan lasten begynner og stopper i spesifikke punkter på bjelken. En triangel-lastefordeling q(x)q(x) kan være representert ved en funksjon som har form av en lineær variasjon mellom to punkter. Dette kan deretter settes inn i de relevante partielle differensialligningene for å integreres flere ganger, og på den måten finne uttrykket for deflasjonen uy(x)u_y(x).

En viktig prosess i beregningen av bjelkens deflasjon er den flertrinnede integrasjonen av differensialligningene, som krever bruk av delvis integrasjon. Dette gir en rekke formler som til slutt kan brukes til å bestemme hvordan bjelken bøyer seg under lasten. Ved å ta hensyn til randbetingelsene som er spesifikert for moment og forskyvning, kan vi finne de nødvendige integrasjonskonstantene og løse for uy(x)u_y(x), som gir deflasjonen ved alle punktene på bjelken.

Når vi har beregnet de forskjellige integrasjonene, får vi et uttrykk som gir oss den vertikale forskyvningen uy(x)u_y(x) langs bjelkens lengde. I noen spesifikke tilfeller kan vi bruke de klassiske formlene som er oppgitt i litteraturen, for eksempel når lasten er jevnt fordelt over hele bjelkens lengde.

I tillegg er det viktig å merke seg at når vi arbeider med bjelker i det elastoplastiske området, må vi også vurdere materialets oppførsel. Bjelken kan gjennomgå både elastiske og plastiske deformasjoner, noe som kan endre stressfordelingen langs bjelken. I det elastoplastiske området vil materialet oppføre seg annerledes enn i det rent elastiske området, spesielt når bjelken er utsatt for høye bøyningsmomenter. Dette krever at vi bruker de relevante materialmodeller for plastisitet og justerer analysen deretter.

For å analysere materialets respons i dette området, kan vi bruke enkle modeller som antar ideell plastisitet, der materialet har et klart definert flytespenning, og der deformasjonene skjer når spenningene overstiger denne grensen. Videre kan vi beskrive spenningsfordelingen over tverrsnittet av bjelken som en funksjon av bøyningsmomentene som virker på bjelken.

For å fullstendig forstå fenomenene som skjer i elastoplastiske bjelker, må man i tillegg til de matematiske formlene også ha en god forståelse av materialets grunnleggende oppførsel under belastning, og hvordan dette påvirker deflasjon og spenningsfordeling. Dette krever en grundig forståelse av plastisitetsmodeller og evnen til å anvende disse i praktiske situasjoner.

Derfor, i tillegg til de teoretiske beregningene som er beskrevet, er det viktig å ha praktiske referanser for elastoplastisk bøyning av bjelker. Klassiske referanseverk som de av Nádai, Prager, og Hodge gir dyptgående beskrivelser av elastoplastisk analyse og er uunnværlige for den som ønsker å mestre disse komplekse problemene.

Med denne forståelsen kan man ikke bare analysere enkle bøyningsproblemer, men også håndtere mer komplekse situasjoner hvor materialet går inn i det plastiske området, og hvor store deformasjonsverdier må tas i betraktning.

Hvordan Gap i Stenger Påvirker Belastning og Reaksjonskrefter i Strukturer

I denne analysen vurderer vi en stangstruktur som er utsatt for forskjellige laster og hvilke effekter et gap i strukturen har på dens respons. Strukturens materialegenskaper er definert ved Youngs modul E, og tverrsnittsarealet A er konstant for visse deler av analysen. Vi ser nærmere på tilfellene med punktbelastninger på forskjellige steder, samt hvordan gapet påvirker både forskyvninger, reaksjonskrefter og stressfordeling i strukturen.

En viktig egenskap ved strukturen er at gapet mellom den ene enden av stangen og den stive veggen introduserer en ikke-lineær respons etter hvert som belastningen øker. I tilfelle det ikke er kontakt mellom stangen og veggen, vil forskyvningen i midten av stangen (u2) avhenge direkte av den påførte belastningen. Når gapet lukkes, skjer det en endring i systemets stivhet, og begge stangsegmentene begynner å bidra til den globale stivhetsmatrisen.

I det første tilfellet, når punktbelastningen påføres i midten av stangen, kan forskyvningen u2 beregnes som F / (EA), hvor F er den påførte belastningen, E er Youngs modulus, A er tverrsnittsarealet, og L er lengden på stangen. Når F øker til et kritisk nivå, som er lik EAδ / L, der δ er avstanden som gapet kan lukkes, vil kontakten mellom stangen og veggen begynne. Denne kontakten endrer det globale systemet av ligninger, og forskyvningene u2 og u3 vil være lik δ, som er gapets størrelse.

For reaksjonskreftene innebærer dette at både R1 og R3 er rettet mot den negative X-retningen etter at gapet er lukket. Reaksjonskreftene kan beskrives ved formler som:

R1=F+EAδ2LR1 = \frac{F + EA\delta}{2L}
R3=FEAδ2LR3 = \frac{ - F - EA\delta}{2L}

Stressfordelingen i stangen vil også endre seg med økningen i belastning. Før gapet lukkes, vil stresset i stangen følge den klassiske formelen for elastisk materialrespons, men etter at gapet er lukket, vil stangen absorbere all ekstra belastning ved støttepunktene, og stresset i stangen vil ikke lenger endres betydelig med økende belastning.

I det andre tilfellet, hvor punktbelastningen påføres på høyre ende, skjer en lignende prosess, men forskyvningene i midten av stangen blir uavhengige av plasseringen av punktbelastningen. Reaksjonskreftene kan igjen beregnes ved hjelp av samme fremgangsmåte som før, og stresset i stangen vil utvikle seg i henhold til materialets elastiske egenskaper.

En annen viktig vurdering er hvordan tverrsnittsarealet påvirker systemets respons. Hvis tverrsnittsarealet endres lineært langs stangen, vil det føre til en variabel stivhet, og reaksjonene vil måtte beregnes for hvert delsegment. Dette kan løses ved å bruke et discretisert finitt elementnett, der hvert element er modellert med et gjennomsnitt av tverrsnittsarealet på venstre og høyre kant.

Denne analysen belyser hvordan gapet i en struktur, samt variasjoner i belastningens plassering og tverrsnittsareal, kan påvirke både de mekaniske responsene i form av forskyvninger og reaksjonskrefter, samt hvordan man kan bruke finitte elementer for å nøyaktig simulere disse effektene.

Det er viktig å merke seg at selv om gapet ikke påvirker strekk- og trykktilstandene i stangen før det lukkes, vil det ha stor betydning når det først skjer. Etter at gapet er lukket, vil all ekstra belastning overføres til støttepunktene, og stangens stressfordeling vil forbli uendret. Dette gjør det mulig å forutsi reaksjoner i strukturer med gap på en mer presis måte, og forståelsen av hvordan gapet påvirker belastningen kan være viktig i mange ingeniørdisipliner, fra byggkonstruksjoner til maskindesign.

Hvordan oppnå en svak formulering av beam-elementer i finitt elementmetoden

I den klassiske tilnærmingen for løsningen av bjelkebøyning, representerer den eksakte løsningen av problemet en løsning som er oppfylt ved hver enkelt posisjon xx på bjelken, og kalles den sterke formuleringen. Ved å erstatte den eksakte løsningen med en tilnærmet løsning uy(x)u_y(x), oppstår et residual r(x)r(x), som kan skrives som:

r(x)=d4uy(x)dx4qy(x)=0.r(x) = \frac{d^4 u_y(x)}{dx^4} - qy(x) = 0.

Som et resultat av introduksjonen av den tilnærmede løsningen, er det generelt ikke lenger mulig å tilfredsstille den differensialligningen på hver enkelt plassering xx på bjelken. Innenfor rammen av den veide residualmetoden, blir det derimot krevd at den differensialligningen oppfylles over et visst lengdeintervall, og den følgende integrale utsagnet blir oppnådd:

0LW(x)(d4uy(x)dx4qy(x))dx=0,\int_0^L W(x) \left( \frac{d^4 u_y(x)}{dx^4} - qy(x) \right) dx = 0,

hvor funksjonen W(x)W(x) kalles vektingsfunksjonen som fordeler feilen eller residualen i det betraktede området. Ved å integrere den første delen av parentesene i ligning (5.210) ved deler, får vi:

0L[d3uydx3W(x)]dx=0L[EIzd2Wdx2]dx.\int_0^L \left[ \frac{d^3 u_y}{dx^3} W(x) \right] dx = \int_0^L \left[ EIz \frac{d^2 W}{dx^2} \right] dx.

Videre fører integrering ved deler på høyre side til:

0L[EIzd2Wdx2]dx=0.\int_0^L \left[ EIz \frac{d^2 W}{dx^2} \right] dx = 0.

Kombinasjonen av de to resultatene gir, etter betraktning av den svake formuleringen, en symmetrisk formulering med hensyn til derivatene av både den tilnærmede løsningen og vektingsfunksjonen. Denne symmetrien med hensyn til derivatene av den tilnærmede løsningen og vektingsfunksjonen garanterer at man får en symmetrisk stivhetsmatrise for bjelkeelementet, akkurat som i tilfelle av stangelementer.

For å kunne fortsette utledningen av den viktigste finitte elementligningen, må forskyvningen uy(x)u_y(x) og vektingsfunksjonen W(x)W(x) uttrykkes innenfor et element, som er nødvendig i form av nodaltilnærmingen. Denne tilnærmingen kan generelt uttrykkes som:

uy(x)=NT(x)up,u_y(x) = N^T(x) u_p,

hvor N(x)N(x) er kolonnematrise for formfunksjonene, og upu_p er kolonnematrise for de nodale ukjente variablene.

Videre kan vi uttrykke spenningsfordelingen innenfor et element på en lignende måte som:

ε(x)=BTup,\varepsilon(x) = B^T u_p,

hvor BB er en matrise som relaterer de deriverte av formfunksjonene til nodale ukjente. Dette gjør det mulig å derivere stivhetsmatrisen KeK_e i henhold til formelen:

Ke=VBTCBdV.K_e = \int_V B^T C B dV.

For bjelkeelementer, kan vi bruke en definisjon for BB-matrisen som i:

B=d2Ndx2.B = \frac{d^2 N}{dx^2}.

Når man deretter setter opp nodaltilnærmingene for uy(x)u_y(x) og W(x)W(x) i den svake formuleringen, får vi den endelige uttrykk for stivhetsmatrisen:

Ke=0L[EIzd2Ndx2N(x)]dx.K_e = \int_0^L \left[ EIz \frac{d^2 N}{dx^2} N(x) \right] dx.

Dette gjør det mulig å formulere den elementære stivhetsmatrisen for bjelkeelementer i et finitt elementnettverk. I Bernoulli-bjelkeelementet består elementet av to noder, som hver har to frihetsgrader: en forskyvning uyu_y og en rotasjon ϕz=duydx\phi_z = \frac{du_y}{dx}, og kan belastes av en kraft FyF_y og et bøyningsmoment MzM_z. Hver node har dermed to ukjente variabler, og det totale antallet ukjente i elementet blir fire.

Formuleringen av nodaltilnærmingene for de fire frihetsgradene kan skrives som:

uy(x)=N1(x)u1+N2(x)u2,u_y(x) = N_1(x) u_1 + N_2(x) u_2,

hvor N1(x)N_1(x) og N2(x)N_2(x) er formfunksjonene for forskyvningene, og u1u_1 og u2u_2 er de nodale forskyvningene for de to nodene. For rotasjonen brukes de samme formfunksjonene:

ϕz(x)=N1(x)ϕ1+N2(x)ϕ2,\phi_z(x) = N_1(x) \phi_1 + N_2(x) \phi_2,

hvor ϕ1\phi_1 og ϕ2\phi_2 er de nodale rotasjonsvinklene.

I en vanlig praksis er formfunksjonene for bjelkeelementet ofte hentet fra Hermite-interpolasjon, som er en kubisk interpolasjon som garanterer både kontinuitet av forskyvning og moment.

For videre utvikling av stivhetsmatrisen og for å anvende metoden i praktiske beregninger, må de eksakte formfunksjonene for N1(x)N_1(x), N2(x)N_2(x) og deres deriverte dN1dx\frac{dN_1}{dx}, dN2dx\frac{dN_2}{dx} bestemmes, og elementets stivhetsmatrise kan deretter beregnes ved integrasjon.

For videre forståelse av hvordan finitt elementmetoden kan anvendes i praktiske ingeniørproblemer, er det viktig å merke seg at valg av formfunksjoner og nøyaktigheten av tilnærmingen har stor innvirkning på kvaliteten av løsningen. Feil i formfunksjonene kan føre til betydelige unøyaktigheter, spesielt når komplekse geometriske eller lastforhold er involvert. Derfor er det avgjørende å bruke passende interpolasjonsteknikker som tilnærmer de fysiske egenskapene nøyaktig, og alltid vurdere feilen i de numeriske resultatene gjennom residualberegninger.

Hvordan bruke den endelige differensmetoden for analyse av stenger med elastiske egenskaper

Når vi undersøker stenger i det elastiske området, er det viktig å forstå hvordan man kan modellere og analysere slike systemer ved hjelp av metoder som den endelige differensmetoden (FDM). Denne metoden gir en numerisk tilnærming til løsningen av partielle differensialligninger som beskriver elastiske problemer, og kan anvendes på en rekke praktiske problemstillinger som involverer stenger utsatt for både fordelt last og påførte krefter.

For et lastet stangsystem med en konstant material- og geometriparameter, som for eksempel et konstant elastisitetsmodul og tverrsnittsområde, kan problemet beskrives av den partielle differensialligningen:

d2udX2=p(X)EA\frac{d^2 u}{dX^2} = - \frac{p(X)}{EA}

hvor u(X)u(X) representerer forskyvningen, p(X)p(X) er den distribuerte lasten, EE er elastisitetsmodulen, og AA er tverrsnittsarealet. I tillegg, for å løse problemet, kreves det tilstrekkelige randbetingelser. For en stang med fast venstre ende og en påført forskyvning eller kraft i høyre ende, kan de relevante randbetingelsene være:

u(X=0)=0ogu(X=L)=uu(X = 0) = 0 \quad \text{og} \quad u(X = L) = u

eller, hvis høyre ende er påført en kraft FF, kan den tilsvarende randbetingelsen være:

dudXX=L=FEA\frac{du}{dX}\Big|_{X=L} = \frac{F}{EA}