Formelsammlung Tensoranalysis

Aus testwiki
Version vom 15. August 2022, 13:53 Uhr von imported>Boehm (typog)
(Unterschied) ← Nächstältere Version | Aktuelle Version (Unterschied) | Nächstjüngere Version → (Unterschied)
Zur Navigation springen Zur Suche springen

Vorlage:Hinweisbaustein

Diese Formelsammlung fasst Formeln und Definitionen der Analysis mit Vektor- und Tensorfeldern zweiter Stufe in der Kontinuumsmechanik zusammen.

Allgemeines

Siehe auch

Formelsammlung Tensoralgebra

Nomenklatur

Kronecker-Delta

Vorlage:Siehe auch

δij=δij=δij=δji={1falls i=j0sonst

Permutationssymbol

Vorlage:Siehe auch

ϵijk=e^i(e^j×e^k)={1falls(i,j,k){(1,2,3),(2,3,1),(3,1,2)}1falls(i,j,k){(1,3,2),(2,1,3),(3,2,1)}0sonst, d. h. bei doppeltem Index

Kreuzprodukt:

aie^i×bje^j=ϵijkaibje^k
ϵijke^k=e^i×e^j

Formelsammlung Tensoralgebra#Kreuzprodukt eines Vektors mit einem Tensor:

(a×𝐀)g:=a×(𝐀g)
b(a×𝐀)=(b×a)𝐀
g(𝐀×a):=(g𝐀)×a
(𝐀×a)b=𝐀(a×b)

Basisvektoren

Kartesische Koordinaten

Vorlage:Siehe auch

x1,x2,x3

mit Basisvektoren

e^1=(100),e^2=(010),e^3=(001)

die Standardbasis oder allgemeiner eine beliebige Orthonormalbasis ist.

Zylinderkoordinaten

Vorlage:Siehe auch

e^ρ=(cos(φ)sin(φ)0),e^φ=(sin(φ)cos(φ)0),e^z=(001)
e^ρ,φ=e^φ,e^φ,φ=e^ρe^z,φ=0

Winkelgeschwindigkeit#Zylinderkoordinaten:

ω=φ˙e^ze^˙ρ/φ/z=ω×e^ρ/φ/z

Kugelkoordinaten

Vorlage:Siehe auch

e^r=(sin(ϑ)cos(φ)sin(ϑ)sin(φ)cos(ϑ)),e^ϑ=(cos(ϑ)cos(φ)cos(ϑ)sin(φ)sin(ϑ)),e^φ=(sin(φ)cos(φ)0)

Winkelgeschwindigkeit#Kugelkoordinaten:

ω=(ϑ˙sin(φ)ϑ˙cos(φ)φ˙)=φ˙cos(ϑ)e^rφ˙sin(ϑ)e^ϑ+ϑ˙e^φe^˙r/ϑ/φ=ω×e^r/ϑ/φ

Krummlinige Koordinaten

Vorlage:Siehe auch

y1,y2,y3
bi=xyi,bi=grad(yi)=yixbibj=δij

Ableitung von Skalar-, Vektor- oder Tensorfunktionen

Gâteaux-Differential

Vorlage:Siehe auch

Df(x)[h]:=ddsf(x+sh)|s=0=lims0f(x+sh)f(x)s

mit s, f,x,h skalar-, vektor- oder tensorwertig aber x und h gleichartig.

Produktregel:

D(f(x)g(x))[h]=Df(x)[h]g(x)+f(x)Dg(x)[h]

Kettenregel:

Df(g(x))[h]=Df(g)[Dg(x)[h]]

Fréchet-Ableitung

Vorlage:Siehe auch

Existiert ein beschränkter linearer Operator 𝒜, sodass

𝒜[h]=Df(x)[h]h

gilt, so wird 𝒜 Fréchet-Ableitung von f nach x genannt. Man schreibt dann auch

fx=𝒜.

Ableitung von Potenzen eines Tensors

(𝐓1)˙=𝐓1𝐓˙𝐓1=(𝐓1𝐓1)23:𝐓˙d𝐓1d𝐓=(𝐓1𝐓1)23(𝐓1)˙=𝐓1𝐓˙𝐓1=(𝐓1𝐓1)24:𝐓˙d𝐓1d𝐓=(𝐓1𝐓1)24

siehe Formelsammlung Tensoralgebra#Spezielle Tensoren vierter Stufe.

Allgemein mit n ∈ ℕ, >0, T0 := 1:

D𝐓n(𝐓)[𝐇]=m=0n1𝐓m𝐇𝐓nm1d𝐓nd𝐓=(m=0n1𝐓m(𝐓nm1))23

#Gâteaux-Differential der Inversen:

𝐓𝐓1=𝟏D𝐓(𝐓)[𝐇]𝐇𝐓1+𝐓D𝐓1(𝐓)[𝐇]=𝟎D𝐓1(𝐓)[𝐇]=𝐓1𝐇𝐓1=(𝐓1𝐓1)23:𝐇D𝐓1(𝐓)[𝐇]=𝐓1𝐇𝐓1=(𝐓1𝐓1)24:𝐇

n ∈ ℕ, >0:

D𝐓n(𝐓)[𝐇]=m=1n0𝐓mD𝐓1(𝐓)[𝐇]𝐓1nm=m=1n0𝐓m1𝐇𝐓nmd𝐓nd𝐓=(m=1n0𝐓m1(𝐓nm))23
D𝐓n(𝐓)[𝐇]=m=1n0(𝐓m1)𝐇(𝐓nm)d𝐓nd𝐓=(m=1n0(𝐓m1)(𝐓nm))24

Orthogonaler Tensor (Q·Q=1):

𝐐˙=𝐐𝐐˙𝐐

Ableitungen nach dem Ort

Nabla-Operator

Vorlage:Siehe auch

#Kartesische Koordinaten x :=e^ixi

#Zylinderkoordinaten: =eρρ+1ρeφφ+ezz

#Kugelkoordinaten: =err+1reϑϑ+1rsin(ϑ)eφφ

#Krummlinige Koordinaten y :=bjyj    mit    bj=yjxie^i.

Gradient

Vorlage:Siehe auch

Definition des Gradienten/Allgemeines

Definierende Eigenschaft bei skalar- oder vektorwertiger Funktion f:[1]

f(y)f(x)=grad(f)(yx)+𝒪(|yx|) wenn yx

Wenn der Gradient existiert, ist er eindeutig. Berechnung bei skalar- oder vektorwertiger Funktion f:

grad(f)h=ddsf(x+sh)|s=0=lims0f(x+sh)f(x)sh𝕍

Integrabilitätsbedingung: Jedes rotationsfreie Vektorfeld ist das Gradientenfeld eines Skalarpotentials:

rot(f)=0g:f=grad(g).

Koordinatenfreie Darstellung als Volumenableitung:

grad(f)=limv0(1vafda)

Skalarfeld f:

grad(f)=f=:fx

Vektorfeld f=fie^i:[2]

grad(f)=(f)=:fx
grad(x)=𝟏

Zusammenhang mit den anderen Differentialoperatoren:

grad(f)=div(f𝟏)=(f𝟏)
grad(f)×c=rot(fc)

Gradient in verschiedenen Koordinatensystemen

#Kartesische Koordinaten:

grad(f)=f,ie^i
grad(f)=f,ie^i=e^igrad(fi)=fi,je^ie^j

#Zylinderkoordinaten:

grad(f)=f,ρe^ρ+f,φρe^φ+f,ze^z
grad(f)=e^ρgrad(fρ)+e^φgrad(fφ)+e^zgrad(fz)+1ρ(fρe^φfφe^ρ)e^φ

#Kugelkoordinaten:

grad(f)=f,re^r+f,ϑre^ϑ+f,φrsin(ϑ)e^φ
grad(f)=e^rgrad(fr)+e^ϑgrad(fϑ)+e^φgrad(fφ)+frr(𝟏e^re^r)e^rfϑe^ϑ+fφe^φr+fϑe^φfφe^ϑrtan(ϑ)e^φ

#Krummlinige Koordinaten: Vorlage:Siehe auch Christoffelsymbole: Γijk=gi,jgk

Vektorfelder:

grad(gi)=Γijkgkgj
grad(gk)=Γijkgigj
grad(figi)=fi|jgigj
grad(figi)=fi|jgigj

Mit den kovarianten Ableitungen

fi|j=f,ji+Γkjifk
fi|j=fi,jΓijkfk

Tensorfelder:

grad(𝐓)[h]=(hgk)𝐓,k=h(gk𝐓,k)=(𝐓,kgk)h

Soll das Argument wie beim Vektorgradient rechts vom Operator stehen, dann lautet der Tensorgradient

grad(𝐓)=𝐓,kgk

Für ein Tensorfeld zweiter Stufe:

grad(Tijgigj)=Tij|kgigjgk,Tij|k=Tij,kΓiklTljΓjklTilgrad(Tijgigj)=Tij|kgigjgk,Tij|k=T,kij+ΓlkiTlj+ΓlkjTilgrad(Ti.jgigj)=Ti.j|kgigjgk,Ti.j|k=Ti,k.jΓiklTl.j+ΓlkjTi.lgrad(T.ji)gigj=T.ji|kgigjgk,T.ji|k=T.j,ki+ΓlkiT.jlΓjklT.li

Produktregel für Gradienten

grad(fg)=(f,ig+fg,i)e^i=grad(f)g+fgrad(g)grad(fg)=(f,ig+fg,i)e^i=ggrad(f)+fgrad(g)grad(fg)=(f,ig+fg,i)e^i=ggrad(f)+fgrad(g)grad(f×g)=(f,i×g+f×g,i)e^i=f×grad(g)g×grad(f)

In drei Dimensionen ist speziell[3]

grad(fg)=grad(f)g+grad(g)f+f×rot(g)+g×rot(f)

Beliebige Basis:

grad(fibi)=bigrad(fi)+figrad(bi)

Divergenz

Vorlage:Siehe auch

Definition der Divergenz/Allgemeines

Vektorfeld f :

div(f)=f=Sp(grad(f))
div(x)=Sp(grad(x))=Sp(𝟏)=3

Klassische Definition für ein Tensorfeld T:[1]

div(𝐓)c=div(𝐓c)c𝕍
div(𝐓)=(𝐓)

Koordinatenfreie Darstellung:

div(f)=limv0(1vafda)

Zusammenhang mit den anderen Differentialoperatoren:

div(f)=f=Sp(grad(f))div(f𝟏)=(f𝟏)=grad(f)

Divergenz in verschiedenen Koordinatensystemen

#Kartesische Koordinaten:

div(f)=f,ie^i=fi,i
div(𝐓)=𝐓,ie^i=Tij,je^i
𝐓=e^i𝐓,i=Tij,ie^j=Tji,je^i

#Zylinderkoordinaten:

div(f)=1ρρ(ρfρ)+1ρfφ,φ+fz,z
div(𝐓)=(Tρρ,ρ+1ρ(Tρφ,φ+TρρTφφ)+Tρz,z)e^ρ+(Tφρ,ρ+1ρ(Tφφ,φ+Tρφ+Tφρ)+Tφz,z)e^φ+(Tzρ,ρ+1ρ(Tzφ,φ+Tzρ)+Tzz,z)e^z

𝐓=div(𝐓) ergibt sich hieraus durch Vertauschen von Tab durch Tba.

#Kugelkoordinaten:

div(f)=fr,r+2fr+fϑ,ϑr+fϑcos(ϑ)+fφ,φrsin(ϑ)div(𝐓)=(Trr,r+2TrrTϑϑTφφ+Trϑ,ϑr+Trφ,φ+Trϑcos(ϑ)rsin(ϑ))e^r(Tϑr,r+2Tϑr+Trϑ+Tϑϑ,ϑr+(TϑϑTφφ)cos(ϑ)+Tϑφ,φrsin(ϑ))e^ϑ(Tφr,r+2Tφr+Trφ+Tφϑ,ϑr+(Tϑφ+Tφϑ)cos(ϑ)+Tφφ,φrsin(ϑ))e^φ

𝐓=div(𝐓) ergibt sich hieraus durch Vertauschen von Tab durch Tba.

Produktregel für Divergenzen

div(fg)=(fg)=(f,ig+fg,i)e^i=grad(f)g+fdiv(g)
div(f×g)=(f×g)=(f,i×g+f×g,i)e^i=grot(f)frot(g)
div(fg)=(f,ig+fg,i)e^i=grad(f)g+div(g)fdiv(f𝐓)=(f,i𝐓+f𝐓,i)e^i=𝐓grad(f)+fdiv(𝐓)div(𝐓f)=(𝐓,if+𝐓f,i)e^i=div(𝐓)f+𝐓:grad(f)div(f×𝐓)=(f,i×𝐓+f×𝐓,i)e^i=i(grad(f)𝐓)+f×div(𝐓)
(fg)=e^i(f,ig+fg,i)=(f)g+(g)f(f𝐓)=e^i(f,i𝐓+f𝐓,i)=(f)𝐓+f𝐓(𝐓f)=e^i(𝐓,if+𝐓f,i)=(𝐓)f+𝐓:(f)(𝐓×f)=e^i(𝐓,i×f+𝐓×f,i)=(𝐓)×fi((f)𝐓)

Beliebige Basis:

div(fibi)=(fibi)=grad(fi)bi+fidiv(bi)
div(Tijbibj)=(grad(Tij)bj)bi+Tij(grad(bi)bj+div(bj)bi)
(Tijbibj)=((Tij)bi)bj+Tij((bi)bj+(bj)bi)

Produkt mit Konstanten:

div(f𝐂)=𝐂grad(f)div(f𝟏)=grad(f)
(f𝐂)=grad(f)𝐂(f𝟏)=f
div(𝐂f)=𝐂:grad(f)div(f)=div(𝟏f)=𝟏:grad(f)=Sp(grad(f))

Rotation

Vorlage:Siehe auch

Definition der Rotation/Allgemeines

Vektorfeld f :

rot(f)=×f

Klassische Definition für ein Tensorfeld T:[1]

rot(𝐓)c=rot(𝐓c)c𝕍
rot(𝐓)=×(𝐓)

Allgemeine Identitäten:

𝐓=𝐓Sp(rot(𝐓))=Sp(×𝐓)=0
rot(x)=0

Integrabilitätsbedingung[4]: Jedes divergenzfreie Vektorfeld ist die Rotation eines Vektorfeldes:

div(f)=0g:f=rot(g).

Siehe auch #Satz über rotationsfreie Felder.

Koordinatenfreie Darstellung:

rot(f)=limv0(1vaf×da)

Zusammenhang mit den anderen Differentialoperatoren:

rot(fc)=grad(f)×crot(f)=i(grad(f))=i(f)=×f

Rotation in verschiedenen Koordinatensystemen

#Kartesische Koordinaten:

rot(f)=fj,ie^i×e^j=ϵijkfj,ie^k=(f3,2f2,3)e^1+(f1,3f3,1)e^2+(f2,1f1,2)e^3
rot(𝐓)=e^i×𝐓,i=e^i×Tlj,ie^je^l=ϵijkTlj,ie^ke^l

#Zylinderkoordinaten:

rot(f)=fz,φρfφ,zρe^ρ+(fρ,zfz,ρ)e^φ+fφ+ρfφ,ρfρ,φρe^z
rot(𝐓)=e^ρ×(𝐓,ρ)+1ρe^φ×(𝐓,φ)+e^z×(𝐓,z)
×𝐓=e^ρ×𝐓,ρ+1ρe^φ×𝐓,φ+e^z×𝐓,z

#Kugelkoordinaten:

rot(f)=fφ,ϑsin(ϑ)+fφcos(ϑ)fϑ,φrsin(ϑ)e^r+(fr,φrsin(ϑ)fφ+rfφ,rr)e^ϑ+fϑ+rfϑ,rfr,ϑre^φ
rot(𝐓)=e^r×(𝐓,r)+1re^ϑ×(𝐓,ϑ)+1rsin(ϑ)e^φ×(𝐓,φ)
×𝐓=e^r×𝐓,r+1re^ϑ×𝐓,ϑ+1rsin(ϑ)e^φ×𝐓,φ

Produktregel für Rotationen

rot(fg)=e^i×(f,ig+fg,i)=grad(f)×g+frot(g)rot(f×g)=e^i×(f,i×g+f×g,i)=(e^ig)f,i(e^if,i)g+(e^ig,i)f(e^if)g,i=grad(f)gdiv(f)g+div(g)fgrad(g)f=div(fg)div(gf)=(gf)(fg)
rot(fg)=e^i×(g,if+gf,i)=rot(g)fg×grad(f)rot(f𝐓)=e^k×(f,k𝐓+f𝐓,k)=grad(f)×(𝐓)+frot(𝐓)

rot(𝐓f)=e^k×(𝐓,kf+𝐓f,k)=rot(𝐓)f+i(e^k𝐓f,k)=rot(𝐓)fi(𝐓grad(f))rot(f×𝐓)=rot((𝐓×f))=×(𝐓×f)=(×𝐓)×f+𝐓#(f)=rot(𝐓)×f+(𝐓#grad(f))

×(fg)=e^i×(f,ig+fg,i)=(×f)gf×((g)×(f𝐓)=e^k×(f,k𝐓+f𝐓,k)=(f)×𝐓+f×𝐓

×(𝐓f)=e^k×(𝐓,kf+𝐓f,k)=(×𝐓)f+i(e^k𝐓f,k)=(×𝐓)fi(𝐓(f))×(𝐓×f)=e^k×(𝐓,k×f+(𝐓e^i)e^i×f,k)=(×𝐓)×f(𝐓e^i)×e^ke^i×f,k=(×𝐓)×f𝐓#(f)

Beliebige Basis:

rot(fibi)=grad(fi)×bi+firot(bi)

Produkt mit Konstanten:

rot(𝐂f)=i(𝐂grad(f))rot(f)=rot(𝟏f)=i(grad(f))rot(f×𝟏)=𝟏#grad(f)=grad(f)div(f)𝟏

In divergenzfreien Feldern ist also: rot(f×𝟏)=grad(f)

Laplace-Operator

Vorlage:Siehe auch

Definition/Allgemeines

Δ:==2

Zusammenhang mit anderen Differentialoperatoren:

Δf=div(grad(f))=(f)Δf=div(grad(f))=(f)

„Vektorieller Laplace-Operator“:

Δf=grad(div(f))rot(rot(f))

Laplace-Operator in verschiedenen Koordinatensystemen

#Kartesische Koordinaten:

Δf=f,kkΔf=Δfie^i=fi,kke^iΔ𝐓=ΔTije^ie^j=Tij,kke^ie^j

#Zylinderkoordinaten:

Δf=f,ρρ+f,ρρ+f,φφρ2+f,zzΔf=(Δfρ2fφ,φ+fρρ2)e^ρ+(Δfφ+2fρ,φfφρ2)e^φ+Δfze^z

#Kugelkoordinaten:

Δf=1r2r(r2fr)+1r2sin(ϑ)ϑ(sin(ϑ)fϑ)+1r2sin2(ϑ)2fφ2=2f,rr+f,rr+f,ϑcos(ϑ)+f,ϑϑsin(ϑ)r2sin(ϑ)+f,φφr2sin2(ϑ)Δf=(Δfr2r2(fr+fϑ,ϑ)2fφ,φ+fϑcos(ϑ)r2sin(ϑ))e^r+(Δfϑ+2fr,ϑr2fϑ+2fφ,φcos(ϑ)r2sin2(ϑ))e^ϑ+(Δfφfφ2fϑ,φcos(ϑ)2fr,φsin(ϑ)r2sin2(ϑ))e^φ

Verknüpfungen

Wegen der in der Literatur teilweise abweichenden Definitionen der Differentialoperatoren kann es in der Literatur zu abweichenden Formeln kommen. Wenn die Definitionen der Literatur hier eingesetzt werden, gehen die hiesigen Formeln in die der Literatur über.

div(rot(f))=(×f)=0rot(grad(f))=×f=0div(grad(f)×grad(g))=(f×g)=g(×f)=0rot(grad(f))=×(f)=𝟎div(rot(𝐓))=(×𝐓)=0
div(grad(f))=(f)=()f=Δfdiv(grad(f))=(f)=()f=Δf
div(grad(f))=(f)=fi,ije^j=grad(div(f))rot(grad(f))=×((f))=×(f,ie^i)=grad(rot(f))
rot(rot(f))=×(×f)=(f)Δf=grad(div(f))Δfrot(rot(𝐓))=(×(×(𝐓)))=(𝐓)()𝐓=grad(div(𝐓))Δ𝐓
rot(rot(𝐓))=Δ𝐓grad(grad(Sp(𝐓)))+grad(div(𝐓))+grad(div(𝐓))+[ΔSp(𝐓)div(div(𝐓))]𝟏

Bei symmetrischem T = T gilt:

rot(rot(𝐓))=Δ𝐓grad(grad(Sp(𝐓)))+grad(div(𝐓))+grad(div(𝐓))+[ΔSp(𝐓)div(div(𝐓))]𝟏


Wenn zusätzlich 𝐓=𝐓=𝐆Sp(𝐆)𝟏 dann ist:

rot(rot(𝐓))=Δ𝐆+grad(div(𝐆))+grad(div(𝐆))div(div(𝐆))𝟏

Der Laplace-Operator kann zwischen den anderen Operatoren wie ein Skalar behandelt werden, also an beliebiger Stelle in die Formeln eingesetzt werden, z. B.:

Δrot(rot(f))=rot(Δrot(f))=rot(rot(Δf))==Δgrad(div(f))ΔΔf=grad(Δdiv(f))ΔΔf=grad(div(Δf))ΔΔf

Grassmann-Entwicklung

f×rot(f)=12grad(ff)grad(f)f=(grad(f)grad(f))f=i(grad(f))×f
grad(f)f=12grad(ff)f×rot(f)

Sätze über Gradient, Divergenz und Rotation

Ein Vektorfeld, dessen Divergenz und Rotation verschwindet, ist harmonisch:

div(f)=0undrot(f)=0Δf=0

Helmholtz-Theorem

Vorlage:Siehe auch

Jedes Vektorfeld lässt sich eindeutig in einen divergenzfreien und einen rotationsfreien Anteil zerlegen. Den Integrabilitätsbedingungen für Rotationen und Gradienten zufolge ist der erste Anteil ein Rotationsfeld und der zweite ein Gradientenfeld.
f=f1+f2:div(f1)=0undrot(f2)=0g,g:f=rot(g)+grad(g)

Satz über rotationsfreie Felder

Vorlage:Siehe auch

I:rot(u):=e^k×u,k=0f:u=grad(f)II:rot(𝐓)=𝟎u:𝐓=grad(u)III:rot(𝐓)=𝟎undSp(𝐓)=0𝐖:𝐓=rot(𝐖)und𝐖=𝐖

oder

II:×(𝐓)=𝟎u:𝐓=grad(u)III:×(𝐓)=𝟎undSp(𝐓)=0𝐖:𝐓=rot(𝐖)und𝐖=𝐖

Gaußscher Integralsatz

Vorlage:Siehe auch

vgrad(f)dv=afdavgrad(f)dv=afdavdiv(f)dv=afdavrot(f)dv=af×davdiv(𝐓)dv=a𝐓dav𝐓dv=a𝐓da

Mit der #Produktregel für Gradienten, #Produktregel für Divergenzen und #Produktregel für Rotationen können Formeln für die partielle Integration im Mehrdimensionalen abgeleitet werden, beispielsweise:

grad(fg)=grad(f)g+fgrad(g)vgrad(f)gdv=afgdavfgrad(g)dv

Klassischer Integralsatz von Stokes

Vorlage:Siehe auch

Gegeben:

Vektorwertige Funktion f(x,t) :

arot(f)da=bfdb

Mit der #Produktregel für Rotationen können Formeln für die partielle Integration im Mehrdimensionalen abgeleitet werden, beispielsweise:

rot(fg)=grad(f)×g+frot(g)a(grad(f)×g)da=bfgdbafrot(g)da

Reynoldscher Transportsatz

Vorlage:Siehe auch Gegeben:

  • Zeit t
  • Zeitabhängiges Volumen v mit Volumenform dv mit
  • Oberfläche des Volumes a und äußerem vektoriellem Oberflächenelement da
  • Ortsvektoren xv
  • Geschwindigkeitsfeld:v(x,t)
  • Eine skalare oder vektorwertige Dichtefunktion pro Volumeneinheit f(x,t), die mit den sich bewegenden Partikeln transportiert wird.
  • Die Integrale Größe für das Volumen:vf(x,t)dv

Skalare Funktion f(x,t) :

ddtvfdv=vftdv+af(vda)=v(ft+div(fv))dv=v(ft+grad(f)v+div(v)f)dv=v(f˙+div(v)f)dv

Vektorwertige Funktion f(x,t) :

ddtvfdv=vftdv+af(vda)=v(ft+div(vf))dv=v(ft+grad(f)v+div(v)f)dv=v(f˙+div(v)f)dv

Transportsatz für Flächenintegrale

Gegeben:

  • Zeit t
  • Ortsvektoren xv
  • Geschwindigkeitsfeld:v(x,t)
  • Zeitabhängige Fläche a:[0,1]2v, die mit dem Geschwindigkeitsfeld transportiert wird und auf der mit räumlichem, vektoriellem Oberflächenelement da im Volumen v integriert wird
  • Eine skalare oder vektorwertige Feldgröße f(x,t), die mit den sich bewegenden Partikeln transportiert wird.
  • Die Integrale Größe auf der Fläche:af(x,t)da

Skalare Funktion f(x,t) :

ddtafda=a[f˙𝟏+fdiv(v)𝟏fgrad(v)]da

Vektorwertige Funktion f(x,t):

ddtafda=a[f˙+fdiv(v)grad(v)f]da

Transportsatz für Kurvenintegrale

Gegeben:

  • Zeit t
  • Ortsvektoren xv
  • Geschwindigkeitsfeld:v(x,t)
  • Zeitabhängige Kurve b:[0,1)v, die mit dem Geschwindigkeitsfeld transportiert wird und entlang derer mit räumlichem, vektoriellem Linienelement db im Volumen v integriert wird
  • Eine skalare oder vektorwertige Feldgröße f(x,t), die mit den sich bewegenden Partikeln transportiert wird.
  • Die Integrale Größe entlang des Weges:bf(x,t)db

Skalare Funktion f(x,t) :

ddtbfdb=b(f˙𝟏+fgradv)db

Vektorwertige Funktion f(x,t):

ddtbfdb=b(f˙+fgradv)db

Kontinuumsmechanik

Kleine Deformationen

Ingenieursdehnungen:

ε=εije^ie^j=12(ui,j+uj,i)e^ie^j

Kompatibilitätsbedingungen:

rot(rot(ε))=×(×ε)=𝟎2ε12,12ε22,11ε11,22=02ε13,13ε33,11ε11,33=02ε23,23ε33,22ε22,33=0ε11,23+ε23,11ε12,13ε13,12=0ε22,13+ε13,22ε12,23ε23,12=0ε12,33+ε33,12ε13,23ε23,13=0

Starrkörperbewegung

Orthogonaler Tensor 𝐐 beschreibt die Drehung.

Ω:=𝐐˙𝐐=(𝐐𝐐˙)=𝐐𝐐˙

Vektorinvariante oder dualer axialer Vektor ω des schiefsymmetrischen Tensors Ω ist die Winkelgeschwindigkeit:

Ωr=ω×rr

Starrkörperbewegung mit r=const. :

x=f+𝐐rr=𝐐(xf)
v=f˙+𝐐˙r=f˙+𝐐˙𝐐(xf)=f˙+Ω(xf)=f˙+ω×(xf)

Ableitungen der Invarianten

I1(𝐓)𝐓=Sp(𝐓)𝐓=𝟏
I2(𝐓)𝐓=I1(𝐓)𝟏𝐓
I3(𝐓)𝐓=det(𝐓)𝐓=det(𝐓)𝐓1=cof(𝐓)=𝐓𝐓I1(𝐓)𝐓+I2(𝐓)𝟏

mit der transponiert inversen T⊤-1 und dem Kofaktor cof(T) des Tensors T.

Funktion f der Invarianten:

f𝐓(I1(𝐓),I2(𝐓),I3(𝐓))=(fI1+I1fI2+I2fI3)𝟏(fI2+I1fI3)𝐓+fI3𝐓𝐓

Ableitung der Frobenius-Norm:

𝐓𝐓=𝐓𝐓

Eigenwerte (aus der impliziten Ableitung des charakteristischen Polynoms):

𝐓v=λvdet(𝐓λ𝟏)=λ3+I1λ2I2λ+I3=0
dλd𝐓=(λ2λI1+I2)𝟏+(λI1)𝐓+𝐓𝐓3λ22I1λ+I2

Eigenwerte symmetrischer Tensoren:

𝐓v=λvλ𝐓=vv

Eigenwerte von 𝐓=i=13λivivi, wo vi dual zu den Eigenvektoren vi sind (vivj=δij):

λi𝐓=vivi (keine Summe)

Die Eigenwerte von 𝐓=cv1v1+a(v2v2+v3v3)+b(v2v3v3v2) sind λ1=c,λ2=a+ib,λ3=aib mit den Eigenvektoren v1,w2=v2+iv3,w3=v2iv3. Hier ist:

λ1𝐓=v1v1,λk𝐓=12wkwk,k=2,3 (keine Summe)

mit w2=v2+iv3,w3=v2iv3 und der Überstrich markiert den konjugiert komplexen Wert.

Konvektive Koordinaten

Vorlage:Siehe auch

Konvektive Koordinaten y1,y2,y3

Kovariante Basisvektoren Bi=dXdyi,    bi=dxdyi

Kontravariante Basisvektoren Bi=dyidX:=GRAD(yi),    bi=dyidx:=grad(yi)

BiBj=bibj=δij

Deformationsgradient 𝐅=biBi

Räumlicher Geschwindigkeitsgradient 𝐥=b˙ibi=bib˙i

Kovarianter Tensor 𝐓=Tijbibj

Kontravarianter Tensor 𝐓=Tijbibj

Geschwindigkeitsgradient

Vorlage:Siehe auch Räumlicher Geschwindigkeitsgradient:𝐥=grad(v)=𝐅˙𝐅1

Divergenz der Geschwindigkeit:div(v)=Sp(𝐥)

Winkelgeschwindigkeit oder Wirbelstärke ist der duale axiale Vektor

ω=𝐥A=12i(𝐥)=12rot(v)
DDtdet(𝐅)=det(𝐅)𝐅1:𝐅˙=det(𝐅)Sp(𝐅˙𝐅1)=det(𝐅)div(v)

Objektive Zeitableitungen

Vorlage:Siehe auch

Bezeichnungen wie in #Konvektive Koordinaten.

Räumlicher Geschwindigkeitsgradient 𝐥=b˙ibi=bib˙i=𝐝+𝐰

Räumliche Verzerrungsgeschwindigkeit 𝐝=12(𝐥+𝐥)

Wirbel- oder Spintensor 𝐰=12(𝐥𝐥)

Objektive Zeitableitungen von Vektoren

Gegeben:v=vibi=vibi:

vΔ=v˙+𝐥v=v˙ibiv=v˙𝐥v=v˙ibiv=v˙𝐰v

Objektive Zeitableitungen von Tensoren

Gegeben:𝐓=Tijbibj=Tijbibj

𝐓Δ=𝐓˙+𝐓𝐥+𝐥𝐓=T˙ijbibj𝐓=𝐓˙𝐥𝐓𝐓𝐥=T˙ijbibj𝐓=𝐓˙+𝐓𝐰𝐰𝐓𝐓=𝐓˙+Sp(𝐥)𝐓𝐥𝐓𝐓𝐥

Materielle Zeitableitung

Vorlage:Siehe auch

f˙(x,t)=DfDt=ft+grad(f)v=ft+(v)f
f˙(x,t)=DfDt=ft+grad(f)v=ft+(v)f

#Kartesische Koordinaten:DfDt:=ft+vxfx+vyfy+vzfz

#Zylinderkoordinaten:DfDt:=ft+vρfρ+vφρfφ+vzfz

#Kugelkoordinaten:DfDt:=ft+vrfr+vφrsin(ϑ)fφ+vϑrfϑ

Materielle Zeitableitungen von Vektoren werden mittels DfDt=DfiDte^i daraus zusammengesetzt.

Fußnoten

  1. 1,0 1,1 1,2 Morton E. Gurtin: „The linear theory of elasticity.“ In: S. Flügge (Hrsg.): Handbuch der Physik. Band VIa/2.: Festkörpermechanik II / C. Truesdell (Bandherausgeber). Springer, Berlin 1972, ISBN 3-540-05535-5, S. 10 ff.
  2. In der Literatur (z. B. Altenbach 2012) wird auch die transponierte Beziehung benutzt:
    grad~(f)=f=e^ifxi=fjxie^ie^j=grad(f)
    Dann muss, um die Formeln zu vergleichen, grad~(f) und grad(f) vertauscht werden.
  3. Vorlage:Literatur
  4. R. Greve (2003), S. 111.

Literatur