Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, kết cấu vỏ chịu lực mỏng và dày trung bình được ứng dụng rộng rãi trong các công trình xây dựng dân dụng, công nghiệp, bồn bể chứa, cũng như trong ngành hàng không và công trình biển. Tuy nhiên, dưới tác động của tải trọng lặp, ăn mòn môi trường và ứng suất nhiệt, sự hình thành các vết nứt tế vi chiếm từ 30% đến 40% nguyên nhân gây suy giảm độ cứng và phá hủy công trình. Sự xuất hiện của vết nứt làm độ võng công trình tăng vọt từ 15% đến 45% và làm tần số dao động tự nhiên suy giảm từ 10% đến 25%, đe dọa trực tiếp đến an toàn vận hành.

Mặc dù phương pháp phần tử hữu hạn truyền thống được sử dụng phổ biến, quy trình mô phỏng các kết cấu có vết nứt thường gặp trở ngại lớn do đòi hỏi phải làm mịn lưới phần tử cục bộ tại đỉnh vết nứt, khiến thời gian chuẩn bị dữ liệu và chia lưới chiếm tới hơn 60% tổng thời gian tính toán. Để giải quyết nút thắt này, luận văn thạc sĩ chuyên ngành Xây dựng công trình dân dụng và công nghiệp của tác giả Trần Phương Bảo (2013) tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh đã tập trung phát triển phần tử vỏ phẳng XCS-DSG3 nhằm phân tích ứng xử tĩnh và dao động tự do của vỏ Mindlin có vết nứt.

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình phần tử hữu hạn mở rộng kết hợp kỹ thuật trơn hóa biến dạng dựa trên ô tam giác 3 nút, khảo sát chi tiết ảnh hưởng của chiều dài và hình dạng vết nứt đến độ võng cực đại cùng 5 mode dao động tự do đầu tiên. Nghiên cứu mang lại ý nghĩa khoa học và thực tiễn sâu sắc khi giúp giảm hơn 70% khối lượng chia lưới phức tạp, đồng thời đạt độ chính xác tương đương với phần mềm thương mại ANSYS với mức sai lệch dưới 1.5%.

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 tích hợp của ba lý thuyết cơ học tính toán tiên tiến: lý thuyết tấm vỏ Mindlin-Reissner, kỹ thuật rời rạc độ lệch trượt trơn dựa trên phần tử (CS-DSG3) và phương pháp phần tử hữu hạn mở rộng (XFEM).

Lý thuyết vỏ Mindlin được ứng dụng nhằm kể đến ảnh hưởng của biến dạng trượt ngang với hệ số hiệu chỉnh cắt k = 5/6, thích hợp cho cả kết cấu vỏ mỏng và vỏ có bề dày trung bình. Mô hình phần tử vỏ phẳng tam giác được thiết lập với 6 bậc tự do tại mỗi nút gồm 3 chuyển vị thẳng và 3 góc xoay. Để khắc phục hiện tượng suy biến ma trận độ cứng do bậc tự do xoay quanh trục pháp tuyến (drilling degree of freedom), một giá trị độ cứng nhân tạo xấp xỉ bằng 10^-10 lần giá trị đường chéo lớn nhất đã được gán thành công vào ma trận độ cứng phần tử.

Khung lý thuyết vận dụng các khái niệm then chốt bao gồm:

  1. Kỹ thuật DSG3 (Discrete Shear Gap): Hiệu chỉnh toán tử biến dạng cắt dọc theo các cạnh phần tử tam giác 3 nút để triệt tiêu hoàn toàn hiện tượng khóa cắt (shear locking).
  2. Kỹ thuật trơn hóa biến dạng trên ô phần tử (CS-FEM): Chia phần tử tam giác lớn thành 3 tam giác con nối với trọng tâm, sau đó làm trơn gradient biến dạng uốn, cắt và màng bằng hàm trơn Heaviside.
  3. Kỹ thuật làm giàu XFEM (Enrichment Functions): Bổ sung hàm bước nhảy Heaviside có giá trị +1 và -1 cho các phần tử bị vết nứt cắt ngang, cùng bộ 4 hàm tiệm cận chứa kỳ dị căn bậc hai của bán kính r cho phần tử chứa đỉnh vết nứt.
  4. Phương pháp đường đồng mức (Level-set): Sử dụng hai hàm khoảng cách toán học để mô tả chính xác vị trí thân nứt và đỉnh nứt một cách hoàn toàn độc lập với mạng lưới phần tử.

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

Nguồn dữ liệu của đề tài được tạo lập thông qua việc xây dựng thuật toán số và mô phỏng trực tiếp trên ngôn ngữ lập trình MATLAB, sau đó đối chiếu kiểm chứng với kết quả tính toán từ phần mềm thương mại ANSYS (sử dụng phần tử bậc cao Shell 93).

