某建筑地基的地层为淤泥质粘土

某建筑地基的地层为淤泥质粘土，固结系数为C h=C v=1.8×10-3cm2/s，粘聚力c=18kPa，内摩擦角φ=5︒，受压土层厚20m，袋装砂井直径d w=700mm，袋装砂井按等边三角形布置，间距l＝1.4m，深度H=20m，砂井底部为不透水层，砂井打穿受压土层。

预压荷载总压力P=100kPa，分两级等速加载，如图所示，计算地基堆载120天后，地层的平均固结度(不考虑竖井井阻和涂抹的影响)，并判断在建筑荷载(假设为条形基础，基础埋深为0m)为150kPa 时，地基中各点是否破坏？
C程序如下：
#include <stdio.h>
#include<math.h>
void main()
{
double a,b,c,f,d,g,Fn,D,l,C1,C2,B,m,n,w,x,temp1,temp2;
int t,H;
double T[20],r[10],Q[10],p[10],U[10],o[10];
double sum=0;
int q[10];
int i,j;
printf("请输入砂井直径、间距、固结系数、深度\n");
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=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("请输入加载次数\n");
scanf("%d",&t);
printf("请输入加载时间\n");
for(i=0;i<20;i++){
scanf("%lf",&T[i]);
if(i>=(2*t-1))
break;
}
printf("请输入加载数\n");
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++;
}
}
加载
次数
每次加载数及各次加载起止时间固结度
1 100
0.9483 0 20
2
30 70
0.9162 0 6 30 44
3
40 30 30
0.9098 0 8 30 36 50 56
4
30 30 30 10
0.8908 0 6 30 36 50 56 70 72
5
30 20 20 20 10
0.8523 0 6 30 34 50 54 70 74 90 92
结论：固结度随着加载次数的增大而减小。

60 40 加载速
率
固结度
0 5 30 33.3 12 0.9393
0 10 30 36.6 6 0.9361
0 15 30 40 4 0.9325
0 20 30 43.3 3 0.9288
0 25 30 46.6 2.4 0.9248 结论：固结度随着加载速率的减小而减小。

Cv Ch 固结度
0.0018 0.0001 0.3047 0.0011 0.8293 0.0021 0.9562 0.0031 0.9883 0.0041 0.9967
结论：固结度随着水平固结系数的增大而增大。

间距固结度
100 0.9968
130 0.9593
150 0.9045
170 0.8355
190 0.7631
结论：固结度随着井间距的增大而减小。

井径固结度
10 0.7976
60 0.9237
110 0.9631
160 0.9819
210 0.9914
结论：固结度随着井径的增大而增大。