Tổng quan nghiên cứu

Kết cấu dàn thép đóng vai trò trọng yếu trong kỹ thuật xây dựng công trình hiện đại, đặc biệt tại các dự án khẩu độ lớn với nhịp vượt từ 180 mét đến 207 mét như sân vận động Astrodome hay nhà thi đấu Superdome. Nhờ khả năng tận dụng tối đa dung lượng chịu tải kéo và nén thuần túy dọc trục thanh, hệ dàn giúp giảm từ 20% đến 30% khối lượng vật liệu so với dầm đặc tương đương. Tuy nhiên, việc tính toán phân tích chính xác nội lực và biến dạng chuyển vị cho các kết cấu dàn phức tạp, đặc biệt là hệ siêu tĩnh bậc cao, luôn đòi hỏi các công cụ toán cơ học mạnh mẽ và tối ưu.

Các phương pháp giải tích cổ điển như phương pháp tách nút, phương pháp mặt cắt đơn giản, mặt cắt phối hợp hay họa đồ Maxwell-Cremona bộc lộ nhiều hạn chế khi chỉ giải quyết thuận lợi các hệ tĩnh định. Đối với hệ siêu tĩnh, phương pháp lực và phương pháp chuyển vị truyền thống thường đòi hỏi khối lượng biến đổi ma trận phức tạp và khó khăn trong việc mở rộng cho các bài toán có liên kết một chiều. Luận văn thạc sĩ kỹ thuật của tác giả Nguyễn Thanh Tuấn, thực hiện năm 2017 tại Trường Đại học Dân lập Hải Phòng dưới sự hướng dẫn khoa học của Tiến sĩ Phạm Văn Đạt, đã đề xuất một hướng tiếp cận mới: ứng dụng Nguyên lý cực trị Gauss để phân tích tuyến tính bài toán nội lực và chuyển vị của kết cấu dàn chịu tải trọng tĩnh.

Nghiên cứu tập trung giải quyết bài toán dàn phẳng và dàn không gian dựa trên 4 giả thiết cơ học kinh điển, xây dựng hàm phiếm hàm năng lượng ràng buộc và sử dụng thuật toán số để tìm nghiệm tối ưu. Đóng góp nổi bật của công trình là thiết lập quy trình tính toán tự động hóa hoàn toàn trên nền tảng máy tính, giúp kiểm soát sai số cân bằng nút ở mức tiệm cận 0 tuyệt đối và mở ra khả năng giải quyết các dạng liên kết phi truyền thống trong kỹ thuật công trình.

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

Khung lý thuyết áp dụng

Cơ sở lý thuyết trọng tâm của đề tài được xây dựng từ Nguyên lý cực tiểu Gauss do nhà toán học Carl Friedrich Gauss đề xuất năm 1829, kết hợp với các phát triển mở rộng của Giáo sư Tiến sĩ Khoa học Hà Huy Cương trong cơ học công trình. Nguyên lý Gauss phát biểu rằng chuyển động hoặc biến dạng thực của một cơ hệ chịu liên kết tùy ý sẽ diễn ra sao cho lượng ràng buộc đạt giá trị cực tiểu so với trạng thái hoàn toàn tự do.

Về mặt toán học, hệ thống liên kết được mô tả thông qua bất đẳng thức Fourier - Gauss - Ostrogradsky nhằm tổng quát hóa điều kiện cân bằng cho cả liên kết hai chiều và liên kết một chiều (liên kết không giữ). Phiếm hàm lượng ràng buộc cơ học đối với hệ thanh chịu biến dạng đàn hồi tuân theo định luật Hooke được thiết lập dưới dạng bình phương tối thiểu:

$$Z = \sum_{i=1}^{n} \frac{(N_i - N_i^0)^2 l_i}{E_i F_i} - 2 \sum_{j=1}^{r} P_j v_j \to \min$$

Trong đó, $N_i$ là nội lực dọc trục thanh cần tìm, $N_i^0$ là nội lực trong hệ so sánh tự do, $E_i F_i$ là độ cứng kéo nén của tiết diện, $l_i$ là chiều dài thanh, $P_j$ là tải trọng nút và $v_j$ là chuyển vị tương ứng tại nút mang tải. Mô hình lý thuyết khảo sát 2 cách tiếp cận biến phân độc lập: chọn ẩn số chính là các thành phần chuyển vị nút $(u, v, w)$ hoặc chọn ẩn số chính là các thành phần nội lực $N_i$ kết hợp phương pháp thừa số Lagrange để thỏa mãn điều kiện liên tục hình học.

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

