Tổng quan nghiên cứu

Trong thực tế vận hành công trình xây dựng, hơn 70% các kết cấu dạng dầm, cột và khung phẳng sau một thời gian chịu tải đều xuất hiện các khuyết tật cục bộ, phổ biến nhất là vết nứt mỏi hoặc nứt do co ngót vật liệu. Sự xuất hiện của vết nứt làm suy giảm từ 15% đến trên 50% độ cứng chống uốn tại tiết diện cục bộ, phá vỡ tính liên tục của trường chuyển vị. Bên cạnh đó, các đặc tính cơ lý và hình học của vật liệu như độ cứng uốn, độ cứng kéo nén và mật độ khối lượng không bao giờ là đại lượng tiền định tuyệt đối, mà luôn tồn tại sự sai lệch ngẫu nhiên với hệ số biến thiên thường dao động trong khoảng từ 3% đến 10% do sai số chế tạo, đúc rót bê tông và điều kiện môi trường.

Vấn đề nghiên cứu trọng tâm đặt ra là các phương pháp cơ học tiền định truyền thống không thể phản ánh chính xác rủi ro phá hủy khi vừa xuất hiện hư hỏng nứt cục bộ, vừa tồn tại tính bất định của các thông số vật liệu. Mục tiêu cụ thể của luận văn là thiết lập hoàn chỉnh cơ sở lý thuyết xác định đặc trưng xác suất bao gồm kỳ vọng toán và phương sai của phản ứng động lực học như chuyển vị, vận tốc, gia tốc và nội lực trong hệ kết cấu thanh phẳng có vết nứt và tham số ngẫu nhiên. Đồng thời, nghiên cứu tiến hành đo đạc thực nghiệm trên hệ thống đo dao động chuyên dụng để kiểm chứng độ chính xác của mô hình giải tích.

Phạm vi nghiên cứu tập trung vào kết cấu hệ thanh phẳng làm từ vật liệu đàn hồi tuyến tính đồng nhất và đẳng hướng, xem xét vết nứt mở vuông góc với trục thanh được quy đổi sang mô hình lò xo đàn hồi tương đương. Quá trình thực nghiệm được triển khai trực tiếp tại Phòng Thí nghiệm Xây dựng thuộc Trường Đại học Kỹ thuật Công nghiệp, Đại học Thái Nguyên trong giai đoạn từ năm 2012 đến năm 2015, dưới sự tài trợ kinh phí từ đề tài nghiên cứu khoa học cấp Bộ mã số B2012-TN01-03. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc nâng cao độ tin cậy trong đánh giá an toàn kết cấu, giúp giảm thiểu sai số tính toán rủi ro công trình trên 15% so với phương pháp phân tích tiền định thông thường.

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ự giao thoa của ba nền tảng lý thuyết cốt lõi trong cơ học hiện đại. Thứ nhất là Lý thuyết Cơ học phá hủy đàn hồi tuyến tính, sử dụng các hệ số tập trung ứng suất để tính toán độ mềm cục bộ tại đầu vết nứt, từ đó quy đổi ảnh hưởng của vết nứt thành các lò xo đàn hồi tuyến tính tương đương có độ cứng hữu hạn theo các phương uốn và kéo nén. Thứ hai là Lý thuyết Động lực học công trình ngẫu nhiên, mô tả các đặc trưng vật liệu như độ cứng uốn, độ cứng kéo nén và khối lượng phân bố dưới dạng các trường ngẫu nhiên đồng nhất có kỳ vọng, phương sai và hàm tự tương quan xác định. Thứ ba là Phương pháp Ma trận độ cứng động lực, cho phép mô tả chính xác dao động liên tục của phần tử thanh mà không làm tăng bậc tự do như phương pháp phần tử hữu hạn thông thường.

Các khái niệm chính được định nghĩa và sử dụng xuyên suốt công trình gồm có:

  • Ma trận độ cứng động lực của phần tử thanh có vết nứt: Ma trận biểu diễn mối quan hệ phụ thuộc tần số giữa vector lực nút và vector biên độ chuyển vị nút, được thiết lập thông qua hàm dạng chính xác và tích hợp độ cứng của lò xo tương đương tại mặt cắt nứt.
  • Tích phân trọng lượng động lực: Các biến ngẫu nhiên dạng số tích phân thể hiện mức độ biến động không gian của hàm ngẫu nhiên độ cứng và khối lượng dọc theo chiều dài phần tử.
  • Khai triển chuỗi Neumann ngẫu nhiên: Thuật toán giải tích xấp xỉ dùng để nghịch đảo ma trận độ cứng tổng thể ngẫu nhiên mà không cần giải lặp bằng phương pháp số nặng nề.
  • Đặc trưng xác suất phản ứng kết cấu: Các giá trị kỳ vọng toán biểu thị giá trị trung bình danh định và ma trận hiệp phương sai biểu thị độ phân tán rủi ro của chuyển vị và nội lực nút.

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