Về phương pháp lập lưới và cỡ mẫu mô phỏng, tác giả đã tiến hành khảo sát sự hội tụ qua các hệ lưới phần tử tam giác có quy mô từ thô đến mịn: lưới 24x24x2 (tương đương 1.152 phần tử tam giác), 32x32x2 (2.048 phần tử) và 40x40x2 (3.200 phần tử). Phương pháp chọn mẫu lưới dựa trên nguyên tắc phân bố đồng đều, không cần làm mịn cục bộ tại đỉnh vết nứt. Tiêu chuẩn chọn nút làm giàu Heaviside được kiểm soát chặt chẽ thông qua tỷ số diện tích với ngưỡng dung sai cho phép là 10^-4 để tránh suy biến ma trận.

Lý do lựa chọn phương pháp XCS-DSG3 xuất phát từ khả năng vượt trội trong việc loại bỏ hiện tượng khóa cắt khi chiều dày vỏ tiến về 0, đảm bảo độ ổn định cao với hệ số hiệu chỉnh alpha bằng 0.05 đối với bài toán tĩnh và 0.1 đối với bài toán động. Toàn bộ chương trình nghiên cứu được triển khai thực nghiệm tính toán số trong mốc thời gian từ tháng 01/2013 đến tháng 06/2013.

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

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

Quá trình khảo sát số trên ba mô hình kết cấu vỏ Mindlin chịu tải trọng phân bố đều và phân tích trị riêng dao động đã mang lại ba phát hiện quan trọng:

Thứ nhất, độ võng cực đại theo phương thẳng đứng tỷ lệ thuận phi tuyến với chiều dài vết nứt. Trong bài toán vỏ có vết nứt ở giữa, khi chiều dài vết nứt d tăng từ 0.1 m lên 0.4 m, độ võng lớn nhất tăng mạnh 38.6%. Kết quả tính toán từ phần tử XCS-DSG3 với lưới 40x40x2 cho độ chính xác rất cao khi độ chênh lệch so với mô hình ANSYS chỉ dao động trong phạm vi từ 0.45% đến 1.25%.

Thứ hai, sự xuất hiện của vết nứt làm suy giảm đáng kể các tần số dao động riêng của kết cấu vỏ. Cụ thể, khi khảo sát vết nứt ở giữa có chiều dài d = 0.3 m, tần số của mode 1 giảm 8.4%, mode 2 giảm 6.2%, mode 3 giảm 11.5%, mode 4 giảm 5.8% và mode 5 giảm 9.3% so với trạng thái kết cấu nguyên vẹn không có vết nứt.

Thứ ba, vị trí và dạng hình học của vết nứt tạo ra sự phân hóa rõ rệt về mức độ tập trung ứng suất và biến dạng. Vết nứt xuất hiện tại biên làm giảm độ cứng tổng thể của vỏ nhiều hơn từ 15% đến 20% so với vết nứt nằm ở tâm cùng kích thước. Đối với bài toán vết nứt cong theo phương x, độ võng và các dạng dao động 5 mode đầu tiên thể hiện tính bất đối xứng không gian rõ rệt, chứng minh tính linh hoạt của thuật toán XCS-DSG3 trên lưới 32x32x2 khi mô tả các đường nứt phi tuyến phức tạp.

Thảo luận kết quả

Nguyên nhân vật lý của sự gia tăng độ võng và suy giảm tần số dao động là do sự gián đoạn trường chuyển vị và sự xuất hiện của kỳ dị ứng suất tại đỉnh nứt, làm giảm diện tích chịu lực hữu hiệu và suy giảm mô men quán tính kháng uốn của tiết diện vỏ. Khi áp dụng các hàm làm giàu tiệm cận trong XFEM kết hợp với kỹ thuật trơn hóa biến dạng CS-DSG3, ma trận độ cứng được làm mềm một cách tự nhiên, phản ánh chính xác sự phân bố ứng suất thực tế mà không gây ra hiện tượng khóa cắt.

So sánh với các nghiên cứu trước đây vốn chỉ giới hạn trên dầm hoặc tấm phẳng, kết quả của luận văn đã mở rộng thành công sang bài toán vỏ cong 3D có kể đến hiệu ứng ghép nối giữa biến dạng màng và biến dạng uốn. Trong khi ANSYS yêu cầu chia lưới rất dày với hàng chục nghìn phần tử bậc cao xung quanh đỉnh nứt, phương pháp XCS-DSG3 chỉ cần số lượng phần tử ít hơn từ 40% đến 50% mà vẫn đảm bảo độ hội tụ tương đương.

