Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Trong cơ học kết cấu và kỹ thuật công trình, kết cấu dạng tấm và vỏ (kết cấu hai chiều - 2D có chiều dày nhỏ so với các kích thước mặt bằng) đóng vai trò bộ phận chịu lực chính trong nhiều lĩnh vực: xây dựng dân dụng và công nghiệp (sàn, tường, mái vòm, silo, bể chứa), cơ khí (tháp giải nhiệt, đường ống dẫn, bình chịu áp lực), kỹ thuật tàu thủy (vỏ tàu ngầm) và hàng không - vũ trụ (thân máy bay, cánh, vỏ thùng nhiên liệu). Dưới tác động của tải trọng cơ học và môi trường làm việc khắc nghiệt, việc tính toán chính xác các trạng thái ứng xử uốn tĩnh, dao động tự do và ổn định của tấm/vỏ là yêu cầu kỹ thuật bắt buộc.

Khi áp dụng phương pháp giải tích, việc giải quyết các kết cấu có hình dạng phức tạp, điều kiện biên hỗn hợp hoặc cấu tạo từ vật liệu không đồng nhất gặp nhiều giới hạn toán học. Phương pháp phần tử hữu hạn (FEM) là công cụ tính toán số được sử dụng phổ biến nhất để giải quyết các bài toán này. Tuy nhiên, FEM truyền thống dành cho tấm và vỏ vẫn bộc lộ những hạn chế kỹ thuật:

  1. Hiện tượng khóa cắt (shear locking) xuất hiện khi phân tích kết cấu tấm/vỏ rất mỏng theo lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin, khiến độ cứng mô hình bị khuếch đại giả tạo và độ võng suy giảm nghiêm trọng.
  2. Hiện tượng khóa màng (membrane locking) xảy ra khi sử dụng lưới phần tử thô hoặc lưới méo để rời rạc hóa kết cấu vỏ cong trong các bài toán mà ứng xử uốn giữ vai trò chủ đạo.
  3. Kỹ thuật tích phân giảm (reduced integration) nhằm triệt tiêu hiện tượng khóa thường kéo theo sự xuất hiện của các dạng năng lượng không giả (spurious zero-energy modes hay hiện tượng hourglass), gây mất ổn định kết quả tính dao động.
  4. Hiện tượng gián đoạn ứng suất và biến dạng qua ranh giới giữa các phần tử lân cận, đòi hỏi công đoạn xử lý hậu kỳ phức tạp.

Do đó, việc nghiên cứu phát triển các kỹ thuật phần tử hữu hạn cải tiến, xây dựng các phần tử tứ giác 4 nút đơn giản nhưng khắc phục được các hiện tượng khóa và nâng cao độ chính xác số là một nhu cầu học thuật và thực tiễn liên tục.

Mục tiêu và nhiệm vụ nghiên cứu

Nhiệm vụ tổng quát của luận án là: "Phát triển các kỹ thuật phần tử hữu hạn dùng để phân tích kết cấu dạng tấm/vỏ."

Các mục tiêu nghiên cứu cụ thể bao gồm:

  1. 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 miền con kết hợp lý thuyết biến dạng cắt bậc cao dạng cải biên $C^0$-HSDT phục vụ phân tích phi tuyến hình học kết cấu tấm phẳng và tấm gấp.
  2. Xây dựng phần tử tứ giác 4 nút SQ4T dựa trên chiến lược nội suy kép (Twice Interpolation Strategy - TIS) nhằm nâng cao bậc xấp xỉ trường chuyển vị và khắc phục tính gián đoạn của trường ứng suất qua biên phần tử trong phân tích tuyến tính và phi tuyến.
  3. Xây dựng phần tử tứ giác 4 nút SQ4C dựa trên kỹ thuật tổ hợp biến dạng (màng, uốn, cắt) kết hợp 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, loại bỏ đồng thời khóa cắt và khóa màng trên lưới chia méo.
  4. Xây dựng phần tử tứ giác 4 nút SQ4P dựa trên chuỗi đa thức trực giao Chebyshev để phân tích kết cấu tấm/vỏ cấu tạo từ vật liệu xốp phân lớp chức năng gia cường hạt tiểu cầu graphene (FGP-GPLs).

