EnvelopeHyperboloids::usage="The equation of a family of hyperboloids is: + - = 1. Find its envelope, its characteristic and its recoil edge.Reference: Personal.";
2
(x-c)
4
2
(y-c)
9
2
z
16
Development:
In[]:=
Solve+-==1,+--1==0,{x,y}
2
(x-c)
4
2
(y-c)
9
2
z
16
∂
c
2
(x-c)
4
2
(y-c)
9
2
z
16
Out[]=
x13c-52c+913c+52c-9
1
13
13
16+
,y2
z
1
52
13
16+
,x2
z
1
13
13
16+
,y2
z
1
52
13
16+
2
z
Envelope:
ζ[t_]:=13c-52c+9[t_]:=13c+52c-9
1
13
13
16+
,2
t
1
52
13
16+
,t;2
t
ζ
1
1
13
13
16+
,2
t
1
52
13
16+
,t;2
t
Characteristic:
For each value of c we have a different characteristic.
Recoil edge :
In[]:=
Solve+-==1,+--1==0,+--1==0==0,{x,y,z}
2
(x-c)
4
2
(y-c)
9
2
z
16
∂
c
2
(x-c)
4
2
(y-c)
9
2
z
16
∂
c,c
2
(x-c)
4
2
(y-c)
9
2
z
16
{};
As it is not a developable surface, it has no recoil edge. But if it has a line of striction we will obtain it from the intersection of the family with the plane z = 0:
ξ[θ_]:={2Cos[θ]+c,3Sin[θ]+c,0};
A family of ellipses.
Animation:
In[]:=
SwatchLegend[{Yellow,Gray,Red},{"Envelope and characteristics","Family elements","Striction lines"}]Manipulateζ[t_]:=13c-52c+9[t_]:=13c+52c-9[t],{t,-20,20},PlotStyle->Yellow]/.Line[pts_,rest___]:>Tube[pts,0.09,rest],{c,-5,t1,1}],Table[ParametricPlot3D[ξ[θ],{θ,0,2π},PlotStyleRed]/.Line[pts_,rest___]:>Tube[pts,0.09,rest],{c,-5,t1,1}],Table[ParametricPlot3D[σ[u,v],{u,-5,5},{v,0,2π},PlotStyle->{Opacity[0.1],Specularity[White,10],Gray},MeshFalse,PerformanceGoal"Quality"],{c,-5,t1,1}],ParametricPlot3D[ζ[t],{t,-20,20},{c,-5,t1},PlotStyle->{Opacity[0.1],Specularity[White,10],Yellow},MeshFalse,PerformanceGoal"Quality"],ParametricPlot3D[[t],{t,-20,20},{c,-5,t1},PlotStyle->{Opacity[0.1],Specularity[White,10],Yellow},MeshFalse,PerformanceGoal"Quality"]},AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},AxesTrue,BoxedFalse,BoxRatiosAutomatic,PlotRangeAll,ViewPoint{8,-2.4,4.},ImageSize900],Style["Envelope of a family of hyperboloids:",Bold,Large],{{t1,-5+0.00002,"Value (t1)"},-5+0.00001,4.9999,1},ControlPlacementTop
1
13
13
16+
,2
t
1
52
13
16+
,t;2
t
ζ
1
1
13
13
16+
,2
t
1
52
13
16+
,t;ξ[θ_]:={2Cos[θ]+c,3Sin[θ]+c,0};σ[u_,v_]:=22
t
1+
Cos[v]+c,32
u
1+
Sin[v]+c,4u;Show[{Table[ParametricPlot3D[ζ[t],{t,-20,20},PlotStyle->Yellow]/.Line[pts_,rest___]:>Tube[pts,0.09,rest],{c,-5,t1,1}],Table[ParametricPlot3D[2
u
ζ
1
ζ
1
Out[]=
Out[]=

