Tổng quan nghiên cứu

Trong xu hướng phát triển mạnh mẽ của ngành kỹ thuật xây dựng công trình dân dụng và công nghiệp, các tòa nhà cao tầng và tổ hợp thương mại quy mô lớn xuất hiện ngày càng nhiều. Những công trình này thường xuyên đòi hỏi các giải pháp kết cấu phức tạp như dầm chuyển, sàn chuyển và hệ cột vách chịu tải trọng tĩnh cực lớn từ hàng chục tầng truyền xuống nền móng. Điểm đặc thù của dầm chuyển hay cột ngắn là tỷ lệ chiều cao tiết diện so với chiều dài nhịp rất lớn, thường đạt ngưỡng tỷ lệ kích thước vượt quá 0.20.

Thực tế tính toán truyền thống dựa trên lý thuyết thanh mảnh Euler – Bernoulli vốn dĩ bỏ qua biến dạng trượt ngang do lực cắt gây ra. Giả thiết này chỉ phản ánh chính xác trạng thái chịu lực khi cấu kiện có tỷ lệ chiều cao trên chiều dài nhịp nhỏ hơn hoặc bằng 0.10. Khi áp dụng vào kết cấu dầm cao hoặc cột chịu uốn lớn, việc bỏ qua biến dạng trượt dẫn đến sai số độ võng và phân bố nội lực từ 15% đến xấp xỉ 70%, gây nguy cơ mất an toàn hoặc lãng phí vật liệu nghiêm trọng trong thiết kế.

Luận văn thạc sĩ kỹ thuật của tác giả Nguyễn Mạnh Hùng, hoàn thành 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 của Giáo sư Hà Huy Cương, tập trung giải quyết bài toán: "Tính toán khung phẳng chịu uốn có xét đến biến dạng trượt ngang". Mục tiêu cốt lõi của đề tài là xây dựng hệ phương trình cân bằng và thuật toán bán giải tích chính xác dựa trên nguyên lý cực trị Gauss, đồng thời phát triển chương trình tính toán số trên máy tính điện tử. Nghiên cứu cung cấp cơ sở định lượng rõ ràng với 5 cấp tỷ số hình học khảo sát, tạo bước đột phá trong việc kiểm soát biến dạng và chuyển vị của hệ khung phẳng chịu tải trọng tĩnh trong xây dựng hiện đại.

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 sự kế thừa và phát triển từ hai lý thuyết cơ bản trong cơ học kết cấu: lý thuyết dầm cổ điển Euler – Bernoulli và lý thuyết biến dạng trượt bậc nhất Timoshenko. Trong lý thuyết dầm thông thường, giả thiết tiết diện phẳng vuông góc với trục dầm sau biến dạng khiến ứng suất tiếp và biến dạng trượt bị triệt tiêu trong công thức tính độ võng. Ngược lại, lý thuyết Timoshenko thừa nhận mặt cắt xoay một góc bổ sung do lực cắt, nhưng cách tiếp cận truyền thống sử dụng hai ẩn chuyển vị độ võng và góc xoay thường gặp hiện tượng khóa trượt khi chiều dày cấu kiện tiến về 0.

Để khắc phục nhược điểm trên, luận văn áp dụng khung lý thuyết cải tiến lấy hàm độ võng và hàm lực cắt làm hai biến phân độc lập. Bốn khái niệm trọng tâm bao gồm:

  • Độ cứng uốn: Tích số giữa mô đun đàn hồi của vật liệu và mô men quán tính của tiết diện.
  • Độ cứng cắt: Tích số giữa mô đun trượt và diện tích mặt cắt ngang.
  • Hệ số phân bố ứng suất trượt: Hệ số tập trung ứng suất cắt tại trục trung hòa, nhận giá trị bằng 1.20 đối với tiết diện hình chữ nhật.
  • Phiếm hàm lượng cưỡng bức Gauss: Biểu thức năng lượng biến dạng toàn phần kết hợp công của ngoại lực phân bố.

