Tổng quan nghiên cứu

Trong kỹ thuật xây dựng và kết cấu công trình hiện đại, sự suy giảm chất lượng vật liệu sau khoảng 10 đến 20 năm vận hành thường dẫn đến sự xuất hiện của các vết nứt tế vi. Các vết nứt này làm suy giảm nghiêm trọng độ cứng tổng thể, đe dọa độ an toàn kết cấu và rút ngắn tuổi thọ công trình khi chịu tác động của tải trọng động lặp. Việc mô phỏng chính xác ứng xử động học của kết cấu tấm có khuyết tật nứt đóng vai trò sống còn trong việc kiểm định an toàn và cảnh báo sớm thảm họa công trình.

Mục tiêu cốt lõi của nghiên cứu là phát triển mô hình phần tử hữu hạn mở rộng sử dụng phần tử tứ giác bốn nút nhằm phân tích dao động riêng và tần số dao động tự do của tấm Mindlin có chứa vết nứt. Nghiên cứu tập trung giải quyết bài toán nứt hai chiều với các dạng vết nứt biên và vết nứt trung tâm, có chiều dài vết nứt biến thiên từ 0 đến 1.0 mét trên tấm vuông kích thước tiêu chuẩn 1.0 x 1.0 mét và chiều dày 0.001 mét.

Nghiên cứu được triển khai thực hiện tại Trường Đại học Bách khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh trong giai đoạn từ tháng 7 năm 2012 đến tháng 6 năm 2013. Luận văn đã đóng góp giải pháp số hóa toàn diện thông qua thuật toán tự lập trình trên nền tảng ngôn ngữ kỹ thuật, cho phép tính toán chính xác 5 dạng dao động tự do đầu tiên của kết cấu. Kết quả tính toán đạt độ tin cậy cao với mức sai số tuyệt đối được kiểm chứng luôn nằm dưới ngưỡng 5 phần trăm so với các công bố quốc tế uy tín, tạo bước đột phá trong việc tối ưu hóa tài nguyên tính toán mà không cần chia lại lưới phần tử.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Khung lý thuyết của công trình dựa trên nền tảng cơ học phá hủy đàn hồi tuyến tính kết hợp lý thuyết tấm biến dạng cắt bậc nhất Mindlin-Reissner. Tấm được xét đến đầy đủ các thành phần biến dạng uốn, biến dạng cắt và biến dạng trượt hình học với hệ số điều chỉnh cắt được chuẩn hóa ở mức 5/6. Mỗi nút của phần tử tấm có 3 bậc tự do độc lập gồm độ võng và hai góc xoay quanh các trục tọa độ.

Để giải quyết tính gián đoạn của trường chuyển vị khi có vết nứt xuất hiện, phương pháp phần tử hữu hạn mở rộng dựa trên nguyên lý phân chia đơn vị được áp dụng. Phương pháp bổ sung hai nhóm hàm làm giàu vào không gian xấp xỉ tiêu chuẩn:

Thứ nhất, hàm bước nhảy Heaviside nhận giá trị từ dương 1 đến âm 1 nhằm mô tả sự gián đoạn của bề mặt vết nứt cắt ngang phần tử. Một dung sai diện tích hỗ trợ bằng 0.0001 được thiết lập để loại trừ các nút có diện tích cắt quá nhỏ, ngăn chặn hiện tượng suy biến ma trận độ cứng hệ thống.

Thứ hai, bốn hàm mở rộng tiệm cận dạng lượng giác phụ thuộc vào hệ tọa độ cực cục bộ tại đầu vết nứt được tích hợp để mô tả trường ứng suất kỳ dị.

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

Nguồn dữ liệu của nghiên cứu được tạo lập từ các bài toán mô phỏng số học với kích thước tấm chuẩn hóa cạnh 1.0 mét, mô-đun đàn hồi Young đạt 210 GPa, hệ số Poisson bằng 0.3 và khối lượng riêng 800 kg/m3. Cỡ mẫu nghiên cứu được thiết lập qua việc khảo sát 4 cấu hình lưới phần tử hữu hạn tăng dần mật độ từ 10x10, 16x16, 20x20 đến 24x24 phần tử tứ giác 4 nút (tương ứng với 576 phần tử và 625 nút trên toàn miền bài toán).

