<< CSPlotter.m

Radial contact structure in R^3

rmax = 6;  dr = rmax/13;
thmax = 2  Pi;  dth = 2  Pi/32;
scale = dr 0.999;

conormalright[x_, y_, z_] = {-y/2, x/2, 
          1};
conormalleft[x_, y_, z_] = {y/2, -x/2, 
          1};

Table[
        contactelement[
{r Cos[θ], 
            r Sin[θ], 0}, 
conormalleft[
            r Cos[θ], r Sin[θ], 0], 
          scale], 
{r, 0, rmax, 
          dr}, 
{θ, 0, Pi 2, 
          dth}
];
rout = 
      Show[%, coords[rmax + 1, rmax + 1, 5.6], 
        Boxed -> False]
Export["Radialleft.eps", rout];
Export["Radialleft.png", rout];
Export["Radialleft.pdf",  
        rout];

Table[
        contactelement[
{r Cos[θ], 
            r Sin[θ], 0}, 
conormalright[
            r Cos[θ], r Sin[θ], 
            0], 
scale], 
{r, 0, rmax, 
          dr}, 
{θ, 0, Pi 2, 
          dth}
];

lout  =  
      Show[%, coords[rmax + 1, rmax + 1, 5.6], 
        Boxed -> False]
Export["Radialright.eps", lout];
Export["Radialright.png", lout];
Export["Radialright.pdf",  lout];