Phương pháp nghiên cứu thực hiện qua việc mô hình hóa số học và kiểm chứng chéo trên các mẫu kết cấu tiêu chuẩn. Cỡ mẫu thực nghiệm bao gồm 2 hệ kết cấu đại diện: mô hình dàn phẳng 31 thanh liên kết tại 14 nút chịu tải trọng nút 40 kN và mô hình hệ dàn vòm tĩnh định, siêu tĩnh chịu tải tập trung 30 kN. Phương pháp chọn mẫu tập trung vào các kết cấu có số bậc tự do từ 25 ẩn số trở lên nhằm đánh giá toàn diện tính hội tụ và độ ổn định của giải thuật.

Quy trình phân tích số được thiết lập tự động hóa qua 5 giai đoạn tuần tự trên phần mềm MATLAB:

  1. Số hóa tọa độ nút, sơ đồ nối thanh và thông số vật liệu $EF = 800.000\text{ kN}$.
  2. Thiết lập biểu thức phiếm hàm lượng ràng buộc $Z$ theo hàm mục tiêu chuyển vị hoặc nội lực.
  3. Lập hệ phương trình đại số tuyến tính từ điều kiện cực trị bậc nhất (đạo hàm riêng theo các ẩn số bằng 0).
  4. Giải hệ phương trình tuyến tính bằng hàm fsolve thuộc hộp công cụ Optimization Toolbox.
  5. Trích xuất vector nghiệm, tính toán nội lực các thanh và kiểm tra cân bằng tĩnh học tại từng nút.

Lý do lựa chọn phương pháp phân tích số trên MATLAB kết hợp thuật toán tối ưu là nhằm loại bỏ các bước lập ma trận độ cứng phần tử phức tạp của phương pháp phần tử hữu hạn truyền thống, đồng thời đảm bảo thời gian xử lý nhanh chóng và độ chính xác số học cao.

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

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

Nghiên cứu mang lại những kết quả định lượng chuẩn xác, chứng minh tính khả thi và ưu thế vượt trội của phương pháp nguyên lý cực trị Gauss trong cơ học kết cấu dàn:

Thứ nhất, trên mô hình dàn phẳng thử nghiệm gồm 31 phần tử thanh và 14 nút với độ cứng $EF = 800.000\text{ kN}$, thuật toán giải quyết chính xác hệ phương trình gồm 25 ẩn số chuyển vị chưa biết. Chuyển vị ngang cực đại ghi nhận tại nút 8 đạt $u = 0,4812\text{ cm}$, trong khi độ võng đứng lớn nhất xuất hiện đối xứng tại nút 11 với $v = -1,3214\text{ cm}$ và nút 5 với $v = -1,3156\text{ cm}$.

Thứ hai, phân tích nội lực ghi nhận lực kéo lớn nhất tại các thanh cánh dưới số 3 và số 4 với giá trị $+229,652\text{ kN}$, trong khi lực nén cực đại tập trung tại các thanh cánh trên số 10 và số 11 với giá trị $-223,681\text{ kN}$. Các thanh bụng xiên chịu lực biến thiên từ $+72,077\text{ kN}$ (thanh 20, 31) đến $-94,590\text{ kN}$ (thanh 25, 26).

Thứ ba, sai số cân bằng tĩnh học tại tất cả 14 nút dàn được triệt tiêu hoàn toàn. Độ lệch tổng hình chiếu lực theo phương ngang $\Delta F_x$ và phương đứng $\Delta F_y$ dao động trong khoảng từ $-4,45 \cdot 10^{-12}\text{ kN}$ đến $+6,09 \cdot 10^{-13}\text{ kN}$. Đây là mức sai số xấp xỉ bằng 0 tuyệt đối, khẳng định sự cân bằng hoàn hảo của hệ thống phương trình nghiệm.

Thứ tư, khi so sánh đối chiếu 100% dữ liệu nội lực của 31 thanh dàn với phần mềm mô phỏng thương mại chuẩn quốc tế SAP2000, độ lệch giữa 2 phương pháp đạt mức $0,00%$. Tất cả các giá trị nội lực từ thanh 1 đến thanh 31 đều trùng khớp chính xác đến 3 chữ số thập phân sau dấu phẩy.