Đố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ỏ 2D (tấm phẳng, tấm gấp, vỏ trụ, vỏ cầu, vỏ bán cầu khoét lỗ, vỏ yên ngựa, tấm/vỏ có sườn gia cường).
  • Vật liệu khảo sát: Vật liệu đẳng hướng đàn hồi; vật liệu composite nhiều lớp cốt sợi (lệch trục, trực hướng, đối xứng/bất đối xứng); vật liệu phân lớp chức năng biến thiên theo quy luật hàm mũ (P-FGM); 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 phân tích cơ học: Phân tích uốn tĩnh (tuyến tính và phi tuyến hình học), phân tích dao động tự do (tần số dao động tự nhiên và các dạng dao động riêng), phân tích ổn định đàn hồi (tải trọng tới hạn đơn trục và đa trục).

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

Các hướng nghiên cứu tiền nhiệm

Trong quá trình tổng quan lịch sử phát triển của các phần tử hữu hạn phân tích tấm và vỏ, luận án đã hệ thống hóa các nhóm công trình tiêu biểu:

  1. Phần tử đẳng tham số $C^0$ và kỹ thuật tích phân giảm: Vào thập niên 1970–1980, Irons và cộng sự cùng Zienkiewicz và cộng sự đã đề xuất phần tử đẳng tham số $C^0$ nội suy độc lập trường chuyển vị và góc xoay dựa trên lý thuyết Reissner-Mindlin. Giới hạn chính của phần tử này là hiện tượng khóa cắt khi chiều dày tấm giảm. Kỹ thuật tích phân giảm được đề xuất nhằm giải tỏa năng lượng biến dạng cắt dư thừa nhưng lại phát sinh hiện tượng mất ổn định dạng đồng hồ cát (hourglass phenomena).

  2. Phương pháp nội suy ten-xơ hỗn hợp (MITC) và phương pháp khe trượt rời rạc (DSG, MIN): Để giải quyết triệt để khóa cắt, Bathe và cộng sự đã xây dựng thành công kỹ thuật MITC thông qua các phần tử tứ giác MITC4, phần tử 8 nút MITC8, 9 nút MITC9 và 16 nút MITC16. Lee và Bathe phát triển biến thể MITC4+ để hạn chế khóa màng. Song song đó, Bletzinger và cộng sự đưa ra phương pháp khe trượt rời rạc DSG (như DSG3), trong khi Tessler và Hughes phát triển họ phần tử MIN (Mindlin) xấp xỉ biến dạng cắt từ các điểm rời rạc trên biên.

  3. Phương pháp phần tử hữu hạn trơn (Smoothed Finite Element Methods - SFEM): Liu và cộng sự, Nguyen-Xuan và cộng sự, Nguyen-Thoi và cộng sự đã phát triển hệ thống các phương pháp làm trơn biến dạng: trơn trên nút (NS-FEM), trơn trên cạnh (ES-FEM) và trơn trên miền con (CS-FEM). Các phần tử tam giác và tứ giác kết hợp như CS-DSG3, ES-DSG3, CSMIN3 chứng minh tốc độ hội tụ cao và độ cứng xấp xỉ mềm hơn FEM truyền thống. Đặc biệt, Nguyen-Van và cộng sự đã kết hợp MITC4 với kỹ thuật trơn biến dạng miền con để tạo nên hai phần tử tứ giác bậc thấp hiệu quả là MISQ20 và MISQ24.

  4. Kỹ thuật nội suy kép (Twice Interpolation Strategy - TIS): Bui và cộng sự, Wu và cộng sự, Zheng và cộng sự đã nghiên cứu chiến lược TIS cho các bài toán cơ học chất rắn 2D phẳng, sử dụng hàm nội suy bậc cao dựa trên giá trị nút và gradient trung bình của các nút lân cận, tạo ra trường ứng suất liên tục qua biên phần tử.

  5. Phương pháp phần tử hữu hạn mở rộng, IGA và đa thức Chebyshev: Hughes và cộng sự khởi xướng phương pháp phân tích đẳng hình học (IGA) sử dụng hàm cơ sở B-spline/NURBS, sau đó được phát triển bởi Bazilevs, Gómez, Nguyen-Thanh, Thai-Hoang, Tran-Van. Zrahia và cộng sự, Sprague và cộng sự khảo sát phần tử phổ tấm (PSE) với các điểm Gauss-Legendre-Lobatto. Ngoài ra, việc ứng dụng đa thức Chebyshev trực giao vào thiết lập hàm dạng xấp xỉ được khảo sát bởi Liu và cộng sự, He và cộng sự, Dang-Trung và cộng sự.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

