公路桥梁荷载受力计算全套

公路桥梁荷载受力计算code------=1考虑抗扭因素=0不考虑抗扭因素It--------每根主梁的抗扭惯性矩II--------每根主梁的抗弯惯性矩RR--------人行道宽度BB--------行车道宽度AA--------主梁间距LL--------主梁跨径N---------主梁个数C0--------车道数M---------荷载类型个数1=<M<=4MI--------存放各种类型荷载序号1----人群荷载2----汽车荷载3----挂车荷载4----履带车荷载单位采用国际单位#include"stdio.h"#include"stdlib.h"#include"alloc.h"#include"process.h"voidcomputeM1(intcode,doubleIt,doubleII,doubleRR,doubleBB,doubleAA,doubleLL,intN,intC0,intM,intMI[4],double*MM);intmain(){double*MM,It,II,RR,BB,AA,LL;inti,j,code,k,N,M,C0,MI[4];FILE*fp1,*fp2;charin[20],out[20],str[30];printf("请输入初始数据文件名:");scanf("%s",in);printf("\n请输入结果数据文件名:");scanf("%s",out);if((fp1=fopen(in,"r"))==NULL){printf("\n\n打开文件%s失败\n",in);exit(0);}if((fp2=fopen(out,"w"))==NULL){printf("创建文件%s失败\n",out);exit(0);}fscanf(fp1,"%s%s%s%s%s",str,str,str,str,str);fscanf(fp1,"%d%lf%lf%lf%lf",&code,&It,&II,&RR,&BB);fscanf(fp1,"%s%s%s%s%s",str,str,str,str,str);fscanf(fp1,"%lf%lf%d%d%d",&AA,&LL,&N,&C0,&M);for(i=0;i<4;i++)MI[i]=0;if(M>4&&M<1){printf("\nM值有误!\n"),exit(0);}if(N>4&&M<1){printf("\nN值有误!\n"),exit(0);}for(i=0;i<M;i++)fscanf(fp1,"%d",&MI[i]);k=(N+1)/2;MM=malloc(k*M*sizeof(double));if(MM==NULL){printf("\nK值太大了\n");return;}for(i=0;i<k*M;i++)MM[i]=0;computeM1(code,It,II,RR,BB,AA,LL,N,C0,M,MI,MM);for(i=0;i<M;i++){switch(MI[i]){case1:fprintf(fp2,"人群荷载");break;case2:fprintf(fp2,"汽车荷载");break;case3:fprintf(fp2,"挂车荷载");break;case4:fprintf(fp2,"履带车荷载");break;}}fprintf(fp2,"\n");for(i=0;i<k;i++){fprintf(fp2,"%4d",i+1);for(j=0;j<M;j++)fprintf(fp2,"%15lf",MM[i*M+j]);fprintf(fp2,"\n");}fclose(fp1);fclose(fp2);}voidcomputeM1(intcode,doubleIt,doubleII,doubleRR,doubleBB,doubleAA,doubleLL,intN,intC0,intM,intMI[4],double*MM){intk,i,j,z;doubleB,M0,MQ,MC,temp=0,E,E1,TC,CC,TG,T1,T2,tt;k=(int)(N+1)/2;if(N<=1||II<=0||AA<=0){printf("输入数据有误!\n");return;}if(code==1){temp=N;temp=temp*(temp*temp-1.0);temp=AA*AA*temp;temp=N*LL*LL/temp;B=1+0.425*temp*It/II;B=1/B;}elseB=1.0;tt=N;E=(N+1)/2.0;tt=tt*(tt*tt-1.0)/12.0;tt=1.0/tt;for(i=0;i<k;i++){T1=1.0/N+tt*(E-i-1.0)*(E-1.0)*B;T2=1.0/N+tt*(E-i-1.0)*(E-N)*B;TC=(T1+T2)/2.0;TG=(T1-T2)/((N-1)*AA);for(j=0;j<M;j++){switch(MI[j]){case1:E1=(RR+BB)/2.0;if(TC>TG*E1)temp=2*TC;elsetemp=TC+TG*E1;break;case2:CC=1.0;temp=0;for(z=1;z<=C0;z++){if(z==3)CC=0.8;if(z==4)CC=0.7;E1=BB/2.0-1.55*(z-1)-1.4;M0=TC+TG*E1;MC=CC*z*