公路橋梁荷載受力計算全套

公路橋梁荷載受力計算code------=1考慮抗扭因素=0不考慮抗扭因素It--------每根主梁的抗扭慣性矩II--------每根主梁的抗彎慣性矩RR--------人行道寬度BB--------行車道寬度AA--------主梁間距LL--------主梁跨徑N---------主梁個數(shù)C0--------車道數(shù)M---------荷載類型個數(shù)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("請輸入初始數(shù)據(jù)文件名:");scanf("%s",in);printf("\n請輸入結(jié)果數(shù)據(jù)文件名:");scanf("%s",out);if((fp1=fopen(in,"r"))==NULL){printf("\n\n打開文件%s失敗\n",in);exit(0);}if((fp2=fopen(out,"w"))==NULL){printf("創(chuàng)建文件%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("輸入數(shù)據(jù)有誤!\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*