Einschaliges Hyperboloid

Prof. Dr. Dörte Haftendorn, MuPAD 4,    Juni 08 Update 29. Juni 08

Web:  https://mathe.web.leuphana.de             www.mathematik-verstehen.de

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Eine Gerade, die windschief zur z-Achse ist, dreht sich um die z-Achse.

Erst ist hier das Ergebnis gezeichnet.

f:=4:c:=2:

hy:=plot::Implicit3d(x^2/f^2+y^2/f^2-z^2/(f^2*c^2)=1,

      x=-10..10,y=-10..10,z=-10..10, Mesh=[20,20,20]):

ger:=plot::Curve3d([f,r,r*c],r=-5..5,

     LineWidth=1,LineColor=[0,1,0]);

math

Rotations-Hyperboloid mit erzeugender Gerade, letztere in Parameterdarstellung

plot(hy,ger)

MuPAD graphics

Nun soll die Gerade sich drehen um die z-Achse

A:=matrix([[cos(t),-sin(t),0],[sin(t),cos(t),0],[0,0,1]]);

math

gera:=plot::Curve3d([f*cos(t)-r*sin(t),f*sin(t)+r*cos(t),r*c],r=-5..5, t=0..2*PI, LineWidth=1, LineColor=[0,1,0]):

plot(hy,gera)

MuPAD graphics

image

Natürlich ist der Einheitsvektor der z-Achse der Eigenvektor der Matrix A

zum Eigenwert 1.

Ev:= linalg::eigenvectors(A)

math

Volumenbrechnung

Für y=0 wird nach x^2 aufgelöst, der dann rechts stehende Term ist der, der über z zu integrieren ist.

hold(x^2/f^2+y^2/f^2-z^2/(f^2*c^2)=1);

hold(x^2=f^2+z^2/c^2)

math

math

2*PI*int(f^2+z^2/c^2, z);

2*PI*int(f^2+z^2/c^2, z=0..10);

Vhyp:=float(%/PI)*PI

math

math

math

Berechnung des oberen Radius bei Höhe 10 von der Mitte aus.

hyp:=x^2/f^2+y^2/f^2-z^2/(f^2*c^2)=1

math

solve((hyp|z=10)|y=0,x);

solve((hyp|z=-10)|y=0,x);

math

math

Das durch Rotationsbilder erzeugte Hyperboloid passt genau.

ob:=plot::ZRotate(c*sqrt(x^2-f^2),x=0..sqrt(41)):

unt:=plot::ZRotate(-c*sqrt(x^2-f^2),x=0..sqrt(41)):

plot(ob,unt);

plot(ob,unt, hy)

MuPAD graphics

MuPAD graphics

Volumen eines umfassenden Zylinders

oder eines innen liegenden Kegels

Vzyl:=PI*41*20;

Vkeg:=2*PI*41*10/3;

2.0*PI*41*10/3;

math

math

math

Berechnung der beiden Kegelstümpfe

kleinerer Radius

solve((hyp|z=0)|y=0,x)

math

Kegelstumpf-Volumenformel

hold(Vst=PI/3*(R^2+r*R+r^2)*h)

math

2 solche Stümpfe

2*float(1/3*(41+4*sqrt(41)+16)*10)*PI;

Vhyp

math

math

Das passt sehr schön so.