Tổng quan luận án

Kết cấu dạng tấm và vỏ hai chiều (2D) là các cấu kiện chịu lực cơ bản, xuất hiện phổ biến trong các công trình dân dụng, cơ khí (mái vòm, tháp giải nhiệt, bể chứa, đường ống, bình chịu áp lực), công nghiệp đóng tàu (vỏ tàu ngầm) và hàng không vũ trụ (thân máy bay, cánh, thùng chứa). Do hình học phức tạp và sự đa dạng về điều kiện biên cũng như vật liệu cấu thành, các phương pháp giải tích giải nghiệm chính xác bị hạn chế lớn về phạm vi ứng dụng. Phương pháp phần tử hữu hạn (FEM) là công cụ tính toán số chủ đạo để mô phỏng và phân tích ứng xử cơ học của các kết cấu này. Tuy nhiên, khi áp dụng phần tử tứ giác bốn nút bậc thấp trong phân tích tấm/vỏ mỏng, FEM truyền thống bộc lộ những khiếm khuyết kỹ thuật điển hình:

  • Hiện tượng khóa cắt (shear locking) phát sinh khi tỷ lệ giữa chiều dày và kích thước kết cấu ($a/h$) trở nên rất nhỏ, khiến năng lượng biến dạng cắt bị thổi phồng giả tạo và làm độ võng tính toán nhỏ hơn thực tế.
  • Hiện tượng khóa màng (membrane locking) xuất hiện khi phân tích kết cấu vỏ có dạng hình học cong hoặc khi lưới phần tử bị méo mó, thô ráp, làm suy giảm nghiêm trọng độ chính xác trong các bài toán uốn nổi trội.
  • Sự bất liên tục của trường biến dạng và ứng suất qua các biên phần tử bậc thấp, gây khó khăn cho công tác đánh giá ứng suất chính xác ở khâu xử lý hậu kỳ.

Luận án tiến sĩ ngành Cơ kỹ thuật của tác giả Tôn Thất Hoàng Lân (năm 2022) tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, dưới sự hướng dẫn khoa học của PGS. Nguyễn Văn Hiếu và PGS. Châu Đình Thành, được thực hiện nhằm giải quyết các hạn chế số nêu trên.

Mục tiêu và nhiệm vụ nghiên cứu của đề tài được xác định cụ thể như sau:

  1. Nghiên cứu và phát triển tập hợp bốn phần tử tứ giác bốn nút phẳng cải tiến đơn giản trong thiết lập công thức toán học nhưng có khả năng triệt tiêu hoặc giảm thiểu tối đa hiện tượng khóa cắt và khóa màng.
  2. Xây dựng phần tử SQ4H kết hợp kỹ thuật làm trơn biến dạng trên miền con (CS-FEM) và lý thuyết biến dạng cắt bậc cao dạng $C^0$ ($C^0$-HSDT) phục vụ phân tích phi tuyến hình học của tấm phẳng và tấm gấp.
  3. Xây dựng phần tử SQ4T dựa trên kỹ thuật nội suy kép (Twice Interpolation Strategy - TIS) để nâng cao bậc xấp xỉ chuyển vị, khử tính bất liên tục của trường ứng suất và biến dạng qua biên phần tử trong phân tích tuyến tính và phi tuyến.
  4. Xây dựng phần tử SQ4C dựa trên kỹ thuật tổ hợp biến dạng (màng, uốn và cắt) kết hợp làm trơn biến dạng miền con để phân tích kết cấu tấm/vỏ có hoặc không có sườn gia cường trên lưới chia méo.
  5. Xây dựng phần tử SQ4P dựa trên chuỗi đa thức Chebyshev bậc cao nhằm mở rộng khả năng điều chỉnh bậc xấp xỉ số và độ chính xác của kết cấu tấm/vỏ theo cả lưới chia lẫn bậc đa thức.
  6. Đánh giá sai số, tính ổn định và thời gian tính toán của các phần tử đề xuất trên các bài toán chuẩn về uốn tĩnh, dao động tự do, mất ổn định và uốn phi tuyến.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Kết cấu dạng tấm và vỏ mỏng đến dày, tấm gấp, tấm/vỏ có sườn gia cường chế tạo từ vật liệu đẳng hướng, composite nhiều lớp (laminated composite), vật liệu phân lớp chức năng (P-FGM) và vật liệu xốp phân lớp chức năng gia cường tiểu cầu graphene (FGP-GPLs).
  • Phạm vi nghiên cứu: Phân tích ứng xử cơ học trong miền đàn hồi tuyến tính (uốn tĩnh, xác định tần số và dạng dao động tự do, xác định tải trọng tới hạn mất ổn định đàn hồi) và phân tích phi tuyến hình học (uốn phi tuyến của tấm phẳng, tấm gấp và tấm FGM).

