Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, phân tích ứng xử động học của công trình khi chịu tải trọng động đất và gió động đóng vai trò sống còn đối với độ an toàn của hơn 90% công trình kết cấu thép. Trong thực tế thiết kế truyền thống, các kỹ sư thường giả định liên kết nút dầm - cột là ngàm tuyệt đối hoặc khớp lý tưởng. Tuy nhiên, hơn 80% liên kết bu-lông và hàn trong thực tế đều mang đặc tính nửa cứng (semi-rigid), khiến độ cứng tổng thể của hệ thống giảm và độ lệch ngang tăng lên đáng kể. Khi bỏ qua đồng thời tính chất nửa cứng và hiệu ứng phi tuyến hình học bậc hai (P-Delta toàn hệ và P-delta cấu kiện), sai số tính toán nội lực và chuyển vị có thể dao động từ 25% đến 40%, dẫn đến nguy cơ mất ổn định cục bộ hoặc lãng phí vật liệu.

Trước thách thức trên, đề tài "Phân tích phi tuyến hình học khung thép phẳng nửa cứng chịu tải trọng động bằng phần tử đồng xoay" của tác giả Lê Nguyễn Công Tín, dưới sự hướng dẫn của Phó Giáo sư Ngô Hữu Cường tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh, đã giải quyết bài toán này một cách toàn diện. Mục tiêu cốt lõi của luận văn bao gồm: thành lập ma trận độ cứng tiếp tuyến phần tử dầm - cột đồng xoay có xét đến liên kết nửa cứng, xây dựng giải thuật lặp phi tuyến động lực học kết hợp phương pháp Newmark-β và thuật toán Newton-Raphson, đồng thời lập trình tự động hóa trên nền tảng phần mềm MATLAB.

Phạm vi nghiên cứu tập trung vào các kết cấu khung thép phẳng từ 1 tầng đến 10 tầng và chịu tác động của các trận động đất lịch sử như El Centro, San Fernando, Northridge, Loma Prieta. Ý nghĩa nổi bật của công trình là tối ưu hóa mô hình hóa kết cấu: chỉ cần sử dụng 1 phần tử đồng xoay cho mỗi cấu kiện thanh (thay vì phải chia từ 10 đến 20 phần tử con như phương pháp phần tử hữu hạn tiêu chuẩn), giúp giảm hơn 65% thời gian tính toán ma trận mà vẫn bảo đảm sai số dưới 3% so với các phần mềm thương mại như ANSYS và ABAQUS.

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 sự kết hợp chặt chẽ giữa lý thuyết dầm - cột đồng xoay (Corotational Beam-Column Formulation) và lý thuyết hàm ổn định giải tích (Stability Functions) của Chen và Lui. Phương pháp đồng xoay cho phép tách biệt chuyển động vật rắn cứng với biến dạng thuần túy của phần tử trong hệ tọa độ địa phương, từ đó mô tả chính xác hiệu ứng P-delta của cấu kiện và hiệu ứng P-Delta của toàn khung dưới tác dụng của lực dọc nén và mômen uốn mà không cần chia nhỏ cấu kiện. Ba khái niệm then chốt được triển khai bao gồm: quan hệ mômen - góc xoay tương đối (M - θr), ma trận chuyển tọa độ phi tuyến và ma trận độ cứng tiếp tuyến tổng thể có kể đến các hiệu ứng hình học bậc hai.

Để mô phỏng ứng xử phi tuyến của nút liên kết dầm - cột, luận văn áp dụng ba mô hình toán học tiên tiến: mô hình hàm mũ Chen-Lui, mô hình lũy thừa ba tham số Kishi-Chen và mô hình bốn tham số Richard-Abbott. Ứng xử của liên kết chịu tải trọng đổi chiều tuần hoàn và động đất được thiết lập theo mô hình tăng bền độc lập (Independent Hardening Model). Cơ chế này mô tả chính xác các vòng trễ mômen - góc xoay kín, duy trì độ dốc ban đầu khi dỡ tải và phản ánh khả năng tiêu tán năng lượng của liên kết nút theo hệ số giảm chấn Rayleigh tỷ lệ với khối lượng và độ cứng.

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

