Plot of the equilibrium shapes of an elastic ring with lines of contact

Case n=9

Contents of this notebook :
Part I. Definitions
Part II. Scale reduction and determination of the contact values of the parameters q, η and s1 for the ring of the reduced length
Part III. Plot of the equilibrium shape with lines of contact
Part IV. Confirmation of the contact conditions and the length

Part I. Definitions

Definitions from Section 4. Analytic expressions for the ring curvatures and slope angles

μ=(1/4) ((α^2)+(β^2)+(γ^2)+α β+α γ+β γ);

γ=-((α+β)/2)+I η ;

λ=(1/4)(Sqrt[A B]); k= (1/2)-((4 (η^2)+(3 α+β) (α+3 β))/(2 A B));

A= (Sqrt[4 (η^2)+(((3 α+β))^2)]); B= (Sqrt[4 (η^2)+ (((α+3 β))^2)]) ;

T=(4/λ)  EllipticK[k];

Left-hand side of the closure condition -- formula (42)

CC= (((A+B) (α-β))/(2 λ (A-B))) EllipticPi[-((((A-B))^2)/(4 A B)),k]+((A β-B α)/(λ (A-B))) EllipticK[k];

New parametrization -- formula (44)

α=((4 σ)/((η^2)+(q^2)))-q; β=((4 σ)/((η^2)+(q^2)))+q;

Input data (contact pressure (Subscript[σ, cn])=464.276)

σ=800; n=9;

Part II. Scale reduction and determination of the contact values of the parameters q, η and s1 for the ring of the reduced length

The contact parameters for a ring of length 2π
 (These can be determined using the notebook Contact.nb. Here, 150 digit precision is used to ensure sufficient accuracy of the computations.)

σ9cn= 464.27591167203005316900723413758081842284745815395936945922213196602376300295145803170315572975892713637252192634430427905552802884788150437492867373612792168518884732`150.;

q9cn=18.50573713223354970890536829136729185743330537365335122913778156232947687087005587591529067931840224337102721193188183947600929108889829261760221151916465850339119657`150.;

η9cn=13.32497064961571646649775179800224354005530665171299945205627047146762956432651842004401591623003515440662170119176698468519944680410475939972825410355494126179869746`150.;

s19cn =0.0933379966441185869355166230068658285389135714993104101029023175155798723723273308540203583274908614560875839215803067667393245315228732799104769;

Scalling

f=(((((σ9cn)/σ)))^(1/3));

q=q9cn/f; η=η9cn/f; s1 = s19cn f;

Part III. Plot of the equilibrium shape with lines of contact

Expressions for the curvature, its derivative and the slope angle

κ[s_]=(((A β+B α)-(A β-B α) JacobiCN[λ s,k])/((A+B)-(A-B) JacobiCN[λ s,k]));

κ1[s_]=(∂_s) κ[s];

∅[s_]= ((A β-B α)/(A-B)) s+(((A+B) (α-β))/(2 λ (A-B))) EllipticPi[-((((A-B))^2)/(4 A B)),JacobiAmplitude[λ s,k],k]+((α-β)/(2 λ (Sqrt[k+((((A-B))^2)/(4 A B))])))  ArcTan[(Sqrt[k+((((A-B))^2)/(4 A B))])((JacobiSN[λ s,k])/(JacobiDN[λ s,k]))];

Parametric equations

X[s_]=(1/σ) κ1[s] Cos[∅[s]]+(( ((κ[s])^2)-2 μ)/(2 σ)) Sin[∅[s]];

Y[s_]=(1/σ) κ1[s] Sin[∅[s]]-(( ((κ[s])^2)-2 μ)/(2 σ)) Cos[∅[s]];

Plot

Calculation of the total length LL of all line segments

LL=((2 π)/n)(1-f);

Plots of the parts of the ring shape

F1=ParametricPlot[{{X[s],Y[s]},{X[s+T],Y[s+T]},{X[s+2 T],Y[s+2 T]},{X[s+3 T],Y[s+3 T]},{X[s+4 T],Y[s+4 T]},{X[s+5T],Y[s+5 T]},{X[s+6T],Y[s+6 T]},{X[s+7T],Y[s+7 T]},{X[s+8T],Y[s+8 T]}},{s,0,s1},AspectRatio->Automatic, Axes -> False,PlotStyle->{{Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick}}];

s2=T-s1;

F2=ParametricPlot[{{X[s],Y[s]},{X[s+T],Y[s+T]},{X[s+2 T],Y[s+2 T]},{X[s+3 T],Y[s+3 T]},{X[s+4 T],Y[s+4 T]},{X[s+5T],Y[s+5 T]},{X[s+6T],Y[s+6 T]},{X[s+7T],Y[s+7 T]},{X[s+8T],Y[s+8 T]}},{s,s2,T},AspectRatio->Automatic, Axes -> False,PlotStyle->{{Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick}}];

ξ1=(Sqrt[1+(LL/(Sqrt[((X[s1])^2)+((Y[s1])^2)]))]);

FL=ParametricPlot[{{X[s1] ξ,Y[s1] ξ},{X[s1+T] ξ,Y[s1+T] ξ},{X[s1+2 T] ξ,Y[s1+2 T] ξ},{X[s1+3 T] ξ,Y[s1+3 T] ξ},{X[s1+4 T] ξ,Y[s1+4 T] ξ},{X[s1+5 T] ξ,Y[s1+5 T] ξ},{X[s1+6 T] ξ,Y[s1+6 T] ξ},{X[s1+7 T] ξ,Y[s1+7 T] ξ},{X[s1+8 T] ξ,Y[s1+8 T] ξ}}, {ξ,1,ξ1},AspectRatio->Automatic, Axes -> False,PlotStyle->{{Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick}}];

F3=ParametricPlot[{{(ξ1-1)  X[s1]+ X[s],(ξ1-1) Y[s1]+Y[s]},{(ξ1-1)  X[s1+T]+ X[s+T],(ξ1-1) Y[s1+T]+Y[s+T]},{(ξ1-1)  X[s1+2 T]+ X[s+2 T],(ξ1-1) Y[s1+2 T]+Y[s+2 T]},{(ξ1-1)  X[s1+3 T]+ X[s+3 T],(ξ1-1) Y[s1+3 T]+Y[s+3 T]},{(ξ1-1)  X[s1+4 T]+ X[s+4 T],(ξ1-1) Y[s1+4 T]+Y[s+4 T]},{(ξ1-1)  X[s1+5 T]+ X[s+5 T],(ξ1-1) Y[s1+5 T]+Y[s+5 T]},{(ξ1-1)  X[s1+6 T]+ X[s+6 T],(ξ1-1) Y[s1+6 T]+Y[s+6 T]},{(ξ1-1)  X[s1+7 T]+ X[s+7 T],(ξ1-1) Y[s1+7 T]+Y[s+7 T]},{(ξ1-1)  X[s1+8 T]+ X[s+8 T],(ξ1-1) Y[s1+8 T]+Y[s+8 T]}},{s,s1,s2},AspectRatio->Automatic, Axes -> False,PlotStyle->{{Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick},{Black,Thick}, {Black,Thick}, {Black,Thick}}];

Combination of these plots gives the ring shape with lines of contact

FC=Show[F3, F1, F2,FL]

Export["C:\\Lines_9_fold.eps", FC];

Part IV. Confirmation of the contact conditions and the length

N[CC-(π/(2 n))]

N[(T/4)-(π/(2 n)) f]

N[∅[s1]-((( π)/n)-(( π)/2))]

N[n T+n LL-2 π]