Mối quan hệ biến dạng trượt được xác định trực tiếp qua lực cắt chia cho độ cứng cắt, giúp công thức toán học tự động hội tụ chính xác về nghiệm lý thuyết dầm cổ điển khi chiều cao tiết diện giảm dần về 0.

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

Nghiên cứu sử dụng phương pháp biến phân dựa trên nguyên lý cực trị Gauss kết hợp kỹ thuật nhân tử Lagrange để giải bài toán cơ học kết cấu phẳng. Quá trình phân tích bao gồm việc thiết lập phiếm hàm dừng cho hệ thanh, đưa các điều kiện biên hình học và điều kiện liên tục tại nút khung vào hệ phương trình ràng buộc.

Về cỡ mẫu và đối tượng khảo sát, tác giả xây dựng 2 mô hình kết cấu tiêu biểu: hệ khung 2 thanh siêu tĩnh bậc 2 và hệ khung phẳng 1 tầng 1 nhịp gồm 3 thanh liên kết ngàm. Trên mỗi mô hình, 5 kịch bản tỷ lệ hình học giữa chiều cao tiết diện và chiều dài nhịp được phân tích chi tiết, bao gồm các tỷ số: 1/100, 1/10, 1/5 và 1/3. Lý do lựa chọn phương pháp bán giải tích thông qua hệ đa thức xấp xỉ bậc 4 là nhằm chuyển đổi hệ phương trình vi phân phức tạp thành hệ phương trình đại số tuyến tính có bậc tự do từ 24 đến 29 phương trình. Timeline nghiên cứu được thực hiện xuyên suốt trong năm 2017, đảm bảo tính chuẩn xác thông qua các thuật toán giải số được lập trình hóa hoàn toàn.

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

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

Kết quả tính toán trên các mô hình khung phẳng đem lại những phát hiện định lượng mang tính bước ngoặt:

Thứ nhất, biến dạng trượt ngang làm tăng đáng kể độ võng của các cấu kiện có tiết diện lớn. Đối với khung hai thanh chịu tải trọng phân bố đều, khi tỷ số chiều cao trên chiều dài nhịp ở mức 1/100, chênh lệch độ võng giữa mô hình có xét và không xét trượt là 0.00%. Tuy nhiên, khi tỷ số này tăng lên 1/10, độ võng tăng thêm 0.37%; ở mức 1/5, chênh lệch độ võng đạt 34.02% và chạm ngưỡng tăng 11.86% tại tỷ số 1/3.

Thứ hai, đối với khung một tầng một nhịp chịu tải trọng thẳng đứng, ảnh hưởng của lực cắt lên chuyển vị giữa nhịp thanh ngang càng thể hiện rõ rệt. Khi tỷ số kích thước đạt 1/10, chênh lệch độ võng ghi nhận mức 0.69%; khi tỷ số đạt 1/5, mức chênh lệch tăng lên 27.77% và đạt con số kỷ lục 69.56% tại tỷ số 1/3.

Thứ ba, sự thay đổi của mô men uốn tại các nút khung chịu ảnh hưởng tương đối nhỏ từ biến dạng trượt. Cụ thể, mô men nút khung chỉ biến động trong khoảng từ 0.50% đến 1.50% khi tỷ số chiều cao trên nhịp biến thiên từ 1/10 đến 1/3. Biến động lực cắt dọc trục thanh dao động cục bộ từ 0.10% đến 1.10%.

Thảo luận kết quả

Các phát hiện trên phản ánh bản chất cơ học của vật liệu đàn hồi khi làm việc ở trạng thái ứng suất phẳng phức tạp. Nguyên nhân chính khiến độ võng tăng vọt lên tới 69.56% ở cấu kiện dầm cao là do năng lượng biến dạng trượt tích lũy dọc theo chiều cao dầm chiếm tỷ trọng lớn trong tổng năng lượng biến dạng của toàn hệ kết cấu. Đối với thanh mảnh thông thường, độ cứng uốn kiểm soát hoàn toàn chuyển vị; nhưng đối với dầm chuyển có tỷ số chiều cao trên chiều dài từ 1/5 đến 1/3, độ cứng cắt trở thành nhân tố quyết định.

