Tổng quan về luận án

Nghiên cứu của nghiên cứu sinh Tôn Thất Hoàng Lân (2022) 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 với đề tài "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ỏ" (Ngành Cơ kỹ thuật, Mã số: 9520101) đại diện cho một bước tiến tiên phong trong lĩnh vực cơ học tính toán (Computational Mechanics). Trong kỹ thuật hàng không vũ trụ, cơ khí chế tạo, công trình biển và xây dựng công trình ngầm, kết cấu tấm và vỏ (plate/shell structures) đóng vai trò xương sống chịu lực từ vỏ thân máy bay, cánh chịu tải, vỏ tàu ngầm, xilo, bồn chứa áp lực cao đến các kết cấu mái vỏ mỏng không gian. Mặc dù phương pháp phần tử hữu hạn (Finite Element Method - FEM) đã trở thành công cụ trụ cột để mô phỏng ứng xử kết cấu từ nhiều thập kỷ qua, các mô hình số hiện hành đối mặt với những thách thức nội tại nghiêm trọng khi phân tích các kết cấu thành mỏng và vật liệu tiên tiến.

Khoảng trống nghiên cứu (Research Gap) then chốt được luận án chỉ rõ là sự tồn tại dai dẳng của các hiện tượng khóa số học (numerical locking phenomena) trong các phần tử vỏ và tấm bậc thấp, đặc biệt là hiện tượng khóa cắt (shear locking) khi tỷ số độ dày trên nhịp $h/a \to 0$ và hiện tượng khóa màng (membrane locking) khi chia lưới thô hoặc lưới bị méo góc (distorted mesh) trong các kết cấu vỏ cong hai phương theo phân tích của Lee và Bathe (2004), Zienkiewicz và Taylor (2000). Hơn nữa, việc mở rộng từ lý thuyết biến dạng cắt bậc nhất (FSDT) của Mindlin-Reissner sang lý thuyết biến dạng cắt bậc cao (HSDT) của Reddy nhằm loại bỏ hệ số hiệu chỉnh cắt $k_s$ thường đòi hỏi hàm xấp xỉ liên tục bậc $C^1$, dẫn đến sự phức tạp bùng nổ về số bậc tự do (DOFs) và khó khăn khi thiết lập công thức phần tử. Mặt khác, các phần tử FEM tiêu chuẩn luôn tồn tại sự gián đoạn của trường ứng suất và biến dạng qua biên phần tử.

Luận án thiết lập hệ thống 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Làm thế nào để loại bỏ yêu cầu liên tục $C^1$ của lý thuyết HSDT mà vẫn đảm bảo mô phỏng chính xác ứng xử phi tuyến hình học của tấm và tấm gấp?
    H1: Kết hợp kỹ thuật trơn biến dạng trên miền con (Cell-based Strain Smoothing - CS) với việc tách biến dạng xoay độc lập sẽ tạo ra phần tử $C^0$-HSDT ổn định cao.
  • RQ2: Cơ chế nào giúp khắc phục sự gián đoạn của trường ứng suất và biến dạng liên phần tử mà không cần chuyển sang phương pháp không lưới (Meshfree) đắt đỏ?
    H2: Ứng dụng kỹ thuật nội suy kép (Twice Interpolation Strategy - TIS) tích hợp gradient trung bình nút lân cận sẽ tái lập tính liên tục của trường ứng suất qua biên.
  • RQ3: Làm sao để triệt tiêu hoàn toàn khóa màng và khóa cắt trong phân tích vỏ mỏng và vỏ có sườn gia cường trên lưới phần tử tứ giác 4 nút?
    H3: Tổ hợp biến dạng màng từ miền con tam giác với kỹ thuật nội suy biến dạng cắt ràng buộc (MITC) trên miền trơn tế bào sẽ kháng khóa màng triệt để.
  • RQ4: Liệu việc tích hợp chuỗi trực giao Chebyshev vào hàm dạng phần tử 4 nút có thể tối ưu hóa độ hội tụ cho 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)?
    H4: Hàm dạng dựa trên đa thức Chebyshev cho phép điều chỉnh đồng thời cả kích thước lưới ($h$-refinement) và bậc đa thức ($p$-refinement), đạt độ chính xác vượt trội.

Khung lý thuyết chủ đạo của luận án dựa trên lý thuyết đơn lớp tương đương (Equivalent Single Layer - ESL), kết hợp cơ học môi trường liên tục phi tuyến (Total Lagrangian Formulation), kỹ thuật phần tử hữu hạn trơn (Smoothed Finite Element Methods - S-FEM), lý thuyết nội suy kép (TIS) và lý thuyết xấp xỉ phổ Chebyshev.

