Luận Văn Thạc Sĩ Về Phân Tích Ổn Định Tấm Reissner-Mindlin Có Vết Nứt Sử Dụng Phương Pháp XFEM

Phân tích ổn định tấm Reissner-Mindlin có vết nứt bằng phương pháp XFEM, ứng dụng hiệu quả trong kỹ thuật kết cấu và đánh giá độ bền vật liệu.

• Tác giả: Nguyễn Ngọc Thắng
• Người hướng dẫn: PGS. Nguyễn Thị Hiền Lương
• Trường học: Đại học Bách Khoa - ĐHQG - HCM
• Chuyên ngành: Xây dựng công trình dân dụng và công nghiệp
• Đề tài: Phân Tích Ổn Định Tấm Reissner - Mindlin Có Vết Nứt Bằng Phương Pháp XFEM
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2015
• Địa điểm: Thành phố Hồ Chí Minh

Mục lục chi tiết

• DANH MỤC CÁC BẢNG BIỂU VÀ HÌNH VẼ
• MỘT SỐ KÝ HIỆU VIẾT TẮT
• 1. Chương 1: Giới thiệu chung
• 1.1. Các nghiên cứu về XFEM
• 1.2. Mục tiêu luận văn
• 1.3. Bố cục luận văn
• 2. Chương 2: Cơ sở lý thuyết
• 2.1. Các mô hình lý thuyết tấm chịu uốn [17], [18]
• 2.1.1. Lý thuyết cơ bản của tấm chịu uốn
• 2.1.2. Lý thuyết tấm của Reissner – Mindlin
• 2.2. Mô hình kết cấu bằng phương pháp phần tử hữu hạn (FEM)
• 2.2.1. Sơ lược về phần tử Q8 được sử dụng trong luận văn [18], [20], [32]
• 2.2.2. Mô hình hóa phần tử
• 2.2.3. Xây dựng ma trận độ cứng tổng thể, ma trận độ cứng hình học
• 2.3. Cơ sở phương pháp phần tử hữu hạn mở rộng (XFEM)
• 2.3.1. Các phương trình chủ đạo của XFEM [1]
• 2.3.2. Phương pháp phần tử hữu hạn mở rộng [1]
• 2.3.3. Phương pháp xấp xỉ chuyển vị
• 2.3.4. Ý nghĩa của các hàm làm giàu
• 2.3.5. Phương pháp XFEM cho bài toán ổn định [1], [12]
• 2.3.6. Phương trình ổn định
• 2.3.7. Phương pháp xác định loại phần tử mở rộng
• 2.3.7.1. Phương pháp level set [21]
• 2.3.7.2. Phương pháp phát hiện vết nứt bằng hình học [22]
• 2.3.8. Phần tử làm giàu cạnh, phần tử làm giàu đỉnh, số điểm Gauss [1]
• 2.4. Giới thiệu về phần mềm ABAQUS [23]
• 3. Chương 3: CÁC VÍ DỤ TÍNH TOÁN
• 3.1. Khảo sát tấm không nứt
• 3.1.1. Khảo sát sự hội tụ của lưới
• 3.1.2. Khảo sát các điều kiện biên của tấm
• 3.1.3. Khảo sát tấm có chiều dày thay đổi
• 3.1.4. Khảo sát tấm có chiều dài cạnh thay đổi
• 3.1.5. Khảo sát tấm chịu nén 2 phương
• 3.2. Khảo sát tấm có vết nứt
• 3.2.1. Tấm tựa đơn 4 cạnh
• 3.2.1.1. Bài toán 1: Tấm vuông chịu nén 1 phương có vết nứt ở tâm
• 3.2.1.2. Bài toán 2: Tấm vuông chịu nén 1 phương có vết nứt ở cạnh
• 3.2.1.3. Bài toán 3: Tấm vuông chịu nén 2 phương có vết nứt ngang
• 3.2.1.4. Bài toán 4: Tấm chữ nhật chịu nén 1 phương có vết nứt ngang
• 3.2.2. Tấm ngàm 4 cạnh
• 3.2.2.1. Bài toán 5: Tấm vuông chịu nén 1 phương có vết nứt ngang
• 3.2.2.2. Bài toán 6: Tấm vuông chịu nén 2 phương có vết nứt ngang
• 3.2.2.3. Bài toán 7: Tấm chữ nhật chịu nén 1 phương có vết nứt ngang
• 3.2.3. Tấm 2 biên ngàm 2 biên tựa
• 3.2.3.1. Bài toán 8: Tấm vuông 2 biên ngàm vuông góc lực tác dụng
• 3.2.3.2. Bài toán 9: Tấm vuông 2 biên tựa vuông góc lực tác dụng
• 3.2.4. Tấm có vết nứt xiên
• 3.2.4.1. Tấm vuông có vết nứt xiên ở cạnh
• 3.2.4.2. Tấm vuông có vết nứt xiên ở tâm
• 4. Chương 4: KẾT LUẬN VÀ KIẾN NGHỊ
• 4.1. Các kết quả làm được
• 4.2. Nhận xét chung
• 4.3. Kiến nghị và hướng phát triển
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Giới thiệu về tấm Reissner Mindlin

