E1 70[GPa] "Block Young's modulus"
E2 210[GPa] "Cylinder Young's modulus"
nu0 0.3 "Poisson's ratio"
Fn 35[kN] "External load"
E_star 2*E1*E2/((E1+E2)*(1-nu0^2)) "Combined Young's modulus"
R 50[mm] "Combined radius"
d 200[mm] "Block width"
th 1[mm] Thickness
lc 10[mm] "Estimated contact length"
a sqrt(8*Fn*R/(pi*E_star*th)) "Analytical contact length"
pmax sqrt(Fn*E_star/(2*pi*R*th)) "Maximum contact pressure"
dist 1[mm] "Initial distance between parts"
k 0 "Spring coefficient"