Wprowadzenie do rachunku tensorowego


Wprowadzenie do rachunku tensorowego 2007-02-05
Wprowadzenie do rachunku tensorowego
Tensorami nazywamy ogólnie takie obiekty jak skalary, wektory i wielkości wyższych
rzedów. Tak więc tensorem rzędu zerowego (np. A, czyli wielkość o walencji 0) są skalary,
pierwszego (np. ji , czyli wielkość o walencji 1) wektory, zaś reprezentacją tensora rzędu
drugiego (np. dij ) jest macierz 2x2. Tensory rzędów (czyli walencji) wyższych (np. 3-go, 4-
go) nie możemy już zapisać na ,,płaszczyznie  .
Operacje na tensorach rządzą się pewnymi zasadami (algebra tensorów):
1. Cij = Aij + Bij - dodawanie tensorów
2. aijbkln = cijkln - iloczyn tensorowy dwóch tensorów jest tensorem o walencji równej
sumie walencji składników
3. Istnieje tzw. tensor jednostkowy  zwany deltą Kroneckera
1 0 0
Ą# ń#
1 dla i = j
ż#
ó#0
ij = 1 0Ą# , czyli ij =
#0 dla i `" j
ó# Ą#
#
ó# Ą#
Ł#0 0 1Ś#
4. Zamianę wskazników dokonujemy poprzez przemnożenie danego tensora przez tensor
jednostkowy
aijik = aik
5. Niezmienniki tensora:
I = Aii
1
Aij ! II = (Aii Ajj - Aij A )
ji
2
III = det(Aij )
Równania tensorowe
Jednym z ,,zastosowań  tensorów jest umożliwienie zapisu relacji funkcyjnej między parą
wektorów, czyli np. dla wektorów ai i bj będzie
ai = Bijbj ,
natomiast dla pary tensorów będzie
Aij = BijklCkl
Konwencja sumacyjna Einsteina
W zapisie równań tensorowych powszechnie wykorzystuje się tzw. notację sumacyjną
Einsteina, pozwalającą zapisać skomplikowane wyrażenia w zwięzłej postaci. Konwencja ta
mówi, że sumowania dokonujemy po powtarzających się wskaznikach, np.
3
aij x a" x = ai1x1 + ai2 x2 + ai3x3
j "aij j
j=1
Te powtarzające się wskazniki nazywa się wskaznikami niemymi. Mają one tą właściwość,
że można je wymienić na dowolne inne, np.:
aij x a" aik xk
j
UWAGA! Przy przekształceniach należy zwrócić uwagę na to, żeby dany wskaznik nie
powtórzył się więcej niż 2 razy po jednej stronie równania.
A. Marynowicz Strona 1
Wprowadzenie do rachunku tensorowego 2007-02-05
Przykłady zadań z równaniami tensorowymi
Poniżej podam kilka przykładów, bez wyjaśnienia podstaw fizycznych  po te odsyłam do
wykładów.
Przykład 1.
Pole przemieszczeń w nieograniczonym ciele sprężystym ma postać:
ui = xi Bj x (1)
j
Należy wyznaczyć tensor odkształceń w postaci:
2ij = ui, j + u (2)
j,i
oraz naprężeń
 = 2Gij + kkij (3)
ij
a następnie sprawdzić równanie równowagi wewnętrznej
 + Fi = 0 (4)
ij, j
Z równania tego należy, przyjmując znane Fi , wyznaczyć składowe wektora Ai.
Rozwiązanie
Aby to zadanie rozwiązać należy w pierwszej kolejności obliczyć pochodne (cząstkowe)
z wektora przemieszczeń ui. Pochodną (cząstkową) w zapisie tensorowym oznacza się
przecinkiem i wskazuje po jakiej współrzędnej się ją liczy  to po przecinku, np.
"xi "2xi
a" xi, j , a" xi, , itp&
"x "x "xk jk
j j
Należy jeszcze pamiętać o pochodnej mieszanej (z iloczynu funkcji). Mamy więc w naszym
przypadku (po zamianie indeksu j na k w iloczynie!  powód podano wyżej):
ui, j = [xi Bk xk ] = [xi ] Bk xk + xi[Bk xk ] = xi, j Bk xk + xi Bk xk , j = ij Bk xk + xi Bk =
kj
, j , j , j
= ij Bk xk + xi Bj
(5)
Wektor Bi jest wektorem stałych, stąd pochodna jego jest równa zero. Wykorzystano tu także
zależności (można to łatwo rozpisać)
xi, j = ij
oraz
Bk = Bj .
kj
Pochodną uj,i otrzymujemy podobnie, zamieniając wskazniki i oraz j, stąd otrzymamy
u = ij Bk xk + x Bi (6)
j,i j
wiedząc, że ij =  .
ji
Podstawiając (5) i (6) do (2) otrzymamy:
2ij = 2ij Bk xk + xi Bj + x Bi . (7)
j
Następnie podstawiamy to do równania fizycznego (3)
 = G(2ij Bk xk + xiBj + x Bi)+ 4ij Bk xk = K =