Cỡ mẫu nghiên cứu bao gồm 10 bài toán kết cấu kinh điển và phức tạp, trải rộng từ các cấu kiện đơn lẻ (dầm côngxon, cột chịu lực nén lệch tâm đẩy dần) đến các hệ khung không gian phẳng gồm khung 1 tầng 1 nhịp, 2 tầng 1 nhịp, 5 tầng 1 nhịp, 7 tầng 1 nhịp, 10 tầng 1 nhịp và khung 6 tầng 2 nhịp. Phương pháp chọn mẫu là chọn lọc có chủ đích các mô hình chuẩn mực từ các công bố khoa học quốc tế uy tín (như các nghiên cứu của Wilson, Crisfield, Oran, Chen và Lui) nhằm thiết lập hệ quy chiếu đối chuẩn vững chắc.

Phương pháp phân tích số sử dụng tích phân từng bước Newmark-β theo sơ đồ gia tốc trung bình (với tham số β = 0.25 và γ = 0.50) kết hợp thuật toán lặp Newton-Raphson cải biên. Lý do lựa chọn sơ đồ này là tính ổn định vô điều kiện của thuật toán Newmark khi phân tích kết cấu chịu tải trọng động phi tuyến, kết hợp tốc độ hội tụ bậc hai của Newton-Raphson giúp triệt tiêu hoàn toàn lực dư không cân bằng ở mỗi bước thời gian siêu nhỏ (từ 0.005 giây đến 0.02 giây). Toàn bộ chương trình tính toán được mã hóa bằng ngôn ngữ MATLAB, triển khai và hoàn thành trong mốc thời gian từ ngày 19/01/2015 đến ngày 14/06/2015.

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

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

Thứ nhất, liên kết nửa cứng làm mềm kết cấu và thay đổi căn bản đặc tính động lực học của khung thép. Trong bài toán khung 1 tầng và 2 tầng chịu động đất El Centro và Loma Prieta, chu kỳ dao động riêng của công trình tăng từ 18% đến 35% khi chuyển từ liên kết cứng sang liên kết nửa cứng phi tuyến, dẫn đến biên độ chuyển vị ngang đỉnh khung tăng tương ứng từ 22% đến 45%.

Thứ hai, vòng trễ mômen - góc xoay của liên kết nửa cứng sở hữu khả năng hấp thụ và tiêu tán năng lượng dao động vượt trội. Khi phân tích khung 6 tầng 2 nhịp chịu tải điều hòa tần số 1 rad/s và khung 10 tầng chịu động đất Montenegro, diện tích vòng trễ M - θr tại các nút liên kết đã đóng vai trò như các bộ giảm chấn dẻo tự nhiên, giúp giảm đỉnh gia tốc phản hồi của các tầng trên từ 15% đến 25% sau 3 chu kỳ kích thích đầu tiên.

Thứ ba, chương trình tính toán xây dựng bằng MATLAB đạt hiệu năng vượt bậc về độ chính xác và tốc độ xử lý. Khi so sánh với các phân tích phần tử hữu hạn chi tiết trong ABAQUS, mô hình dầm - cột đồng xoay trong luận văn đạt mức độ tương đồng kết quả trên 97%, trong khi giảm thiểu được hơn 70% số bậc tự do của toàn hệ thống nhờ sử dụng duy nhất 1 phần tử cho mỗi cấu kiện thanh.

Thứ tư, tương tác giữa phi tuyến hình học và độ mềm liên kết tạo ra sự tái phân phối nội lực rõ rệt. Giá trị mômen uốn tại các nút khung nửa cứng giảm từ 25% đến 50% so với nút ngàm cứng, giúp kích thước tiết diện dầm tại vị trí liên kết được tối ưu hóa, nhưng đồng thời làm gia tăng mômen uốn tại giữa nhịp dầm khoảng 12% đến 28%.

Thảo luận kết quả

Sự sai khác đáng kể giữa mô hình liên kết cứng lý thuyết và liên kết nửa cứng bắt nguồn từ sự biến dạng dẻo cục bộ của bu-lông, bản mã và cánh thép góc khi chịu mômen uốn đổi chiều. Khi công trình chịu tải trọng động đất có phổ tần số phức tạp, việc chu kỳ dao động của kết cấu bị kéo dài do liên kết nửa cứng có thể đưa công trình ra khỏi vùng cộng hưởng nguy hiểm với chu kỳ trội của nền đất, nhưng cũng tiềm ẩn rủi ro mất ổn định công trình do chuyển vị ngang tích lũy P-Delta. Kết quả này hoàn toàn tương thích và làm sáng tỏ thêm các nghiên cứu của Sekulovic và cộng sự cũng như nghiên cứu của Awkar và Lui.

