Phân Tích Kết Cấu Tấm Vỏ Bằng Phần Tử Vỏ Khối Có Trường Biến Dạng Trơn (CSn-S8)

Tóm tắt nghiên cứu (200-250 từ)

Nghiên cứu này giới thiệu công thức phần tử hữu hạn (PTHH) vỏ khối 8 nút cải tiến hoàn toàn mới mang tên CSn-S8 ($n = 1, 2, 3, 4$), được phát triển bởi nhóm nghiên cứu tại Trường Đại học Sư phạm Kỹ thuật TP.HCM do PGS. Châu Đình Thành chủ nhiệm (Mã số đề tài: T2019-73TĐ). Mục tiêu trọng tâm là giải quyết triệt để các hiện tượng khóa số học kinh điển gồm khóa cắt (shear locking), khóa màng (membrane locking) và khóa hình thang (trapezoidal locking) khi phân tích tĩnh các kết cấu tấm và vỏ mỏng trong miền đàn hồi tuyến tính.

Về mặt phương pháp, phần tử vỏ khối CSn-S8 tận dụng hàm dạng liên tục $C^0$ chỉ chứa các bậc tự do chuyển vị thẳng (tịnh tiến). Lần đầu tiên, kỹ thuật xấp xỉ biến dạng uốn giả định (assumed bending strains) vốn dùng cho phần tử vỏ suy biến được tích hợp thành công vào phần tử vỏ khối 8 nút. Đồng thời, kỹ thuật làm trơn biến dạng trên miền phần tử (Cell-based Strain Smoothing - CS) được áp dụng cho biến dạng màng, chuyển đổi tích phân diện tích phức tạp thành tích phân đường biên thông qua định lý phân kỳ Gauss.

Kết quả thử nghiệm số trên 7 bài toán chuẩn (benchmark) kinh điển khẳng định phần tử CS4-S8 có khả năng khắc phục hoàn toàn hiện tượng khóa cắt ngay cả khi tỉ số độ mảnh $h/L$ đạt ngưỡng cực hạn $10^{-8}$. Sai số tính toán độ võng giảm xuống mức ấn tượng (chỉ 0,11% đối với bài toán tấm hình thoi Morley cực nhọn). Nghiên cứu mang lại một giải pháp tính toán bền vững, chính xác cao và mở ra tiềm năng lớn cho mô phỏng kết cấu composite hay nứt 3D.

Bối cảnh và tầm quan trọng (300-350 từ)

Kết cấu tấm và vỏ là giải pháp hình học và chịu lực tối ưu trong nhiều ngành công nghiệp mũi nhọn như xây dựng dân dụng & công nghiệp, hàng không vũ trụ, công nghệ đóng tàu và sản xuất ô tô. Do khả năng chịu tải trọng lớn trên một trọng lượng bản thân nhẹ, việc tính toán chính xác ứng xử biến dạng và phân bố nội lực của tấm vỏ đóng vai trò sống còn đối với độ an toàn và tuổi thọ công trình.

Trong cơ học tính toán, Phương pháp Phần tử Hữu hạn (FEM) là công cụ chủ lực. Các công thức phần tử vỏ truyền thống thường được chia thành 3 nhóm:

  1. Phần tử vỏ phẳng (Flat shell): Kết hợp giữa phần tử tấm chịu uốn và màng chịu kéo nén, bỏ qua sự tương tác phi tuyến uốn - màng phức tạp.
  2. Phần tử vỏ 3 chiều suy biến (Degenerated continuum shell): Giả thiết biến dạng pháp vuông góc bằng 0, đòi hỏi cả bậc tự do chuyển vị và xoay tại nút, gây khó khăn khi ghép nối với các phần tử khối 3D thông thường.
  3. Phần tử vỏ khối (Solid-shell): Được xây dựng hoàn toàn từ lý thuyết đàn hồi 3 chiều và chỉ sử dụng các bậc tự do chuyển vị thẳng. Nhờ đó, vỏ khối dễ dàng ghép nối với các phần tử cơ học khác và thuận lợi mô phỏng sự biến thiên vật liệu theo chiều dày.

