Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình hiện đại, hơn 75% các cấu trúc mái nhịp lớn, cầu vượt khẩu độ dài và tháp truyền không gian sử dụng hệ kết cấu dàn thép nhờ ưu điểm trọng lượng nhẹ cùng khả năng chịu lực tối ưu. Tuy nhiên, các phương pháp phân tích truyền thống thường giả thiết mắt dàn là các liên kết khớp lý tưởng và xem chiều dài tính toán của thanh bằng đúng khoảng cách giữa hai tâm nút. Trên thực tế thi công, các phần tử thép liên kết với nhau qua các bản mã hàn hoặc bu-lông tạo nên vùng giao nhau có độ cứng rất lớn, được định nghĩa là vùng cứng với kích thước chiếm khoảng 5% đến 10% chiều dài mỗi thanh. Đồng thời, dưới tác dụng của tải trọng lớn, các thanh chịu nén xuất hiện mômen uốn phụ do chuyển vị ngang gây ra bởi hiệu ứng phi tuyến hình học, làm suy giảm đáng kể độ cứng tổng thể.

Mục tiêu cốt lõi của đề tài là xây dựng mô hình toán học và thuật toán số hoàn chỉnh để phân tích dao động tự nhiên của hệ dàn phẳng, tích hợp đồng thời hai yếu tố then chốt: kích thước hình học thực của vùng cứng và ảnh hưởng của phi tuyến hình học bậc hai. Nghiên cứu tập trung giải quyết bài toán xác định chính xác các tần số dao động riêng, chu kỳ dao động và các dạng dao động tương ứng khi kết cấu chịu tải trọng thay đổi từ trạng thái ban đầu $0 \text{ kN}$ cho đến tải trọng tới hạn $1829,89 \text{ kN}$. Phạm vi đề tài được triển khai trong chuyên ngành Xây dựng Dân dụng và Công nghiệp tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh. Kết quả phân tích giúp nâng cao độ chính xác dự báo động lực học, giảm độ lệch chuẩn chuyển vị từ 8% đến 15% so với phương pháp cổ điển, đảm bảo tính ổn định và an toàn chịu lực cho công trình nhịp lớn.

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 hai lý thuyết cơ học kết cấu: Lý thuyết phương pháp phần tử hữu hạn mở rộng cho phần tử thanh dầm phẳng và Lý thuyết động lực học công trình kết hợp phân tích ổn định phi tuyến bậc hai. Mô hình tính toán tích hợp phần tử thanh có xét độ cứng dọc trục $EA$, độ cứng chống uốn $EI$, cùng các ma trận chuyển đổi hình học tại hai đầu liên kết có kích thước vùng cứng $L_{ei}$ và $L_{ej}$.

Ba khái niệm lý thuyết trung tâm chi phối toàn bộ mô hình tính bao gồm:

  1. Vùng cứng (Rigid Zone): Miền liên kết cục bộ tại các mắt dàn có độ cứng uốn và độ cứng dọc trục vô cùng lớn so với thân phần tử, làm giảm chiều dài tự do thực tế của thanh dầm.
  2. Phi tuyến hình học (Geometric Nonlinearity): Hiện tượng ma trận độ cứng của kết cấu thay đổi phi tuyến phụ thuộc vào trạng thái biến dạng và sự xuất hiện của lực dọc trục $N$, kích hoạt hiệu ứng mômen phụ làm thay đổi ứng xử động lực học.
  3. Dao động tự nhiên và tần số riêng ($\omega, f$): Trạng thái dao động tự do không cản của hệ nhiều bậc tự do, phản ánh đặc tính động lực học nội tại của kết cấu khi chịu kích động ngoại vi.

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

Nghiên cứu sử dụng phương pháp mô hình hóa số phần tử hữu hạn trực tiếp. Kết cấu dàn phẳng khảo sát gồm hệ thống 13 phần tử thanh với các thông số vật liệu chuẩn: mô đun đàn hồi $E = 2,15 \times 10^{11} \text{ kN/m}^2$, diện tích mặt cắt ngang $A = 0,0001 \text{ m}^2$, mômen quán tính $I = 8,333 \times 10^{-10} \text{ m}^4$ và khối lượng riêng vật liệu $\rho = 1800 \text{ kN/m}^3$. Kích thước vùng cứng $L_e$ được khảo sát biến thiên từ $0 \text{ m}$ đến $0,15 \text{ m}$.

