1 Wyznaczenie naprężeń cieplnych w rurze, przez którą przepływa medium o temperaturze 400 C Zadanie 4-5 Cel: Zapoznanie studentów z modelowaniem zjawisk przepływu ciepła i określaniem warunków brzegowych i początkowych, kontakt cieplny i mechaniczny. Wykorzystanie zmiennego w czasie pola temperatury do określenia zmian stanu naprężenia w elemencie konstrukcji. Literatura: 1. R. Grądzki: Wprowadzenie do metody elementów skończonych, Politechnika Łódzka, 2002 2. W. Śródka: Trzy lekcje metody elementów skończonych, Politechnika Wrocławska, 2004 3. A. Skrzat: Modelowanie liniowych i nieliniowych problemów mechaniki ciała stałego i przepływów ciepła w programie Abaqus, Rzeszów 2010 4. Z. Orłoś: Naprężenia cieplne PWN 1991 5. Wykłady ! Komputerowa optymalizacja konstrukcji odlewu pod względem wytrzymałościowym
2 Właściwości cieplne i mechaniczne stopu aluminium Współczynnik przewodzenia ciepła = 206 W/m o K Ciepło właściwe c p = 909 J/kg o K Gęstość = 2700 kg/m 3 Współczynnik rozszerzalności liniowej = 23,8*10 -6 o K -1 Moduł sprężystości E = 7,17*10 10 Pa Liczba Poissona = 0,33
3 (0,0) (0.04,0 ) (0.06,0) (0.04,0.2) (0.06,0.2) 400 C /10000 20/500 Model temperaturowy 400 C 20 C Rys 1. Model zjawiska Rys. 2. Model geometryczny
4 Property Thermal - Conductivity = 206 Thermal – Specific Heat = 909 General – Mass Density = 2700 Name = Aluminium Step Heat Transfer Transient Time Period = 100 Incrementation – Maximum Number of Increments = 100 Type = Fixed Increment Size = 1 Ciąg wprowadzanych instrukcji dla modelu temperaturowego
5 Interaction Create – Step-1 – Surface Film Condition Film coefficient = 500 Sink temperature = 20 Create – Step-1 – Surface Film Condition Film coefficient = 10000 Sink temperature = 400
6 Load Create Predefined Field Step – Initial Category - Other – Temperature = 20 Seed = 0.005 Element Type – Quadratic – Heat Transfer Mesh Job = Temperatura
7 Model naprężeniowy Rys. 3. Model geometryczny – osiowosymetryczny z określonym mocowaniem końców rury (0,0) U2 = 0
8 File – Save As = Naprezenia W drzewie zadania rozwinąć zakładkę Materials (1) i kliknąć 2X materiał o nazwie Aluminium Elastic E=7.17e10, v=0.33 Expansion a= 23.8e-6 Edycja właściwości materiałowych Edycja rodzaju analizy W drzewie zadania rozwinąć zakładkę Steps (2) i skasować Step – 1, Następnie utworzyć nowy, klikając 2X zakładkę Step (1), wybrać : Static, General Ciąg wprowadzanych instrukcji dla modelu naprężeniowego
9 Definiowanie zbioru wyników W drzewie zadania rozwinąć zakładkę Field Output Requests (1) i 2X kliknąć F-Output-1, odznaczyć Forces/Reactions i Contact. Zaznaczyć Thermal - NT Warunki brzegowe Odebrać stopnie swobody U2 końcom rury
10 W drzewie zadania rozwinąć zakładkę Predefined Fields (1), skasować istniejące Pole Predefined Fields – 1. W jego miejsce utworzyć nowe pole. 2X kliknąć Predefined Fields. W oknie dialogowym wprowadzić Step – Step-1 Category – Other Types for Selected Step – Temperature Distribution – From results or output database file File name – Temperatura Begin step = 0 Begin increment = 1 End step = 1 End increment = 100 Interpolation – Mesh compatibility = Compatible
11 Minimum teorii naprężenia cieplne Rys. 4. Stany naprężenia w belce grzanej jednostronnie
12 Nieliniowe pole temperatury T=T(y) rys. 4a powoduje powstanie pola naprężeń, które w przypadku braku więzów na bocznych krawędziach belki, wywołuje jej swobodne ugięcie rys 4b; naprężenie jest w tym przypadku opisane zależnością : W przypadku braku możliwości swobodnego przemieszczania się końców belki naprężenie jest opisane wzorem (2) rozbudowanym o człony związanie z hamowaniem ugięcia rys.4c i hamowaniem rozszerzalności rys. 4d przez narzucone więzy: (1) (2)