Tổng quan tài liệu và vị trí của luận án

Luận án hệ thống hóa các mốc phát triển chính của kỹ thuật phần tử hữu hạn phân tích tấm/vỏ từ thập niên 1970 đến nay:

  • Các phần tử đẳng tham số $C^0$ ban đầu: Nhóm tác giả Irons và cộng sự, Zienkiewicz và cộng sự (1970-1980) đã thiết lập các phần tử dựa trên động học Reissner-Mindlin với trường chuyển vị và góc xoay độc lập. Nhược điểm chính là khóa cắt ở tấm mỏng. Kỹ thuật tích phân giảm (reduced integration) tuy giảm được độ cứng cắt nhưng lại sinh ra các mode năng lượng không (spurious zero energy modes / hourglass phenomena) trong phân tích động lực học.
  • Phương pháp nội suy ten-xơ hỗn hợp và sai lệch cắt rời rạc: Bathe và cộng sự đã giải quyết cơ bản khóa cắt với họ phần tử MITC (MITC4, MITC8, MITC9, MITC16). Nhóm tác giả Nguyen-Xuan, Nguyen-Thoi và cộng sự phát triển họ phần tử tam giác trơn dựa trên kỹ thuật sai lệch cắt rời rạc (CS-DSG, ES-DSG, NS-DSG). Tessler và cộng sự đề xuất các phần tử MIN (Mindlin plate elements) xử lý biến dạng cắt ngang.
  • Kỹ thuật khử khóa màng: Lee và cộng sự đề xuất phân chia miền tứ giác thành các tam giác con để tính toán biến dạng màng và chuyển về các điểm buộc biên (hình thành phần tử MITC4+), hỗ trợ tính toán kết cấu vỏ khi dùng lưới méo.
  • Kỹ thuật phần tử hữu hạn trơn (SFEM): Liu và cộng sự, Nguyen-Xuan và cộng sự, Nguyen-Thoi và cộng sự đưa ra SFEM trên nút (NS-FEM), trên cạnh (ES-FEM) và trên miền con (CS-FEM). Nguyen-Van và cộng sự phát triển thành công các phần tử tứ giác trơn MISQ20 và MISQ24 kết hợp MITC4 và trơn biến dạng màng/uốn cho kết cấu tấm và vỏ.
  • Kỹ thuật phần tử phổ (PSE) và đa thức Chebyshev: Zrahia và cộng sự, Sprague và cộng sự ứng dụng hàm Lagrangian qua các điểm Gauss-Legendre-Lobatto. Liu, He, Dang-Trung và cộng sự khai thác tính trực giao của đa thức Chebyshev trong xấp xỉ phần tử hữu hạn.
  • Kỹ thuật nội suy kép (TIS) và phân tích đẳng hình học (IGA): Kỹ thuật TIS được Bui, Wu, Zheng phát triển nhằm tạo tính liên tục cho trường ứng suất trong bài toán phẳng 2D. Phương pháp IGA do Hughes, Bazilevs, Gómez, Nguyen-Thanh, Thai-Hoang, Tran-Van phát triển dựa trên hàm cơ sở B-spline/NURBS tích hợp CAD/CAE.