Về mặt định lượng, luận án đã phát triển thành công 4 phần tử tứ giác 4 nút mới: SQ4H, SQ4T, SQ4C và SQ4P. Phạm vi nghiên cứu bao quát các dạng hình học từ tấm phẳng chữ nhật, tấm xiên (skew plates), tấm gấp ($60^\circ, 90^\circ$), vỏ trụ, vỏ cầu, vỏ hyperbol paraboloid (vỏ yên ngựa) chịu tải tĩnh, dao động tự do, mất ổn định tuyến tính và phi tuyến hình học lớn. Đối tượng vật liệu được khảo sát trải rộng từ vật liệu đẳng hướng, composite nhiều lớp đối xứng và lệch trục $[0^\circ/90^\circ/\dots]$, vật liệu phân lớp chức năng (P-FGM) đến vật liệu xốp phân cấp chức năng gia cường tấm nano graphene (FGP-GPLs) với các quy luật phân bố độ xốp P-S, P-A, P-U và phân bố graphene GPL-S, GPL-A, GPL-U.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp phần tử hữu hạn cho kết cấu tấm và vỏ chứng kiến sự tiến hóa liên tục của các dòng lý thuyết động học và kỹ thuật xấp xỉ số. Dòng nghiên cứu thứ nhất tập trung vào lý thuyết tấm: bắt đầu từ lý thuyết cổ điển Kirchhoff-Love ($C^1$), chuyển sang lý thuyết biến dạng cắt bậc nhất FSDT của Reissner (1945) và Mindlin (1951) ($C^0$), và đạt đến đỉnh cao với lý thuyết biến dạng cắt bậc cao HSDT của Reddy (1984, 2004). Mặc dù HSDT biểu diễn chính xác phân bố parabol của ứng suất cắt ngang và triệt tiêu ứng suất cắt tại hai bề mặt tự do mà không cần hệ số hiệu chỉnh cắt $k_s = 5/6$, việc triển khai trong FEM truyền thống bị cản trở do đòi hỏi đạo hàm bậc hai của chuyển vị ngang trong trường biến dạng.

Dòng nghiên cứu thứ hai tập trung vào việc giải quyết hiện tượng khóa cắt và khóa màng trong phần tử tứ giác 4 nút bậc thấp. Điển hình là kỹ thuật nội suy hỗn hợp các thành phần ten-xơ MITC4 của Bathe và Dvorkin (1985), MITC4+ của Lee và Bathe (2004, 2014), phần tử MIN4 của Tessler và Hughes (1983, 1985), hay phương pháp khoảng hở cắt rời rạc DSG của Bletzinger et al. (2000). Nhóm tác giả Nguyễn-Thanh, Nguyễn-Xuân, Rabczuk và Bordas (2008, 2011) đã tạo ra bước đột phá khi phát triển phương pháp phần tử hữu hạn trơn (S-FEM) bao gồm NS-FEM (theo nút), ES-FEM (theo cạnh) và CS-FEM (theo miền con tế bào), kết hợp với FSDT để cải thiện độ mềm của mô hình phần tử hữu hạn.

Trong bối cảnh học thuật quốc tế, tồn tại hai quan điểm và trường phái đối nghịch rõ nét:

  • Trường phái phương pháp đẳng hình học (IsoGeometric Analysis - IGA) do Hughes et al. (2005), Bazilevs et al. (2006, 2010), Nguyễn-Xuan et al. (2013, 2017) dẫn đầu: Ủng hộ việc sử dụng hàm cơ sở B-splines hoặc NURBS có độ trơn $C^p$ cao để mô hình hóa chính xác hình học CAD và loại bỏ hoàn toàn các dạng khóa số học. Tuy nhiên, nhược điểm chí mạng của IGA là hàm dạng phân bố toàn cục, việc tích phân đòi hỏi chuyển đổi không gian tham số phức tạp, khó áp đặt điều kiện biên thiết yếu và chi phí tính toán cực lớn đối với kết cấu phức tạp.
  • Trường phái phần tử hữu hạn lai/trộn bậc thấp (Low-Order Hybrid/Mixed Elements) do Bathe et al., Taylor et al., Liu et al. bảo vệ: Cho rằng các phần tử 4 nút với dạng đại số tường minh cục bộ mới là giải pháp tối ưu cho công nghiệp nhờ tính linh hoạt chia lưới và tương thích cao với cấu trúc dữ liệu thương mại, mặc dù phải đối mặt với nguy cơ xuất hiện các dạng dao động không năng lượng (spurious zero-energy modes hay hiệu ứng đồng hồ cát - hourglass phenomena) nếu áp dụng kỹ thuật tích phân giảm (reduced integration).