Mô hình được rời rạc hóa thành các nút có 3 bậc tự do (chuyển vị ngang, chuyển vị đứng và góc xoay). Thuật toán số được thiết lập tự động hóa hoàn toàn trên môi trường tính toán Matlab để xây dựng ma trận độ cứng tổng thể $[K]$ và ma trận khối lượng nhất quán $[M]$. Lý do lựa chọn phương pháp này là nhằm kiểm soát trực tiếp việc hiệu chỉnh các ma trận phần tử chứa tham số phi tuyến và kích thước vùng cứng, giải quyết bài toán trị riêng tổng quát $[K] - \omega^2 [M] = 0$. Đồng thời, toàn bộ kết quả phân tích được đối chứng chéo với mô hình mô phỏng trên phần mềm SAP2000 nhằm kiểm định độ tin cậy và sai số tính toán. Quá trình thu thập dữ liệu và xử lý mô phỏng được thực hiện liên tục trong khoảng thời gian 12 tháng.

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

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

Phân tích số trên mô hình dàn phẳng mang lại bốn phát hiện quan trọng có ý nghĩa thực tiễn lớn trong cơ học công trình:

  • Tác động gia cường độ cứng của kích thước vùng cứng: Khi kích thước vùng cứng $L_e$ tăng dần từ $0 \text{ m}$ lên $0,15 \text{ m}$, chiều dài tính toán thực của phần tử giảm xuống, làm tăng độ cứng tổng thể của toàn hệ kết cấu. Dưới tác dụng của trọng lượng bản thân, chuyển vị đứng tại nút dàn giảm mạnh từ $0,00402 \text{ mm}$ (trường hợp khớp lý tưởng) về các giá trị nhỏ hơn tương ứng với độ lớn vùng cứng.
  • Quy luật suy giảm tần số dao động theo tải trọng: Khi tải trọng nút $P$ tăng dần từ $0 \text{ kN}$ qua mức $651,34 \text{ kN}$, $1305 \text{ kN}$ đến mức tải trọng tới hạn $1829,89 \text{ kN}$ (tại hệ số biến dạng $\delta = 1,4$), tần số dao động tự nhiên $f$ của hệ dàn suy giảm liên tục. Ở dạng dao động thứ 1 (mode cơ bản), tần số dao động giảm rõ rệt nhất với mức giảm trên 60%, phản ánh sự mềm hóa của kết cấu dưới ảnh hưởng của lực nén dọc trục.
  • Tính bất biến của các dạng dao động bậc cao: Khác với các mode dao động đầu, các dạng dao động bậc cao (từ mode thứ 5 trở lên) có tần số dao động biến thiên không đáng kể (mức thay đổi dưới 3%) khi tăng tải trọng ngoại lực, chứng minh hiệu ứng phi tuyến hình học tập trung chi phối chủ yếu các dạng mất ổn định uốn tổng thể ban đầu.
  • Ngưỡng phá hoại mất ổn định kết cấu: Tại hệ số vượt tải $\delta = 1,5$, định thức của ma trận độ cứng tiếp cận giá trị suy biến, hệ dàn hoàn toàn mất ổn định động lực học. Tải trọng tới hạn an toàn xác định được là $1829,89 \text{ kN}$.

Thảo luận kết quả

Nguyên nhân căn bản của các hiện tượng trên bắt nguồn từ cơ chế làm việc cơ học tại liên kết nút. Trong thực tế, sự hiện diện của bản mã dày tạo ra một vùng cứng không biến dạng ở hai đầu thanh, làm triệt tiêu một phần chiều dài uốn tự do. Khi tính toán kết cấu theo mô hình khớp truyền thống, độ võng và chu kỳ dao động bị phóng đại quá mức so với thực tế làm việc của công trình. Ngược lại, việc bỏ qua phi tuyến hình học khi thanh chịu nén lớn sẽ gây nguy hiểm nghiêm trọng vì bỏ sót mômen uốn bậc hai $P\text{-}\Delta$.

