Tổng quan về luận án

Trong cơ học công trình và kỹ thuật kết cấu hiện đại, kết cấu dạng tấm và vỏ (plate/shell structures) giữ vai trò then chốt tại các hệ thống chịu lực trọng yếu như vỏ tàu ngầm, thân cánh máy bay, tấm chắn nhiệt hàng không vũ trụ, vỏ lò phản ứng, tháp giải nhiệt, si-lô và bồn chứa áp lực cao. Về mặt phân loại hình học, tấm và vỏ là các kết cấu 2 chiều (2D) có kích thước bề dày $h$ rất nhỏ so với các kích thước mặt phẳng $a, b$. Do sự phức tạp về mặt hình học cong, tính dị hướng của vật liệu composite đa lớp, tính biến thiên liên tục của vật liệu phân lớp chức năng (FGM) và vật liệu xốp nano (FGP-GPLs), các phương pháp giải tích thuần túy (analytical methods) như nghiệm chuỗi Fourier của Navier hay Levy chỉ giải quyết được các bài toán có dạng hình học và điều kiện biên lý tưởng. Phương pháp phần tử hữu hạn (Finite Element Method - FEM) là công cụ số vượt trội và phổ biến nhất để giải quyết các bài toán biên phức tạp này.

Tuy nhiên, FEM truyền thống dựa trên phép nội suy chuyển vị chuẩn $C^0$ liên tục đối mặt với các rào cản tính toán nan giải: hiện tượng khóa cắt (shear locking) khi độ mảnh $a/h \ge 100$ trong lý thuyết tấm Reissner-Mindlin; hiện tượng khóa màng (membrane locking) khi mô hình hóa vỏ mỏng chịu uốn chủ đạo bằng phần tử bậc thấp hoặc lưới biến dạng (distorted mesh); sự gián đoạn ứng suất và biến dạng liên phần tử qua các biên chia; và chi phí tính toán tăng vọt khi áp dụng lý thuyết biến dạng cắt bậc cao (HSDT) đòi hỏi tính liên tục $C^1$. Luận án tiến sĩ ngành Cơ kỹ thuật (mã số 9520101) của nghiên cứu sinh Tôn Thất Hoàng Lân, 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 tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh (2022), mang tên: "Phát triển các kỹ thuật phần tử hữu hạn cho phân tích kết cấu dạng tấm và vỏ", đã trực tiếp giải quyết các khoảng trống học thuật này.

Nghiên cứu tập trung giải quyết 4 câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Làm thế nào để xây dựng một phần tử tứ giác 4 nút bậc thấp khử hoàn toàn khóa cắt và thỏa mãn lý thuyết cắt bậc cao HSDT mà chỉ yêu cầu tính liên tục $C^0$ phục vụ phân tích phi tuyến hình học?
  2. RQ2: Kỹ thuật nội suy kép (Twice Interpolation Strategy - TIS) có thể khắc phục triệt để 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ử tấm/vỏ mà không làm tăng bậc tự do toàn cục hay không?
  3. RQ3: Làm thế nào để loại bỏ đồng thời khóa cắt và khóa màng cho kết cấu tấm/vỏ có sườn gia cường phức tạp trên lưới méo thông qua kỹ thuật tổ hợp biến dạng và làm trơn miền con?
  4. RQ4: Đa thức trực giao Chebyshev có thể được tích hợp vào công thức phần tử tấm/vỏ như thế nào để đạt khả năng tinh chỉnh bậc $p$ linh hoạt và phân tích hiệu quả vật liệu tiên tiến FGP gia cường tiểu cầu graphene (GPLs)?

Phạm vi nghiên cứu bao quát các bài toán phân tích uốn tĩnh tuyến tính, dao động tự do, phân tích mất ổn định (buckling), và uốn phi tuyến hình học (theo giả thuyết biến dạng lớn von Kármán) cho tấm phẳng, tấm gấp, vỏ trụ, vỏ cầu, vỏ yên ngựa và tấm/vỏ có sườn gia cường. Đối tượng vật liệu khảo sát trải rộng từ vật liệu đẳng hướng, composite nhiều lớp đối xứng/lệch trục, vật liệu chức năng gradient P-FGM, đến vật liệu xốp FGP với 3 dạng phân bố lỗ rỗng (P-S, P-A, P-U) kết hợp 3 mô hình phân bố tiểu cầu graphene (GPL-S, GPL-A, GPL-U).

Literature Review và Positioning

Khung lý thuyết cho phân tích tấm/vỏ đã trải qua nhiều giai đoạn phát triển với các cuộc tranh luận học thuật sâu sắc:

                            [Trường phái Lý thuyết Tấm & Vỏ]
                                          │
            ┌─────────────────────────────┴─────────────────────────────┐
            ▼                                                           ▼
  Lý thuyết Cổ điển & Cắt bậc 1                               Lý thuyết Biến dạng Cắt bậc cao
  • Kirchhoff-Love (CPT): C¹                                  • Reddy (1984): HSDT (C¹)
  • Reissner (1945), Mindlin (1951): FSDT (C⁰)                • C⁰-HSDT Formulation
    └─ Khóa cắt (Shear Locking) khi a/h lớn                     └─ Khử hệ số hiệu chỉnh cắt k_s
            │                                                           │
            └─────────────────────────────┬─────────────────────────────┘
                                          ▼
                         [Chiến lược Phần tử Hữu hạn Hiện đại]
                                          │
      ┌───────────────────┬───────────────┴───────────────┬───────────────────┐
      ▼                   ▼                               ▼                   ▼
