Tổng quan nghiên cứu

Trong kỹ thuật kết cấu công trình, hàng không vũ trụ và công nghiệp đóng tàu, kết cấu vỏ chịu lực chiếm hơn 70% các cấu kiện chịu tải trọng phức tạp. Tuy nhiên, sự xuất hiện của các vết nứt tế vi hoặc vết nứt mỏi làm suy giảm từ 20% đến 45% khả năng chịu lực tổng thể, thay đổi đáng kể tần số dao động tự nhiên và gia tăng độ võng cục bộ. Việc mô phỏng chính xác ứng xử cơ học của vỏ có vết nứt là một thách thức lớn trong cơ học tính toán do hiện tượng tập trung ứng suất cao tại đỉnh nứt và hiện tượng khóa cắt (shear locking) khi bề dày vỏ giảm dần.

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 tập trung giải quyết bài toán phân tích ứng xử tĩnh và dao động tự do của vỏ Mindlin có vết nứt thông qua việc phát triển phần tử vỏ phẳng mở rộng XCS-DSG3. Mục tiêu cụ thể của nghiên cứu là xây dựng công thức phần tử tam giác 3 nút kết hợp giữa kỹ thuật trơn hóa biến dạng trượt (CS-DSG3) và phương pháp phần tử hữu hạn mở rộng (XFEM), từ đó khảo sát định lượng ảnh hưởng của chiều dài vết nứt từ 0,1 m đến 0,4 m đối với độ võng lớn nhất và 5 dạng dao động tự nhiên đầu tiên.

Nghiên cứu được triển khai tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong thời gian 5 tháng (từ tháng 01/2013 đến tháng 06/2013). Giá trị học thuật và thực tiễn của công trình thể hiện qua việc giảm hơn 60% thời gian tái tạo lưới phần tử so với phương pháp phần tử hữu hạn truyền thống, đồng thời duy trì sai số tương đối dưới 1,5% khi đối sánh với phần mềm thương mại ANSYS.

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 lý thuyết tấm và vỏ Mindlin-Reissner, cho phép 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 bằng 5/6 và mô đun đàn hồi trượt G phụ thuộc vào mô đun Young E cùng hệ số Poisson qua công thức G = E / [2(1 + ν)]. Phần tử vỏ phẳng được hình thành thông qua sự kết hợp giữa phần tử tấm Mindlin chịu uốn (3 bậc tự do w, βx, βy tại mỗi nút) và phần tử màng ứng suất phẳng (3 bậc tự do u, v, βz tại mỗi nút), tạo nên phần tử vỏ hoàn chỉnh với 6 bậc tự do trên mỗi nút sau khi chuyển đổi từ hệ tọa độ địa phương sang hệ tọa độ tổng thể.

Để loại bỏ hiện tượng khóa cắt, nghiên cứu tích hợp phương pháp rời rạc độ lệch trượt trơn dựa trên phần tử (Cell-based smoothed Discrete Shear Gap - CS-DSG3). Mỗi phần tử tam giác lớn được chia thành 3 tam giác con nối từ trọng tâm đến 3 đỉnh, sau đó áp dụng hàm trơn Heaviside để tính toán gradient biến dạng uốn, cắt và màng.

Đối với bài toán nứt, phương pháp phần tử hữu hạn mở rộng (XFEM) được áp dụng dựa trên nguyên lý phân hoạch đơn vị (Partition of Unity). Trường chuyển vị xấp xỉ được làm giàu bởi 2 nhóm hàm:

  1. Hàm gián đoạn Heaviside H(x) nhận giá trị +1 hoặc -1 để mô tả bước nhảy chuyển vị qua bề mặt vết nứt.
  2. Bộ 4 hàm tiệm cận B(r, θ) chứa thành phần căn bậc hai của bán kính r kết hợp hàm lượng giác góc θ để nắm bắt kỳ dị ứng suất bậc 1/sqrt(r) tại đỉnh vết nứt.

Vị trí vết nứt được định vị độc lập với lưới nhờ phương pháp đường đồng mức (Level-set) qua hai hàm khoảng cách φ(x) và ψ(x). Nhằm đảm bảo tính ổn định số, nghiên cứu thiết lập hệ số ổn định dao động cắt α bằng 0,05 cho bài toán tĩnh, α bằng 0,1 cho bài toán động, và bổ sung giá trị độ cứng nhân tạo xấp xỉ 10^-10 lần giá trị đường chéo lớn nhất cho bậc tự do xoay quanh trục pháp tuyến.

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

