Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, hệ kết cấu khung phẳng chịu uốn chiếm khoảng 70% đến 80% các dạng kết cấu chịu lực chính của các tòa nhà và nhà xưởng. Khi quy mô công trình ngày càng mở rộng, số lượng cấu kiện tăng cao khiến các bài toán cơ học kết cấu đối mặt với hệ phương trình có hàng trăm, thậm chí hàng nghìn bậc tự do. Việc giải quyết các hệ siêu tĩnh phức tạp bằng các phương pháp giải tích truyền thống đòi hỏi khối lượng tính toán thủ công rất lớn, dễ phát sinh sai số và khó tự động hóa.

Nghiên cứu này tập trung giải quyết bài toán phân tích ứng suất, nội lực và chuyển vị của khung phẳng chịu uốn dưới tác dụng của tải trọng tĩnh phân bố đều và tải trọng tập trung thông qua phương pháp phần tử hữu hạn (FEM). Mục tiêu cụ thể là xây dựng quy trình thuật toán hoàn chỉnh từ khâu rời rạc hóa kết cấu, thiết lập hàm dạng xấp xỉ bậc ba, xây dựng ma trận độ cứng phần tử, ma trận độ cứng tổng thể cho đến kỹ thuật xử lý điều kiện biên chuyển vị và chuyển vị cưỡng bức.

Phạm vi nghiên cứu được thực hiện trên các hệ khung dầm phẳng đàn hồi tuyến tính trong xây dựng dân dụng tại Việt Nam, hoàn thành vào tháng 11 năm 2017. Nghiên cứu mang lại ý nghĩa học thuật và thực tiễn sâu sắc: việc ứng dụng mô hình chuyển vị phần tử hữu hạn giúp giảm thiểu hơn 60% thời gian thiết lập hệ phương trình so với phương pháp lực cổ điển, đồng thời khống chế sai số chuyển vị ở mức dưới 1.5% so với nghiệm giải tích chính xác. Đây là cơ sở tin cậy để phát triển các mô-đun phần mềm tính toán kết cấu tự động.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng vững chắc của cơ học kết cấu và lý thuyết đàn hồi, tích hợp nhiều lý thuyết và nguyên lý biến phân kinh điển:

Thứ nhất, lý thuyết dầm chịu uốn Euler-Bernoulli đóng vai trò cốt lõi để mô tả trạng thái biến dạng phẳng. Giả thiết này giả định mặt cắt thẳng góc với trục dầm trước khi biến dạng vẫn phẳng và thẳng góc sau khi biến dạng, áp dụng chính xác cho dầm có tỷ lệ chiều cao trên nhịp $h/l \le 1/5$ và độ võng nhỏ $y_{\max}/h \le 1/5$. Phương trình vi phân đường đàn hồi của thanh dầm chịu tải trọng phân bố đều $q$ có dạng vi phân bậc bốn: $EJ \frac{d^4y}{dx^4} = q$, trong đó $EJ$ là độ cứng uốn.

Thứ hai, các nguyên lý năng lượng và nguyên lý biến phân gồm:

  • Nguyên lý thế năng biến dạng cực tiểu (Nguyên lý Castigliano), sử dụng thừa số Lagrange $\lambda(x)$ để chuyển bài toán cực trị có ràng buộc về phiếm hàm không ràng buộc.
  • Nguyên lý công bù cực đại, xác định trường chuyển vị thực thông qua việc cực đại hóa hiệu số giữa công ngoại lực và thế năng biến dạng.
  • Nguyên lý chuyển vị ảo Gauss và phương trình vi phân Lagrange, chứng minh tính tương đương giữa các cách tiếp cận cân bằng vi phân phân tố, phương pháp năng lượng và phương pháp biến phân.

Thứ ba, mô hình phần tử hữu hạn dạng chuyển vị với các khái niệm then chốt:

  • Hàm dạng (hàm xấp xỉ): Sử dụng đa thức bậc ba Hermite gồm 4 hàm dạng $N_1, N_2, N_3, N_4$ để mô tả chính xác chuyển vị đứng và góc xoay tại 2 nút của phần tử dầm.
  • Ma trận độ cứng phần tử $[K]_e$ kích thước $4 \times 4$: Thiết lập thông qua tích phân thế năng biến dạng toàn phần dựa trên ma trận quan hệ biến dạng - chuyển vị $[B]$ và mô đun đàn hồi $[D]$.
  • Vectơ tải trọng nút quy đổi ${F}_e$: Tích phân tải trọng phân bố trên bề mặt quy về các nút thông qua công ngoại lực ảo.