MITC & DSG (Bathe)   SFEM (Liu, Nguyen-Xuan)     TIS (Bui, Wu)         IGA (Hughes)
MITC4/MITC8          CS-FEM / ES-FEM             Hàm dạng bậc cao      NURBS cơ sở
Trói buộc ten-xơ     Làm trơn biến dạng          Liên tục ứng suất     Biên chính xác

Trường phái thứ nhất bắt nguồn từ lý thuyết tấm cổ điển Kirchhoff-Love (Classical Plate Theory - CPT) bỏ qua biến dạng cắt ngang, đòi hỏi hàm dạng liên tục đạo hàm bậc nhất $C^1$, gây khó khăn khi xây dựng phần tử tứ giác 4 nút. Để khắc phục, lý thuyết biến dạng cắt bậc nhất FSDT của Reissner (1945) và Mindlin (1951) chỉ yêu cầu tính liên tục $C^0$, chuyển vị thẳng và góc xoay độc lập. Tuy nhiên, FSDT lại dẫn đến hiện tượng khóa cắt trầm trọng khi bề dày $h \to 0$ và phụ thuộc vào hệ số hiệu chỉnh cắt $k_s = 5/6$. Trường phái thứ hai phát triển lý thuyết cắt bậc cao HSDT (Reddy, 1984, 2004), mô tả phân bố ứng suất cắt dạng parabol qua chiều dày, triệt tiêu ứng suất cắt tại hai bề mặt đáy/đỉnh mà không cần hệ số $k_s$. Dẫu vậy, dạng thiết lập nguyên bản của HSDT đòi hỏi các đạo hàm bậc cao, khiến việc xây dựng phần tử hữu hạn 4 nút $C^0$ truyền thống trở nên vô cùng phức tạp.

Về mặt công nghệ phần tử số, các nghiên cứu quốc tế đã đề xuất nhiều giải pháp:

  1. Phương pháp nội suy ten-xơ hỗn hợp (Mixed Interpolation Tensorial Components - MITC): Bathe và Dvorkin (1985, 1986) giới thiệu phần tử MITC4, sau đó mở rộng thành MITC8, MITC9, MITC16. Gần đây, Lee và Bathe (2004, 2014) phát triển MITC4+ để tăng cường độ hội tụ và giảm khóa màng trên lưới méo.
  2. Phương pháp phần tử hữu hạn trơn (Smoothed Finite Element Methods - SFEM): Khởi xướng bởi Liu và cộng sự (2007), phát triển mạnh mẽ qua các công trình của Nguyen-Xuan và cộng sự (2008, 2010), Nguyen-Thoi và cộng sự (2012, 2013) với các biến thể 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/CSMIN3, CS-DSG3, MISQ20, MISQ24). SFEM làm mềm độ cứng quá mức (overly stiff) của FEM chuẩn, cải thiện tốc độ hội tụ và độ chính xác của trường ứng suất.
  3. Phương pháp phân tích đẳng hình học (IsoGeometric Analysis - IGA): Hughes và cộng sự (2005), Bazilevs và cộng sự (2010), Trần Văn và cộng sự ứng dụng hàm cơ sở NURBS để đồng nhất hóa mô hình hình học CAD và mô hình tính toán. Nhược điểm của IGA nằm ở tính toàn cục của hàm cơ sở, khiến việc tích phân Gauss và xử lý điều kiện biên cục bộ tiêu tốn nhiều tài nguyên tính toán.
  4. Phương pháp phần tử phổ và đa thức trực giao: Zrahia và cộng sự (2008), Sprague và Geers (2003) nghiên cứu phần tử phổ PSE dựa trên điểm Gauss-Legendre-Lobatto; Liu, He, Đặng Trung và cộng sự (2014-2018) khai thác tính trực giao của chuỗi Chebyshev để giải phương trình vi phân kết cấu.
  5. Kỹ thuật nội suy kép (Twice Interpolation Strategy - TIS): Bùi và cộng sự (2011), Wu và cộng sự (2014), Zheng và cộng sự (2016) áp dụng TIS cho bài toán đàn hồi 2D phẳng để đạt trường ứng suất liên tục qua biên phần tử.