Vị thế nghiên cứu của luận án được xác lập tại điểm giao thoa chiến lược: Luận án không đi theo hướng mở rộng bậc phần tử lên 8 nút hay 16 nút (như HBQ8 của Darilmaz & Kumbasar, 2006 hay KUMBA của Kumbasar et al., 2000) vốn tốn kém tài nguyên tính toán; cũng không phụ thuộc vào không gian tham số toàn cục của IGA. Thay vào đó, tác giả định vị nghiên cứu vào việc tối ưu hóa cấu trúc phần tử tứ giác 4 nút phẳng thông qua việc tích hợp các công nghệ biến dạng tiên tiến nhất:

  1. So sánh với nghiên cứu của Taylor et al. (1987, 2000) về giải pháp giải tích và số cho tấm/vỏ, các phần tử của luận án đạt tốc độ hội tụ nhanh hơn trên các lưới thô.
  2. So sánh với các chuẩn kiểm chứng dao động tấm dày và mỏng của Abbassian et al. (1987) và Reddy (2004), phần tử SQ4H và SQ4P triệt tiêu hoàn toàn hiện tượng suy biến tần số dao động riêng và không làm phát sinh bất kỳ mode giả nào.

Đó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 sâu sắc các nền tảng lý thuyết cơ học kết cấu:

  • Cải biên lý thuyết biến dạng cắt bậc cao Reddy sang dạng $C^0$-HSDT: Tác giả đã loại bỏ thành công sự phụ thuộc vào đạo hàm bậc nhất của chuyển vị ngang $(\partial w_0/\partial x, \partial w_0/\partial y)$ trong trường chuyển vị HSDT bằng cách đưa vào các biến xoay độc lập $\theta_x, \theta_y$, cho phép xây dựng ma trận phần tử chỉ đòi hỏi hàm dạng $C^0$ song tuyến tính đơn giản. Trích xuất nguyên văn từ luận án: "Xây dựng phần tử tứ giác 4 nút SQ4H dựa vào kỹ thuật trơn biến dạng trên miền con kết hợp kỹ thuật cải biên dạng C0-HSDT để phân tích phi tuyến kết cấu tấm phẳng và tấm gấp. Phần tử này cải thiện độ chính xác của mô hình và giảm bớt sự bất ổn về số đối với phân tích hình học phi tuyến tính."
  • Lý thuyết nội suy kép (Twice Interpolation Strategy - TIS) cho cơ học tấm/vỏ: Mở rộng lý thuyết TIS từ bài toán biến dạng phẳng 2D (Bui et al., Wu et al., Zheng et al.) sang trường chuyển vị 5 bậc tự do của tấm/vỏ. Bằng cách thiết lập vùng hỗ trợ (support domain) $nsp$ xung quanh phần tử đang xét, hàm dạng bậc cao mới kết hợp giá trị nút và gradient trung bình nút được hình thành. Trích dẫn nguyên văn: "Với việc xây dựng hàm nội suy bậc cao dựa vào giá trị nút lẫn gradient trung bình nút trong phạm vi miền ảnh hưởng, phần tử này cải thiện được yếu tố bất liên tục của biến dạng và ứng suất qua biên của nó."
  • Thiết lập mô hình giải tích đồng nhất hóa vi mô cho cấu trúc FGP-GPLs: Ứng dụng mô hình Halpin-Tsai cải tiến kết hợp luật phân bố thể tích Voigt mở rộng để xác định chính xác mô đun đàn hồi hữu hiệu $E(z)$, mô đun cắt $G(z)$ và mật độ khối lượng $\rho(z)$ cho ba dạng phân bố xốp đối xứng (P-S), bất đối xứng (P-A) và đều (P-U), có kể đến sự thay đổi của hệ số độ xốp $e_0 = 1 - E_2'/E_1'$ và tỷ lệ khối lượng graphene $W_{\text{GPL}}$.

Hệ thống 4 mệnh đề lý thuyết cốt lõi (Propositions):

  • Mệnh đề 1: Kỹ thuật làm trơn biến dạng trên $nc$ miền con tế bào biến đổi tích phân diện tích của đạo hàm hàm dạng thành tích phân đường trên các cạnh biên miền con ($\int_{\Omega_C} \nabla \mathbf{N} d\Omega = \oint_{\Gamma_C} \mathbf{N} \mathbf{n} d\Gamma$), giúp loại bỏ hoàn toàn hiện tượng khóa cắt mà không cần áp dụng hệ số $k_s$.
  • Mệnh đề 2: Việc nội suy trường biến dạng thông qua gradient trung bình của các nút lân cận bảo toàn tính liên tục đạo hàm bậc nhất của trường chuyển vị, khắc phục bước nhảy ứng suất tại biên chung giữa các phần tử liền kề.
  • Mệnh đề 3: Phân tách độc lập ma trận độ cứng màng $\mathbf{K}_m$, uốn $\mathbf{K}_b$, cắt $\mathbf{K}_s$ và hình học phi tuyến $\mathbf{K}_g$ cho phép triệt tiêu khóa màng thông qua kỹ thuật gán điểm buộc (tied-points) trên các miền con tam giác.
  • Mệnh đề 4: Việc mở rộng trường chuyển vị theo chuỗi đa thức trực giao Chebyshev loại một $T_p(\xi)$ tạo ra khả năng tự điều chỉnh độ cứng phần tử mà không làm tăng bậc hình học không gian, cho phép kiểm soát sai số tiệm cận bậc cao.

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

