Tor einer Lagerhalle Fa. Ernesting in Coesfeld (Architekt S. Calatrava)

[Maple OLE 2.0 Object]

Herleitung der Fächengleichung

[Maple OLE 2.0 Object]

Definitionen: Torhöhe h

Torhöhe offen h1

Parabel im geschlossenen Zustand z=p y^2 +s

Da zz+(h1-zz)=h gilt ist der Kreuzriss der von den Punkten P beschriebenen Leitkurve eine Ellipse.

> restart;
with(plots):

Warning, the name changecoords has been redefined

> k1:=x^2+z^2-zz^2;k2:= x^2+(z-h1)^2-(-zz-h)^2;

k1 := x^2+z^2-zz^2

k2 := x^2+(z-h1)^2-(-zz-h)^2

> solve({k1,k2},{x,z});

{z = -1/2*(-h1^2+2*zz*h+h^2)/h1, x = 1/2*RootOf(h1^...

> allvalues(%);

{z = -1/2*(-h1^2+2*zz*h+h^2)/h1, x = 1/2*(-h1^4+4*h...
{z = -1/2*(-h1^2+2*zz*h+h^2)/h1, x = 1/2*(-h1^4+4*h...

> L:=%[1];

L := {z = -1/2*(-h1^2+2*zz*h+h^2)/h1, x = 1/2*(-h1^...

> L1:=(subs(zz=p*y^2+s,L));

L1 := {z = -1/2*(-h1^2+2*(p*y^2+s)*h+h^2)/h1, x = 1...

> L1[1];

z = -1/2*(-h1^2+2*(p*y^2+s)*h+h^2)/h1

> lz:=op(2,%);

lz := -1/2*(-h1^2+2*(p*y^2+s)*h+h^2)/h1

> L1[2];

x = 1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(p*y...

> lx:=op(2,%);

lx := 1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(p...

Damit haben wir die Leitkurve konstruiert:

> LK:=[lx,y,lz];

>

LK := [1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(...

> LK1:=subs(h=5,h1=-3,s=-3,p=.1,LK);

LK1 := [-1/6*sqrt(704-32.0*y^2-64*(.1*y^2-3)^2), y,...

> spacecurve(LK1,y=-5..5,color=black,thickness=3,axes=framed,orientation=[-130,90]);

[Maple Plot]

FL ist nun die obere Regelfläche und FU ist die untere, sie haben die Leitkurve gemeinsam und sind gerade Konoide mit der jeweiligen Torkante als eigentlicher Leitgeraden und der Ferngeraden der xz-Ebene als uneigentlicher Leitgeraden

> FL:=[lx-lambda*lx,y,lz-lambda*lz];

FL := [1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(...
FL := [1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(...

> FU:=[lx-lambda*lx,y,lz-lambda*(lz-h1)];

FU := [1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(...
FU := [1/2*(-h1^4+4*h1^2*(p*y^2+s)*h+2*h^2*h1^2-4*(...

> FL1:=subs(h=5,h1=-2,s=-3,p=.05,FL);

FL1 := [-1/4*sqrt(819-21.00*y^2-84*(.5e-1*y^2-3)^2)...

> Fu1:=subs(h=5,h1=-2,s=-3,p=.05,FU);

Fu1 := [-1/4*sqrt(819-21.00*y^2-84*(.5e-1*y^2-3)^2)...

> OB:=plot3d(FL1,y=-5..5,lambda=0..1):

> UN:=plot3d(Fu1,y=-5..5,lambda=0..1):

> display3d({OB,UN},orientation=[180,90],axes=framed,scaling=constrained);

[Maple Plot]

> display3d({OB,UN},orientation=[-130,90],axes=framed,scaling=unconstrained);

[Maple Plot]

>