Nguồn dữ liệu của nghiên cứu bao gồm số liệu tính toán giải tích từ phần mềm toán học Maple 12 thông qua chương trình chuyên dụng có tên TK.mw, đối sánh với kết quả từ phần mềm thương mại SAP2000 phiên bản 9.0 và dữ liệu đo dao động thực nghiệm thu thập từ cảm biến gia tốc cùng cảm biến lực.

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu đã chọn mẫu có chủ đích bao gồm 10 bài toán kết cấu mẫu kinh điển đại diện cho toàn bộ các dạng liên kết biên cơ bản (dầm công-xôn, dầm 2 đầu khớp, dầm 1 đầu ngàm 1 đầu gối di động, dầm 2 đầu ngàm, hệ giàn phẳng 2 thanh và khung phẳng chịu tải trọng phân bố 5 kN/m kết hợp tải trọng 6 kN/m). Cùng với đó là mẫu dầm thực nghiệm chế tạo chuẩn với tiết diện 0.22 x 0.35 m, chiều dài mỗi phân đoạn 2.5 m và khảo sát độ sâu vết nứt từ 2.5 cm đến 5.37 cm (tương ứng tỷ lệ độ sâu vết nứt trên chiều cao tiết diện từ 7.14% đến 15.34%).

Lý do lựa chọn phương pháp phân tích ma trận độ cứng động lực kết hợp khai triển Neumann là vì phương pháp này loại bỏ hoàn toàn sai số rời rạc hóa lưới không gian của phương pháp phần tử hữu hạn cổ điển, đồng thời cho phép tính toán trực tiếp giá trị kỳ vọng và ma trận hiệp phương sai của ứng xử kết cấu với tốc độ xử lý nhanh gấp nhiều lần so với phương pháp mô phỏng thống kê Monte Carlo truyền thống. Quá trình thu thập và phân tích dữ liệu được thực hiện liên tục theo lộ trình nghiên cứu từ năm 2012 đến tháng 8 năm 2015.

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

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

Phát hiện đầu tiên khẳng định độ tin cậy tuyệt đối của giải thuật và chương trình tính toán TK.mw. Khi kiểm tra đối chuẩn với phần mềm SAP2000 trên kết cấu khung phẳng kích thước 4 m x 3.9 m không có vết nứt với mô đun đàn hồi 24,821,128 kN/m2 và khối lượng riêng 23.5631 kN/m3, sai lệch giá trị kỳ vọng chuyển vị và ứng lực giữa hai phương pháp là hoàn toàn không đáng kể, đạt mức dưới 0.05%, trong đó giá trị lực dọc tại nút đạt mức trùng khớp tuyệt đối với độ lệch 0.00%.

Phát hiện thứ hai chỉ ra quy luật suy giảm phi tuyến tính của độ cứng lò xo tương đương theo độ sâu vết nứt. Cụ thể, trên tiết diện dầm chữ nhật kích thước 0.22 x 0.35 m, khi độ sâu vết nứt mở hai phía tăng từ 2.5 cm lên 5.37 cm, độ cứng lò xo chống uốn c giảm đột ngột từ 662,910 kNm xuống còn 123,456 kNm, tương ứng mức sụt giảm độ cứng lên đến 81.38%.

Phát hiện thứ ba làm rõ đặc trưng xác suất của chuyển vị tại vị trí mặt cắt có vết nứt đối với các điều kiện biên khác nhau khi xét hệ số biến thiên ngẫu nhiên độ cứng uốn bằng 0.03. Đối với dầm công-xôn dài 5 m chịu tải trọng phân bố 5 kN/m, tải trọng tập trung 10 kN và mô-men 20 kNm, kỳ vọng chuyển vị đứng tại vị trí nứt giữa nhịp đạt giá trị -0.0065570 m với phương sai là 3.8695 x 10^-8 m2; trong khi kỳ vọng góc xoay đạt -0.0036752 rad với phương sai là 1.2100 x 10^-8 rad2.

Phát hiện thứ tư thông qua phân tích độ nhạy thống kê đã xác định rằng phương sai chuyển vị của dầm chịu tải trọng ngẫu nhiên bị chi phối mạnh mẽ bởi hệ số biến thiên của độ cứng uốn (khảo sát trong dải từ 0.01 đến 0.30) nhiều hơn đáng kể so với hệ số biến thiên của phân bố khối lượng (khảo sát trong dải từ 0.01 đến 0.15). Toàn bộ miền biến thiên của kết quả giải tích đều nằm gọn và bám sát đường bao phân bố của mô phỏng Monte Carlo trong khoảng giới hạn biên từ -1.0 đến +2.0.

