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