Để phục vụ công tác đối sánh trực quan, toàn bộ dữ liệu phân tích được hệ thống hóa qua hệ thống bảng biểu và đồ thị kỹ thuật. Đồ thị hàm truyền biểu diễn mối tương quan giữa tải trọng $P$ và tần số riêng $f$ cho thấy đường cong suy giảm dốc đứng ở mode 1 và tiệm cận nằm ngang ở các mode cao. Biểu đồ cột so sánh sai số lực dọc giữa mô hình vùng cứng và mô hình khớp lý tưởng chỉ ra mức chênh lệch nội lực dao động từ 5% đến 18% tại các thanh giàn chịu nén chính. Kết quả tính toán từ thuật toán Matlab hoàn toàn tương thích và tiệm cận chặt chẽ với kết quả phân tích từ phần mềm chuyên dụng SAP2000, khẳng định độ tin cậy tuyệt đối của giải thuật đã xây dựng.

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

Dựa trên các kết quả số và phân tích cơ học chuyên sâu, bốn khuyến nghị kỹ thuật cụ thể được đề xuất nhằm nâng cao chất lượng thiết kế và thi công kết cấu dàn thép:

  1. Hiệu chỉnh mô hình tính toán trong thiết kế dàn thép: Các kỹ sư kết cấu và đơn vị tư vấn thiết kế cần tích hợp kích thước vùng cứng từ 5% đến 10% chiều dài phần tử vào sơ đồ tính của các hệ dàn khẩu độ lớn trên $24 \text{ m}$, nhằm giảm thiểu từ 10% đến 15% khối lượng thép dư thừa do giả định khớp quá an toàn, áp dụng ngay trong giai đoạn 2025 - 2026.
  2. Bắt buộc thực hiện phân tích phi tuyến bậc hai cho cấu kiện nén cao: Các cơ quan quản lý tiêu chuẩn xây dựng cần bổ sung quy định tính toán kiểm tra hiệu ứng phi tuyến hình học $P\text{-}\Delta$ đối với tất cả các thanh dàn chịu lực nén vượt quá 50% tải trọng tới hạn Euler, nhằm loại trừ triệt để nguy cơ sụp đổ đột ngột do mất ổn định đàn hồi trước năm 2027.
  3. Ứng dụng công cụ tự động hóa tính toán liên kết: Các viện nghiên cứu và doanh nghiệp phần mềm xây dựng nên phát triển các plugin tích hợp thuật toán ma trận vùng cứng trên các phần mềm thiết kế thông dụng, giúp rút ngắn 30% thời gian thiết lập mô hình và nâng cao độ chính xác kiểm tra dao động công trình trong vòng 18 tháng tới.
  4. Tăng cường kiểm định thực nghiệm dao động công trình: Các nhà thầu thi công và đơn vị thẩm tra cần triển khai đo đạc thực tế tần số dao động riêng của hệ dàn mái sau khi lắp dựng tại hiện trường, duy trì sai số cho phép dưới 5% so với mô hình mô phỏng số để nghiệm thu chất lượng liên kết nút.

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

Công trình nghiên cứu này mang lại giá trị học thuật và thực tiễn cao cho nhiều nhóm đối tượng chuyên môn:

  • Kỹ sư thiết kế kết cấu công trình: Nắm vững phương pháp thiết lập sơ đồ tính chính xác, lựa chọn kích thước bản mã và kiểm soát chuyển vị động học cho các công trình nhà xưởng nhịp lớn, giàn không gian và kết cấu cầu thép.
  • Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Xây dựng: Tiếp cận hệ thống công thức giải tích hoàn chỉnh về ma trận độ cứng và khối lượng mở rộng, sử dụng làm tài liệu tham khảo chất lượng để phát triển các thuật toán mô phỏng kết cấu phức tạp hơn.
  • Giảng viên và nhà nghiên cứu cơ học kết cấu: Sử dụng mô hình toán học và các ví dụ số chi tiết trong luận văn làm học liệu giảng dạy chuyên đề Động lực học công trình nâng cao và Phương pháp phần tử hữu hạn ứng dụng.
  • Chuyên gia thẩm tra và quản lý chất lượng dự án: Khai thác các chỉ số tới hạn về tải trọng và tần số riêng để đánh giá mức độ an toàn chịu lực của các công trình công nghiệp trọng điểm trước tải trọng gió động và động đất.

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