Nghiên cứu sử dụng phương pháp mô phỏng số và phân tích phần tử hữu hạn nâng cao trên môi trường lập trình MATLAB, đối chứng trực tiếp với mô hình phần tử vỏ bậc cao Shell93 trong phần mềm ANSYS.

Dữ liệu tính toán được khảo sát trên 3 mô hình kết cấu vỏ Mindlin tiêu chuẩn:

  • Vỏ có vết nứt thẳng ở tâm dọc theo trục y.
  • Vỏ có vết nứt thẳng ở biên dọc theo trục y.
  • Vỏ có vết nứt cong ở tâm theo phương trục x.

Quy trình chọn mẫu phân tích số sử dụng hệ thống lưới tam giác có cấu trúc với các mật độ lưới khảo sát lần lượt là 24x24x2 phần tử (1.152 phần tử), 32x32x2 phần tử (2.048 phần tử) và 40x40x2 phần tử (3.200 phần tử tương ứng 1.681 nút). Lý do lựa chọn phương pháp phân tích số rời rạc không phụ thuộc vết nứt là vì phương pháp này loại bỏ hoàn toàn yêu cầu tái tạo lưới phức tạp tại đỉnh nứt, giải quyết triệt để sự suy biến ma trận độ cứng nhờ tiêu chuẩn diện tích hỗ trợ với ngưỡng dung sai cắt r_ab, r_be lớn hơn hoặc bằng 10^-4. Tích phân số trên các tam giác con được thực hiện thông qua sơ đồ cầu phương Gauss với số lượng từ 3 đến 5 điểm Gauss, bảo đảm nghiệm hội tụ chính xác trong suốt timeline nghiên cứu 5 tháng.

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ố đã mang lại 4 phát hiện cốt lõi về mặt cơ học:

Thứ nhất, đối với bài toán phân tích tĩnh, độ võng lớn nhất theo phương z của vỏ Mindlin tỷ lệ thuận phi tuyến với chiều dài vết nứt. Khi chiều dài vết nứt ở giữa vỏ tăng từ 0,1 m lên 0,4 m, độ võng cực đại tại tâm vỏ tăng trưởng liên tục. Trên lưới 40x40x2 phần tử, giá trị chuyển vị tính bằng phần tử XCS-DSG3 bám sát kết quả phần mềm thương mại ANSYS với độ lệch cực đại luôn duy trì dưới 1,2%.

Thứ hai, đối với bài toán dao động tự do, tần số dao động tự nhiên của cả 5 mode đầu tiên đều suy giảm khi chiều dài vết nứt tăng lên. Cụ thể, khi khảo sát vết nứt ở tâm có chiều dài 0,3 m, tần số dao động của mode 1 giảm khoảng 4,5% so với trạng thái nguyên vẹn, sai số tương đối giữa XCS-DSG3 và ANSYS cho cả 5 mode dao động dao động trong khoảng hẹp từ 0,25% đến 0,85%.

Thứ ba, vị trí vết nứt đóng vai trò quyết định đến mức độ suy giảm độ cứng kết cấu. Vết nứt nằm ở biên kết cấu gây ra sự tập trung biến dạng cục bộ mạnh hơn so với vết nứt ở giữa tấm vỏ cùng chiều dài 0,4 m, làm giảm tần số dao động tự nhiên của các mode bậc cao từ 3% đến 8% so với trường hợp nứt tâm.

Thứ tư, thuật toán thể hiện tính ưu việt vượt trội khi xử lý vết nứt cong phức tạp dọc trục x trên lưới 32x32x2 phần tử. Phần tử XCS-DSG3 mô tả chính xác trường ứng suất và biến dạng uốn mà không cần phần tử kỳ dị đẳng tham số bậc hai hay kỹ thuật chia nhỏ lưới bán kính quanh đỉnh nứt.

Thảo luận kết quả

Nguyên nhân vật lý của các hiện tượng trên bắt nguồn từ việc xuất hiện vết nứt làm gián đoạn trường ứng suất màng và ứng suất uốn, dẫn đến sự suy giảm độ cứng tổng thể của ma trận K trong phương trình dao động (K - ω^2 M)d = 0. Do ma trận khối lượng M gần như không đổi, sự suy giảm của ma trận độ cứng K tất yếu dẫn đến sự sụt giảm của các giá trị riêng ω (tần số dao động tự nhiên).