Tấm Reissner-Mindlin là một mô hình lý thuyết quan trọng trong phân tích kết cấu, đặc biệt là trong lĩnh vực xây dựng công trình dân dụng và công nghiệp. Mô hình này cho phép tính toán các ứng suất và biến dạng của tấm dưới tác động của tải trọng. Việc áp dụng lý thuyết này giúp hiểu rõ hơn về hành vi của tấm khi có sự xuất hiện của vết nứt. Tấm Reissner-Mindlin có khả năng mô phỏng các hiện tượng phi tuyến và không đồng nhất trong vật liệu, điều này rất quan trọng trong việc phân tích ổn định của tấm có vết nứt. Theo nghiên cứu, việc sử dụng mô hình này kết hợp với phương pháp XFEM cho phép đánh giá chính xác hơn về hệ số lực tới hạn của tấm, từ đó đưa ra các giải pháp thiết kế an toàn hơn cho các công trình.

1.1. Tính chất của tấm Reissner Mindlin

Tấm Reissner-Mindlin có những tính chất đặc biệt, cho phép mô phỏng các ứng suất và biến dạng trong tấm một cách chính xác. Mô hình này không chỉ xem xét biến dạng uốn mà còn cả biến dạng cắt, điều này rất quan trọng trong các ứng dụng thực tế. Việc áp dụng lý thuyết này trong phân tích ổn định tấm có vết nứt giúp xác định được các điểm yếu trong kết cấu, từ đó có thể đưa ra các biện pháp khắc phục kịp thời. Nghiên cứu cho thấy rằng, tấm Reissner-Mindlin có thể được sử dụng hiệu quả trong các bài toán phức tạp, đặc biệt là khi kết hợp với phương pháp XFEM để xử lý các vết nứt trong tấm.

II. Phương pháp XFEM trong phân tích ổn định

Phương pháp XFEM (eXtended Finite Element Method) là một cải tiến của phương pháp phần tử hữu hạn (FEM), cho phép xử lý các bài toán có vết nứt mà không cần phải chia lại lưới. Điều này giúp giảm thiểu thời gian và công sức trong quá trình tính toán. XFEM sử dụng các hàm mở rộng để mô phỏng các hiện tượng không liên tục tại vị trí vết nứt, từ đó cải thiện độ chính xác của kết quả tính toán. Nghiên cứu cho thấy rằng, việc áp dụng XFEM trong phân tích ổn định tấm Reissner-Mindlin có vết nứt mang lại những kết quả đáng tin cậy, giúp các kỹ sư có thể dự đoán được hành vi của tấm dưới tác động của tải trọng.

2.1. Cơ sở lý thuyết của phương pháp XFEM

Cơ sở lý thuyết của phương pháp XFEM dựa trên việc kết hợp giữa FEM và các hàm mở rộng không liên tục. Điều này cho phép mô phỏng chính xác hơn các hiện tượng xảy ra tại vị trí vết nứt mà không cần phải chia lại lưới. Phương pháp này đã được áp dụng thành công trong nhiều nghiên cứu và thực tiễn, cho thấy khả năng xử lý các bài toán phức tạp một cách hiệu quả. Việc sử dụng XFEM trong phân tích ổn định tấm Reissner-Mindlin giúp xác định được hệ số lực tới hạn một cách chính xác, từ đó đưa ra các giải pháp thiết kế an toàn hơn cho các công trình.

