2、>#includevoidmain(){doublea,b,c,f,d,g,Fn,D,l,C1,C2,B,m,n,w,x,temp1,temp2;intt,H;doubleT[20],r[10],Q[10],p[10],U[10],o[10];doublesum=0;intq[10];inti,j;printf("请输入砂井直径、间距、固结系数、深度");scanf("%lf,%lf,%lf,%lf,%d",&w,&l,&C1,&C2,&H);D=1.05*l;x=D/w;a=x*x;b=x*x-1;c
3、=a/b;f=log(x);d=c*f;g=(3*a-1)/(4*a);Fn=d-g;temp1=(3.14*3.14*C1)/(4*H*H);temp2=(8*C2)/(Fn*D*D);B=(temp1+temp2)*86400;m=exp(-120*B);n=0.81/B*m;printf("请输入加载次数");scanf("%d",&t);printf("请输入加载时间");for(i=0;i<20;i++){scanf("%lf",&T[i]);if(i>=(2*t-1))break;}printf("请输
4、入加载数");for(j=0;j<10;j++){scanf("%d",&q[j]);if(j>=(t-1))break;}for(i=0,j=0;i<20,j<10;i++,j++){if(j>=t&&i>=2*t)break;r[j]=T[i+1]-T[i];Q[j]=q[j]/r[j];o[j]=exp(B*r[j]);p[j]=exp(B*T[i])*(o[j]-1);U[j]=Q[j]/100*(r[j]-n*p[j]);sum=sum+U[j];printf("%lf",sum);i++;}}加载次数每次加
5、载数及各次加载起止时间固结度11000.9483020230700.916206304434030300.909808303650564303030100.890806303650567072530202020100.8523063034505470749092结论:固结度随着加载次数的增大而减小。6040加载速率固结度053033.3120.93930103036.660.9361015304040.93250203043.330.92880253046.62.40.9248结论:固结度随着加载速率的减小而减小。CvCh
6、固结度0.00180.00010.30470.00110.82930.00210.95620.00310.98830.00410.9967结论:固结度随着水平固结系数的增大而增大。间距固结度1000.99681300.95931500.90451700.83551900.7631结论:固结度随着井间距的增大而减小。井径固结度100.7976600.92371100.96311600.98192100.9914结论:固结度随着井径的增大而增大。