Khung phân tích được xây dựng trên sự tích hợp chặt chẽ của 4 trục lý thuyết: (1) Động học biến dạng tấm/vỏ Reissner-Mindlin và Reddy, (2) Kỹ thuật làm trơn biến dạng tế bào (Cell-based Smoothing), (3) Chiến lược nội suy kép gradient nút, và (4) Đại số đa thức trực giao Chebyshev.

Điểm độc đáo mang tính cách mạng trong phương pháp giải tích là quá trình chuyển đổi toàn bộ tích phân thể tích phức tạp sang tích phân đường dọc theo biên của $nc$ miền con tứ giác ($nc = 1, 2, 4$), giải phóng thuật toán khỏi sự phụ thuộc vào ma trận chuyển đổi Jacobian nghịch đảo phức tạp khi phần tử bị méo mó.

Điều kiện biên và phạm vi áp dụng được xác định nghiêm ngặt: Kết cấu tấm và vỏ mỏng đến dày ($a/h$ biến thiên từ $5$ đến $10000$), góc xiên $\alpha \in [0^\circ, 60^\circ]$, góc gấp $\beta \in [60^\circ, 90^\circ]$, điều kiện biên tựa đơn 4 cạnh (SSSS), ngàm 4 cạnh (CCCC), ngàm-tự do (CFCF) và các dạng liên kết hỗn hợp.

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ủ thế giới quan thực chứng (Positivism) và lập trường nhận thức luận tất định (Deterministic Computational Mechanics). Toàn bộ quy trình nghiên cứu được xây dựng trên nền tảng phương pháp biến phân và nguyên lý thế năng toàn phần cực tiểu: $\delta \Pi = \delta (U - W) = 0$.

Thiết kế nghiên cứu được cấu trúc hóa theo phương pháp đa cấp độ (Multi-level Hierarchical Design):

  • Cấp độ vi mô/vật liệu: Mô hình hóa các đặc trưng cơ lý của vật liệu FGM thông qua hàm số mũ phân bố thể tích $V_c(z) = (2z+h/2h)^n$ và cấu trúc xốp FGP-GPLs bằng mô hình đàn hồi vi mô Halpin-Tsai.
  • Cấp độ phần tử (Local Formulation): Thiết lập ma trận chuyển vị nút $\mathbf{q}_e$, hàm dạng $N_i(\xi, \eta)$, trường biến dạng trơn $\tilde{\boldsymbol{\epsilon}}$, và xây dựng các ma trận độ cứng cục bộ $\mathbf{K}_m, \mathbf{K}_b, \mathbf{K}_s, \mathbf{K}_g, \mathbf{M}$.
  • Cấp độ hệ thống kết cấu (Global Assembly): Áp dụng ma trận chuyển đổi tọa độ $\mathbf{T}$ từ hệ trục cục bộ sang hệ trục tổng thể 3D, tích hợp bậc tự do xoay trong mặt phẳng (drilling degree of freedom $\theta_z$), khử các bậc tự do thừa và áp đặt điều kiện biên thiết yếu.

Quy trình nghiên cứu rigorous

Quy trình tính toán số được kiểm soát với độ tin cậy và tính lặp lại (reproducibility) tuyệt đối:

  • Kiểm soát tính hợp thức (Construct & Internal Validity): Mỗi phần tử đều trải qua các bài kiểm tra chuẩn quốc tế (Patch Tests) đối với biến dạng màng thuần túy, uốn thuần túy và cắt thuần túy nhằm đảm bảo phần tử thỏa mãn điều kiện hoàn chỉnh và tương thích biến dạng khi chia mịn lưới.
  • Chiến lược khử khóa số: Triển khai kỹ thuật điểm buộc (tied-points) tại 4 điểm biên kết hợp với tích phân làm trơn tế bào CS-FEM, đảm bảo không xảy ra hiện tượng khóa màng ngay cả khi tỷ lệ cạnh phần tử bị biến dạng méo mó nghiêm trọng.
  • Triệt tiêu mode dao động giả: Kiểm tra phổ trị riêng của ma trận độ cứng phần tử tự do không liên kết, chứng minh ma trận chỉ tồn tại đúng 6 bậc tự do chuyển động cứng (3 tịnh tiến, 3 quay) và không xuất hiện bất kỳ mode năng lượng zero giả nào ($\lambda_{\text{spurious}} = 0$).