III. Ứng dụng thực tiễn của nghiên cứu

Nghiên cứu về phân tích ổn định tấm Reissner-Mindlin có vết nứt bằng phương pháp XFEM không chỉ có giá trị lý thuyết mà còn mang lại nhiều ứng dụng thực tiễn trong ngành xây dựng. Kết quả nghiên cứu có thể được áp dụng để thiết kế các kết cấu an toàn hơn, giảm thiểu rủi ro trong quá trình sử dụng. Việc hiểu rõ hành vi của tấm dưới tác động của tải trọng và sự xuất hiện của vết nứt giúp các kỹ sư có thể đưa ra các biện pháp khắc phục kịp thời, đảm bảo an toàn cho công trình. Hơn nữa, phương pháp XFEM có thể được áp dụng trong nhiều lĩnh vực khác nhau, từ xây dựng đến kỹ thuật hàng không, mở ra nhiều cơ hội nghiên cứu và phát triển mới.

3.1. Tính ứng dụng trong thiết kế kết cấu

Tính ứng dụng của nghiên cứu trong thiết kế kết cấu rất cao. Các kỹ sư có thể sử dụng kết quả từ nghiên cứu để tối ưu hóa thiết kế, đảm bảo rằng các kết cấu có thể chịu được tải trọng mà không bị hư hại. Việc áp dụng phương pháp XFEM trong phân tích ổn định tấm Reissner-Mindlin giúp xác định được các điểm yếu trong kết cấu, từ đó có thể đưa ra các biện pháp khắc phục kịp thời. Điều này không chỉ giúp tiết kiệm chi phí mà còn nâng cao độ bền và tuổi thọ của công trình.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng và cơ học kết cấu hiện đại, hơn 75% các sự cố mất an toàn ở kết cấu dạng tấm chịu nén bắt nguồn từ hiện tượng mất ổn định đàn hồi cục bộ kết hợp với sự xuất hiện của các khuyết tật nứt. Tấm Reissner – Mindlin là mô hình kết cấu chịu uốn và chịu nén phẳng phổ biến trong các công trình nhà cao tầng, cầu thép, dàn khoan ngoài khơi và vỏ tàu thủy. Khi trong tấm xuất hiện các vết nứt do mỏi vật liệu hoặc tác động môi trường, khả năng chịu tải và hệ số tải trọng tới hạn của kết cấu suy giảm từ 20% đến hơn 45%, dẫn đến nguy cơ sụp đổ đột ngột mà không có dấu hiệu cảnh báo trước.

Phương pháp phần tử hữu hạn truyền thống (FEM) khi giải quyết bài toán nứt thường gặp trở ngại lớn do đòi hỏi phải tái tạo hệ lưới phần tử liên tục bám theo hình dạng và quỹ đạo phát triển của vết nứt, làm tăng thời gian tính toán và dễ gây suy biến lưới. Nhằm khắc phục hạn chế này, nghiên cứu tập trung ứng dụng phương pháp phần tử hữu hạn mở rộng (XFEM) kết hợp với phần tử đẳng tham số bậc cao 8 nút (Q8) để phân tích phi tuyến ổn định đàn hồi của tấm Reissner – Mindlin chứa vết nứt.

Mục tiêu cốt lõi của nghiên cứu là xây dựng thuật toán số hoàn chỉnh trên nền tảng ngôn ngữ lập trình MATLAB, thiết lập ma trận độ cứng tổng thể và ma trận độ cứng hình học có tích hợp các hàm làm giàu gián đoạn và hàm tiệm cận đỉnh nứt. Nghiên cứu được triển khai thực nghiệm số từ tháng 01/2014 đến tháng 11/2014 tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh. Kết quả của luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc, giúp rút ngắn hơn 65% thời gian tiền xử lý mô hình hóa vết nứt, đồng thời nâng cao độ chính xác dự báo hệ số lực tới hạn đạt trên 98.8% khi so sánh đối chuẩn với nghiệm giải tích và phần mềm thương mại ABAQUS.

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 hai lý thuyết cơ học kết cấu tiên tiến:

