RadiusCurvatureSpiral::usage="Demonstrate that the radius of curvature of a conical spiral a Cos[φ ],a Sin[φ ],b isproportional to the distance between a point of the spiral and the axis of the cone.References: A.S. Fedenko, Problemas de Geometría Diferencial. Editorial Mir: Moscú, p.68, 1981.";
hφ
hφ
hφ
Development:
In[]:=
r[φ_]:=aCos[φ],aSin[φ],b;rc=
hφ
hφ
hφ
3
(Norm[r'[φ]])
Norm[r'[φ]r''[φ]]
In[]:=
Simplify++
3/2
++
2Re[hφ]
2
(bh)
2
ahCos[φ]-aSin[φ]
hφ
hφ
2
aCos[φ]+ahSin[φ]
hφ
hφ
2
-abCos[φ]+abhSin[φ]
2hφ
2
h
2hφ
2
-abhCos[φ]-abSin[φ]
2hφ
2hφ
2
h
2
+++
2
a
2hφ
2
Cos[φ]
2
a
2hφ
2
h
2
Cos[φ]
2
a
2hφ
2
Sin[φ]
2
a
2hφ
2
h
2
Sin[φ]
3/2
(+(1+))
2
b
2Re[hφ]
2
h
2
a
2hφ
2
h
2
a
4hφ
2
h
2
b
2
h
2
a
2
h
3hφ
3/2
(+(1+))
2
b
2
h
2
a
2
h
a
2hφ
1+
2
h
hφ
2
b
2
h
2
a
2
h
a
1+
2
h
Distance between a point of the spiral and the axis of the cone:
In[]:=
d=
Norm[aCos[φ],aSin[φ],b-{0,0,0}{0,0,1}]
hφ
hφ
hφ
Norm[{0,0,1}]
2Re[hφ]
2
(aCos[φ])
2Re[hφ]
2
(aSin[φ])
hφ
If we divide both quantities:
In[]:=
hφ
2
b
2
h
2
a
2
h
a
1+
2
h
a
hφ
Out[]=
(+(1+))
2
b
2
h
2
a
2
h
2
a
1+
2
h
It is a constant amount, which means that they are proportional amounts.
Gráfica, sean a = 1, b = , h = :
1
1
4
In[]:=
ManipulateClearAll;r[φ_]:=φCos[φ],φSin[φ],φ;a=1;b=1;h=;d1[t1_]:=(+(1+));d2[t1_]:=a;p=Grid[{{CurvatureRadius,ConeAxisSpiralDistance,DistancesRatio},{d1[t1],d2[t1],d1[t1]/d2[t1]}},FrameAll,ItemSize->15,ItemStyle->15];Show[{ParametricPlot3D[r[φ],{φ,-4Log[t1],4Log[t1]},PlotStyle{Thick,Blue}],ParametricPlot3D[-r[φ],{φ,-4Log[t1],4Log[t1]},PlotStyle{Thick,Blue}],ParametricPlot3D[{vCos[u],vSin[u],v},{u,0,2π},{v,-20,20},PlotStyleDirective[Green,Specularity[White,20],Opacity[0.1]],MeshFalse,PerformanceGoal"Quality"],Graphics3D[{Black,Dashed,Thick,Line[{{0,0,t1},r[4Log[t1]]}]}],Graphics3D[{Black,Dashed,Thick,Line[{{0,0,-t1},-r[4Log[t1]]}]}],Graphics3D[{Red,Ball[{0,0,t1},0.43]}],Graphics3D[{Red,Ball[r[4Log[t1]],0.43]}],Graphics3D[{Red,Ball[-r[4Log[t1]],0.43]}],Graphics3D[{Red,Ball[{0,0,-t1},0.43]}]},AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},AxesTrue,BoxedFalse,BoxRatiosAutomatic,PlotRangeAll,ViewPoint{-3,3,2},ImageSize500,PerformanceGoal"Quality"],Style["Conical Spiral:",Bold,Large],{{t1,0.1,"Value (t1)"},0.01,20,0.01},Delimiter,{{p,2,"Measurable parameters"}},ControlPlacementTop
1
4
1
4
1
4
1
4
h(4Log[t1])
2
b
2
h
2
a
2
h
a
1+
2
h
h(4Log[t1])
Out[]=