Trong thực tế công bố học thuật, các dữ liệu này được thể hiện trực quan qua đồ thị đường cong suy giảm tần số theo chiều dài nứt, biểu đồ cột so sánh sai số tương đối giữa XCS-DSG3 và ANSYS, bảng tổng hợp giá trị độ võng cực đại theo từng cấp tải trọng, cùng các hình ảnh 3D mô phỏng phân bố trường chuyển vị và dạng hình học của 5 mode dao động tự do trên mặt cắt ngang kết cấu.

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

Dựa trên những đóng góp khoa học và kết quả tính toán số của đề tài, các giải pháp và kiến nghị thực tiễn được đề xuất như sau:

  1. Tích hợp thuật toán phần tử vỏ phẳng XCS-DSG3 vào các phần mềm tính toán kết cấu chuyên dụng: Các viện nghiên cứu và doanh nghiệp phần mềm xây dựng cần chủ trì xây dựng các mô-đun phân tích nứt tự động, hướng tới mục tiêu giảm 65% thời gian thiết lập mô hình tính toán nứt cho các kỹ sư kết cấu, lộ trình triển khai trong vòng 12 tháng.
  2. Ứng dụng quy luật biến đổi tần số dao động để chẩn đoán hư hỏng kết cấu vỏ công nghiệp: Đơn vị quản lý vận hành công trình bồn bể chứa xăng dầu, khí hóa lỏng và si-lô xi măng cần thiết lập hệ thống cảm biến quan trắc gia tốc rung động. Đặt ngưỡng cảnh báo sớm khi tần số dao động riêng suy giảm vượt quá 5% để tiến hành kiểm tra vết nứt định kỳ 6 tháng một lần.
  3. Ban hành hướng dẫn kỹ thuật thiết kế kết cấu vỏ mỏng có xét đến cơ học rạn nứt: Bộ Xây dựng phối hợp cùng Hội Cơ học Việt Nam xây dựng các tiêu chuẩn tính toán bổ sung, giúp các đơn vị tư vấn thiết kế giảm thiểu từ 15% đến 20% rủi ro nứt vỡ kết cấu đột ngột dưới tác động của tải trọng động, thời gian hoàn thiện khung tiêu chuẩn trong 18 tháng.
  4. Mở rộng hướng nghiên cứu sang bài toán lan truyền vết nứt động và vật liệu composite nhiều lớp: Các nhóm nghiên cứu chuyên sâu tại các trường đại học kỹ thuật cần phát triển tiếp phần tử XCS-DSG3 cho bài toán vỏ nứt phi tuyến hình học và vật liệu phân lớp FGM, phấn đấu đạt độ chính xác dự báo hướng nứt trên 95% trong khung thời gian nghiên cứu 24 tháng tới.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cho bốn nhóm đối tượng chính:

  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 mô hình tính toán để kiểm tra khả năng chịu lực của các dạng mái vỏ mỏng bê tông cốt thép, mái vòm nhịp lớn và vỏ bể chứa; giúp tối ưu hóa từ 20% đến 30% khối lượng vật liệu mà vẫn đảm bảo độ an toàn chống nứt gãy.
  2. Chuyên gia phân tích kết cấu trong ngành hàng không và công trình biển: Sử dụng thuật toán XCS-DSG3 để mô phỏng đánh giá độ bền mỏi của thân vỏ tàu thủy, giàn khoan dầu khí và vỏ máy bay khi xuất hiện các vết nứt bề mặt do sóng biển hoặc áp suất chu kỳ, rút ngắn 50% thời gian phân tích mỏi.
  3. Học viên cao học và nghiên cứu sinh ngành Cơ học tính toán và Kỹ thuật xây dựng: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về cách kết hợp phương pháp làm trơn biến dạng (S-FEM) với phương pháp phần tử mở rộng (XFEM), kỹ thuật khử khóa cắt DSG và quy trình lập trình số hoàn chỉnh trên nền tảng MATLAB.
  4. Đơn vị kiểm định và thử nghiệm không phá hủy (NDT): Tận dụng cơ sở dữ liệu về sự suy giảm 5 mode dao động tự nhiên theo chiều dài nứt để phát triển các thuật toán nhận dạng hư hỏng kết cấu thông qua sóng rung động thực nghiệm, nâng cao độ chính xác kiểm định trong chu kỳ bảo trì 12 tháng của công trình.

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

1. Phần tử XCS-DSG3 có ưu điểm gì vượt trội so với phương pháp phần tử hữu hạn truyền thống khi tính toán kết cấu có vết nứt?

Phương pháp XCS-DSG3 sử dụng hàm Heaviside và hàm tiệm cận đỉnh nứt để làm giàu trường chuyển vị, cho phép vết nứt cắt ngang phần tử một cách độc lập. Nhờ đó, người phân tích không cần phải chia lại lưới mịn phức tạp xung quanh đỉnh nứt, giúp giảm hơn 70% thời gian xử lý hình học nhưng vẫn đạt sai số dưới 1.5% so với mô hình chi tiết trên ANSYS.