Tại sao cần phải xét kích thước vùng cứng thay vì xem nút dàn là khớp lý tưởng?
Trong thực tế thi công, các thanh dàn liên kết với nhau qua bản mã hàn hoặc bu-lông có độ cứng lớn, tạo ra vùng chồng lấn không biến dạng chiếm khoảng 5% đến 10% chiều dài thanh. Việc xét vùng cứng giúp phản ánh chính xác chiều dài làm việc tự do, giảm sai số tính toán chuyển vị đứng từ $0,00402 \text{ mm}$ về mức chuẩn xác hơn.

Phi tuyến hình học ảnh hưởng như thế nào đến tần số dao động riêng của hệ dàn?
Phi tuyến hình học tạo ra lực dọc trục nén lớn, sinh ra mômen uốn phụ làm suy giảm độ cứng chống uốn tổng thể của hệ thanh. Điều này làm tần số dao động riêng ở mode 1 suy giảm mạnh hơn 60% khi tải trọng tăng từ $0 \text{ kN}$ lên mức tới hạn $1829,89 \text{ kN}$.

Tải trọng tới hạn của hệ dàn được xác định bằng phương pháp nào trong nghiên cứu?
Tải trọng tới hạn được xác định thông qua phương pháp quét tham số tải trọng trên mô hình phần tử hữu hạn Matlab kết hợp tiêu chuẩn suy biến ma trận độ cứng tiếp tuyến. Kết quả chỉ ra hệ dàn mất ổn định hoàn toàn tại hệ số vượt tải $\delta = 1,5$, đạt tải trọng tới hạn an toàn ở mức $1829,89 \text{ kN}$ tại $\delta = 1,4$.

Kích thước vùng cứng làm thay đổi nội lực trong thanh dàn như thế nào?
Khi xét đến kích thước vùng cứng, tại hai đầu thanh xuất hiện mômen uốn cục bộ thay vì chỉ có lực dọc thuần túy như giả định khớp. Chênh lệch nội lực giữa hai mô hình đạt từ 5% đến 18% tùy thuộc vào vị trí thanh biên hay thanh bụng trong kết cấu.

Thuật toán trong luận văn có thể tích hợp vào các phần mềm kết cấu thương mại không?
Thuật toán được xây dựng hoàn toàn dựa trên phương pháp phần tử hữu hạn chuẩn tắc nên rất thuận tiện để lập trình nhúng vào các phần mềm như SAP2000, ETABS hoặc ANSYS thông qua các giao diện lập trình ứng dụng mở rộng nhằm tự động hóa quy trình phân tích động lực học kết cấu thực tế.

Kết luận

  • Đã thiết lập thành công hệ phương trình phần tử hữu hạn và ma trận độ cứng, ma trận khối lượng nhất quán tích hợp đồng thời kích thước vùng cứng và phi tuyến hình học cho hệ dàn phẳng.
  • Xây dựng hoàn chỉnh chương trình tính toán tự động trên Matlab, giải quyết bài toán dao động riêng với độ chính xác cao và đã được kiểm chứng tương đương với kết quả trên phần mềm SAP2000.
  • Chứng minh định lượng sự suy giảm tần số riêng ở các mode cơ bản (trên 60%) khi tải trọng tiến sát tải tới hạn $1829,89 \text{ kN}$, trong khi các mode bậc cao ít bị ảnh hưởng.
  • Khẳng định vai trò của việc mô hình hóa vùng cứng giúp tối ưu hóa thiết kế, phản ánh sát thực tế ứng xử chịu lực và nâng cao độ an toàn cho các công trình dàn thép vượt nhịp.
  • Hướng phát triển tiếp theo là mở rộng nghiên cứu cho hệ dàn không gian 3D chịu tác động của tải trọng động lực học ngẫu nhiên như động đất và gió giật trong giai đoạn tới; các đơn vị tư vấn thiết kế nên áp dụng ngay quy trình tính toán này vào các dự án kết cấu thép nhịp lớn hiện nay.