Tổng quan nghiên cứu

Trong quá trình khai thác và vận hành công trình xây dựng, sự suy giảm cơ lý tính và xuất hiện các dạng hư hỏng cục bộ như vết nứt là hiện tượng tất yếu. Theo các khảo sát kỹ thuật, phần lớn các kết cấu dầm, cột và bản sàn đều tiềm ẩn khuyết tật nứt, làm gián đoạn tính liên tục và giảm cục bộ độ cứng chống uốn, chống cắt của cấu kiện. Đa số các mô hình tính toán truyền thống đều giả định vật liệu hoàn toàn đồng nhất và tải trọng có tính tiền định, điều này dễ dẫn đến việc đánh giá sai lệch biên độ an toàn thực tế. Trên thực tế, các thông số vật liệu như độ cứng uốn $EI_0 = 24.821.128 \text{ kN/m}^2$, khối lượng phân bố $m_0$ với khối lượng riêng $23{,}5631 \text{ kN/m}^3$, cũng như kích thước tiết diện luôn tồn tại độ phân tán ngẫu nhiên dưới tác động của điều kiện thi công và môi trường.

Nghiên cứu tập trung giải quyết bài toán phân tích ứng xử động lực học ngẫu nhiên trong hệ kết cấu thanh phẳng có vết nứt, được thực hiện trong giai đoạn từ năm 2012 đến 2015 tại Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên theo khuôn khổ đề tài nghiên cứu mã số B2012-TN01-03. Mục tiêu cụ thể là xác định đặc trưng xác suất của các phản ứng kết cấu bao gồm giá trị kỳ vọng toán và phương sai của chuyển vị, vận tốc, gia tốc và nội lực khi kể đến sự biến thiên ngẫu nhiên của độ cứng uốn $EI(x)$, độ cứng kéo nén $AE(x)$, phân bố khối lượng $m(x)$ và tải trọng tác dụng.

Ý nghĩa của nghiên cứu thể hiện qua việc cung cấp công cụ toán học tin cậy giúp giảm thiểu sai số dự báo chuyển vị xuống dưới $0{,}01%$ so với nghiệm giải tích, đồng thời định lượng hóa mức độ rủi ro hư hỏng của các vết nứt có chiều sâu từ $2{,}5 \text{ cm}$ đến $5{,}37 \text{ cm}$ trên cấu kiện thự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 sự kết hợp giữa cơ học phá hủy đàn hồi tuyến tính, động lực học công trình và lý thuyết xác suất thống kê nâng cao:

  • Mô hình vết nứt tương đương: Sử dụng lý thuyết cơ học phá hủy của Thomson và Aliabadi để quy đổi trạng thái tập trung ứng suất tại vết nứt ngang thành hệ lò xo đàn hồi tuyến tính cục bộ. Độ mềm cục bộ $c$ được biểu diễn qua tích phân hệ số cường độ ứng suất $K_I$, liên kết trực tiếp với độ sâu vết nứt $a$ và chiều cao tiết diện $h$.
  • Phương pháp ma trận độ cứng động lực (Dynamic Stiffness Matrix - MTĐCĐL): Dựa trên phương trình vi phân dao động uốn và kéo nén đồng thời của thanh Bernoulli-Euler, thiết lập hàm dạng chính xác chứa các hàm lượng giác và hàm hyperbol thay vì các hàm đa thức xấp xỉ bậc thấp.
  • Lý thuyết trường ngẫu nhiên và khai triển Neumann: Mô hình hóa các đại lượng cơ lý tính dưới dạng hàm ngẫu nhiên đồng nhất $EI(x) = EI_0(1 + \epsilon_1 g_1(x))$ và $m(x) = m_0(1 + \epsilon_2 g_2(x))$ với kỳ vọng bằng $0$ và hàm tự tương quan $R_{ii}(\tau) = \exp(-\tau^2 / d^2)$. Phương pháp khai triển Neumann được áp dụng để nghịch đảo ma trận độ cứng ngẫu nhiên bậc cao mà không cần giải trực tiếp hệ phương trình vi phân ngẫu nhiên phi tuyến phức tạp.

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