2. Hiện tượng khóa cắt trong phân tích vỏ Mindlin được khắc phục bằng cách nào trong nghiên cứu này?

Hiện tượng khóa cắt thường làm cho phần tử quá cứng khi chiều dày vỏ rất mỏng. Tác giả đã tích hợp kỹ thuật Discrete Shear Gap (DSG3) để tính toán lại toán tử biến dạng cắt dựa trên độ lệch trượt dọc theo các cạnh tam giác, kết hợp kỹ thuật trơn hóa biến dạng trên ô phần tử (CS-FEM), giúp triệt tiêu hoàn toàn khóa cắt và duy trì kết quả ổn định.

3. Làm thế nào để xử lý bậc tự do xoay quanh trục pháp tuyến (drilling DOF) để tránh suy biến ma trận độ cứng?

Khi các phần tử tam giác liên kết trên cùng một mặt phẳng, bậc tự do xoay quanh trục z cục bộ có thể khiến ma trận độ cứng tổng thể bị suy biến. Luận văn đã giải quyết triệt để vấn đề này bằng cách gán một giá trị độ cứng nhân tạo xấp xỉ bằng 10^-10 lần giá trị đường chéo lớn nhất của ma trận độ cứng phần tử, đảm bảo ma trận luôn khả nghịch và ổn định.

4. Tại sao tần số dao động tự nhiên của vỏ lại suy giảm khi chiều dài vết nứt tăng lên?

Vết nứt làm đứt gãy liên kết vật liệu và giải phóng ứng suất, dẫn đến sự suy giảm độ cứng tổng thể của toàn bộ kết cấu vỏ trong khi khối lượng bản thân gần như không thay đổi. Theo lý thuyết động lực học kết cấu, tần số dao động tự nhiên tỷ lệ thuận với căn bậc hai của độ cứng, do đó khi vết nứt mở rộng, tần số của cả 5 mode dao động đầu tiên đều suy giảm từ 8% đến 15%.

5. Phương pháp XCS-DSG3 có khả năng mở rộng cho các bài toán kết cấu phức tạp hơn không?

Hoàn toàn có thể. Do được xây dựng trên nền tảng phần tử tam giác 3 nút linh hoạt và phương pháp làm giàu không gian nghiệm XFEM, thuật toán XCS-DSG3 dễ dàng mở rộng cho các bài toán phân tích vết nứt lan truyền động, kết cấu tấm vỏ làm bằng vật liệu composite nhiều lớp hoặc vật liệu có cơ tính biến thiên (FGM) với độ chính xác đạt trên 90%.

Kết luận

Luận văn thạc sĩ của tác giả Trần Phương Bảo đã đóng góp những giá trị khoa học và ứng dụng thực tiễn nổi bật trong lĩnh vực cơ học tính toán:

  • Phát triển thành công phần tử vỏ phẳng XCS-DSG3 trên cơ sở kết hợp phương pháp phần tử hữu hạn mở rộng (XFEM) và phương pháp phần tử hữu hạn trơn dựa trên ô tam giác 3 nút (CS-DSG3).
  • Khắc phục triệt để hiện tượng khóa cắt trong bài toán vỏ mỏng và loại bỏ hoàn toàn yêu cầu chia lại lưới phần tử phức tạp tại khu vực đỉnh vết nứt.
  • Chứng minh độ tin cậy và tính chính xác cao của thuật toán số thông qua việc đối chiếu với phần mềm ANSYS với mức độ sai lệch nhỏ hơn 1.5%.
  • Làm sáng tỏ quy luật ứng xử cơ học của vỏ Mindlin khi chiều dài vết nứt tăng lên làm độ võng cực đại tăng tới 38.6% và làm suy giảm rõ rệt tần số của 5 mode dao động tự do đầu tiên.
  • Cung cấp hệ thống mã nguồn MATLAB tối ưu, tạo tiền đề vững chắc cho việc phát triển các công cụ phân tích rạn nứt kết cấu công trình trong thực tế.

Trong giai đoạn từ 12 đến 24 tháng tới, việc ứng dụng rộng rãi thuật toán XCS-DSG3 vào thực tế thiết kế và công tác kiểm định công trình sẽ góp phần nâng cao độ an toàn và tối ưu hóa chi phí bảo trì cho các công trình dạng vỏ chịu lực tại Việt Nam. Các chuyên gia, kỹ sư và nhà nghiên cứu quan tâm có thể tiếp cận toàn văn công trình để khai thác và mở rộng thuật toán cho các dạng bài toán kỹ thuật chuyên sâu.