Dữ liệu tính toán từ luận văn có thể được biểu diễn trực quan qua các bảng tổng hợp số liệu so sánh và biểu đồ biến thiên độ võng theo từng tỷ lệ kích thước. Khi trực quan hóa đường độ võng trên đồ thị hai chiều, đường cong võng có xét biến dạng trượt tách xa rõ rệt khỏi đường đàn hồi Euler – Bernoulli truyền thống tại khu vực giữa nhịp. So sánh với các nghiên cứu trước đây vốn phụ thuộc vào phương pháp phần tử hữu hạn thương mại, phương pháp bán giải tích của luận văn giúp loại bỏ hoàn toàn sai số do hiện tượng khóa trượt, đồng thời bảo toàn tính liên tục của hàm mô men uốn và hàm lực cắt tại các điểm nút liên kết.

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

Dựa trên các kết luận khoa học vững chắc, luận văn đưa ra 4 nhóm đề xuất và khuyến nghị kỹ thuật có tính ứng dụng cao:

  1. Cập nhật quy trình tính toán thiết kế: Các tổ chức tư vấn thiết kế cần bắt buộc áp dụng bài toán uốn có xét biến dạng trượt đối với mọi cấu kiện dầm chuyển, dầm vách hoặc cột ngắn có tỷ số kích thước lớn hơn 0.20, nhằm kiểm soát độ võng thực tế chính xác hơn 30% so với phương pháp cũ.
  2. Bổ sung hệ số hiệu chỉnh vào tiêu chuẩn xây dựng: Khuyến nghị Viện Khoa học Công nghệ Xây dựng và các cơ quan quản lý ban hành hướng dẫn bổ sung trong vòng 12 đến 24 tháng tới, quy định cụ thể hệ số phân bố ứng suất trượt bằng 1.20 cho dầm bê tông cốt thép tiết diện chữ nhật.
  3. Tích hợp thuật toán vào phần mềm kết cấu: Các kỹ sư phát triển phần mềm mô phỏng xây dựng cần lập trình tích hợp hệ phương trình giải tích hai ẩn độ võng và lực cắt vào các mô-đun phân tích khung phẳng, đặt mục tiêu hoàn thiện gói thuật toán trong lộ trình 6 tháng.
  4. Mở rộng phạm vi nghiên cứu thực nghiệm: Các trường đại học khối kỹ thuật và viện nghiên cứu cần tiếp tục thực hiện các chương trình thử nghiệm tải trọng động và nghiên cứu ứng xử phi tuyến vật liệu của khung phẳng dầm cao trong giai đoạn 2026 - 2028.

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

Tài liệu này là nguồn tư liệu học thuật và thực tiễn chuyên sâu, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Kỹ sư kết cấu công trình: Nắm vững phương pháp tính toán chuẩn xác cho các hệ kết cấu đặc biệt như dầm chuyển nhà cao tầng, móng bè sườn, sàn chuyển ứng suất trước, giúp tối ưu hóa việc bố trí cốt thép và kiểm tra võng nứt.
  2. Giảng viên và nghiên cứu sinh ngành Xây dựng: Sử dụng làm tài liệu tham khảo giảng dạy chuyên đề Cơ học kết cấu nâng cao, Lý thuyết biến phân và Phương pháp năng lượng trong cơ học công trình.
  3. Chuyên viên nghiên cứu và phát triển phần mềm kết cấu: Ứng dụng thuật toán ma trận và hệ phương trình đại số tuyến tính 24 đến 29 ẩn số để xây dựng các công cụ phân tích phần tử thanh tổng quát không bị khóa trượt.
  4. Kỹ sư thẩm tra và quản lý dự án xây dựng: Có thêm căn cứ khoa học định lượng để đánh giá mức độ an toàn của hồ sơ thiết kế các công trình có kết cấu cột ngắn và dầm tiết diện lớn.

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