Thảo luận kết quả

Các kết quả tính toán có thể được trực quan hóa rõ nét thông qua biểu đồ biến dạng hình học của dàn trước và sau khi chịu tải. Sơ đồ biến dạng phản ánh sự dịch chuyển trơn tru của các mắt dàn về phía dưới theo hình parabol đối xứng, hoàn toàn phù hợp với quy luật làm việc đàn hồi tuyến tính của vật liệu thép. Bảng dữ liệu đối chiếu giữa phương pháp Gauss và SAP2000 chứng minh thuật toán tối ưu hóa phiếm hàm ràng buộc đạt độ tin cậy ngang hàng với phương pháp phần tử hữu hạn.

So với phương pháp chuyển vị kinh điển, việc tiếp cận từ phiếm hàm Gauss giúp việc thiết lập phương trình trở nên trực quan hơn. Khi giải bài toán theo ẩn số chuyển vị, điều kiện liên tục tại các nút được tự động thỏa mãn mà không cần bổ sung các phương trình tương thích biến dạng phức tạp. Khi giải theo ẩn số nội lực, việc tích hợp phương pháp thừa số Lagrange cho phép trích xuất đồng thời cả trường ứng lực và trường chuyển vị trong cùng một chu trình tính toán, mang lại sự linh hoạt tối đa cho người làm công tác phân tích kết cấu.

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

Dựa trên những phát hiện khoa học từ đề tài, các khuyến nghị kỹ thuật và ứng dụng thực tiễn được đề xuất như sau:

Thứ nhất, các viện nghiên cứu và doanh nghiệp phần mềm xây dựng cần tích hợp thuật toán nguyên lý cực trị Gauss vào các bộ giải (solver) của phần mềm phân tích kết cấu nội địa. Mục tiêu nâng cao tốc độ xử lý bài toán dàn không gian khẩu độ lớn từ 15% đến 20% trong giai đoạn 2026 - 2027.

Thứ hai, nhóm nghiên cứu và các kỹ sư kết cấu nên tiếp tục mở rộng phạm vi áp dụng của phương pháp sang bài toán phi tuyến hình học (biến dạng lớn) và phi tuyến vật liệu. Đích ngắm là hoàn thiện mô hình tính toán cho các kết cấu dàn vòm thép không gian có khẩu độ vượt nhịp trên 200 mét trong vòng 12 đến 24 tháng tới.

Thứ ba, phát triển các module chuyên dụng xử lý bài toán tiếp xúc và gối tựa một chiều (liên kết không giữ). Tận dụng tối đa ưu thế của bất đẳng thức Gauss nhằm giải quyết bài toán dàn làm việc trên nền đàn hồi hoặc kết cấu dàn chịu tải trọng động đất có vùng tiêu tán năng lượng dẻo (hệ khung STMFs).

Thứ tư, các trường đại học khối kỹ thuật công trình cần đưa phương pháp cực trị Gauss vào chương trình giảng dạy chuyên đề Cơ học kết cấu nâng cao và Tin học ứng dụng công trình, đặt chỉ tiêu 100% học viên cao học nắm vững quy trình lập trình tự động hóa giải thuật trên MATLAB.

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

Công trình nghiên cứu mang lại giá trị học thuật và ứng dụng thiết thực cho 4 nhóm đối tượng trọng tâm:

  1. Kỹ sư thiết kế kết cấu công trình: Ứng dụng phương pháp để xây dựng công cụ kiểm tra độc lập (independent check), đối soát nhanh kết quả nội lực và chuyển vị của các hệ dàn khẩu độ lớn trong các dự án nhà thi đấu, sân vận động, nhà xưởng công nghiệp.
  2. Học viên cao học và nghiên cứu sinh ngành Xây dựng - Cơ học: Sử dụng tài liệu làm cơ sở lý thuyết tham khảo về nguyên lý biến phân cơ học, phương pháp lượng ràng buộc tối thiểu và kỹ thuật áp dụng toán tối ưu trong phân tích kết cấu.
  3. Giảng viên các trường đại học kỹ thuật: Bổ sung nguồn học liệu tham khảo phong phú cho các học phần Cơ học kết cấu, Phương pháp tính trong xây dựng, phân tích hệ thanh và kết cấu công trình dân dụng - công nghiệp.
  4. Các lập trình viên và chuyên gia phát triển phần mềm CAE: Tham khảo sơ đồ khối giải thuật và mã nguồn MATLAB mẫu để phát triển các thuật toán phần tử hữu hạn thế hệ mới, tối ưu hóa ma trận độ cứng cho các bài toán kết cấu siêu tĩnh.

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

