ConchoidSurface::usage="The definition of a conchoid surface is:​
r
c
​[u, v] = r[u, v] ± k
r[u,v]-P
r[u,v]-P
​, P is the pole, k distance and r[u, v] the base surface.In this notebook we will obtain the conchoid surface of the cone: ​{v Cos[u]​, v Sin[u]​, v} with the distance k = 1.5 and as pole P the {0, 0, 1}.Reference: Peternell, M. et al.Conchoid surfaces of rational ruled surfaces. Vienna University of Technology; Personal. ";
Graph:
In[]:=
SwatchLegend[{Orange,Blue,Yellow},{"Positive conchoid surface","Negative conchoid surface","Base surface"}]​​Showr[u_,v_]:={vCos[u],vSin[u],v};k=1.5;P={0,0,1};r1[u_,v_]:=r[u,v]+k
r[u,v]-P
Norm[r[u,v]-P]
;​​r2[u_,v_]:=r[u,v]-k
r[u,v]-P
Norm[r[u,v]-P]
;{ParametricPlot3D[r1[u,v],{u,0,2π},{v,-2,2},PerformanceGoal"Quality",PlotStyle{Opacity[0.3],Orange},PlotPoints60],​​ParametricPlot3D[r2[u,v],{u,0,2π},{v,-2,2},PerformanceGoal"Quality",PlotStyle{Opacity[0.3],Blue},PlotPoints60],ParametricPlot3D[r[u,v],{u,0,2π},{v,-3,3},MeshFalse,PlotPoints60,PlotStyle->{Yellow,Opacity[0.3]}]},AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},AxesTrue,BoxedFalse,BoxRatiosAutomatic,PlotLabelStyle[Framed["Conchoid surface of a cone"],16,Blue,BackgroundYellow],PlotRange{{-3,3},{-3,3},{-5,5}},ViewPoint{-7,-2.4,2.},ImageSize700
Out[]=
Positive conchoid surface
Negative conchoid surface
Base surface
Out[]=
Explanation in 2D:
In[]:=
SwatchLegend[{Orange,Blue,DarkYellow,Red,Black},{"Positive Conchoid","Negative Conchoid","Original function","Pole","Generatrix line"}]​​Manipulatef[t_]:=t;
f
1
[t_]:=-t;g[t_]:=t;
x
0
=0;
y
0
=1;k=1.5;
ζ
+
[t_]:=f[t]+
k(f[t]-
x
0
)
2
(f[t]-
x
0
)
+
2
(g[t]-
y
0
)
,g[t]+
k(g[t]-
y
0
)
2
(f[t]-
x
0
)
+
2
(g[t]-
y
0
)
;​​
ζ
-
[t_]:=f[t]-
k(f[t]-
x
0
)
2
(f[t]-
x
0
)
+
2
(g[t]-
y
0
)
,g[t]-
k(g[t]-
y
0
)
2
(f[t]-
x
0
)
+
2
(g[t]-
y
0
)
;​​
ζ
1
+
[t_]:=
f
1
[t]+
k(
f
1
[t]-
x
0
)
2
(
f
1
[t]-
x
0
)
+
2
(g[t]-
y
0
)
,g[t]+
k(g[t]-
y
0
)
2
(
f
1
[t]-
x
0
)
+
2
(g[t]-
y
0
)
;​​
ζ
1
-
[t_]:=
f
1
[t]-
k(
f
1
[t]-
x
0
)
2
(
f
1
[t]-
x
0
)
+
2
(g[t]-
y
0
)
,g[t]-
k(g[t]-
y
0
)
2
(
f
1
[t]-
x
0
)
+
2
(g[t]-
y
0
)
;​​Show[{ParametricPlot[
ζ
+
[t],{t,-3,t1},PlotStyle{Orange,Thick},PerformanceGoal->"Quality"],ParametricPlot[
ζ
-
[t],{t,-3,t1},PlotStyle{Blue,Thick},PerformanceGoal->"Quality"],​​ParametricPlot[
ζ
1
+
[t],{t,-3,t1},PlotStyle{Orange,Thick},PerformanceGoal->"Quality"],ParametricPlot[
ζ
1
-
[t],{t,-3,t1},PlotStyle{Blue,Thick},PerformanceGoal->"Quality"],ParametricPlot[{f[t],g[t]},{t,-3,3},PlotStyle{DarkYellow,Thick}],ParametricPlot[{
f
1
[t],g[t]},{t,-3,3},PlotStyle{DarkYellow,Thick}],Graphics[{Orange,PointSize[0.015],Point[
ζ
+
[t1]]}],Graphics[{Orange,PointSize[0.015],Point[
ζ
1
+
[t1]]}],Graphics[{Blue,PointSize[0.015],Point[
ζ
-
[t1]]}],Graphics[{Blue,PointSize[0.015],Point[
ζ
1
-
[t1]]}],Graphics[{Black,Thick,Line[{
ζ
-
[t1],{
x
0
,
y
0
},
ζ
+
[t1]}]}],Graphics[{Black,Thick,Line[{
ζ
1
-
[t1],{
x
0
,
y
0
},
ζ
1
+
[t1]}]}],Graphics[{Red,PointSize[0.015],Point[{
x
0
,
y
0
}]}],Graphics[{DarkYellow,PointSize[0.015],Point[{f[t1],g[t1]}]}],Graphics[{DarkYellow,PointSize[0.015],Point[{
f
1
[t1],g[t1]}]}]},AxesTrue,AxesOrigin{0,0},AspectRatioAutomatic,ImageSize700,PlotRange{{-4,4},{-5,6}},PlotLabelStyle[Framed["Conchoid"],16,Blue,BackgroundYellow]],{{t1,-3+0.00002,"Value (t1)"},-3+0.00001,3,0.00001},ControlPlacement->Top​​
Out[]=
Positive Conchoid
Negative Conchoid
Original function
Pole
Generatrix line
Out[]=
​
Value (t1)