Nghiên cứu kết hợp chặt chẽ giữa mô phỏng số giải tích và thực nghiệm đối chứng:

  • Phương pháp chọn mẫu và cỡ mẫu: Thu thập và phân tích dữ liệu thực nghiệm trên $12$ kịch bản dao động khác nhau của mẫu dầm thép hình học chữ nhật $(0{,}22 \times 0{,}35 \text{ m})$ với chiều dài nhịp $L = 5 \text{ m}$, bao gồm mẫu dầm nguyên vẹn và mẫu dầm có vết nứt nhân tạo với tỷ số chiều sâu vết nứt $a/h$ biến thiên từ $0{,}07$ đến $0{,}15$. Phương pháp chọn mẫu có chủ đích được áp dụng nhằm tập trung vào các trạng thái giới hạn nguy hiểm nhất khi kết cấu chịu kích động điều hòa.
  • Phương pháp phân tích số: Xây dựng thuật toán chuyên dụng TK.mw trên môi trường Maple 12, giải hệ phương trình động lực ngẫu nhiên với $9$ bậc tự do nút. Lý do lựa chọn Maple là khả năng xử lý đại số ký hiệu chính xác các tích phân trọng lượng động lực $W_{kr}^u$ và cho phép so sánh kiểm chứng trực tiếp với phần mềm thương mại SAP2000 v9 cùng phương pháp mô phỏng Monte-Carlo với $10.000$ mẫu thử.
  • Quy trình thực nghiệm: Sử dụng hệ thống thiết bị đo động gồm các cảm biến gia tốc trở kháng cao, cảm biến lực Loadcell gắn trực tiếp tại nguồn kích rung mô tơ biến tần, và bộ chuyển đổi tín hiệu số kết nối máy tính để ghi nhận phản ứng dao động tại các dải tần số kích thích $10 \text{ Hz}$, $13 \text{ Hz}$ và $14 \text{ Hz}$.

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

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

  1. Quy luật suy giảm độ cứng theo kích thước vết nứt: Khi chiều sâu vết nứt $a$ tăng từ $2{,}5 \text{ cm}$ lên $5{,}37 \text{ cm}$ trên dầm tiết diện $(0{,}22 \times 0{,}35 \text{ m})$, độ cứng lò xo xoay tương đương $c$ giảm mạnh từ $662.910 \text{ kNm}$ xuống còn $123.456 \text{ kNm}$, tương đương mức suy giảm độ cứng cục bộ lên tới $81{,}37%$. Sự sụt giảm này làm thay đổi cục bộ phân bố góc xoay và biên độ dao động tại nút lân cận.
  2. Ảnh hưởng rõ rệt của điều kiện biên đến đặc trưng thống kê: Kỳ vọng và phương sai chuyển vị võng tại giữa nhịp chịu tác động sâu sắc từ dạng liên kết. Đối với dầm conson, kỳ vọng chuyển vị đạt $-0{,}006557 \text{ m}$ với phương sai $3{,}8695 \times 10^{-8} \text{ m}^2$; trong khi ở dầm hai đầu ngàm, kỳ vọng chuyển vị chỉ còn $-0{,}0008561 \text{ m}$ với phương sai $4{,}9795 \times 10^{-10} \text{ m}^2$, tức giảm hơn $86{,}94%$ về độ võng và giảm gần $98{,}71%$ về phương sai độ phân tán.
  3. Mức độ nhạy cảm của chuyển vị đối với các tham số ngẫu nhiên: Kết quả phân tích độ nhạy chứng minh rằng phương sai chuyển vị phụ thuộc chủ yếu vào hệ số biến thiên của mô đun đàn hồi ($\epsilon_1$) hơn là hệ số biến thiên của tải trọng ($\epsilon_2$). Khi $\epsilon_1$ thay đổi từ $0{,}01$ đến $0{,}30$, biên độ phương sai chuyển vị tiệm cận đường biên trên $+2{,}0$ của mô phỏng Monte-Carlo, trong khi sự biến thiên của $\epsilon_2$ chỉ tạo ra độ lệch giới hạn trong khoảng $[-0{,}4; +0{,}4]$.
  4. Độ tương đồng cao giữa thực nghiệm và mô hình lý thuyết: Dữ liệu gia tốc và vận tốc đo tại các tần số kích thích $10 \text{ Hz}$, $13 \text{ Hz}$ và $14 \text{ Hz}$ cho thấy dạng sóng dao động đo được bằng Loadcell và gia tốc kế trùng khớp chặt chẽ với đường cong đáp ứng tần số tính toán lý thuyết, sai số biên độ đỉnh dao động duy trì dưới $5{,}8%$.