Tuy nhiên, rào cản lớn nhất của phần tử vỏ khối 8 nút $C^0$ là hiện tượng khóa số học (numerical locking) nghiêm trọng khi kết cấu có chiều dày rất mỏng hoặc độ cong lớn, khiến mô hình tính toán trở nên quá cứng (over-stiff) so với thực tế. Mặc dù các kỹ thuật làm trơn biến dạng (Smoothed Finite Element Method - S-FEM) đã phát triển mạnh mẽ cho phần tử vỏ phẳng tại Việt Nam và thế giới, việc tích hợp kỹ thuật làm trơn lên phần tử vỏ khối 8 nút vẫn còn là một khoảng trống nghiên cứu lớn.

Đề tài T2019-73TĐ được thực hiện nhằm lấp đầy khoảng trống này, mang lại một công cụ số tiên tiến, giúp các kỹ sư và nhà nghiên cứu tiếp cận phương pháp tính toán kết cấu tấm vỏ mỏng và dày có độ tin cậy vượt bậc.

Methodology và approach (350-400 từ)

Nghiên cứu áp dụng quy trình thiết lập giải thuật toán học chặt chẽ kết hợp lập trình mô phỏng số trên nền tảng MATLAB, đối soát trực tiếp với các nghiệm giải tích và nghiệm số tham khảo uy tín trên thế giới.

                          Quy trình công thức hóa phần tử CSn-S8

1. Mô hình hình học và trường chuyển vị

Phần tử vỏ khối 8 nút được định nghĩa qua hệ tọa độ tự nhiên $rst$, trong đó mặt trung bình phần tử được tham số hóa bởi các hàm dạng đẳng tham số $C^0$ bậc nhất ($N_I, I=1..4$). Bậc tự do tại mỗi nút chỉ bao gồm 3 chuyển vị thẳng toàn cục $[U, V, W]^T$, loại bỏ hoàn toàn các bậc tự do góc xoay phức tạp.

2. Kỹ thuật xấp xỉ trường biến dạng và xử lý khóa

  • Khử khóa cắt (Shear Locking): Ứng dụng kỹ thuật nội suy ten-xơ hỗn hợp MITC4 (Bathe & Dvorkin), tính toán biến dạng cắt ngoài mặt phẳng tại 4 điểm buộc độc lập $A_1(0,-1,0)$, $A_2(1,0,0)$, $A_3(0,1,0)$, $A_4(-1,0,0)$.
  • Khử khóa hình thang (Trapezoidal Locking): Xấp xỉ lại trường biến dạng thẳng vuông góc với mặt trung bình $\varepsilon_{zz}$ tại 4 điểm buộc $B_1(-1,-1,0)$, $B_2(1,-1,0)$, $B_3(1,1,0)$, $B_4(-1,1,0)$ để khử độ cứng ảo sinh ra do độ cong của vỏ.
  • Xấp xỉ biến dạng uốn giả định: Kế thừa và cải tiến kỹ thuật của Choi và Paik (1994), tái lập trường biến dạng uốn thông qua các điểm buộc, tăng độ chính xác của các bài toán chịu uốn thuần túy.

3. Kỹ thuật làm trơn biến dạng màng trên miền phần tử (CS)

Mặt trung bình của phần tử được chia thành $n$ miền con làm trơn ($n = 1, 2, 3, 4$). Áp dụng định lý phân kỳ Gauss: $$\int_{\Omega_{cs}} N_{I,x} , d\Omega = \oint_{\Gamma_{cs}} N_I n_x , d\Gamma$$ Tích phân diện tích của gradient biến dạng màng được chuyển hóa thành tích phân đường trên các cạnh biên $\Gamma_{cs}$ của từng miền con với 1 điểm Gauss trên mỗi cạnh. Kỹ thuật này giúp ma trận độ cứng màng trở nên mềm dẻo hơn, loại bỏ hiện tượng khóa màng một cách tự nhiên.

Phát hiện chính (400-450 từ)

Phần tử đề xuất CSn-S8 đã được kiểm chứng thông qua chuỗi bài toán kiểm tra độ chính xác, tốc độ hội tụ và độ ổn định:

Bài toán kiểm chuẩn Đặc điểm hình học & Tải trọng Kết quả của CS4-S8 so với các phần tử khác
1. Tấm kéo nén Cook Chịu lực màng thuần túy, lưới méo Cho kết quả chuyển vị tiệm cận nghiệm chính xác tốt hơn phần tử vỏ khối Xsolid85 và RH8s-4.
2. Tấm vuông ngàm chịu uốn Thử nghiệm độ mảnh $h/L = 10^{-3} \to 10^{-8}$ Không bị khóa cắt ngay cả khi $h/L = 10^{-8}$ (độ võng chuẩn hóa đạt 0,1074, trong khi MITC3+, MITC4 bị khóa cứng hoàn toàn tại $10^{-6}$).
3. Tấm chữ nhật chịu lực tập trung Tấm ngàm 4 cạnh, tỉ lệ cạnh $1:5$ Lưới $16 \times 16$ đạt sai số tương đối $< 0,1%$ (tấm vuông) và $< 1,0%$ (tấm chữ nhật).
4. Tấm thoi Morley ($30^\circ$) Góc nhọn méo lưới nghiêm trọng Sai số độ võng tâm chỉ 0,11% (lưới $20 \times 20$), vượt trội hoàn toàn so với MITC4 (7,99%) và Xsolid85 (5,21%).
5. Vỏ mái Scordellis-Lo Tải trọng lượng bản thân CS4-S8 đạt sai số chỉ 0,53% ở lưới $16 \times 16$, hội tụ mượt mà không xảy ra biến dạng giả.
6. Vỏ trụ chịu tải tập trung Bài toán Pinched Cylinder Sai số chỉ 0,94% ở lưới $32 \times 32$; thời gian giải tương đương phần tử vỏ phẳng 4 nút.
7. Vỏ bán cầu có lỗ $18^\circ$ Trạng thái uốn và màng kết hợp phức tạp Đạt độ chính xác cao nhất với sai số chỉ 0,108% tại lưới $32 \times 32$, vượt trội hơn hẳn RH8s-4 và Xsolid85.

Các điểm phát hiện khoa học đột phá:

  • Ảnh hưởng của số miền làm trơn ($n$): Khi tăng số miền làm trơn từ $n=1$ lên $n=4$, độ cứng màng của phần tử được cân bằng tối ưu. Phiên bản CS4-S8 cho kết quả ổn định và chính xác nhất trong mọi trường hợp tải trọng phức hợp.
  • Miễn nhiễm với khóa cắt cực hạn: Phần tử duy trì độ chính xác cao ở dải chiều dày cực mỏng ($h/L = 10^{-8}$), nơi mà hầu hết các phần tử vỏ phẳng chuẩn công nghiệp (như MITC3+, MITC4) đều mất tính hội tụ do khóa cắt.
  • Khả năng chịu méo lưới xuất sắc: Nhờ kỹ thuật làm trơn màng trên biên miền con kết hợp biến dạng uốn giả định, phần tử CSn-S8 thể hiện độ trơ tuyệt vời trước sự suy giảm chất lượng lưới (mesh distortion), đặc biệt là bài toán tấm Morley góc nhọn $30^\circ$.

Đóng góp khoa học (250-300 từ)

Nghiên cứu mang lại những đóng góp học thuật và kỹ thuật mang tính nền tảng:

1. Đóng góp về mặt lý thuyết và phương pháp luận

  • Tiên phong kết hợp lý thuyết: Lần đầu tiên tích hợp thành công kỹ thuật xấp xỉ biến dạng uốn Choi-Paik vào công thức phần tử vỏ khối 8 nút, mở ra hướng tiếp cận mới trong mô hình hóa biến dạng uốn tách biệt mà không cần hàm bậc cao.
  • Hoàn thiện công thức S-FEM cho Solid-Shell: Mở rộng thành công kỹ thuật làm trơn biến dạng ô cơ sở (Cell-based Smoothing) cho phần tử vỏ khối 8 nút, chứng minh tính ưu việt của việc tích phân biên thay cho tích phân diện tích truyền thống.

2. Ứng dụng thực tiễn và chuyển giao

  • Công cụ tính toán đáng tin cậy: Tạo ra bộ giải thuật PTHH có độ chính xác cao, thời gian tính toán tối ưu, phù hợp tích hợp vào các phần mềm phân tích kết cấu công trình xây dựng, cơ khí, hàng không.
  • Tiền đề phát triển vật liệu tiên tiến: Nhờ chỉ sử dụng bậc tự do tịnh tiến và mô phỏng được ứng xử theo phương chiều dày, công thức CSn-S8 là nền tảng hoàn hảo để mở rộng cho:
    • Phân tích tấm vỏ composite nhiều lớp (Laminated composites).
    • Kết cấu vật liệu biến tính chức năng (Functionally Graded Materials - FGMs).
    • Mô phỏng sự lan truyền vết nứt 3D và phân tích phi tuyến hình học lớn.