Mặc dù các kỹ thuật trên đã đạt nhiều kết quả, việc thiết lập các phần tử tứ giác 4 nút phẳng có cấu trúc thuật toán đơn giản nhưng đồng thời xử lý được khóa cắt, triệt tiêu khóa màng khi lưới chia méo, khử tính gián đoạn ứng suất liên phần tử và ứng dụng được cho các bài toán phi tuyến hình học của vật liệu composite nhiều lớp, FGM và FGP-GPLs vẫn là một khoảng trống học thuật cần tiếp tục hoàn thiện.

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

Các lý thuyết biến dạng cơ học của tấm và vỏ

Luận án sử dụng hai lý thuyết biến dạng cắt chính:

  1. Lý thuyết biến dạng cắt bậc nhất (FSDT): Dựa trên giả thiết Reissner-Mindlin, đoạn thẳng vuông góc với mặt trung hòa trước biến dạng vẫn thẳng nhưng không nhất thiết vuông góc với mặt trung hòa sau biến dạng: $$u(x,y,z) = u_0(x,y) + z\beta_x(x,y)$$ $$v(x,y,z) = v_0(x,y) + z\beta_y(x,y)$$ $$w(x,y,z) = w_0(x,y)$$ Trong đó $u_0, v_0, w_0$ là các chuyển vị tại mặt trung hòa; $\beta_x, \beta_y$ là các góc xoay quanh trục $y$ và $x$. Hệ số hiệu chỉnh biến dạng cắt ngang được lấy bằng $k_s = 5/6$.

  2. Lý thuyết biến dạng cắt bậc cao ($C^0$-HSDT): Trường chuyển vị uốn phi tuyến theo tọa độ chiều dày $z$ (TSDT - Reddy) triệt tiêu ứng suất cắt ngang trên hai bề mặt tự do $z = \pm h/2$, không cần sử dụng hệ số hiệu chỉnh cắt $k_s$. Bằng cách phân rã biến dạng góc xoay thành các biến độc lập, kỹ thuật $C^0$-HSDT loại bỏ yêu cầu tính liên tục $C^1$ của đạo hàm chuyển vị võng, cho phép sử dụng các hàm dạng $C^0$ chuẩn của phần tử hữu hạn 4 nút.

  3. Mô hình phần tử vỏ phẳng (Flat Shell Element): Phần tử vỏ phẳng được thiết lập thông qua sự kết hợp giữa trường chuyển vị màng trong mặt phẳng (bậc tự do $u, v, \theta_z$) và trường chuyển vị uốn ngoài mặt phẳng (bậc tự do $w, \theta_x, \theta_y$), sử dụng ma trận chuyển tọa độ cục bộ sang tọa độ tổng thể.