Thảo luận kết quả

Các kết quả nghiên cứu chỉ ra rằng vết nứt đóng vai trò như một bộ lọc làm suy giảm tần số dao động riêng và gia tăng độ bất định của phản ứng kết cấu. Trên các đồ thị so sánh phân bố xác suất chuyển vị, sự hiện diện của vết nứt làm giãn rộng hàm mật độ xác suất, đồng nghĩa với việc gia tăng rủi ro vượt ngưỡng chuyển vị cho phép trong thiết kế.

Khi thể hiện dữ liệu qua các biểu đồ phân tán và bảng so sánh chuyển vị nút, các giá trị kỳ vọng tính từ chương trình TK.mw hoàn toàn trùng khớp với kết quả từ phần mềm SAP2000 (độ lệch bằng $0{,}00%$ ở hầu hết các chuyển vị thẳng). Điều này khẳng định ma trận độ cứng động lực kết hợp khai triển Neumann bậc nhất hoàn toàn đủ độ tin cậy để thay thế các phương pháp mô phỏng lặp tốn kém tài nguyên tính toán như Monte-Carlo trong các bài toán cơ học kết cấu ngẫu nhiên.

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

Dựa trên các kết quả tính toán giải tích và kiểm chứng thực nghiệm, nghiên cứu đưa ra 4 khuyến nghị ứng dụng thực tế:

  • Hoàn thiện quy trình kiểm định kết cấu có khuyết tật: Yêu cầu các đơn vị kiểm định chất lượng công trình áp dụng mô hình liên kết mềm lò xo tương đương để đánh giá độ cứng còn lại của dầm nứt, đặt mục tiêu giảm thiểu ít nhất $25%$ rủi ro đánh giá sai khả năng chịu lực trong vòng $6$ tháng kể từ khi áp dụng quy trình mới.
  • Tích hợp giải thuật khai triển Neumann vào phần mềm thiết kế: Đề xuất các viện nghiên cứu và doanh nghiệp phần mềm kết cấu bổ sung module tính toán tham số ngẫu nhiên vào các công cụ phân tích động lực học, hướng tới nâng cao độ chính xác phân tích độ tin cậy lên $15\text{ - }20%$, hoàn thành trong lộ trình $12$ tháng.
  • Thiết lập hệ thống quan trắc dao động thời gian thực: Ban quản lý các công trình cầu và nhà cao tầng cần lắp đặt cụm cảm biến gia tốc và cảm biến đo lực tại các vị trí xung yếu để phát hiện sớm các vết nứt bề mặt có độ sâu từ $2{,}5 \text{ cm}$ trở lên thông qua sự dịch chuyển dải tần số $10\text{ - }14 \text{ Hz}$, triển khai trong vòng $18$ tháng.
  • Cập nhật các tiêu chuẩn thiết kế kháng chấn: Kiến nghị Bộ Xây dựng và các cơ quan quản lý chuyên ngành bổ sung các chỉ dẫn kỹ thuật về hệ số an toàn ngẫu nhiên cho vật liệu và tải trọng động vào hệ thống tiêu chuẩn quốc gia, lộ trình thực hiện trong $24$ tháng.

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

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng chính:

  • Kỹ sư kết cấu và chuyên gia thiết kế công trình: Giúp nắm vững phương pháp tính toán chuyển vị và nội lực chính xác khi công trình làm việc trong điều kiện vật liệu biến tính và chịu tải trọng ngẫu nhiên.
  • Cán bộ kiểm định và giám định chất lượng: Cung cấp cơ sở khoa học để đánh giá mức độ an toàn của cấu kiện xuất hiện vết nứt cục bộ thông qua các thông số đo đạc dao động thực tế.
  • Giảng viên và nhà nghiên cứu cơ học công trình: Sử dụng làm tài liệu giảng dạy chuyên đề động lực học công trình ngẫu nhiên, cơ học phá hủy và phương pháp phần tử hữu hạn nâng cao.
  • Học viên cao học và nghiên cứu sinh ngành Cơ kỹ thuật, Xây dựng: Cung cấp khung phương pháp luận kết hợp giữa tính toán giải tích trên Maple, mô phỏng số và kỹ thuật thực nghiệm đo dao động bằng cảm biến.

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

