Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình hiện đại, kết cấu khung thép đóng vai trò chủ đạo khi chiếm khoảng 65% đến 70% các dự án nhà xưởng công nghiệp và nhà nhịp lớn trên thế giới. Dẫu vậy, đa số các tiêu chuẩn thiết kế hiện hành vẫn áp dụng phương pháp phân tích tuyến tính đàn hồi bậc một. Cách tiếp cận này bỏ qua biến dạng lớn và sự phân bố lại nội lực dẻo, gây nên sai lệch từ 15% đến 25% trong việc đánh giá khả năng chịu lực thực tế của công trình. Ngược lại, phân tích phi tuyến hoàn chỉnh bao gồm phi tuyến hình học và phi tuyến vật liệu theo phương pháp vùng dẻo truyền thống lại đòi hỏi chia nhỏ cấu kiện thành hàng trăm thớ sợi vi mô, làm bùng nổ khối lượng tính toán trên các hệ kết cấu quy mô lớn.

Nhằm khắc phục triệt để rào cản tính toán này, luận văn thạc sĩ chuyên ngành Xây dựng công trình dân dụng và công nghiệp với mã số ngành 605820 tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh đã tập trung giải quyết bài toán: thiết lập hàm nội suy đa thức bậc 5 mô phỏng ứng xử bậc hai của cấu kiện dầm-cột, kết hợp phương pháp khớp dẻo hiệu chỉnh để phân tích phi tuyến khung thép phẳng chịu tải trọng tĩnh với quy mô chỉ duy nhất 1 phần tử cho mỗi cấu kiện. Công trình được thực hiện từ tháng 09 năm 2012 đến tháng 12 năm 2013 và bảo vệ thành công vào tháng 01 năm 2014. Ý nghĩa thực tiễn nổi bật của nghiên cứu là cắt giảm từ 50% đến 70% số bậc tự do của hệ phương trình kết cấu so với phương pháp chia lưới phần tử hữu hạn thông thường, đồng thời vẫn kiểm soát sai số tải trọng cực hạn dưới mức 3.0% so với các nghiệm chuẩn quốc tế.

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 dựa trên nền tảng lý thuyết dầm-cột Timoshenko chịu uốn và nén đồng thời, kết hợp lời giải phương trình vi phân cân bằng bậc hai cho cấu kiện chịu tải trọng đầu mút. Hàm chuyển vị ngang được giả thiết có dạng đa thức bậc 5 tổng quát: y(x) = a0 + a1x + a2x^2 + a3x^3 + a4x^4 + a5*x^5. Sáu hệ số ẩn từ a0 đến a5 được xác định tường minh thông qua việc thỏa mãn 3 điều kiện tương thích chuyển vị, góc xoay tại các mút và 3 phương trình cân bằng tĩnh học mô-men uốn cùng lực cắt dọc trục. Khi lực dọc trục tiến về 0, hàm dạng bậc 5 này tự động suy biến hoàn toàn về hàm dạng Hermite bậc 3 truyền thống.

Hiện tượng phi tuyến vật liệu được mô hình hóa bằng lý thuyết khớp dẻo hiệu chỉnh (Refined Plastic Hinge) nhằm loại bỏ sự chuyển trạng thái tức thời phi thực tế của mô hình khớp dẻo cứng cổ điển. Sự suy giảm độ cứng chống uốn được kiểm soát liên tục qua hệ số dẻo e = 4alpha(1 - alpha) khi thông số nội lực alpha di chuyển trong khoảng từ 0.5 đến 1.0. Các khái niệm then chốt xuyên suốt mô hình gồm: hiệu ứng hình học bậc hai P-Delta và P-delta, mô-đun tiếp tuyến Et phản ánh ứng suất dư nhiệt ban đầu theo tiêu chuẩn CRC và Liew, cùng 3 mô hình đường mặt dẻo giới hạn bao gồm Orbison, AISC-LRFD và Balling.

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