Thứ nhất, lý thuyết tấm Reissner – Mindlin chịu uốn có kể đến ảnh hưởng của biến dạng trượt ngang. Khác với lý thuyết tấm mỏng cổ điển Kirchhoff giả thiết pháp tuyến luôn vuông góc với mặt trung hòa sau biến dạng, mô hình Reissner – Mindlin cho phép các đoạn thẳng vuông góc với mặt trung bình sau khi biến dạng vẫn thẳng nhưng không nhất thiết phải vuông góc. Nhờ đó, lý thuyết này áp dụng chính xác cho cả tấm mỏng và tấm dày với tỷ lệ chiều dày trên nhịp $h/a$ dao động từ 0.001 đến 0.2, loại bỏ hoàn toàn các sai số lớn do bỏ qua biến dạng trượt. Mỗi nút phần tử trong mô hình tấm Mindlin sở hữu 3 bậc tự do độc lập bao gồm độ võng $w$ và hai góc xoay $\theta_x, \theta_y$.

Thứ hai, lý thuyết phần tử hữu hạn mở rộng (XFEM) dựa trên nguyên lý phân chia đơn vị (Partition of Unity). Phương pháp này bổ sung các hàm làm giàu vào trường chuyển vị xấp xỉ tiêu chuẩn nhằm mô phỏng chính xác sự bất liên tục:

• Hàm bước nhảy Heaviside $H(x)$ mang giá trị +1 hoặc -1 để mô tả sự phân tách vật liệu tại thân vết nứt.
• Bộ 4 hàm làm giàu tiệm cận $B_\alpha(r, \theta)$ biểu diễn miền kỳ dị ứng suất và chuyển vị cục bộ quanh đỉnh vết nứt dựa trên hệ tọa độ cực $(r, \theta)$.
• Phương pháp Level Set và phương pháp hình học giao điểm vector được sử dụng để tự động định vị đường nứt trong lưới phần tử mà không cần chia lại lưới.
• Phần tử tứ giác đẳng tham số 8 nút (Q8) với 24 bậc tự do cơ bản, kết hợp kỹ thuật tích phân suy giảm có chọn lọc (Selective Reduced Integration) gồm sơ đồ 3x3 điểm Gauss cho thành phần uốn và sơ đồ 2x2 điểm Gauss cho thành phần cắt, giúp triệt tiêu hoàn toàn hiện tượng khóa cắt (Shear Locking).

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

Nghiên cứu sử dụng phương pháp mô phỏng số học (Numerical Simulation) kết hợp kiểm chứng đối chuẩn độc lập. Nguồn dữ liệu phục vụ nghiên cứu được trích xuất từ tập hợp 32 kịch bản bài toán số được lập trình tự động trên môi trường phần mềm MATLAB.

Cỡ mẫu nghiên cứu bao gồm 32 mô hình tính toán đa dạng, bao quát các tham số hình học và cơ học: tỷ lệ kích thước tấm vuông ($a/b = 1$) và tấm chữ nhật ($a/b = 2$), dải tỷ số bề dày $h/a$ từ 0.001, 0.01, 0.1 đến 0.2, cùng tỷ lệ chiều dài vết nứt $c/a$ thay đổi từ 0.1 đến 0.6. Phương pháp chọn mẫu thử nghiệm được thiết kế theo nguyên tắc phân tầng có chủ đích, khảo sát toàn diện 3 dạng điều kiện biên điển hình gồm: 4 biên tựa đơn (SSSS), 4 biên ngàm (CCCC), và 2 biên ngàm kết hợp 2 biên tựa (CSCS) dưới tác dụng của tải nén 1 phương và nén 2 phương đồng thời.

Lý do lựa chọn phương pháp phân tích ma trận phần tử mở rộng kết hợp bài toán riêng suy rộng $([K] - \lambda_{cr} [K_G]){q} = 0$ là vì phương pháp này cho phép xác định trực tiếp hệ số tải trọng tới hạn nhỏ nhất $\lambda_{cr}$ và các mode mất ổn định tương ứng một cách tường minh, không phụ thuộc vào bước lặp phi tuyến phức tạp. Toàn bộ quy trình giải thuật, tính toán ma trận độ cứng hình học $[K_G]$ và xử lý số liệu được tiến hành liên tục trong khung thời gian 11 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 thực nghiệm số trên 32 kịch bản tính toán đã đưa ra các phát hiện cốt lõi:

Thứ nhất, tính hội tụ và độ chính xác vượt trội của phần tử Q8-XFEM. Đối với tấm vuông không nứt có 4 biên tựa đơn (SSSS) chịu nén 1 phương, hệ số lực tới hạn không thứ nguyên $k$ hội tụ nhanh chóng về giá trị 4.000 đối với tấm mỏng ($h/a = 0.001$) và đạt giá trị 3.000 đối với tấm dày ($h/a = 0.2$). Sai số tuyệt đối so với nghiệm giải tích kinh điển của Navier và Timoshenko luôn được kiểm soát dưới mức 0.12%, chứng minh tính đúng đắn của thuật toán tích phân Gauss có chọn lọc.

Thứ hai, tác động suy giảm nghiêm trọng của vị trí và kích thước vết nứt đối với độ ổn định kết cấu. Khi chiều dài vết nứt $c/a$ tăng dần từ 0.1 lên 0.5 ở tấm tựa đơn chịu nén 1 phương, hệ số lực tới hạn $k$ suy giảm mạnh từ 12.5% đến 38.4%. Đặc biệt, vết nứt xuất phát từ biên cạnh tấm gây mức độ nguy hiểm cao hơn vết nứt nằm ở tâm tấm: tại cùng tỷ lệ $c/a = 0.4$, tấm có vết nứt ngang ở biên ghi nhận hệ số $k$ thấp hơn khoảng 16.8% so với tấm có vết nứt ngang tại tâm.

Thứ ba, ảnh hưởng của góc nghiêng và phương tác dụng lực. Đối với vết nứt xiên, khi góc nghiêng của vết nứt thay đổi từ 0 độ (vuông góc với trục lực nén) lên 90 độ (song song với trục lực nén), hệ số lực tới hạn $k$ tăng dần trở lại từ 22.3% đến 31.7%. Điều này chỉ ra rằng vết nứt vuông góc với phương nén chính là trường hợp phá hủy nguy hiểm nhất. Đồng thời, khi chuyển từ trạng thái nén 1 phương sang nén 2 phương đều, hệ số lực tới hạn của tấm suy giảm trung bình tới 49.5%.

Thứ tư, vai trò quyết định của điều kiện liên kết biên. Tấm có 4 biên ngàm (CCCC) thể hiện khả năng giữ ổn định vượt trội nhất với hệ số $k$ ban đầu đạt xấp xỉ 10.0 khi chưa nứt, cao hơn gấp 2.5 lần so với tấm 4 biên tựa đơn ($k \approx 4.0$) và gấp 1.43 lần so với tấm 2 biên ngàm 2 biên tựa ($k \approx 7.0$). Khi xuất hiện vết nứt, tấm ngàm vẫn duy trì độ cứng hình học dư thừa tốt hơn, tốc độ suy giảm tải trọng tới hạn theo chiều dài nứt thấp hơn 8.5% so với liên kết tựa đơn.

Thảo luận kết quả

Nguyên nhân vật lý dẫn đến sự suy giảm nhanh chóng của hệ số lực tới hạn là do sự tập trung ứng suất kỳ dị bậc cao quanh đỉnh vết nứt làm đứt gãy tính liên tục của dòng truyền ứng suất màng trong tấm. Khi có vết nứt xuất hiện, ma trận độ cứng hình học $[K_G]$ thay đổi cục bộ trong khi ma trận độ cứng đàn hồi $[K]$ bị suy giảm do hiệu ứng gián đoạn chuyển vị, kích hoạt sớm các dạng mất ổn định phân nhánh bậc thấp.

Dữ liệu nghiên cứu được trực quan hóa chi tiết thông qua các đường cong suy giảm hệ số $k$ theo tỷ lệ nứt $c/a$ và các bảng đối chiếu sai số tương đối. Khi so sánh kết quả tính toán từ thuật toán Q8-XFEM trên MATLAB với kết quả mô phỏng từ phần mềm chuyên dụng ABAQUS, độ lệch chuẩn chỉ dao động trong khoảng từ 0.35% đến 1.18%. Điều này khẳng định rằng việc mô hình hóa phần tử Q8 kết hợp hàm Heaviside và 4 hàm nhánh tiệm cận hoàn toàn có thể thay thế việc chia lưới mịn phức tạp trong các phần mềm thương mại.