Phương pháp chọn lưới có cấu trúc đồng nhất được áp dụng độc lập với đường nứt, cho phép vết nứt đi qua bất kỳ vị trí nào trong phần tử mà không đòi hỏi tái lập lưới cục bộ. Lý do lựa chọn phương pháp phân tích số này là nhằm tận dụng khả năng tính toán ma trận độ cứng và ma trận khối lượng giải tích kết hợp tích phân số Gauss hai chiều. Trong đó, tích phân thành phần uốn sử dụng sơ đồ 2x2 điểm Gauss và thành phần trượt cắt sử dụng sơ đồ thu gọn 1x1 điểm Gauss nhằm triệt tiêu triệt để hiện tượng khóa trượt. Thuật toán giải bài toán trị riêng tuyến tính được thực hiện hoàn chỉnh trong chu trình 12 tháng nghiên cứu.

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

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

Quá trình phân tích số trên mô hình lưới 24x24 phần tử cho thấy sự thay đổi rõ rệt của tần số dao động riêng theo 5 mode đầu tiên khi tỷ lệ chiều dài vết nứt c/a tăng từ 0 đến 1.0:

Thứ nhất, đối với vết nứt ở cạnh tấm vuông tựa đơn, chiều dài vết nứt tăng làm suy giảm tần số ở tất cả các mode dao động nhưng với mức độ phân hóa rất lớn. Tần số mode 1 giảm 18.39% từ mức ban đầu khi chưa nứt xuống khi nứt hoàn toàn, trong khi tần số mode 3 chỉ suy giảm nhẹ 5.3%.

Thứ hai, các mode dao động bậc chẵn thể hiện độ nhạy cảm vượt trội đối với khuyết tật cạnh. Tần số dao động của mode 2 sụt giảm mạnh nhất lên tới 67.5%, mode 4 giảm 40.86% và mode 5 giảm 24.5%. Điều này chứng minh vết nứt cạnh cắt đứt trực tiếp các đường truyền ứng suất uốn chính của các dạng dao động xoắn và uốn bậc cao.

Thứ ba, sự kiểm chứng độ chính xác khi so sánh với phương pháp của Bachene cho thấy sai số lớn nhất chỉ là 4.10% ở mode 5 và sai số nhỏ nhất đạt 0.84% ở mode 2 khi tỷ lệ vết nứt c/a bằng 0.5.

Thứ tư, đối với bài toán vết nứt ở giữa tấm, kết quả so sánh với dữ liệu của Liew và Stahl ghi nhận mức sai số cực nhỏ từ 0.83% đến 4.98%, khẳng định tính ổn định hội tụ cao của phần tử tứ giác mở rộng.

Thảo luận kết quả

Nguyên nhân vật lý sâu xa của hiện tượng sụt giảm tần số dao động là do sự xuất hiện và mở rộng của vết nứt làm đứt gãy tính liên tục của trường chuyển vị, dẫn đến việc suy giảm trực tiếp ma trận độ cứng tổng thể của kết cấu tấm. Trong khi đó, khối lượng bản thân của tấm gần như không đổi, làm cho tỷ số độ cứng trên khối lượng giảm xuống, kéo theo tần số góc riêng giảm theo hàm căn bậc hai.

Dữ liệu nghiên cứu được biểu diễn trực quan và sinh động thông qua các bảng đối chiếu tham số sai số tương quan và hệ thống đồ thị đường cong phân bố tần số. Bảng tổng hợp số liệu thể hiện sự hội tụ của 5 mode dao động từ lưới 10x10 đến 24x24 giúp quan sát rõ tốc độ ổn định nghiệm của phương pháp mở rộng. Đồ thị biểu diễn mối quan hệ giữa chiều dài vết nứt từ 0 đến 1.0 mét và tần số dao động phản ánh đường cong suy giảm phi tuyến tính rõ nét, giúp các kỹ sư dễ dàng nhận biết ngưỡng nguy hiểm của cấu kiện. Kết quả này hoàn toàn tương thích và tiệm cận với các công bố của những nhà khoa học hàng đầu thế giới về cơ học tính toán.

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