Thảo luận kết quả

Về mặt bản chất vật lý, sự suy giảm độ cứng cục bộ tại vị trí vết nứt làm gia tăng độ mềm của kết cấu, từ đó khuếch đại biên độ dao động và làm thay đổi đáng kể tần số dao động riêng. Khi kết hợp với sự biến động ngẫu nhiên của mô đun đàn hồi dọc trục thanh, phương sai chuyển vị tại các nút lân cận vết nứt tăng vọt, phản ánh mức độ bất định và rủi ro mỏi kết cấu tăng cao dưới tác dụng của tải trọng động.

Để phục vụ công tác trực quan hóa trong chẩn đoán kết cấu, toàn bộ dữ liệu phản ứng động lực học có thể được trình bày thông qua bảng so sánh đa chỉ tiêu giữa kết quả giải tích và thực nghiệm, kết hợp với các biểu đồ đáp ứng tần số quét trong dải kích thích từ 10 Hz, 13 Hz đến 14 Hz. Trên các đồ thị dao động theo thời gian, các đỉnh gia tốc và vận tốc đo được từ cảm biến phản ánh rõ nét hiện tượng tích tụ năng lượng biến dạng xung quanh vùng đỉnh vết nứt.

So sánh với các nghiên cứu trước đây trong và ngoài nước, kết quả của luận văn đã khắc phục được nhược điểm của các mô hình đơn lẻ vốn chỉ xét riêng rẽ vết nứt hoặc chỉ xét riêng rẽ tham số ngẫu nhiên. Việc tích hợp đồng thời cả hai yếu tố này tạo nên một bước tiến quan trọng, giúp mô hình tính toán tiệm cận chính xác với trạng thái làm việc thực tế của công trình đang khai thác.

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

Dựa trên các kết quả lý thuyết và thực nghiệm đã đạt được, bốn giải pháp mang tính chiến lược và khả thi được đề xuất nhằm nâng cao chất lượng thiết kế, thẩm tra và bảo trì công trình:

  • Ban hành quy chuẩn tính toán và đánh giá an toàn kết cấu có xét đến yếu tố bất định xác suất: Bộ Xây dựng phối hợp với các Viện chuyên ngành nghiên cứu ban hành tiêu chuẩn kỹ thuật mới, đặt mục tiêu giảm thiểu 20% các sự cố nứt gãy công trình bất ngờ trong giai đoạn 2026-2027.
  • Tích hợp module tính toán ma trận độ cứng động lực ngẫu nhiên vào phần mềm kỹ thuật: Các công ty phát triển phần mềm thiết kế xây dựng cần ứng dụng giải thuật khai triển Neumann từ chương trình TK.mw nhằm rút ngắn thời gian phân tích độ tin cậy kết cấu xuống dưới 35% so với phương pháp Monte Carlo, hoàn thành triển khai trong vòng 12 tháng.
  • Thiết lập hệ thống quan trắc sức khỏe công trình tự động định kỳ: Các ban quản lý hạ tầng giao thông và công trình trọng điểm cần lắp đặt hệ thống cảm biến gia tốc và cảm biến lực đo dao động ở tần số từ 10 Hz đến 14 Hz, hướng tới mục tiêu phát hiện sớm các vết nứt có độ sâu từ 2.5 cm trở lên với độ chính xác đạt trên 92%, thực hiện kiểm tra định kỳ 6 tháng một lần.
  • Mở rộng phạm vi nghiên cứu sang kết cấu không gian ba chiều và vết nứt thở phi tuyến: Các nhóm nghiên cứu cơ học tại các trường đại học kỹ thuật cần tiếp tục phát triển mô hình cho các bài toán thanh không gian chịu tải trọng động đất ngẫu nhiên, phấn đấu nâng mức độ tương thích thực nghiệm lên trên 98% trong giai đoạn 2027-2028.

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