Nguồn dữ liệu kiểm chứng của đề tài bao gồm 6 hệ kết cấu chuẩn (benchmark problems) kinh điển được công bố trong các tài liệu khoa học quốc tế uy tín từ năm 1983 đến năm 2010: cột hai đầu khớp W8x31, khung Vogel 1 tầng 1 nhịp, khung Balling 2 tầng 1 nhịp không giằng, khung Zhou 4 tầng 2 nhịp, khung Kassimali 4 tầng 1 nhịp và khung McGuire 1 tầng 2 nhịp. Cỡ mẫu nghiên cứu gồm 6 bài toán số đại diện đầy đủ cho các chế độ làm việc từ mất ổn định nén đúng tâm đàn hồi đến mất ổn định dẻo phi đàn hồi toàn hệ thống.

Phương pháp chọn mẫu là phương pháp chọn mẫu có chủ đích (purposive sampling). Tác giả chủ động chọn lọc các mô hình chuẩn quốc tế đã có sẵn nghiệm giải tích chính xác hoặc nghiệm số từ phương pháp vùng dẻo (Plastic Zone) vốn được xem là chuẩn mực đo lường cao nhất trong cơ học kết cấu.

Phương pháp phân tích được lựa chọn là thuật toán gia số lặp Newton-Raphson đầy đủ, được tác giả lập trình độc lập trên môi trường tính toán khoa học MATLAB. Lý do lựa chọn thuật toán Newton-Raphson đầy đủ thay vì phương pháp Euler hay Newton-Raphson hiệu chỉnh là vì thuật toán này cập nhật ma trận độ cứng tiếp tuyến sau từng vòng lặp con. Nhờ đó, thuật toán khử triệt để lực dư không cân bằng ở mỗi bước tăng tải, đạt tốc độ hội tụ bậc hai và ngăn chặn sự tích lũy sai số khi các khớp dẻo hình thành liên tiếp. Toàn bộ quy trình tính toán được thực hiện trong thời gian 16 tháng với tệp dữ liệu đầu vào chuẩn hóa gồm 25 khối thông số hình học, tiết diện và điều kiện biên.

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

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

Phát hiện thứ nhất: Trên mô hình cột hai đầu khớp tiết diện W8x31 với mô-đun đàn hồi E = 29000 ksi và ứng suất chảy dẻo 36 ksi, khi có xét đến ứng suất dư thông qua mô-đun tiếp tuyến Et, đường cong tải trọng tới hạn phân tích bám sát tuyệt đối đường chuẩn CRC với sai số tối đa chỉ 2.1%. Trong trường hợp không xét ứng suất dư, sai số tải trọng tới hạn so với kết quả của tác giả Kim năm 2009 trên dải độ mảnh từ 40 đến 160 được khống chế ở mức xấp xỉ 2.5%.

Phát hiện thứ hai: Đối với khung Vogel 1 tầng 1 nhịp dùng tiết diện HEA340 chịu tải trọng đứng P = 2800 kN và tải trọng ngang H = 35 kN cùng độ lệch lạc ban đầu 1/400, thuật toán đề xuất xác định hệ số tải trọng cực hạn đạt mức 1.02. Giá trị này hoàn toàn trùng khớp với kết quả từ phương pháp vùng dẻo của Vogel năm 1985, chứng minh khả năng dự báo chính xác sức chịu tải giới hạn mà không cần chia nhỏ cấu kiện.

Phát hiện thứ ba: Ở bài toán khung không giằng Balling 2 tầng 1 nhịp sử dụng thép W14x30 với E = 29000 ksi và ứng suất chảy dẻo 50 ksi, hệ số tải cực hạn phân tích đạt 128.00. So sánh với hệ số tải chuẩn 129.79 của tác giả Balling năm 2010, sai số chỉ dừng ở mức 1.39%.

Phát hiện thứ tư: Tại khung Zhou 4 tầng 2 nhịp sử dụng thép W16x40 có mô-đun đàn hồi 200 GPa và giới hạn chảy 250 MPa, hệ số tải trọng cực hạn đạt chính xác 1.83, sai số tải cực hạn so với công bố của Kukreti và Zhou năm 2006 đạt mức 0%. Đặc biệt, chuyển vị ngang tại đỉnh khung đạt 235.87 mm, lớn hơn giá trị 207.23 mm của mô hình tham chiếu khoảng 13.82%, phản ánh đúng độ suy giảm độ cứng tổng thể của khung khi biến dạng uốn lớn phát triển.

