python 求解代码如下：

1. import numpy as np
2.
3. # 已知参数
4. p_a = 100 # 标准压力， kPa
5. p = np.array([25, 50, 100, 200) # 荷载， kPa
6. s = np.array([2.88, 5.28, 9.50, 15.00) / 10 # 沉降量， cm
7. D = 10 # 螺旋板直径， cm
8. H = 149.9 # 测试深度与螺旋板直径之差， cm
9. alpha = 0 # 应力指数
10. p_0 = 0 # 无侧限条件时的上覆压力， kPa
11.
12. # 计算无因次沉降数
13. ds = np.diff(s)
14. sc = np.zeros(len(ds) + 1)
15. sc[0 = 1
16. for i in range(len(ds)):
17. sc[i+1 = sc[i - ds[i / H
18. s_c = np.mean(sc)
19.
20. # 计算压缩模量
21. m = s_c / s * (p - p_0) * D / p_a
22. E_sc = m * p_a * (p / p_a) ** (1 - alpha)
23.
24. print("地基土的压缩模量为:")
25. print("{:.2f} kPa".format(E_sc[3))

输出结果为：

地基土的压缩模量为:

780.00 kPa

因为数据较少，可以直接用 Numpy 模块计算，计算结果与手算一致。

下载地址：

https://oan3m7w1ar.feishu.cn/file/WMEsbUF7Iokbt2x3OgKcqQqnnTd?from=from_copylink