附录D 此程序用于计算受
best china phone椭圆分布压力的半空间体内的应力σx，σy
主程序名：stress.f90
parameter (n=41)
dimension p(n)
china phone !光滑的data (p(i),i=1,41)/5.9700012E-04,1.6415119E-03,2.1514073E-03,2.5552623E-03, 2.8793626E-03,3.1586811E-03,3.3939332E-03,3.6049113E-03,3.7837699E-03,3.9469935E-03,
flip phone ttk 4.0850602E-03,4.2105112E-03,4.3150913E-03,4.4099651E-03,4.4862051E-03,4.5492928E-03, 4.6069678E-03,4.6492573E-03,4.6739280E-03,4.6903510E-03,4.6991799E-03,4.6903323E-03, 4.6739057E-03,4.6492890E-03,4.6069380E-03,4.5493543E-03,4.4861548E-03,4.4099726E-03, 4.3150689E-03,4.2105401E-03,4.0850714E-03,3.9469711E-03,3.7838016E-03,3.6049038E-03, 3.3939546E-03,3.1586601E-03,2.8793511E-03,2.5552886E-03,2.1514120E-03,1.6415853E-03, 5.9686293E-04/
!粗糙的
china phones free shippingdata (p(i),i=1,41)/0.0,0.2175937,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.1650926, 0.0,0.0,0.0,0.0, 0.1806793, 0.0, 0.0, 1.6479224E-02,9.3706310E-02, 0.0,0.0,0.0,0.1879733, 0.0,0.0,0.0,0.0, 7.0323713E-02,0.0, 0.0,0.0,3.8797136E-02,0.0,0.0,0.0,0.0,3.4677655E-02,0.0,0.0,0.0/
n0=301 !轮廓接触线长度为301个条元宽度
pai=3.14159265359 !圆周率
thickness_layer1=1.0 !
chinese unlocked cell phones 涂层厚度
e_halfspace_1=0.21*2.0 !涂层1的弹性模量
v_halfspace_1=0.3 !涂层1的泊松比
e_halfspace_2=0.21 !基体的弹性模量
v_halfspace_2=0.3 !基体的泊松比
e_cylinder=0.21 !
iphone 4gs assembled in china压头的弹性模量
v_cylinder=0.3 !压头的泊松比
r_e1e2=e_halfspace_1/e_halfspace_2 !涂层与基体的弹性模量之比
r=7.8248776E-02 !选择适当大小的压头
eps=1.0e-5 !误差参数