Mô hình hóa các loại vật liệu

  1. Vật liệu đẳng hướng: Đặc trưng bởi mô đun Young ($E$), hệ số Poisson ($\mu$), mật độ khối lượng ($\rho$), mô đun đàn hồi cắt ($G = \frac{E}{2(1+\mu)}$) và hệ số giãn nở nhiệt.
  2. Vật liệu composite phân lớp: Cấu tạo từ các lớp trực hướng xếp chồng với góc định hướng sợi $\theta$. Ma trận độ cứng đàn hồi của từng lớp được biến đổi từ hệ trục chính $(1,2)$ sang hệ trục tọa độ kết cấu $(x,y)$ thông qua ma trận quay ten-xơ ứng suất và biến dạng.
  3. Vật liệu phân lớp chức năng (P-FGM): Tỷ phần thể tích gốm ($V_c$) và kim loại ($V_m$) biến thiên liên tục theo chiều dày $z$ theo quy luật hàm lũy thừa: $$V_c(z) = \left(\frac{2z+h}{2h}\right)^n, \quad V_m(z) = 1 - V_c(z)$$ Đặc tính đàn hồi hiệu dụng $P(z)$ (mô đun $E$, mật độ $\rho$) được tính toán qua: $$P(z) = (P_c - P_m)V_c(z) + P_m$$
  4. Vật liệu xốp phân lớp chức năng gia cường tiểu cầu graphene (FGP-GPLs):
    • Các dạng phân bố độ xốp theo chiều dày: Đối xứng (P-S), Bất đối xứng (P-A), và Đồng đều (P-U). Hệ số độ xốp danh định $e_0 = 1 - E_2'/E_1'$.
    • Các dạng phân bố tiểu cầu graphene: Đối xứng (GPL-S), Bất đối xứng (GPL-A), và Đồng đều (GPL-U).
    • Đặc tính vật liệu hiệu dụng được xác định dựa trên mô hình cơ học vi mô Halpin-Tsai và phân phối chuẩn Gauss.

Phương pháp phân tích và công thức biến phân

Hệ phương trình phần tử hữu hạn được thiết lập từ nguyên lý thế năng toàn phần dừng:

  • Phân tích uốn tĩnh tuyến tính: $[K]{q} = {F}$
  • Phân tích uốn phi tuyến hình học: Thiết lập theo công thức von Kármán, giải hệ phương trình phi tuyến $[K_{NL}({q})]{q} = {F}$ bằng thuật toán lặp Newton-Raphson.
  • Phân tích dao động tự do: $([K] - \omega^2[M]){q} = {0}$, xác định tần số dao động tự nhiên $\omega$ và dạng dao động.
  • Phân tích ổn định: $([K] - \lambda[K_g]){q} = {0}$, xác định hệ số tải trọng tới hạn $\lambda$ và dạng mất ổn định.

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

Chương 1: Tổng quan

Chương 1 trình bày khái quát về hình học và trạng thái ứng xử cơ học của kết cấu tấm/vỏ; phân tích đặc điểm làm việc của phần tử màng, phần tử tấm uốn và phần tử vỏ phẳng. Chương đánh giá chi tiết nguyên nhân phát sinh hiện tượng khóa cắt, khóa màng, hiện tượng đồng hồ cát và lược khảo các giải pháp khắc phục trong FEM (MITC, DSG, MIN, SFEM, TIS, IGA, PSE, đa thức Chebyshev). Từ đó, chương xác định mục tiêu xây dựng 4 phần tử tứ giác 4 nút mới (SQ4H, SQ4T, SQ4C, SQ4P) cùng cấu trúc tổng thể của luận án.

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

Chương 2 thiết lập hệ thống phương trình cấu nạp và hình học cho 4 nhóm vật liệu: vật liệu đẳng hướng, composite nhiều lớp, vật liệu phân lớp chức năng (P-FGM, S-FGM, E-FGM) và vật liệu xốp FGP-GPLs. Trình bày chi tiết cơ sở toán học của lý thuyết FSDT, lý thuyết HSDT và kỹ thuật làm trơn $C^0$-HSDT. Thiết lập công thức phần tử màng tứ giác đẳng tham số, phần tử tấm uốn có kể biến dạng cắt, phần tử vỏ phẳng và các ma trận hệ thống ($[K], [M], [K_g]$).

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 sự kết hợp giữa lý thuyết biến dạng cắt bậc cao dạng $C^0$-HSDT và kỹ thuật làm trơn biến dạng trên miền con ($n_c = 1, 2$).

  • Luận điểm và giải pháp: Việc kết hợp $C^0$-HSDT loại bỏ hệ số hiệu chỉnh cắt $k_s$, trong khi kỹ thuật trơn biến dạng trên miền con làm giảm độ cứng giả tạo của mô hình, giúp quá trình phân tích phi tuyến hình học hội tụ ổn định.
  • Kiểm chứng số: Phân tích uốn phi tuyến cho tấm vuông nhiều lớp ngàm chịu tải đều, tấm hình bình hành ngàm, tấm hình thang ngàm, và kết cấu tấm gấp góc $90^\circ$ và $60^\circ$ chịu tải phân bố đều vuông góc với bề mặt. Kết quả cho thấy đường cong tải trọng - độ võng của phần tử SQ4H phản ánh chính xác ứng xử phi tuyến mà không bị suy biến số.

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