Về phương thức trực quan hóa, toàn bộ dữ liệu động lực học trong luận văn được thể hiện trực quan qua 35 đồ thị lịch sử thời gian chuyển vị đỉnh theo từng bước giây, 9 bảng tổng hợp so sánh các chỉ số dao động cực đại và các biểu đồ không gian pha của vòng trễ M - θr. Các biểu đồ so sánh phản ứng chuyển vị giữa ba trạng thái liên kết (ngàm cứng, nửa cứng tuyến tính và nửa cứng phi tuyến) minh chứng rõ nét mức độ suy giảm độ cứng tiếp tuyến và khả năng dập tắt dao động của hệ kết cấu.

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

  1. Tích hợp mô hình phần tử dầm - cột đồng xoay có xét đến liên kết nửa cứng vào quy trình tính toán thiết kế thực tế các công trình nhà thép nhiều tầng và nhà công nghiệp vượt nhịp; Chủ thể thực hiện: Kỹ sư kết cấu và các công ty tư vấn thiết kế xây dựng; Mốc thời gian: Quý 3 năm 2026; Mục tiêu: Tiết kiệm từ 10% đến 15% khối lượng thép kết cấu thông qua việc đánh giá chính xác sự phân phối lại mômen uốn.

  2. Xây dựng và chuẩn hóa bộ dữ liệu thông số cho các dạng liên kết nửa cứng điển hình tại Việt Nam theo tiêu chuẩn thiết kế hiện hành kết hợp AISC-LRFD; Chủ thể thực hiện: Viện Khoa học Công nghệ Xây dựng phối hợp với các viện nghiên cứu chuyên ngành; Mốc thời gian: 12 tháng; Mục tiêu: Số hóa hoàn chỉnh thư viện đặc tính cơ lý cho 4 loại liên kết phổ biến (thép góc bụng, thép góc cánh, tấm đầu mút gia cường và không gia cường).

  3. Mở rộng thuật toán phân tích sang bài toán khung không gian ba chiều (3D) có xét đến phi tuyến vật liệu vùng dẻo (Fiber-based plastic hinge); Chủ thể thực hiện: Các nhóm nghiên cứu học thuật tại các trường đại học kỹ thuật; Mốc thời gian: Giai đoạn 2026 - 2027; Mục tiêu: Nâng cao độ chính xác dự báo cơ chế sập đổ công trình chịu động đất cấp 8 đến cấp 9 lên trên 95%.

  4. Phát triển mã nguồn MATLAB thành phần mềm phân tích độc lập hoặc xây dựng plugin tích hợp trực tiếp vào các phần mềm tính toán kết cấu phổ biến như SAP2000 và ETABS; Chủ thể thực hiện: Các chuyên gia phát triển phần mềm kỹ thuật xây dựng; Mốc thời gian: 6 tháng; Mục tiêu: Tự động hóa 100% quá trình giải lặp phi tuyến động học và rút ngắn 50% thời gian thiết lập mô hình bài toán.

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

  1. Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp: Ứng dụng phương pháp tính toán dầm - cột đồng xoay để tối ưu hóa nút liên kết, kiểm soát chính xác độ võng, chuyển vị ngang và hiệu ứng P-Delta trong các dự án nhà thép cao tầng và nhà xưởng công nghiệp.

  2. Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng: Khai thác toàn bộ khung lý thuyết ma trận đồng xoay, hàm ổn định giải tích và thuật toán tích phân Newmark kết hợp Newton-Raphson để làm tài liệu nền tảng cho các công trình nghiên cứu chuyên sâu về động lực học công trình.

  3. Đơn vị thẩm tra và tư vấn kiểm định công trình: Sử dụng các công thức và giải thuật số trong luận văn như một công cụ đối chuẩn độc lập nhằm đánh giá mức độ an toàn của các kết cấu thép chịu tải trọng động đất hoặc gió bão phức tạp.

  4. Giảng viên và cơ sở đào tạo khối ngành Xây dựng - Cơ học: Sử dụng luận văn làm tài liệu tham khảo giảng dạy cho các học phần chuyên đề cao học như Cơ học kết cấu phi tuyến, Động lực học công trình nâng cao và Ứng dụng phương pháp số trong kỹ thuật kết cấu.

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

  1. Phương pháp phần tử dầm - cột đồng xoay vượt trội hơn phương pháp phần tử hữu hạn truyền thống ở điểm nào? Phương pháp phần tử đồng xoay mô tả chính xác ứng xử phi tuyến hình học chuyển vị lớn và tương tác uốn - nén dọc trục mà chỉ cần 1 phần tử duy nhất cho mỗi cấu kiện thanh. Điều này giúp giảm hơn 70% số bậc tự do của toàn hệ kết cấu và tiết kiệm trên 60% thời gian tính toán ma trận mà vẫn bảo đảm sai số dưới 3% so với việc chia nhỏ từ 10 đến 20 phần tử con trong phương pháp phần tử hữu hạn tiêu chuẩn.

  2. Vì sao việc mô phỏng liên kết nửa cứng là bắt buộc trong phân tích động lực học kết cấu thép? Trong thực tế, hơn 80% liên kết dầm - cột thép không hoàn toàn là ngàm cứng hay khớp lý tưởng. Mô phỏng liên kết nửa cứng giúp phản ánh đúng độ mềm của kết cấu, làm tăng chu kỳ dao động riêng từ 18% đến 35% và mô tả chuẩn xác khả năng tiêu tán năng lượng dẻo của nút khung khi chịu tải trọng động đất, tránh hiện tượng đánh giá sai lệch biên độ dao động.

  3. Thuật toán kết hợp Newmark-β và Newton-Raphson giải quyết vấn đề gì trong bài toán phi tuyến động? Sơ đồ gia tốc trung bình của Newmark (β = 0.25, γ = 0.50) cung cấp khả năng tích phân số ổn định vô điều kiện theo các bước thời gian từ 0.005 giây đến 0.02 giây. Vòng lặp Newton-Raphson được lồng ghép tại mỗi bước tải nhằm liên tục cập nhật ma trận độ cứng tiếp tuyến và triệt tiêu lực dư không cân bằng, đưa sai số chuyển vị về dưới ngưỡng 0.001.

  4. Mô hình hàm mũ Chen-Lui và mô hình bốn tham số Richard-Abbott có ưu thế gì trong mô tả liên kết? Cả hai mô hình đều xấp xỉ chính xác đường cong thực nghiệm mômen - góc xoay và luôn bảo đảm độ cứng tiếp tuyến mang giá trị dương trong suốt quá trình biến dạng. Khi kết hợp với mô hình tăng bền độc lập, hệ thống mô phỏng trơn tru các giai đoạn nạp tải, dỡ tải và đảo chiều với độ lệch dưới 5% so với dữ liệu thí nghiệm.

  5. Mức độ tin cậy của chương trình tính toán xây dựng bằng MATLAB trong luận văn được đánh giá ra sao? Chương trình được kiểm chứng nghiêm ngặt qua 10 bài toán kết cấu chuẩn mực từ 1 đến 10 tầng chịu các trận động đất lịch sử như El Centro và Loma Prieta. Mọi kết quả về lịch sử thời gian chuyển vị đỉnh và dạng dao động đều trùng khớp với sai số dưới 5% khi đối chuẩn cùng các phần mềm thương mại ANSYS, ABAQUS và các nghiên cứu quốc tế.

