Définition d'un Pseudo-Tore "épycycloîdal"






Jean-François COLONNA
www.lactamme.polytechnique.fr
jean-francois.colonna@polytechnique.edu
CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641, Ecole Polytechnique, CNRS, 91128 Palaiseau Cedex, France

[Site Map, Help and Search [Plan du Site, Aide et Recherche]]
[The Y2K bug [Le bug de l'an 2000]]
[Croyez-vous que les Nombres Réels existent dans un ordinateur et que les calculs flottants sont sûrs ?]
[N'oubliez pas de visiter Une Machine Virtuelle à Explorer l'Espace-Temps où vous trouverez plus de 6140 images à la frontière de l'Art et de la Science]
(Site WWW CMAP28 : cette page a été créée le 12/12/2016 et mise à jour le 12/12/2016 13:11:38 -CET-)



/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D E F I N I T I O N   D E S   T R O I S   F O N C T I O N S   ' F '  :                                                     */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Definition :                                                                                                               */
/*                                                                                                                                   */
/*                    Le tore est defini parametriquement                                                                            */
/*                  en fonction des deux parametres 'u'                                                                              */
/*                  (appele aussi 'distance polaire' ou                                                                              */
/*                  'teta' ou encore 'latitude') et 'v'                                                                              */
/*                  (appele aussi 'longitude' ou 'phi') :                                                                            */
/*                                                                                                                                   */
/*                                      XPlan(u,v) = B.cos(C,v) - D.cos(E.v)                                                         */
/*                                      YPlan(u,v) = F.sin(G.v) - H.sin(I.v)                                                         */
/*                                                                              Definition d'une courbe plane "epycycloidale".       */
/*                                                                                                                                   */
/*                                      dXPlan(u,v) = d#XPlan(u,v)                                                                   */
/*                                      dYPlan(u,v) = d#YPlan(u,v)                                                                   */
/*                                                                              La tangente T est donc {+dXPlan(u,v),+dYPlan(u,v).   */
/*                                                                                                                                   */
/*                                      Psi(u,v) = arctg(-dXPlan(u,v),+dYPlan(u,v))                                                  */
/*                                                                              La normale N est perpendiculaire a la normale T      */
/*                                                                              et est donc {+dYPlan(u,v),-dXPlan(u,v)}.             */
/*                                                                                                                                   */
/*                                      Xellipse(u,v) = Ra.cos(u)                                                                    */
/*                                      Yellipse(u,v) = 0                                                                            */
/*                                      Zellipse(u,v) = Rb.sin(u)                                                                    */
/*                                                                              Definition d'une ellipse dans le plan {OX,OZ}.       */
/*                                                                                                                                   */
/*                                      F (u,v) = [Xellipse(u,v)*cos(psi)] - [Yellipse(u,v)*sin(psi)] + XPlan(u,v)                   */
/*                                       x                                                                                           */
/*                                                                                                                                   */
/*                                      F (u,v) = [Xellipse(u,v)*sin(psi)] + [Yellipse(u,v)*cos(psi)] + YPlan(u,v)                   */
/*                                       y                                                                                           */
/*                                                                                                                                   */
/*                                      F (u,v) = Zellipse(u,v)                                                                      */
/*                                       z                                                                                           */
/*                                                                              Definition d'une ellipse dont le centre est le point */
/*                                                                              courant {XPlan(u,v),YPlan(u,v),0} et situee dans le  */
/*                                                                              plan {N,OZ}.                                         */
/*                                                                                                                                   */
/*                  et (parametres par defaut) :                                                                                     */
/*                                                                                                                                   */
/*                                      A=8                                                                                          */
/*                                      a=1                                                                                          */
/*                                      L=4                                                                                          */
/*                                                                                                                                   */
/*                                      B=A+a                                                                                        */
/*                                      C=1                                                                                          */
/*                                      D=L.a                                                                                        */
/*                                      E=(A+a)/a                                                                                    */
/*                                                                                                                                   */
/*                                      F=A+a                                                                                        */
/*                                      G=1                                                                                          */
/*                                      H=L.a                                                                                        */
/*                                      I=(A+a)/a                                                                                    */
/*                                                                                                                                   */
/*                  avec :                                                                                                           */
/*                                                                                                                                   */
/*                                      u E [ 0 , 2.p ]                                                                              */
/*                                                                                                                                   */
/*                                      v E [ 0 , 2.p ]                                                                              */
/*                                                                                                                                   */
/*                  (ou 'p' designe 'pi').                                                                                           */
/*                                                                                                                                   */
/*                    Tout ceci est tres inspire de 'v $xtc/epycycloide.03$c'                                                        */
/*                  en faisant :                                                                                                     */
/*                                                                                                                                   */
/*                                      u = phi  E [0,2.p]                                                                           */
/*                                      v = teta E [0,2.p]                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/


(Nota : les lignes d'explications qui précèdent sont des commentaires extraits des programmes ayant été utilisés pour calculer les images correspondantes. Ce programme en est un exemple parmi des centaines.)


Copyright (c) Jean-François Colonna, 2016.
Copyright (c) CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2016.