Chương 4 thiết lập phần tử tứ giác 4 nút SQ4T dựa trên chiến lược nội suy kép (TIS - Twice Interpolation Strategy).

  • Luận điểm và giải pháp: Hàm dạng bậc cao được xây dựng thông qua tập hợp các nút lân cận ($n_{sp}$ nút) nằm trong miền ảnh hưởng của phần tử, kết hợp giá trị nút và gradient trung bình của nút.
  • Kiểm chứng số:
    1. Bài toán màng Cook (Cook's membrane): đánh giá chuyển vị đỉnh $w_A$ và biểu đồ ứng suất $\sigma_x, \tau_{xy}$; phần tử SQ4T cho trường ứng suất liên tục qua biên phần tử tốt hơn so với MITC4.
    2. Phân tích uốn tĩnh tấm vuông composite nhiều lớp chịu tải đều và tải hình sin.
    3. Phân tích dao động tự do tấm và vỏ trụ tựa đơn.
    4. Phân tích ổn định tấm vuông $[0^\circ/90^\circ/90^\circ/0^\circ]$ chịu nén đơn trục và vỏ trụ 5 lớp.
    5. Phân tích uốn phi tuyến tấm hình bình hành làm bằng vật liệu P-FGM với các liên kết ngàm (CCCC) và tựa đơn (SSSS).

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

Chương 5 phát triển phần tử tứ giác 4 nút SQ4C dựa trên kỹ thuật tổ hợp biến dạng (Combined strain strategy).

  • Luận điểm và giải pháp: Tích hợp 4 điểm buộc (tying points) trên các cạnh phần tử để nội suy biến dạng cắt ngang (khử khóa cắt), kết hợp kỹ thuật chia miền con tam giác với 4 điểm buộc để tính toán biến dạng màng (khử khóa màng) và làm trơn biến dạng trên các miền con tứ giác. Mô hình hóa phần tử tấm/vỏ có sườn gia cường ($n_{st}$ sườn) bằng cách tích hợp trực tiếp độ cứng của sườn vào ma trận độ cứng phần tử vỏ.
  • Kiểm chứng số: Phân tích vỏ trụ có màng cứng ở hai đầu (Scordelis-Lo roof), vỏ bán cầu khoét lỗ $18^\circ$, vỏ cầu chịu tải phân bố đều; phân tích tấm vuông tựa đơn có gia cường sườn ngang chịu nén đơn trục. Phần tử SQ4C duy trì độ chính xác cao trên cả lưới vuông chuẩn và lưới chia méo (distorted mesh), triệt tiêu hiện tượng khóa màng.

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 hàm xấp xỉ chuỗi đa thức trực giao Chebyshev 1D/2D với bậc đa thức tùy chọn ($p_1, p_2$).

  • Luận điểm và giải pháp: Tận dụng tính chất trực giao và hội tụ nhanh của đa thức Chebyshev để nâng cao bậc xấp xỉ trường chuyển vị trong phần tử mà không cần tăng thêm số nút hình học ngoài 4 nút góc.
  • Kiểm chứng số:
    1. Khảo sát sự hội tụ của độ võng chuẩn hóa $w^$ và ứng suất chuẩn hóa $\sigma_x^$ khi thay đổi bậc đa thức ($p_1 = p_2 = 3, 4, 5, 6$) và tỷ số độ mảnh $a/h$ từ $10$ đến $100000$.
    2. Kiểm tra độ nhạy trên bài toán vỏ bán cầu khoét lỗ $18^\circ$ và vỏ nửa yên ngựa (hyperbolic paraboloid).
    3. Phân tích uốn tĩnh, tần số dao động tự do mode 1 và tải trọng tới hạn mất ổn định đơn trục/đa trục của tấm vuông làm bằng vật liệu FGP-GPLs dưới ảnh hưởng của hệ số độ xốp $e_0$ và tỷ lệ khối lượng graphene $W_{GPL}$ ($1.0%$).

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

Chương 7 tổng hợp và so sánh hiệu năng tính toán của cả 4 phần tử đề xuất (SQ4H, SQ4T, SQ4C, SQ4P) trên cùng một tập hợp bài toán chuẩn:

  1. Tấm vuông đẳng hướng chịu tải đều liên kết tựa đơn (SSSS) và ngàm (CCCC) với các tỷ số độ mảnh $a/h = 10, 100, 1000, 10000$.
  2. Tấm vuông đẳng hướng dao động tự do: so sánh 4 tần số dao động đầu tiên.
  3. Vỏ cầu đẳng hướng tựa đơn chịu tải phân bố đều.
  4. Đánh giá thời gian tiêu tốn CPU (tính bằng giây) giữa các phần tử, phân tích tương quan giữa chi phí tính toán và độ chính xác hội tụ.

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 về mặt phương pháp luận, công thức phần tử và mô phỏng số; đồng thời chỉ ra các hạn chế hiện tại và định hướng các nghiên cứu tiếp theo.

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

Các đóng góp mới của luận án

  1. 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 lý thuyết biến dạng cắt bậc cao dạng $C^0$-HSDT, cải thiện độ chính xác và giảm thiểu sự mất ổn định số trong phân tích phi tuyến hình học kết cấu tấm phẳng và tấm gấp.
  2. Phần tử SQ4T: Đề xuất phần tử tứ giác 4 nút ứng dụng kỹ thuật nội suy kép (TIS), thiết lập các hàm dạng bậc cao kể đến ảnh hưởng của nhóm nút lân cận và gradient trung bình nút, khắc phục hiện tượng gián đoạn ứng suất và biến dạng qua biên phần tử trong phân tích tĩnh và phi tuyến.
  3. Phần tử SQ4C: Phát triển phần tử tứ giác 4 nút dựa trên chiến lược tổ hợp biến dạng (màng, uốn, cắt) và làm trơn miền con, loại bỏ hiện tượng khóa màng và khóa cắt khi phân tích kết cấu tấm/vỏ có hoặc không có sườn gia cường, đặc biệt hiệu quả khi sử dụng lưới chia méo.
  4. Phần tử SQ4P: Ứng dụng chuỗi đa thức Chebyshev để xây dựng phần tử tứ giác 4 nút có khả năng nâng cao độ chính xác linh hoạt theo cả việc chia mịn lưới ($h$-refinement) lẫn tăng bậc đa thức ($p$-refinement), áp dụng hiệu quả cho bài toán kết cấu tấm/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).