Khoảng trống nghiên cứu luận án lựa chọn: Mặc dù các kỹ thuật tiên tiến trên đã mang lại nhiều cải tiến, các phần tử tứ giác 4 nút phẳng bậc thấp dùng cho vỏ vẫn gặp khó khăn trong việc cân bằng giữa cấu trúc đại số đơn giản, chi phí tính toán thấp và khả năng khắc phục đồng thời cả hiện tượng khóa cắt lẫn khóa màng trên lưới chia méo. Kỹ thuật TIS trước đây chủ yếu chỉ ứng dụng cho bài toán đàn hồi phẳng 2D; việc mở rộng TIS sang phân tích tấm/vỏ và phân tích phi tuyến chưa được thực hiện đầy đủ. Việc kết hợp lý thuyết biến dạng cắt bậc cao dạng $C^0$ với kỹ thuật trơn biến dạng miền con để phân tích tấm gấp phi tuyến, cũng như việc ứng dụng trực tiếp đa thức Chebyshev xây dựng phần tử tứ giác 4 nút bậc cao linh hoạt cho kết cấu FGP-GPLs, là những vấn đề cần được nghiên cứu và kiểm chứng hệ thống.

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

Khung lý thuyết vật liệu

Luận án thiết lập và áp dụng mô hình cơ học cho bốn nhóm vật liệu:

  1. Vật liệu đẳng hướng: Đặc trưng bởi mật độ khối lượng $\rho$, mô đun Young $E$, hệ số Poisson $\mu$, hệ số giãn nở nhiệt và mô đun đàn hồi cắt $G = E / [2(1+\mu)]$.
  2. Vật liệu composite nhiều lớp: Xét trạng thái ứng suất phẳng cho từng lớp trực hướng lệch trục. Quan hệ ứng suất - biến dạng trong hệ tọa độ quy chiếu $(x, y)$ được chuyển đổi từ hệ trục chính $(1, 2)$ qua góc hướng sợi $\theta$ bằng ma trận độ cứng biến đổi $Q'_{ij}$ ($i, j = 1, 2, 6$).
  3. Vật liệu phân lớp chức năng (P-FGM): Cấu tạo từ hỗn hợp gốm (ceramic) và kim loại (metal). Tỷ phần thể tích gốm $V_c(z)$ biến thiên liên tục theo chiều dày $z \in [-h/2, h/2]$ theo quy luật hàm mũ: $$V_c(z) = \left(\frac{2z + h}{2h}\right)^n, \quad V_m(z) = 1 - V_c(z)$$ với $n$ là chỉ số tỷ lệ thể tích ($n \ge 0$).
  4. Vật liệu xốp phân lớp chức năng gia cường graphene (FGP-GPLs): Mô hình hóa độ xốp theo ba dạng phân bố (P-S: đối xứng, P-A: bất đối xứng, P-U: phân bố đều) kết hợp ba dạng phân bố tiểu cầu graphene (GPL-S, GPL-A, GPL-U). Các đặc trưng cơ lý $E(z), G(z), \rho(z)$ phụ thuộc vào hệ số xốp $e_0 = 1 - E_2'/E_1'$, hàm phân bố khối lượng $e_m$, và tỷ phần khối lượng graphene $W_{F(GPL)}$.

Lý thuyết động học tấm/vỏ và công thức phần tử hữu hạn

  • Động học: Luận án áp dụng lý thuyết biến dạng cắt bậc nhất FSDT (Reissner-Mindlin) với 5 bậc tự do tại mỗi nút $(u, v, w, \beta_x, \beta_y)$ hoặc 6 bậc tự do khi bổ sung góc xoay khoan $\beta_z$ trong mặt phẳng; đồng thời sử dụng lý thuyết biến dạng cắt bậc cao HSDT được cải biên về dạng liên tục $C^0$ ($C^0$-HSDT) để triệt tiêu đòi hỏi về đạo hàm bậc hai của hàm dạng và không cần dùng hệ số hiệu chỉnh cắt $k_s$.
  • Phần tử vỏ phẳng: Xây dựng bằng cách kết hợp ma trận độ cứng của phần tử màng tứ giác đẳng tham số (chịu kéo/nén trong mặt phẳng) và phần tử tấm uốn tứ giác (chịu uốn và cắt ngoài mặt phẳng), sau đó chuyển từ hệ tọa độ địa phương sang hệ tọa độ tổng thể qua ma trận chuyển trục $T$.
  • Phương trình chủ đạo trong tính toán:
    • Uốn tĩnh: $K q = F$
    • Dao động tự do: $(K - \omega^2 M) q = 0$
    • Ổn định đàn hồi: $(K - P^* K_g) q = 0$
    • Uốn phi tuyến hình học: Áp dụng thuật toán lặp giải hệ phi tuyến $K_{NL}(q) q = F$.