So với các nghiên cứu quốc tế kể trên, luận án đã định vị một hướng tiếp cận độc đáo: tích hợp các ưu điểm của kỹ thuật làm trơn miền con CS-FEM, kỹ thuật nội suy kép TIS, chiến lược tổ hợp biến dạng (combined strain), và chuỗi đa thức Chebyshev vào họ 4 phần tử tứ giác 4 nút bậc thấp (SQ4H, SQ4T, SQ4C, SQ4P). Hệ phần tử này giữ nguyên tính đơn giản của cấu trúc 4 nút trực quan, tương thích tuyệt đối với các phần mềm FEM thương mại hiện hành, đồng thời giải quyết triệt để các dạng khóa số và hiện tượng gián đoạn ứng suất.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và tái cấu trúc các lý thuyết cơ học kết cấu kinh điển thông qua 4 đóng góp đột phá:

Thứ nhất, mở rộng lý thuyết biến dạng cắt bậc cao HSDT của Reddy thành dạng $C^0$-HSDT phi tuyến: Bằng cách biến đổi trường chuyển vị để cô lập các thành phần đạo hàm bậc cao vào các bậc tự do xoay ảo, luận án đã thiết lập thành công mô hình phần tử tứ giác 4 nút $C^0$ thỏa mãn điều kiện triệt tiêu ứng suất tiếp tuyến tại mặt phẳng trên và dưới ($z = \pm h/2$) mà không cần hàm dạng liên tục $C^1$. Khi tích hợp phương trình phi tuyến hình học von Kármán: $$\varepsilon_{xx} = \varepsilon_{xx}^0 + z \kappa_{xx}^0 + z^3 \kappa_{xx}^* + \frac{1}{2}\left(\frac{\partial w}{\partial x}\right)^2$$ mô hình đã phản ánh chính xác hiệu ứng tăng cứng màng (membrane stiffening) ở miền chuyển vị lớn của tấm phẳng và tấm gấp.

Thứ hai, phát triển lý thuyết trường liên tục đa phần tử thông qua Twice Interpolation Strategy (TIS): Luận án chứng minh rằng việc xây dựng hàm nội suy bậc hai dựa trên giá trị nút kết hợp gradient trung bình của các nút lân cận trong miền ảnh hưởng ($nsp$) tạo ra trường chuyển vị có đạo hàm liên tục bậc cao. Điều này phá vỡ hạn chế cố hữu của phần tử $C^0$ tiêu chuẩn—vốn có đạo hàm bước nhảy tại biên phần tử gây gián đoạn trường ứng suất $\sigma_{xx}, \sigma_{yy}, \tau_{xy}$.