Mô hình lò xo đàn hồi tương đương đại diện cho vết nứt hoạt động như thế nào? Mô hình này chuyển đổi sự gián đoạn hình học và tập trung ứng suất tại vết nứt thành các lò xo đàn hồi cục bộ có độ cứng xoay và độ cứng trượt hữu hạn. Chẳng hạn, khi chiều sâu vết nứt là $5{,}37 \text{ cm}$, độ cứng lò xo xoay tương đương được xác định bằng $123.456 \text{ kNm}$, phản ánh độ mềm tăng lên tại mặt cắt nứt.

Tại sao nên sử dụng phương pháp khai triển Neumann thay vì mô phỏng Monte-Carlo? Khai triển Neumann cho phép biểu diễn ma trận nghịch đảo ngẫu nhiên thành chuỗi chu kỳ giải tích, giúp tính trực tiếp kỳ vọng và phương sai chuyển vị mà chỉ mất vài giây xử lý. Phương pháp này giảm hơn $90%$ thời gian tính toán so với việc chạy hàng chục nghìn vòng lặp ngẫu nhiên của mô phỏng Monte-Carlo.

Yếu tố nào ảnh hưởng mạnh nhất đến phương sai chuyển vị của dầm có vết nứt? Phân tích độ nhạy khẳng định hệ số biến thiên của mô đun đàn hồi ($\epsilon_1$) tác động vượt trội đến phương sai chuyển vị so với hệ số biến thiên tải trọng ($\epsilon_2$). Khi $\epsilon_1$ tăng từ $0{,}01$ lên $0{,}30$, biên độ phương sai chuyển vị biến thiên mạnh mẽ và bao quát toàn bộ dải giá trị tới hạn của kết cấu.

Điều kiện biên liên kết ảnh hưởng thế nào đến độ võng của kết cấu có vết nứt? Liên kết càng cứng thì kỳ vọng chuyển vị và phương sai ngẫu nhiên càng giảm mạnh. Chuyển vị dầm conson ngàm một đầu đạt $-0{,}006557 \text{ m}$ (phương sai $3{,}8695 \times 10^{-8} \text{ m}^2$), trong khi dầm hai đầu ngàm chỉ võng $-0{,}0008561 \text{ m}$ (phương sai $4{,}9795 \times 10^{-10} \text{ m}^2$), giúp giảm thiểu tối đa độ bất định dao động.

Các tần số thí nghiệm 10 Hz, 13 Hz và 14 Hz có ý nghĩa gì trong việc kiểm chứng lý thuyết? Các tần số này quét qua các chế độ kích thích dao động cưỡng bức gần miền cộng hưởng của dầm thí nghiệm. Kết quả đo gia tốc và lực kích thích tại $10\text{ - }14 \text{ Hz}$ đã xác nhận sự suy giảm tần số riêng khi dầm có vết nứt, đồng thời khẳng định độ chính xác của mô hình lý thuyết với sai số dưới $5{,}8%$.

Kết luận

  • Đã thiết lập thành công mô hình ma trận độ cứng động lực ngẫu nhiên cho phần tử thanh phẳng có vết nứt mở ngang chịu uốn và kéo nén đồng thời.
  • Ứng dụng hiệu quả phương pháp khai triển Neumann để xác định tường minh biểu thức kỳ vọng và phương sai của chuyển vị và nội lực kết cấu.
  • Giải quyết hoàn chỉnh $10$ bài toán ứng dụng trên dầm và hệ thanh với các điều kiện biên khác nhau, làm sáng tỏ mức độ nhạy cảm của đáp ứng kết cấu đối với các biến ngẫu nhiên vật liệu.
  • Xây dựng hệ thống thí nghiệm dao động dầm hoàn chỉnh với cảm biến gia tốc và Loadcell, cung cấp bộ dữ liệu thực nghiệm chuẩn mực xác nhận độ tin cậy của thuật toán tính toán.
  • Mở ra hướng nghiên cứu tiếp theo về mô hình vết nứt thở phi tuyến và mở rộng cho kết cấu không gian chịu tải trọng động đất và sóng biển ngẫu nhiên.

Quý độc giả, kỹ sư kết cấu và các nhà nghiên cứu quan tâm có thể ứng dụng các giải thuật và bảng số liệu của luận văn để nâng cao chất lượng tính toán và kiểm định an toàn công trình.