Tổng quan luận án
Văn bản thuộc nhóm công trình khoa học chuyên ngành Kỹ thuật Xây dựng Công trình Giao thông (Địa kỹ thuật công trình), được thực hiện tại Trường Đại học Giao thông Vận tải (Hà Nội, năm 2023). Bản văn bản được cung cấp là tập phụ lục các chương trình tính toán số phục vụ cho luận án tiến sĩ kỹ thuật, tập trung giải quyết bài toán phân tích ổn định mái dốc và ổn định trượt sâu của mố cầu đặt trên móng nông.
Tính cấp thiết và khoảng trống nghiên cứu được xác định xuất phát từ yêu cầu kiểm soát an toàn cho mố cầu móng nông chịu tác động đồng thời của tĩnh tải kết cấu, áp lực đất đắp sau lưng mố với chiều cao thay đổi, các tải trọng thi công và khai thác bất lợi, cũng như các yếu tố tự nhiên như mực nước ngầm và lực quán tính động đất. Trong thực tế tính toán thiết kế công trình giao thông, việc xác định hệ số an toàn trượt sâu nhỏ nhất ($F_{s\min}$) và vị trí mặt trượt nguy hiểm nhất thường gặp khó khăn do cấu trúc hình học phức tạp của hệ tương tác mố - đất và tính phi tuyến của các phương trình cân bằng giới hạn.
Mục tiêu nghiên cứu của công trình bao gồm hai nội dung chính:
- Xây dựng và thiết lập mô hình ma trận giải tích dựa trên phương pháp cân bằng giới hạn tổng quát (GLEM - Generalized Limit Equilibrium Method) nhằm xác định trạng thái ứng suất, lực tương tác giữa các phân thỏi và hệ số an toàn ổn định của mái đất dính đồng nhất và hệ mố cầu trên móng nông.
- Ứng dụng công cụ tối ưu hóa Solver trong môi trường bảng tính Microsoft Excel để tự động hóa quá trình tìm kiếm mặt trượt nguy hiểm nhất tương ứng với hệ số an toàn trượt sâu nhỏ nhất trong các giai đoạn thi công và khai thác.
Văn bản đặt ra các câu hỏi nghiên cứu cụ thể:
- Phương pháp cân bằng giới hạn tổng quát GLEM giải quyết bài toán cân bằng tĩnh học lực pháp tuyến và tiếp tuyến giữa các thỏi trượt đất và khối mố cầu như thế nào thông qua hệ phương trình ma trận $[A]{K} = {B}$?
- Chiều cao đất đắp thay đổi ($H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$) cùng với sự tham gia của các tổ hợp tải trọng thi công (xe 3 trục, xe 2 trục, xe lu), hoạt tải khai thác và tải trọng động đất tác động ra sao đến vị trí mặt trượt tối ưu và giá trị hệ số an toàn trượt sâu của mố cầu?
Đối tượng nghiên cứu là ổn định tổng thể của mái dốc đất dính đồng nhất và ổn định trượt sâu của mố cầu trên móng nông. Phạm vi nghiên cứu về không gian là mô hình hình học phẳng 2D của mái dốc và hệ mố - đất đắp; phạm vi về điều kiện làm việc khảo sát hai giai đoạn kỹ thuật: giai đoạn thi công (xét chiều cao đắp $H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$, tải trọng xe 3 trục, xe 2 trục, xe lu, động đất) và giai đoạn khai thác (xét hoạt tải khai thác, mực nước ngầm, động đất theo các trạng thái giới hạn cường độ 2 và cường độ 3).
Tổng quan tài liệu và vị trí của luận án
Trong phạm vi văn bản phụ lục tính toán được cung cấp, cơ sở tài liệu được xây dựng trực tiếp trên nền tảng phương pháp cân bằng giới hạn tổng quát GLEM (Generalized Limit Equilibrium Method) trong cơ học đất và địa kỹ thuật công trình. Phương pháp GLEM là sự mở rộng của các phương pháp phân chia thỏi cổ điển, cho phép thiết lập đầy đủ các phương trình cân bằng tĩnh học đối với từng phân thỏi cũng như xét đến lực tương tác tiếp xúc trên các mặt biên giữa các thỏi liền kề.
Bên cạnh lý thuyết GLEM, luận án tích hợp triết lý thiết kế hệ số tải trọng và sức kháng (LRFD) theo định dạng tiêu chuẩn thiết kế cầu đường bộ hiện hành, thể hiện qua việc phân loại và áp dụng các tổ hợp tải trọng:
- Trạng thái giới hạn cường độ I (Strength I Limit State): áp dụng cho giai đoạn thi công khi xét tĩnh tải kết cấu mố và tĩnh tải đất đắp trong lòng mố với các hệ số tải trọng quy định.
- Trạng thái giới hạn cường độ 2 và cường độ 3 (Strength II, Strength III Limit States): áp dụng cho giai đoạn khai thác để đánh giá ảnh hưởng của hoạt tải xe cơ giới, mực nước ngầm bão hòa và lực quán tính động đất.
Khoảng trống nghiên cứu mà tác giả lựa chọn giải quyết trong phần tính toán này là chuyển hóa toàn bộ hệ phương trình vi phân và phương trình cân bằng tĩnh học phức tạp của GLEM thành dạng đại số ma trận cấp cao $[A]{K} = {B}$, từ đó kết nối trực tiếp với công cụ Solver trong Excel để tìm nghiệm tối ưu mặt trượt và hệ số an toàn mà không cần phụ thuộc vào các phần mềm địa kỹ thuật chuyên dụng đóng kín.
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng khung lý thuyết cân bằng giới hạn tổng quát (GLEM) kết hợp định luật chống cắt Coulomb. Mô hình phân chia khối trượt giả định thành $n$ phân thỏi đất thẳng đứng và một khối mố kết cấu ký hiệu là khối $M$. Mỗi phân thỏi $i$ chịu tác dụng của trọng lượng bản thân $W_i$, phản lực pháp tuyến đáy $N_i$, lực cắt tiếp tuyến đáy $T_i$, cùng các lực tương tác tại hai mặt biên tiếp xúc giữa thỏi thứ $i$ và thỏi thứ $i+1$ gồm lực pháp tuyến $H_i, H_{i+1}$ và lực tiếp tuyến $V_i, V_{i+1}$.
Phương trình cân bằng tĩnh học của phân thỏi thứ $i$ được thiết lập bằng cách chiếu toàn bộ hệ lực lên phương pháp tuyến và phương tiếp tuyến của mặt đáy khối trượt:
- Chiếu lên phương pháp tuyến đáy:
$$N_i - W_i \cos \alpha_i - H_i \cos(\beta_i - \alpha_i) + V_i \sin(\beta_i - \alpha_i) + H_{i+1} \cos(\beta_{i+1} - \alpha_i) - V_{i+1} \sin(\beta_{i+1} - \alpha_i) = 0$$
- Chiếu lên phương tiếp tuyến đáy:
$$T_i - W_i \sin \alpha_i + H_i \sin(\beta_i - \alpha_i) + V_i \cos(\beta_i - \alpha_i) - H_{i+1} \sin(\beta_{i+1} - \alpha_i) - V_{i+1} \cos(\beta_{i+1} - \alpha_i) = 0$$
Đối với khối mố $M$, phương trình cân bằng được bổ sung thêm thành phần tải trọng thẳng đứng của bản thân kết cấu mố $Q_M$ và tĩnh tải đất đắp trong lòng mố $Q_{soil}$. Điều kiện trượt trên các mặt phẳng đáy và mặt biên giữa các khối tuân theo biểu thức sức kháng cắt Mohr - Coulomb được chuẩn hóa qua hệ số an toàn $F_s$:
- Tại đáy khối thứ $i$: $T_i = \frac{k_i (c_1 S_i B + N_i \tan \varphi_1)}{F_s}$
- Tại mặt biên giữa khối thứ $i$ và $i+1$: $V_i = \frac{m_i (c_1 R_i B + H_i \tan \varphi_1)}{F_{si}}$
Toàn bộ hệ thức cân bằng và điều kiện biên được cấu trúc thành phương trình ma trận tuyến tính dạng:
$$[A] \cdot {K} = {B}$$
Trong đó:
- ${K}$ là véc-tơ chứa 24 ẩn số lực: ${K} = [N_1, T_1, H_2, V_2, N_2, T_2, H_3, V_3, N_3, T_3, H_4, V_4, N_4, T_4, H_5, V_5, N_5, T_5, H_6, V_6, N_6, T_6, H_7, V_7]^T$.
- $[A]$ là ma trận vuông hệ số kích thước $24 \times 24$ chứa các hàm lượng giác của góc nghiêng mặt đáy $\alpha_i$, góc nghiêng mặt biên $\beta_i$, góc ma sát trong $\varphi$, chiều dài đáy $S_i$, chiều dài biên $R_i$ và các hệ số điều hướng $k, m$.
- ${B}$ là véc-tơ tải trọng ngoại lực phụ thuộc vào trọng lượng khối trượt $W_i$, tải trọng mố $Q_M$, tải trọng đất đắp $Q_{soil}$, lực dính đơn vị $c_1, c_2$ và bề rộng tính toán $B$.
Hệ thống ký hiệu quy ước trong thuật toán được tổng hợp trong các bảng danh mục ký hiệu của tài liệu:
| Ký hiệu |
Đơn vị |
Bản chất đại lượng và mô tả |
| $Q_M$ |
kN |
Tĩnh tải mố cầu |
| $Q_{soil}$ |
kN |
Tĩnh tải đất đắp trong lòng mố |
| $W, W_M$ |
kN |
Trọng lượng khối trượt đất và trọng lượng khối trượt mố |
| $W'_i, W''_i$ |
kN |
Trọng lượng phần đất bị bão hòa (ngập nước) và không bão hòa |
| $N, T$ |
kN |
Lực pháp tuyến và lực tiếp tuyến tác dụng lên mặt phẳng đáy khối trượt |
| $H, V$ |
kN |
Lực pháp tuyến và lực tiếp tuyến tác dụng lên mặt phẳng giữa hai khối |
| $U, U'$ |
kN |
Tổng lực đẩy nổi trên một đơn vị chiều dài theo phương pháp tuyến |
| $Q_{DL}, DD$ |
kN |
Tĩnh tải kết cấu nhịp và tải trọng kéo xuống do ma sát âm |
| $DC, DW$ |
kN |
Tĩnh tải bộ phận kết cấu, phụ trợ và tĩnh tải lớp phủ mặt cầu |
| $v, h$ |
$\text{m/s}^2$ |
Gia tốc quán tính động đất theo phương thẳng đứng và nằm ngang |
| $BR, FR, IM$ |
$\text{kN, %}$ |
Lực hãm xe, lực ma sát và độ gia tăng lực do xung kích của xe |
| $LS, PL$ |
kN |
Hoạt tải chất thêm và tải trọng người đi bộ |
| $WL, WS$ |
kN |
Tải trọng gió tác dụng lên hoạt tải và tải trọng gió tác dụng lên kết cấu |
| $F_s, F_{si}$ |
- |
Hệ số an toàn trên mặt phẳng đáy và mặt phẳng giữa các khối |
| $\alpha_i, \beta_i, \theta_e$ |
Độ |
Góc nghiêng mặt đáy, góc nghiêng mặt biên và mặt biên phần ngập nước |
| $S_i, R_i$ |
m |
Chiều dài mặt đáy và chiều dài mặt phẳng giữa các khối |
| $k, m$ |
- |
Thông số xác định hướng của lực cắt đáy $T$ và lực cắt biên $V$ |
| $c_1, \varphi_1, \gamma_1$ |
$\text{kN/m}^2, \text{kN/m}^3$ |
Lực dính, góc ma sát trong và trọng lượng riêng của lớp đất tự nhiên |
| $c_2, \varphi_2, \gamma_2$ |
$\text{kN/m}^2, \text{kN/m}^3$ |
Lực dính, góc ma sát trong và trọng lượng riêng của lớp đất đắp |
| $g, \gamma_c$ |
$\text{m/s}^2, \text{kN/m}^3$ |
Gia tốc trọng trường ($g = 10\text{ m/s}^2$) và trọng lượng riêng bê tông |
| $H_đ$ |
m |
Chiều cao đất đắp sau lưng mố |
Phương pháp phân tích số được thực hiện thông qua công cụ Solver trong Microsoft Excel. Hàm mục tiêu được thiết lập là cực tiểu hóa hệ số an toàn ($F_s \rightarrow \min$), với các biến số thay đổi là tọa độ hình học $(X_i, Y_i)$ của các điểm nút mặt trượt và bán kính $R_i$, thỏa mãn hệ ràng buộc cân bằng ma trận $[A]{K} = {B}$ và các điều kiện biên hình học của mái dốc.
Nội dung chính theo từng chương
Thí dụ tính ổn định mái dốc bằng phương pháp giới hạn tổng quát GLEM khi mái dốc là mái đất dính đồng nhất (trường hợp trượt qua chân dốc)
Nội dung phần này trình bày bài toán phân tích ổn định mái dốc đất dính đồng nhất với cung trượt đi qua chân dốc. Khối trượt được chia thành 6 phân thỏi chính với 7 mặt biên tiếp xúc. Sau khi thiết lập ma trận cân bằng $24 \times 24$ và thực hiện tối ưu hóa bằng Solver trong Excel, kết quả hình học và cơ học của mặt trượt được xác định chi tiết trong bảng thông số:
| STT |
$\beta_i$ (Độ) |
$R_i$ (m) |
$XP_i$ (m) |
$YP_i$ (m) |
$X_i$ (m) |
$Y_i$ (m) |
$\alpha_i$ (Độ) |
Diện tích $A_i$ ($\text{m}^2$) |
Trọng lượng $W_i$ (kN) |
Chiều dài đáy $S_i$ (m) |
| 1 |
33,69 |
2,42 |
-7,14 |
-4,76 |
-9,15 |
-6,10 |
-10,71 |
1,24 |
24,74 |
1,46 |
| 2 |
70,47 |
1,71 |
-7,14 |
-4,76 |
-7,71 |
-6,37 |
10,81 |
3,74 |
74,71 |
2,23 |
| 3 |
84,74 |
2,43 |
-5,30 |
-3,53 |
-5,52 |
-5,95 |
28,04 |
6,66 |
133,12 |
3,65 |
| 4 |
110,77 |
2,31 |
-3,12 |
-2,08 |
-2,30 |
-4,24 |
41,55 |
5,31 |
106,17 |
3,00 |
| 5 |
130,26 |
1,85 |
-1,25 |
-0,83 |
-0,05 |
-2,25 |
63,98 |
2,31 |
46,22 |
1,98 |
| 6 |
150,38 |
0,94 |
0,00 |
0,00 |
0,82 |
-0,46 |
90,00 |
0,19 |
3,79 |
0,46 |
| 7 |
180,00 |
0,82 |
0,00 |
0,00 |
0,82 |
0,00 |
- |
- |
- |
- |
Kết quả tính toán hệ số an toàn sau tối ưu hóa:
- Khi xét hệ số tương tác biên $FSI = \infty$ (tương đương điều kiện biên tự do), hệ số an toàn nhỏ nhất hội tụ đạt giá trị $FS_{\min_tu} = 1,380$.
- Khi xét trường hợp phân bố ứng suất trung gian với giả định $FSI(i) = FS(i) = 1,7$, hệ số an toàn sau tối ưu hóa đạt $FS_{\min_tu} = 1,510$.
Thí dụ tính ổn định mái dốc bằng phương pháp giới hạn tổng quát GLEM khi mái dốc là mái đất dính đồng nhất (trường hợp trượt sâu)
Trường hợp trượt sâu phân tích cung trượt phát triển sâu xuống phía dưới chân dốc, vượt qua cao độ chân dốc tự nhiên. Quá trình giải hệ phương trình ma trận $[A]{K} = {B}$ kết hợp hàm Solver mang lại bảng thông số hình học và tải trọng phân thỏi sau tối ưu:
| STT |
$\beta_i$ (Độ) |
$R_i$ (m) |
$XP_i$ (m) |
$YP_i$ (m) |
$X_i$ (m) |
$Y_i$ (m) |
$\alpha_i$ (Độ) |
Diện tích $A_i$ ($\text{m}^2$) |
Trọng lượng $W_i$ (kN) |
Chiều dài đáy $S_i$ (m) |
| 1 |
0,00 |
1,85 |
-9,15 |
-6,10 |
-11,00 |
-6,10 |
-20,06 |
0,52 |
10,35 |
1,63 |
| 2 |
60,45 |
0,64 |
-9,15 |
-6,10 |
-9,47 |
-6,66 |
-0,09 |
2,74 |
54,76 |
2,51 |
| 3 |
95,46 |
1,91 |
-7,14 |
-4,76 |
-6,96 |
-6,66 |
15,72 |
5,29 |
105,73 |
2,67 |
| 4 |
110,76 |
2,57 |
-5,30 |
-3,53 |
-4,39 |
-5,94 |
29,92 |
8,06 |
161,20 |
3,52 |
| 5 |
130,25 |
2,76 |
-3,12 |
-2,08 |
-1,34 |
-4,18 |
49,96 |
8,78 |
175,65 |
4,37 |
| 6 |
150,38 |
1,70 |
0,00 |
0,00 |
1,47 |
-0,84 |
49,96 |
0,91 |
18,26 |
1,09 |
| 7 |
180,00 |
2,18 |
0,00 |
0,00 |
2,18 |
0,00 |
- |
- |
- |
- |
Kết quả hệ số an toàn ổn định trượt sâu:
- Ở trường hợp $FSI = \infty$, hệ số an toàn nhỏ nhất sau tối ưu hóa đạt $FS_{\min_tu} = 1,469$.
- Ở trường hợp hiệu chỉnh ứng suất biên $FSI(i) = FS(i) = 1,7$, hệ số an toàn sau tối ưu hóa đạt $FS_{\min_tu} = 1,654$.
Thí dụ tính ổn định trượt sâu mố cầu trên móng nông trong giai đoạn thi công khi xét áp lực đất đắp sau lưng mố với chiều cao đất đắp thay đổi ($H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$)
Phần này khảo sát bài toán thực tế của mố cầu trên móng nông trong giai đoạn thi công xây dựng. Tĩnh tải các bộ phận kết cấu mố cầu được xác định cụ thể theo hồ sơ thiết kế:
- Tĩnh tải bệ mố ($Q_{M1}$): $2829,75\text{ kN}$
- Tĩnh tải tường thân ($Q_{M2}$): $2263,80\text{ kN}$
- Tĩnh tải tường cánh ($Q_{M3}$): $902,40\text{ kN}$
- Tĩnh tải tường đỉnh ($Q_{M4}$): $294,43\text{ kN}$
- Tĩnh tải vai kê ($Q_{M5}$): $36,38\text{ kN}$
- Tổng tĩnh tải kết cấu mố ($Q_M$): $6326,76\text{ kN}$ (tổng trọng lượng chưa nhân hệ số tải trọng). Khi tính toán tổ hợp tải trọng mố cầu $Q_M$, tổng lực thẳng đứng chưa kể đất đắp là $7908,45\text{ kN}$.
Sự biến thiên chiều cao đất đắp sau lưng mố làm thay đổi tĩnh tải đất đắp trong lòng mố $Q_{soil}$ và tổng tĩnh tải khối mố $W_M$:
- Với trường hợp $H_đ = 2,05\text{ m}$: Tĩnh tải đất đắp $Q_{soil} = 6096,87\text{ kN}$, tổng tĩnh tải phần mố $W_M = 12423,63\text{ kN}$. Theo tổ hợp Trạng thái giới hạn Cường độ I (hệ số tải trọng mố $\gamma_P = 1,25$, hệ số đất đắp $\gamma_P = 1,35$), tổng lực thẳng đứng tính toán đạt $16139,22\text{ kN}$.
- Với trường hợp $H_đ = 4,09\text{ m}$: Tĩnh tải đất đắp $Q_{soil} = 3223,94\text{ kN}$, tổng tĩnh tải phần mố $W_M = 9550,69\text{ kN}$. Theo tổ hợp Trạng thái giới hạn Cường độ I, tổng lực đứng tính toán đạt $12260,76\text{ kN}$.
- Với trường hợp $H_đ = 8,19\text{ m}$: Tĩnh tải đất đắp $Q_{soil} = 6096,87\text{ kN}$, tổng tĩnh tải phần mố $W_M = 12423,63\text{ kN}$, tổng lực đứng tính toán theo Cường độ I đạt $16139,22\text{ kN}$.
Bảng thông số hình học mặt trượt dự kiến ban đầu cho trường hợp $H_đ = 2,05\text{ m}$ và $H_đ = 4,09\text{ m}$:
| STT |
$H_đ = 2,05\text{ m}$: $\beta_i (^\circ)$ |
$H_đ = 2,05\text{ m}$: $W_i$ (kN) |
$H_đ = 2,05\text{ m}$: $S_i$ (m) |
$H_đ = 4,09\text{ m}$: $\beta_i (^\circ)$ |
$H_đ = 4,09\text{ m}$: $W_i$ (kN) |
$H_đ = 4,09\text{ m}$: $S_i$ (m) |
| 1 |
0,00 |
1724,22 |
2,31 |
0,00 |
939,48 |
1,72 |
| 2 |
15,00 |
2351,25 |
3,78 |
12,36 |
1352,71 |
3,40 |
| 3 |
45,00 |
5309,22 |
9,97 |
45,00 |
2881,55 |
8,64 |
| 4 |
130,00 |
11782,10 |
13,15 |
145,00 |
5839,00 |
12,57 |
| 5 |
128,00 |
7250,69 |
7,52 |
150,00 |
2568,07 |
3,45 |
| 6 |
165,00 |
3395,73 |
3,37 |
171,00 |
1320,73 |
1,69 |
| 7 |
180,00 |
- |
- |
180,00 |
- |
- |
Hệ thống ma trận hệ số $[A]$ ($24 \times 24$) và véc-tơ ngoại lực ${B}$ được thiết lập riêng cho từng trường hợp chiều cao đắp, làm đầu vào để công cụ Solver giải lặp tìm nghiệm $F_{s\min}$.
Các trường hợp mở rộng trong giai đoạn thi công và giai đoạn khai thác
Mục lục văn bản chỉ rõ hệ thống các chương trình tính toán mở rộng đánh giá toàn diện các yếu tố tác động khác:
- Giai đoạn thi công kết hợp hoạt tải thiết bị thi công:
- Tác động của xe 3 trục thiết kế tương ứng với 3 cấp chiều cao đất đắp ($H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$).
- Tác động của xe 2 trục thiết kế tương ứng với 3 cấp chiều cao đất đắp ($H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$).
- Tác động của tải trọng xe lu phục vụ đầm nén đất sau mố tương ứng với 3 cấp chiều cao đất đắp ($H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$).
- Giai đoạn thi công khi có động đất xảy ra: tính toán ổn định trượt sâu mố cầu dưới tác động đồng thời của áp lực đất đắp ($H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$) và gia tốc quán tính động đất theo phương thẳng đứng ($v$) và nằm ngang ($h$).
- Giai đoạn khai thác công trình:
- Xét hoạt tải khai thác: tính toán xác định $F_{s\min}$ theo hai trạng thái giới hạn là Trạng thái giới hạn cường độ 2 và Trạng thái giới hạn cường độ 3.
- Xét tác động của mực nước ngầm: đánh giá lực đẩy nổi $U, U'$ và sự suy giảm sức kháng cắt của các khối đất bão hòa nước ($W'_i$) theo Trạng thái giới hạn cường độ 2 và cường độ 3.
- Xét tổ hợp tải trọng khi có động đất xảy ra trong giai đoạn khai thác: tính toán hệ số an toàn nhỏ nhất $F_{s\min}$ theo Trạng thái giới hạn cường độ 2 và cường độ 3.
Kết quả và những đóng góp mới
Dựa trên các khối lượng tính toán được trình bày trong tập phụ lục, công trình đạt được các kết quả và đóng góp kỹ thuật:
- Về mặt phương pháp luận tính toán: Luận án đã cụ thể hóa phương pháp cân bằng giới hạn tổng quát GLEM thành hệ phương trình ma trận tường minh $[A]{K} = {B}$ có kích thước $24 \times 24$ để giải bài toán trượt sâu mố cầu móng nông. Mô hình ma trận này xử lý đồng thời liên kết cơ học giữa khối mố kết cấu cứng và các phân thỏi đất đắp, đất nền xung quanh.
- Về mặt công nghệ tính toán: Ứng dụng giải pháp tính toán trên bảng tính Microsoft Excel kết hợp công cụ tối ưu hóa Solver. Phương pháp này cho phép tìm nghiệm hệ số an toàn nhỏ nhất ($F_{s\min}$) và bề mặt trượt nguy hiểm nhất trực tiếp mà không cần lập trình mã nguồn phức tạp trên các ngôn ngữ bậc cao hay sử dụng phần mềm thương mại đắt tiền.
- Về mặt dữ liệu thực nghiệm số: Cung cấp tập số liệu chi tiết về tọa độ mặt trượt, diện tích, trọng lượng phân thỏi và hệ số an toàn cho các kịch bản điển hình: mái dốc đồng nhất trượt chân dốc ($F_{s\min} = 1,380 - 1,510$), trượt sâu ($F_{s\min} = 1,469 - 1,654$), và hệ mố cầu chịu tải trọng phân tầng với các chiều cao đắp $H_đ = 2,05\text{ m}$; $4,09\text{ m}$; $8,19\text{ m}$.
Hạn chế và hướng nghiên cứu tiếp
Văn bản phụ lục bộc lộ một số giới hạn kỹ thuật:
- Thuật toán ma trận và sơ đồ cân bằng lực được xây dựng hoàn toàn trong không gian hai chiều (bài toán biến dạng phẳng 2D), chưa xét đến hiệu ứng không gian ba chiều (3D) của mố cầu và tường cánh góc chéo.
- Các thí dụ kiểm chứng cơ sở giả định môi trường đất dính đồng nhất; trường hợp đất nền phân lớp phức tạp hoặc có sự thay đổi dị hướng của lực dính $c$ và góc ma sát $\varphi$ chưa được thể hiện cụ thể trong tập ma trận mẫu.
- Việc tìm kiếm nghiệm tối ưu bằng Solver trong Excel phụ thuộc vào việc lựa chọn vị trí mặt trượt giả định ban đầu, có nguy cơ rơi vào điểm cực trị cục bộ nếu hình học ban đầu sai lệch lớn so với mặt trượt thực tế.
Hướng nghiên cứu tiếp được xác định gồm việc phát triển mô hình ma trận GLEM cho bài toán trượt không gian ba chiều (3D), tích hợp các thuật toán tối ưu hóa toàn cục (như thuật toán di truyền GA hoặc bầy đàn PSO) để tự động hóa tìm kiếm mặt trượt trong địa tầng nhiều lớp phi đồng nhất.
Giá trị tham khảo
Tài liệu là nguồn tham khảo kỹ thuật cho:
- Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Địa kỹ thuật, Công trình Giao thông, Xây dựng Cầu - Hầm quan tâm đến phương pháp giải tích số và mô hình hóa cân bằng giới hạn tổng quát.
- Các kỹ sư thiết kế kết cấu công trình giao thông đang thực hiện tính toán kiểm toán ổn định trượt sâu mố cầu trên móng nông theo tiêu chuẩn LRFD, cần tài liệu đối soát cấu trúc ma trận lực và tổ hợp tải trọng thi công/khai thác.
- Phần đáng tham khảo nhất là quy trình thiết lập ma trận $[A]$ kích thước $24 \times 24$ và cách thức thiết lập bảng tính Excel kết hợp Solver để giải bài toán cơ học đất phi tuyến.
Câu hỏi thường gặp
1. Hệ phương trình ma trận trong luận án có kích thước bao nhiêu và chứa các ẩn số lực nào?
Hệ phương trình ma trận có dạng $[A]{K} = {B}$ với kích thước vuông $24 \times 24$. Véc-tơ ẩn số ${K}$ gồm 24 đại lượng lực nội tại và tiếp xúc đáy của 6 phân thỏi và 7 mặt biên: các lực pháp tuyến đáy $N_1$ đến $N_6$, các lực tiếp tuyến đáy $T_1$ đến $T_6$, các lực pháp tuyến mặt biên giữa các khối $H_2$ đến $H_7$, và các lực tiếp tuyến mặt biên $V_2$ đến $V_7$.
2. Hệ số an toàn tối ưu ($F_{s\min_tu}$) của mái dốc đất dính đồng nhất được xác định là bao nhiêu?
Trong trường hợp trượt qua chân dốc, hệ số an toàn nhỏ nhất sau tối ưu hóa đạt $1,380$ (khi $FSI = \infty$) và $1,510$ (khi $FSI = 1,7$). Trong trường hợp trượt sâu, hệ số an toàn nhỏ nhất đạt $1,469$ (khi $FSI = \infty$) và $1,654$ (khi $FSI = 1,7$).
3. Tải trọng kết cấu mố cầu và đất đắp trong lòng mố được tổ hợp như thế nào trong bài toán giai đoạn thi công?
Tải trọng được tổ hợp theo Trạng thái giới hạn Cường độ I với hệ số tải trọng $\gamma_P = 1,25$ cho các bộ phận kết cấu mố (bệ mố $Q_{M1} = 2829,75\text{ kN}$, tường thân $Q_{M2} = 2263,80\text{ kN}$, tường cánh $Q_{M3} = 902,40\text{ kN}$, tường đỉnh $Q_{M4} = 294,43\text{ kN}$, vai kê $Q_{M5} = 36,38\text{ kN}$) và hệ số tải trọng $\gamma_P = 1,35$ cho tĩnh tải đất đắp trong lòng mố ($Q_{soil} = 6096,87\text{ kN}$ với $H_đ = 2,05\text{ m}$ và $8,19\text{ m}$; $Q_{soil} = 3223,94\text{ kN}$ với $H_đ = 4,09\text{ m}$). Tổng lực đứng tính toán tương ứng là $16139,22\text{ kN}$ và $12260,76\text{ kN}$.
4. Luận án giải bài toán tìm mặt trượt nguy hiểm nhất bằng công cụ tính toán nào?
Tác giả sử dụng trực tiếp công cụ Solver tích hợp trong Microsoft Excel để thực hiện bài toán tối ưu hóa phi tuyến, với hàm mục tiêu là tìm giá trị nhỏ nhất của hệ số an toàn $F_s$ thông qua việc điều chỉnh các tọa độ điểm nút hình học $(X_i, Y_i)$ của mặt trượt.
5. Các tham số $k$ và $m$ trong bảng dữ liệu biểu diễn điều gì?
Các tham số $k$ và $m$ là các đại lượng quy ước không thứ nguyên được sử dụng để xác định chiều và hướng của lực tiếp tuyến (lực cắt) tác dụng trên mặt phẳng đáy khối trượt ($T_i$) và trên mặt phẳng tiếp xúc giữa các khối trượt ($V_i$).
Kết luận
Bản văn bản là tập phụ lục tính toán số chuyên sâu của luận án tiến sĩ kỹ thuật thuộc Trường Đại học Giao thông Vận tải (năm 2023), tập trung vào bài toán ổn định mái dốc và ổn định trượt sâu mố cầu trên móng nông. Công trình đã chuẩn hóa và ma trận hóa thành công phương pháp cân bằng giới hạn tổng quát GLEM thành hệ ma trận $[A]{K} = {B}$ kích thước $24 \times 24$ giải bằng công cụ Solver trong Microsoft Excel. Các kết quả tính toán số cung cấp dữ liệu định lượng cụ thể về hệ số an toàn và tọa độ mặt trượt nguy hiểm nhất dưới tác động biến thiên của chiều cao đất đắp, thiết bị thi công, hoạt tải khai thác, mực nước ngầm và động đất theo tiêu chuẩn trạng thái giới hạn.