Bảng tổng hợp đặc tính kỹ thuật của các phần tử đề xuất

Tên phần tử Kỹ thuật nền tảng Bậc tự do / Nút Khắc phục hiện tượng khóa Ứng dụng phân tích tiêu biểu
SQ4H Trơn biến dạng miền con ($n_c=1,2$) + $C^0$-HSDT $C^0$ cải biên Khóa cắt (không cần $k_s$) Uốn phi tuyến hình học tấm phẳng, tấm gấp
SQ4T Kỹ thuật nội suy kép (TIS) + Miền ảnh hưởng ($n_{sp}$) 5–6 DOFs Gián đoạn ứng suất qua biên Uốn tĩnh, dao động, ổn định, phi tuyến tấm/vỏ
SQ4C Tổ hợp biến dạng + Điểm buộc + Trơn miền con 6 DOFs Khóa cắt và Khóa màng (lưới méo) Tấm/vỏ có hoặc không có sườn gia cường
SQ4P Chuỗi đa thức Chebyshev ($p_1, p_2$) Biến thiên theo bậc $p$ Khóa cắt trên tấm cực mỏng Tấm/vỏ vật liệu phức hợp FGP-GPLs

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

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

  • Phân tích kết cấu vỏ chủ yếu dựa trên mô hình phần tử vỏ phẳng (flat shell element), chưa mở rộng trực tiếp cho các phần tử vỏ cong liên tục bậc cao (curved shell elements).
  • Đối với vật liệu FGM, luận án tập trung khảo sát dạng phân bố theo quy luật hàm lũy thừa (P-FGM), chưa đi sâu phân tích toàn diện cho các quy luật phân bố Sigmoid (S-FGM) hay hàm mũ siêu việt (E-FGM) trong toàn bộ các chương số.
  • Các khảo sát số tập trung vào ứng xử tĩnh tuyến tính/phi tuyến hình học, dao động tự do và ổn định tĩnh; chưa xét đến bài toán động lực học phi tuyến theo thời gian thực (transient dynamic) hay bài toán phá hủy vật liệu.

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

  1. Mở rộng các kỹ thuật phần tử SQ4H, SQ4T, SQ4C, SQ4P cho bài toán phân tích động lực học phi tuyến, phân tích va chạm và tải trọng nổ tác dụng lên kết cấu tấm/vỏ.
  2. Nghiên cứu tích hợp các kỹ thuật phần tử đề xuất vào bài toán phân tích đa trường vật lý tương tác (nhiệt - áp điện - từ giảo - cơ học).
  3. Phát triển các module tính toán dựa trên các phần tử này để tích hợp vào 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 kỹ thuật thực tế.

