stepx=1.0e-3
unlocked chinese cellular phones !选择适当的条元宽度
open(31,file='stress_x.txt',status='unknown',form='formatted',access='append')
open(32,file='stress_y.txt',
video door phone nikon schematics status='unknown',form='formatted',access='append')
!do i_x=1,30,3
!#############################
tiphone 3g a8 a520 i_x=7
!#############################
do i_y=0,60,3
sum_value_stress_x=0.0
sum_
china unlocked cell phones value_stress_y=0.0
sum_value_stress_t=0.0
x_kexi=i_x*stepx !when change the coeffience of 0.1 here, change the respond co in the function sub too.
do i_n=1,n
if (i_n.lt.(n+1)/2) then
y_gita=(stepx*i_y-stepx*(i_n-(n+1)/2))/thickness_layer1
else
phones/a1900 y_gita=(stepx*i_y+stepx*(i_n-(n+1)/2))/thickness_layer1
endif
!计算x方向的正应力
call stress_x(x_kexi,y_gita,eps,pai,stepx,thickness_layer1,e_halfspace_1,v_halfspace_1,e_halfspace_2,v_halfspace_2,e_cylinder,v_cylinder,value_stress_x)
sum_value_stress_x=sum_value_stress_x+value_stress_x*p(i_n)
!计算y方向的正应力
hiphone c9 call stress_y(x_kexi,y_gita,eps,pai,stepx,thickness_layer1,e_halfspace_1,v_halfspace_1,e_halfspace_2,v_halfspace_2,e_cylinder,v_cylinder,value_stress_y)
sum_value_stress_y=sum_value_stress_y+value_stress_y*p(i_n)
enddo
!输出到屏幕
phone se450 print *,' i_x=',i_x,' i_y=',i_y
print *,' stress_x=',sum_value_stress_x,' stress_y=',sum_value_stress_y
!输出到
phone no a1221文件
if (i_y.eq.0) then
write(31,*) i_y,i_x,sum_value_stress_x
write(32,*) i_y,i_x,sum_value_stress_y
