公路桥梁荷载受力计算全套.docx
公路桥梁荷载受力计算code=1考虑抗扭因素=0不考虑抗扭因素It每根主梁的抗扭惯性矩II每根主梁的抗弯惯性矩RR人行道宽度BB行车道宽度AA主梁间距1.L主梁跨径N主梁个数CO车道数M荷载类型个数1=<M<=4MI存放各种类型荷载序号I-人群荷载2汽车荷载3-挂车荷载4一一履带车荷载单位采用国际单位*include,stdio.hninclude,stdlib.h#include,alloc.h,'#includenprocess.h"voidcomputeMl(intcodezdoubleItzdouble!,doubleRRzdoubleBBzdoubleAA,doubleLLJntNJntCOJntMJntMI4rdouble*MM);intmain()(double*MM,It,11,RR,BB,AA,LL;intiJcode,k,N,M,C0,MI4;FILE*fplz*fp2;charin20zout20zstr30;Primfe请输入初始数据文件名：")；scanf(',%s,zin);printf("n请输入结果数据文件名：");scanf("%s,zout);if(fpl=fopen(inzr,)=NULL)printf(,nn打开文件s失败nn,in);exit(0);if(fp2=fopen(outr,w)=NLL)printf(仓U建文件s失败n"zout);exit(0);fscanf(fpl,"%s%s%s%s%s"zstr,str,str,str,str);fscanf(fpL"d%lf%lf%lf%lf”,&code,&It,&n,&RR,&BB);fscanf(fpl,"%s%s%s%s%snzstr,strzstr,str,str);fscanf(fpl,”lf%lf%d%d%d”,&AA,&LL,&N,&CO,&M);for(i=0;i<4;i+)MIi=0;if(M>4&&Ml)printf(”nM值有误!n)exit(0)if(N>4&&M<l)printf(”nN值有误!n)exit(0)for(i=0;i<M;i+)fscanf(fpl%dMIi);k=(N+l)2;MM=malloc(k*M*sizeof(double);if(MM=NULL)printf(nK值太大了n");retum;for(i=0;i<k*M;i+)MMi=0;ComputeMl(CodezItzIIzRRzBBzAAzLLzNzCOzMzMIzMM);for(i=0;i<M;i+)switch(MIi)case1:fprintf(fp2z人群荷载)break;case 2:fprintf(fp2/汽车荷载)break;case 3: fprintf(fp2,”挂车荷载)break;case 4: fprintf(fp2zn履带车荷载)break;)fprintf(fp2,nn");for(i=0;i<k;i+)fpntf(fp2,%4d,J+l);for(j=0;j<M;j+)fprintf(fp2,"%15lf",MMi*M+j);fprintf(fp2znn");)fclose(fpl);fclose(fp2);voidcomputeMl(intcodezdoubleItzdoubleIIzdoubleRRdoubIeBBzdoubleAAzdoubleLLJntNJntCOJntMJntMI4zdouble*MM)(intkjj,z;doubleB,MO,MQMCtemP=O,E,ELTC,CC,TG,Tl,T2,tt;k=(int)(N+l)2;if(N<=:LIlnC=OIlAA<=O)printf(,输入数据有误!n");return;if(code=l)temp=N;temp=temp*(temp*temp-1.0);temp=AA*AA*temp;temp=N*LL*LLtemp;B=l+0.425*temp*ItII;B=1B;elseB=1.0;tt=N;E=(N+l)/2.0;tt=tt*(tt*tt-1.0)12.0;tt=1.0tt;for(i=0;i<k;i+)T1=1.0/N+tt*(E-i-1.0)*(E-1.0)*B;T2=1.0N+tt*(E-i-1.0)*(E-N)*B;TC=(T1+T2)2.O;TG=(T1-T2)(N-1)*AA);forQ=0;j<M;j+)(switch(MIj)El=(RR+BB)2.0;if(TC>TG*El)temp=2*TC;elsetemp=TC+TG*El;break;case2:CC=1.0;temp=0;for(z=l;z<=CO;z+)if(z=3)CC=0.8;if(z=4)CC=0.7;El=BB2.0-1.55*(Z-I)-1.4;MO=TC+TG*E1;MC=CC*z*M0;if(MC>temp)temp=MC;)break;temp=TC+TG*(BB-4.7)2.0;break;case4:temp=TC+TG*(BB-4.5)2.0;break;)MMi*M+j=temp;)