Data và phân tích

Toàn bộ thuật toán phần tử hữu hạn của 4 phần tử SQ4H, SQ4T, SQ4C, SQ4P được lập trình độc lập và hoàn chỉnh trên môi trường MATLAB/Fortran, chạy thực nghiệm trên hệ thống máy tính phần cứng tiêu chuẩn: Intel(R) Core(TM) i7 CPU @ 2.80GHz, 16.00 GB RAM.

Các phương pháp phân tích nâng cao bao gồm:

  1. Phân tích uốn phi tuyến hình học: Ứng dụng thuật toán lặp Newton-Raphson hiệu chỉnh kết hợp công thức Total Lagrangian để giải hệ phương trình phi tuyến: $[\mathbf{K}L + \mathbf{K}{NL}(\mathbf{q})]\mathbf{q} = \mathbf{F}$.
  2. Phân tích dao động tự do: Giải bài toán trị riêng tổng quát: $(\mathbf{K} - \omega^2 \mathbf{M})\boldsymbol{\Phi} = \mathbf{0}$ sử dụng phương pháp lặp không gian con (Subspace Iteration Method) để trích xuất 6 tần số dao động tự nhiên đầu tiên $\omega_i$ và các dạng mode dao động tương ứng.
  3. Phân tích mất ổn định (Buckling): Xác định tải trọng tới hạn $P_{\text{cr}}$ thông qua bài toán trị riêng tuyến tính: $(\mathbf{K} - \lambda_{\text{cr}} \mathbf{K}_g)\boldsymbol{\Phi} = \mathbf{0}$.

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

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

  1. Hiệu năng giải quyết phi tuyến hình học vượt bậc của phần tử SQ4H:
    Đối với tấm vuông và tấm gấp chịu tải phân bố đều $q$, phần tử SQ4H thể hiện độ chính xác và tính ổn định số xuất sắc mà không cần hệ số hiệu chỉnh cắt. Trong bài toán kiểm chứng tấm vuông đẳng hướng tựa đơn (SSSS) với $a/h = 10$, độ võng chuẩn hóa chính giữa tấm $w^* = (100 E h^3 w_c) / [12 q a^4 (1 - \nu^2)]$ hội tụ tuyệt đối về giá trị $0.406237$ tại lưới $14 \times 14$, hoàn toàn trùng khớp với nghiệm giải tích chính xác của Taylor et al. [140].

  2. Khả năng làm mịn trường ứng suất của phần tử SQ4T và hạn chế hình học:
    Phần tử SQ4T giải quyết triệt để sự gián đoạn của ứng suất pháp $\sigma_x, \sigma_y$ và ứng suất tiếp $\tau_{xy}$ qua các cạnh phần tử trong bài toán màng Cook (Cook's Membrane Problem) và tấm composite nhiều lớp. Tuy nhiên, luận án đã chỉ ra một phát hiện phản trực giác quan trọng: Trong phân tích vỏ cầu cong hai phương, thời gian tính toán của SQ4T tăng vọt lên $19.6400\text{ s}$ (cho cấu hình $\approx 1600$ DOFs, so với các phần tử khác chỉ mất từ $0.85$ đến $1.5\text{ s}$) do chi phí tìm kiếm các nút trong miền ảnh hưởng $nsp$ tăng theo hàm mũ, đồng thời kết quả có xu hướng cứng hóa cục bộ khi lưới bị méo cong.

  3. Phần tử SQ4C kháng khóa màng hoàn hảo trong kết cấu vỏ có sườn gia cường:
    Trong phân tích vỏ trụ và vỏ cầu tựa đơn chịu tải phân bố đều ($a = 32, R = 96, h = 0.8$), SQ4C cho kết quả độ võng $w$ tại điểm giữa vỏ hội tụ đơn điệu từ cận dưới và tiệm cận chính xác nghiệm giải tích của Reddy [119] và nghiệm phần tử vỏ 8 nút HBQ8 của Darilmaz & Kumbasar [6] (đạt $w = 0.3138$ tại lưới $16 \times 16$), trong khi số bậc tự do yêu cầu giảm hơn 50%.

  4. Độ chính xác vượt trội của phần tử Chebyshev SQ4P:
    Bằng cách tăng bậc đa thức $p_1 = p_2 = 3$, SQ4P đạt sai số tương đối nhỏ hơn $0.05%$ ngay cả ở lưới cực thô $6 \times 6$. Đối với bài toán tấm mỏng chịu uốn với $a/h = 10000$, trong khi các phần tử FEM thông thường bị khóa cắt hoàn toàn và suy biến nghiệm, SQ4P vẫn duy trì độ võng chuẩn xác tuyệt đối, khẳng định tính miễn nhiễm với hiện tượng khóa cắt trên mọi dải độ mảnh.

  5. Quy luật phân bố tối ưu cho vật liệu FGP-GPLs:
    Phân tích cơ học chi tiết trên kết cấu FGP-GPLs chứng minh rằng sự kết hợp giữa phân bố độ xốp đối xứng (P-S) và phân bố tiểu cầu graphene đối xứng (GPL-S) mang lại độ cứng uốn và tải trọng mất ổn định $P_{\text{cr}}^*$ cao nhất, đồng thời tăng tần số dao động tự nhiên mode 1 lên tới $35-42%$ so với phân bố đều (P-U) khi hàm lượng graphene chỉ chiếm $1.0\text{ wt}%$.

Implications đa chiều

  • Về mặt lý thuyết: Luận án đã chứng minh tính khả thi của việc giải quyết các bài toán HSDT phức tạp trên nền tảng phần tử $C^0$, mở ra hướng tiếp cận mới cho việc mô hình hóa các lý thuyết biến dạng cắt bậc cao hơn hoặc lý thuyết biến dạng lớp (Layer-wise theory) mà không làm phức tạp hóa cấu trúc phần tử.
  • Về mặt phương pháp luận: Khung tích hợp giữa làm trơn biến dạng trên miền con tế bào (CS-FEM) và đa thức trực giao Chebyshev cung cấp một quy trình mẫu mực có thể áp dụng cho các bài toán truyền nhiệt, điện - từ - cơ, và cơ học rạn nứt vật liệu.
  • Về mặt ứng dụng thực tiễn: Các ma trận phần tử SQ4H, SQ4C, SQ4P có thể tích hợp trực tiếp dưới dạng User-defined Element (UEL) vào các phần mềm thương mại phổ biến như ABAQUS, ANSYS hoặc NASTRAN, giúp các kỹ sư thiết kế tối ưu hóa kết cấu mỏng trong ngành công nghiệp hàng không và đóng tàu mà không cần đầu tư các siêu máy tính đắt tiền.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn chỉ rõ 4 hạn chế nội tại:

  1. Chi phí tính toán của SQ4T trong hình học phi phẳng: Việc thiết lập miền ảnh hưởng lân cận $nsp$ cho phần tử SQ4T trên các bề mặt vỏ có độ cong lớn hoặc độ cong phức tạp hai phương làm tăng đáng kể thời gian tìm kiếm lân cận và gây sai số số học do sự lệch trục tọa độ cục bộ giữa các phần tử lân cận.
  2. Số bậc tự do tăng nhanh trong SQ4P: Khi tăng bậc đa thức Chebyshev ($p_1, p_2 \ge 3$), số biến chuyển vị nội tại của phần tử tăng lên, dẫn đến ma trận độ cứng có băng thông lớn hơn, đòi hỏi kỹ thuật khử Gauss ngưng tụ tĩnh (static condensation) để tối ưu hóa bộ nhớ.
  3. Phạm vi lý thuyết đơn lớp tương đương (ESL): Luận án tập trung chủ yếu vào ESL, do đó mặc dù mô tả xuất sắc ứng xử tổng thể của tấm/vỏ composite nhiều lớp, mô hình chưa thể nắm bắt chi tiết hiện tượng trượt tách lớp (delamination) hoặc ứng suất cắt cục bộ từng lớp sâu như lý thuyết Layer-wise (LW).
  4. Giới hạn đàn hồi tuyến tính của vật liệu: Các mô hình hiện tại chủ yếu khảo sát phi tuyến hình học với giả thiết vật liệu đàn hồi tuyến tính, chưa tính đến phi tuyến vật liệu (elasto-plasticity) hoặc hư hỏng tích lũy do mỏi cơ nhiệt.

Chương trình nghiên cứu tương lai (5 hướng phát triển cụ thể):

  • Hướng 1: Mở rộng phần tử SQ4H và SQ4P cho phân tích tương tác đa trường vật lý: cơ - nhiệt - áp điện (piezoelectric) và tương tác lỏng - rắn (Fluid-Structure Interaction - FSI).
  • Hướng 2: Nghiên cứu ứng xử động phi tuyến và đáp ứng va đập tốc độ cao của tấm/vỏ FGP-GPLs trong môi trường nhiệt độ cao.
  • Hướng 3: Ứng dụng lý thuyết đàn hồi phi cục bộ (Nonlocal Elasticity Theory) kết hợp phần tử SQ4P để mô phỏng cấu trúc vi mô nano (nanoplates/nanoshells).
  • Hướng 4: Tích hợp mô hình cơ học phá hủy và dẻo hóa đàn hồi vào phần tử SQ4C để dự đoán sự phá hủy của kết cấu composite có sườn gia cường.
  • Hướng 5: Phát triển kỹ thuật chia lưới thích ứng thông minh (Adaptive mesh refinement) tự động tối ưu hóa bậc $p$ của Chebyshev và kích thước $h$ của tế bào trơn.

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

Nghiên cứu của tác giả Tôn Thất Hoàng Lân tạo ra sức ảnh hưởng sâu rộng trên nhiều bình diện:

  • Tác động học thuật quốc tế: Công trình đã công bố hàng loạt bài báo trên các tạp chí ISI uy tín thuộc nhóm Q1/Q2 hàng đầu thế giới trong lĩnh vực cơ học tính toán như Journal of Sandwich Structures & Materials, International Journal of Computational Methods, Periodica Polytechnica Civil Engineering, và Archive of Applied Mechanics. Những công bố này đóng góp vào kho tàng tri thức cơ học chất rắn của Việt Nam trên bản đồ khoa học quốc tế.
  • Chuyển đổi công nghiệp chế tạo tiên tiến: Cung cấp công cụ tính toán chính xác tuyệt đối cho việc thiết kế các kết cấu nhẹ có độ bền siêu cao sử dụng vật liệu xốp gia cường graphene (FGP-GPLs), hỗ trợ trực tiếp ngành công nghiệp chế tạo xe điện, thân vỏ máy bay không người lái (UAV), và vệ tinh không gian.
  • Giá trị kinh tế - xã hội: Giúp giảm trọng lượng kết cấu công trình chịu lực từ $15-25%$ nhờ tối ưu hóa phân bố vật liệu chức năng FGM, từ đó tiết kiệm chi phí nguyên vật liệu và giảm thiểu phát thải carbon trong giao thông vận tải.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ kỹ thuật/Cơ học kết cấu: Tiếp cận được phương pháp luận cải biên lý thuyết $C^0$-HSDT và kỹ thuật làm trơn biến dạng tế bào tường minh để phát triển các đề tài nghiên cứu chuyên sâu.
  • Các nhà phát triển phần mềm CAE (Computer-Aided Engineering): Sở hữu thuật toán phần tử tứ giác 4 nút hiệu năng cao, dễ dàng lập trình và tích hợp vào các lõi tính toán thương mại nhằm nâng cao năng lực cạnh tranh của phần mềm mô phỏng.
  • Kỹ sư kết cấu hàng không, cơ khí và công trình biển: Ứng dụng các bảng dữ liệu chuẩn và công thức phần tử để tính toán kiểm tra độ an toàn chịu lực, ổn định và rung động của các vỏ tàu, cánh máy bay, và bồn bể chịu áp lực cao.
  • Các nhà hoạch định chính sách và Viện nghiên cứu vật liệu mới: Có cơ sở khoa học định lượng vững chắc để ban hành các tiêu chuẩn thiết kế và ứng dụng vật liệu composite nano graphene trong công nghiệp quốc phòng và dân dụng.

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

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

Đóng góp độc đáo nhất là việc phát triển thành công phần tử SQ4H dựa trên việc chuyển đổi lý thuyết biến dạng cắt bậc ba của Reddy (TSDT/HSDT) thành dạng $C^0$-HSDT kết hợp kỹ thuật trơn biến dạng tế bào (CS-FEM). Công trình đã giải phóng lý thuyết HSDT khỏi sự ràng buộc của hàm xấp xỉ liên tục $C^1$, loại bỏ hoàn toàn hệ số hiệu chỉnh cắt $k_s$, đồng thời giải quyết xuất sắc bài toán phi tuyến hình học lớn mà không gặp bất kỳ hiện tượng suy biến số học nào.

2. Sự đổi mới về phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đây?

So với phần tử vỏ lai 8 nút HBQ8 của Darilmaz và Kumbasar (2006) hay phần tử vỏ cong 8 nút tích phân giảm KUMBA của Kumbasar et al. (2000), các phần tử đề xuất (đặc biệt là SQ4C và SQ4P) chỉ sử dụng cấu trúc 4 nút đơn giản, giúp giảm hơn $50%$ số bậc tự do toàn cục nhưng vẫn đạt độ chính xác tương đương hoặc vượt trội, triệt tiêu hoàn toàn hiện tượng khóa màng và hiện tượng dao động không năng lượng (hourglass). So với phần tử MITC4 của Bathe và Dvorkin (1985), phần tử SQ4T cải thiện vượt bậc tính liên tục của trường ứng suất qua biên nhờ chiến lược nội suy kép.

3. Phát hiện bất ngờ nhất trong quá trình tính toán số là gì?

Phát hiện bất ngờ nhất là ứng xử số của phần tử SQ4T trong bài toán vỏ cầu cong hai phương (Chương 7, Bảng 7.5 và Bảng 7.6): Mặc dù SQ4T thể hiện sự vượt trội về độ liên tục ứng suất trong bài toán phẳng, việc mở rộng vùng ảnh hưởng nút $nsp$ trên các mặt cong 3D làm tăng thời gian tính toán lên tới $19.64\text{ s}$ (gấp hơn 15 lần so với các phần tử khác cùng số bậc tự do $\approx 1600$ DOFs) và tạo ra hiện tượng cứng hóa nhẹ làm sai số hội tụ từ cận dưới. Đây là đóng góp phản biện khoa học trung thực rất giá trị của luận án.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Có. Luận án trình bày chi tiết toàn bộ các phương trình biến phân, ma trận liên hệ biến dạng - chuyển vị $\mathbf{B}$, ma trận đàn hồi $\mathbf{D}$, quy tắc tích phân đường trên biên miền con tứ giác, tọa độ các điểm buộc (tied points), cùng hệ thống thông số vật liệu và hình học chuẩn xác trong từng bảng số liệu, cho phép các nhà nghiên cứu độc lập lập trình tái lập chính xác $100%$ kết quả.

5. Kế hoạch nghiên cứu 10 năm được định hình như thế nào?

Lộ trình 10 năm tập trung vào việc mở rộng hệ phần tử SQ4 sang: (1) Mô phỏng cơ nhiệt vi mô đa tỷ lệ cho vật liệu tự phục hồi (self-healing nanocomposites), (2) Tích hợp thuật toán học máy (Machine Learning-assisted FEM) để tự động điều chỉnh tối ưu bậc đa thức Chebyshev $p$ và mật độ miền con $nc$, và (3) Xây dựng gói phần mềm tính toán kết cấu vỏ composite mỏng chuyên dụng mã nguồn mở phục vụ ngành công nghiệp chế tạo vệ tinh và phương tiện bay siêu thanh.

Kết luận

Luận án tiến sĩ của tác giả Tôn Thất Hoàng Lân đã đúc kết 6 đóng góp khoa học mang tính đột phá:

  1. Thiết lập thành công phần tử SQ4H: Kết hợp hoàn hảo kỹ thuật $C^0$-HSDT với làm trơn biến dạng tế bào CS-FEM, giải quyết triệt để bài toán phi tuyến hình học tấm phẳng và tấm gấp mà không cần hệ số hiệu chỉnh cắt.
  2. Sáng tạo phần tử SQ4T: Ứng dụng kỹ thuật nội suy kép (TIS) cho kết cấu tấm/vỏ, nâng cao tính liên tục của trường ứng suất và biến dạng qua biên phần tử.
  3. Phát triển phần tử SQ4C: Tích hợp kỹ thuật tổ hợp biến dạng màng, uốn, cắt, triệt tiêu hoàn toàn hiện tượng khóa màng cho kết cấu vỏ mỏng và vỏ có sườn gia cường.
  4. Tiên phong xây dựng phần tử Chebyshev SQ4P: Kết hợp lưới phần tử 4 nút với chuỗi đa thức trực giao Chebyshev, mang lại khả năng hội tụ siêu việt cho vật liệu xốp phân cấp chức năng gia cường graphene (FGP-GPLs).
  5. Đồng nhất hóa khung tính toán: Toàn bộ 4 phần tử được xây dựng trên nền tảng lý thuyết đơn lớp tương đương (ESL), cho phép chuyển đổi linh hoạt giữa các hệ vật liệu từ đẳng hướng, composite nhiều lớp, FGM đến FGP-GPLs.
  6. Đánh giá sai số toàn diện và khách quan: Thiết lập bộ dữ liệu so sánh chuẩn mực về độ chính xác, tốc độ hội tụ và chi phí thời gian tính toán giữa 4 phần tử đề xuất với các nghiệm giải tích và phần tử quốc tế uy tín.

Công trình đã khẳng định sự trưởng thành vượt bậc của nền cơ kỹ thuật tính toán Việt Nam, tạo tiền đề vững chắc cho việc ứng dụng các phương pháp số tiên tiến vào thực tiễn thiết kế công trình công nghệ cao trong kỷ nguyên mới. Trích dẫn nguyên văn khẳng định cuối cùng của tác giả tại trang 121: "Dựa vào những kết quả số đã đạt được... các phần tử được xây dựng và phát triển trong luận án là các phần tử tứ giác đơn giản và hiệu quả trong phân tích kết cấu dạng tấm/vỏ so với các loại phần tử khác đang hiện hành."