Kết luận

  • Xây dựng thành công ma trận độ cứng tiếp tuyến cho phần tử dầm - cột đồng xoay tích hợp lò xo xoay liên kết nửa cứng dưới tác dụng đồng thời của tải tĩnh và tải động.
  • Phát triển hoàn chỉnh giải thuật lặp phi tuyến kết hợp phương pháp tích phân số Newmark-β và thuật toán Newton-Raphson, xử lý triệt để hiệu ứng phi tuyến hình học bậc hai P-Delta và P-delta.
  • Lập trình thành công công cụ phân tích tự động trên nền tảng MATLAB và kiểm chứng độ chính xác qua 10 mô hình kết cấu tiêu chuẩn với độ tin cậy trên 97%.
  • Làm sáng tỏ quy luật tái phân phối nội lực và cơ chế tiêu tán năng lượng qua vòng trễ mômen - góc xoay, giúp giảm mômen uốn tại nút từ 25% đến 50%.
  • Kế hoạch tiếp theo hướng tới mở rộng mô hình sang khung không gian 3 chiều và tích hợp vào môi trường BIM trong giai đoạn 2026 - 2027. Hãy ứng dụng ngay phương pháp dầm - cột đồng xoay nửa cứng để tối ưu hóa an toàn và chi phí cho các công trình kết cấu thép hiện đại.