Đại lượng / Ký hiệu Ý nghĩa cơ học trong luận án
$u, v, w$ Chuyển vị thẳng theo ba trục tọa độ $x, y, z$
$\beta_x, \beta_y, \beta_z$ Góc xoay quanh ba trục tọa độ $x, y, z$
$w^*$ Độ võng không thứ nguyên chuẩn hóa
$\sigma_x^, \tau^$ Ứng suất pháp và ứng suất tiếp chuẩn hóa
$\omega^*$ Tần số dao động tự do chuẩn hóa
$P^*$ Tải trọng tới hạn mất ổn định chuẩn hóa
$K, M, K_g$ Ma trận độ cứng đàn hồi, ma trận khối lượng, ma trận độ cứng hình học
$n_c, n_{sp}$ Số miền con trơn, số nút trong miền ảnh hưởng nội suy
$p, p_1, p_2$ Bậc của chuỗi đa thức Chebyshev

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 khái quát về mặt cơ học các dạng kết cấu thanh (1D), tấm/vỏ (2D) và khối (3D). Tác giả phân tích chi tiết ưu điểm và hạn chế của ba loại phần tử vỏ (vỏ khối, vỏ cong, vỏ phẳng), khẳng định tính ưu việt của phần tử vỏ phẳng 4 nút khi được bổ sung bậc tự do xoay trong mặt phẳng và kỹ thuật xử lý cắt ngang. Chương này phân tích nguồn gốc hiện tượng khóa cắt, khóa màng, tổng thuật các kỹ thuật phần tử hữu hạn hiện đại (MITC, SFEM, TIS, PSE, đa thức Chebyshev, IGA), từ đó xác lập mục tiêu, đối tượng, nhiệm vụ và cấu trúc 8 chương của luận án.

Chương 2: Cơ sở lý thuyết

