Formelsammlung Tensoranalysis

Formelsammlung Tensoranalysis


xn Dieser Artikel ist eine Formelsammlung zum Thema Tensoranalysis. Es werden mathematische Symbole verwendet, die im Artikel Liste mathematischer Symbole erläutert werden.

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

Allgemeines

Siehe auch

Formelsammlung Tensoralgebra

Nomenklatur

  • Operatoren wie „grad“ werden nicht kursiv geschrieben.
  • Buchstaben in der Mitte des Alphabets werden als Indizes benutzt: i,j,k,l{1,2,3}
  • Es gilt die Einsteinsche Summenkonvention ohne Beachtung der Indexstellung.
    • Kommt in einer Formel in einem Produkt ein Index doppelt vor wie in c=aibi wird über diesen Index von eins bis drei summiert:
      c=aibi=i=13aibi.
    • Kommen mehrere Indizes doppelt vor wie in c=AijBji wird über diese summiert:
      c=AijBji=i=13j=13AijBji.
    • Ein Index, der nur einfach vorkommt wie i in vi=Aijbj, ist ein freier Index. Die Formel gilt dann für alle Werte der freien Indizes:
      vi=Aijbjvi=j=13Aijbji{1,2,3}.
  • Vektoren:
    • Alle hier verwendeten Vektoren sind geometrische Vektoren im dreidimensionalen euklidischen Vektorraum 𝕍.
    • Vektoren werden mit Kleinbuchstaben bezeichnet.
    • Einheitsvektoren mit Länge eins werden wie in e^ mit einem Hut versehen.
    • Vektoren mit unbestimmter Länge werden wie in a mit einem Pfeil versehen.
    • Standardbasis e^1,e^2,e^3
    • Beliebige Basis b1,b2,b3 mit dualer Basis b1,b2,b3
    • Der Vektor x=xie^i wird durchgängig Ortsvektor genannt.
  • Tensoren zweiter Stufe werden wie in 𝐓 mit fetten Großbuchstaben notiert.
  • Koordinaten:
  • Konstanten: c,c,𝐂
  • Zeit t
  • Variablen: skalar r,s oder vektorwertig r,s𝕍3
  • Funktionen:
    • Skalar f,g oder vektorwertig f,g𝕍3
    • Tensorwertig: 𝐓(x,t) oder 𝐓(y,t)
  • Operatoren:
  • Differentialoperatoren:
  • Kontinuumsmechanik:

Kronecker-Delta

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

Basisvektoren

Kartesische Koordinaten: e^1,e^2,e^3

Zylinderkoordinaten:

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

Kugelkoordinaten:

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

Krummlinige Koordinaten: y1,y2,y3

bi=xyi,bi=grad(yi)=yixbibj=δij

Ableitungen nach dem Ort

Gâteaux-Differential

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:

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

Fréchet-Ableitung

Existiert ein beschränkter linearer Operator 𝒜, so dass

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

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

fx=𝒜.

Nabla Operator

Kartesische Koordinaten x : =e^ixi

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

Gradient

Skalarfeld f :

grad(f)=f=fxie^i=f,ie^i

Vektorfeld f=fie^i:[1]

grad(f)=(f)=fxje^j=fixje^ie^j=fi,je^ie^j
grad(x)=𝐈

Skalar- oder vektorwertige Funktion f:

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

Zylinderkoordinaten:

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

Kugelkoordinaten:

grad(f)=fre^r+1rfϑe^ϑ+1rsinϑfφe^φ
grad(f)=e^rgrad(fr)+frr𝐈frre^re^r+e^ϑgrad(fϑ)+fϑrtanϑe^φe^φfϑre^re^ϑ+e^φgrad(fφ)fφrtanϑe^ϑe^φfφre^re^φ

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

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

Koordinatenfreie Darstellung als Volumenableitung:

grad(f)=limv0(1vafda)

Zusammenhang mit den anderen Differentialoperatoren:

grad(f)=div(f𝐈)

Produktregel:

grad(fg)=e^if,ig+fe^ig,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=grad(f)g+grad(g)fgrad(f×g)=(f,i×g+f×g,i)e^i=g×grad(f)+f×grad(g)

Beliebige Basis:

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

Produkt mit Konstanten:

grad(cf)=cgrad(f)
grad(fc)=cgrad(f)
grad(cf)=cgrad(f)
grad(cf)=grad(f)c
grad(c×f)=c×grad(f)=grad(f×c)
grad(𝐂f)=𝐂grad(f)

Divergenz

Vektorfeld f=fie^i :

div(f)=f=fxie^i=fixi=Sp(grad(f))
divx=3

Tensorfeld 𝐓=Tije^ie^j:[2]

div(𝐓)=𝐓=e^i𝐓xi=Tijxie^j

Zylinderkoordinaten:

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

Kugelkoordinaten:

div(f)=1r2r(r2fr)+1rsinϑϑ(fϑsinϑ)+1rsinϑfφφdiv(𝐓)=(1r2r(r2Trr)+1rsinϑTφr,φ+1rsinϑϑ(Tϑrsinϑ)1r(Tϑϑ+Tφφ))e^r+(1r3r(r3Trϑ)+1rsinϑTφϑ,φ+1rsinϑϑ(Tϑϑsinϑ)+1r(TϑrTrϑTφφcotϑ))e^ϑ+(1r3r(r3Trφ)+1rsinϑTφφ,φ+1rsinϑϑ(Tϑφsinϑ)+1r(TφrTrφ+Tφϑcotϑ))e^φ

Koordinatenfreie Darstellung:

div(f)=limv0(1vafda).

Zusammenhang mit den anderen Differentialoperatoren:

div(f)=Sp(grad(f))div(f𝐈)=grad(f)div(fie^i)=grad(fi)e^idiv(f×e^i)=rot(f)e^i

Produktregel:

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

Beliebige Basis:

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

Produkt mit Konstanten:

div(cf)=cdiv(f)
div(fc)=grad(f)c
div(c𝐓)=cdiv(𝐓)
div(f𝐂)=grad(f)𝐂div(f𝐈)=grad(f)
div(𝐂f)=𝐂:grad(f)div(𝐈f)=div(f)=𝐈:grad(f)=Sp(grad(f))
div(𝐓c)=div(𝐓)c
div(cf)=grad(f)c
div(fc)=div(f)c

Rotation

Vektorfeld f=fie^i :

rot(f)=×f=fjxie^i×e^j=εijkfjxie^k=(f3x2f2x3)e^1+(f1x3f3x1)e^2+(f2x1f1x2)e^3
rotx=0

Tensorfeld 𝐓=Tije^ie^j:[3]

rot(𝐓)=×𝐓=e^i×𝐓xi=e^i×Tjl,ie^je^l=εijkTjl,ie^ke^l
𝐓=𝐓Sp(rot(𝐓))=0

Zylinderkoordinaten:

rotf=[1ρfzφfφz]e^ρ+[fρzfzρ]e^φ+1ρ[ρ(ρfφ)fρφ]e^z

Kugelkoordinaten:

rot(f)=1rsinϑ[ϑ(fφsinϑ)fϑφ]e^r+[1rsinϑfrφ1rr(rfφ)]e^ϑ+1r[r(rfϑ)frϑ]e^φ

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

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

Koordinatenfreie Darstellung:

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

Zusammenhang mit den anderen Differentialoperatoren:

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