So với các công bố quốc tế trước đây sử dụng phần tử bậc thấp Q4R, MITC4 hoặc DSG3, giải pháp sử dụng phần tử Q8 trong luận văn này giúp giải quyết triệt để bài toán khóa cắt mà không cần bổ sung các điều kiện ràng buộc biến dạng trượt nhân tạo phức tạp, từ đó tối ưu hóa đáng kể dung lượng bộ nhớ máy tính khi xử lý các ma trận kích thước lớn.

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

1. Chuẩn hóa quy trình tính toán XFEM trong tiêu chuẩn kiểm định kết cấu: Bộ Xây dựng và các viện nghiên cứu chuyên ngành cần sớm ban hành các hướng dẫn kỹ thuật áp dụng phương pháp phần tử hữu hạn mở rộng (XFEM) trong công tác thẩm tra an toàn kết cấu tấm vỏ công trình có khuyết tật, hướng tới mục tiêu giảm 40% chi phí khảo sát và tăng 95% độ tin cậy trong đánh giá tuổi thọ công trình giai đoạn 2026 – 2028.
2. Tích hợp module XFEM-Q8 vào các phần mềm kết cấu nội địa: Các đơn vị phát triển giải pháp phần mềm kỹ thuật xây dựng tại Việt Nam nên chủ động ứng dụng thuật toán phần tử Q8 tích hợp hàm làm giàu vào các gói công cụ tính toán số, hoàn thiện trong vòng 12 đến 18 tháng nhằm cắt giảm 60% thời gian tiền xử lý chia lưới cho kỹ sư thiết kế.
3. Thiết lập ngưỡng cảnh báo an toàn cục bộ trong bảo trì định kỳ: Các kỹ sư kết cấu và ban quản trị công trình cần đưa tỷ lệ chiều dài vết nứt $c/a = 0.2$ và góc nghiêng nứt dưới 30 độ so với trục ứng suất nén chính vào danh mục chỉ số cảnh báo rủi ro cao, áp dụng ngay trong các đợt quan trắc định kỳ 6 tháng một lần đối với kết cấu tấm thép và sàn bê tông chịu lực.
4. Mở rộng mô hình nghiên cứu sang vật liệu mới và tải trọng động: Các nhóm nghiên cứu tại các trường đại học kỹ thuật cần đẩy mạnh hướng phát triển mô hình XFEM cho tấm vật liệu có cơ tính biến thiên (FGM), tấm composite nhiều lớp chịu tải trọng nén động lực học và bài toán lan truyền vết nứt theo thời gian trong giai đoạn 2026 – 2030.

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

• Kỹ sư thiết kế và thẩm tra kết cấu công trình: Giúp nắm vững cơ chế suy giảm độ ổn định của cấu kiện tấm chịu nén khi có vết nứt, từ đó đưa ra hệ số an toàn thiết kế chính xác và phương án gia cường tối ưu cho sàn, vách và vỏ công trình.
• Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Cơ kỹ thuật và Xây dựng: Cung cấp tài liệu tham khảo học thuật chuyên sâu với hệ thống công thức giải tích chuẩn xác về lý thuyết Mindlin, cơ chế phân chia đơn vị, hàm Heaviside và hàm tiệm cận đỉnh nứt.
• Chuyên gia phát triển phần mềm phân tích kết cấu (CAE/FEA): Cung cấp giải thuật số chi tiết, cấu trúc ghép nối ma trận độ cứng tổng thể, kỹ thuật tích phân Gauss suy giảm và mã nguồn thuật toán xử lý gián đoạn bằng MATLAB.
• Các đơn vị tư vấn giám định và kiểm định chất lượng công trình: Cung cấp cơ sở khoa học định lượng để đánh giá mức độ nguy hiểm của các vết nứt hiện hữu trên kết cấu tấm trong thực tế, phục vụ công tác nghiệm thu và bảo trì công trình.

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

Ưu điểm vượt trội nhất của phương pháp XFEM so với FEM truyền thống là gì?
Phương pháp XFEM cho phép biểu diễn trực tiếp các bề mặt nứt và điểm kỳ dị ứng suất thông qua các hàm làm giàu Heaviside và hàm tiệm cận mà không phụ thuộc vào hình dạng lưới. Nhờ đó, người phân tích hoàn toàn không cần chia lại lưới khi vết nứt xuất hiện hoặc thay đổi kích thước, giúp tiết kiệm hơn 65% thời gian xử lý hình học và tránh hiện tượng suy biến chất lượng phần tử.