ij j
(8)
= G(xiBj + x Bi)+ ij Bk xk (2G + 4)
j
A. Marynowicz Strona 2
Wprowadzenie do rachunku tensorowego 2007-02-05
Występujące w (3) wyrażenie kk otrzymano z przekształcenia (7):
2kk = 2 Bl xl + xk Bk + xk Bk = 2 " 3" Bl xl + 2xk Bk = 6Bk xk + 2xk Bk = 8Bk xk ,
kk
czyli
kk = 4Bk xk .
Chcąc otrzymać ostatni element, czy wyliczyć siłę masową z równania równowagi, musimy
obliczyć pochodną z równania (8), czyli
ij, j = [G(xiBj + xjBi)] +[ijBk xk(2G + 4)] =
, j , j
 
ij # kj
3
# ś#
} } }ś#
ź#
= Gś# xi, B + xj, j Bi ź# + (2G + 4)ś#ij Bk x =
ś#
ś# ź#
1j3j 1 3j ź#
2 2k , ź#
ś# ź#
ś#
Bi B
# # j
# #
= 4GBi + (2G + 4)Bi = Bi(6G + 4)
Wykorzystano tu fakt, że pochodna z ij jest równa zero. Mając tak policzoną pochodną
otrzymamy równanie równowagi
Fi = -Bi(6G + 4) , (9)
z którego, przy znanym Fi , wyznaczymy składowe wektora Bi :

B1 = - F1
6G + 4

B2 = - F2 (10)
6G + 4

B3 = - F3
6G + 4
Równania (7), (8), (9) i (10) stanowią rozwiązanie zadania.
Przykład 2.
Pole przemieszczeń w nieograniczonym ciele sprężystym ma postać:
ui = (xi - Bi )xj Bj (11)
Wyznaczyć wielkości z przykładu 1 (tylko zamiast Ai wyznaczyć składowe Bi).
Rozwiązanie
Możemy na początku przekształcić (zamieniając indeksy ,,nieme  (czyli nie występujące po
lewej stronie równania) j na k) (11)
ui = (xi - Bi )x Bj = xi xj Bj - xiBiBj = xi xk Bk - xk BiBk
j
Mamy więc:
k , j ij kj
# ś#
} } }ź#
ś#
ui, j = [xi xk ] Bk - xk j Bk Bi = xi, j xk + xi xk , j ź#Bk - Bj Bi = ij xk Bk + xikjBk - Bj Bi =
, j
ś#
1, 3 13
2 2
ś# ź#
Bj # # Bj
= ij xk Bk + Bj(xi - Bi )
oraz
u = ij xk Bk + Bi(xj - Bj).
j,i
Tak więc tensor odkształcenia ma postać
2ij = 2ij xk Bk + 2BiBj + Bj xi + Bi x (12)
j
A. Marynowicz Strona 3
Wprowadzenie do rachunku tensorowego 2007-02-05
Tensor kk :
2
kk = 4Bk xk - Bk
Stąd tensor naprężeń (równanie fizyczne):
2
 = (2G + 4)ij xk Bk + GBi(xi + x )- ij Bk - 2GBiBj (13)
ij j
Pochodna z wyrażenia (13) ma postać (uwaga na xj,j=3)
2
 = (2G + 4)[ij xk Bk] + GBi[(xi + xj)] -[ij Bk ] -[2GBiBj] = (2G + 4)Bi+G(Bj - Bi)
,
ij, j
1 3 14243 1 3 1424, j
424, j 4 4, j 424, j 3
Bi G(Bj -3Bi ) 0 0
skąd otrzymamy wyrażenie na składowe wektora siły masowej
Bi
ij
}
Fi = -Bi(- G + 4)- G Bj = -Bi[G(-1+ ij)+ 4]= -4Bi (14)
Ostatecznie składowe wektora Bi wyliczymy przekształcając (14), czyli otrzymamy

-1
Bi = -[G(-1+ ij)+ 4] Fi = - Fi (15)
4
czyli po rozpisaniu:
-1
B1 = -(4) F1
-1
B2 = -(4) F2 (16)
-1
B3 = -(4) F3
co kończy zadanie.
A. Marynowicz Strona 4


Wyszukiwarka

Podobne podstrony:
Marynowicz A Wprowadzenie do rachunku tensorowego v4
WYKŁAD 1 Wprowadzenie do biotechnologii farmaceutycznej
Medycyna manualna Wprowadzenie do teorii, rozpoznawanie i leczenie
01 Wprowadzenie do programowania w jezyku C
wprowadzenie do buddyzmu z islamskiego punktu widzenia
1 wprowadzenie do statystyki statystyka opisowa
Informatyka Wprowadzenie Do Informatyki Ver 0 95
Wprowadzenie do psychologii wykł UG
645 Informacja dodatkowa wprowadzenie do sprawozdania finasowego

więcej podobnych podstron