Produktregel:

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(gf)div(fg)rot(fg)=e^i×(f,ig+fg,i)=rot(f)gf×grad(g)rot(𝐓f)=e^k×(𝐓,kf+(𝐓e^i)e^if,k)=rot(𝐓)f(e^if,k)(𝐓e^i)×e^k=rot(𝐓)f𝐓×grad(f)rot(f𝐓)=e^k×(f,k𝐓+f𝐓,k)=grad(f)×𝐓+frot(𝐓)rot(𝐓×f)=e^k×(𝐓,k×f+(𝐓e^i)e^i×f,k)=rot(𝐓)×fe^k×(𝐓e^i)f,k×e^i=rot(𝐓)×fgrad(f)××𝐓rot(f×𝐂)=e^k×(f,k×(𝐂e^i)e^i)=(e^k𝐂e^i)f,ke^i(e^kf,k)𝐂e^ie^i=grad(f)𝐂div(f)𝐂

Beliebige Basis:

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

Produkt mit Konstanten:

rot(cf)=crot(f)rot(fc)=grad(f)×crot(c×f)=rot(f×c)=div(f)cgrad(f)c=div(fc)div(cf)rot(cf)=c×grad(f)rot(fc)=rot(f)crot(𝐂f)=𝐂×grad(f)rot(f)=rot(𝐈f)=𝐈×grad(f)rot(𝐓c)=rot(𝐓)crot(c𝐓)=crot(𝐓)rot(f𝐂)=grad(f)×𝐂rot(𝐂×f)=grad(f)××𝐂rot(𝐓×c)=rot(𝐓)×crot(f×𝐈)=grad(f)div(f)𝐈

Satz über rotationsfreie Felder

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

Laplace-Operator

Δ:==2

Kartesische Koordinaten:

Δf=f2x12+f2x22+f2x32=f,kkΔf=2f=()f=Δfie^i=fi,kke^iΔ𝐓=2𝐓=()𝐓=ΔTije^ie^j=Tij,kke^ie^j

Zylinderkoordinaten:

Δf=1ρρ(ρfρ)+1ρ22fϕ2+2fz2Δf=(Δfρ1ρ2fρ2ρ2fφ,φ)e^ρ+(Δfφ1ρ2fφ+2ρ2fρ,φ)e^φ+Δfze^z

Kugelkoordinaten:

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

Zusammenhang mit anderen Differentialoperatoren:

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

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)=0div(rot(𝐓))=×𝐓=0div(rot(rot(𝐓)))=(×(×𝐓))=0rot(grad(f))=×f=0rot(grad(f))=×(f)=𝟎div(grad(f)×grad(g))=(f×g)=g(×f)=0div(grad(f))=(f)=()f=Δfdiv(grad(f))=(f)=()f=Δfdiv(grad(f))=(f)=(f)=grad(div(f))rot(grad(f))=×(f)=(×f)=grad(rot(f))rot(rot(f))=×(×f)=(f)Δf=grad(div(f))Δfrot(rot(𝐓))=×[×(𝐓e^i)e^i]=[div(𝐓)e^i]e^iΔ𝐓=(div(𝐓))e^ie^iΔ𝐓=grad(div(𝐓))Δ𝐓
rot(rot(𝐓))=Δ𝐓grad(grad(Sp(𝐓)))+grad(div(𝐓))+grad(div(𝐓))+[ΔSp(𝐓)div(div(𝐓))]𝐈

Bei symmetrischem 𝐓=𝐆Sp(𝐆)𝐈 gilt außerdem:

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

Der #Laplace-Operator kann 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
grad(f)f=12grad(ff)f×rot(f)

Sätze über Gradient, Divergenz und Rotation

  • Helmholtz-Theorem: 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)=0rot(f2)=0g,g:f=rot(g)+grad(g)
  • Ein Vektorfeld, dessen Divergenz und Rotation verschwindet, ist harmonisch:
    div(f)=0rot(f)=0Δf=0.
  • Siehe auch Satz über rotationsfreie Felder

Integralsätze

Gaußscher Integralsatz

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

Klassischer Integralsatz von Stokes

Gegeben:

Vektorwertige Funktion f(x,t) :

arot(f)da=bfdb

Kontinuumsmechanik

Kleine Deformationen