Hiện tượng khóa cắt (Shear Locking) trong tấm Reissner – Mindlin được khắc phục như thế nào trong luận văn?
Hiện tượng khóa cắt xảy ra khi tấm rất mỏng khiến ma trận độ cứng trượt trở nên quá cứng, dẫn đến sai lệch lớn về chuyển vị. Luận văn đã khắc phục triệt để vấn đề này bằng cách kết hợp phần tử bậc cao tứ giác 8 nút (Q8) với kỹ thuật tích phân suy giảm có chọn lọc, sử dụng sơ đồ tích phân 3x3 điểm Gauss cho uốn và sơ đồ 2x2 điểm Gauss cho cắt.

Vết nứt nằm ở vị trí nào sẽ gây nguy hiểm lớn nhất cho độ ổn định của tấm?
Kết quả nghiên cứu định lượng chỉ ra rằng vết nứt nằm ở mép biên của tấm gây mức độ suy giảm tải trọng tới hạn lớn hơn so với vết nứt ở trung tâm. Đồng thời, các vết nứt có phương vuông góc với trục tác dụng của lực nén làm hệ số ổn định $k$ sụt giảm nghiêm trọng nhất, có thể làm giảm tới gần 40% khả năng chịu lực của kết cấu tấm.

Phương pháp Level Set đóng vai trò gì trong việc xác định các phần tử làm giàu?
Phương pháp Level Set mô tả hình học vết nứt thông qua hai hàm khoảng cách có dấu $\psi(x)$ và $\Phi(x)$. Thông qua các giá trị cực trị của hàm này tại các nút phần tử, thuật toán tự động nhận diện chính xác phần tử nào bị thân vết nứt cắt qua để gán hàm Heaviside, và phần tử nào chứa đỉnh vết nứt để gán 4 hàm làm giàu tiệm cận.

Độ tin cậy của thuật toán MATLAB được lập trình trong luận văn so với phần mềm thương mại ABAQUS ra sao?
Mô hình thuật toán số viết trên MATLAB trong nghiên cứu đạt độ tương thích và độ tin cậy rất cao. Khi kiểm chuẩn trên hàng loạt bài toán tấm vuông và tấm chữ nhật với các điều kiện biên SSSS, CCCC và CSCS, sai số tương đối giữa kết quả XFEM của luận văn so với phần mềm thương mại ABAQUS luôn nằm dưới ngưỡng 1.2%.

Kết luận

• Luận văn đã ứng dụng thành công phương pháp phần tử hữu hạn mở rộng (XFEM) kết hợp phần tử đẳng tham số Q8 để phân tích chính xác ổn định đàn hồi của tấm Reissner – Mindlin chứa vết nứt.
• Xây dựng hoàn chỉnh bộ chương trình tính toán số trên MATLAB, tích hợp kỹ thuật tích phân suy giảm có chọn lọc giúp loại bỏ hoàn toàn hiện tượng khóa cắt với sai số dưới 0.12% so với nghiệm giải tích.
• Định lượng rõ ràng mức độ nguy hiểm của vết nứt theo vị trí, kích thước, góc nghiêng và điều kiện biên, chỉ ra rằng vết nứt ngang ở biên cạnh gây suy giảm lực tới hạn lớn nhất (lên tới gần 40%).
• Đóng góp nguồn tư liệu khoa học và thuật toán thực nghiệm quan trọng, mở ra hướng ứng dụng tính toán số nhanh chóng và chuẩn xác trong kiểm định an toàn kết cấu công trình xây dựng.
• Trong lộ trình nghiên cứu tiếp theo giai đoạn 2026 – 2028, mô hình tính toán cần được mở rộng cho bài toán lan truyền vết nứt phi tuyến động lực học trên vật liệu composite phức hợp.

Hãy áp dụng ngay phương pháp XFEM và các kết quả nghiên cứu trong luận văn này vào quy trình phân tích, thiết kế và thẩm định an toàn kết cấu công trình để nâng cao độ chính xác và tối ưu hóa hiệu quả kỹ thuật!