研究生隧道工程學(xué)課程考核格式及作業(yè)-資料下載頁
2025-06-28 05:17本頁面
【正文】 in=K kmin=K。 x=x./b。 vector_kmin=x。 end endendvector_kmax,kmax,vector_kmin,kminfunction K=funx4(x)h=*x(3)*tan(x(4))。Hq=x(6)。if h=HqG=*x(2)*x(3)^2*tan(x(4))。F=*x(1)*(1+tan(x(5)))*x(3)^2/x(6)/(cos(x(4)))^2。else G=*x(3)*h*(1((hHq)/h)^2)。 F=*x(1)*(1+tan(x(5)))*x(3)^2/x(6)/(cos(x(4)))^2*(Hq/h)^2+x(1)*(1Hq/h)*x(3)。end。K=F/G。程序3 model3下計(jì)算K值的范圍x2=20。 x3=1。x5=1*pi/180。side2=5。kmax=0。kmin=100000。vector_kmax=[0 0 0 0 0 0 0]。vector_kmin=[0 0 0 0 0 0 0]。b=[1 1 1 pi/180 pi/180 1 1]。side=[x3 side2]。 Hq=fun2(side)。 x1=1 for x4=pi/360:pi/180:85*pi/180 x=[x1 x2 x3 x4 x5 Hq side2]。 K=funx4(x)。 if kmax=K kmax=K。 x=x./b。 vector_kmax=x。 end if kmin=K kmin=K。 x=x./b。 vector_kmin=x。 end endvector_kmax,kmax,vector_kmin,kminfunction[Hq]=fun2(side)if side(1)=5 a1=。else a1=。endw=1+a1*(side(1)5)。H1=w**2^(side(2)1)。 if side(2)=3 a=2。else a=。endHq=a*H1。程序4 model4大跨度圍巖分級(jí)x2=20。%209。210。205。193。204。229。214。216。182。200。x3=100。%191。231。182。200。x1=50。%213。179。190。219。193。166。x5=10*pi/180。%196。218。196。166。178。193。189。199。side2=5。%178。224。195。230。206。167。209。210。188。182。177。240。if side2==1 x4=*pi/180。elseif side2==2 x4=*pi/180。elseif side2=3amp。side2=5 x4=*pi/180。else sprintf(39。a wrong surrounding rock classification39。)endif side2==1amp。x5=45*pi/180|side2==2amp。x5=45*pi/180amp。x5=35*pi/180.... |side2==3amp。x5=35*pi/180amp。x5=25*pi/180|side2==4amp。x5=25*pi/180.... amp。x5=15*pi/180|side2==5amp。x5=15*pi/180amp。x5=1*pi/180 returnelse sprintf(39。a wrong angle of internal friction39。)endif side2==1amp。x1=400|side2==2amp。x1=400amp。x1=300.... |side2==3amp。x1=300amp。x1=200|side2==4amp。x1=200.... amp。x1=100|side2==5amp。x1=100amp。x1=1 returnelse sprintf(39。a wrong cohesion39。)endb=[1 1 1 pi/180 pi/180 1 1]。side=[x3 side2]。%208。206。179。201。sideHq=fun2(side)。%211。176。207。236。201。238。182。200。x=[x1 x2 x3 x4 x5 Hq side2]。h=*x(3)*tan(x(4))。Hq=x(6)。if h=HqG=*x(2)*x(3)^2*tan(x(4))。 F=*x(1)*(1+tan(x(5)))*x(3)^2/x(6)/(cos(x(4)))^2。 else G=*x(3)*h*(1((hHq)/h)^2)*x(2)。 F=*x(1)*(1+tan(x(5)))*x(3)^2/x(6)/(cos(x(4)))^2*(Hq/h)^2+x(1)*(1Hq/h)*x(3)。end。K=F/G,x3if K=0amp。K= sprintf(39。162。245。188。182。206。167。209。210。39。)elseif K=amp。K=8 sprintf(39。162。244。188。182。206。167。209。210。39。) elseif K=8amp。K=23 sprintf(39。162。243。188。182。206。167。209。210。39。) elseif K=23amp。K=105 sprintf(39。162。242。188。182。206。167。209。210。39。) elseif K=105 sprintf(39。162。241。188。182。206。167。209。210。39。)else sprintf(39。a wrong result39。) endXXXVII /
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1