Dựa trên kết quả nghiên cứu toàn diện, tác giả đưa ra 4 khuyến nghị then chốt nhằm thúc đẩy ứng dụng thực tiễn:

Thứ nhất, tích hợp ngay thuật toán phần tử tứ giác mở rộng vào các phần mềm tính toán kết cấu thương mại trong vòng 12 tháng tới. Các viện nghiên cứu và đơn vị phát triển phần mềm cần triển khai giải pháp này nhằm giảm ít nhất 40% thời gian thiết lập mô hình tính toán so với phương pháp chia lưới truyền thống.

Thứ hai, thiết lập hệ thống quan trắc sức khỏe công trình dựa trên độ nhạy tần số dao động cho các kết cấu tấm sàn nhịp lớn và vỏ công nghiệp. Các kỹ sư bảo trì cần thực hiện đo đạc tần số dao động định kỳ 6 tháng một lần, lấy ngưỡng sụt giảm tần số mode 2 vượt quá 30% làm chỉ số cảnh báo nguy cấp về sự tồn tại của vết nứt biên có chiều dài trên 0.4 mét.

Thứ ba, mở rộng phạm vi nghiên cứu sang bài toán vết nứt cong phức tạp và vết nứt không gian ba chiều trong thời gian 24 tháng. Nhóm nghiên cứu cần kết hợp phương pháp điều chỉnh mức độ cao để tự động hóa hoàn toàn việc cập nhật quỹ đạo lan truyền nứt động dưới tác động của tải trọng mỏi lặp chu kỳ lớn.

Thứ tư, xây dựng bộ tiêu chuẩn hướng dẫn đánh giá độ an toàn kết cấu có khuyết tật cục bộ trong vòng 18 tháng. Cơ quan quản lý ngành xây dựng cần quy định hệ số an toàn tăng thêm từ 1.15 đến 1.30 đối với các cấu kiện tấm mỏng làm việc trong môi trường rung động liên tục khi phát hiện vết nứt cục bộ.

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

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng cụ thể:

Nhóm thứ nhất là các kỹ sư kết cấu công trình dân dụng và công nghiệp. Luận văn cung cấp phương pháp luận chuẩn xác để đánh giá khả năng chịu lực còn lại và mức độ an toàn động học của các tấm sàn, vách cứng, bản mặt cầu bê tông hoặc thép đang xuất hiện vết nứt từ 0.1 đến 1.0 mét trong quá trình khai thác.

Nhóm thứ hai là các nhà khoa học và nghiên cứu sinh chuyên ngành Cơ học tính toán và Cơ học phá hủy. Công trình cung cấp hệ phương trình toán học chi tiết, các công thức ma trận độ cứng, ma trận khối lượng và kỹ thuật làm giàu nút, làm nền tảng vững chắc để phát triển các phần tử tấm vỏ bậc cao.

Nhóm thứ ba là giảng viên và học viên cao học tại các trường đại học kỹ thuật. Đây là tài liệu tham khảo chất lượng cao cho các môn học chuyên đề như Phương pháp phần tử hữu hạn nâng cao, Động lực học công trình và Cơ học rạn nứt với 6 ví dụ số minh họa rõ ràng.

Nhóm thứ tư là các chuyên gia giám định chất lượng và kiểm định công trình. Tài liệu giúp xây dựng quy trình chẩn đoán không phá hủy dựa trên phổ tần số dao động riêng thu thập từ hiện trường.

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

Phương pháp phần tử hữu hạn mở rộng mang lại ưu điểm vượt trội gì so với phương pháp phần tử hữu hạn truyền thống khi giải bài toán nứt?

Phương pháp mở rộng cho phép biểu diễn sự gián đoạn của vết nứt hoàn toàn độc lập với lưới phần tử thông qua các hàm làm giàu Heaviside và hàm tiệm cận đầu vết nứt. Nhờ đó, người tính toán không cần phải chia lại lưới mịn phức tạp quanh đầu vết nứt, tiết kiệm hơn 50% khối lượng lao động tính toán và giải quyết triệt để vấn đề biến dạng phần tử.