Nguyên lý cực trị Gauss có điểm gì khác biệt so với phương pháp phần tử hữu hạn thông thường? Phương pháp phần tử hữu hạn truyền thống thiết lập ma trận độ cứng qua tích phân trên từng phần tử rồi ghép nối ma trận tổng thể. Trong khi đó, nguyên lý Gauss tiếp cận trực tiếp từ phiếm hàm lượng ràng buộc năng lượng toàn hệ thống, giúp biểu diễn bài toán cơ học dưới dạng bài toán tối ưu hóa thuần túy và dễ dàng tích hợp các điều kiện ràng buộc bất đẳng thức.

Vì sao luận văn lại xây dựng hai cách chọn ẩn số chuyển vị và nội lực riêng biệt? Cách chọn ẩn chuyển vị giúp tự động thỏa mãn tính liên tục hình học tại các nút, phù hợp cho lập trình giải nhanh bài toán đàn hồi. Cách chọn ẩn nội lực kết hợp nhân tử Lagrange lại cho phép xác định trực tiếp ứng lực thanh mà không qua bước vi phân trung gian, rất hữu ích khi cần kiểm tra giới hạn chịu lực của thanh.

Độ chính xác của phương pháp khi áp dụng cho dàn siêu tĩnh bậc cao như thế nào? Phương pháp đạt độ chính xác tương đương 100% so với phần mềm SAP2000. Thực nghiệm trên mô hình dàn 31 thanh cho thấy sai số cân bằng tĩnh học tại các nút chỉ ở mức $10^{-12}$ đến $10^{-14}\text{ kN}$, triệt tiêu hoàn toàn sự sai lệch số học trong quá trình hội tụ.

Thuật toán trong luận văn được lập trình trên công cụ nào và có khó triển khai không? Thuật toán được xây dựng hoàn toàn trên phần mềm MATLAB 7.0 bằng hàm giải fsolve thuộc Optimization Toolbox. Cấu trúc chương trình đơn giản, gồm 5 bước tuần tự từ nhập tọa độ đến xuất kết quả, cho phép kỹ sư dễ dàng tùy biến mã nguồn theo từng bài toán cụ thể.

Phương pháp này có thể mở rộng cho các bài toán kết cấu phức tạp hơn không? Hoàn toàn có thể. Do được xây dựng trên nền tảng bất đẳng thức Gauss, phương pháp có tiềm năng to lớn để mở rộng giải quyết các bài toán phi tuyến hình học, phân tích biến dạng dẻo giảm chấn động đất STMFs, và phân tích các hệ kết cấu có liên kết một chiều không giữ.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp tiếp cận mới phân tích nội lực và chuyển vị kết cấu dàn tuyến tính dựa trên nguyên lý cực trị Gauss và bất đẳng thức Fourier.
  • Thiết lập hoàn chỉnh 2 giải thuật toán học độc lập: giải theo ẩn số chuyển vị nút và giải theo ẩn số nội lực thanh có bổ sung điều kiện liên tục bằng thừa số Lagrange.
  • Ứng dụng thành công phần mềm MATLAB 7.0 để tự động hóa quá trình giải số, kiểm soát sai số cân bằng lực tại các nút ở mức cực nhỏ $10^{-12}\text{ kN}$.
  • Kiểm chứng kết quả trên mô hình thực nghiệm 31 thanh, đạt độ trùng khớp $100%$ so với phần mềm SAP2000, khẳng định độ tin cậy khoa học tuyệt đối.
  • Định hướng phát triển giai đoạn 2026 - 2027 tập trung mở rộng giải thuật cho dàn không gian phi tuyến và tích hợp vào các phần mềm tính toán kết cấu thực tế. Hãy tải ngay toàn văn tài liệu để làm chủ phương pháp phân tích kết cấu tiên tiến này!