495发动机动力计算vb实验答案

以下为代码：PrivateSubCommand1_Click()Picture1.ClsPicture1.PrintPicture1.PrintPicture1.Print初始数据Picture1.PrintPicture1.Print机型495Q型柴油机Picture1.Print主要原始参数Picture1.Print气缸直径D(mm)95Picture1.Print活塞行程S(mm)115Picture1.Print连杆长度L(mm)210Picture1.Print压缩比20Picture1.Print供油提前角:18+2Picture1.Print进气门开:前17Picture1.Print进气门关:后47Picture1.Print排气门开:前47Picture1.Print燃烧室型式涡流式Picture1.Print发动机型式四冲程、水冷、直列Picture1.Print发火顺序1-3-4-2Picture1.Print曲柄销Dk直径m0.062Picture1.PrintLs长度m0.028Picture1.Print主颈Dk直径m0.07Picture1.PrintLs长度m0.04fh=π*d^2/4'求运动值ForS1=0To360Step1α=S1*π/180SIB=λ*Sin(α)COB=Sqr(1-λ^2*(Sin(α))^2)EX(S1)=r*1000*(1-Cos(α))+l*1000*(1-COB)FV(S1)=r*w*(Sin(α)+Cos(α)*SIB/COB)BB=Cos(α)-Sin(α)*SIB/COB+λ*(Cos(α))^2/COB^3GA(S1)=r*w^2*BBJIAO(S1)=S1NextS1fh=π*d^2/4ForS1=0To360Step1α=S1*π/180SIB=λ*Sin(α)COB=Sqr(1-λ^2*(Sin(α))^2)EX(S1)=r*1000*(1-Cos(α))+l*1000*(1-COB)FV(S1)=r*w*(Sin(α)+Cos(α)*SIB/COB)BB=Cos(α)-Sin(α)*SIB/COB+λ*(Cos(α))^2/COB^3GA(S1)=r*w^2*BBJIAO(S1)=S1NextS1pa=0.9*po:HC=s/(ε-1)'求气体力值Fori=0To180Step10a1=i/10:PG(a1)=(0.9-1)*poNextiForI1=190To350Step10n1=1.41-100/n:A2=I1/10PG(A2)=pa*((s+HC)/(EX(I1)/1000+HC))^n1-poNextI1pz=7.5:pc=pa*ε^n1PG(36)=pc+0.45*(pz-pc)-po:PG(37)=7.5pb=0.39:n2=1.18Fori3=380To540Step10A3=i3/10:AA3=i3-360:n2=1.18PG(A3)=pb*((s+HC)/(EX(AA3)/1000+HC))^n2-poNexti3Fori4=550To720Step10P5=(1.15-1)*po:A4=i4/10PG(A4)=P5Nexti4ForAC=0To36A3=AC*10:MJ=1.838PJJ(AC)=-MJ*GA(A3)/fh*10^-6PJJ(AC+36)=-MJ*GA(A3)/fh*10^-6NextACForii=0To720Step10aa=ii*π/180SIB=λ*Sin(aa):COB=Sqr(1-λ^2*(Sin(aa)^2))A2=ii/10:PE=PG(A2)+PJJ(A2):PP(A2)=PEPHH(A2)=PP(A2)*SIB/COBPCC(A2)=PP(A2)/COBPTT(A2)=PP(A2)*(Sin(aa)+Cos(aa)*SIB/COB)PNN(A2)=PP(A2)*(Cos(aa)-Sin(aa)*SIB/COB)NextiiForI1=10To720Step10'求矩aa=I1/10M(1,aa)=PTT(aa)*r*fh*10^6NextI1Foraa=1To72Ifaa+1872ThenC1=aa-54M(2,aa)=M(1,C1)ElseC1=aa+18M(2,aa)=M(1,C1)EndIfIfaa+3672ThenC2=aa-36M(4,aa)=M(1,C2)ElseC2=aa+36M(4,aa)=M(1,C2)EndIfIfaa+5472ThenC3=aa-18M(3,aa)=M(1,C3)ElseC3=aa+54M(3,aa)=M(1,C3)EndIfNextaaForii=0To72'求合成扭矩ForII1=1To5Mm=0ForII2=0ToII1-1Mm=Mm+M(II2,ii)NextII2Z(II1,ii)=MmNextII1NextiiForff=0To72ForFF1=1To4Q(FF1,ff)=Z(FF1,ff)+M(FF1,ff)/2NextFF1NextffForDD=0To72PRB(1)=4*0.000001*MCB*r*w^2/(π*d^2)RBB(DD)=Sqr(PTT(DD)^2+(PNN(DD)-PRB(1))^2)Zz=PNN(DD)-PRB(1)SETA1=Atn(Abs(PTT(DD)/Zz))SETA(DD)=(1-Sgn(Zz))*π/4+Sgn(Zz)*Sgn(PTT(DD))*SETA1NextDDForI2=1To36'STEP3'曲柄销值RBBb=0Forj=1To13k=I2+(j-7)Ifk36Thenb=k-36ElseIfk1Thenb=36+kElseb=kEndIfRBBb=RBBb+RBB(b)+RBB(36+b)NextjRB1(I2)=RBBbNextI2RB1(0)=RB1(36)fh=π*d^2/4:MK2=1.149:MK1=0.7833:MCA=1.395ForS1=0To72IfS1+1872ThenA5=S1-54ElseA5=S1+18IfS1+3672ThenA6=S1-36ElseA6=S1+36IfS1+5472ThenA7=S1-18ElseA7=S1+54PN2(S1)=PNN(A5):PN4(S1)=PNN(A6):PN3(S1)=PNN(A7)PT2(S1)=PTT(A5):PT3(S1)=PTT(A7):PT4(S1)=PTT(A6)NextS1PL=32.3/1000:MK2=1.49:MCA=1.359:MK1=0.7833'主轴颈值及轴承值RP=55.9/1000:MP=1.616Fori=0To720Step10MR=MK1+2*PL/r*MK2+MCApr=MR*r*w^2*10^-6/fhKP=MP*RP*w^2*0.00001/fhKP1=KP*92.5/124:A=i/10PT(0,A)=-PTT(A)/2:PN(0,A)=-PNN(A)/2PT(1,A)=-(PTT(A)-PT2(A))/2:PN(1,A)=-(PNN(A)-PN2(A))/2PT(2,A)=(PT2(A)+PT3(A))/2:PN(2,A)=(PN2(A)+PN3(A))/2PT(3,A)=-(PT3(A)-PT4(A))/2:PN(3,A)=-(PN3(A)-PN4(A))/2PZZ(1,A)=(Sqr((2*KP1+PNN(A)-pr)^2+PTT(A)))/2PZZ(2,A)=(Sqr((PTT(A)-PT2(A))^2+(PNN(A)-PN2(A))^2))/2PZZ(3,A)=(Sqr((PT2(A)+PT3(A))^2+(2*pr-4*KP1-(PN2(A)+PN3(A)))^2))/2PZZ(4,A)=(Sqr((PT3(A)-PT4(A))^2+(PN3(A)-PN4(A))^2))/2PZZ(4,A)=(Sqr((PN4(A)+2*KP1-pr)^2+PT4(A)^2))/2NextiMK1=0.7833:MK2=1.149'求主轴承负荷MR=MK1+2*32.2/57.5*MK2+MCBForIE=0To72ALFA=IE*π/18ROH(0,IE)=-PT(0,IE)*Cos(α)+(PN(0,IE)-pr)*Sin(α)ROV(0,IE)=PT(0,IE)*Sin(α)+(PN(0,IE)-pr)*Cos(α)ROH(1,IE)=-PT(1,IE)*Cos(α)+PN(1,IE)*Sin(α)ROV(1,IE)=PT(1,IE)*Sin(α)+PN(1,IE)*Cos(α)ROH(2,IE)=-PT(2,IE)*Cos(α)+(PN(2,IE)-pr)*Sin(α)ROV(2,IE)=PT(2,IE)*Sin(α)+(PN(2,IE)-pr)*Cos(α)ROH(3,IE)=-PT(3,IE)*Cos(α)+(PN(3,IE)-pr)*Sin(α)ROV(3,IE)=PT(3,IE)*Sin(α)+(PN(3,IE)-pr)*Cos(α)NextIEEndSubPrivateSubCommand10_Click()Picture3.ClsPicture3.DrawMode=13blndrawnow=FalsePicture3.Scale(0,400)-(650,-20)'绘制坐标系E=15:X=-20:Y=0Picture3.Line(120+X,200+Y)-(420+X,200+Y),SEPicture3.Line(270+X,30+Y)-(270+X,380+Y),SEPicture3.Line(267+X,40+Y)-(270+X,30+Y),SEPicture3.Line(270+X,30+Y)-(273+X,40+Y),SEPicture3.Line(410+X,197+Y)-(420+X,200+Y),SEPicture3.Line(420+X,200+Y)-(410+X,203+Y),SEForI1=1To71Step1Picture3.Line(250+ROH(0,I1-1)*20,200+ROV(0,I1-1)*20)-(250+ROH(0,I1)*20,200+ROV(0,I1)*20),vbGreenPicture3.Line(250+ROH(0,I1)*20,200+ROV(0,I1)*20)-(250+ROH(0,I1+1)*20,200+ROV(0,I1+1)*20),vbGreenNextI1Picture3.PSet(300,300),vbBluePicture3.Print第一,第五主轴承负荷图EndSubPrivateSubCommand11_Click()Picture3.ClsPicture3.DrawMode=13blndrawnow=FalsePicture3.Scale(0,600)-(650,-20)'绘制坐标系SE=15:X=-20:Y=0Picture3.Line(120+X,330+Y)-(420+X,330+Y),SEPicture3.Line(270+X,30+Y)-(270+X,400+Y),SEPicture3.Line(267+X,40+Y)-(270+X,30+Y),SEPicture3.Line(270+X,30+Y)-(273+X,40+Y),SEPicture3.Line(410+X,327+Y)-(420+X,330+Y),SEPicture3.Line(420+X,330+Y)-(410+X,333+Y),SEForI1=1To71Step1Picture3.Line(250+ROH(1,I1-1)*70,330+ROV(1,I1-1)*70)-(250+ROH(1,I1)*70,330+ROV(1,I1)*70),vbGreenPicture3.Line(250+ROH(1,I1)*70,330+ROV(1,I1)*70)-(250+ROH(1,I1+1)