WykÅ‚ad 1 Zagadnienie brzegowe liniowej teorii sprężystoÅ›ci. Metody rozwiÄ…zywania, metody wytrzymaÅ‚oÅ›ci materiałów. Zestawienie wzorów i okreÅ›leÅ„. UkÅ‚ad współrzÄ™dnych KartezjaÅ„ski, prostokÄ…tny. Osie x y z oznaczono odpowiednio przez x1, x2, x3. t x3 St x2 P &! &! &! &! x1 Su Rozpatrujemy obszar &! z brzegiem S podzielonym na dwie rozÅ‚Ä…czne części Su i St. Obszar ten wypeÅ‚niony jest ciÄ…gÅ‚ym oÅ›rodkiem materialnym M. Dowolny punkt P w obszarze &! ma współrzÄ™dne xP. Wektor poÅ‚ożenia punktu P zapisywany jest również przy użyciu nastÄ™pujÄ…cej konwencji: P P P P x = () x1 ,x2 ,x3 lub w zapisie indeksowym: xiP , gdzie i=1, 2, 3. Na części brzegu St. zadany jest wektor naprężenia t0 o skÅ‚adowych t0i (to znaczy: t01,. t02, t03). Pod dziaÅ‚aniem tego obciążenia punkty ciaÅ‚a M zmienia swoje poÅ‚ożenie w przestrzeni. Obszar &! zmienia ksztaÅ‚t. Nowe poÅ‚ożenia P' punktów ciaÅ‚a M nazywamy jego konfiguracjÄ… odksztaÅ‚conÄ…. Dowolny punkt P zajmuje w konfiguracji odksztaÅ‚conej nowe poÅ‚ożenie xP': P' P' P' P' x = () x1 ,x2 ,x3 Różnica wektorów tych dwóch poÅ‚ożeÅ„ punktu P opisana jest przez wektor przemieszczenia uP(xP): uP(xP)=xP'-xP W zdeformowanym obszarze &! zdefiniowane jest tym samym pole wektora przemieszczenia u(x) o skÅ‚adowych skalarnych (u1(x), u2(x), u3(x)). Pole to powinno być "zgodne" z danym przemieszczeniem u0(x) zadanym na części Su brzegu &!. To znaczy, że: dla x‚"Su u(x)=u0(x) gdzie u0(x) jest dane OdksztaÅ‚cenie nieskoÅ„czenie maÅ‚ego elementu d&! opisane jest polem tensora maÅ‚ego odksztaÅ‚cenia µ µ: µ µ 1 ëÅ‚ "u (x)öÅ‚ "ui (x) 1 j 1 ìÅ‚ ÷Å‚ µij (x)= + lub zapisujÄ…c skrótowo: µij = ("ui, j + "u ) j,i ìÅ‚ 2 "x "xi ÷Å‚ 2 j íÅ‚ Å‚Å‚ SkÅ‚adowe tensora maÅ‚ego odksztaÅ‚cenia reprezentowane sÄ… przez macierz pól skalarnych: îÅ‚µ11 µ12 µ13 Å‚Å‚ ïÅ‚µ śł µ = µ µ 21 22 23 ïÅ‚ śł ïÅ‚µ31 µ32 µ33 śł ðÅ‚ ûÅ‚ MyÅ›lowo wyodrÄ™bniony z obszaru &! element d&! oddziaÅ‚uje z resztÄ… obszaru, poprzez element powierzchni dS, za poÅ›rednictwem wektora naprężeÅ„ t. SkÅ‚adowe wektora naprężeÅ„ t zależą nie tylko od punktu M w poÅ‚ożeniu odksztaÅ‚conym P' ale także od orientacji przestrzennej element powierzchni dS poprowadzonej przez P'. Ta orientacja zadana jest za poÅ›rednictwem n - jednostkowego wektora normalnego do dS. Wektor ten ma współrzÄ™dne nastÄ™pujÄ…ce: n={cos(n,Ox1), cos(n,Ox2), cos(n,Ox3)} WielkoÅ›ciÄ… charakteryzujÄ…cÄ… stan naprężenia w punkcie ciaÅ‚a M w jego poÅ‚ożeniu P' jest tensor naprężenia Ä. Jest on zwiÄ…zany z wektorem naprężenia w punkcie materialnym M Ä Ä Ä zależnoÅ›ciÄ…: 3 "Ä ij (x)n j (x)= ti (x,n) j=1 Interpretacja i umowa o znakach dla macierzy skÅ‚adowych tensora naprężenia znana jest z dotychczasowego kursu wytrzymaÅ‚oÅ›ci materiałów. îÅ‚Ä11 Ä12 Ä13 Å‚Å‚ ïÅ‚Ä Å›Å‚ Ä = Ä Ä 21 22 23 ïÅ‚ śł ïÅ‚Ä31 Ä 32 Ä33 śł ðÅ‚ ûÅ‚ Warunki równowagi zapisane dla element d&! prowadzÄ… do szeÅ›ciu równaÅ„ (b jest wektorem siÅ‚ masowych): Ä = Ä (symetria tensora naprężenia) ij ji 3 "Ä "Ä "Ä i1 i2 i3 dla i=1,2,3: + + + bi = 0 lub skrótowo: + bi = 0 "Ä ij, j "x1 "x2 "x3 j=1 ZwiÄ…zek pomiÄ™dzy tensorem odksztaÅ‚cenia a tensorem naprężenia jest wÅ‚aÅ›ciwoÅ›ciÄ… materiaÅ‚u i powinien być okreÅ›lony doÅ›wiadczalnie. Liniowa teoria sprężystoÅ›ci zakÅ‚ada, że zwiÄ…zek ten opisany jest uogólnionym prawem Hooke'a. Zapis tego prawa konstytutywnego, dla przypadku materiaÅ‚u izotropowego podany jest poniżej: 3 1 +½ ½ µij = Ä + ´ "Ä ij ij kk E E k =1 2 1 dla i = j Å„Å‚ pomocniczy symbol ´ okreÅ›la siÄ™ nastepujÄ…co: ´ = ´ ´ ´ òÅ‚ ij ół0 dla i `" j Zależność µ Ä) można Å‚atwo "odwrócić", obliczajÄ…c Ä µ: µ(Ä Ä w funkcji µ µ Ä Ä µ µ Ä Ä µ 3 E ëÅ‚ ½ öÅ‚ Äij = ìÅ‚µij + ´ "µ ÷Å‚ ìÅ‚ ij kk ÷Å‚ 1 +½ 1 - 2½ íÅ‚ k=1 Å‚Å‚ Ramka nr 1 zawiera sformulowanie zadania liniowej teorii spreżystoÅ›ci: Ramka 1. SformuÅ‚owanie liniowego problemu teorii sprężystoÅ›ci: Znalezć wektorowe pole przemieszczeÅ„ u(x), oraz tensorowe pole naprężeÅ„ Ä Ä(x), takie Ä Ä Å¼e: u(x)=u0(x) na części brzegu Su naprężenia Ä(x) speÅ‚niajÄ… naprężeniowe warunki brzegowe (równowaga elementu w Ä Ä Ä sÄ…siedztwie brzegu St): 3 na części brzegu St "Ä ij (u)n j = ti0 j=1 oraz warunki równowagi w caÅ‚ym obszarze &!: 3 Ä = Ä "Ä ij, j + bi = 0 ij ji j=1 Ponadto speÅ‚nione jest prawo Hooke'a, które wiąże tensor naprężenia z wektorem przemieszczenia za poÅ›rednictwem tensora maÅ‚ych odksztaÅ‚ceÅ„: 3 E ëÅ‚ ½ 1 Ä = ìÅ‚µij (u)+ ´ (u)öÅ‚ gdzie µij (u)= ("ui, j + "u ) "µ ÷Å‚ ij ìÅ‚ ij kk ÷Å‚ j,i 1 +½ 1 - 2½ 2 íÅ‚ k =1 Å‚Å‚ Prawo Hooke'a może być równie dobrze zadane wzorem: 3 1 +½ ½ µij = Ä + ´ "Ä ij ij kk E E k =1 (we wszystkich równaniach i,j,k=1,2,3) Dwa sformuÅ‚owania pochodne: przemieszczeniowe i naprężeniowe - pozwalajÄ… zredukować liczbÄ™ niewiadomych pól skalarnych. Zapis µ(u) w formuÅ‚ach w ramce nr 2. oznacza, że naprężenia zostaÅ‚y wyrażone przez µ µ µ przemieszczenia u. Podobnie zapis Ä(u) oznacza, że naprężenia zostaÅ‚y wyrażone przez Ä Ä Ä przemieszczenia u i w powyższym sformuÅ‚owaniu sÄ… tylko trzy niewiadome funkcje skalarne: u1(x), u2(x), u3(x), nie zaÅ› 9, jak poprzednio. 3 Ramka 2. Przemieszczeniowe sformuÅ‚owanie liniowego problemu teorii sprężystoÅ›ci: Znalezć wektorowe pole przemieszczeÅ„ u(x), takie że: u(x)=u0(x) na części brzegu Su i takie, że skÅ‚adowe tensora naprężenia obliczone za poÅ›rednictwem zwiÄ…zku konstytutywnego: E 1 ëÅ‚µ ½ ´ 3 µ gdzie µij Ä (u)= (u)+ (u)öÅ‚ (u)= ("ui, j + "u ) ìÅ‚ " ÷Å‚ ij ij ij kk j,i 2 1+½ 1- 2½ íÅ‚ k =1 Å‚Å‚ speÅ‚niajÄ… warunki równowagi: 3 w caÅ‚ym obszarze &!: Ä = Ä "Ä ij, j + bi = 0 ij ji j=1 oraz naprężeniowe warunki brzegowe (równowaga elementu w sÄ…siedztwie brzegu St): 3 na części brzegu St "Ä ij (u)n j = ti0 j=1 Poniżej przedstawiono przykÅ‚ad sformuÅ‚owania przemieszczeniowego. Jest to ukÅ‚ad równaÅ„ Naviera. Jego wyprowadzenie wymaga wykonania nastÄ™pujÄ…cych operacji: " u µ Ä=Ä µ µ(u) prawo Hooke'a Ä Ä( µ µ Ä Ä µ(u)) µ Ä Ä µ " Podstawienie Ä( µ Ä µ(u)) do równaÅ„ równowagi. Ä µ Ä µ Otrzyma siÄ™ w ten sposób równania równowagi wyrażone przez przemieszczenia: 3 "Ä ij, j (u)+ bi = 0 j=1 3 3 E ½ (u)+ (u)öÅ‚ + bi = 0 ìÅ‚ "1 +½ ëÅ‚µij, j 1 - 2½ ´ ij "µ kk , j ÷Å‚ j=1 íÅ‚ k =1 Å‚Å‚ 3 3 E 1 E ½ (ui, + u )+ uk ,ki + bi = 0 "1+½ 2 jj j,ij 1+½ 1- 2½ " j=1 k =1 3 3 E 1 E ½ E 1 ui, + + uk + bi = 0 "1+½ 2 jj ëÅ‚1+½ 1- 2½ 1+½ 2 öÅ‚ ,ki ìÅ‚ ÷Å‚" j=1 íÅ‚ Å‚Å‚ k=1 (zauważmy, że nie ma to znaczenia czy wskaznik sumowania od 1 do 3 nazywa siÄ™ j czy k) 3 3 1 2(1 +½ )b = 0 ui, jj + uk ,ki + "" i (1 - 2½ ) E j=1 k=1 Jest to ukÅ‚ad trzech równaÅ„ różniczkowych (kolejno dla i=1, 2, 3) z drugimi pochodnymi czÄ…stkowymi, okreÅ›lony w obszarze &! i z warunkami brzegowymi na Su. SformuÅ‚owanie to wygodnie jest stosować jeÅ›li warunki brzegowe zadane sÄ… jedynie na Su. Ewentualne warunki brzegowe naprężeniowe powinny być wyrażone w przemieszczeniach w nastÄ™pujÄ…cy sposób: 3 E ½ ëÅ‚ (ui, j + ui, j )n + uk ,k ni öÅ‚ = pi ìÅ‚ " ÷Å‚ j 1+½ 1- 2½ íÅ‚ k =1 Å‚Å‚ Problem liniowej teorii sprężystoÅ›ci w ujÄ™ciu przemieszczeniowym można teraz sformuÅ‚ować tak: 4 Ramka 2.1. Przemieszczeniowe sformuÅ‚owanie liniowego problemu teorii sprężystoÅ›ci - równania Naviera: Znalezć wektorowe pole przemieszczeÅ„ u(x), takie że: ui=u0i na części brzegu Su 3 E ½ ëÅ‚ oraz (ui, j + ui, j )n + uk ,k ni öÅ‚ = pi na części brzegu St ìÅ‚ " ÷Å‚ j 1+½ 1- 2½ íÅ‚ k =1 Å‚Å‚ i takie, że skÅ‚adowe wektora przemieszczenia speÅ‚niajÄ… ukÅ‚ad równaÅ„ różniczkowych o pochodnych czÄ…stkowych (równania Naviera): 3 3 1 2(1 +½ )b = 0 ui, jj + uk ,ki + "" i (1 - 2½ ) E j=1 k=1 SkÅ‚adowe tensora naprężenia można obliczyć za poÅ›rednictwem zwiÄ…zku konstytutywnego: 3 E ëÅ‚ ½ Ä (u)= ìÅ‚µij (u)+ ´ (u)öÅ‚ "µ ÷Å‚ ij ìÅ‚ ij kk ÷Å‚ 1 +½ 1 - 2½ íÅ‚ k=1 Å‚Å‚ 1 gdzie µij (u)= ("ui, j + "u ) zaÅ› u jest rozwiÄ…zaniem równaÅ„ różniczkowych Naviera. j,i 2 (OczywiÅ›cie naprężenia takie automatycznie speÅ‚niajÄ… równania równowagi oraz naprężeniowe warunki brzegowe.) Aby uzyskać sformuÅ‚owanie naprężeniowe liniowego problemu teorii sprężystoÅ›ci, należy wyeliminować przemieszczenia ze sformuÅ‚owania wyjÅ›ciowego. W tym celu należy zróżniczkować dwukrotnie µ11 wzglÄ™dem x2 zaÅ› µ22 wzglÄ™dem x1 nastÄ™pnie dodać dodać stronami oba rownania: µ11,22 + µ = u1,122 + u2,211 = Å‚ 22,11 12,12 Podobnie można uzyskać dwa kolejne wzory: µ11,33 + µ33,11 = Å‚ 13,13 µ + µ33,22 = Å‚ 22,33 23,23 Trzy nastÄ™pne otrzymuje sie przez wykonanie podobnych przeksztaÅ‚ceÅ„ na równaniach definiujÄ…cych Å‚ Å‚: Å‚ Å‚ Å‚ = u1,23 + u2,13 12,3 Å‚ = u1,32 + u3,12 13,2 Å‚ = u2,31 + u3,21 23,1 Å‚12,3 + Å‚ - Å‚ = u1,23 + u1,32 13,2 23,1 Å‚12,31 + Å‚ - Å‚ = u1,123 + u1,132 = 2µ11,23 13,21 23,11 5 podobnie otrzymuje siÄ™ dwa poniższe wzory: Å‚ + Å‚ - Å‚ = 2µ 21,32 23,12 13,22 22,13 Å‚ + Å‚ - Å‚ = 2µ33,12 32,13 31,23 12,33 Równania powyższe nazywane sÄ… warunkami nierozdzielnoÅ›ci odksztaÅ‚ceÅ„. Ich speÅ‚nienie daje gwarancje znalezienia ciÄ…gÅ‚ego i różnowartoÅ›ciowego pola przemieszczeÅ„ u, takiego, że: 1 µij (u)= (ui, j + u ) j,i 2 Ramka nr 3. Naprężeniowe sformuÅ‚owanie liniowego problemu teorii sprężystoÅ›ci: Znalezć tensorowe pole naprężeÅ„ Ä Ä(x), takie że: Ä Ä speÅ‚nione sÄ… warunki równowagi: 3 w caÅ‚ym obszarze &!: Ä = Ä "Ä ij, j + bi = 0 ij ji j=1 oraz naprężeniowe warunki brzegowe 3 na części brzegu St "Ä ij (u)n j = ti0 j=1 OdksztaÅ‚cenia wyznaczone dla tych naprężeÅ„ za poÅ›rednictwem prawa Hooke'a: 3 1 +½ ½ µij (Ä)= Ä + ´ "Ä ij ij kk E E k =1 powinny speÅ‚niać równania nierozdzielnoÅ›ci w caÅ‚ym obszarze &!: µ11,22 + µ = Å‚ 22,11 12,12 µ11,33 + µ33,11 = Å‚ 13,13 µ + µ33,22 = Å‚ 22,33 23,23 Å‚12,31 + Å‚13,21 - Å‚ = 2µ11,23 23,11 Å‚ + Å‚ - Å‚ = 2µ 21,32 23,12 13,22 22,13 Å‚ + Å‚ - Å‚ = 2µ33,12 32,13 31,23 12,33 Zauważmy, że ponieważ µ µ Ä), w równaniach nierozdzielnoÅ›ci w ramce nr 3. wystÄ™pujÄ… µ= µ Ä µ µ(Ä µ µ Ä tylko pola naprężeÅ„ speÅ‚niajÄ…ce równania równowagi i warunki naprężeniowe na części brzegu St (takie pola naprężeÅ„ nazywamy statycznie dopuszczalnymi). Cwiczenie: Przedstawić wszystkie równania wystÄ™pujÄ…ce w wykÅ‚adzie tak, aby wystÄ™powaÅ‚y w nich jawnie wielkoÅ›ci Ãx, Ãy, Ãz, Äxy, Äxz, Äzy, µx, µy, µz, Å‚xy, Å‚xz, Å‚zy, u, v, w. Porównać otrzymany zapis z tym jaki jest podany w rozdziale 17 (Zadania i podstawowe równania teorii sprężystoÅ›ci) podrÄ™cznika "WytrzymaÅ‚ość Materiałów" (A. Jakubowicz, Z. OrÅ‚oÅ›). 6