Giá trị tham khảo

  • Đối với nghiên cứu sinh, học viên cao học và nhà nghiên cứu ngành Cơ kỹ thuật, Cơ học tính toán: Cung cấp tài liệu tham khảo chi tiết về quy trình toán học thiết lập phần tử hữu hạn trơn (SFEM), kỹ thuật nội suy kép (TIS), kỹ thuật tổ hợp biến dạng xử lý khóa màng và kỹ thuật xấp xỉ Chebyshev.
  • Đối với giảng viên và kỹ sư kết cấu công trình, hàng không, tàu thủy: Cung cấp công thức phần tử vỏ phẳng có sườn gia cường và các thuật toán tính toán ứng xử kết cấu làm từ vật liệu composite, FGM và FGP-GPLs.
  • Phần đáng tham khảo nhất:
    • Chương 4: Quy trình thiết lập hàm dạng theo chiến lược nội suy kép TIS.
    • Chương 5: Kỹ thuật chọn điểm buộc và tổ hợp biến dạng để khử khóa màng trên lưới méo.
    • Chương 7: Hệ thống bảng số liệu kiểm chứng và so sánh sai số giữa các phần tử với nghiệm giải tích chuẩn.

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

1. Hiện tượng khóa cắt và khóa màng trong phân tích tấm/vỏ là gì và luận án giải quyết bằng cách nào?

Khóa cắt (shear locking) là hiện tượng năng lượng biến dạng cắt bị phóng đại quá mức khi chiều dày kết cấu rất mỏng trong lý thuyết Reissner-Mindlin, dẫn đến chuyển vị tính toán nhỏ hơn thực tế. Khóa màng (membrane locking) là hiện tượng độ cứng màng bị giả tạo khi phân tích vỏ cong bằng lưới chia thô hoặc méo dưới tải trọng uốn. Luận án giải quyết khóa cắt bằng cách sử dụng lý thuyết biến dạng cắt bậc cao $C^0$-HSDT (phần tử SQ4H), kỹ thuật điểm buộc ten-xơ biến dạng cắt rời rạc (phần tử SQ4C), và giải quyết khóa màng bằng kỹ thuật chia miền con tam giác để tính biến dạng màng chuyển về 4 điểm buộc biên kết hợp làm trơn biến dạng (phần tử SQ4C).

2. Sự khác biệt căn bản giữa phần tử SQ4H và phần tử SQ4T là gì?