Các dữ liệu số trong nghiên cứu được trực quan hóa thông qua bảng so sánh sai số và hệ thống đồ thị phân tán 2D biểu diễn biến thiên độ võng, biến thiên tần số 5 mode đầu theo chiều dài nứt d. Dạng dao động riêng được mô tả sinh động qua các biểu đồ phân bố biên độ chuyển vị 3D và các đường cắt mặt ngang dọc theo bán kính và chiều dài của vỏ.

So sánh với các nghiên cứu sử dụng phần tử vỏ truyền thống hoặc phương pháp nút ảo (PNM), kỹ thuật trơn hóa biến dạng CS-DSG3 giúp loại bỏ hoàn toàn độ cứng giả tạo do khóa cắt, trong khi các hàm làm giàu của XFEM tái hiện chính xác điểm kỳ dị ứng suất tại đỉnh vết nứt. Điều này chứng minh rằng phần tử tam giác 3 nút XCS-DSG3 vừa tiết kiệm tài nguyên tính toán vừa đạt độ tin cậy tương đương các phần tử bậc cao 8 nút hay 9 nút phức tạp.

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

Dựa trên các kết quả nghiên cứu đạt được, 4 giải pháp hành động cụ thể được đề xuất nhằm ứng dụng và phát triển phương pháp:

Thứ nhất, tích hợp phần tử vỏ phẳng XCS-DSG3 vào các phần mềm tính toán kết cấu mã nguồn mở và thương mại. Nhóm phát triển phần mềm kỹ thuật cần xây dựng module phần tử tam giác 3 nút tích hợp hàm làm giàu Heaviside và hàm tiệm cận, đặt mục tiêu giảm từ 40% đến 50% thời gian xử lý tiền kỳ và chia lưới trong vòng 6 đến 12 tháng tới.

Thứ hai, ứng dụng thuật toán vào quy trình kiểm định và chẩn đoán hư hỏng kết cấu công trình thực tế. Các kỹ sư kết cấu tại các doanh nghiệp tư vấn và thi công cần sử dụng dữ liệu suy giảm tần số dao động riêng 5 mode đầu để phát hiện sớm các vết nứt có kích thước từ 0,05 m trên bồn bể chứa áp lực, vỏ tàu thủy và silo công nghiệp, triển khai thí điểm trong thời gian từ 3 đến 6 tháng.

Thứ ba, chuẩn hóa quy trình đánh giá an toàn động lực học công trình dạng vỏ mỏng. Viện nghiên cứu cơ học và các cơ quan quản lý chuyên ngành xây dựng cần xây dựng bộ hướng dẫn tính toán dao động vỏ có vết nứt dựa trên phương pháp số mở rộng, hướng đến mục tiêu kiểm soát sai số tần số dao động dưới ngưỡng 1,0% trong vòng 12 tháng.

Thứ tư, mở rộng hướng nghiên cứu thuật toán 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à nghiên cứu và nghiên cứu sinh tại các trường đại học kỹ thuật cần phối hợp thực hiện đề tài phát triển phần tử XCS-DSG3 cho vỏ composite lớp và vật liệu FGM có cơ tính biến thiên, nâng cao hiệu năng phân tích phi tuyến hình học thêm 30% trong lộ trình từ 18 đến 24 tháng.

Đố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ị thiết thực cho 4 nhóm đối tượng cụ thể:

  1. Kỹ sư thiết kế và thẩm tra kết cấu công trình: Tiếp cận công cụ tính toán kiểm tra độ võng và dao động của kết cấu mái vỏ mỏng, vòm bê tông cốt thép, bồn chứa khí hóa lỏng và thân tàu biển có khuyết tật nứt, ứng dụng trực tiếp vào việc đánh giá an toàn chịu lực của công trình mà không cần phần mềm thương mại đắt tiền.
  2. Học viên cao học và nghiên cứu sinh ngành Cơ kỹ thuật - Xây dựng: Khai thác tài liệu tham khảo chất lượng cao về phương pháp phần tử hữu hạn trơn (S-FEM) kết hợp phương pháp phần tử hữu hạn mở rộng (XFEM), hỗ trợ xây dựng thuật toán và mã nguồn MATLAB phục vụ các đề tài nghiên cứu chuyên sâu.
  3. Giảng viên và nhà khoa học tại các trường đại học khối kỹ thuật: Sử dụng làm học liệu giảng dạy chuyên đề Cơ học rạn nứt công trình, Phương pháp số nâng cao trong cơ học kết cấu và Lý thuyết tấm vỏ Mindlin cho sinh viên và học viên sau đại học.
  4. Chuyên gia phát triển phần mềm mô phỏng (CAE): Tham khảo cấu trúc thuật toán ma trận độ cứng trơn, phương pháp khử khóa cắt và kỹ thuật xử lý suy biến bậc tự do xoay để tối ưu hóa bộ giải (solver) trong các phần mềm tính toán kết cấu thế hệ mới.

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