Thứ ba, thiết lập mô hình cơ học vi mô cho vật liệu xốp FGP gia cường tiểu cầu graphene (GPLs): Dựa trên mô hình Halpin-Tsai kết hợp định luật hỗn hợp mở rộng (extended rule of mixture), luận án đã thiết lập hệ phương trình liên hệ ứng suất - biến dạng biểu diễn chính xác ảnh hưởng đồng thời của hệ số xốp $e_0$ và tỷ phần khối lượng hạt nano $WF_{GPL}$: $$E(z) = E_s [1 - e_0 \lambda(z)], \quad \rho(z) = \rho_s [1 - e_m \lambda(z)]$$ trong đó $\lambda(z)$ tương ứng với các phân bố đối xứng (P-S), bất đối xứng (P-A), và đều (P-U).

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đa lý thuyết và phương pháp tính toán:

                  ┌─────────────────────────────────────────────────────────┐
                  │                   THIẾT KẾ PHẦN TỬ                      │
                  └────────────────────────────┬────────────────────────────┘
                                               │
         ┌───────────────────┬─────────────────┴─────────────────┬───────────────────┐
         ▼                   ▼                                   ▼                   ▼
     [ SQ4H ]            [ SQ4T ]                            [ SQ4C ]            [ SQ4P ]
  ┌──────────────┐    ┌──────────────┐                    ┌──────────────┐    ┌──────────────┐
  │ • C⁰-HSDT    │    │ • TIS Scheme │                    │ • Combined   │    │ • Chebyshev  │
  │ • CS-FEM     │    │ • Nodal Grad │                    │   Strain     │    │   Polynomial │
  │   nc = 1, 2  │    │ • Local Sub- │                    │ • CS-FEM     │    │ • p-Order    │
  │ • Geometric  │    │   domain     │                    │ • Tying Pts  │    │   Refinement │
  │   Nonlinear  │    │ • Continuous │                    │ • Stiffened  │    │ • FGP-GPLs   │
  │   Analysis   │    │   Stress     │                    │   Shells     │    │   Modeling   │
  └──────────────┘    └──────────────┘                    └──────────────┘    └──────────────┘
         │                   │                                   │                   │
         └───────────────────┼───────────────────────────────────┴───────────────────┘
                             ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │                   PHÂN TÍCH KẾT CẤU                     │
                  │  • Uốn tĩnh (Static Bending)  • Dao động tự do (Free    │
                  │    Linear/Nonlinear             Vibration Analysis)     │
                  │  • Mất ổn định (Buckling)     • Vỏ mỏng, Tấm gấp,       │
                  │    Uniaxial/Biaxial             Kết cấu có sườn         │
                  └─────────────────────────────────────────────────────────┘

  1. Phần tử SQ4H (Smoothed Q4 based on $C^0$-HSDT): Kết hợp kỹ thuật trơn biến dạng trên miền con tứ giác ($nc = 1, 2$) với phương trình động học $C^0$-HSDT. Ma trận độ cứng tiếp tuyến $K_T$ được làm mềm thông qua toán tử trơn biến dạng biến phân: $$\tilde{\boldsymbol{\varepsilon}} = \int_{\Omega_k} \boldsymbol{\varepsilon}(\mathbf{x}) \Phi(\mathbf{x}) d\Omega$$ giúp khử hoàn toàn khóa cắt mà không cần bất kỳ hệ số hiệu chỉnh nhân tạo nào.

  2. Phần tử SQ4T (Twice Interpolation Strategy Q4): Sử dụng chiến lược nội suy 2 bước. Bước 1: Tính toán gradient nút từ các phần tử lân cận. Bước 2: Bổ sung các số hạng gradient vào hàm dạng cơ sở để nâng bậc xấp xỉ trường chuyển vị trong lòng phần tử tứ giác 4 nút, tạo nên phân bố ứng suất trơn tru trên toàn bộ kết cấu màng Cook và tấm chịu uốn.

  3. Phần tử SQ4C (Combined Strain Q4): Sử dụng chiến lược tách biệt và tái tổ hợp 3 trạng thái biến dạng: màng ($\boldsymbol{\varepsilon}_m$), uốn ($\boldsymbol{\varepsilon}_b$) và cắt ngang ($\boldsymbol{\gamma}_s$). Kỹ thuật 4 điểm buộc (tying points) được tích hợp trên các biên phần tử để triệt tiêu các mode năng lượng giả (spurious zero-energy modes), khử đồng thời hiện tượng khóa màng khi chia lưới thô/méo trên vỏ cong và khóa cắt khi tấm chuyển sang giới hạn siêu mỏng ($a/h = 10^4$).

  4. Phần tử SQ4P (Chebyshev Polynomial Q4): Xây dựng hàm dạng 2D từ tích tensor của các đa thức trực giao Chebyshev loại một $T_p(\xi)$: $$N_i(\xi, \eta) = \phi_m(\xi) \psi_n(\eta)$$ Cho phép nâng bậc đa thức $p_1, p_2$ từ bậc 3 lên bậc 6 mà không cần thay đổi cấu trúc lưới hình học, đạt độ hội tụ dạng phổ (spectral convergence) cho các bài toán phân tích tần số dao động riêng bậc cao và tải trọng tới hạn đa trục.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý Thực chứng (Positivism) kết hợp Chủ nghĩa hiện thực có tính phản biện (Critical Realism) trong cơ học tính toán. Mọi mô hình toán học và phần tử số đều được kiểm chứng chéo thông qua hệ thống phương trình biến phân năng lượng toàn phần $\Pi$: $$\delta \Pi = \delta U - \delta W = 0$$ trong đó $U$ là nội năng biến dạng đàn hồi và $W$ là công ngoại lực.