Phần tử SQ4H tập trung vào việc tích hợp lý thuyết biến dạng cắt bậc cao cải biên ($C^0$-HSDT) với kỹ thuật trơn biến dạng miền con ($n_c = 1, 2$) nhằm loại bỏ hệ số hiệu chỉnh cắt và ổn định hóa phân tích uốn phi tuyến hình học. Trong khi đó, phần tử SQ4T sử dụng chiến lược nội suy kép (TIS), xây dựng hàm dạng bậc cao thông qua việc xét thêm gradient trung bình và giá trị của nhóm $n_{sp}$ nút lân cận trong miền ảnh hưởng, nhằm khắc phục sự gián đoạn của trường ứng suất và biến dạng qua các biên phần tử.

3. Phần tử SQ4C xử lý bài toán kết cấu tấm/vỏ có sườn gia cường như thế nào?

Phần tử SQ4C tích hợp độ cứng của sườn gia cường trực tiếp vào phần tử vỏ phẳng thông qua việc ghép nối ma trận độ cứng của phần tử dầm/sườn với ma trận độ cứng màng và uốn của phần tử tấm/vỏ tại các bậc tự do tương ứng, cho phép phân tích chính xác ứng xử tĩnh, dao động và ổn định của tấm/vỏ có sườn ($n_{st}$ sườn) mà không đòi hỏi phải chia lưới quá mịn tại vị trí tiếp giáp giữa sườn và bản vỏ.

4. Kết cấu vật liệu FGP-GPLs trong luận án được khảo sát với các mô hình phân bố nào?

Vật liệu xốp phân lớp chức năng gia cường tiểu cầu graphene (FGP-GPLs) trong luận án được khảo sát với:

  • 3 dạng phân bố độ xốp theo chiều dày: Đối xứng (P-S), Bất đối xứng (P-A) và Đồng đều (P-U).
  • 3 dạng phân bố hàm lượng graphene: Đối xứng (GPL-S), Bất đối xứng (GPL-A) và Đồng đều (GPL-U). Các đặc tính cơ học hiệu dụng được xác định thông qua mô hình vi cơ học Halpin-Tsai kết hợp quy luật phân phối Gauss.

5. Vì sao luận án tập trung phát triển các phần tử tứ giác 4 nút thay vì các phần tử bậc cao 8 nút hay 9 nút?

Phần tử tứ giác 4 nút bậc thấp có ưu điểm lớn về sự đơn giản trong thiết lập công thức toán học, kích thước ma trận độ cứng nhỏ, thuận tiện trong việc chia lưới tự động và tiết kiệm chi phí tính toán (thời gian CPU). Tuy nhiên, phần tử 4 nút truyền thống dễ bị khóa cắt và khóa màng. Việc phát triển các kỹ thuật mới (trơn biến dạng, nội suy kép, tổ hợp biến dạng, đa thức Chebyshev) trên nền phần tử tứ giác 4 nút giúp đạt được độ chính xác tương đương các phần tử bậc cao nhưng vẫn duy trì tính hiệu quả và kinh tế về mặt tài nguyên tính toán.

Kết luận

Luận án đã nghiên cứu và phát triển thành công 4 kỹ thuật phần tử hữu hạn tứ giác 4 nút (SQ4H, SQ4T, SQ4C, SQ4P) ứng dụng cho phân tích ứng xử cơ học của kết cấu dạng tấm và vỏ. Các phần tử đề xuất đã giải quyết hiệu quả các hạn chế cố hữu của FEM truyền thống như hiện tượng khóa cắt, khóa màng, sự mất ổn định số trong phân tích phi tuyến và tính gián đoạn ứng suất liên phần tử. Hệ thống kết quả số đã kiểm chứng độ chính xác và độ tin cậy của các mô hình trên nhiều dạng kết cấu tấm phẳng, tấm gấp, vỏ cong phức tạp cấu tạo từ vật liệu đẳng hướng, composite nhiều lớp, FGM và FGP-GPLs. Đây là công trình nghiên cứu hoàn chỉnh, đóng góp các công cụ số hữu ích cho lĩnh vực cơ học tính toán và kỹ thuật công trình.