Chương 2 trình bày chi tiết công thức toán học và cơ sở vật lý cho các loại vật liệu: vật liệu đẳng hướng, composite nhiều lớp lệch trục (với ma trận độ cứng quay trục $Q'_{ij}$), vật liệu FGM (với ba mô hình P-FGM, S-FGM, E-FGM, trong đó tập trung vào P-FGM), và vật liệu xốp FGP-GPLs với đầy đủ các hàm phân bố độ xốp và tiểu cầu graphene. Tiếp đó, chương thiết lập động học tấm theo FSDT và $C^0$-HSDT, xây dựng công thức phần tử màng tứ giác, phần tử tấm uốn, phần tử vỏ phẳng và hệ phương trình phần tử hữu hạn giải bài toán uốn tĩnh, dao động tự do và phân tích ổn định.

Chương 3: Phần tử SQ4H

Chương 3 trình bày việc xây dựng phần tử tứ giác 4 nút SQ4H dựa trên kỹ thuật làm trơn biến dạng trên các miền con ($n_c = 1, 2$) kết hợp lý thuyết biến dạng cắt bậc cao cải biên $C^0$-HSDT. Nhờ trường chuyển vị $C^0$, phần tử dễ dàng cài đặt vào mô hình số mà không đòi hỏi tính liên tục bậc cao của đạo hàm chuyển vị. Kết quả kiểm chứng số được thực hiện qua các bài toán phân tích uốn phi tuyến hình học:

  • Tấm vuông ngàm composite nhiều lớp chịu tải phân bố đều.
  • Tấm hình bình hành và hình thang ngàm composite nhiều lớp.
  • Tấm gấp $90^\circ$ và $60^\circ$ liên kết ngàm chịu tải phân bố đều vuông góc với bề mặt tấm. Kết quả cho thấy phần tử SQ4H triệt tiêu hiện tượng khóa cắt, cho độ hội tụ nhanh và giảm thiểu sự bất ổn định số trong phân tích phi tuyến hình học so với các phần tử truyền thống.

Chương 4: Phần tử SQ4T

Chương 4 phát triển phần tử tứ giác 4 nút SQ4T dựa trên chiến lược nội suy kép (TIS). Bằng cách mở rộng miền ảnh hưởng ra các nút lân cận ($n_{sp}$ nút) và tích hợp gradient trung bình của trường chuyển vị tại các nút, hàm dạng của phần tử đạt bậc cao hơn nhưng số bậc tự do tại từng nút không đổi. Kết quả số được kiểm chứng qua:

  • Bài toán dầm màng Cook: Đánh giá độ chính xác chuyển vị và tính liên tục của trường ứng suất qua biên phần tử so với phần tử MITC4.
  • Tấm vuông composite nhiều lớp $[0^\circ/90^\circ/90^\circ/0^\circ]$ và $[\theta^\circ/-\theta^\circ]$ chịu tải đều và tải hình sin: Khảo sát độ võng chuẩn hóa $w^*$ và sai số với tỷ số $a/h = 20$ và $100$.
  • Phân tích dao động tự do tấm composite đa lớp ($a/h = 5$) và vỏ trụ tựa đơn.
  • Phân tích mất ổn định tấm composite chịu nén đơn trục và vỏ trụ 5 lớp $[0^\circ/90^\circ/0^\circ/90^\circ/0^\circ]$.
  • Phân tích uốn phi tuyến hình học tấm hình bình hành FGM với điều kiện biên ngàm và tựa đơn.

Chương 5: Phần tử SQ4C

Chương 5 thiết lập phần tử tứ giác 4 nút SQ4C ứng dụng kỹ thuật tổ hợp biến dạng:

  1. Khử khóa cắt: Tính toán biến dạng cắt ngang thông qua 4 điểm buộc rời rạc trên các cạnh phần tử.
  2. Khử khóa màng: Xác định biến dạng màng thông qua việc chia phần tử thành các miền tam giác con với sự hỗ trợ của nút ảo số 5 ở tâm phần tử, sau đó chiếu về 4 điểm buộc màng trên biên.
  3. Làm trơn biến dạng: Áp dụng kỹ thuật trơn trên miền con (Cell-based smoothing) cho cả thành phần uốn và màng.

Phần tử SQ4C được áp dụng phân tích tĩnh, dao động tự do và ổn định cho:

  • Tấm vuông composite ngàm và tựa đơn chịu tải phân bố đều và hình sin.
  • Kết cấu vỏ kinh điển: Vỏ trụ có màng cứng ở hai đầu (Scordelis-Lo roof) trên cả hai loại lưới chia vuông đều và lưới méo; vỏ cầu chịu tải phân bố đều.
  • Kết cấu gia cường: Tấm vuông tựa đơn có một sườn ngang hoặc $n_{st}$ sườn ngang cách đều chịu uốn tĩnh, dao động tự do (khảo sát 6 dạng dao động đầu tiên) và chịu nén đơn trục mất ổn định.
Đặc tính kỹ thuật Phần tử SQ4H Phần tử SQ4T Phần tử SQ4C Phần tử SQ4P
Nền tảng lý thuyết $C^0$-HSDT & CS-FEM Twice Interpolation (TIS) Combined Strain & CS-FEM Đa thức Chebyshev ($p_1, p_2$)
Bậc tự do mỗi nút 5 hoặc 6 5 hoặc 6 5 hoặc 6 5 hoặc 6
Cơ chế khử khóa Trơn miền con $n_c=1,2$ TIS qua miền ảnh hưởng nút Điểm buộc cắt và màng Tăng bậc đa thức $p$
Khắc phục méo lưới Tốt trên tấm phẳng/gấp Rất tốt cho ứng suất biên Tối ưu cho vỏ cong trên lưới méo Hội tụ tốt theo bậc $p$
Ứng dụng chính Phi tuyến tấm phẳng, tấm gấp Tuyến tính/phi tuyến tấm, vỏ, FGM Tấm/vỏ có sườn gia cường Tấm, vỏ, vật liệu FGP-GPLs

Chương 6: Phần tử SQ4P

Chương 6 xây dựng phần tử tứ giác 4 nút SQ4P dựa trên chuỗi đa thức Chebyshev trực giao trên đoạn $[-1, 1]$. Hàm dạng 2D được xây dựng từ tích phân bố của các đa thức Chebyshev 1D bậc $p_1$ và $p_2$. Các khảo sát số bao gồm:

  • Khảo sát sự hội tụ của độ võng $w^$ và ứng suất $\sigma_x^$ tại tâm tấm vuông tựa đơn và ngàm khi thay đổi bậc $p_1, p_2$ từ bậc 3 đến bậc 6.
  • Khảo sát hiện tượng khóa cắt trên dải độ mảnh cực rộng ($a/h$ biến thiên từ 10 đến 100.000), chứng minh phần tử SQ4P hoàn toàn không bị khóa cắt ngay cả ở tỷ số $a/h = 10^5$.
  • Phân tích các kết cấu vỏ phức tạp: Vỏ bán cầu có khoét lỗ $18^\circ$ chịu tải tập trung và vỏ nửa yên ngựa (hyperbolic paraboloid shell).
  • Phân tích uốn tĩnh, tần số dao động tự do mode 1 và lực tới hạn mất ổn định đơn trục/đa trục của tấm FGP-GPLs theo các tổ hợp phân bố độ xốp (P-S, P-A) và phân bố graphene (GPL-S, GPL-A, GPL-U) khi biến thiên $e_0$ và $W_{F(GPL)}$.

Chương 7: Đánh giá sai số chung giữa các phần tử

Chương 7 thực hiện đối sánh trực tiếp hiệu năng tính toán giữa các phần tử đề xuất (SQ4H, SQ4T, SQ4C, SQ4P) và các phần tử trong y văn thông qua ba bài toán chuẩn:

  1. Tấm vuông đẳng hướng chịu tải phân bố đều: So sánh sai số độ võng tâm tấm với các tỷ số $a/h = 10, 100, 1000, 10000$ dưới điều kiện biên tựa đơn (SSSS) và ngàm (CCCC).
  2. Tấm vuông đẳng hướng dao động tự do: So sánh độ chính xác của bốn tần số dao động đầu tiên và thời gian tính toán (giây).
  3. Vỏ cầu đẳng hướng tựa đơn chịu tải phân bố đều: So sánh độ võng tâm và mức độ hội tụ lưới.

Kết quả đánh giá tổng hợp cho thấy phần tử SQ4C đạt hiệu năng cao và ổn định nhất đối với các kết cấu vỏ cong có lưới chia méo; phần tử SQ4T vượt trội về độ mịn của trường ứng suất; phần tử SQ4P cho tốc độ hội tụ nhanh khi tăng bậc đa thức mà không cần làm mịn lưới quá mức; phần tử SQ4H ổn định trong phân tích phi tuyến hình học của tấm và tấm gấp.

Chương 8: Kết luận và hướng phát triển

Chương 8 tổng kết toàn bộ các kết quả đạt được của luận án, khẳng định các mục tiêu đề ra đã được hoàn thành, đồng thời chỉ rõ các hướng mở rộng nghiên cứu trong tương lai.

Kết quả và những đóng góp mới

Luận án tự xác định các đóng góp khoa học chính bao gồm:

  1. Phát triển phần tử SQ4H: Xây dựng thành công phần tử tứ giác 4 nút dựa trên kỹ thuật trơn biến dạng trên miền con kết hợp với lý thuyết biến dạng cắt bậc cao cải biên dạng $C^0$ ($C^0$-HSDT). Phần tử này cải thiện độ chính xác và giảm thiểu sự bất ổn định số trong phân tích uốn phi tuyến hình học của kết cấu tấm phẳng và tấm gấp composite nhiều lớp.
  2. Phát triển phần tử SQ4T: Thiết lập phần tử tứ giác 4 nút dựa trên chiến lược nội suy kép (TIS). Bằng cách thiết lập các hàm dạng bậc cao kể đến ảnh hưởng của tập hợp các nút lân cận và gradient trung bình nút, phần tử cải thiện đáng kể sự bất liên tục của trường biến dạng và trường ứng suất qua các biên phần tử.
  3. Phát triển phần tử SQ4C: Xây dựng phần tử tứ giác 4 nút dựa trên kỹ thuật tổ hợp biến dạng (màng, uốn và cắt) kết hợp trơn biến dạng miền con. Phần tử này loại bỏ đồng thời hiện tượng khóa cắt và khóa màng, đảm bảo độ chính xác cao khi phân tích kết cấu tấm/vỏ có hoặc không có sườn gia cường ngay cả khi sử dụng lưới phần tử bị méo mó.
  4. Phát triển phần tử SQ4P: Đề xuất phần tử tứ giác 4 nút dựa trên chuỗi đa thức Chebyshev. Phần tử cho phép nâng cao độ chính xác số linh hoạt thông qua việc thay đổi mật độ lưới chia hoặc bậc của đa thức Chebyshev ($p_1, p_2$), hoàn toàn không bị khóa cắt ở dải tỷ số độ mảnh cực lớn ($a/h = 10^5$).
  5. Ứng dụng và kiến nghị giải pháp: Cung cấp công cụ tính toán số hiệu quả để phân tích các dạng vật liệu phức tạp (composite nhiều lớp, P-FGM, FGP-GPLs) và kết cấu có sườn. Các công thức phần tử đề xuất có thể tích hợp trực tiếp vào các mô-đun giải bài toán của các phần mềm phần tử hữu hạn thương mại hoặc mã nguồn mở phục vụ tính toán công trình.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án:

  • Nghiên cứu tập trung vào mô hình ứng xử vật liệu đàn hồi tuyến tính; chưa xét đến ứng xử phi tuyến vật liệu (chảy dẻo, từ biến, tổn thương tích lũy hoặc nứt gãy).
  • Phân tích kết cấu vỏ được thực hiện trên cơ sở phần tử vỏ phẳng (flat shell element); chưa mở rộng trực tiếp sang các phần tử vỏ cong liên tục bậc cao hoặc phần tử vỏ khối 3D.
  • Chưa khảo sát các bài toán phân tích động lực học phi tuyến theo lịch sử thời gian (time-history dynamic analysis) và bài toán tương tác đa trường vật lý phức tạp (nhiệt - cơ - điện).

Hướng nghiên cứu tiếp theo:

  • Mở rộng các kỹ thuật phần tử SQ4H, SQ4T, SQ4C, SQ4P cho phân tích phi tuyến vật liệu và phân tích phá hủy kết cấu tấm/vỏ.
  • Phát triển kỹ thuật tính toán cho các bài toán phân tích động lực học phi tuyến và tương tác đa môi trường vật lý.
  • Nghiên cứu kết hợp các kỹ thuật đề xuất với thuật toán tinh chỉnh lưới thích ứng (adaptive mesh refinement) để tự động tối ưu hóa lưới chia trong các khu vực tập trung ứng suất cao.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị chuyên môn cao cho:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Cơ kỹ thuật, Cơ học tính toán, Xây dựng công trình và Kỹ thuật Cơ khí: Tham khảo cơ sở lý thuyết chi tiết về các kỹ thuật phần tử hữu hạn hiện đại (CS-FEM, TIS, Combined strain, Chebyshev polynomial FEM), các kỹ thuật khử khóa số (locking-free techniques), và các mô hình cơ học vật liệu mới (FGM, FGP-GPLs).
  • Các nhà nghiên cứu phát triển phần mềm mô phỏng (CAE): Cung cấp thuật toán chi tiết và các ví dụ kiểm chứng số (benchmark tests) để lập trình, đóng gói thành các mô-đun phần tử vỏ phẳng 4 nút có độ chính xác cao và khả năng chống méo lưới tốt.
  • Kỹ sư thiết kế kết cấu công trình dân dụng, cơ khí, tàu thủy và hàng không: Tham khảo phương pháp đánh giá ứng xử uốn, tần số dao động tự do, tải trọng mất ổn định và đáp ứng phi tuyến của tấm/vỏ có sườn gia cường hoặc làm từ vật liệu composite tiên tiến.

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

1. Bốn phần tử tứ giác đề xuất trong luận án có tên là gì và dựa trên những kỹ thuật cốt lõi nào? Bốn phần tử gồm: SQ4H (dựa trên kỹ thuật trơn biến dạng miền con kết hợp lý thuyết cắt bậc cao $C^0$-HSDT); SQ4T (dựa trên chiến lược nội suy kép TIS); SQ4C (dựa trên kỹ thuật tổ hợp biến dạng màng, uốn, cắt và trơn biến dạng miền con); và SQ4P (dựa trên chuỗi đa thức Chebyshev bậc cao).

2. Phần tử SQ4C khử đồng thời hiện tượng khóa cắt và khóa màng bằng phương pháp nào? SQ4C tính toán biến dạng cắt ngang qua 4 điểm buộc biên để khử khóa cắt; tính toán biến dạng màng thông qua việc chia phần tử thành các tam giác con kết hợp nút ảo trung tâm rồi đưa về 4 điểm buộc màng trên biên để khử khóa màng; đồng thời áp dụng kỹ thuật làm trơn biến dạng trên các miền con cho cả thành phần màng và uốn.

3. Ưu điểm nổi bật nhất của kỹ thuật nội suy kép (TIS) trong phần tử SQ4T là gì? Kỹ thuật TIS trong SQ4T thiết lập các hàm dạng bậc cao dựa trên giá trị nút và gradient trung bình nút trong miền ảnh hưởng lân cận, giúp khắc phục nhược điểm bất liên tục của trường ứng suất và biến dạng qua các biên phần tử bậc thấp mà không làm tăng số bậc tự do tại các nút của phần tử.

4. Luận án đã khảo sát những mô hình vật liệu tiên tiến nào cho kết cấu tấm và vỏ? Luận án khảo sát vật liệu đẳng hướng, vật liệu composite nhiều lớp lệch trục, vật liệu phân lớp chức năng P-FGM (với tỷ phần thể tích gốm/kim loại biến thiên theo quy luật hàm mũ qua chiều dày), và vật liệu xốp phân lớp chức năng gia cường tiểu cầu graphene FGP-GPLs (với ba kiểu phân bố độ xốp P-S, P-A, P-U và ba kiểu phân bố graphene GPL-S, GPL-A, GPL-U).

5. Phần tử SQ4P thể hiện khả năng chống khóa cắt như thế nào khi chiều dày tấm cực mỏng? Kết quả kiểm chứng số của phần tử SQ4P cho thấy khi thay đổi tỷ số chiều dài trên chiều dày $a/h$ từ 10 đến 100.000, nghiệm độ võng chuẩn hóa và ứng suất chuẩn hóa vẫn hội tụ chính xác và ổn định về nghiệm giải tích, chứng minh phần tử hoàn toàn không bị hiện tượng khóa cắt ngay cả ở tỷ số độ mảnh $a/h = 10^5$.

Kết luận

Luận án tiến sĩ của tác giả Tôn Thất Hoàng Lân đã phát triển và hoàn thiện thành công bốn kỹ thuật phần tử hữu hạn tứ giác bốn nút phẳng cải tiến (SQ4H, SQ4T, SQ4C, SQ4P) dùng cho phân tích cơ học kết cấu tấm và vỏ. Các phần tử đề xuất đã giải quyết hiệu quả các hạn chế kinh điển của phương pháp phần tử hữu hạn bậc thấp như hiện tượng khóa cắt, khóa màng, sự bất liên tục của trường ứng suất qua biên phần tử và sự nhạy cảm với lưới chia méo. Hệ thống ví dụ số phong phú trên các bài toán uốn tĩnh, dao động tự do, mất ổn định đàn hồi và uốn phi tuyến hình học đã chứng minh tính chính xác, độ tin cậy và khả năng ứng dụng thực tiễn của các mô hình số đối với nhiều dạng kết cấu và vật liệu kỹ thuật phức tạp.