二次近似曲線の式を定義します。
//**********************************************************************************
#macro Contact_curve_function(Sp,E1,E2,X0,Z0,Xc,Zc,X1,Z1,Px0,Pz0)
#local Px0=(X0*pow(1-Sp,2)+(2+E1)*Xc*Sp*(1-Sp)+
(1+E1+E2)*X1*pow(Sp,2))/(1+E1*Sp+E2*pow(Sp,2));
#local Pz0=(Z0*pow(1-Sp,2)+(2+E1)*Zc*Sp*(1-Sp)+
(1+E1+E2)*Z1*pow(Sp,2))/(1+E1*Sp+E2*pow(Sp,2));
#end
//**********************************************************************************