Thiết kế nghiên cứu thực hiện theo quy trình 4 giai đoạn khép kín:

  • Giai đoạn 1: Thiết lập phương trình vi phân chủ đạo và chuyển đổi trường chuyển vị sang dạng biến phân rời rạc hóa.
  • Giai đoạn 2: Lập trình thuật toán phần tử hữu hạn bằng ngôn ngữ MATLAB/Fortran, xây dựng các hàm tính ma trận độ cứng $\mathbf{K}$, ma trận khối lượng nhất quán $\mathbf{M}$, và ma trận độ cứng hình học $\mathbf{K}_g$.
  • Giai đoạn 3: Mô phỏng số các bài toán chuẩn mực (Standard Benchmarks) quốc tế: bài toán màng Cook (Cook's membrane), vỏ trụ Scordelis-Lo (Scordelis-Lo roof), vỏ bán cầu có lỗ khoét $18^\circ$ (Hemispherical shell), vỏ trụ màng cứng (Pinched cylinder).
  • Giai đoạn 4: Khảo sát tham số chuyên sâu trên các vật liệu tiên tiến FGM và FGP-GPLs, phân tích sai số tương đối và độ nhạy của lưới chia.

Quy trình nghiên cứu rigorous

Quy trình giải quyết các bài toán cơ học phi tuyến và đặc trưng động học được thực hiện nghiêm ngặt qua các thuật toán giải số chuẩn:

  1. Bài toán uốn phi tuyến hình học: Áp dụng thuật toán lặp Newton-Raphson hiệu chỉnh kết hợp kỹ thuật kiểm soát độ dài cung (arc-length method) để vượt qua các điểm giới hạn (limit points) và điểm rẽ nhánh (bifurcation points): $$\mathbf{K}T(\mathbf{q}i) \Delta \mathbf{q}{i+1} = \mathbf{F}{ext} - \mathbf{F}_{int}(\mathbf{q}_i)$$ Tiêu chuẩn hội tụ lực và chuyển vị được kiểm soát chặt chẽ với sai số dư $\epsilon \le 10^{-4}$.
  2. Bài toán dao động tự do: Giải bài toán trị riêng tổng quát: $$(\mathbf{K} - \omega^2 \mathbf{M})\mathbf{q} = \mathbf{0}$$ sử dụng thuật toán phân rã subspace hoặc Lanczos để xác định chính xác 6 tần số tự nhiên đầu tiên $\omega_i$ và các dạng dao động tương ứng.
  3. Bài toán ổn định (Buckling): Xác định tải trọng tới hạn $P_{cr}$ và hệ số tải trọng $\lambda_{cr}$ qua phương trình đặc trưng: $$(\mathbf{K} - \lambda_{cr} \mathbf{K}_g)\mathbf{q} = \mathbf{0}$$

Data và phân tích

Hệ thống kết quả số trong luận án được kiểm chứng dựa trên hàng loạt bộ dữ liệu đối chuẩn:

  • Phân tích tấm vuông composite nhiều lớp: Khảo sát các cấu hình xếp lớp đối xứng và phản đối xứng $[0^\circ/90^\circ/90^\circ/0^\circ]$, $[\theta^\circ/-\theta^\circ]$, $[0^\circ/90^\circ/0^\circ/90^\circ/0^\circ]$ với tỷ số mô-đun đàn hồi sợi $E_1/E_2 = 10, 20, 30, 40$; tỷ số kích thước/chiều dày biến thiên từ tấm rất dày ($a/h = 5, 10$) đến siêu mỏng ($a/h = 100, 1000, 10000$).
  • Phân tích kết cấu tấm/vỏ có sườn gia cường: Mô hình hóa sườn gia cường dọc và ngang bằng kỹ thuật phần tử dầm lệch tâm tương thích, đánh giá ảnh hưởng của số lượng sườn $n_{st}$, chiều cao sườn $t_{st}$, và độ cứng uốn tương đối của sườn lên tải trọng mất ổn định đơn trục và hai trục.
  • Phân tích vật liệu FGP-GPLs: Đánh giá sự phụ thuộc của độ võng không thứ nguyên $w^$, tần số dao động chuẩn hóa $\bar{\omega} = \omega (a^2/h)\sqrt{\rho_m/E_m}$, và lực tới hạn chuẩn hóa $P^ = P_{cr} a^2 / (E_m h^3)$ theo hệ số xốp $e_0 \in [0.1, 0.7]$ và hàm lượng khối lượng hạt nano graphene $WF_{GPL} \in [0.0%, 1.0%]$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Khử triệt để hiện tượng khóa cắt và khóa màng trên mọi dải chiều dày: Phần tử SQ4C và SQ4H duy trì độ chính xác tuyệt đối khi tỷ số độ mảnh $a/h$ tăng từ $10$ lên $10000$. Tại bài toán tấm vuông chịu tải phân bố đều, sai số độ võng của SQ4C so với nghiệm giải tích chính xác của Pagano và Reddy duy trì dưới $0.15%$ mà không xuất hiện bất kỳ dấu hiệu cứng hóa số học (numerical locking) hay dao động áp lực giả.

  2. Nâng cao vượt bậc tính liên tục của trường ứng suất nhờ SQ4T: Trong bài toán kinh điển màng Cook (Cook's membrane problem)—vốn chịu ứng xử uốn kết hợp cắt cao làm biến dạng lưới tứ giác nghiêm trọng—phần tử SQ4T cho thấy trường ứng suất tương đương von Mises hoàn toàn liên tục và trơn tru qua các đường ranh giới phần tử ngay từ lưới thô ($4 \times 4$), trong khi phần tử tiêu chuẩn MITC4 thể hiện sự gián đoạn rõ rệt giữa các biên phần tử lân cận.

  3. Khả năng hội tụ phổ linh hoạt của SQ4P mà không cần chia lại lưới: Bằng cách tăng bậc đa thức Chebyshev từ $p=3$ lên $p=5$, độ võng chuẩn hóa và 6 mode dao động đầu tiên của vỏ bán cầu có lỗ khoét $18^\circ$ hội tụ tiệm cận về nghiệm tham chiếu chính xác cao với số lượng phần tử giảm hơn $60%$ so với việc tinh chỉnh lưới theo phương pháp $h$-refinement của FEM truyền thống.

  4. Phát hiện cơ chế tối ưu hóa độ cứng của vật liệu xốp FGP-GPLs: Kết quả số chỉ ra rằng dạng phân bố lỗ xốp đối xứng (P-S) kết hợp với cấu trúc phân bố hạt nano graphene tập trung ở hai bề mặt ngoài (GPL-S) đem lại hiệu quả tăng cứng uốn và kháng mất ổn định cao nhất. Cụ thể, khi bổ sung chỉ $1.0\ \text{wt.%}$ GPLs vào nền kim loại xốp với hệ số rỗng $e_0 = 0.5$, tải trọng tới hạn $P^*$ của tấm tăng lên hơn $185%$ so với trường hợp không có hạt gia cường.

  5. Hiệu ứng bất đối xứng hình học và phân bố phi đối xứng: Cấu hình phân bố lỗ xốp P-A và GPL-A tạo ra sự lệch trục trung hòa và gây xuất hiện hiệu ứng liên kết uốn-màng (bending-extension coupling), làm giảm tần số dao động cơ bản của kết cấu từ $15%$ đến $28%$ so với cấu hình phân bố đối xứng P-S.

+---------------------------------------------------------------------------------------------------+
|               BẢNG SO SÁNH TỔNG HỢP NĂNG LỰC HỌC THUẬT CỦA 4 PHẦN TỬ ĐỀ XUẤT                      |
+==========+======================+=========================+====================+==================+
| Phần tử  | Nền tảng Lý thuyết   | Khắc phục Lỗi Số Học    | Loại Phân tích     | Ứng dụng Vật liệu|
+----------+----------------------+-------------------------+--------------------+------------------+
| SQ4H     | C⁰-HSDT + CS-FEM     | Shear locking,          | Phi tuyến hình học | Composite,       |
|          | (Làm trơn nc=1, 2)   | Triệt tiêu bậc C¹       | Tấm phẳng, Tấm gấp | Tấm gấp FGM      |
+----------+----------------------+-------------------------+--------------------+------------------+
| SQ4T     | Twice Interpolation  | Inter-element stress    | Tuyến tính &       | Đẳng hướng, FGM, |
|          | Strategy (TIS)       | discontinuity           | Phi tuyến uốn/vỏ   | Composite        |
+----------+----------------------+-------------------------+--------------------+------------------+
| SQ4C     | Combined Strain      | Membrane locking &      | Uốn tĩnh, Dao động,| Tấm/Vỏ có sườn,  |
|          | Strategy + CS-FEM    | Shear locking           | Buckling           | Vỏ cong phức tạp |
+----------+----------------------+-------------------------+--------------------+------------------+
| SQ4P     | Chebyshev Orthogonal | Locking-free via        | Dao động bậc cao,  | FGP-GPLs xốp,    |
|          | Polynomial Series    | p-refinement            | Buckling đa trục   | Composite nano   |
+----------+----------------------+-------------------------+--------------------+------------------+

Implications đa chiều

  • Về mặt học thuật: Luận án khẳng định tính khả thi vượt trội của việc chuyển đổi lý thuyết biến dạng cắt bậc cao HSDT sang dạng tương thích $C^0$ thông qua kỹ thuật làm trơn biến dạng trên miền con (CS-FEM). Mô hình này mở ra cánh cửa mới để giải các bài toán phi tuyến phức tạp mà không phải chịu gánh nặng tính toán của phần tử $C^1$ hay sự phức tạp lưới của IGA.
  • Về mặt phương pháp luận: Sự kết hợp giữa phép nội suy kép (TIS) và đa thức phổ Chebyshev cung cấp một công cụ mạnh mẽ để xử lý các bài toán cơ học có gradient ứng suất cục bộ cao (như xung quanh lỗ khoét, góc nhọn, hoặc vùng tập trung ứng suất tiếp xúc).
  • Về mặt thực tiễn kỹ thuật: Các kỹ sư kết cấu có thể ứng dụng trực tiếp các ma trận phần tử SQ4H, SQ4T, SQ4C vào các phần mềm tính toán kết cấu chuyên dụng (như SAP2000, ANSYS, ABAQUS User-defined Element - UEL) để thiết kế các kết cấu vỏ mỏng hàng không vũ trụ và bồn bể công nghiệp với độ tin cậy cao và chi phí phần cứng tối thiểu.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Mô hình vật lý: Nghiên cứu chỉ tập trung vào phân tích phi tuyến hình học (biến dạng lớn von Kármán) trong miền đàn hồi tuyến tính của vật liệu, chưa xét đến tính phi tuyến phi đàn hồi (elastoplasticity), hiện tượng phá hủy (damage mechanics), hay hiện tượng tách lớp (delamination) trong vật liệu composite nhiều lớp.
  2. Điều kiện môi trường: Ảnh hưởng tương tác của môi trường nhiệt - ẩm phức tạp (hygrothermal effects) và hiện tượng rão (creep) theo thời gian chưa được tích hợp đồng thời vào mô hình phần tử.
  3. Phân tích tương tác động học: Các bài toán tương tác phức tạp như tương tác chất lỏng - kết cấu (Fluid-Structure Interaction - FSI) trong bể chứa dao động tự do có chất lỏng sóng sánh hay bài toán va đập tốc độ cao chưa nằm trong phạm vi khảo sát.

Chương trình nghiên cứu tương lai (Future Research Agenda) bao gồm:

  • Hướng 1: Mở rộng hệ phần tử SQ4H và SQ4C sang phân tích phi tuyến vật liệu dẻo và mô hình phá hủy liên tục (Continuum Damage Mechanics) cho composite nhiều lớp.
  • Hướng 2: Tích hợp thuật toán tối ưu hóa hình học (Topology Optimization) dựa trên phần tử SQ4P để thiết kế phân bố lỗ xốp và hàm lượng hạt nano graphene tối ưu cho các cấu kiện hàng không vũ trụ.
  • Hướng 3: Phát triển các phần tử vỏ 3D liên tục (solid-shell elements) dựa trên kỹ thuật nội suy kép TIS để mô phỏng chính xác ứng suất pháp vuông góc với mặt phẳng tấm ($\sigma_{zz}$) trong các kết cấu rất dày.
  • Hướng 4: Mở rộng thuật toán sang bài toán tương tác đa trường vật lý (Nhiệt - Điện - Từ - Cơ) cho vật liệu áp điện thông minh (Smart Piezoelectric Materials).

Tác động và ảnh hưởng

Luận án của tác giả Tôn Thất Hoàng Lân đã tạo ra tác động học thuật và ứng dụng sâu rộng:

  • Chất lượng công bố quốc tế: Công trình đã công bố 5 bài báo khoa học trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (Q1/Q2), bao gồm:
    • Journal of Sandwich Structures & Materials (ISI)
    • International Journal of Computational Methods (ISI)
    • Archive of Applied Mechanics (ISI)
    • Periodica Polytechnica Civil Engineering (ISI)
    • Journal of Applied and Computational Mechanics (Scopus) Cùng nhiều bài báo trên kỷ yếu hội nghị cơ học toàn quốc và Springer quốc tế.
  • Tác động ngành và công nghệ: Cung cấp giải pháp phần mềm lõi giúp các viện nghiên cứu thiết kế tàu thủy, hàng không vũ trụ và các tập đoàn xây dựng tối ưu hóa quá trình tính toán kết cấu tấm/vỏ phức tạp, giảm thiểu sai số tính toán từ $10-15%$ của các phần tử thương mại thông thường xuống dưới $1%$.
  • Đóng góp xã hội: Giúp giảm khối lượng kết cấu công trình thông qua việc sử dụng hiệu quả vật liệu tiên tiến FGM/FGP-GPLs, thúc đẩy phát triển bền vững và tiết kiệm tài nguyên vật liệu trong công nghiệp chế tạo.

Đối tượng hưởng lợi

                                [HỆ SINH THÁI HƯỞNG LỢI TỪ LUẬN ÁN]
                                                 │
          ┌──────────────────────────┬───────────┴───────────┬──────────────────────────┐
          ▼                          ▼                       ▼                          ▼
  [Nghiên cứu sinh &          [Chuyên gia &           [Bộ phận R&D             [Cơ quan Quản lý &
     Học viên Cao học]           Giảng viên Đại học]     Doanh nghiệp Kỹ thuật]   Hiệp hội Tiêu chuẩn]
  • Khung lý thuyết C⁰-HSDT   • Tài liệu giảng dạy    • Module UEL tích hợp    • Cơ sở dữ liệu khoa học
  • Mã nguồn giải thuật       • Nền tảng mở rộng        vào phần mềm CAE         cho tiêu chuẩn thiết kế
    phần tử hữu hạn nâng cao    đa vật lý               thương mại               vật liệu FGM mới
  • Nghiên cứu sinh và học viên cao học chuyên ngành Cơ kỹ thuật/Cơ học tính toán: Tiếp cận các thuật toán phần tử hữu hạn cải tiến, mã nguồn tính toán biến dạng trơn và phương pháp giải phi tuyến hình học mẫu mực.
  • Giảng viên và nhà khoa học: Sử dụng luận án như một tài liệu tham khảo chuyên sâu về kỹ thuật khử khóa cắt, khóa màng và mô hình hóa vật liệu cấu trúc nano tiên tiến.
  • Các kỹ sư R&D trong ngành Hàng không vũ trụ, Đóng tàu và Dầu khí: Khai thác công thức phần tử SQ4C và SQ4H để xây dựng các chương trình tính toán nội bộ (in-house codes) hoặc viết các module UEL cho ABAQUS/ANSYS phục vụ phân tích độ bền bồn áp lực và tấm gia cường.
  • Các nhà hoạch định chính sách và cơ quan tiêu chuẩn kỹ thuật: Có thêm cơ sở khoa học định lượng về ứng xử của vật liệu mới FGM và nanocomposite để xây dựng các quy chuẩn thiết kế kết cấu tiên tiến.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết nào của luận án mang tính độc đáo nhất và đã mở rộng lý thuyết kinh điển nào?

Đóng góp độc đáo nhất là việc thiết lập công thức 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 ($C^0$-HSDT) và kỹ thuật làm trơn biến dạng trên miền con (CS-FEM). Công trình đã mở rộng lý thuyết HSDT của J.N. Reddy (1984)—vốn đòi hỏi tính liên tục $C^1$ phức tạp—thành dạng biến phân $C^0$ khả thi cho phần tử 4 nút, loại bỏ hoàn toàn nhu cầu sử dụng hệ số hiệu chỉnh cắt $k_s$, đồng thời giải quyết bài toán phi tuyến hình học von Kármán với độ ổn định số học vượt bậc.

2. Sự cách tân về phương pháp luận của họ phần tử đề xuất so với các nghiên cứu quốc tế trước đây thể hiện ở điểm nào?

So với phần tử MITC4/MITC4+ của Bathe và cộng sự (chủ yếu dựa vào trói buộc trường biến dạng cắt tại các điểm nút) hay CSMIN3/MISQ24 của Liu và Nguyen-Van (chỉ áp dụng kỹ thuật làm trơn trên FSDT), phương pháp luận của luận án mang tính tổng hòa:

  • Phần tử SQ4C kết hợp tổ hợp cả 3 trường biến dạng (màng, uốn, cắt) với 4 điểm buộc và làm trơn miền con, khắc phục đồng thời cả hiện tượng khóa màng lẫn khóa cắt trên lưới méo.
  • Phần tử SQ4T đưa khái niệm nội suy kép TIS từ bài toán 2D phẳng sang cơ học kết cấu tấm/vỏ, đạt được tính liên tục của ứng suất qua biên phần tử mà các phần tử đẳng tham số thông thường không thể đạt được.

3. Phát hiện số học nào gây bất ngờ nhất trong quá trình phân tích kết cấu và được giải thích bằng cơ chế cơ học nào?

Phát hiện bất ngờ nhất xuất hiện ở phân tích tấm xốp FGP-GPLs phân bố bất đối xứng (P-A và GPL-A): Mặc dù hàm lượng vật liệu gia cường graphene được giữ nguyên so với dạng đối xứng (P-S), tần số dao động tự nhiên và tải trọng tới hạn lại sụt giảm tới $28%$. Nguyên nhân cơ học là do phân bố bất đối xứng tạo ra sự lệch tâm giữa mặt phẳng trung hòa hình học và mặt phẳng trung hòa cơ học, sinh ra hiệu ứng ghép nối màng - uốn (coupling effect) làm suy giảm độ cứng chống uốn tổng thể của tấm.

4. Quy trình lặp và tái lập kết quả (Replication Protocol) có được đảm bảo trong công trình không?

Luận án cung cấp đầy đủ và chi tiết toàn bộ các biểu thức giải tích của ma trận hàm dạng, ma trận liên hệ biến dạng - chuyển vị ($\mathbf{B}_m, \mathbf{B}_b, \mathbf{B}_s$), tọa độ các điểm tích phân Gauss, sơ đồ chia miền con ($nc=1, 2, 4$), và các tham số cơ lý của vật liệu. Bất kỳ nhà nghiên cứu nào cũng có thể lập trình tái lập chính xác $100%$ các kết quả uốn tĩnh, tần số dao động và tải trọng mất ổn định được công bố.

5. Lộ trình nghiên cứu 10 năm (10-Year Research Agenda) mở ra từ công trình này là gì?

Lộ trình 10 năm phát triển từ luận án tập trung vào: (1) Mở rộng hệ phần tử sang cơ học phá hủy và mô phỏng vết nứt lan truyền (XFEM-CSFEM); (2) Tích hợp học máy (Machine Learning / PINNs) để tối ưu hóa siêu tham số lưới và phân bố vật liệu nano; (3) Phát triển phần mềm tính toán kết cấu tấm/vỏ chuyên dụng phục vụ công nghiệp hàng không vũ trụ và công trình biển của Việt Nam.

Kết luận

Luận án tiến sĩ của NCS. Tôn Thất Hoàng Lân đã đóng góp toàn diện cho ngành Cơ kỹ thuật thông qua các kết quả cụ thể:

  1. Xây dựng thành công 4 phần tử tứ giác 4 nút mới (SQ4H, SQ4T, SQ4C, SQ4P) cho phân tích tuyến tính và phi tuyến kết cấu tấm/vỏ, đạt độ chính xác cao và triệt tiêu các hiện tượng khóa cắt, khóa màng.
  2. Hiện thực hóa lý thuyết biến dạng cắt bậc cao $C^0$-HSDT kết hợp CS-FEM trong phân tích phi tuyến hình học tấm phẳng và tấm gấp, nâng cao độ ổn định số học.
  3. Ứng dụng đột phá kỹ thuật nội suy kép (TIS) trong phân tích tấm/vỏ, giải quyết triệt để bài toán gián đoạn trường ứng suất và biến dạng liên phần tử.
  4. Thiết lập chiến lược tổ hợp biến dạng và chuỗi đa thức trực giao Chebyshev, mang lại giải pháp phân tích hiệu quả cho kết cấu tấm/vỏ có sườn gia cường và tấm xốp nano FGP-GPLs.
  5. Công bố 5 công trình ISI/Scopus chất lượng cao, khẳng định tầm vóc học thuật và khả năng ứng dụng thực tiễn của nghiên cứu trong kỷ nguyên số hóa cơ học kết cấu.