Ingenieursdehnungen:

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

Kompatibilitätsbedingungen:

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 axialer Vektor ω des schiefsymmetrischen Tensors Ω :

Ω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(𝐓)𝐓T1

mit der transponiert inversen 𝐓T1 des Tensors 𝐓.

𝐓𝐓=𝐓𝐓

Eigenwerte (keine Summe über i):

𝐓vi=λiviλi𝐓=vivi

Funktion f der Invarianten:

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

Zeitableitungen der Invarianten

DI1(𝐓)Dt=Sp(𝐓˙)
DI2(𝐓)Dt=Sp(𝐓)Sp(𝐓˙)Sp(𝐓𝐓˙)
DI3(𝐓)Dt=Ddet(𝐓)Dt=det(𝐓)Sp(𝐓˙𝐓1)
D𝐓Dt=𝐓𝐓˙𝐓

Zeitableitung von inversen Tensoren

DDt(𝐓1)=𝐓1𝐓˙𝐓1
DDt(𝐓1)=𝐓1𝐓˙𝐓1

Orthogonale Tensoren:

𝐐˙=𝐐𝐐˙𝐐

Konvektive Koordinaten

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

Räumlicher Geschwindigkeitsgradient: 𝐥=grad(v)

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

DDtdet(𝐅)=det(𝐅)div(v)
𝐥=𝐅˙𝐅1=𝐅DDt(𝐅1)

Objektive Zeitableitungen

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

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

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.

Transportsätze

Reynoldscher Transportsatz

Gegeben:

  • Zeit t
  • Zeitabhängiges Volumen v mit Volumenform dv mit
  • Die 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 Punkten 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 Kurvenintegrale

Gegeben:

  • Zeit t
  • Zeitabhängige Kurve b:[0,1)v, entlang derer mit räumlichem, vektoriellem Linienelement db im Volumen v integriert wird
  • Ortsvektoren xv
  • Geschwindigkeitsfeld: v(x,t)
  • Eine skalare oder vektorwertige Feldgröße f(x,t), die mit den sich bewegenden Punkten 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

Fußnoten

  1. 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.
  2. In der Literatur (z. B. Truesdell 1972) wird auch die transponierte Beziehung benutzt:
    div~(𝐓)=𝐓=𝐓xke^k=Tik,ke^i=div(𝐓)
    Dann muss, um die Formeln zu vergleichen, div~(𝐓) und div(𝐓) vertauscht werden.
  3. In der Literatur (z. B. Truesdell 1972) wird auch die transponierte Beziehung benutzt:
    rot~(𝐓)=×𝐓=e^k×𝐓xk=e^k×Tij,ke^je^i=rot(𝐓)
    Dann muss, um die Formeln zu vergleichen, rot~(𝐓) und rot(𝐓) vertauscht werden.
  4. R. Greve (2003), S. 111

Literatur

  • H. Altenbach: Kontinuumsmechanik. Springer, 2012, ISBN 978-3-642-24118-5.
  • M. Bestehorn: Hydrodynamik und Strukturbildung. Springer, 2006, ISBN 978-3-540-33796-6.
  • Adolf J. Schwab: Begriffswelt der Feldtheorie. Praxisnahe, anschauliche Einführung. Elektromagnetische Felder, Maxwellsche Gleichungen, Gradient, Rotation, Divergenz. 6., unveränderte Auflage. Springer, Berlin u. a. 2002, ISBN 3-540-42018-5.
  • Konrad Königsberger: Analysis. überarbeitete Auflage. Band 2. 4. Springer, Berlin u. a. 2000, ISBN 3-540-43580-8.
  • Ralf Greve: Kontinuumsmechanik. Springer, 2003, ISBN 3-540-00760-1.
  • C. Truesdell: Festkörpermechanik II. In: S. Flügge (Hrsg.): Handbuch der Physik. Band VIa/2. Springer, 1972, ISBN 3-540-05535-5.