Tại sao cần phải xét đến biến dạng trượt ngang khi tính toán dầm cao và cột ngắn? Khi cấu kiện có chiều cao tiết diện lớn với tỷ số trên nhịp vượt quá 1/5, biến dạng do ứng suất tiếp gây ra chiếm tỷ trọng rất lớn trong tổng năng lượng biến dạng. Bỏ qua yếu tố này sẽ làm sai lệch độ võng thực tế lên tới xấp xỉ 70%, dẫn đến đánh giá sai trạng thái giới hạn sử dụng của công trình.

Phương pháp hai hàm chưa biết của tác giả có ưu điểm gì so với lý thuyết dầm Timoshenko truyền thống? Lý thuyết truyền thống dùng độ võng và góc xoay làm ẩn số nên thường gặp hiện tượng khóa trượt và không tự động hội tụ về dầm Euler – Bernoulli khi mô đun trượt tiến ra vô cùng. Phương pháp của luận văn dùng độ võng và lực cắt làm ẩn, giúp bài toán hội tụ trơn tru và phản ánh chính xác phân bố ứng suất tiếp.

Hệ số phân bố ứng suất trượt trong luận văn được xác định bằng bao nhiêu? Đối với tiết diện hình chữ nhật chịu uốn ngang phẳng, tác giả sử dụng hệ số phân bố ứng suất trượt tại trục trung hòa bằng 1.20. Giá trị này được xác lập dựa trên nguyên lý cân bằng công của ứng suất cắt thực hiện trên toàn bộ chiều cao tiết diện thanh.

Nội lực mô men uốn tại nút khung có thay đổi nhiều khi kể đến biến dạng trượt không? Các kết quả tính toán trên mô hình thử nghiệm cho thấy mô men uốn tại các nút khung chịu ảnh hưởng tương đối nhỏ, với mức chênh lệch chỉ dao động từ 0.50% đến 1.50% khi tỷ số kích thước thay đổi từ 1/10 đến 1/3. Biến dạng trượt chủ yếu tác động trực tiếp lên độ võng và lực cắt.

Phương pháp bán giải tích trong luận văn giải quyết hệ phương trình như thế nào? Tác giả thiết lập hàm xấp xỉ đa thức bậc 4 cho đường độ võng và đường lực cắt, sau đó cực trị hóa phiếm hàm Gauss kết hợp nhân tử Lagrange. Quy trình này dẫn đến hệ phương trình đại số tuyến tính gồm 24 phương trình cho khung hai thanh và 29 phương trình cho khung một tầng một nhịp.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp tính toán khung phẳng chịu uốn có xét đến biến dạng trượt ngang dựa trên nguyên lý cực trị Gauss và hệ hai hàm ẩn độ võng kết hợp lực cắt.
  • Khẳng định định lượng mức độ ảnh hưởng của lực cắt: làm tăng độ võng từ 27.77% đến 69.56% đối với các cấu kiện có tỷ số chiều cao trên nhịp từ 1/5 đến 1/3.
  • Chứng minh tính ổn định của mô men uốn tại nút khung khi mức độ biến thiên chỉ duy trì trong khoảng hẹp từ 0.50% đến 1.50%.
  • Cung cấp chương trình tính toán số hóa hoàn chỉnh trên máy tính điện tử, sẵn sàng chuyển giao cho công tác thiết kế thực tế các kết cấu dầm chuyển và cột vách.
  • Mở ra định hướng phát triển bài toán cho hệ khung không gian 3 chiều và bài toán dao động động lực học trong giai đoạn tiếp theo; các kỹ sư và nhà nghiên cứu quan tâm nên ứng dụng ngay giải pháp này vào các dự án công trình quy mô lớn hiện nay.