Phần tử XCS-DSG3 khắc phục hiện tượng khóa cắt trong vỏ Mindlin như thế nào? Hiện tượng khóa cắt thường xảy ra khi tỷ số chiều dày trên nhịp giảm xuống dưới 0,01. XCS-DSG3 chia mỗi tam giác thành 3 tam giác con và áp dụng kỹ thuật trơn hóa biến dạng CS-FEM kết hợp toán tử độ lệch trượt DSG3, giúp ma trận biến dạng cắt ổn định và không bị khóa cứng.

Phương pháp XFEM kết hợp Level-set có ưu điểm gì nổi bật so với FEM truyền thống? FEM truyền thống đòi hỏi lưới phần tử phải trùng khớp với biên vết nứt và cần làm mịn cục bộ tại đỉnh nứt với hàng ngàn phần tử nhỏ. XFEM sử dụng hàm Level-set cho phép mô phỏng vết nứt độc lập hoàn toàn với lưới, giảm hơn 60% thời gian phân chia lại phần tử.

Bậc tự do xoay quanh trục pháp tuyến được xử lý như thế nào để tránh suy biến ma trận? Do phần tử tấm phẳng ban đầu không có độ cứng chống xoay quanh trục z, ma trận độ cứng tổng thể có thể bị suy biến khi các phần tử đồng phẳng. Nghiên cứu đã bổ sung độ cứng nhân tạo bằng 10^-10 lần giá trị đường chéo lớn nhất, giúp hệ phương trình giải ổn định tuyệt đối.

Sự hiện diện của vết nứt làm thay đổi tần số dao động tự nhiên của vỏ ra sao? Vết nứt làm suy giảm độ cứng uốn và độ cứng cắt của kết cấu, dẫn đến việc tần số dao động của 5 mode đầu tiên giảm từ 1,5% đến 8% tùy thuộc vào chiều dài vết nứt từ 0,1 m đến 0,4 m và vị trí nứt ở biên hay ở giữa vỏ.

Mức độ chính xác của phương pháp XCS-DSG3 so với phần mềm ANSYS đạt bao nhiêu? Trên các mô hình lưới chuẩn 32x32x2 và 40x40x2 phần tử, kết quả độ võng tĩnh và tần số dao động riêng 5 mode đầu của XCS-DSG3 đạt độ tương đồng rất cao với ANSYS Shell93, sai số tương đối luôn nhỏ hơn 1,5%.

Kết luận

  • Luận văn đã xây dựng thành công phần tử vỏ phẳng XCS-DSG3 kết hợp giữa phương pháp phần tử hữu hạn mở rộng XFEM và kỹ thuật trơn hóa biến dạng phần tử CS-DSG3.
  • Xử lý triệt để hiện tượng khóa cắt và mô tả chính xác trường kỳ dị ứng suất tại đỉnh vết nứt trên phần tử tam giác 3 nút đơn giản.
  • Khảo sát định lượng toàn diện ảnh hưởng của vết nứt thẳng, vết nứt biên và vết nứt cong với chiều dài từ 0,1 m đến 0,4 m đến độ võng và tần số dao động 5 mode đầu của vỏ Mindlin.
  • Kết quả mô phỏng trên nền tảng MATLAB đạt độ chính xác cao với sai số dưới 1,5% so với phần mềm thương mại ANSYS nhưng tiết kiệm đáng kể thời gian chia lưới.
  • Đóng góp giải pháp học thuật mới cho cơ học tính toán và mở ra hướng ứng dụng hiệu quả trong chẩn đoán rạn nứt công trình dạng vỏ.

Trong giai đoạn 6 đến 12 tháng tiếp theo, thuật toán cần được phát triển mở rộng cho các bài toán rạn nứt 3D và vật liệu composite nhiều lớp. Các kỹ sư và nhà nghiên cứu quan tâm có thể khai thác mã nguồn và khung lý thuyết của luận văn để áp dụng vào các dự án kiểm định kết cấu thực tế ngay hôm nay.