Đối tượng quan tâm (200-250 từ)

  • Các nhà nghiên cứu và học giả: Chuyên gia trong lĩnh vực Cơ học tính toán (Computational Mechanics), Phương pháp phần tử hữu hạn nâng cao, lý thuyết tấm vỏ và vật liệu mới.
  • Kỹ sư kết cấu & Nhà phát triển phần mềm CAE: Kỹ sư thiết kế trong ngành xây dựng, cơ khí chính xác, hàng không vũ trụ và ô tô đang tìm kiếm các phần tử vỏ khối có độ chính xác cao để mô phỏng vỏ mỏng mà không lo ngại lỗi khóa số học.
  • Giảng viên, Học viên Cao học & Nghiên cứu sinh: Là tài liệu tham khảo giá trị phục vụ công tác giảng dạy chuyên đề PTHH nâng cao, luận văn thạc sĩ và luận án tiến sĩ kỹ thuật xây dựng công trình.

FAQ (Câu hỏi thường gặp) (250-300 từ)

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Phần tử vỏ khối CS4-S8 có khả năng khắc phục hoàn toàn hiện tượng khóa cắt ở dải chiều dày cực mỏng ($h/L = 10^{-8}$) và kháng méo lưới vượt trội trên bài toán tấm thoi Morley với sai số chỉ 0,11%.

2. Phương pháp của phần tử CSn-S8 có gì đặc biệt so với phần tử thông thường?

CSn-S8 chỉ dùng bậc tự do tịnh tiến, kết hợp kỹ thuật nội suy biến dạng cắt MITC4, khử khóa hình thang, tích hợp xấp xỉ biến dạng uốn giả định và làm trơn biến dạng màng bằng tích phân biên theo định lý phân kỳ Gauss.

3. Kết quả của nghiên cứu có thể tổng quát hóa (generalize) không?

Có. Mặc dù báo cáo khảo sát vật liệu đàn hồi đồng nhất tuyến tính, cấu trúc toán học của phần tử vỏ khối 8 nút này hoàn toàn sẵn sàng để mở rộng cho các bài toán phi tuyến hình học, vật liệu composite đa lớp hoặc vật liệu FGM.

4. Thời gian tính toán của CSn-S8 có lâu hơn phần tử vỏ phẳng không?

Thời gian tính toán của CS4-S8 chỉ chênh lệch nhẹ so với phần tử vỏ phẳng 4 nút và nhanh hơn đáng kể so với phần tử vỏ tam giác MITC3+ khi mô phỏng cùng một mức độ chính xác của lưới.

5. Tại sao nên dùng phần tử vỏ khối thay vì phần tử vỏ suy biến truyền thống?

Phần tử vỏ khối không có bậc tự do góc xoay, giúp dễ dàng ghép nối trực tiếp với các phần tử 3D solid khác mà không cần phương trình liên kết phụ, đồng thời mô phỏng trực tiếp trường ứng suất biến dạng theo chiều dày kết cấu.

Kết luận (150 từ)

Đề tài nghiên cứu khoa học cấp trường trọng điểm T2019-73TĐ do PGS. Châu Đình Thành chủ nhiệm đã phát triển thành công phần tử vỏ khối 8 nút làm trơn CSn-S8. Bằng việc kết hợp hài hòa giữa các kỹ thuật khử khóa tiên tiến, xấp xỉ biến dạng uốn giả định và làm trơn biến dạng màng, công thức đề xuất đã chứng minh tính ưu việt về độ chính xác, tốc độ hội tụ và khả năng khử khóa số học trên hàng loạt bài toán kết cấu tấm vỏ phức tạp.

Công trình đã công bố quốc tế thành công trên tạp chí Journal of Applied and Computational Mechanics (ESCI) năm 2020. Hướng phát triển tiếp theo của nghiên cứu tập trung vào bài toán phi tuyến hình học chuyển vị lớn, tấm vỏ phân lớp chức năng FGM và mô phỏng cơ học rạn nứt 3D.