Phát hiện thứ năm: Khi kiểm chứng trên khung Kassimali 4 tầng 1 nhịp chịu tải trọng đứng P = 1334 kN và tỷ lệ tải ngang alpha = 0.24, hệ số tải cực hạn tính toán đạt 1.534. Kết quả này tiệm cận sát sao nghiệm 1.502 của tác giả Kassimali năm 1983 với độ lệch nhỏ chỉ khoảng 2.13%.

Thảo luận kết quả

Nguyên nhân căn bản giúp mô hình đề xuất đạt được độ chuẩn xác cao là nhờ hàm nội suy đa thức bậc 5 thỏa mãn đồng thời các điều kiện cân bằng vi phân dọc theo chiều dài phần tử và phản ánh chính xác hiệu ứng co ngắn thanh do uốn. Khi vật liệu bắt đầu chảy dẻo cục bộ, hệ số suy giảm độ cứng e chuyển tiếp mượt mà từ 1.0 về 0, giúp cấu kiện thích ứng dần với sự suy thoái độ cứng mà không bị gián đoạn số học như mô hình khớp dẻo cứng.

So với các nghiên cứu sử dụng phần tử bậc 3 Hermite thông thường vốn bắt buộc phải chia mỗi thanh thành 4 đến 8 phần tử con để bắt được hiệu ứng P-delta, mô hình của luận văn chỉ cần dùng 1 phần tử cho mỗi cấu kiện. Cách tiếp cận này giúp giảm đáng kể kích thước ma trận độ cứng tổng thể mà vẫn đạt độ tin cậy tương đương với phương pháp vùng dẻo (Plastic Zone) vốn tiêu tốn gấp hàng chục lần thời gian chạy máy.

Trong thực tế trình bày kết quả phân tích kết cấu, ứng xử phi tuyến của công trình được mô tả trực quan và sinh động thông qua biểu đồ quan hệ lực - chuyển vị. Đồ thị thể hiện rõ 3 giai đoạn làm việc: đoạn dốc thẳng ban đầu thể hiện sự làm việc đàn hồi tuyến tính, đoạn cong thoải dần thể hiện sự hình thành khớp dẻo tiệm tiến từ ngưỡng alpha = 0.5, và điểm uốn nằm ngang thể hiện tải trọng tới hạn trước khi kết cấu chuyển sang nhánh dỡ tải mất ổn định. Đi kèm biểu đồ là bảng tổng hợp các chỉ số cơ lý bao gồm hệ số tải cực hạn, chuyển vị đỉnh tại từng tầng và độ lệch phần trăm so với tiêu chuẩn, mang lại cái nhìn minh bạch cho người thiết kế.

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

Tích hợp thuật toán phần tử bậc 5 vào các phần mềm kết cấu thương mại: Triển khai phát triển các module mở rộng hoặc plugin cho các nền tảng tính toán như SAP2000, ETABS và OpenSees thông qua giao diện lập trình ứng dụng. Mục tiêu giảm thiểu từ 50% đến 60% thời gian phân tích phi tuyến bậc hai cho các công trình thép cao trên 15 tầng. Lộ trình thực hiện kéo dài 18 tháng, do Viện Tin học Xây dựng phối hợp cùng các kỹ sư phát triển phần mềm kết cấu đảm nhiệm.

Chuẩn hóa quy trình thiết kế trực tiếp theo tiêu chuẩn xây dựng: Biên soạn tài liệu kỹ thuật hướng dẫn áp dụng phương pháp thiết kế trực tiếp (Direct Analysis Method) sử dụng phân tích phi tuyến bậc hai thay thế cho phương pháp tra hệ số chiều dài tính toán K truyền thống. Mục tiêu giảm từ 8% đến 12% trọng lượng thép tiêu hao trên mỗi công trình mà vẫn đảm bảo hệ số an toàn. Kế hoạch triển khai trong 24 tháng, do Hội Kết cấu và Công nghệ Xây dựng Việt Nam phối hợp với các cơ quan quản lý chuyên ngành ban hành.