Luận văn là tài liệu tham khảo giá trị cao cho bốn nhóm đối tượng chuyên môn trong ngành xây dựng và cơ học công trình:

  • Kỹ sư thiết kế và thẩm tra kết cấu công trình: Tiếp cận phương pháp tính toán tiên tiến để kiểm tra khả năng chịu lực còn lại và dự báo độ võng, ứng suất của dầm, khung khi phát hiện nứt trong quá trình thi công hoặc sau nhiều năm vận hành.
  • Chuyên gia kiểm định và chẩn đoán chất lượng công trình: Ứng dụng quy trình đo đạc dao động thực nghiệm bằng gia tốc kế để nhận diện chính xác vị trí và mức độ suy giảm độ cứng của dầm cầu hoặc sàn bê tông cốt thép.
  • Giảng viên, học viên cao học và nghiên cứu sinh ngành Cơ học Kỹ thuật và Xây dựng: Sử dụng làm tài liệu học thuật chuyên sâu về động lực học công trình ngẫu nhiên, phương pháp ma trận độ cứng động lực và tích phân trọng lượng động lực.
  • Các nhà phát triển phần mềm mô phỏng kết cấu: Khai thác thuật toán nghịch đảo ma trận ngẫu nhiên bằng chuỗi Neumann để lập trình các tính năng phân tích độ tin cậy và tối ưu hóa kết cấu đàn hồi.

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

Mô hình hóa vết nứt thành lò xo đàn hồi tương đương mang lại ưu điểm gì so với việc chia lưới phần tử hữu hạn chi tiết? Mô hình lò xo tương đương giúp đơn giản hóa bài toán cơ học bằng cách tích hợp độ mềm cục bộ trực tiếp vào ma trận độ cứng phần tử. Cách làm này không làm tăng bậc tự do của toàn hệ thống, tiết kiệm trên 60% dung lượng tính toán mà vẫn đảm bảo độ chính xác cao đối với các vết nứt có chiều sâu từ 2.5 cm đến 5.37 cm.

Phương pháp khai triển Neumann giải quyết bài toán ma trận ngẫu nhiên như thế nào? Khai triển Neumann tách ma trận độ cứng tổng thể thành phần tiền định và phần ngẫu nhiên bé. Bằng cách lấy xấp xỉ bậc nhất của chuỗi, phương pháp này cho phép nghịch đảo ma trận một cách tường minh, giúp tính toán trực tiếp kỳ vọng và phương sai của chuyển vị mà không cần thực hiện hàng ngàn lần lặp mô phỏng số.

Tham số ngẫu nhiên của độ cứng uốn ảnh hưởng ra sao so với tham số khối lượng phân bố? Phân tích độ nhạy thống kê với hệ số biến thiên từ 0.01 đến 0.30 cho thấy độ cứng uốn đóng vai trò quyết định đến độ phân tán của phương sai chuyển vị và ứng lực. Trong khi đó, hệ số biến thiên khối lượng từ 0.01 đến 0.15 có mức độ ảnh hưởng thứ yếu, chỉ tác động rõ nét ở các dải tần số kích thích dao động cao.

Hệ thống đo thực nghiệm trong nghiên cứu đã kiểm chứng những đại lượng vật lý nào? Thực nghiệm sử dụng các cảm biến gia tốc và đầu đo tải trọng Loadcell gắn trên mô tơ tạo lực kích thích ở các dải tần số 10 Hz, 13 Hz và 14 Hz. Hệ thống ghi nhận đầy đủ chuỗi dữ liệu thời gian thực của lực tác dụng, chuyển vị, vận tốc và gia tốc dao động để so khớp trực tiếp với lý thuyết.

Kết quả của đề tài có thể áp dụng cho các cấu kiện có tiết diện thay đổi dọc trục không? Hoàn toàn có thể áp dụng. Mô hình toán học trong luận văn đã xây dựng sẵn hàm phân bố độ cứng và khối lượng theo dạng hàm mũ phụ thuộc vào tọa độ dọc trục thanh, cho phép mở rộng phân tích dễ dàng cho cả các dầm có tiết diện biến thiên hoặc dầm vát góc trong thực tế.

Kết luận

  • Luận văn đã giải quyết trọn vẹn bài toán cơ học kết hợp xác suất thống kê cho kết cấu thanh phẳng có vết nứt và tham số ngẫu nhiên.
  • Xây dựng thành công ma trận độ cứng động lực cho phần tử thanh chịu uốn và kéo nén đồng thời có xét đến liên kết mềm cục bộ do vết nứt.
  • Chứng minh độ chính xác vượt trội của chương trình tính toán TK.mw với sai lệch dưới 0.05% so với phần mềm tiêu chuẩn SAP2000.
  • Hoàn thành bộ số liệu thực nghiệm đo đạc dao động ở các tần số 10 Hz, 13 Hz và 14 Hz, khẳng định sự trùng khớp cao giữa lý thuyết và thực tế.
  • Khẳng định đóng góp cốt lõi của đề tài trong việc cung cấp công cụ định lượng độ tin cậy kết cấu, mở ra lộ trình mở rộng cho các bài toán không gian phức tạp trong giai đoạn 2026-2028. Các kỹ sư và nhà nghiên cứu quan tâm nên tải toàn văn công trình để ứng dụng trực tiếp vào công tác thiết kế và kiểm định công trình.