i
Förord
Denna rapport utgör en vidare fördjupning inom ett av de kursmoment som behandlar stålbalkars kapacitet vid böjknäckning i kursen ”AF2130 – Dimensionering av
stålkonstruktioner enligt Eurocode 3”.
I samband med ovan nämnda kurs har jag fått ovärderlig hjälp av kursens examinator och tillika kursansvarige, universitetslektor Bert Norlin, som ägnat mycket tid utanför lektionerna till att ge mig fördjupad kunskap inom områdena stålkonstruktion och Eurocode 3.
Jag vill också rikta ett särskilt tack till Maija Engberg, studievägledare för Samhällsbyggnadsprogrammet på KTH, som gjorde denna specialinriktning på kandidatarbetet möjlig.
Stockholm, maj 2018 Jonas Gerland
ii
iii
Sammanfattning
I den mest grundläggande delen av Eurocode 3, SS-EN 1993-1-1, finns olika sätt att
kontrollera bärförmågan för en balkpelare som utsätts för böjknäckning, d.v.s. en kombination av transversell och axiell last. Genom att använda formler där det dimensionerande värdet på andra ordningens moment ingår kan man kontrollera om balkpelaren håller för den aktuella belastningen.
Beräkningar av andra ordningens moment ger differentialekvationer som kan vara besvärliga och tidskrävande att lösa, men beräkningarna kan förenklas om man känner till den
momentfaktor, kallad Cm-faktor, som tar hänsyn till momentfördelningens form i det aktuella lastfallet.
Denna rapport visar hur man kan härleda både exakta och approximativa värden på Cm- faktorer genom att studera två lastfall som skulle kunna vara av intresse för praktiskt bruk då konsolbalkar utsätts för böjknäckning.
Resultaten visar att den approximativa metoden, som innebär avsevärt förenklade beräkningar, ger fullgoda resultat för de två fall som behandlas i rapporten.
Cm-faktorerna presenteras som grafer för att kunna läsas av då normalkraften som angriper i konsolbalkens fria ände varierar. I rapporten ges även bakgrunden till begreppet
momentfaktor genom förklaring och beräkning av förstoringsfaktorn, ofta kallad f.
En sammanfattning av de i rapporten härledda Cm-faktorerna redovisas i Tabell 1.
Tabell 1: Cm- faktorer för två olika lastfall.
Fall Upplagsvillkor
& laster
Momentfaktorer (Cm)
”Exakt1”
x (1 − N N⁄ cr)
”Dimensio- neringsvärde2”
”Approx3”
1 x =tan(kL)kL (1 + 1
1− π2 k2L2
) 1 − 0,446 N
Ncr 1 − 0,428 N Ncr
2 x =βL1 [sin[k L (β−1)]+sin(k L)
k cos(k L) ]
Beror av
värdet på β 1 + [−π2β(β−3)
24 − 1] N
Ncr
1 Motsvarar den exakta lösningen till den styrande differentialekvationen.
2 Linjär approximation av den exakta Cm- faktorn.
3 Cm-faktor helt baserad på en approximativ beräkning av andra ordningens moment via traditionell förstoringsfaktor.
iv
v
Abstract
In the most basic part of Eurocode 3, SS-EN 1993-1-1, there are various ways to verify the load carrying capacity of a beam column exposed to flexural column buckling, i.e. a combination of transverse and axial loads, by using different formulas.
These formulas include the value of the second order moment that is provided by differential equations that can be tedious and difficult to solve. By using a correction factor Cm that takes the shape of the moment distribution for the current load case into account, the calculations can be considerably simplified.
This report shows how to derive both exact and approximate Cm-factors by studying two load cases that might be of practical use when a beam column is subjected to flexural column buckling.
The results show that the approximate method, which involves considerably simplified calculations, gives good results for the two load cases discussed in the report.
The Cm-factors are presented graphically for varying values of the normal force acting at the free end of the beam column. The report also gives the background to the concept of Cm- factors by explaining and deriving the magnification factor, often called f.
A summary of the Cm-factors derived in the report is shown in Table 2.
Table 2: Cm-factors for two different load cases.
Case Boundary conditions and load cases
Moment factors (Cm)
”Exact4”
x (1 − N N⁄ cr)
”Dimesion value5”
”Approx6”
1 x =tan(kL)kL (1 + 1
1−π2 k2L2
) 1 − 0,446 N
Ncr 1 − 0,428 N Ncr
2 x =βL1 [sin[k L (β−1)]+sin(k L)
k cos(k L) ]
Varies with β
1 + [−π2β(β−3)
24 − 1] N
Ncr
4Corresponds to the exact solution of the governing differential equation.
5 Linear approximation of the exact Cm-factor.
6 Cm-factor entirely based on an approximate calculation of second order moment via traditional magnification factor.
vi
vii
Innehållsförteckning
Förord ... i
Sammanfattning ... iii
Abstract ... v
Kapitel 1 - Inledning ... 1
1.1 Bakgrund ... 1
1.2 Syfte ... 2
1.3 Beskrivning av rapporten ... 2
1.4 Avgränsningar ... 2
1.5 Rapportens struktur och innehåll ... 2
Kapitel 2 - Analys enligt Eurocode 3 ... 3
2.1 Olika dimensioneringsmetoder där Cm-faktorer kan användas i Eurocode 3 ... 3
2.1.1 Snittkontroll enligt SS-EN 1993-1-1 6.2 ... 3
2.1.2 Interaktionssamband enligt SS-EN 1993-1-1 6.3.3 ... 5
Kapitel 3 - Andra ordningens teori och momentfaktorer ... 7
3.1 Teoretiska förutsättningar ... 7
3.2 Grundläggande materialsamband ... 7
3.3 Andra ordningens teori ... 8
3.4 Bakgrund till begreppet förstoringsfaktor... 8
3.5 Ekvivalent momentfaktor ... 11
Kapitel 4 - Fallstudier ... 13
4.1 Konsol med sinusformad last och normalkraft ... 13
4.2 Konsolbalk med punktlast i varierande läge samt normalkraft ... 17
Kapitel 5 - Resultat och slutsats ... 23
5.1 Resultat ... 23
5.2 Slutsats ... 23
Källförteckning ... 25
Bilaga - Beräkningsexempel ... 27
1
Kapitel 1 - Inledning
1.1 Bakgrund
Böjknäckning innebär att en balk eller pelare samtidigt belastas av tryckande normalkraft och böjande moment, och då den fungerar som både balk och pelare samtidigt benämns den alltså som ”balkpelare”.
Det böjande momentet orsakas av yttre transversallast och/eller böjande ändmoment (t.ex.
orsakade av excentriskt angripande normalkraft).
Böjknäckning är ett plant problem vilket innebär att balkpelaren endast böjer ut i samma plan som transversallasten och/eller ändmomenten verkar.
Om normalkrafterna är stora kan de få stor inverkan på balkböjningen, vilket gör att man måste ta hänsyn till dem vid dimensionering.
Då man vill veta exempelvis momentfördelningen för en balkpelare som utsätts för samtidig transversell last och stor normalkraft måste därför en jämviktsekvation ställas upp i
balkpelarens deformerade läge, vilket ger en differentialekvation av andra ordningen som är besvärligare och mer tidskrävande att lösa än om inte hänsyn hade tagits till den verkande normalkraften. Detta kallas för beräkningar enligt andra ordningens teori.
I de formler som anges i den första delen av Eurocode 3, SS-EN 1993-1-1, behöver man känna till värdet på det dimensionerande moment som verkar i något snitt av balkpelaren för att kontrollera bärförmågan vid den aktuella belastningen.
Vid böjknäckning beräknas detta dimensionerande moment utifrån andra ordningens teori.
För att slippa att använda andra ordningens teori kan man använda momentfaktorer och med hjälp av ex. tabellfall utifrån första ordningens teori få fram både exakta och approximativa värden på andra ordningens moment för att kontrollera balkpelarens bärförmåga. Dessa momentfaktorer kallas för Cm- faktorer och tar hänsyn till momentfördelningens form för det aktuella belastningsfallet.
I denna rapport görs den typen av beräkningar som i Sarvell (2012), men för två elementarfall som inte behandlas i rapporten eller i Appendix A i SS-EN 1993-1-1. För en djupare
beskrivning av de begränsningar som formelpaketet för bärförmågan vid böjknäckning i Eurocode 3 har hänvisas till Sarvell (2012).
2
1.2 Syfte
Syftet med denna rapport är härleda Cm-faktorer för två elementarfall då balkpelare utsätts för böjknäckning.
1.3 Beskrivning av rapporten
Denna rapport beskriver teorin bakom momentfaktorer och visar genom två beräkningsexempel hur Cm-faktorer kan tas fram för balkpelare.
1.4 Avgränsningar
Cm-faktorerna begränsas till att användas endast i det plana böjknäckningsfallet.
1.5 Rapportens struktur och innehåll
I kapitel 2 beskrivs hur bärförmågan vid böjknäckning kontrolleras i SS-EN 1993-1-1.
I kapitel 3 ges bakgrunden till andra ordningens teori och härledning av förstoringsfaktorn, samt hur denna är relaterad till Cm-faktorn.
I kapitel 4 härleds exakta och approximativa Cm-faktorer för två fall då balkpelare utsätts för böjknäckning.
I kapitel 5 redovisas resultat och rapportens slutsats.
I bilagan ges ett beräkningsexempel där det visas hur man kan tillämpa Cm-faktorer i ett praktiskt dimensioneringsfall med andra ordningens beräkning och snittkontroll enligt SS-EN 1993-1-1 som redovisas i kapitel 2.
3
Kapitel 2 - Analys enligt Eurocode 3
2.1 Olika dimensioneringsmetoder där C
m-faktorer kan användas i Eurocode 3
Då man vill kontrollera bärförmågan hos en balkpelare som utsätts för böjknäckning utifrån kriterierna i Eurocode 3, SS-EN 1993-1-1, kan man använda snittkontroll eller
interaktionskontroll.
Snittkontroll måste genomföras om det inte är uppenbart att snittkontroll är gynnsammare än stabilitetskontroll med interaktionssamband.
Notera att de flesta ekvationerna i detta kapitel har getts samma nummer som i SS-EN 1993- 1-1. Redogörelsen av kapitlets innehåll är hämtad ur Norlin (2003).
2.1.1 Snittkontroll enligt SS-EN 1993-1-1 avsnitt 6.2
Balktvärsnitten förutsätts tillhöra tvärsnittsklass (TK) 1, 2 eller 3. Balkpelarens
initialkrokighet, egenspänningar samt inverkan av plasticering för TK 1 och 2 måste beaktas genom att använda tabell 5.1 i SS-EN 1993-1-1.
Vid beräkningar enligt andra ordningens teori måste längs en given struktur alla balktvärsnitt uppfylla snittkontrollen. Andra ordningens moment kan enkelt beräknas exakt eller
approximativt om man känner till Cm-faktorn för det aktuella böjknäckningsfallet.
Det exakta eller approximativa värdet på andra ordningens moment är det dimensionerande värde för momentet som ska användas i brottkriterierna (6.1), (6.2) och (6.31) då balkpelaren utsätts för böjknäckning.
Då man utför snittkontroll utifrån en elastisk analys med ett balktvärsnitt tillhörande tvärsnittsklass 3 måste den mest utsatta delen på balkpelaren uppfylla brottkriteriet (6.1)
( σx,Ed
fy⁄γM0)
2
+ ( σz,Ed
fy⁄γM0)
2
− ( σx,Ed
fy⁄γM0) ( σz,Ed
fy⁄γM0) + 3 ( τEd
fy⁄γM0)
2
≤ 1 (6.1) där σx,Ed är det dimensionerande värdet av den lokala longitudinella
spänningen vid bedömningspunkten.
σz,Ed är det dimensionerande värdet av den lokala tvärspänningen vid bedömningspunkten.
τEd är dimensioneringsvärdet av den lokala skjuvspänningen vid bedömningspunkten.
Genom att summera kvoten mellan dimensioneringsvärde och kapacitet för varje spänningsresultant fås
NEd
NRd+My,Ed
My,Rd+Mz,Ed
Mz,Rd ≤ 1 (6.2)
4
där NEd, My,Ed och Mz,Ed är dimensioneringsvärdena för normalkraft och böjande moment kring y- respektive z-axeln, och där NRd, My,Rd och Mz,Rd är normalkrafts- respektive
momentkapaciteten för motsvarande. (6.1) och (6.2) – speciellt (6.2) – får användas för att på säkra sidan kontrollera tvärsnitt i klass 1 och 2.
För tvärsnitt tillhörande klass 1 eller 2 måste följande villkor vara uppfyllt
MEd ≤ MN,Rd (6.31)
Vilket betyder att det dimensionerande momentet måste vara mindre än eller lika stort som tvärsnittets momentkapacitet då detta samtidigt påverkas av en normalkraft.
För dubbelsymmetriska I-och H-tvärsnitt tillhörande klass 1 eller 2 och med böjning kring y- axeln kan momentkapaciteten kontrolleras genom
MN,y,Rd = Mpl,y,Rd 1−n
1−0,5a men MN,y,Rd< Mpl,y,Rd (6.36) där MN,y,Rd är momentkapaciteten vid y-axelböjning.
Mpl,y,Rd är tvärsnittets plastiska momentkapacitet vid y-axelböjning.
n = NEd
Npl,Rd är kvoten mellan yttre normalkraft och plastisk normalkraftskapacitet.
a =A−2 b tf
A men a ≤ 0,5 kvoten mellan livarea och total tvärsnittsarea.
För dubbelsymmetriska I-och H-tvärsnitt tillhörande klass 1 eller 2 och med böjning kring z- axeln kan momentkapaciteten kontrolleras genom
MN,z,Rd = Mpl,z,Rd om n ≤ a (6.37)
MN,z,Rd = Mpl,z,Rd[1 − (n−a
1−a)2] om n > a (6.38) där MN,z,Rd är momentkapaciteten vid z-axelböjning.
Mpl,z,Rd är tvärsnittets plastiska momentkapacitet vid z-axelböjning.
5
För rektangulära ihåliga tvärsnitt med konstant tjocklek samt för sammansvetsade tvärsnitt med lika flänsar och liv kan momentkapaciteten för y- samt z-axelböjning kontrolleras genom
MN,y,Rd = Mpl,y,Rd 1−n
1−0,5aw men MN,y,Rd< Mpl,y,Rd (6.39) MN,z,Rd = Mpl,z,Rd 1−n
1−0,5af men MN,z,Rd < Mpl,z,Rd (6.40) där aw =(A−2 b t)
A men aw ≤ 0,5 för ihåliga tvärsnitt.
aw = (A−2 b tf)
A men aw ≤ 0,5 för sammansvetsade lådtvärsnitt.
af =(A−2 h t)
A men af≤ 0,5 för ihåliga tvärsnitt.
af =(A−2 b tw)
A men af≤ 0,5 för sammansvetsade lådtvärsnitt.
2.1.2 Interaktionssamband enligt SS-EN 1993-1-1 6.3.3 Begränsningar
För att tillämpa de interaktionssamband för konstruktionselement som är föremål för böjknäckning i SS-EN 1993-1-1 avsnitt 6.6.3 måste hänsyn tas till följande begränsningar:
• Endast dubbelsymmetriska tvärsnitt kan hanteras,
• Tvärsnittsformen får inte ändras vid belastning av t.ex. koncentrerade laster, men tvärsnittsklass 4 är okej.
• Interaktionsformlerna är framtagna och kalibrerade för en vid båda ändar fritt upplagd balk.
Interaktionskontroll
Interaktionssambanden i SS-EN 1993-1-1 avsnitt 6.6.3 används för konstruktionselement som utsätts för böjknäckning med syftet att andra ordningens teori ej behöver tillämpas.
Interaktionssambanden kan användas på två sätt, antingen genom att beräkna de ingående parametrarna enligt Appendix A eller Appendix B. I den nationella bilagan till SS-EN 1993- 1-1 anges att Appendix A bör användas.
Enligt SS-EN 1993-1-1 avsnitt 6.3.3 måste följande villkor vara uppfyllda för en balkpelare som utsätts för böjknäckning:
NEd χyNRk
γM1
+ kyyMy,Ed+∆My,Ed
χLTMy,Rk γM1
+ kyzMz,EdMz,Rk+∆Mz,Ed
γM1
≤ 1 (6.61)
NEd χzNRk
γM1
+ kzyMy,Ed+∆My,Ed
χLTMy,Rk γM1
+ kzz Mz,EdMz,Rk+∆Mz,Ed
γM1
≤ 1 (6.62)
6
där NEd, My,Ed, Mz,Ed är maximala dimensioneringsvärdet för yttre normalkraft och böjmoment kring y-respektive z-axeln längs med hela balkpelaren.
∆My,Ed, ∆Mz,Ed är tillskottsmoment kring y-respektive z-axeln som uppkommer p.g.a. att neutrala lagrets läge kan ändras något vid areareducering för profiler i tvärsnittsklass 4.
NRk, My,Rk, Mz,Rk är normalkraftskapaciteten samt momentkapaciteten kring y- respektive z- axeln. Dessa värden är beroende av tvärsnittsklass.
χy, χz är reduktionsfaktorer för knäckning av centriskt tryckt stång.
χLT är en reduktionsfaktor för inverkan av vippning.
kyy, kyz, kzy, kzz är interaktionsfaktorer enligt tabellerna A.1 och A.2 i Appendix A i SS-EN 1993-1-1.
Då böjknäckning är ett plant problem kan endast böjning ske kring antingen y- eller z-axeln vilket förenklar (6.61) till
NEd χyNRk
γM1
+ kyyMy,Ed+∆My,Ed
χLTMy,Rk γM1
≤ 1 (2.1)
och (6.62) till
NEd χzNRk
γM1
+ kzz Mz,EdMz,Rk+∆Mz,Ed
γM1
≤ 1 (2.2)
där kyy= Cmy
(1−χyNEd Ncr,y)Cyy
och kzz = Cmz
(1−χzNEd Ncr,z)Czz
för tvärsnittsklass 1 och 2
samt kyy= Cmy
(1−χyNEd
Ncr,y) och kzz = Cmz
(1−χzNEd
Ncr,z) för tvärsnittsklass 3 och 4
Cyy och Czz i formlerna för kyy och kzz i fallen med tvärsnittsklass 1 eller 2 ovan är en korrektionsfaktor som beaktar inverkan av den plastiska responsen.
Ncr,y är balkpelarens kritiska knäckningslast med avseende på knäckning i den styva riktningen.
Cmy är den Cm-faktor som beaktar momentfördelningens form enligt tabellerna A.1 och A.2 i SS-EN 1993-1-1 Appendix A och som härleds för två andra lastfall i denna rapport.
7
Kapitel 3 - Andra ordningens teori och momentfaktorer
3.1 Teoretiska förutsättningar
Förutsättningar för matematisk behandling av elastisk instabilitet som de beskrivs i Höglund (2006):
• Linjärt elastiskt material, d.v.s. töjning och spänning förutsätts vara linjärt relaterade oavsett vilken spänningsnivå som gäller.
• Små töjningar och därmed små böjdeformationer.
• Vid bestämning av kritiska laster anses konstruktionselementen som imperfektionsfria, d.v.s.
ingen initialkrokighet och inga egenspänningar anses förekomma.
3.2 Grundläggande materialsamband
Den grundläggande teorin i detta kapitel är hämtad ur Heyden m.fl. (2011).
Vid ren böjbelastning utgörs materialsambandet av elastiska linjens differentialekvation som beskriver relationen mellan deformation, böjmomentet och böjstyvheten:
d2w
dx2 = −M
EI (3.1) som kan skrivas som
M = −EId2w
dx2 (3.2) Om elastiska linjens differentialekvation löses för en ledat infäst pelare som belastas med en tryckande normalkraft N i båda ändar fås den kritiska knäckningslasten Ncr =π2EI
L2 . Den allmänna lösningen kan skrivas som
Ncr =π2EI
Lcr2 (3.3) där Lcr är knäckningslängden, en fiktiv längd som används för att knäckningslasten skall kunna uttryckas med samma ekvation oberoende av belastningens form, stångens geometri samt upplagsförhållandena.
För pelare med olika upplagsvillkor är de vanligaste de fyra fall som benämns som Eulers knäckningsfall, då Lcr = 2L för en i ena änden fast inspänd pelare, Lcr= 1L för en ledat infäst pelare, Lcr = 0,7L för en i ena änden fast inspänd och i den andra änden ledat infäst pelare och Lcr= 0,5L för en i båda ändar fast inspänd pelare.
8
3.3 Andra ordningens teori
För balkpelare, alltså en balk som utsätts för både transversell och axiell last, måste jämvikten ställas upp i det deformerade läget då normalkrafterna är stora. I momentjämvikten ingår då en okänd utböjning w(x) som insatt i elastiska linjens differentialekvation ger en
differentialekvation av andra ordningen.
Beräkningar enligt andra ordningens teori visar att en stor normalkraft ger större utböjningar och moment i en balkpelare än vad en beräkning enligt första ordningen gör. Ju större normalkraft som balkpelaren utsätts för, desto mindre transversallast klarar den av.
Då normalkraften närmar sig värdet av den kritiska knäckningslasten kommer balkpelaren att kollapsa utan att klara någon transversell belastning alls.
3.4 Bakgrund till begreppet förstoringsfaktor
Om andra ordningens moment skall kunna beräknas måste det finnas en initialkrokighet, som antingen finns där från början eller orsakas av en transversallast eller ändmoment.
Nedan följer en beräkning av utböjningen för en ledat infäst pelare som utsätts för en normalkraft och som är krökt redan då den är obelastad. Pelaren antas ha en initialutböjning med formen av en halv sinusvåg som visas i Figur 3.1.
a(x) = a0sin (πx
L) (3.4)
Figur 3.1: Pelare med sinusformad initialutböjning belastad med normalkraft i ena änden.
9
En momentjämvikt enligt Figur 3.2 ställs upp i pelarens deformerade läge då normalkraften verkar.
M(x) − N(a(x) + w) = 0, där w är tillskottsutböjningen som ges av normalkraften N.
Figur 3.2: Momentjämvikt i snitt för normalkraftsbelastad pelare med initialutböjning.
Detta ger momentekvationen
M(x) = N(a(x) + w) (3.5) som insatt (3.1) ger
d2w
dx2 = −N
EIa0sin (πx
L) −Nw
EI (3.6) Genom att definiera k = √N
EI blir k2 = N
EI och (3.6) kan skrivas om som
d2w
dx2 + k2w = −a0 sin (πx
L) (3.7) Den allmänna lösningen till (3.7) blir
w(x) = A cos(kx) + B sin(kx) +π2k2
L2−k2a0sin (πx
L) (3.8) Randvillkoren w(0) = 0 och w(L) = 0 för (3.8) ger A = B = 0
Om k2 skrivs om genom att utnyttja (3.3) med Lcr= 1L för en ledat infäst pelare blir k2 = Nπ2
NcrL2 och tillskottsutböjningen orsakad av normalkraften N kan skrivas som:
w(x) = N N⁄ cr
1−N N⁄ cra0sin (πx
L) (3.9)
10
I mitten av pelaren blir denna tillskottsutböjning, när x =L
2 sätts in i (3.9) w (L
2) = N N⁄ cr
1−N N⁄ cra0 (3.10) Genom att i (3.4) sätta x =L
2 fås pelarens initialutböjning i mitten a (L
2) = a0 (3.11) Den maximala utböjningen, som alltså inträffar i pelarens mitt och som fås genom addition av (3.10) och (3.11), blir enligt andra ordningens teori
w2 = w (L
2) + a (L
2) = N N⁄ cr
1−N N⁄ cra0+ a0 = a0
1−N N⁄ cr (3.12) Initialutböjningen vid pelarens mitt, L
2 , som utan normalkraftspåverkan var a0 har nu multiplicerats med faktorn
f = 1
1−N N⁄ cr (3.13) Faktorn (3.13) brukar kallas för en förstoringsfaktor då uttryckets nämnare aldrig kan bli större än 1 och förstorar alltså täljarens värde. Om normalkraften N närmar sig värdet av den kritiska knäckningslasten Ncr kommer uttrycket att bli oändligt stort.
Om man ställer upp momentjämvikten och bortser från den tillskottsutböjning w som ges av normalkraften N hade momentet blivit M = Na0
Om (3.10) och (3.11) sätts in i (3.5) fås värdet på andra ordningens moment M2 = N a0
1−N N⁄ cr (3.14) Även här har momentet M = Na0 multiplicerats med faktorn f = 1
1−N N⁄ cr.
För att erhålla en mycket bra approximation av andra ordningens maximala utböjning kan alltså första ordningens utböjning multipliceras med förstoringsfaktorn f = 1
1−N N⁄ cr. Detta förutsätter dock att första ordningens utböjning har en form som liknar formen hos aktuell knäckningsmod. I varianten med enbart sinusformig initialkrokighet har både knäckningsmoden och initialkrokigheten exakt samma form och då ger förstoringsfaktorn exakt lösning. Om formerna avviker lite blir approximationen mycket bra, om formerna är helt olika blir approximationen också helt fel.
11
Då första ordningens moment är formad som eller liknar formen av den kurva som knäckningsmoden motsvarar kan även första ordningens moment multipliceras med
förstoringsfaktorn f för att få en approximation av andra ordningens moment. I de flesta fall är approximationen av andra ordningens utböjning betydligt bättre än approximationen av andra ordningens moment.
3.5 Ekvivalent momentfaktor
Man kan alltså uppnå ett bra värde på andra ordningens utböjning och moment, dock med bättre precision på utböjningen, genom att multiplicera första ordningens utböjning och moment med förstoringsfaktorn f under antagande att de båda i formen liknar den sinuskurva som knäckningsmoden motsvarar.
w2 = w1 1
1−N N⁄ cr (3.15) M2 = M1 1
1−N N⁄ cr (3.16) Om så inte är fallet måste ettan i förstoringsfaktorns täljare ersättas med en korrektionsfaktor Cm som kan ge det exakta eller approximativa värdet på andra ordningens moment, M2. Ekvation (3.16) kan då skrivas som
M2 = M1 Cm
1−N N⁄ cr (3.17) och om man ur denna ekvation löser ut korrektionsfaktorn Cm, som härefter benämns som momentfaktorn eller Cm-faktorn, samt skriver första och andra ordningens moment som funktioner av x-koordinaten längs balkpelaren kan momentfaktorn anges som en funktion av x och kvoten N N⁄ cr:
Cm(x, N
Ncr) =M2(x)
M1(x)(1 − N
Ncr) (3.18) Denna momentfaktor kommer att ha olika värden för varje värde på x-koordinaten längs balkpelaren.
För en balkpelare som belastas med transversallaster eller transversallaster i kombination med koncentrerade moment kan ett approximativt värde på andra ordningens moment beräknas enligt
M2,approx(x) = M1(x) + N[w1(x) − wm(x)]f (3.19)
12
Jämför med ekvation (3.14), där i ekvation (3.19) första ordningens moment M1(x) adderats till följd av de transversallaster eller ändmoment som också kan påverka pelaren.
Faktorn wm är en korrektion som gör att skillnaden w1(x) − wm(x) får en form som så nära som möjligt liknar den på första ordningens moment. Denna faktor används inte i några av de beräkningar som görs i denna rapport, men för en mer ingående förklaring hänvisas till Sarvell (2012).
Om andra ordningens moment, M2 i ekvation (3.18), syftar till det exakta värdet på M2 som fås genom att lösa den styrande differentialekvationen (3.7) (dock med högerledet utbytt mot första ordningens moment M1(x)), så betecknas också Cm som exakt.
Cm,exakt(x, N
Ncr) =M2,exakt(x)
M1(x) (1 − N
Ncr) (3.20) Om M2(x) löses ur (3.18) och sätts lika med (3.19) kan den approximativa Cm- faktorn skrivas som
Cm,approx(x, N
Ncr) = 1 + [Ncr(w1(x)−wm(x))
M1(x) − 1] N
Ncr (3.21) Denna formel återfinns som fall 2 i SS-EN 1993-1-1 Appendix A tabell A.2.
I denna rapport kommer faktorn wm alltid vara lika med noll.
Det bör påpekas att det finns ett alternativt sätt att uttrycka (3.17) i litteratur med koppling till Eurocode 3, där Cm = 1 + δ N
Ncr , och där faktorn δ kallas för ”Dischinger faktorn”, se Lindner m.fl. (2017). Ekvation (3.17) skrivs då som
M2 = M11+δ
N Ncr 1−N
Ncr
(3.22)
13
Kapitel 4 - Fallstudier
Nedan följer två härledningar av Cm-faktorer för fast inspända balkpelare som utsätts för både transversell och axiell last.
För båda fallen gäller att maximalt första ordningens moment M1, maximalt andra ordningens moment M2, och maximal första ordningens utböjning w1 inträffar i den fast inspända änden, dvs där x = L . Utböjningen sker i endast det plan som lasterna verkar.
Cm-faktorerna kommer att redovisas med grafer för olika värden på kvoten N N⁄ cr i intervallet 0 ≤ N N⁄ cr≤ 1. En approximation med ett förstagradspolynom kommer även att göras till kurvan för Cm,exakt.
4.1 Konsol med sinusformad last och normalkraft
Figur 4.1: Konsol med sinusformad last och normalkraft.
Momentjämvikt enligt andra ordningens teori ger enligt Figur 4.1
M2(x) = N w(x) + M1(x) (4.1) där första ordningens moment är
M1(x) =q0Lx
π −q0L2
π2 sin (πx
L) (4.2) För beräkning av (3.21) behövs även värdet på första ordningens böjdeformation
w1(x) = 1
EI[−Lq0(π
3x3+6L3sin(πxL)+6πL2x−3π3L2x)
6π4 ] (4.3)
14
Uttrycket (4.2) insatt i (4.1) ger M2(x) = N w(x) +q0Lx
π −q0L2
π2 sin (πx
L) (4.4) och (4.4) insatt i elastiska linjens differentialekvation, (3.1), ger krökningen enligt
d2w
dx2 = −M2(x)
EI = −1
EI(N w(x) +q0Lx
π −q0L2
π2 sin (πx
L)) (4.5) Detta är en linjär differentialekvation av andra ordningen havande konstanta koefficienter.
Genom att definiera k = √N
EI blir k2 = N
EI, och andra ordningens böjdeformation - lösningen till (4.5) - kan skrivas som
w(x) = A cos(kx) + B sin(kx) + q0L2
Nπ2(1− π2
k2L2)sin (πx
L) −q0Lx
Nπ (4.6) Randvillkoret w(0) = 0 → A = 0
Andra ordningens vinkeländring blir w′(x) = B k cos(kx) + q0L
Nπ(1− π2
k2L2)cos (πx
L) −q0L
Nπ (4.7)
Randvillkoret w′(L) = 0 → B = q0L
N π k cos(kL)(1 + 1
1− π2 k2L2
)
Andra ordningens nedböjning kan nu skrivas som
w(x) = q0L
N π k cos(kL)(1 + 1
1− π2 k2L2
) sin(kx) + q0L2
Nπ2(1− π2
k2L2)sin (πx
L) −q0Lx
Nπ (4.8) Andra ordningens moment M2 erhålles genom att derivera utböjningen två gånger och
multiplicera med −EI.
M2(x) = q0L
π k cos(kL)(1 + 1
1−π2 k2L2
) sin(kx) + q0
k2(1− π2
k2L2)sin (πx
L) (4.9) Genom insättning av x = L i (4.9) fås maximalt andra ordningens moment i inspänningssnittet enligt
15
M2(L) = M2,max= q0L tan(kL)
πk (1 + 1
1− π2 k2L2
) (4.10)
Även maximalt första ordningens moment inträffar i inspänningssnittet, x = L sätts in i ekvation (4.2)
M1(L) = M1,max =q0L2
π (4.11) Genom insättning av (4.10) och (4.11) i (3.20) fås
Cm,exakt( N
Ncr) = tan(kL)
kL (1 + 1
1− π2 k2L2
) (1 − N
Ncr) (4.12) Det maximala värdet på första ordningens utböjning behövs för beräkningen av Cm,approx. Insättning av x = L i (4.3) ger
w1(L) =q0L4
EI (1
3π− 1
π3) (4.13) Insättning av (4.11), (4.13) och (3.3) med knäckningslängden Lcr= 2L, samt wm = 0 i (3.21) ger
Cm,approx( N
Ncr) = 1 + (π2
12−5
4) N
Ncr≈ 1 − 0,4275 N
Ncr (4.14) Graferna i Figur 4.2 visar sambandet mellan Cm-faktorn och kvoten N N⁄ cr.
Den första som benämns ”Exact” visar Cm-värden enligt ekvation (4.12).
Den andra, benämnd ”Exact design curve”, är en kurvanpassning med ett förstagradspolynom till kurvan ”Exact” i intervallet 0 ≤ N N⁄ cr≤ 1.
Den tredje, ”Approx”, visar approximativa värden enligt ekvation (4.14). Knäckningslängden har satts till Lcr= 2,0L enligt Eulers första knäckningsfall.
Graferna visar att approximativ lösning enligt (3.21) eller (4.14) stämmer mycket bra med exakt lösning enligt (3.20) eller (4.12).
16
Figur 4.2: Momentfaktorer för en konsol med sinusformad last.
17
4.2 Konsolbalk med punktlast i varierande läge samt normalkraft
Figur 4.3: Konsolbalk med punktlast i varierande läge samt normalkraft.
Beteckningen wh i följande beräkningar syftar på nedböjningen för den del av balken som ligger till höger om punktlastens verkningspunkt och har ingenting med riktningen på x-axeln att göra. Nedböjningen wv berör den del av balken som ligger mellan punktlastens
verkningspunkt och den fasta inspänningen. Dessa beteckningar används då skarvvillkor förekommer i beräkningen. βL är sträckan mellan den fasta inspänningen och
transversallastens verkningspunkt.
För delen 0 ≤ x ≤ L − βL:
Momentjämvikt enligt Figur 4.3 ger
M2,h(x) = N wh(x) (4.15) Funktionen (4.15) insatt i (3.1) ger krökningen
d2wh
dx2 = −M2,h(x)
EI = −N wh(x)
EI (4.16) Genom definition
k = √N
EI och k2 = N
EI (4.17) skrivs (4.16) som
d2wh
dx2 + k2wh(x) = 0 (4.18)
18
Den allmänna lösningen till (4.18) är
wh(x) = A cos(kx) + B sin(kx) (4.19) Det randvillkor som kan användas är wh(0) = 0 → A = 0
wh(x) = B sin(kx) (4.20) Derivering av (4.20) ger lutningen
wh′(x) = B k cos(kx) (4.21) För delen L − βL ≤ x ≤ L:
Momentjämvikt ger
M2,v = N wv(x) + P(x − L + βL) (4.22) Uttrycket (4.22) insatt i (3.1) ger krökningen
d2wv
dx2 = −M2,v(x)
EI = −1
EI(N w (x) + Px − PL + PβL) (4.23) och genom att använda definitionen (4.17) kan (4.23) skrivas som
d2wv
dx2 + k2wv(x) =k2P
N [−x + L − βL] (4.24) Differentialekvationen (4.24) har lösningen
wv(x) = C cos(kx) + D sin(kx) +P
N[−x + L − βL] (4.25) Derivering ger
wv′(x) = −C k sin(kx) + D k cos(kx) −PN (4.26) Randvillkoret wv′(L) = 0 ger
D = C sin (kL)
cos (kL)+ P
k N cos (kL) (4.27) Ekvationerna (4.25) och (4.26) kan nu skrivas om innan skarvvillkoren används genom att använda (4.27)
19
wv(x) = C cos(kx) + [Csin (kL)
cos (kL)+ P
k N cos (kL)] sin(kx) +P
N[−x + L − βL] (4.28)
wv′(x) = −C k sin(kx) + k [Ccos(kL)sin(kL)+ P
k N cos(kL)] cos(kx) −P
N (4.29) En balk som deformeras elastiskt får inga knyckar vilket innebär att lutning och nedböjning i
den punkt där punktlasten verkar måste vara lika. Skarvvillkoren wh(L − βL) = wv(L − βL) och wh′(L − βL) = wv′(L − βL) ger de hittills okända konstanterna B, C och D:
B = P
k N[sin2k(L−βL)
cosk(L−βL) −sink(L−βL)
cos(kL) sin(kL) + 1
cos(kL)− 1
cos k(L−βL)] (4.30) C = − P
k Nsink(L − βL) (4.31) D = P [ 1
k N cos (kL)−sink(L−βL)
k N
sin(kL)
cos(kL)] (4.32) Genom att använda (4.30), (4.31) och (4.32) skrivs (4.28) och (4.29) som
wv(x) = −P sin k(L−βL)
k N cos(kx) + [ P
k N cos(kL)−P sin k(L−βL) k N
sin(kL)
cos(kL)] sin(kx) +P
N[−x + L − β] (4.33) wv′(x) =Psin k(L−βL)
N sin(kx) + [ P
N cos(kL)−P sin k(L−βL) N
sin(kL)
cos(kL)] cos(kx) −P
N (4.34)
Den del av balken som är betydelsefull då Cm-faktorerna beräknas är den vänstra delen, mellan den fasta inspänningen och punktlasten.
Nedan följer de ekvationer som behövs för att ta fram Cm-faktorerna för den delen, d.v.s. för:
L − βL ≤ x ≤ L
Första ordningens moment:
M1(x) = P(x − L + βL) (4.35) Första ordningens nedböjning:
w1(x) = 1
EI[−P(L3β3−3L3β2 +3L3β−L3−6L2βx+3L2x+3Lβx2−3Lx2+x3)
6 ] (4.36) Andra ordningens moment fås genom att derivera (4.34) vilket ger andra ordningens
krökning. Efter multiplikation med −EI fås det exakta värdet på andra ordningens moment.
M2(x) = −P sin k(L−βL)
k cos(kx) + [ P
k cos(kL)−P sin k(L−βL) k
sin(kL)
cos(kL)] sin (kx) (4.37)
20
Genom insättning av x = L i (4.37) fås maximalt andra ordningens moment i inspänningssnittet enligt
M2(L) = M2,max= Psin[kL(β−1)+sin(kL)]
k cos(kL) (4.38) Även maximalt första ordningens moment inträffar i inspänningssnittet, x = L sätts
in i (4.35)
M1(L) = M1,max = PβL (4.39) Genom insättning av (4.38) och (4.39) i (3.20) fås
Cm,exakt( N
Ncr, β) = 1
βL
sin[k L (β−1)]+sin(k L)
k cos(k L) (1 − N
Ncr) (4.40) Det maximala värdet på första ordningens utböjning behövs för beräkningen av Cm,approx. Insättning av x = L i (4.36) ger
w1(L) =PL3β2(β−3)
6 (4.41) Insättning av (4.39), (4.41) och (3.3) med knäckningslängden Lcr= 2L, samt wm = 0 i (3.21) ger
Cm,approx( N
Ncr, β) = 1 + [−π2β(β−3)
24 − 1] N
Ncr (4.42) Det är viktigt att notera att (4.40) och (4.42) endast är giltiga för den vänstra sidan av balken.
Om β = 0 försvinner den vänstra sidan vilket gör att (4.40) och (4.42) förlorar sin innebörd.
Det som är på den högra sidan är helt enkelt Eulers första knäckningsfall då det inte finns någon initialkrokighet eller någon initiallutning hos stången.
Figur 4.4 och Figur 4.5 visar sambandet mellan Cm-faktorer och kvoten mellan normalkraft och knäckningslast för fallet då β = 1, d.v.s. då punktlasten angriper i konsolens fria ände och fallet då β = 0,5, då punktlasten angriper i konsolens mitt.
När β minskar ökar lutningen på graferna vilket betyder att andra ordningens effekter avtar och att punktlasten P gör mindre skada ju närmre den flyttas mot balkpelarens fasta
inspänning.
21
Figur 4.4: Momentfaktorer för en konsol med punktlast i fria änden, 𝛃 = 𝟏.
Figur 4.5: Momentfaktorer för en konsol med punktlast i balkens mitt, 𝛃 = 𝟎, 𝟓.
22
Figur 4.6 visar för tre olika kvoter på N N⁄ cr hur de approximativa Cm-faktorerna varierar beroende av var på balkpelaren punktlasten angriper mätt från den fasta inspänningen. Man ser återigen tydligt att lutningen på graferna ökar då β minskar, vilket är lätt att se på formel (4.42) som är den Cm-faktor som visas i Figur 4.6. Punktlasten P gör alltså mer skada ju närmre balkpelarens fria ände den angriper.
Figur 4.6: Den i inspänningssnittet approximativa momentfaktorns variation beroende av punktlastens läge.
23
Kapitel 5 - Resultat och slutsats
5.1 Resultat
Resultaten för de i föregående kapitel beräknade Cm-faktorerna presenteras i Tabell 5.1. De tre fallen motsvarar ”Exakt” - beräknad utifrån andra ordningens teori,
”Dimensioneringsvärde” - en kurvanpassning med ett förstagradspolynom till det exakta värdet, samt ”Approx” – beräknad med den approximativa metoden. Notera att det är svårt att skriva ett enkelt uttryck för kurvanpassningen ”Dimensioneringsvärde” för fall 2 som beror av värdet på β utmed balkpelaren.
Tabell 5.1: Cm- faktorer för två olika lastfall.
Fall Upplagsvillkor
& laster
Momentfaktorer (Cm)
”Exakt7”
x (1 − N N⁄ cr)
”Dimensio- neringsvärde8”
”Approx9”
1
x =tan(kL)kL (1 +1−1π2k2L2
) 1 − 0,446 N
Ncr 1 − 0,428 N Ncr
2
x =βL1 [sin[k L (β−1)]+sin(k L)k cos(k L) ]
Beror av
värdet på β 1 + [−π2β(β−3)
24 − 1] N
Ncr
5.2 Slutsats
I denna rapport har jag noga redovisat beräkningar för två fall av exakta och approximativa värden på Cm-faktorer som jag ej har funnit i litteratur som behandlar dylika fall. Resultaten i rapporten visar att det i dessa två fall är mycket liten skillnad på de exakta och approximativa värdena.
7 Motsvarar den exakta lösningen till den styrande differentialekvationen.
8 Linjär approximation av den exakta Cm- faktorn.
9 Cm-faktor helt baserad på en approximativ beräkning av andra ordningens moment via traditionell förstoringsfaktor.
24
25
Källförteckning
Skriftliga källor
EN 1993-1-1. Eurocode 3: Design of steel structures - part 1-1: General rules and rules for buildings. (2005). European Committee for Standartisation.
Heyden, S., Dahlblom, O., Olsson, A., & Sandberg, G. (2011). Introduktion till strukturmekaniken. Lund: Studentlitteratur.
Höglund, T. (2006). Modul 6 Stabilitet för balkar och stänger. Stockholm: Kunliga Tekniska Högskolan.
Lindner, J., Kuhlmann, U., & Jörg, F. (2017). Initial bow imperfections e_0 for the verification of Flexural Buckling According to Eurocode 3 Part 1-1- additional considerations.
Muratagic, D. (2010). Verifiering av interaktionsformlers gränser i Eurokod 3.
Norlin, B. (2003). Dimensionering av stålkonstruktioner enligt Eurocode 3, föreläsningar &
tal.
Sarvell, F. (2012). Verifiering av interaktionsformlers gränser för balkpelare i Eurokod 3 En fallstudie av balkpelare med godtyckliga last- och randvillkor. Stockholm.
Muntliga källor
Norlin, Bert. Universitetslektor på Institutionen för byggvetenskap, bro- och stålbyggnad, återkommande samtal och granskningar under vårterminen 2018.
26
27
Bilaga - Beräkningsexempel
I följande exempel beräknas normalkraftsbärförmågan för tvärsnittet HEA 280 då ett elastiskt brottkriterium enligt (6.2) och ett plastiskt brottkriterium enligt (6.31) i SS-EN-1993-1-1 används. Balkpelaren kan bara böjknäcka i veka riktningen och är fast inspänd med en normalkraft verkande i den fria änden.
Tvärsnittets sträckgräns och balkpelarens längd är olika för det elastiska och det plastiska fallet. Detta är för att samma tvärsnittsprofil skall kunna användas tillsammans med båda brottkriterierna samt för att anpassa så att slankhetstalet λ, som är ett mått på en tryckt pelares benägenhet att knäcka, i båda fallen blir lika med ett.
I denna bilaga används Eurokodens numrering av ekvationer hämtade från denna standard.
Tvärsnittsdata HEA 280
Stålets E-modul E = 210 GPa Partialkoefficient γM1 = 1,0
Profilbredd b = 280 mm Profilhöjd h = 270 mm Livtjocklek tw = 8,0 mm Flänstjocklek tf= 13,0 mm Hålkälsradie R = 24 mm Tvärsnittsarea A = 9726 mm2 Tröghetsmoment Iz = 47,63 ∗ 106 mm4 Elastiskt böjmotstånd Wz,el = 340 ∗ 103 mm3 Plastiskt böjmotstånd Wz,pl = 518 ∗ 103 mm3
Elastiskt brottkriterium enligt SS-EN-1993-1-1 (6.2)
Tvärsnittsklassning
Flytspänning fy = 620 MPa
28
Roten ur kvoten mellan referensflytspänningen och aktuell flytspänning
ε = √235 MPa
fy = 0,616
Slankhet hos flänsplåten som utsätts för tryck
cf= b 2 −
tw 2 − R
𝑡𝑓 = 8,62
Gränser för tryckt fläns enligt Tabell 5.2 SS-EN-1993-1-1
( 9ε 10ε 14ε
) = ( 5,54 6,16 8,62
)
Klass 1 Tryckflänsen Klass 2 tillhör klass 3 Klass 3
Slankhet hos livplåten under antagande att hela livplåten är tryckt
cw =h − 2 tf− 2R
𝑡𝑤 = 24,5
Gränser för tryckt livplåt enligt Tabell 5.2 SS-EN-1993-1-1
( 33ε 38ε 42ε
) = ( 20,3 23,4 25,9
)
Klass 1 Livplåten Klass 2 tillhör klass 3 Klass 3
Tvärsnittet tillhör alltså tvärsnittsklass 3.
Bärförmåga med hänsyn till normalkraft och momentkapacitet
Tvärsnittets normalkraftsbärförmåga NRd = Afy
γM1 = 6,03 ∗ 103 kN (6.10) Tvärsnittets elastiska momentbärförmåga
Mz,Rd =Wz,el fy
γM1 = 211 kNm (6.14)
29
Slankheten
Då slankheten λ har bestämts till värdet 1 för både den elastiska och plastiska analysen uttrycks längden L1 via parametrarna: slankheten λ, elasticitetsmodulen E, tröghetsmomentet Iz och normalkraftsbärförmågan NRd.
L1 =π
2λ√NEIz
Rk = 2,02 m Kritisk knäckninglast
Den kritiska knäckningslasten beräknas Ncr,z =π2EIz
Lcr2 = 6,03 ∗ 103 kN där Lcr= 2L1
Bärförmåga för knäckning enligt SS-EN 1993-1-1 6.3.1.1 & 6.3.1.2
Brottvillkoret
NEd
Nb,Rd≤ 1,0 (6.46) där normalkraftskapaciteten för knäckning är
Nb,Rd =χAfy
γM1 (6.47) Reduktionsfaktorn χ kan antingen avläsas ur Figur 6.4 eller beräknas enligt formeln
χ = 1
Φ+√Φ2−λ2 (6.49) där hjälpparametern Φ skrivs som
Φ = 0,5[1 + α(λ − 0,2) + λ2]
30
Då (6.49) används avläses imperfektionsparametern α i Tabell 6.1 utifrån villkoren i Tabell 6.2.
31
Då flytspänningen fy = 620 MPa inte finns med i Tabell 6.2 används värdet för S 460 och knäckningskurva a används vilket ger α = 0,21.
Φ = 1,08 insatt i (6.49) ger χ = 0,67
Normalkraftskapaciteten för knäckning (6.47) blir Nb,Rd = 4,01 ∗ 103 kN