Mở rộng mô hình giải thuật cho kết cấu khung không gian 3 chiều: Tiếp tục nghiên cứu bổ sung các bậc tự do xoắn vênh và hiện tượng mất ổn định uốn - xoắn ngoài mặt phẳng vào ma trận độ cứng tiếp tuyến của phần tử bậc 5. Đặt mục tiêu thử nghiệm thành công trên tối thiểu 10 hệ khung không gian nhịp lớn trong thời hạn 12 tháng, do các nhóm nghiên cứu tại các trường đại học khối kỹ thuật chủ trì thực hiện.

Đào tạo chuyển giao công nghệ cho đội ngũ kỹ sư kết cấu: Tổ chức các chương trình tập huấn chuyên sâu về lý thuyết ổn định phi tuyến và kỹ năng phân tích ứng xử khớp dẻo trên mô hình số. Mục tiêu phổ cập kiến thức chuyên môn cho hơn 400 kỹ sư tư vấn thiết kế tại các doanh nghiệp xây dựng trong vòng 12 tháng, do Tổng hội Xây dựng Việt Nam chủ trì điều phối.

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

Kỹ sư tư vấn thiết kế kết cấu công trình thép: Tiếp cận phương pháp phân tích phi tuyến trực tiếp với định dạng 1 phần tử cho mỗi cấu kiện để thẩm tra khả năng chịu tải cực hạn và ổn định tổng thể của các dự án nhà xưởng vượt nhịp lớn, giàn không gian và khung thép nhiều tầng một cách chuẩn xác, loại bỏ hoàn toàn các sai số do giả định chiều dài tính toán K gây ra.

Học viên cao học và nghiên cứu sinh ngành Xây dựng: Sử dụng tài liệu như một nguồn tham khảo học thuật chỉn chu về phương pháp giải phương trình vi phân thanh dầm-cột, kỹ thuật thiết lập ma trận độ cứng tiếp tuyến phi tuyến hình học bậc cao và quy trình lặp Newton-Raphson đa biến trong cơ học kết cấu.

Giảng viên và nhà nghiên cứu tại các viện, trường đại học: Vận dụng cơ sở lý thuyết hàm nội suy bậc 5 và các mô hình mặt dẻo Orbison, AISC-LRFD, Balling để biên soạn giáo trình chuyên đề nâng cao về Cơ học kết cấu phi tuyến, đồng thời mở rộng hướng phát triển cho kết cấu thép liên kết nửa cứng chịu tải trọng động đất.

Chuyên gia phát triển phần mềm tính toán kỹ thuật: Tham khảo cấu trúc tệp dữ liệu đầu vào chuẩn hóa gồm 25 khối thông tin cùng mã nguồn MATLAB để tích hợp các module phân tích phi tuyến tĩnh bậc cao vào các giải pháp phần mềm cơ khí và xây dựng nội địa.

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

Tại sao nên sử dụng hàm nội suy đa thức bậc 5 thay vì hàm bậc 3 Hermite thông thường? Hàm bậc 3 Hermite vốn chỉ giải chính xác bài toán uốn thuần túy. Khi cấu kiện chịu lực nén dọc trục lớn, hiệu ứng P-delta làm gia tăng độ võng phụ, đòi hỏi phải chia nhỏ cấu kiện thành 4 đến 8 phần tử bậc 3 để đạt độ chính xác. Ngược lại, hàm bậc 5 là nghiệm của phương trình vi phân dầm-cột chịu tải đầu mút, cho phép mô phỏng chính xác ứng xử uốn nén kết hợp chỉ với duy nhất 1 phần tử cho mỗi cấu kiện.

