Numerisk prediktering av uppkomst av solkurvor — modellering ...
Numerisk prediktering av uppkomst av solkurvor — modellering ...
Numerisk prediktering av uppkomst av solkurvor — modellering ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
lateral styvhet<br />
[MN/m]<br />
<strong>Numerisk</strong> <strong>prediktering</strong> <strong>av</strong> <strong>uppkomst</strong> <strong>av</strong> <strong>solkurvor</strong><br />
vertikal styvhet<br />
[MN/m]<br />
longitudinell<br />
styvhet [MN/m]<br />
skalfaktor för<br />
flytgräns<br />
5.0 100.0 3.5 1<br />
Tabell 1: Materialparametrar i den användardefinierade ballastmodellen<br />
Den enda last som analyserats är temperaturbelastning, vilken relaterar till en axialkraft<br />
i rälen via temperaturutvidgningskoefficienten α = 11.7 · 10 −6 [C −1 ].<br />
3 Analys <strong>av</strong> resultat<br />
En linjär knäckningsanalys ger en kritisk temperatur <strong>av</strong> cirka 132 °C för referensfallet.<br />
Detta är naturligtvis ett resultat som är oanvändbart i praktiken. Anledningen till den<br />
höga kritiska temperaturen är den ideala spårgeometrin. Matematiskt relaterar en linjär<br />
knäckningsanalys till egenvärdet hos strukturens styvhetsmatris vilket förhindrar en<br />
ordentlig analys <strong>av</strong> inverkan <strong>av</strong> initiala deformationer, ytterligare laster (t.ex. fordonslaster),<br />
m.m. För att komma runt detta måste en icke-linjär gränslastanalys genomföras.<br />
I denna studie används en så kallad Riks-analys. Här behandlas magnituden <strong>av</strong> temperaturlasten<br />
som en obekant och är alltså inte den variabel som primärt styr analysen.<br />
Iterationen för att finna den numeriska lösningen styrs istället <strong>av</strong> en hjälpvariabel som<br />
i Abaqus benämns båglängd (arc length). Resultaten i form <strong>av</strong> förskjutningar, krafter,<br />
m.m. ges som funktioner <strong>av</strong> denna båglängd. För att relatera till en given temperatur<br />
måste därför båglängden kopplas till temperaturen i en båglängds–temperaturplot.<br />
Det resultat som söks i undersökningarna är gränslasten då de laterala deformationerna<br />
börjar att accelerera under ytterligare temperaturlast. Genom att plotta ett antal<br />
storheter visade det sig att denna gränslast var tydligast då sektionskrafterna i den mest<br />
förskjutna sliperns yttersta nod plottades. Man bör notera att denna sektionskrafts relation<br />
till t.ex. ballasttrycket mot sliperns kortända är oklar eftersom den beräknas i en<br />
nod till vilken en fjäder är kopplad.<br />
I figur 2 visas ett exempel på sektionskraften i den yttersta noden hos den mest förskjutna<br />
slipern (nod nummer 2383), samt den yttersta noden i grannslipern (nod nummer<br />
2372). Man ser tydligt en knyck som motsvarar ett mjuknande i sidomotståndet. Detta<br />
inträffar dock vid något olika temperaturer för de två sliprarna. Vår tolkning är att detta<br />
beror på att ballasten vid den mest förskjutna slipern mjuknar först. Detta leder till<br />
att den närliggande slipern får ta mer last. Då lasten uppgår till ballastens “flytgräns”<br />
kommer även denna slipers sidostyvhet att minska.<br />
Genom att plotta rältemperaturen mot båglängden, se figur 3 kan man relatera “knycken”<br />
i kurvan till en temperatur då sidostyvheten minskar. I det exempel som visas i<br />
figur 2 är dessa temperaturer 38.5 °C (nod 2383), respektive 40.4 °C (nod 2372). I denna<br />
rapport har den temperatur som motsvarar mjuknande i sidomotståndet för den mest<br />
deformerade slipern, betecknats kritisk temperatur.<br />
En graf <strong>av</strong> sliperförskjutningen mot båglängden visas i figur 3. Som nämndes ovan<br />
är “knycken” i förskjutningskurvorna inte lika markant som den för sektionskrafterna, se<br />
3