EXAMENSARBETE Skruvförbands inverkan på rotorer i vattenkraftverk Jimmy Norberg Westling 2015 Civilingenjörsexamen Maskinteknik Luleå tekniska universitet Institutionen för teknikvetenskap och matematik Förord Detta examensarbete är den avslutande delen för civilingenjörsutbildningen Maskinteknik med inriktning Konstruktion vid Luleå Tekniska Universitet. Jag vill tacka mina handledare Jan-Olov Aidanpää på Luleå Tekniska Universitet och Rolf Gustavsson på Vattenfall för all vägledning och råd under projektets gång. Jimmy Norberg Westling Luleå Sammanfattning I vattenkraftverk finns det stora turbiner som roterar på grund av vattnet som forsar förbi dem. Turbinerna är kopplade till generatorer som generar en stor del av elen i Sverige. Mellan turbinen och generatorn finns det flera skruvförband som håller alla delar fixerade. På grund av vattnets rörelse utsätts skruvförbanden för stora belastningar som gör att de böjs, vrids och sträcks ut. Belastningarna kan göra att flänsarna börjar separera ifrån varandra. Detta beteende gör att ett skruvförband inte är lika styvt som en solid motsvarighet. Trots denna svaghet används de flitigt i industrin framförallt för att de gör det möjligt att kunna byta ut och reparera eller uppgradera felande delar. Detta är en viktig egenskap för att underlätta underhåll och för att skapa lång livslängd på systemet. I denna rapport undersöks hur skruvförbandens böjstyvhet, vridstyvhet och axiella styvhet minskar med högre belastningar. Styvheten i skruvförbanden minskar när det belastas för att flänsarna separerar ifrån varandra. Man kan minska den här effekten genom att spänna åt skruvarna. Förspänningskraften gör att flänsen pressas ihop ytterligare vilket gör det svårare för skruvförbandets flänsar att separera ifrån varandra. I två av de studerade skruvförbanden, A och B, används skruvar som har en upphöjd yta som presspassas in vid flänsarnas skarv. Detta för att ta upp de stora vridmoment som uppstår när turbinen roterar. Skruvförband C har en annorlunda form, det sitter mellan vattenkraftverkets övre röraxel och turbinaxel. Den övre röraxelns diameter är dubbelt så stor som turbinaxeln vilket har gjort att skruvarna skruvas in på insidan av röraxeln. Skruvar med M90 gängor används för att hålla ihop detta skruvförband, för att ta upp vridmomenten används en hylsa istället som pressas emellan flänsskarven. För att ta reda på hur styvheten förändras med ökande moment skapades CAD-modeller ur ritningar, därefter gjordes finita element analyser på skruvförbanden. Eftersom det är flera delar som kan ta i varandra behövdes kontaktanalyser göras. De är betydligt mer tidskrävande analyser än vanliga FEM analyser vilket gav motivation att finna förenklade modeller som ger liknande svar som kontakanalyserna. Framförallt har tre förenklade modeller studerats: Solid modell där flänsarna har avfasats med 30° och 45°. Solid modell där flänsarna har tagits bort så att modellen liknar en cylinder. Solid modell där flänsarnas E-modul har minskats med en konstant. Kontaktanalyserna visade att skruvförbanden har en konstant böjstyvhet tills de böjande momenten överstiger förspänningskraften. Därefter börjar böjstyvheten att minska, med extrema böjande moment, förutsatt att skruvarna inte redan har gått av, börjar de inre ytorna att komma i kontakt med varandra vilket gör att styvheten ökar igen och konvergerar mot ett värde eftersom kontaktytorna slutar att förändras. Vridstyvheten är konstant för små vridande moment på grund av friktionen men när belastningarna höjs börjar flänsytorna att glida relativt varandra vilket gör att vridstyvheten sjunker. Skruvarna kommer därefter att komma i kontakt med skruvhålen vilket gör att vridstyvheten ökar en aning och konvergerar till en vridstyvhet där skruvarna tar upp majoriteten av de vridande momenten. Den axiella styvheten är konstant tills de axiella krafterna överstiger förspänningskraften i skruvarna. Därefter börjar de inre ytorna att separera vilket gör att den axiella styvheten minskar tills en stor del flänsarna har helt separerat ifrån varandra. Att skruvförbandens styvhet är konstant vid små belastningar styrker användandet av förenklade modeller eftersom deras styvheter inte förändras med högre belastningar. Men att fasa av flänsarna gav styvheter som var anmärkningsvärt högre än kontaktanalysen vid små belastningar. Detta berodde till stor del på att det fanns mellanrum i flänsarna som minskade styvheten. Detta innebär att de förenklade modellerna också måste ha dessa mellanrum för att få jämförbara resultat som kontaktanalyserna. Att minska E-modulen i den förenklade modellens flänsar med en konstant har gett jämförbara resultat som kontaktanalysen förutsatt att modellen stämmer bra överrens med kontaktanalysen vid små belastningar. Detta öppnar upp möjligheten att variera E-modulen i modellens flänsar för att simulera styvhetsminskningen med betydligt lägre beräkningstider som resultat. Abstract Inside hydropower plants there are large turbines that rotate due to the flow of the water. The turbines are connected to generators, between the turbine and the generator there are several bolted joints that keeps all the parts fixed. Due to the water flow the bolted joints are subjected to large loads which deform the bolted joints. If the loads are large enough they may force the flanges to separate from each other. This behavior implies that bolted joints aren’t as stiff as a solid counterpart. Still, the bolted joints are frequently used in the industry, mainly because they make it possible to repair and switch faulty parts with ease. This is an important feature to generate a long life-cycle of hydropower plants. In this report the bolted joints varying stiffness is studied as the bolted joints are subjected to bending moments, torques and axial forces. The bolted joints stiffness will decrease as it is subjected to larger loads due to the flange separation. The decrease can be halted by adding a pre-load to the bolts. The pre-load further deforms the flanges which means larger external loads are need to force the flanges to separate. Three different bolted joints have been studied in this report. The bolted joints in two of them, Bolted joint A and B, have a similar type of bolted joint. The bolts used in these bolted joints have an uncommon feature; an elevated surface that is fitted inside the bolt holes. The reason for this elevated surface is to make sure that flange surfaces are sliding due to the large torques from the rotating turbine. Bolted joint C has a different shape because the upper tube shafts diameter is considerably larger than the turbine shafts diameter, forcing the bolts to be fastened on the inside of the upper tube shaft. Bolts without elevated surfaces with screw thread size M90 are used in this bolted joint however to account for the torque loads a case is fitted inside bolt holes. To determine how the stiffness varies with higher external moments several models were created in a CAEprogram based on drawings from the Bolted joint A, B and C followed by finite element analysis. Since there are several parts that can come in contact with each other contact analysis was required. Contact analysis takes considerably longer time to solve compared to regular FEM operations. Therefore simplified models were also created that gives similar results to the contact analysis. Mainly three models have been studied: Solid model where the flanges have been chamfered by 30° and 45°. Solid model where the flanges have been removed with a similar form as a cylinder. Solid model where the flanges’ Young’ modulus has been a constant. From the contact analysis it was shown that the bolted joints have a constant bending stiffness until the bending moments are larger than the pre-load. As the bending moment gets larger the bending stiffness becomes lower. Eventually with extreme bending moments, granted the bolts have not already broken, the inner surfaces will come in contact with each other which leads to an increase of the bending stiffness and will converge to a value because the surfaces in contact stops varying. The torsional stiffness is constant at small torques due to friction between the flanges but as the torque increases the flange surfaces starts to slide relative to each other results in the torsional stiffness to decrease. Eventually the bolts will come in contact with the bolt holes leading to an increase in torsional stiffness and will converge to torsional stiffness where the bolts are subjected to the majority of the torque. The axial stiffness is also constant until the external axial forces exceed the pre-load in the screws. The axial stiffness decreases as the inner surfaces separates due to the increasing axial forces until a large part of contact surfaces have separated. Since the stiffness varies very little at lower loads it further strengthens the use of simplified models because their stiffness does not vary. Howerer, the simplified analysis with a beveled solid flange resulted in stiffness values considerably larger than the contact analysis. The main part of the reason is the spaces that are between the flanges. This means it is important to take these spaces into account when modelling a simplified model to get accurate results. By lowering the Young’s modulus in the simplified analysis by a constant has given comparable results with the contact analysis, given the model has a similar stiffness with the original Young’s modulus as the contact analysis at smaller loads. This means it’s possible to simulate the stiffness variations by varying the Young’s modulus, with this method the computation time can be decrease considerably. Innehållsförteckning 1 Inledning ........................................................................................................................................................ 12 1.1 Bakgrund ............................................................................................................................................... 12 1.2 Problembeskrivning .............................................................................................................................. 13 1.3 Syfte ...................................................................................................................................................... 13 1.4 Avgränsningar ....................................................................................................................................... 13 2 Metod ............................................................................................................................................................. 14 3 Teori............................................................................................................................................................... 16 3.1 Deformation .......................................................................................................................................... 16 3.2 Materials elastiska och plastiska egenskaper ........................................................................................ 16 3.2.1 3.3 Förspänning i ett skruvförband ............................................................................................................. 17 3.4 Böjning av ett skruvförband .................................................................................................................. 18 3.4.1 Elementarfall för fritt upplagd balk .............................................................................................. 18 3.4.2 Böjning av ett symmetriskt skruvförband med elastisk axel ........................................................ 19 3.5 Finita Element Metoden ........................................................................................................................ 21 3.5.1 4 Styvhet .......................................................................................................................................... 17 Kontaktanalys ............................................................................................................................... 21 Resultat .......................................................................................................................................................... 22 4.1 Skruvförband A ..................................................................................................................................... 22 4.1.1 Böjstyvhet för skruvförband A ..................................................................................................... 23 4.1.2 Vridstyvhet för skruvförband A ................................................................................................... 31 4.1.3 Axiell styvhet i skruvförband A ................................................................................................... 33 4.2 Skruvförband B ..................................................................................................................................... 35 4.2.1 Böjstyvhet för skruvförband B ..................................................................................................... 36 4.2.2 Vridstyvhet för skruvförband B.................................................................................................... 42 4.2.3 Axiell styvhet i skruvförband B ................................................................................................... 44 4.3 Skruvförband C ..................................................................................................................................... 46 4.3.1 Böjstyvhet för skruvförband C ..................................................................................................... 47 4.3.2 Vridstyvhet för skruvförband C.................................................................................................... 52 4.3.3 Axiell styvhet i skruvförband C ................................................................................................... 54 5 Slutsatser ........................................................................................................................................................ 56 6 Diskussion...................................................................................................................................................... 58 7 Referenser ...................................................................................................................................................... 59 8 Appendix........................................................................................................................................................ 61 Tabell 1. Nomenklatur Symbol 𝑨 𝑨𝒌 𝑨𝒔 𝑫 𝑬 𝑬𝒌 𝑬𝒔 𝑭 𝑭𝒑𝒓𝒆 𝑭𝒔 𝑰 𝑴 𝑴𝒄𝒚𝒍 𝑷 𝒅 𝒅𝒉 𝒅𝒘 𝒌 𝒌𝒌 𝒌𝒔 𝒌𝒕𝒐𝒕 𝒍 𝒍𝟎 𝒍𝒄𝒚𝒍 𝒍𝒉 𝒍𝒌 𝒍𝒔 𝒏 𝒓 𝒕 𝜶 𝜷 𝜷𝟏 𝜷𝟐 𝜷𝒌 𝜸 𝜹 𝜹𝒃 𝝐 𝜽𝒃 𝝈 Förklaring Area Arean för de klämda delarna Nominella spänningsarean för skruven Ytterdiametern på axeln E-modul E-modul för de klämda delarna E-modul för skruven Axiell kraft Förspänningskraft Skruvkraft Böjtröghetsmoment för en tjockväggig cylinder Moment vid skruvförbandets ändar Moment för elementarfallet för en fritt upplagd cylinder Längd från flänscentrum till separationspunkten Innerdiametern på axeln Hålets diameter Diameter för den plana ytan under skruvhuvdet Styvhet Styvhet för de klämda delarna Styvhet för skruven Ekvivalent styvhet för 𝑘1 , 𝑘2 , … , 𝑘𝑛 Totala längden för ena halvan av skruvförbandet Begynnelselängd Längden för cylindern Längden från flänscentrum till skruvhålen i y-led Längden för de klämda delarna Längden för skruven Totalt antal skruvar Skruvhålsradie Flänsens tjocklek Vinkel mellan skruvhålen Böjningsvinkel Vinkel för vänstra flänsen Vinkel för högra flänsen Halva konvinkeln, sätts till 31° Flänsens separationsvinkel Materialdeformation Deformation för en konsolbalk Materialtöjning Böjningsvinkeln för en konsolbalk Spänning Enhet m 𝑚2 𝑚2 m Pa Pa Pa N N N 𝑚4 Nm Nm m m 𝑚2 m N/m N/m N/m N/m m m m m m m m m ° ° ° ° ° ° m m ° Pa 1 Inledning Inom vattenkraften förekommer ett stort antal varianter av skruvförband där skruvarna har en låg förspänning. Ett skruvförband består av en skruv, mutter och en eventuell bricka som håller fast delarna i förbandet, förspänningen uppstår när muttern spänner åt skruven. En skruv är definitionsmässigt en utvändigt gängad cylindrisk kropp [1]. Skruvförbanden som studeras i denna rapport kopplar ihop axlarna mellan turbinen och generatorn eller den övre röraxeln. Flänsarna har en diameter på ca 1,2 − 4 𝑚 och skruvarna som håller ihop skruvförbanden väger ca 35 𝑘𝑔. Typiska belastningar som dessa skruvförband utsätts för under drift är böjande, vridande och axiella belastningar. Detta innebär att om skruvarna inte är tillräckligt förspända kan flänsarna börja glappa. Glappar skruvförbandet kommer det att kunna deformera sig mer än om skruvförbandet inte varit där vilket kan leda till haveri om inte konstruktören har haft detta i åtanke när vattenkraftverket designades. 1.1 Bakgrund Skruvförband är i sin enklaste form två delar som hålls ihop av skruvar och har stor användning i dagens samhälle. Användningen av gängor kan spåras tillbaka till 400 f. kr. men de hade ofta andra funktioner än bara att fixera delar. Skruvar återfinns i rustningar från 1400-talet, de hölls ihop muttrar. Den första maskinen för att tillverka skruvar gjordes av Besson i Frankrike 1568. Senare vid 1641 förbättrade Hindley of York-maskinen och den började användas allt mer. Under första och andra världskriget var bristen på standardiserade gängor ett stort bekymmer vilket gjorde att Storbritannien, USA och Canada beslutade att använda en universell gängstandard för alla länder som använder det brittiska måttsystemet [2]. I Sverige skapades nationella standarder på 1920-talet bland annat på gängor, skruvar muttrar och fästelement. Under 1970-talet hade Sverige ett stort inflytande över gängutvecklingen och lyckades förena tyska, ryska och amerikanska särintressen till internationellt system av metriska gängor och gängtoleranser. Resultatet är den ISO-standard som är vida använd i dagens samhälle. Metrisk gänga är numera en av de vanligaste, men tumgängor används fortfarande i hög utsträckning [3]. Skruvförbands beteenden har diskuterats i olika artiklar och böcker. Waters et. Als [4] forskning och Schneiders fortsatta forskning i artikeln ”Flat Face Flanges With Metal-to-Metal Contact Beyond the Bolt Circle” [5] ligger till grunden för många av The American Society of Mechanical Engineers standarder när det gäller skruvförbandsdesign. I LTUs Maskinelement bok [1] diskuteras hur förspänningen i ett skruvförband kan bestämmas utifrån en känd axiell kraft, även sättning nämns och hur man kan kompensera för det. Bouzid har skrivit flera artiklar angående skruvförband. I artikeln ”Analytical Modeling of Flat Face Flanges With Metal-to-Metal Contact Beyond the Bolt Circle” [6] har två analytiska modeller tagits fram för skruvförband som har metall-metall kontakt utanför skruvcirkeln. Modellen antar dock att det är bara ett symmetriskt innertryck som påverkar själva deformationen. I artikeln ”On the Effect of External Bending Loads in Bolted Flange Joints” [7] har Bouzid föreslagit hur externa moment ska behandlas i tätade skruvförband och visar hur externa moment kan göras om till ekvivalenta axiella krafter som föreslogs av Koves i artikeln ”Analysis of Flange Joints Under External Loads” [8], dessa axiella krafter kan då behandlas på ett liknade sätt som innertrycket i [6]. Men detta är för tätade skruvförband, detta antagande gör stor skillnad för kontaktpunkten antas vara på själva tätningen som ligger innanför bultcirkeln. I artikeln ”Unimportance of geometric nonlinearity in analysis of flanged joints with metal-to-metal contact” [9] visade Joshi et. Al. att skruvförbands olinjära beteende inte har så stor inverkan när belastningarna är små. De gjorde undersökningen genom att utsätta skruvförbandsmodellen för ett innertryck och en axiell kraft. De visade att mycket tid kan sparas genom att göra finita element analyserna utan kontaktvillkor, dock undersöktes bara hur axiella krafter påverkade flänsseparationen. 12 1.2 Problembeskrivning Inom hållfasthetsläran används många olika elementarfall för att beskriva hur balkar, stänger och cylindrar. beter sig under belastning [10]. Dessa modeller är väl beprövade och används till många applikationer. Skruvförband används i liknande miljöer men eftersom det är två delar som kopplas ihop får skruvförbanden en lägre styvhet jämfört med en solid kropp när den böjs. Vattenfall har i sina vattenkraftsrotorer stora skruvförband som kopplar ihop deras generatorer med turbinerna och vill därför veta hur skruvförbanden påverkar styvheten i rotorn och hur väl olika förenklade modeller som används stämmer överens med varandra. 1.3 Syfte Om skruvförbanden sitter löst kommer deras egenskaper att ändras betydligt eftersom skruvförbanden kan dela sig vid kopplingen. Detta fenomen är därför intressant att analysera för att se hur stor skillnad lösa och hårt förspända skruvförband påverkar styvheten. Inom vattenkraften förekommer olika förenklade modeller för att minska styvheten i skruvförbanden vilket bör utvärderas. Exempelvis finns det förenklade solida modeller där flänsarna har fasats av med 30° och 45° på grund av geometriförändringen som sker vid flänsarna. En annan förenklad modell är att skapa en solid modell och sedan sänka E-modulen i flänsarna för att efterlikna flänsseparationen. Detta gav följande forskningsfrågor att undersöka: Hur ändras styvheten i ett skruvförband när det utsätts för böjande moment, vridande moment och axiella krafter? Hur väl stämmer olika förenklade modeller överens med varandra? 1.4 Avgränsningar Projektet påbörjades 2015-01-20, en muntlig presentation skedde 2015-05-29. Projektet kommer i detalj att studera styvhetsförändring när skruvförbanden utsätts för böjande moment. Analytiska modeller kommer bara att avse stela skruvförband eller med hjälp av kända elementarfall. För beräkningar används MATLAB och CAE-modeller kommer att göras i Siemens NX 10. Modellerna innehåller enbart nödvändiga geometrier för att uppfylla skruvförbandets funktion, hela axlarna som skruvförbanden sitter i modelleras inte, istället kapas axlarna 640 mm från kopplingen detta för att minska beräkningstider. Skruvar och muttrar modelleras som en solid, utan gängor och andra avancerade geometrier. FEM-analyserna tar inte plasticering i åtanke, detta för att kunna se om styvheten konvergerar till ett värde och om det går att se ett tydligt samband. Kunskapen kommer däremot bara att kunna appliceras där alla delar i skruvförbandet beter sig elastiskt. 13 2 Metod I Figur 1 visas projektets tillvägagångsätt. Figur 1. Examensarbetets arbetsgång I början av projektet hölls ett möte tillsammans med handledare från LTU och Vattenfall där projektets problembeskrivning förklarades. Syfte, mål och avgränsningar formulerades utifrån problembeskrivningen. Dessa delar har ändrats och förbättrats med tiden med nya erfarenheter. Därefter gjordes en övergripande planering av projektet med tidsplan. Utifrån problembeskrivningen måste ett beslut tas hur undersökningen skall göras. De vanligaste undersökningsansatserna är survey-undersökning, fallstudie och experiment [11]. Survey innebär att undersökningen görs på en större avgränsad grupp med hjälp av frågeformulär eller intervjuer. Fallstudie innebär att undersökningen görs på en mindre avgränsad grupp. Experiment är en beteckning på en undersökningsuppläggning där några enstaka variabler undersöks. Survey och fallstudie är dock svåra att applicera på naturvetenskapliga problem därmed valdes en experimentell undersökningsansats. Den beroende variabeln i detta projekt är styvheten och de oberoende variablerna är externa belastningar och 14 förspänningskrafter, när förenklade modeller har testats har även deras geometri varierats för att se hur avancerad geometri som krävs för att ge liknande styvheter. En opponering har även gjorts på ett tidigare examensarbete, detta för att få en förståelse av vad som ska vara med i ett examensarbete. Opponeringen gjordes på Lisa Kvarnströms examensarbete ”Förspänningsförluster i maskinfundamentens skruvförband - en jämförelse mellan praktik och teori” [12]. Rapporten handlar om hur förspänningen ändras med tiden på ett skruvförband som sitter på ett maskinfundament, maskinfundament är oftast tillverkade i betong. Att rapporten handlar om skruvförband och hur de förlängs med tiden hade en viss relevans för examensarbetet. Litteraturstudier har utförts för att få en bättre förståelse av skruvförbands beteende. Under litteraturstudierna har ritningar på diverse skruvförband som finns i vattenkraftverk runt om i Sverige studerats. Även böcker och artiklar angående skruvförband, skruvar och rotorer har studerats. I vissa artiklar som har studerats har olika slarvfel i härledningar framkommit vilket visar hur viktigt det är att bekräfta referenser och beräkningar. I detta projekt har framförallt en kvantitativ analysmetod använts, det innebär att numeriska värden från mätningar har bearbetats och analyserats [11]. Utifrån ritningarna har modeller på skruvförbanden skapats med hjälp av CAE-programmet NX 10 [13] för att kunna studera hur olika styvhetsvärden förändras med varierande förspänningar, belastningar och geometrier. Eftersom det är flera delar som kan komma i kontakt med varandra har kontaktanalyser behövts göras och denna typ av finita element analyser är betydligt mer tidskrävande. Olika förenklade modeller togs även fram för att jämföra med kontaktanalyserna. Kontaktanalyser tar längre tid än vanliga därför fanns behovet att hitta förenklade modeller som ger liknande svar som kontaktanalysen. Analytiska modeller som beskriver böjningen på ett acceptabelt sätt togs även fram. Om analytiska modeller stämmer väl överens förkortas framtagningstiden än mer om man jämför med förenklade FEM-modeller. Intiellt antogs att skruvförbandet var stelt men detta gav inte tillräcklig noggrannhet jämfört med de finita element analyserna. Istället antogs att axlarna var flexibla vilket gav jämförbara resultat med flänsar som inte har mellanrum i flänsarna. Ur resultaten från de finita element analyserna och analytiska analyserna kunde slutsatser dras. Det finns tre angreppssätt att relatera teori och empiri och att kunna dra slutsatser: deduktion, induktion, och abduktion [11]. Att arbeta deduktivt innebär att man utifrån allmänna principer och teorier drar slutsatser. En fara med detta arbetssätt är att befintliga teorier kan rikta forskningen så att så att nya intressanta rön inte upptäcks. Att arbeta induktivt innebär att man följer empirin och formulerar utifrån de slutsatser utan någon tidigare vedertagen teori som stöd. Risken finns att man inte vet hur lång räckvidd slutsatserna har eftersom det sällan finns empiriska värden för alla tänkbara situationer. Abduktion är en kombination av deduktion och induktion. Det innebär att man utifrån enstaka fall formulerar teorier, vilket kännetecknas av att vara induktiv. Nästa steg testas teorin på ett nytt fall och utvecklas vidare, vilket kännetecknas att vara deduktiv. Risken finns att nya fall väljs utifrån tidigare erfarenheter och formulerar hypotetiska teorier som utesluter andra tolkningar. Inom naturvetenskapen används framförallt ett deduktivt angreppssätt eftersom de består av redan sammanhängande lagar och är empiriskt prövade. Detta projekt har framförallt haft ett deduktivt arbetssätt. Data har tagits fram från etablerade finita element metoder och analyserats, förenklade modeller har tagits fram som jämförs med FEM-modellerna. För att öka validiteten har flera skruvförband undersöks med olika geometrier för att se kunna se liknande beteenden och för att öka reliabiliteten i detta projekt har värdena jämförts med väl etablerade elementarfall för att säkerställa att styvheten är inom rimliga gränser Examensarbetet avslutades med en presentation på Luleå Tekniska Universitet. 15 3 Teori I detta kapitel presenteras olika teorier som används i rapporten såsom deformation, elasticitet och styvhet. 3.1 Deformation Om ett material utsätts för belastningar kommer materialet att börja deformera. Detta beteende är lätt att se och testa med gummi men kan appliceras på de flesta material. Sambandet mellan spänning och deformation kan beskrivas med Hooke’s lag [14] med följande samband (3.1) 𝜎 = 𝐸𝜖 Där 𝜎 är spänningen, 𝐸 är materialets E-modul och 𝜖 är materialets töjning vilket kan skrivas som 𝛿 𝜖= 𝑙0 (3.2) Där 𝛿 är deformationen och 𝑙0 är den ursprungliga längden. 3.2 Materials elastiska och plastiska egenskaper De vanligt förkommande materialen som används i industrier har elastiska egenskaper, till exempel stål, gjutjärn eller aluminium [10]. Ett elastiskt material innebär efter att materialet har deformerats så återgår materialet till sitt ursprungliga tillstånd. Ekvation (3.1) visade att det är ett linjärt samband mellan spänning och töjning. Detta samband gäller dock bara upp till materialets sträckgräns, 𝜎𝑠 . I Figur 2 kan man se sambandet mellan spänning och töjning [14]. Figur 2. Spänning- och töjningsdiagram Om 𝜎 > 𝜎𝑠 kommer materialet att plasticera vilket innebär att materialet inte återgår till sitt ursprungliga tillstånd efter avlastning. När spänningen når brottgränsen, 𝜎 = 𝜎𝑏 , går materialet av. De olika faserna elasticitet, plasticering och brott är enkelt att exemplifiera och testa med mjukplast eftersom de har en låg E-modul som gör att det inte krävs en stor belastning för att töja ut plasten. Belastas plasten lite töjs den ut en aning men återgår till sin ursprungliga form, med högre belastningar upptäcker man att den fortfarande har deformeringar efter belastning och med ännu högre belastningar går plasten sönder. 16 3.2.1 Styvhet Styvhet beskriver hur motståndskraftigt ett material är mot deformation när den utsätts för laster. Styvhet påminner mycket om E-modulen men är även beroende av kroppens geometri [14]. Styvheten kan härledas ur Hooke’s lag 𝐸𝐴 𝐹 = 𝜎𝐴 = 𝐸𝜖𝐴 = 𝛿 = 𝑘𝛿 (3.3) 𝑙0 Där 𝑘 är styvheten, 𝛿 är deformationen. Det är värt att notera att (3.3) gäller för enaxliga fall. Om fallet har fler frihetsgrader krävs fler samband. 3.3 Förspänning i ett skruvförband Skruvförband förspänns med en viss kraft för att säkerställa att de inte går lös, läcker eller separerar. Förspänningen uppstår när skruven och muttern dras åt. Kraftflödet som uppstår i förbandet när skruven förspänns kan ses i Figur 3. Figur 3. Kraftflöde i förbandet, tagen ur [1] I Figur 3 ser man att när klämkrafterna ökar fördelas spänningszonen i form av en dubbelkon. Detta gäller om de hopklämda delarna liknar två plåtar som ligger dikt an mot varandra. Därför görs ofta beräkningar med elastiska kroppar som liknar dubbelkoner [1]. Förspänningskraften, 𝐹𝑝𝑟𝑒 , går att bestämma om man antar skruvförbandet inte får glappa när den utsätts för en axiell kraft 𝐹. Styvheten för skruven och de klämda delarna är: 𝑘𝑠 = 𝐸𝑠 𝐴𝑠 𝑙𝑠 , 𝑘𝑘 = 𝐸𝑘 𝐴𝑘 𝑙𝑘 (3.4) Där 𝐸𝑠 är skruvens E-modul, 𝑙𝑠 är den töjbara delen av skruvstammen vilket är längden mellan mutter och skruvhuvud. 𝐴𝑠 är skruvens nominella spänningsarea och finns ofta tabellerad. 𝐸𝑘 är de klämda delarnas Emodul, 𝐿𝑘 är de bägge flänsarnas tjocklek. 17 För två hopfogade plåtar kan den klämda zonens area, 𝐴𝑘 , beräknas enligt följande: 𝜋𝑑ℎ 𝑙𝑘 𝑡𝑎𝑛𝛽𝑘 𝐴𝑘 = 𝑙 𝑡𝑎𝑛𝛽𝑘 1+ 𝑘 𝑑𝑤 − 𝑑ℎ 2 ln ( ) 𝑙 𝑡𝑎𝑛𝛽𝑘 1+ 𝑘 𝑑𝑤 + 𝑑ℎ (3.5) Där 𝑑ℎ är hålets diameter, 𝑑𝑤 är diametern för den plana ytan under skruvhuvudet som klämmer mot flänsen och 𝛽𝑘 är halva konvinkeln, denna sätts enligt [1] till 31°. Figur 4. Friläggning av ett skruvförband tagen ur [1] Utifrån Figur 4 kan förspänningskraften som krävs för att ingen separation skall uppstå bestämmas med 𝐹 𝐹𝑝𝑟𝑒 = (3.6) 𝑘 1+ 𝑠 𝑘𝑘 3.4 Böjning av ett skruvförband När ett skruvförband böjs kommer flänsarna tillslut att separera ifrån varandra. Detta leder till att förbandet blir lättare att böja ju mer flänsarna separerar ifrån varandra. Att finna förenklade modeller som beskriver böjningen och styvheten under vanliga belastningar är därmed intressant att undersöka. Nedan presenteras två analytiska modeller som ger liknande böjstyvheter som ett skruvförband. 3.4.1 Elementarfall för fritt upplagd balk Ett sätt att minska böjstyvheten är att helt ta bort flänsarna så att denna förenklade modell liknar en cylinder. För ett elementarfall för fritt upplagda balkar kan sambandet mellan böjningsmoment och böjningsvinkel vid ändarna beskrivas med följande samband [10]: 2𝐸𝐼 𝑀𝑐𝑦𝑙 = 𝛽 (3.7) 𝑙𝑐𝑦𝑙 Där 𝐸 är materialets E-modul, 𝑙 är balkens längd, 𝛽 är böjningsvinkeln vid ändarna och för en tjockväggig cylinder är böjtröghetsmomentet [15] 𝜋 𝐼= (𝐷4 − 𝑑 4 ) (3.8) 64 18 3.4.2 Böjning av ett symmetriskt skruvförband med elastisk axel Elementarfallet i 3.4.1 tar inte förspänningskrafter i åtanke. Därför gjordes det till en början en modell för ett stelt skruvförband, härledningarna kan ses i Appendix C. Det visade sig dock att en styv axel inte är tillräckligt för att beskriva böjningen på ett acceptabelt sätt. Istället antogs att axeln är elastisk, flänsen är styv och att den endra halvan har samma dimensioner, se Figur 5. Figur 5. Friläggning för ena halvan av flänsen För en konsolbalk kan deformationen vid änden när den utsätts för ett moment skrivas med följande samband [10] 𝑀(𝑙 − 𝑡)2 (3.9) 𝛿𝑏 = 2𝐸𝐼 Och vinkeln vid änden beräknas med 𝜃𝑏 = 𝑀𝑙 𝐸𝐼 Längden från flänsens centrum till skruvhålen är 𝑙ℎ𝑖 = 𝑟 𝑐𝑜𝑠𝛼 (3.10) (3.11) Skruvarna har skruvats åt med en förspänning 𝐹𝑝𝑟𝑒 . När flänsarna separerar i punkten 𝑃 kommer skruven att deformeras vilket ökar skruvkraften betydligt. 19 Figur 6. Skruvkrafter I Figur 6 har antagandet att skruvkrafterna alltid är vågräta gjorts, detta stämmer bra om flänsarna är symmetriska. Antagandet faller om flänsen skulle vara osymmetrisk men eftersom separationen är generellt väldigt små så påverkar det resultatet väldigt lite i osymmetriska fall också. Skruvkrafterna blir 𝐸𝑠 𝐴𝑠 (3.12) 𝐹𝑠 𝑖 = 𝐹𝑝𝑟𝑒 + 2 (𝑃 − 𝑙ℎ i ) sin 𝛾 𝑡 Krafterna som uppstår från deformationen multipliceras med 2 för att skruven deformeras lika mycket åt andra hållet på grund av att flänsarna är symmetriska. 𝑃 är längden från flänscentrum till punkten där flänsarna separerar. I Figur 6 är 𝑃 = 𝑅. Om separationspunkten inte är längst ut på flänsen, utan är till exempel vid 𝑃 = 𝑅/2. Då finns risken att några av skruvarna är utanför separationspunken, 𝑃 < 𝑙ℎ 𝑖 . De uteliggande skruvarna kommer inte att deformeras alls, vilket innebär att 𝐹𝑠 𝑖 = 𝐹𝑝𝑟𝑒 vid de platserna. Momentjämvikten runt punkt 1 för Figur 5 blir ↶ 1 𝑛 0 = 𝑀 + ∑ 𝐹𝑠 𝑖 (𝑙ℎ 𝑖 𝑐𝑜𝑠𝛾 − (𝑙𝑠𝑖𝑛𝛾 + 𝛿𝑐𝑜𝑠𝛾)) (3.13) 𝑖=1 Substituerar (3.9) i (3.13) ↶ 1 𝑛 0 = 𝑀 + ∑ 𝐹𝑠 𝑖 (𝑙ℎ 𝑖 𝑐𝑜𝑠𝛾 − (𝑙𝑠𝑖𝑛𝛾 + 𝑖=1 𝑀𝑙 2 𝑐𝑜𝑠𝛾)) 2𝐸𝐼 (3.14) och löser för 𝑀. 𝑀= − ∑𝑛𝑖=1 𝐹𝑠 𝑖 (𝑙ℎ 𝑖 𝑐𝑜𝑠𝛾 − 𝑙𝑠𝑖𝑛𝛾) (𝑙 − 𝑡)2 1 − ∑𝑛𝑖=1 𝐹𝑠 𝑖 ( 𝑐𝑜𝑠𝛾) 2𝐸𝐼 Vinkeln som uppstår i ändarna, 𝛽, beräknas genom att addera 𝛾 och (3.10) vilket blir 𝑀(𝑙 − 𝑡) 𝛽=𝛾+ 𝐸𝐼 20 (3.15) (3.16) 3.5 Finita Element Metoden Finita element metoden används för avancerade analyser inom hållfasthetsläran. Metoden går ut på att man delar upp modellen i mindre delområden, finita element, sedan använder hållfasthetslärans ekvationer för att lösa ut deformeringar, spänningar och andra storheter för varje element. Det finns många typer av element som stångelement, balkelement och 3D-element. Elementen består av noder, hur många den består av beror på vad för typ elementet är. Till exempel har ett stångelement två noder. Varje nod kan röra sig i tre riktningar x, y, z och även rotera runt x-, y- och z-axlarna. Detta innebär att noderna har 6 stycken frihetsgrader, 3 translationsrörelser och 3 rotationer [14]. Elementen kan låsa nodernas olika frihetsgrader beroende vad för typ elementet är. Under definieringen av modellen måste frihetsgrader låsas för att kunna lösa simuleringen. Anledningen är enkel att förstå om man tänker att en stångmodell ligger svävande i rymden. Ponera att stången utsätts för en axiell kraft för att se hur mycket stången deformeras. Om inte någon nod har låsts i samma riktning som kraften kommer stången att försvinna iväg i intet, i samma riktning som kraften för att det inte finns något som håller stången på plats. Detta exempel visar hur viktigt det är att låsa frihetsgrader för att få korrekta resultat. 3.5.1 Kontaktanalys Om det finns en risk att olika ytor kan komma i kontakt med varandra när finita element analyser pågår måste kontaktvillkor definieras mellan ytorna. Kontaktvillkoren förhindrar ytorna från penetrering och tillåter glidning med valfria friktionsförluster [16]. Till exempel i skruvförband finns det ett flertal ställen där delar kan komma eller redan är i kontakt med varandra: Skruvhuvudet mot flänsen Skruvgängor mot flänshålen Skruvgängor mot mutter Mutter mot fläns Fläns mot fläns FEM-programmet upptäcker att beräkningsnätet kommer i kontakt med ett annat beräkningsnät, reaktionskrafter uppstår på grund av kontakten. Lösningen itererar tills kontaktkraftsförändringarna mellan kropparna har börjat konvergera [16]. 21 4 Resultat Detta kapitel redovisar resultaten som framkommit ur de finita element analyser och analytiska analyser som gjorts under projektets gång. Finita element analyser har gjorts i NX 10 för att se hur skruvförbandet böjs när det utsätts för externa belastningar. Eftersom det är flera delar som inverkar på varandra har kontaktanalyser gjorts för att se hur styvheten ändras. Därefter har olika förenklade modeller gjorts för att jämföra med kontaktanalysen. Alla delar har ståls materialegenskaper med E-modul 206,94 𝑀𝑃𝑎, Poissons konstant 0,3 och friktionskoefficienten har satts till 0,1 om inte annat anges. Skruvförbanden som redovisas har kapats 640 𝑚𝑚 från koppling för att underlätta analysen. 4.1 Skruvförband A Detta är ett skruvförband som har skruvarna på yttersidan om axeln. Se Appendix A.1 för översiktsritningar. 14 stycken skruvar håller skruvförbandet fixerat och är designade för att kunna ta upp vridmomenten som uppstår. Skruvarna har en upphöjd yta vid kopplingen som presspassas in i hålen på flänsarna, se Figur 7. Den upphöjda ytan hamnar vid flänsskarven när skruven är i rätt position, detta för att flänsarna inte ska kunna vridas om friktionskrafterna mellan flänsarna inte är tillräckligt stor för att hålla flänsarna på plats. Skruvarna är 460 𝑚𝑚 lång och diametern vid gängorna är 105 𝑚𝑚 och skruven väger ungefär 35 𝑘𝑔. Figur 7. Skruv som används i skruvförband A Skruvarna ska förspännas till 1200 𝑘𝑝/𝑐𝑚2 enligt ritningen vilket motsvarar 117 𝑀𝑃𝑎. Modellen som har gjorts utifrån skruvförbandets ritningar kan ses i Figur 8. Figur 8. CAD-modell av skruvförband A Flänsarnas diameter är 1170 mm och axlarna 780 mm i diameter på den vänstra axeln och 720 mm i diameter på den högra axeln. Eftersom geometrierna skiljer sig från varandra, är ena sidan enklare att deformera jämfört med den andra. 22 4.1.1 Böjstyvhet för skruvförband A Eftersom det är två flänsar som hålls ihop av presspassade skruvar behövdes en kontaktanalys göras för alla delar som kan eller tar i varandra. Ytorna mellan flänsarna, skruv och fläns definierades som kontaktytor. Muttrarna och skruvhuvudet ska inte kunna rör sig därför limmades (surface-to-surface glue) de fast på flänsytorna. Den generella arbetsgången i NX 10 förklaras detaljerat i Appendix B. Elementtypen CTETRA(10) som har använts under analyserna är ett solidelement, det används ofta för att modellera komplicerade system med många avrundningar och avancerade geometrier [16]. Detta solidelement tillåter inte att noderna roterar vilket innebär att nya punkter behövde införas i ändarnas centrumpunkter för att kunna använda momentbelastningarna som finns i NX 10. Sedan kopplades skruvförbandets ändar ihop med centrumpunkterna med hjälp av stela element. Dessa punkter användes också till att låsa frihetsgraderna för translationsrörelserna som kan ses i Figur 9. Figur 9. Beräkningsnät för skruvförband A Skruvförbandet utsattes för två lika stora böjande moment runt z-axeln vid ändarna på skruvförbandet vilket gör att skruvförbandet böjer sig enligt Figur 10. Figur 10. Böjningsvinklarna som studerades när skruvförbanden utsattes för ett böjande moment 23 I Figur 11 har momentbelastningen varierats för att se hur vinklarna 𝛽1 och 𝛽2 ändrades med ökat moment. Figur 11. Förändring av böjningsvinklarna 𝜷𝟏 och 𝜷𝟐 för skruvförband A Figur 11 visar att böjningsvinkelns lutning blir högre när förspänningen minskar. Det är inget överraskande resultat eftersom lägre förspänning gör att flänsarna lättare separerar och böjs därmed lättare. Det är däremot svårt att avgöra om lutningen är konstant eller om den varierar. Böjstyvheten går att få fram genom att ta inversen på kurvornas lutning i Figur 11 enligt Δ𝑀 Δ𝑀 𝑘1 = ,𝑘 = (4.1) Δ𝛽1 2 Δ𝛽2 Figur 12 illustrerar styvheten som beräknades enligt (4.1). Figur 12. Böjstyvhetsförändring för skruvförband A Figur 12 visar att böjstyvheten är konstant för små böjande moment. När belastningarna ökar kommer flänsytorna att börja separera ifrån varandra vilket gör att böjstyvheten minskar och konvergerar mot kurvan där skruvarna inte förspänns, en 15 % minskning av böjstyvheten. Med extrema moment börjar böjstyvheten att öka igen, detta beror på att de inre ytorna på flänsarna har deformerats så mycket att de börjar ta i varandra. Notera att styvheten inte är lika vid ändarna, detta beror på att skruvförbandet inte är symmetrisk utformad, se Figur 9. 24 Figur 13. Böjstyvhetsförändring för skruvförband A, förstorad i x-led I Figur 13 har böjstyvheten förstorats för att illustrera att båda förspänningarna ger en nästintill konstant böjstyvhet. Förspänningen som rekommenderas i ritningen är 117 𝑀𝑃𝑎 vilket är en relativt låg förspänningskraft. En ofta rekommenderad förspänning som till exempel Schneider [5] föreslår är att förspänna skruven till 80 % av dess sträckgräns. För en skruv med hållfasthetsklass 8.8 är sträckgränsen 640 𝑀𝑃𝑎 [17] och för en sådan skruv blir förspänningskraften 512 𝑀𝑃𝑎. Figur 13 visar att böjstyvheten kan hållas konstant för betydligt högre böjningsvinklar om skruvarna förspänns till 512 𝑀𝑃𝑎, vilket är ca 4 gånger större än vad ritningen rekommenderar. Att höja förspänningskraften skulle vara ett alternativ för att få mindre variationer i styvheten. Men man måste ta i åtanke att skruvarna är designad för att också utsättas för skjuvspänningar, det är därför inte säkert att skruven tål lika mycket som vanliga skruvar. 25 Kontaktkrafter Figur 14 illustrerar hur kontaktkrafterna i skruvförband A med skruvarna förspända till 117 𝑀𝑃𝑎 förändrades med ökat böjmoment enligt Tabell 2. Figur 14. Kontaktkraftsförändring [N] med högre belastning för en fläns förspänd till 𝟏𝟏𝟕 𝑴𝑷𝒂 Tabell 2. Momentbelastningsfall för Figur 14 [Nm] 𝟎 𝟓𝟎 ∙ 𝟏𝟎𝟑 𝟐𝟓𝟎 ∙ 𝟏𝟎𝟑 6 𝟔 2 ∙ 10 2,5 ∙ 106 𝟏 ∙ 𝟏𝟎 6 𝟔 25 ∙ 10 500 ∙ 106 𝟓 ∙ 𝟏𝟎 Figur 14 visar att kontaktkrafterna börjar förändras vid 1 𝑀𝑁𝑚, vid lägre böjande moment kontaktförändringen relativt liten. Med högre böjande moment fokuseras kontaktkrafterna på utsticken som används för att rikta in skruvförbanden vid montering. Utformningen på kontaktkrafterna förändras mindre med moment mellan 5 − 25 𝑀𝑁𝑚, vilket också bekräftas i Figur 12 där styvheten har konvergerat vid dessa moment. Som kuriosa börjar de inre ytorna att kommat i kontakt med varandra med extrema böjande moment som 500 𝑀𝑁𝑚 vilket gör att böjstyvheten skulle öka om skruvförbandet varit elastiskt i detta område. Dock har materialen plasticerat och gått av innan detta sker. 26 Förenklade modeller Kontaktanalyser tar betydligt längre tid än vanliga statiska FEM-analyser därför var det intressant att se vad förenklade modellers styvhet blir jämfört med kontaktanalyserna som gjordes ovan. Solid fläns Att skapa en solid modell där avrundningar och skruvhål har prioriterats bort för att få modellen så enkel som möjligt bör ge en styvhet som liknar kontaktanalysen vid små belastningar. Figur 15 illustrerar modellen. Figur 15. Förenklad solid modell där skruvhålen och avrundningar har tagits bort för skruvförband A Simpel fläns Ett sätt att minska styvheten är att helt ta bort flänsarna så att modellen liknar en cylinder istället. Figur 16. Förenklad solid modell där flänsarna har tagits bort så att den liknar en simpel cylinder för skruvförband A 𝟑𝟎° avfasning I denna modell har flänsarna bytts ut mot 30° avfasningar, detta används ofta för att beskriva geometriövergångar eftersom spänningszonens utformning påminner om en dubbelkon [1]. Därför är det intressant att se hur denna beter sig jämfört med kontaktanalysen. Modellen kan ses i Figur 17. Figur 17. Förenklad solid modell där flänsarna har fasats av med 𝟑𝟎° lutning 27 𝟒𝟓° avfasning Figur 18 visar en modell där flänsarna har fasats av med 45°. Denna modell kommer att få en lägre styvhet än modellen med 30° avfasade flänsar eftersom mer material har tagits bort från modellen. Figur 18. Förenklad solid modell där flänsarna har fasats av med 𝟒𝟓° lutning för skruvförband A Minskning av E-modul i flänsarna Om man delar upp skruvförbandet så att den består av fyra delar, två stycken axlar och två flänsar. Eftersom det är flänsseparationen som gör att styvheten minskar borde det gå att göra en förenklad modell genom att limma ihop flänsarna till en solid och sedan sänka dess E-Modul. Figur 19 illustrerar modellen. Figur 19. Förenklad solid modell där E-modulen i flänsarna har minskats med en faktor 𝒄 för skruvförband A Solid med mellanrum Figur 20 visar en helt solid modell utan skruvhål men har ett mellanrum i skarven, om mellanrummet är stort som det är i skruvförband A måste detta mellanrum tas i åtanke. Figur 20. Förenklad solid med ett mellanrum vid skarven för skruvförband A 28 Finita element analyser gjordes på dessa förenklade modeller, resultaten kan ses i Figur 21. Figur 21. Förenklade modeller av böjstyvheten för skruvförband A Figur 21 visar att den solida versionen är den styvaste modellen. Modellerna med avfasade flänsar har en aning lägre styvhet. Kontaktanalysen har ungefär 20 % lägre styvhet jämfört med den helt solida modellen. Skillnaden beror framförallt på att skruvförbandet har ett mellanrum där flänsarna inte är i kontakt med varandra vilket gör att böjstyvheten blir lägre jämfört med en helt solid version. Om mellanrummet tas i åtanke minskar böjstyvheten till ett värde som är aningen lägre än kontaktanalysen. Modellen med mellanrum har fortfarande en lägre böjstyvhet vilket beror framförallt på att avrundningarna som finns på originalet inte tas med i denna förenklade modell. Kontaktanalysen har minskat till 86 % av ursprungligt värde i Figur 21, genom att minska flänsarnas Emodul i den solida modellen med mellanrum till 86 % ger det ett liknande svar som det konvergerade värdet. Det skulle gå att få ett ännu noggrannare resultat genom att ta med avrundningarna som skruvförbandet har. Dessa förenklade modeller tog under 1 minut att lösa, jämför det med kontaktanalysen som tar 4-6 timmar att lösa, man kan få bort en stor flaskhals om man kan få en bra förenklad modell av böjstyvhetens förändring. 29 Analytisk modell för skruvförband A I Figur 22 jämförs den analytiska modellen som togs fram i kapitel 3.4.2 med FEM-modellen där skruvarna var förspända till 283 𝑀𝑃𝑎. Ekvation (3.15) och (3.16) har lösts med dimensioner från översiktsritningarna i Appendix A.1. Ett elementarfall för en fritt upplagd cylinder visas även här och dess böjning kan beskrivas med ekvation (3.7). Denna modell liknar FEM-modellen Simpel cylinder men den avancerade geometrin på insidan har inte tagits i åtanke, det är bara ett medelvärde mellan den största och minsta innerdiametern. Figur 22. Jämförelse mellan de analytiska modellerna och FE modellen för skruvförband A Böjstyvheten för de olika analytiska modellerna kan ses i Figur 23. Figur 23. Böjstyvheten för de analytiska modellerna för skruvförband A Figur 23 visar att den analytiska modellen har en liknande styvhet som FEM-modellen där flänsarna har fasats av med 30°. Men den analytiska modellen är för styv för att vara konkurrenskraftig med elementarfallet för en fritt upplagd cylinder. Den stora differensen beror på att skruvförbandet har ett mellanrum vid flänsytorna som gör att bara en del av flänsen är i kontakt med den andra flänsen. Om mellanrummen i flänsarna tas bort hade böjstyvheten ökat till liknande nivåer som den helt solida modellen och då hade den analytiska modellen varit mer relevant. Det finns dock flera nackdelar med det: Böjstyvheten minskar omedelbart istället för att vara konstant under en period eftersom kontaktytorna som är långt ifrån skruvarna inte påverkas nämnvärt av förspänningskraften och separerar därmed lättare. Vid tillverkning ökar kostnaderna betydligt om hela flänsens yta ska få en bra tolerans. 30 4.1.2 Vridstyvhet för skruvförband A Turbinerna i vattenkraftverken roterar på grund av att det forsande vattnet. Därmed uppstår även vridande moment i skruvförbanden. Att studera hur vridstyvheten i skruvförband A ändras med ökande vridmoment är därför också intressant. I det här fallet har friktionen mellan flänsytorna betydligt större roll än i böjfallen som har studerats i rapporten, friktionskoefficienten som använts har värdet 0.1. Figur 24 visar hur analysen har gjorts, den vänstra flänsens frihetsgrader låstes fullständigt, vilket innebär att 𝛽1 = 0° i detta belastningsfall. Skruvarna förspändes och ett ökande vridmoment applicerades vid den högra flänsens ände. Med hjälp av NX 10 löstes 𝛽2 ut. Figur 24. Vridningsvinkeln, 𝜷𝟐 , som uppstår ifrån det applicerade vridmomentet för skruvförband A Styvhetsförändringen beräknades med hjälp av 𝑘𝑣𝑟𝑖𝑑 = Δ𝑀 Δ𝛽2 Resultaten illustreras i Figur 25. Figur 25. Vridstyvhetsförändring för skruvförband A 31 (4.2) Vridstyvheten sjunker och konvergerar till 88 % av ursprungligt värde. Men vridstyvheten varierar på ett intressant sätt. Förspänningskrafterna klämmer åt flänsarna så att styvheten liknar en solid modell till en början. När momentet ökar börjar flänsytorna att glida mer och mer relativt varandra, det är den abrupta nedgången i kurvorna. Därefter tar skruvarna upp vridningen vilket visas i figuren med att styvheten ökar och konvergerar mot kurvan utan förspänning. Förenklade modeller Förenklade modeller har även gjorts för denna typ av belastningsfall. Samma typ av förenklade modeller har använts som när böjstyvheten undersöktes. Figur 26. Förenklade modeller av vridstyvheten för skruvförband A Figur 26 visar vridstyvheten för de förenklade modellerna. Från kontaktanalysen kan man se att vridstyvheten minskar till 88 % av ursprungligt värde. När E-modulen i flänsarna sänks till 88 % av ursprungligt av värde för den solida modellen med mellanrum blir vridstyvheten för hög jämfört med de konvergerade värdena, skillnaden beror på att skruvhålen inte finns med i de förenklade modellerna vilket gör att den får en högre vridstyvhet. Den förenklade modell som fungerar bäst för små vridmoment är den 45° avfasade modellen. Den är dock lite för styv, det beror på att skruvförbandet har fräst ut ytorna på flänsen vilket gör att styvheten minskar. 32 4.1.3 Axiell styvhet i skruvförband A Skruvförband A sträcks även ut i axiell riktning under drift, till exempel på grund av egenvikten, det strömmande vattnet och vattentrycket. Därför är det också intressant att undersöka hur den axiella styvheten förändras när skruvförbandet utsätts för axiella krafter. I Figur 27 har den vänstra änden fixerats, därefter belastas den högra änden med en axiell kraft. Figur 27. Skruvförband A utsätts för en axiell kraft Den axiella styvheten tas fram med följande samband 𝑘𝑎𝑥𝑖𝑒𝑙𝑙 = Δ𝐹 Δ𝑥 (4.3) I Figur 28 visas hur den axiella styvheten minskar när skruvförbandet deformeras. Figur 28. Axiell styvhetsförändring för skruvförband A Figur 28 visar att styvheten är konstant till en början, notera att den högra flänsen har i detta fall förflyttats i motsatt riktning vilket beror på att förspänningskraften deformerar flänsen i den riktningen. Den axiella styvheten har en tydligt konstant utformning i detta fall. Det beror på att när en axiell kraft appliceras deformeras hela skruvförbandet på samma sätt, flänsarna separerar därmed samtidigt. Den axiella styvheten minskar i detta fall med ungefär 52 % och då har flänsarna helt separerat ifrån varandra. Det krävs dock mycket kraft för att få flänsarna att separera ifrån varandra. Figur 28 visar det inte men 33 när skruvförbandet har deformerats 3 𝑚𝑚 har även 512 𝑀𝑃𝑎-kurvan konvergerat mot 52 % av sitt ursprungliga värde. Förenklade modeller I Figur 29 visas styvheten för de olika förenklade modellerna och hur de jämför med kontaktanalysen. Figur 29. Förenklade modeller av den axiella styvheten för skruvförband A Figur 29 visar att det inte går att få rimliga axiella styvheter genom att använda en simpel cylinder, 30° och 45° avfasade flänsar som förenklade modeller. Att göra en solid modell med ett mellanrum ger liknande vridstyvhet som kontaktanalysens högsta vridstyvhet vilket visar hur viktigt det är att ta med mellanrummen i de förenklade modellerna. Eftersom den solida modellen med mellanrum stämmer bra överrens med kontaktanalysen gör det möjligt att sänka dess flänsar med 52 % och få liknande svar som det konvergerade värdet. Mellanrummet har diametern är 810 𝑚𝑚 vilket motsvarar ungefär 70 % av flänsens diameter och uppstår som en bieffekt av att den vänstra flänsen har ett utstick som används för att rikta in den högra flänsen vid montering, utsticket har ett 4 𝑚𝑚 mellanrum mellan flänsarna för att garantera att de inte kommer i kontakt när de skruvas fast. 4 𝑚𝑚 är däremot ett stort mellanrum och går att minska betydligt utan att försvåra tillverkning eller montering. 34 4.2 Skruvförband B Detta skruvförband är en större version av skruvförband A med liknande typ av skruvar med en upphöjd yta vid skruvens mitt. Skruvarna har inget skruvhuvud, istället hålls den ihop av två muttrar. Skruvarna förspänns till 2688 𝑘𝑁 vilket motsvarar 283 𝑀𝑃𝑎 för dessa skruvar. Figur 30 visar hur modellen ser ut och översiktsritningar finnes i Appendix 0. En kontaktanalys gjordes även på denna med samma metoder som beskrivs i kapitel 4.1.1, förspänningen sattes vid muttern och skruvförbandet utsattes för ökande böjmoment vid ändarna av axlarna som kan ses i Figur 31. Figur 30. CAD-modell av skruvförband B Figur 31. Beräkningsnät för skruvförband B Båda flänsarna har liknande dimensioner men den största skillnaden kan ses i Figur 32 och Figur 33. På generatoraxeln har innerytorna och ytorna längs skruvcirkeln frästs ned. Detta innebär att det är de två upphöjda ytorna mellan skruvhålen som endast kommer i kontakt med den andra flänsen när skruvarna förspänns. Figur 32. Turbinfläns för skruvförband B Figur 33. Generatorfläns för skruvförband B 35 4.2.1 Böjstyvhet för skruvförband B Axlarna som skruvförband B håller fast kan utsättas för böjande moment under drift. Figur 34 visar vad böjningsvinkeln blev vid ändarna när axlarna utsattes för böjande moment. Figur 34. Förändring av böjningsvinklarna 𝜷𝟏 och 𝜷𝟐 för skruvförband B Med ekvation (4.1) fås böjstyvheten fram och resultaten visas i Figur 35. Figur 35. Böjstyvhetsförändring för skruvförband B Figur 35 visar att böjstyvheten är relativt konstant upp till 0,005° − 0,01°. Därefter börjar böjstyvheten minska. Böjstyvheten för de förspända skruvarna hinner inte konvergera mot kurvan med 0 𝑀𝑃𝑎 förspänning innan de inre ytorna kommer i kontakt med varandra. Detta beror på att det bara är 0,3 𝑚𝑚 mellanrum i detta skruvförband. De inre ytorna påverkas inte nämnvärt av förspänningen vilket gör att det böjande moment som krävs för att de ska börja komma i kontakt med varandra inte förändras när förspänningen ökar. Böjstyvheten har näst intill samma värde vid 𝛽1 och 𝛽2 , anledningen är att skruvförband B har en mer symmetrisk design jämfört med skruvförband A. Figur 35 visar att böjstyvheten inte minskar lika mycket om skruvarna förspänns till 512 𝑀𝑃𝑎, det är dock inte lika stor skillnad som det var i skruvförband A. Det beror på att förspänningen är dubbelt så stor i skruvförband B jämfört med skruvförband A där förspänningen är fyra gånger större. 36 Kontaktkrafter Figur 36 visar hur kontakkrafterna förändras för ett skruvförband som är förspänd till 283 𝑀𝑃𝑎 när det böjande momentet ökar enligt Tabell 3. Figur 36. Kontaktkrafternas förändring [N] för skruvförband B med 𝟐𝟖𝟑 𝑴𝑷𝒂 förspända skruvar Tabell 3. Momentbelastning för Figur 36 [Nm] 𝟎 𝟏𝟐𝟓 ∙ 𝟏𝟎𝟑 𝟔𝟐𝟓 ∙ 𝟏𝟎𝟑 6,25 ∙ 106 8,75 ∙ 106 𝟐, 𝟓 ∙ 𝟏𝟎𝟔 6 𝟔 62,5 ∙ 10 125 ∙ 106 𝟏𝟐, 𝟓 ∙ 𝟏𝟎 Kontaktkrafterna har en liknande utformning som skruvförband A i Figur 14, dock tar inte ytorna närmast skruvcirkeln i varandra på grund av att generatoraxelns yta har frästs ned runt skruvcirkeln, därför ser man inga kontaktrafter längs skruvcirkeln. Upp till 2,5 𝑀𝑁𝑚 förändras ytorna ytterst lite vilket gör att böjstyheten är nästan konstant. Med högre böjande moment börjar den inre kontaktytan att separera ifrån varandra vid 6 𝑀𝑁𝑚. Vid böjmoment högre än 60 𝑀𝑁𝑚 kommer de inre ytorna i kontakt med varandra som gör att böjstyvheten börjar ökar. 37 Förenklade modeller Förenklade modeller gjordes också för skruvförband B. Eftersom detta skruvförband liknar skruvförband A kunde samma typer av förenklingar göras. En solid fläns gjordes även för skruvförband B som kan ses i Figur 37. Varken Solid fläns avrundningarna eller skruvhålen som finns i skruvförbandet har inte tagits med i denna för att gör modellen så enkel som möjligt. Figur 37. Förenklad solid modell där skruvhålen och avrundningarna har tagits bort för skruvförband B Simpel fläns Figur 38 visar en modell där flänsarna har tagits bort som gör att modellen liknar en cylinder. På detta sätt kan böjstyvheten minskas betydligt. Figur 38. Förenklad simpel fläns för skruvförband B 𝟑𝟎° avfasning Figur 39 illustrerar en 30° avfasad modell av skruvförband B. Den avfasade delen har en aning högre densitet för att motsvara materialet som tagits bort. Figur 39. Förenklad solid modell där flänsarna har fasats av med 𝟑𝟎° lutning 38 𝟒𝟓° avfasning Figur 40 visar en modell med flänsarna avfasade 45°. Den avfasade delen har en aning högre densitet för att motsvara materialet som tagits bort. Figur 40. Förenklad solid modell där flänsarna har fasats av med 𝟒𝟓° lutning Minskning av E-modul i flänsarna En solid modell som har E-modulen sänkt i flänsarna gjordes även för skruvförband B. Detta gör att modellen blir betydligt mjukare i flänsarna enbart. Den har sänkts med en faktor 𝑐 eftersom skruvförbandets styvhete minskade och konvergerade mot ett värde. Därmed borde denna modell ge samma böjstyvhet som det konvergerade värdet. Figur 41. Förenklad solid modell där E-modulen i flänsarna har minskats med en faktor 𝒄 för skruvförband B Solid med tre mellanrum Figur 42 visar en helt solid modell utan skruvhål men har tre mellanrum i skarven för att likna utformningen som skruvförband B har. Figur 42. Förenklad solid med ett mellanrum vid skarven för skruvförband B Variationer på denna finns med avrundade kanter och lika breda kontaktytor för att simulera olika styvheter vid olika belastningar. 39 I Figur 43 har de förenklade modellerna lösts med hjälp av NX 10, därefter har ekvation (4.1) använts för att ta fram böjstyvheten. Figur 43. Förenklade modeller av böjstyvheten för skruvförband B Figur 43 visar att den styvaste modellen är den helt solida modellen. Därefter kommer modellerna med avfasade flänsar. I detta skruvförband har kontaktanalysen en högre böjstyvhet än den simpla cylindern. Böjstyvheten är konstant upp till 0,01° och minskar därefter tills de inre ytorna kommer i kontakt med varandra och börjar därför öka igen. För detta skruvförband krävs det att både utfräsningarna och avrundade kanter mellan flänsen och axeln för att ge ett bra jämförbart resultat med kontaktanalysen. Att sänka flänsarnas E-modul för den till 86 % av dess ursprungliga värde stämmer väl överrens med de konvergerade värdena i kontaktanalysen. Det inte finns dock ingen likhet mellan böjstyvheterna för skruvförband A och B. I skruvförband B är böjstyvheten för kontaktanalysen betydligt högre relativt de förenklade modellerna, 86 % av den helt solida modellens böjstyvhet. I skruvförband A är böjstyvheten 80 % av dess helt solida modell. Det som skiljer dessa skruvförband åt är framförallt hur mellanrummen mellan flänsarna är designade. Skruvförband A har ett mellanrum på 4 𝑚𝑚 som är 810 𝑚𝑚 i diameter på grund av utsticket som används för att rikta in axeln vid montering. Skruvförband B har tre mellanrum med olika gap, utsticket har diametern 540 𝑚𝑚 och där är mellanrummet mellan flänsarna 2 𝑚𝑚. Mellan diametrarna 540 − 900 𝑚𝑚 finns ett mellanrum som är 0,3 𝑚𝑚, eftersom mellanrummet är så litet blir dess inverkan på böjstyvheten liten. Det finns även ett mellanrum vid skruvcirkeln som också är 0,3 𝑚𝑚. Detta innebär att för att få noggranna förenklade modeller måste de avancerade geometrierna också tas i åtanke. 40 Analytisk modell för skruvförband B I Figur 44 visas en jämförelse med den analytiska modellen som togs fram i kapitel 3.4.2 och FEM-modellen för skruvförband B med skruvarna förspända till 283 𝑀𝑃𝑎. Ekvation (3.15) och (3.16) har lösts med dimensioner från översiktsritningarna i Appendix 0. Ett elementarfall för en fritt upplagd cylinder visas även här där ekvation (3.7) löstes för detta skruvförband. Figur 44. Jämförelse mellan de analytiska modellerna och FE modellen för skruvförband B Styvheten för de olika analytiska modellerna kan ses i Figur 45 Figur 45. Böjstyvheten för de analytiska modellerna för skruvförband B Figur 45 visar att den analytiska modellen är jämförbar med de avfasade FEM-modellerna men är fortfarande för styv för att kunna användas för att noggrant approximera böjstyvheten i detta skruvförband. Notera att den analytiska modellen är styvare än modellen med 30° avfasade flänsar, för skruvförband A är det tvärtom. Anledningen är att skruvförband A har en avancerad geometri på insidan, den analytiska modellen tar bara ett medelvärde på innerdiametern och i (3.8) upphöjs innerdiametern med fyra. Små variationer i innerdiametern ger därmed betydligt större variationer i böjstyvheten. I detta fall är elementarfallet ett alternativ eftersom den kommer att fortsätta vara relevant även när böjstyvheten minskar. 41 4.2.2 Vridstyvhet för skruvförband B Skruvförband B utsätts för vridande moment när turbinen roterar hur vridstyvheten förändras är därför också intressant att studera. Figur 46 visar hur skruvförbandet belastades, vänstra flänsen fixerades vid axelns ände och ett vridande moment applicerades på den högra axelns ände. Figur 46. Vridningsvinkeln, 𝜷𝟐 , som uppstår ifrån det applicerade vridmomentet för skruvförband B När NX har löst detta belastningsfall används ekvation (4.2) för att ta fram vridstyvheten och resultaten illustreras i Figur 47. Figur 47. Vridstyvhetsförändring för skruvförband B Förändringen av vridstyvheten i Figur 25 påminner om skruvförband As förändring i Figur 25. Vridstyvheten är hög till en början men minskar med högre vridande moment på grund av att flänsytorna börjar glida relativt varandra. Vid en viss punkt börjar styvheten att öka igen, detta sker när skruvarna börjar klämmas fast flänsarna. Vilket gör att flänsytorna slutar glida relativt varandra och vridstyvheten börjar konvergera mot 89 % av dess ursprungliga värde. 42 Förenklade modeller Figur 48 visar en jämförelse mellan kontaktanalysen och de förenklade modellerna som presenterades i kapitel 4.2.1. Figur 48. Förenklade modeller av vridstyvheten för skruvförband B Figur 48 visar att en solid modell med flera mellanrum ger liknande vridstyvhet som kontaktanalysen har vid små vridmoment. Men vridstyvheten minskar fort när vridmomentet ökar vilket gör att det bara gäller upp till vridningsvinkeln 0,05°. Kontaktanalysen minskar till 89 % av ursprungligt värde, när detta värde används för att sänka flänsarna i den solida modellen med flera mellanrum ger det ett bra approximerat värde på den konvergerade vridstyvheten. 43 4.2.3 Axiell styvhet i skruvförband B Skruvförband B kan utsättas för stora krafter längs dess axel, bland annat från vattnet som forsar förbi turbinen. Den axiella styvheten har undersökts genom att fixera den högra axelns ände och applicera en axiell kraft på den högra axelns ände och därefter mäta deformationen i x-led som i Figur 49. Figur 49. Skruvförband B utsätts för en axiell kraft Den axiella styvheten kan beskådas i Figur 50 som togs fram med ekvation (4.3). Figur 50. Axiella styvhetsförändringen i skruvförband B Figur 50 visar att detta skruvförband har två platåer som den axiella styvheten är konstant på. Detta beror på att detta skruvförband har fräst bort ytorna runt skruvcirkeln, se Figur 33. Vid första platån är det två ytor som är i kontakt med den andra flänsen. Den innersta ytan separerar först vilket illustreras i Figur 50 som den första styvhetsfallet. Den andra platån innebär att den inre ytan har separerat tills den inre kontaktytan har ungefär samma bredd som den yttre kontaktytan, den axiella styvheten har då sjunkit till 76 % av ursprungligt värde. Den axiella styvheten förändras inte för att ytorna som fortfaranade är i kontakt med varandra känner av en stor 44 del av förspänningen. Tillslut separerar den inre kontaktytan helt och då minskar den axiella styvheten ytterligare och konvergerar tillslut till 25 % av dess ursprungliga värde vid 2 𝑚𝑚. Förenklade modeller I Figur 51 visas de förenklade modeller som använts för att beskriva den axiella styvheten för skruvförband B. Figur 51. Förenklade modeller av den axiella styvheten för skruvförband B Figur 51 visar att den helt solida modellen, de avfasade modellerna och den cylindermodellen är för styva för att kunna ge bra approximation på den axiella styvheten. Men när man tog med de tre mellanrummen som finns i skruvförband B är det möjligt att få väldigt jämförbara värden som med kontaktanalysen. Det är därmed viktigt att ta med mellanrummen som finns i skruvförbanden för att kunna få bra resultat, speciellt för axiella styvheter. Eftersom modellen med tre mellanrum stämmer bra överens med kontaktanalysens största axiella styvheter gör det möjligt att kunna minska E-modulen i flänsarna för att kunna replikera styvhetsminskningen. Vid andra platån har kontaktanalysen sjunkit till 76 % av sitt ursprungliga värde. Därefter minskar styvheten snabbt till 31 % av dess ursprungliga värde och konvergerar långsamt till kurvan med inga förspända skruvar som är ungefär 25 % av kontaktanalysens ursprungliga värde. Som Figur 51 visar korrelerar E-modulsminkningarna väldigt bra med kontaktanalysen. Det går att få fram styvhetsminskningarna med förenklade modeller genom att ändra hur stor kontaktytorna är. För den andra platån är den inre kontakytan lika bred som den yttre kontakytan. Genom att minska den inre kontaktytan för soliden med tre mellanrum så att den blir lika bred som den yttre kontaktytan går det att få ett jämförbart värde som kontaktanalysen. För den tredje platån tas den inre kontaktytan bort helt, denna har dock inte samma noggrannhet vilket beror framförallt på att skruvarna inte har modellerats i denna modell därför blir den axiella styvheten lägre. 45 4.3 Skruvförband C Skruvförband C har en stor övre röraxel som kopplas till en turbinaxel. Skruvarna är vanliga utvändigt gängade skruvar i storlek M90. För att ta upp vridmomenten som uppstår när turbinen roterar används en hylsa som presspassas in, denna tar upp vridmomenten istället för skruven. Detta gör att man kan förspänna skruvarna mer eftersom de inte ska ta upp skjuvande spänningar. Dessa skruvar förspänns enligt ritningarna till 4200 𝑘𝑝/ 𝑐𝑚2 vilket motsvarar ca 412 𝑀𝑃𝑎, översiktsritningar för modellen kan studeras i Appendix 0. Figur 52 illustrerar hur det modellerade skruvförbandet ser ut. Figur 52. CAD-modell av skruvförband C Figur 53 visar beräkningsnätet för modellen och var förspänningskraften har placerats, den lasten har definierats där skruvhålsgängorna startar. Kontaktytorna har definierats mellan flänsarna, hylsan och flänsarna skruven och hylsan. Skruvarna limmades fast vid skruvhuvudet, hylsan limmades fast i turbinaxeln så att den sitter kvar i skruvhålet. Figur 53. Beräkningsnät för skruvförband C 46 4.3.1 Böjstyvhet för skruvförband C Axlarna som skruvförbandet sitter på kan böjas på grund av bland annat vattnets rörelse, när skruvförbandet böjs kommer flänsarna att börja separera vilket gör att böjstyvheten i skruvförbandet minskar. Figur 54 visar hur flänsrotationen ändrades med högre böjande moment. Skruvarna förspändes till 412 𝑀𝑃𝑎, 512 𝑀𝑃𝑎, 720 𝑀𝑃𝑎 och 0 𝑀𝑃𝑎. De högre förspänningskrafterna motsvarar 80 % av sträckgränsen för en 8.8 och 10.9 skruv. Figur 54. Förändring av böjningsvinklarna 𝜷𝟏 och 𝜷𝟐 för skruvförband C Med (4.1) beräknas böjstyvheten och visas Figur 55. Figur 55. Böjstyvhetsförändring i skruvförband C I Figur 55 illustreras hur böjstyvheten minskar för skruvförband C, böjstyvheten minskar till ca 60 − 65 % av sitt ursprungliga värde. Detta beror framförallt på röraxelns utformning, den har betydligt större diameter jämfört med turbinaxeln och relativt liten cylindertjocklek. Att förspänna skruvarna med högre förspänningar påverkade inte böjstyvhetens förändring lika mycket som i skruvförband A vilket beror att förspänningen enligt ritning redan är hög. 512 𝑀𝑃𝑎 och 720 𝑀𝑃𝑎 motsvarar en 20 % respektive 70 % ökning av förspänningen jämfört med den rekommenderade förspänningen 412 𝑀𝑃𝑎, vilket förklarar varför styvheten inte påverkas lika mycket som för till exempel skruvförband A där det var en fyrfaldig ökning. Kontaktkrafter Figur 56 visar hur kontaktkrafterna mellan flänsarna ändras när momentet höjs för skruvförband C med 412 𝑀𝑃𝑎 förspända skruvar. De här visar också att när momentet stiger över 1 − 4 𝑀𝑁𝑚 börjar tydliga kraftförändringar att uppträda. 47 Figur 56. Kontaktkraftsförändring med högre belastning med skruvar förspända till 𝟒𝟏𝟐 𝑴𝑷𝒂 Tabell 4. Momentbelastningsfall för Figur 56 [Nm] 𝟖𝟎𝟎 ∙ 𝟏𝟎𝟑 16 ∙ 106 56 ∙ 106 400 ∙ 106 𝟎 𝟖 ∙ 𝟏𝟎𝟔 𝟑𝟐 ∙ 𝟏𝟎𝟔 𝟖𝟎 ∙ 𝟏𝟎𝟔 48 Figur 56 visar en annorlunda utformning av kontaktkrafterna. Det är en liten förändring upp till 8 𝑀𝑁𝑚, därefter börjar kontaktrafterna på den övre halvan att minska. Utformningen slutar att förändras runt 30 – 40 𝑀𝑁𝑚. Notera i detta skruvförband att det är en liten yta som är kontakt vid 30 − 40 𝑀𝑁𝑚 jämfört med skruvförband A och B. Detta beror på att den övre röraxeln är betydligt större än turbinaxeln, detta gör det möjligt för den övre röraxeln att deformeras över turbinaxeln och därmed blir de största kontaktkrafterna vid turbinaxelns kanter, se Figur 57 Figur 57. Deformation av skruvförband C vid högre moment Figur 57 visar överdrivet deformationen av skruvförband C, jämför med skruvförband A i Figur 10 för att se hur de olika typerna av skruvförband deformerar. 49 Förenklade modeller Skruvförband C har en annorlunda utformning jämfört med skruvförband A och B. Detta gjorde det svårt att hitta förenklade modeller. Att göra en solid av den är fullt möjligt men dessa skulle fortfarande vara för styva som resultaten från skruvförband A och B visade. Att göra avfasningar för denna modell fungerar inte på grund av att ytterdiametrarna är så olika varandra. Att sänka E-modulen på modellen är fortfarande ett alternativ. Eftersom styvheten sjönk till ca 60 % − 65 % av ursprungligt värde för skruvförband C beroende på vilken ände som studeras är det ett bra värde att testa med. Denna modell har inte heller tydliga flänsar som man kan sänka E-modulen på men om man definierar soliderna som skruvarna sitter på som flänsar och sänker Emodulen går det att använda samma metod som i skruvförband A och B. I Figur 58 har väggen som skruvarna sitter på och delen av turbinaxeln som skruven är gängad i sänkts till 60 % och 65 % av ursprungligt värde, 124 𝑀𝑃𝑎 respektive 135 𝑀𝑃𝑎. Figur 58. E-modulsförändring i skruvförband C Figur 58 visar att hålen har tagits med i både den solida versionen och när E-modulen sänks. Behovet att ta med hålen har inte behövts i de förenklade modeller som gjordes med skruvförband C. Men med detta skruvförband behövdes hålen vara med för att få acceptabla värden. Visserligen har detta skruvförband fler skruvhål än både skruvförband A och B, 18 st. skruvhål, volymen som tas bort från hålen verkar ha en signifikant påverkan på styvheten för denna modell. 50 Figur 59. Förenklade modeller av böjstyvheten för skruvförband C Figur 59 visar hur styvheten för de olika förenklade modellerna jämför med kontaktanalysen. Den solida modellens styvhet är lika stor som kontakanalysen i början. Detta är ett bra beteende eftersom kontaktanalysens styvhet bör vara lik en solids om den inte utsätts för större externa belastningar. Att sänka E-modulen fungerar inte lika bra som den har gjort i skruvförband A och B. Men det är en relativt enkel förenkling som ger ett ungefärligt svar på under 1 min jämfört med kontaktanalysen som tog 5 timmar att lösa. 51 4.3.2 Vridstyvhet för skruvförband C Skruvförband C utsätts för vridande moment när turbinen roterar detta gör det intressant att undersöka hur vridstyvheten förändras. Vridstyvheten undersöktes genom att fixera den vänstra röraxelns ände och ett vridande moment applicerades på den högra turbinaxelns ände och illustreras i Figur 60. Figur 60. Vridningsvinkeln, 𝜷𝟐 , som uppstår ifrån det applicerade vridmomentet för skruvförband C Resultaten från FEM-analyserna illustreras i Figur 61. Figur 61. Vridstyvhetsförändring för skruvförband C Figur 61 visar att styvheten minskar till 60 % av ursprungligt värde. Vilket är betydligt lägre än skruvförband A och B, detta beror framförallt på dess utformning. Det är dock ett styvare skruvförband, vridstyvheten är 10 − 20 gånger större än både skruvförband A och B. 52 Förenklade modeller I Figur 62 visas hur bra modellen med minskad E-modul, se Figur 58, stämmer med vridstyvheten för kontaktanalysen. Figur 62. Förenklade modeller av vridstyvheten för skruvförband C Figur 62 visar att den solida modellen med skruvhål ger ett liknande värde som kontaktanalysen men om dess E-modul minskas med 40 % blir vridstyvheten för styv jämfört med kontaktanalysens konvergerade värde. 53 4.3.3 Axiell styvhet i skruvförband C Skruvförband kan även utsättas för axiella krafter. Genom att fixera den vänstra flänsen och applicera en axiell kraft går det att se hur den axiella styvheten förändras med ökad kraft. Figur 63. Skruvförband C utsätts för en axiell kraft och deformeras längden x I Figur 64 visas resultaten hur den axiella styvheten minskas när skruvförband C utsätts för en axiell kraft. Figur 64. Axiell styvhetsförändring för skruvförband C Figur 64 visar att skruvförband C behöver deformeras mer än 6 mm för att konvergera mot samma axiella styvhet som skruvförbandet har ytan någon förspänningskraft. För axiella krafter sjunker den axiella styvheten till ungefär 45 % av det ursprungliga värdet. Figur 65 visar hur den axiella styvheten förändras vid små deformationer. 54 Figur 65. Axiell styvhetsförändring för skruvförband C, förstorad i x-led Förenklade modeller Figur 64 visade att den axiella styvheten sjönk till 45 % av ursprungligt värde. I Figur 66 har den solida modellen med skruvhål jämförts med kontaktanalysen, solidens E-modul har också sänkts till 45 % av ursprungligt värde för att undersöka om det ger det konvergerade värdet. Figur 66. Förenklade modeller av den axiella styvheten för skruvförband C Figur 66 visar att sänka E-modulen på den solida modellen med skruvhål till 45 % ger liknande värde som kontaktanalysens konvergerade värde. Dock är den solida modellen inte tillräckligt styv vid små deformationer, vilket nog beror mest på skruvhålen. 55 5 Slutsatser I detta avsnitt besvaras forskningsfrågorna som ställdes i 1.3 Syfte. Hur ändras styvheten i ett skruvförband när det utsätts för olika belastningsfall? För böjande belastningar är böjstyvheten konstant till en början, detta beror på att flänsytorna som är i kontakt med varandra bara är i närheten av skruvcirkeln. Förspänningskraften gör att ytorna som är närheten av skruvcirkeln pressas ihop med liknande krafter vilket innebär att ytorna inte kommer att separera förrän det externa momentet överstiger förspänningskraften i skruvarna. Förspänningskraften sprids dock inte ut jämnt på alla ytor som är i kontakt med varandra utan minskar ju längre ifrån kontaktytan är vilket gör att dessa ytor separerar för små belastningar. Detta innebär att det inte är nödvändigt att ha en så stor yta som möjligt i kontakt med varandra om det är viktigt att få en konstant styvhet, visserligen blir böjstyvheten högre tillfälligt men minskar omedelbart istället för att vara konstant för vissa belastningar. Fördelen med konstant styvhet är att man enklare kan göra vidare beräkningar för att till exempel se vad skruvförbandets egenfrekvenser blir. Böjstyvheten börjar tillslut att minska. Hur mycket den minskar har varierat för alla skruvförband, för skruvförband som liknar A och B har också flänsskarven utformning betydelse. Skruvförand A och B liknar varandra medan C är en helt annan typ av skruvförband. Som väntat var det tydlig skillnad på hur böjstyvheten såg ut mellan de olika typerna. Men trots att A och B har en liknande geometri är det skillnad på deras böjstyvheter. Om en helt solid modell jämförs med kontaktanalysen ser man att skruvförband A har en 20 % lägre böjstyvhet jämfört med en helt solid modell. För skruvförband B är det 15 % skillnad mellan en helt solid modell. Detta var ett överraskande resultat med tanke på att skruvförband B har en mindre area som är i kontakt med varandra. Anledningen är mellanrummen som finns i flänsarna, skruvförband A har ett större mellanrum och det är det som minskar dess böjstyvhet. När skruvförband utsätts för vridande moment har friktionen en större påverkan, för små vridande moment är styvheten hög men sjunker med högre moment på grund av att flänsytorna börjar glida mer och mer relativt varandra. Skruvförbandet kommer tillslut att ha vridits så mycket att skruvarna och hylsorna fixerar flänsarna, detta gör att styvheten ökar en aning och konvergerar mot ett värde eftersom flänsytorna inte glider relativt varandra längre. För axiella belastningar kan man se att den axiella styvheten är nästintill konstant till en början sedan minskar den axiella styvheten och konvergerar mot 52 % av ursprungligt värde för skruvförband A och för skruvförband B sjunker den ned till 22 % av sitt ursprungliga värde. Skruvförband C konvergerar mot 45 % av ursprunglig axiell styvhet. Alla skruvförband har en tydlig platå med axiella belastningar, detta beteende beror på att en axiell belastning deformerar hela skruvförbandet symmetriskt. I framförallt skruvförband A är förspänningskraften som rekommenderas väldigt låg. Att höja förspänningskraften till 80 % av skruvens sträckgräns gör att det tar längre tid för de olika styvheterna att konvergera. Det innebär också att styvhetsförändringen är liten vid små vinklar vilket förbättrar de förenklade modellernas noggrannhet. 56 Hur väl stämmer olika förenklade modeller överrens med varandra? Tiden det tog för att lösa kontaktanalyserna varierade mellan 4-6 timmar. De förenklade modeller som har visats i rapporten tar mindre än 1 minut att lösa vilket innebär att en stor flaskhals försvinner. På grund av skruvförbandens utformning är det inte möjligt att göra en enkel modell, såsom att fasa av flänsarna med 30° eller 45°, för de får ofta för hög styvhet. En stor del beror på att flänsarna ofta har ett mellanrum, både för att minska ytorna som behöver ha hög tolerans och för att göra det möjligt att enkelt rikta in axlarna vid montering. Detta innebär att man måste ta med mellanrummen som finns i skruvförbanden i de förenklade modellerna. Med mellanrummen kan man få jämförbara styvheter men detta gäller bara vid små belastningar då styvheten är näst intill konstant. De förenklade modeller som har presenterats i resultaten har konstant styvhet, vilket gör det svårt att få jämförbar styvhet vid högre belastningar eftersom kontaktanalysens styvhet minskar. Det är möjligt att få fram hur mycket styvheten sjunker genom att först skapa en modell som stämmer bra överrens med kontakanalysen vid små belastningar. Därefter förändra hur stor del av solidens flänsar som är i kontakt med varandra som i Figur 51, kontaktytorna bör ha liknande som kontaktanalysen vid de belastningarna. Utmaningen är att korrekt förutspå hur kontaktytorna har förändrats vilket kräver en del erfarenhet för att med säkerhet kunna använda resultaten utan att först ha gjort en kontaktanalys. Att sänka E-modulen på en solid modells flänsar med samma procentsats som de konvergerade värdena i kontaktanalyserna har visat sig ge väldigt liknande resultat. I skruvförband A och B stämmer sänkningen väldigt bra för böj-, vrid- och axiell styvhet, i skruvförband C är precisionen inte lika hög vilket beror framförallt på svårigheten i att avgöra vilka delar som E-modulen ska minskas. Det är dock viktigt att den solida modellen stämmer väl överrens med kontaktanalysen vid små belastningar. För att få den noggrannheten måste mellanrummen som finns i flänsarna tas i åtanke, även de stora avrundningarna är viktiga att ta med i de förenklade modellerna, speciellt vid böjande moment. Det som är intressant med den här metoden är att man kan simulera styvhetsminskningen i en solid modell genom att ändra E-modulen i flänsarna när belastningarna ökar. På detta sätt kan man kringgå den långdragna kontaktanalysen men ändå få jämförbara styvhetsminskningar. 57 6 Diskussion Undersökningen som har gjorts på tre stycken olika skruvförband har visat att för små belastningar är kontaktanalysens styvhet nästan konstant vilket gör det möjligt att använda förenklade modeller eftersom deras styvheter förblir konstant för alla belastningar. Det är dock viktigt att ta med de avancerade geometrierna i de förenklade modellerna, speciellt insidorna på flänsarna har haft stor påverkan för att få jämförbara resultat. Att göra kontaktanalyser visade sig vara svårare än planerat. Kunskapströskeln var relativt hög för att bara kunna lösa analysen. Att felsöka modellen kräver en hel del kunskap och erfarenhet, för felmeddelandena är ofta för generella i NX vilket gör det svårt att förstå felet. När kontaktanalyserna ska lösas tar varje kontaktanalys generellt 4-6 timmar att lösa vilket har varit en stor flaskhals i detta projekt. När analysen är klar ska värdena kopieras över till Matlab. Många värden måste kopieras över manuellt vilket gör att risken för fel ökar. Detta har inte varit ett stort problem för det syns tydligt i styvhetsgraferna om det är något enstaka fel, men det är enkelt att råka välja fel fil och skriva över den. Lagringen av värdena hade kunnat automatiseras genom att programmera in förfarandet men den kunskapen har inte funnits. Under litteraturstudien var det svårt att hitta relevant information, det finns litteratur som diskuterade flänsars beteende men det var ingen som gick in på hur externa moment påverkar styvheten. Detta har inte varit något hinder för finita element analyserna. Men har varit ett större bekymmer för de analytiska analyserna. Det finns artiklar som snuddar vid ämnet som [7] där tätade flänsar utsätts för externa moment. Men när artiklarna studeras i detalj finner man olika slarvfel eller så refererar de viktiga härledningar vidare som är svårfunna vilket har gjort det svårt att reproducera deras resultat och kunna applicera det på detta projekt. I [9] appliceras axiella krafter för att se hur mycket flänsarna separerar ifrån varandra och de visade att flänsseparationen är linjär. Men de undersökte inte hur böjande och vridande moment påverkas vilket också är relevanta belastningar som måste tas i åtanke. Experimentella tester hade varit intressant att göra för att bekräfta resultaten men det hade krävts tidigare erfarenheter med kontaktanalyser och testriggen hade redan behövt vara färdigställd för att hinna med det. Förslag till fortsatt arbete Utöver att göra experiment för att bekräfta styvhetsförändringarna ligger en av förbättringspunkter i att förbättra den analytiska modellen för böjning. Den stämmer bra om skruvförbanden inte har mellanrum, tyvärr är detta inte realistiskt att anta eftersom mellanrummen i flänsarna har visats sig att ha en stor påverkan på styvheten. Det är nog möjligt att få en bättre modell genom att anta att flänsen också är elastisk, därefter ta med mellanrummen i flänsarna. Att minska E-modulen har kunnat beskriva styvhetsminskningen på ett övertygande sätt. Detta fungerar speciellt bra för skruvförband A och B där det är enkelt att bestämma vilken del som E-modulen ska sänkas i. Detta förutsätter dock att man vet när ett skruvförband har konvergerat och hur mycket styvheten minskar. Det är möjligt att ta reda på hur mycket styvhetsminskningen är genom att förutspå hur kontaktytorna ser ut, därefter göra en förenklad modell som liknar den utformningen och jämföra styvheterna. Vilket leder till nästa förslag att arbeta vidare med; att ta reda på när skruvförbanden har konvergerat och finna en bra kurvanpassning för en sådan förenklad modell. 58 7 Referenser [1] E. Eriksson, E. Kassfeldt, E. Höglund, O. Isaksson, S. Lagerkrans och J. Lundberg, Maskinelement, Luleå: Luleå Tekniska Universitet, 1993. [2] A. Eames, ”Nord-Lock,” 04 Februari 2012. [Online]. Available: http://www.nord-lock.com/wpcontent/uploads/2012/10/Bolted_2_2012_EN.pdf. [Använd 18 April 2015]. [3] Swedish Fasteners Network, ”Swedish Fasteners Network,” 23 07 2010. [Online]. Available: http://handbok.sfnskruv.se/template.asp?lank=150. [Använd 18 04 2015]. [4] E. O. Waters, D. B. Rossheim, D. B. Wesstrom och F. S. G. Williams, Development of General Formulas for Bolted Flanges, Chicago: Taylor Forge and Pipe Works, 1949. [5] R. W. Schneider, ”Flat Face Flanges With Metal-to-Metal Contact Beyond the Bolt Circle,” ASME, J. Eng. Gas Turbines Power, vol. 90, nr 1, pp. 82-88, 1968. [6] A.-H. Bouzid och H. Galai, ”Analytical Modeling of Flat Face Flanges With Metal-to-Metal Contact Beyond the Bolt Circle,” ASME,Journal of Pressure Vessel Technology, vol. 132, nr 6, pp. 061207-061207-8, 2010. [7] A.-H. Bouzid, ”On the Effect of External Bending Loads in Bolted Flange Joints,” ASME, Journal of Pressure Vessel Technology, vol. 131, nr 2, pp. 021201-1 - 021201-8, 2009. [8] W. J. Koves, ”Analysis of Flange Joints Under External Loads,” ASME, Journal of Pressure Vessel Technology, vol. 118, pp. 59-63, 1996. [9] D. Joshi, P. Mahadevan, A. Marathe och A. Chatterjee, ”Unimportance of geometric nonlinearity in analysis of flanged joints with metal-to-metal contact,” Internatiional journal of pressure vessels and piping, vol. 84, nr 7, pp. 405-411, 2007. [10] B. Sundström, Handbok och formelsamling i Hållfasthetslära, Stockholm: Institutionen för hållfasthetslära, KTH, 2008. [11] R. Patel och B. Davidsson, Forskningsmetodikens grunder; Att planera, genomföra och rapportera en undersökning, Lund: Studentlitteratur, 2011. [12] L. Kvarnström, ”Förspänningsförluster i maskinfundamentens skruvförband,” 2008. [Online]. Available: http://pure.ltu.se/portal/sv/studentthesis/foerspanningsfoerluster-i-maskinfundamentens-skruvfoerband-en-jamfoerelsemellan-praktik-och-teori(a39a2bcf-b193-4e73-b2fe-7040a147fb15).html. [Använd 23 Januari 2015]. [13] Siemens, ”NX 10,” Siemens, 2015. [Online]. Available: http://www.plm.automation.siemens.com/en_us/products/nx/10/. [Använd 30 04 2015]. [14] H. Lundh, Grundläggande hållfasthetslära, Stockholm: Institutionen för hållfasthetslära, KTH, 200. [15] S. Sjöström, ”TABELLSAMLING I HÅLLFASTHETSLÄRA,” 15 3 2012. [Online]. Available: http://www.solid.iei.liu.se/Education/TMHL09/tables.pdf. [Använd 3 05 2015]. [16] Siemens, ”NX 10 Manual,” München , 2014. [17] N.-O. Eriksson och B. Karlsson, Verkstadshandboken, Stockholm: Liber AB, 2012. 59 60 8 Appendix A Översiktsritningar .................................................................................................................................. A-1 A.1 Skruvförband A ............................................................................................................................. A-1 A.2 Skruvförband B ............................................................................................................................. A-3 A.3 Skruvförband C ............................................................................................................................. A-5 B Arbetsflöde för kontaktanalyser i NX 10 ............................................................................................... B-7 C Analytisk modell för ett stelt skruvförband ........................................................................................... C-9 C.1 Kontaktpunktens rörelse i y-led .................................................................................................... C-9 C.2 Skruvkrafter................................................................................................................................. C-10 C.3 Momentjämvikt ........................................................................................................................... C-11 61 A Översiktsritningar Här presenteras översiktsritningar för CAD-modellerna som använts. Dessa ritningar saknar toleranser och detaljmått som en detaljritning hade haft. Dessa är därmed bara till för att få en överblick av skruvförbandens geometri och storlek. A.1 Skruvförband A A-1 A-2 A.2 Skruvförband B A-3 A-4 A.3 Skruvförband C A-5 A-6 B Arbetsflöde för kontaktanalyser i NX 10 Finita element analyser i NX 10 kan göras ytan större svårigheter, ofta fås korrekt svar utan några ändringar i solvern. Men om kontaktanalyser ska göras krävs det förändringar i solvern om man vill kunna använda vissa funktioner. Nedan presenteras hur arbetsflödet ser ut och några tips hur man undviker orimliga svar och fel. Texten är inriktad på skruvförband och använder de som exempel men metodiken går att applicera på andra modeller också. 1. 2. 3. 4. Skapa part-modellerna i NX. Detta gjordes genom en gemensam master-sketch som sedan WAVElänkades till olika parter. Skruven behöver splittas där förspänningskraften skall appliceras, vilket valdes till mutterns kontaktyta mot flänsytan. Assemblera parterna och constraina så att alla frihetsgrader är låsta. FEM: a. Byt till Advanced Simulations. b. Skapa en FEM and Simulation. Använd inte Assembly FEM då den ej stödjer kontaktanalyser. Välj solver, NX Nastran har använts i denna rapport. c. De länkade soliderna kan man t.ex. splitta om man vill använda symmetrier för att förkorta beräkningstider, gör sådana ändringar i den idealiserade filen så att assembly-filen kan användas till annat också. d. Skapa en Mesh Collector för varje part som ska ha olika storlekar och egenskaper på mesherna. Här kan man enkelt välja/skapa/ändra material också. e. Skapa ett 3D meshnät för soliderna. Klicka på 3D Tetrahedral Mesh med CTETRA(10). Börja med den rekommenderade elementstorleken för soliden. Bekräfta att meshen läggs i rätt Mesh Collector. f. Eftersom skruvarna är splittade för att få en yta för förspänningskraften måste de constrainas så att de sitter ihop igen. Mesh Mating Condition kan användas för detta, den kopplar ihop solider och deras meshnät med varandra. Men MMC är väldigt känslig för förändringar i partfilen, en ändring i partfilen kan göra att MMC slutar fungera. Använd inte denna för saker som ändras ofta, använd istället Surface-to-Surface Glue i Simulations istället. g. 3D-meshnätstypen CTETRA(10) tillåter inte nodrotationer. För att kunna använda momentlaster måste man därför skapa en punkt man kan applicera momentet på istället som 3D-meshen är beroende av. Välj 1D Connection med Point to face, som Source välj arcentrummet där momentet skall belastas. Som Target välj ytan som är i samma plan som punkten. Dessa element skall vara stela så i Element Properties välj RBE2. I Mesh Control välj samma yta och använd samma elementstorlek. Simulations: a. Ändringar i Solvern. i. För att visa kontaktresultat: Edit Solution -> Output requests -> Contact Result ->Enable BCRESULTS Request ii. Visa applicerade belastningar: Edit Solution -> Output requests -> Contact Result ->Enable OLOAD Request iii. Initial penetration är en viktig parameter, den bestämmer hur solvern skall bete sig om ytor penetrerar varandra initiellt. Det finns tre parametrar att välja mellan: Calculate from geometry, Ignore Penetration/Use Gaps, Set to zero. Dessa kan ge väldigt varierande lösningar vilket gör att det är viktigt att vara säker att rätt parameter är vald. Jämför aldrig med en uppskalade version av modellen då den kan ge vilseledande svar. I detta projekt har Initial Penetration/Use gaps använts: Edit Solution -> Global Contact Parameters -> Initial Penetration/gap -> Ignore Penetration/Use gaps. iv. Om man vill använda Enforced Displacements med förspänningskrafter, för att t.ex. böja en fläns med 30 grader för att se vad för moment det blir måste System Cell 581 ändras till 1 för inte få ett felmeddelande: B-7 b. c. d. e. Edit Solution -> System Cells -> Hitta 581 och ändra till 1. Skapa lasterna som skall användas i analysen. T.ex. Moment, Bolt pre-load. Moment kräver att det finns en nod som tillåter rotation. Därför gjordes I Bolt pre-load måste Force on 3D Elements användas för att förspänna 3D modeller, denna kräver en yta som är normal mot förspänningskraften. Därför bör skruven splittas där majoriteten av förspänningskrafterna uppstår. Notera att Bolt pre-load är känslig, om det uppstår fel kan det mycket väl vara den belastningen som är inblandad. Lås frihetsgrader så att modellen bara kan röra sig i de korrekta riktningarna, använd med fördel punkterna som skapades i 3.g. Definiera ytor mellan flänsarna och skruvarna. Använd Surface-to-Surface Contact för ytor som kan glida och ta i varandra, t.ex. fläns mot fläns och skruv mot fläns. Använd Surface-to-Surface Glue för ytor som ska ta i varandra och inte glida, t.ex. ytorna mellan skruvhuvud, mutter och fläns kommer inte glida mellan varandra nämnvärt. I Surface-to-Surface Contact/Glue måste en Source och Target yta defineras. Ytan som har finaste meshnät skall vara Source yta. Om de är lika stora gör det ingen skillnad vad som är Source och Target region. I Surface-to-Surface Contact måste en Max/Min Search Distance definieras så att NX vet hur långt bort ytorna kan vara, om värdena är för små kan ytorna gå igenom varandra för att NX inte hade tillräcklig räckvidd att se alla ytor. Ett enkelt sätt är skapa en variabel 𝑋 i Expressions och välja i Max/Min fälten 𝑋 och (−𝑋) respektive. Om flera olika lastfall med en förspänningslast ska testas är det viktigt att tänka på vilken ordning lasterna läggs på, detta eftersom förspänningslasterna är beroende av resultaten från tidigare subcase. Detta betyder att man bör lägga lasterna i ordning så att de ökar succesivt i storlek. Att använda WAVE-länkningar fungerar i allmänhet väldigt bra i Part- och Assembly-delen i NX 10, det är ett intuitivt sätt att kunna uppdatera dimensioner som är beroende av varandra i olika parter. Men i Advanced Simulations-modulen kan de konstanta uppdateringarna ställa till det. Speciellt är Mesh Mating Conditionfunktionen känslig för uppdateringar, denna slutar att fungera för små uppdateringar som inte hade något med parten alls att göra. Fördelen med den är att den limmar ihop splittade delar och de 2 ytorna som limmades ihop sammanfogas till 1 yta, detta gjorde det lätt att välja ytan som ska förspännas i Bolt pre-load. Man kan använda Surface-to-Surface Glue för detta också men då finns det ofta två ytor som man måste välja mellan vilket inte är helt intuitivt och detta tar en aning längre tid. Initial penetration är en viktig parameter i solvern. Var noga med att välja rätt. När de olika parametrarna skall jämföras måste det göras med en modell som inte är uppskalad. En uppskalad version kan i vissa fall se väldigt bra ut men fortfarande vara fel! Till exempel med Set to Zero ser det ofta väldigt bra ut med en uppskalad men det är fel för att ytor som inte ska vara i kontakt med varandra, det som bestämmer vad som är i kontakt för Set to Zero är Max/Min Search Distance. B-8 C Analytisk modell för ett stelt skruvförband Ett stelt skruvförband böjs så att skruvförbandet börjar öppna sig vid kopplingen. Skruvarna är förspända med kraften 𝐹𝑝𝑟𝑒 . Om antagandet att skruvförbandet är stelt är det enbart skruvarna som kommer att deformeras och skapa reaktionskrafter. Figur 67. Stelt skruvförband utsatt för moment i ändarna. Längden från flänsens centrum till skruvhålen är 𝑙ℎ = 𝑟 𝑐𝑜𝑠𝛼 (8.1) Längden från kontaktpunkten P till hålen kan skrivas som 𝑙𝑃−ℎ = 𝑃 − 𝑙ℎ (8.2) Där 𝑃 är längden till kontaktpunkten och 𝑙ℎ = 𝑟𝑐𝑜𝑠𝛼 som är längden mellan hålen i y-led. C.1 Kontaktpunktens rörelse i y-led Om flänsarna alltid tar i kontakt vid kontaktpunkten P kan man beskriva rörelsen i y-led som ↑ 𝑙1 sin(𝛽1 ) − 𝑅𝑐𝑜𝑠(𝛽1 ) = 𝑙2 sin(𝛽2 ) − 𝑅𝑐𝑜𝑠(𝛽2 ) (8.3) Den enda variabeln som är okänd är 𝛽2 , löser ut den genom att förenkla vänstra sidan till 𝑙1 sin(𝛽1 ) − 𝑅𝑐𝑜𝑠(𝛽1 ) = 𝑢 (8.4) Cos och sin-termerna på den högra sidan kan skrivas om med de trigonometriska halva vinkeln-sambanden 1 − tan2 (𝛽2 /2 ) (8.5) cos(𝛽2 ) = 1 + tan2 (𝛽2 /2 ) sin(𝛽2 ) = 2 tan(𝛽2 /2 ) 1 + tan2 (𝛽2 /2 ) (8.6) (8.6), (8.5) och (8.4) i (8.3) ger 𝑢 = 𝑙2 2 tan(𝛽2 /2 ) 1 − tan2 (𝛽2 /2 ) − 𝑅 1 + tan2 (𝛽2 /2 ) 1 + tan2 (𝛽2 /2 ) Vilket kan skrivas om som 𝑢(1 + tan2 (𝛽2 /2 )) = 2 𝑙2 tan(𝛽2 /2 ) − 𝑅(1 − tan2 (𝛽2 /2 )) 𝑢 + (𝑢 − 𝑅) tan2 (𝛽2 /2 ) = 2 𝑙2 tan(𝛽2 /2 ) − 𝑅 Sätter 𝑥 = tan(𝛽2 /2) i (8.9) ger en andragradsekvation 𝑢 + (𝑢 − 𝑅)𝑥 2 = 2𝑙2 𝑥 − 𝑅 C-9 (8.7) (8.8) (8.9) (8.10) 2𝑙2 𝑅+𝑢 𝑥− (𝑢 − 𝑅) 𝑢−𝑅 (8.11) 2𝑙2 2𝑙2 𝑅+𝑢 ±√ − 2(𝑢 − 𝑅) 2(𝑢 − 𝑅) 𝑢 − 𝑅 (8.12) 𝑥2 = 𝑥= Lutningen för den andra flänsen, 𝛽2 , blir 𝛽2 = 2arctan(𝑥) (8.13) C.2 Skruvkrafter Skruvarna har skruvats åt med en förspänning 𝐹𝑝𝑟𝑒 . När flänsarna separerar i punkten 𝑃 kommer skruven att deformeras vilket ökar skruvkraften betydligt. Figur 68. Skruvkrafter I Figur 68 kan man se att skruvkrafterna blir 𝐸𝑠 𝐴𝑠 𝐸𝑠 𝐴𝑠 𝐹𝑠 = 𝐹𝑝𝑟𝑒 + 𝑙ℎ sin 𝛽1 + 𝑙 sin 𝛽2 𝑡1 𝑡2 𝑃−𝐻𝑜𝑙𝑒𝑠 C-10 (8.14) C.3 Momentjämvikt I Figur 69 har en friläggning gjorts på de styva flänsarna. Figur 69. Friläggning av flänsar En kraftjämvikt på den vänstra flänsen i y-led ger 𝑜 = 𝑅𝑦1 + 𝑅𝑦2 ↑ (8.15) Och i x-led → (8.16) 𝑜 = 𝑅𝑥1 + 𝑅𝑥2 En momentjämvikt runt punkt 1 ger 𝑜 = 𝑀1 + 𝑅𝑥2 (𝑃 − 𝑙1 sin 𝛽1 ) + 𝑅𝑦2 (𝑙1 𝑐𝑜𝑠𝛽1 + 𝑃𝑠𝑖𝑛𝛽1 ) ↶ 1 + ∑ 𝐹𝑠1 (𝑙ℎ 𝑐𝑜𝑠𝛽1 − 𝑙1 𝑠𝑖𝑛𝛽1 ) (8.17) Löser ut reaktionsmomenten i (8.17) −𝑅𝑥2 (𝑃 − 𝑙1 sin 𝛽1 ) − 𝑅𝑦2 (𝑙1 𝑐𝑜𝑠𝛽1 + 𝑃𝑠𝑖𝑛𝛽1 ) = 𝑀1 + ∑ 𝐹𝑠1 (𝑙ℎ𝑜𝑙𝑒𝑠 𝑐𝑜𝑠𝛽1 − 𝑙1 𝑠𝑖𝑛𝛽1 ) Den högra flänsens kraftjämvikt: ↑ → (8.18) 𝑜 = −𝑅𝑦2 + 𝑅𝑦3 (8.19) 𝑜 = −𝑅𝑥2 + 𝑅𝑥3 (8.20) En momentjämvikt runt punkt 3 ger 𝑜 = 𝑀3 + 𝑅𝑥2 (𝑃 − 𝑙2 sin 𝛽2 ) + 𝑅𝑦2 (𝑙2 𝑐𝑜𝑠𝛽2 + 𝑃𝑠𝑖𝑛𝛽2 ) ↶ 3 − ∑ 𝐹𝑠2 (𝑙ℎ 𝑐𝑜𝑠𝛽2 − 𝑙2 𝑠𝑖𝑛𝛽2 ) (8.21) Löser ut reaktionsmomenten i (8.21) −𝑅𝑥2 (𝑃 − 𝑙1 sin 𝛽1 ) − 𝑅𝑦2 (𝑙1 𝑐𝑜𝑠𝛽1 + 𝑃𝑠𝑖𝑛𝛽1 ) = 𝑀1 + ∑ 𝐹𝑠1 (𝑙ℎ 𝑐𝑜𝑠𝛽1 − 𝑙1 𝑠𝑖𝑛𝛽1 ) (8.22) Reaktionsmomenten måste vara lika stora för flänsarna vilket gör det möjligt att lösa ut 𝑀1 om 𝑀1 = 𝑀3 − ∑ 𝐹𝑠2 (𝑙ℎ 𝑐𝑜𝑠𝛽2 − 𝑙2 𝑠𝑖𝑛𝛽2 ) − ∑ 𝐹𝑠1 (𝑙ℎ 𝑐𝑜𝑠𝛽1 − 𝑙1 𝑠𝑖𝑛𝛽1 ) (8.23) 𝑀1 = 2 C-11
© Copyright 2024