Tại sao tần số dao động của mode thứ 2 lại có mức sụt giảm nghiêm trọng lên tới 67.5% khi vết nứt cạnh phát triển?

Dạng dao động thứ 2 có đường nút dao động và trường ứng suất uốn cực đại nằm trùng hoặc rất gần với vị trí vết nứt biên ở giữa cạnh tấm. Khi vết nứt xuất hiện với chiều dài tăng dần đến 1.0 mét, liên kết truyền mô-men tại vị trí này bị triệt tiêu hoàn toàn, làm suy giảm gần như toàn bộ độ cứng kháng uốn cục bộ của mode dao động này.

Kỹ thuật tích phân Gauss hai chiều được thiết lập ra sao để tránh hiện tượng khóa trượt đối với tấm mỏng?

Nghiên cứu áp dụng kỹ thuật tích phân Gauss có chọn lọc với sơ đồ 2x2 điểm Gauss cho ma trận độ cứng uốn và sơ đồ tích phân thu gọn 1x1 điểm Gauss cho ma trận độ cứng cắt. Việc thu gọn điểm tích phân đối với thành phần cắt triệt tiêu triệt để hiện tượng khóa trượt giả tạo khi tỷ lệ chiều dày trên nhịp tấm đạt mức 1/1000.

Luận văn đã xử lý nguy cơ suy biến ma trận độ cứng khi vết nứt nằm quá sát nút phần tử như thế nào?

Tác giả áp dụng thuật toán kiểm soát diện tích hỗ trợ nút với tỷ lệ diện tích bị cắt phân định trên và dưới vết nứt. Nếu tỷ lệ diện tích nhỏ hơn ngưỡng dung sai 0.0001, nút đó sẽ tự động không được gán hàm mở rộng Heaviside, từ đó loại bỏ hoàn toàn tính phụ thuộc tuyến tính và giữ cho ma trận độ cứng luôn dương vô cùng định.

Tính chính xác của thuật toán lập trình trên ngôn ngữ kỹ thuật được kiểm chứng qua những cơ sở nào?

Thuật toán tự phát triển được kiểm chứng chéo với các công bố quốc tế kinh điển của Bachene, Liew và Stahl qua 5 mode dao động tự do. Mức sai số tương đối luôn duy trì ổn định trong khoảng từ 0.83% đến 4.98% trên các cấu hình lưới từ 10x10 đến 24x24 phần tử, khẳng định mã nguồn có độ chuẩn xác và độ tin cậy tuyệt đối.

Kết luận

Luận văn thạc sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp khoa học nổi bật:

  • Xây dựng thành công thuật toán phần tử tứ giác 4 nút mở rộng cho bài toán động lực học tấm Mindlin có xét đầy đủ biến dạng cắt và uốn.
  • Thiết lập quy trình làm giàu hàm Heaviside và hàm tiệm cận đầu vết nứt chính xác, loại trừ hiện tượng suy biến ma trận độ cứng qua ngưỡng dung sai diện tích 0.0001.
  • Chứng minh tính hội tụ và độ chính xác vượt trội của phương pháp với sai số luôn dưới 5 phần trăm so với các lời giải chuẩn mực quốc tế.
  • Định lượng chi tiết quy luật suy giảm tần số dao động của 5 mode đầu tiên khi chiều dài vết nứt thay đổi từ 0 đến 1.0 mét, phát hiện mức giảm sâu nhất 67.5% ở mode 2.
  • Cung cấp bộ mã nguồn lập trình hoàn chỉnh, mở ra công cụ phân tích động học tin cậy cho ngành xây dựng công trình.

Trong giai đoạn 6 đến 12 tháng tiếp theo, nhóm nghiên cứu dự kiến sẽ phát triển mô hình cho bài toán tấm composite nhiều lớp và phân tích lan truyền vết nứt động theo thời gian. Các đơn vị tư vấn thiết kế và kiểm định kết cấu hãy chủ động liên hệ và ứng dụng ngay giải pháp tính toán này vào các dự án thực tế để nâng cao chất lượng công trình và tối ưu hóa chi phí bảo trì.