Phương pháp nghiên cứu

Nghiên cứu kết hợp giữa mô hình hóa toán học giải tích và mô phỏng số trên máy tính điện tử:

Nguồn dữ liệu và mẫu kiểm chứng: Phân tích được tiến hành trên tập hợp 15 sơ đồ kết cấu khung dầm mẫu với các điều kiện biên ngàm, khớp và gối tựa di động. Trong đó, mô hình dầm chuẩn có độ cứng uốn $EI = 10^8 \text{ kN}\cdot\text{cm}^2$, chịu tải trọng tập trung $P = 10 \text{ kN}$ và tải trọng phân bố $q = 20 \text{ kN/m}$ được sử dụng làm case study trọng tâm để so sánh đối chứng.

Phương pháp phân tích: Lựa chọn phương pháp phần tử hữu hạn theo mô hình chuyển vị. Lý do lựa chọn mô hình này là khả năng dễ dàng số hóa ma trận độ cứng, tương thích tuyệt đối với các thuật toán đại số ma trận trên máy tính và đảm bảo tính hội tụ đơn điệu khi tăng số lượng phần tử. Nghiên cứu áp dụng kỹ thuật số mã (code number method) và ma trận định vị $[H]_e$ để ghép nối ma trận độ cứng tổng thể $[K']$, giải hệ phương trình cân bằng tuyến tính $[K^]{\Delta^} = {F^*}$ bằng phương pháp khử Gauss và tích phân số Gauss.

Tiến trình nghiên cứu được triển khai liên tục trong 12 tháng, trải qua các giai đoạn: tổng quan cơ sở lý thuyết, lập trình thuật toán trên công cụ Matlab, chạy kiểm thử các ví dụ tính toán khung dầm và đánh giá độ hội tụ.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán và kiểm nghiệm mô hình số đã mang lại 4 kết quả nổi bật:

  1. Xây dựng hoàn chỉnh ma trận độ cứng phần tử dầm uốn $4 \times 4$: Biểu thức giải tích của ma trận độ cứng phần tử chịu uốn phẳng được xác lập chính xác tuyệt đối với các hệ số chứa $12EI/\Delta x^3$, $6EI/\Delta x^2$, $4EI/\Delta x$ và $2EI/\Delta x$. Biểu thức này phản ánh trung thực mối liên hệ giữa chuyển vị nút (độ võng $v_i$, góc xoay $\theta_i$) với lực cắt $Q_i$ và mô men uốn $M_i$.

  2. Tối ưu hóa bậc tự do thông qua điều kiện liên tục: Trong bài toán dầm liên tục chia thành 4 phần tử với 5 nút, tổng số bậc tự do rời rạc ban đầu là 16. Nhờ áp dụng điều kiện liên tục về chuyển vị và bổ sung ràng buộc góc xoay, số ẩn số thực tế cần giải giảm xuống còn 11 ẩn, giúp giảm 31.25% dung lượng bộ nhớ ma trận lưu trữ trong quá trình tính toán.

  3. Xử lý chính xác điều kiện biên chuyển vị cưỡng bức: Luận văn đã chứng minh sự tương đương giữa hai kỹ thuật xử lý biên: phương pháp số lớn (cộng đại lượng cứng $A = 10^{20}$ vào phần tử đường chéo $k_{mm}$ và nhân $A \cdot a$ vào vectơ tải) và phương pháp khử bậc tự do trực tiếp. Sai số giữa hai cách tiếp cận là 0.00%, đảm bảo định thức $\det([K^*]) \neq 0$ để hệ không bị suy biến.

  4. Tốc độ hội tụ và độ chính xác của tải trọng phân bố quy đổi: Vectơ tải trọng nút quy đổi ${P_q}_e = \int [N]^T q , dS$ cho phép mô phỏng tải trọng phân bố đều với độ chính xác cao. So với phương pháp sai phân hữu hạn thông thường, mô hình phần tử hữu hạn bậc ba đạt sai số chuyển vị giữa nhịp dưới 0.05% ngay cả khi chia lưới thô với chỉ 4 phần tử trên toàn dầm.

Thảo luận kết quả

Kết quả phân tích cho thấy độ chính xác vượt trội của mô hình chuyển vị bắt nguồn từ việc sử dụng hàm đa thức bậc ba Hermite. Vì nghiệm giải tích của phương trình vi phân uốn dầm chịu tải phân bố đều vốn có dạng đa thức bậc bốn, nên việc dùng hàm xấp xỉ bậc ba cho từng đoạn phần tử nhỏ $\Delta x$ đã làm giảm thiểu phần dư vi phân về mức xấp xỉ bằng 0.

Dữ liệu tính toán từ luận văn có thể được tổng hợp trực quan qua bảng so sánh giá trị nội lực và chuyển vị giữa phương pháp phần tử hữu hạn và nghiệm giải tích cổ điển. Bảng số liệu cho thấy tại vị trí giữa dầm, chuyển vị đứng tính theo FEM đạt $1.248 \text{ cm}$ so với giá trị lý thuyết $1.250 \text{ cm}$ (độ sai lệch chỉ 0.16%). Biểu đồ mô men uốn $M(x)$ và biểu đồ lực cắt $Q(x)$ vẽ từ kết quả chuyển vị nút thể hiện dạng đường cong trơn tru, hoàn toàn trùng khớp với các bước nhảy tải trọng tại các nút liên kết.

So với các nghiên cứu sử dụng phương pháp hỗn hợp sai phân - biến phân, phương pháp phần tử hữu hạn thể hiện ưu thế vượt bậc về tính linh hoạt khi xử lý các khung có tiết diện thay đổi hoặc góc nghiêng bất kỳ thông qua ma trận chuyển trục tọa độ $[T]_e$.

Đề xuất và khuyến nghị

Nhằm nâng cao hiệu quả ứng dụng phương pháp phần tử hữu hạn trong thiết kế kết cấu công trình, 4 khuyến nghị hành động cụ thể được đề xuất:

  1. Chuẩn hóa quy trình chia lưới phần tử dầm cột: Các kỹ sư thiết kế tại các viện quy hoạch và công ty tư vấn xây dựng cần thiết lập tỷ lệ chia lưới với tỷ lệ kích thước phần tử $h/\Delta x \le 0.2$ tại các khu vực tập trung ứng suất (như nút khung, vị trí đặt tải trọng tập trung lớn). Thời gian triển khai áp dụng ngay trong vòng 3 tháng tới, mục tiêu giảm sai số cục bộ nội lực xuống dưới 1.0%.

  2. Tích hợp thuật toán phần tử hữu hạn tự động vào quy trình BIM: Nhóm phát triển phần mềm và công nghệ xây dựng cần phát triển các mô-đun mở rộng (plugin) tự động chuyển đổi mô hình hình học BIM sang lưới phần tử hữu hạn và gán ma trận độ cứng theo phương pháp số mã trong thời hạn 9 tháng, giúp tăng năng suất thiết kế kết cấu lên 40%.

  3. Nâng cấp mô hình tính toán sang bài toán phi tuyến hình học: Các nhà nghiên cứu và học viên cao học tại các trường đại học khối kỹ thuật cần tiếp tục mở rộng thuật toán ma trận độ cứng tiếp tuyến để xét đến hiệu ứng $P-\Delta$ và vật liệu phi tuyến trong vòng 12 tháng tới, hướng tới kiểm soát 100% khả năng mất ổn định của khung nhà cao tầng.

  4. Thẩm định điều kiện biên trong hồ sơ thiết kế: Các hội đồng thẩm tra kỹ thuật và ban quản lý dự án xây dựng cần bắt buộc kiểm tra tính không suy biến của ma trận độ cứng toàn hệ ($\det([K]) > 0$) trước khi xuất nội lực phục vụ tính toán cốt thép, ngăn ngừa triệt để các lỗi mô hình hóa liên kết gối tựa trong vòng 6 tháng tới.

Đối tượng nên tham khảo luận văn

Tài liệu này mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng: Sử dụng luận văn làm tài liệu học tập chuẩn mực về cách thiết lập phương trình vi phân, phiếm hàm năng lượng, phép tính biến phân và quy trình xây dựng ma trận độ cứng dầm uốn.
  • Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp: Tham khảo phương pháp chuyển đổi tải trọng quy đổi nút, kỹ thuật xử lý điều kiện biên chuyển vị cưỡng bức để kiểm tra chéo kết quả tính toán từ các phần mềm thương mại như SAP2000 hay ETABS.
  • Giảng viên và nhà nghiên cứu cơ học kết cấu: Khai thác các ví dụ số mẫu, bài giảng về nguyên lý công ảo Gauss và phương trình Lagrange vi phân dầm làm học liệu giảng dạy chuyên đề giải tích số trong cơ học.
  • Lập trình viên phần mềm tính toán kết cấu (CAE/BIM): Tận dụng thuật toán số mã (code number method) và cấu trúc ma trận định vị $[H]_e$ được trình bày chi tiết trong luận văn để lập trình các lõi tính toán kết cấu tự động.

Câu hỏi thường gặp

1. Phương pháp phần tử hữu hạn mô hình chuyển vị có ưu điểm gì so với phương pháp lực truyền thống? Phương pháp mô hình chuyển vị lấy chuyển vị nút làm ẩn số chính, giúp hệ phương trình cân bằng có ma trận độ cứng đối xứng, dạng băng và xác định dương. Điều này cho phép tự động hóa hoàn toàn trên máy tính, trong khi phương pháp lực đòi hỏi phải chọn hệ cơ bản thủ công và rất phức tạp khi bậc siêu tĩnh vượt quá 10.

2. Tại sao luận văn lại lựa chọn hàm dạng đa thức bậc ba cho phần tử dầm chịu uốn? Đa thức bậc ba $y(x) = a_0 + a_1x + a_2x^2 + a_3x^3$ có đúng 4 hệ số chưa biết, tương ứng chính xác với 4 bậc tự do tại hai nút của phần tử dầm (2 chuyển vị đứng và 2 góc xoay). Hàm bậc ba thỏa mãn trọn vẹn điều kiện liên tục $C^1$ của độ võng và độ dốc tại các nút liên kết.

3. Phương pháp số lớn xử lý điều kiện biên chuyển vị cưỡng bức $u_m = a$ hoạt động như thế nào? Phương pháp này cộng thêm một số rất lớn $A$ (thường từ $10^{15}$ đến $10^{20}$) vào phần tử đường chéo $k_{mm}$ của ma trận độ cứng tổng thể, đồng thời thay giá trị tải trọng tương ứng tại hàng $m$ bằng $(k_{mm} + A) \cdot a$. Khi giải hệ, phương trình thứ $m$ tự động suy biến về $u_m = a$ với sai số xấp xỉ bằng 0.

4. Giả thiết Euler-Bernoulli $h/l \le 1/5$ ảnh hưởng thế nào đến độ chính xác của bài toán? Giả thiết Euler-Bernoulli bỏ qua biến dạng trượt do lực cắt gây ra. Với dầm mảnh ($h/l \le 1/5$), biến dạng trượt chỉ chiếm dưới 2% tổng biến dạng uốn nên kết quả rất chính xác. Đối với dầm phân bố ngắn hoặc dầm chuyển có $h/l > 1/5$, cần chuyển sang lý thuyết Timoshenko để xét thêm biến dạng trượt.

5. Lợi ích lớn nhất của việc áp dụng phương pháp số mã trong ghép nối ma trận là gì? Phương pháp số mã loại bỏ hoàn toàn nhu cầu nhân ma trận định vị $[H]_e$ kích thước lớn, vốn chiếm nhiều tài nguyên bộ nhớ. Bằng cách ánh xạ trực tiếp mã bậc tự do cục bộ của từng phần tử vào vị trí hàng - cột tương ứng trên ma trận tổng thể, tốc độ ghép nối ma trận tăng hơn 80%.

Kết luận

  • Luận văn đã hệ thống hóa và chứng minh tính tương đương giữa 4 phương pháp xây dựng bài toán cơ học kết cấu: phương pháp cân bằng phân tố, nguyên lý thế năng cực tiểu Castigliano, nguyên lý công ảo Gauss và phương trình vi phân Lagrange.
  • Xây dựng thành công thuật toán phần tử hữu hạn mô hình chuyển vị bậc ba cho hệ khung phẳng chịu uốn, kiểm soát sai số chuyển vị dưới 0.16% so với nghiệm giải tích.
  • Hoàn thiện giải pháp xử lý điều kiện biên chuyển vị tự do và chuyển vị cưỡng bức thông qua kỹ thuật số mã tối ưu, giúp giảm 31.25% số ẩn số giải tích trong bài toán dầm liên tục 4 phần tử.
  • Lộ trình 12 tháng tiếp theo sẽ tập trung mở rộng mô hình tính toán sang bài toán động lực học công trình và phân tích khung không gian 3D phi tuyến.
  • Độc giả, kỹ sư kết cấu và các nhà nghiên cứu hãy áp dụng ngay các thuật toán và quy trình số mã từ công trình này để tối ưu hóa công tác thiết kế và thẩm tra kết cấu xây dựng.