Mô hình khớp dẻo hiệu chỉnh khắc phục nhược điểm gì của mô hình khớp dẻo cứng? Mô hình khớp dẻo cứng giả định tiết diện chuyển đột ngột từ trạng thái hoàn toàn đàn hồi sang chảy dẻo hoàn toàn, gây sai lệch lớn và khó khăn cho sự hội tụ số. Khớp dẻo hiệu chỉnh sử dụng hệ số suy giảm e = 4alpha(1 - alpha) trên dải thông số alpha từ 0.5 đến 1.0, giúp mô phỏng chân thực sự chảy dẻo dần dần của mặt cắt dưới tác động kết hợp của lực dọc, mô-men uốn và ứng suất dư.

Việc dùng 1 phần tử cho mỗi cấu kiện mang lại lợi ích tài nguyên cụ thể ra sao? Trong các hệ khung thép phức hợp có hàng trăm cấu kiện, việc dùng 1 phần tử giúp giảm quy mô ma trận độ cứng tổng thể từ 50% đến 70% so với phương pháp chia lưới vi mô. Điều này giúp giảm đáng kể dung lượng bộ nhớ RAM và rút ngắn hàng chục lần thời gian giải hệ phương trình tuyến tính ở mỗi bước tăng tải.

Thuật toán Newton-Raphson xử lý bài toán hội tụ như thế nào khi kết cấu đạt tải tới hạn? Thuật toán tính toán ma trận độ cứng tiếp tuyến tại mỗi bước lặp con nhằm xác định chính xác gia số chuyển vị, qua đó khử dần véc-tơ lực dư giữa ngoại lực tác dụng và nội lực phần tử. Khi kết cấu tiến gần tải cực hạn và xuất hiện các khớp dẻo, việc liên tục cập nhật ma trận tiếp tuyến đảm bảo tốc độ hội tụ bậc hai và kiểm soát sai số dưới ngưỡng cho phép.

Kết quả của luận văn đã được đối chứng và kiểm định qua những cơ sở nào? Chương trình tính toán MATLAB đã được kiểm định chéo với 6 mô hình chuẩn quốc tế kinh điển như cột Euler, khung Vogel năm 1985, khung Balling năm 2010 và khung Zhou năm 2006. Mọi kết quả về tải trọng cực hạn đều tương thích chặt chẽ với phương pháp vùng dẻo và các nghiệm chuẩn với sai số tối đa chỉ dao động từ 0% đến 3.0%.

Kết luận

  • Xây dựng thành công hàm nội suy đa thức bậc 5 mô tả chuẩn xác ứng xử bậc hai của cấu kiện dầm-cột thép phẳng chịu tải trọng đầu mút.
  • Ứng dụng hiệu quả lý thuyết khớp dẻo hiệu chỉnh kết hợp mô-đun tiếp tuyến Et, mô phỏng chân thực sự chảy dẻo tiệm tiến và ảnh hưởng của ứng suất dư nhiệt.
  • Phát triển hoàn chỉnh chương trình phân tích phi tuyến tĩnh trên nền tảng ngôn ngữ MATLAB dựa trên thuật toán gia số lặp Newton-Raphson đầy đủ với độ tin cậy hội tụ cao.
  • Kiểm chứng thành công trên 6 hệ kết cấu chuẩn quốc tế với sai số tải trọng cực hạn khống chế chặt chẽ trong khoảng từ 0% đến 3.0%.
  • Giảm thiểu tối đa khối lượng tính toán khi chỉ cần sử dụng 1 phần tử duy nhất cho mỗi thanh kết cấu, tiết kiệm đáng kể thời gian xử lý dữ liệu.

Đóng góp cốt lõi của luận văn là đã đề xuất một công cụ giải tích số tin cậy, dung hòa tối ưu giữa độ chính xác cơ học nghiêm ngặt của phân tích dẻo bậc hai và tính tinh gọn, thực dụng trong thiết kế công trình xây dựng. Trong giai đoạn 6 đến 12 tháng tiếp theo, nhóm nghiên cứu dự kiến mở rộng mô hình cho bài toán khung không gian 3D chịu tác động của tải trọng động gió và động đất. Các kỹ sư thiết kế và chuyên gia học thuật hãy tham khảo và ứng dụng giải pháp giải tích bậc cao này để nâng cao chất lượng an toàn và tối ưu hóa chi phí cho các công trình kết cấu thép hiện đại.