Tổng quan về luận án

Nghiên cứu của tác giả Nguyễn Đình Dư thuộc chuyên ngành Cơ kỹ thuật (Mã số: 9520101.01) tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của GS.TSKH Nguyễn Đình Đức và PGS.TS Bùi Quốc Tính, đại diện cho một bước đột phá trong lĩnh vực cơ học tính toán (computational mechanics) và mô phỏng số môi trường liên tục. Kể từ khi Phương pháp Phần tử hữu hạn (Finite Element Method - FEM) được khởi xướng bởi Courant (1940) và hệ thống hóa toán học chặt chẽ qua công trình kinh điển của Gilbert Strang và George Fix (1973), FEM đã trở thành trụ cột mô phỏng kết cấu then chốt. Tuy nhiên, FEM tiêu chuẩn vẫn bộc lộ các khiếm khuyết nội tại nghiêm trọng: tính gián đoạn vật lý của trường đạo hàm (gradient trường ứng suất và biến dạng) tại các nút phần tử, độ nhạy cao dẫn tới suy giảm độ chính xác khi lưới bị méo mó, cùng hiện tượng "khóa số" (locking phenomena) điển hình như khóa thể tích (volumetric locking) khi hệ số Poisson tiệm cận giới hạn không nén được ($\nu \to 0.5$) và khóa cắt (shear locking) trong kết cấu tấm/vỏ mỏng.

Khoảng trống nghiên cứu (research gap) mà luận án giải quyết tập trung vào việc khắc phục sự đánh đổi giữa độ chính xác và chi phí tính toán trong phân tích kết cấu phức tạp. Trong khi các giải pháp khắc phục truyền thống như nâng cao bậc phần tử (p-refinement với phần tử Q8, HH20), phương pháp đẳng hình học IGA (Hughes et al., 2005) hay phương pháp không lưới (Meshfree methods như EFG, RKPM) đòi hỏi thuật toán xử lý biên phức tạp do đánh mất thuộc tính Kronecker-delta hoặc làm phình to số bậc tự do (Degrees of Freedom - DOF), luận án phát triển toàn diện Phương pháp Phần tử hữu hạn nội suy kép (Consecutive Interpolation Finite Element Method - CFEM) dựa trên Thủ tục nội suy liên tiếp (Consecutive Interpolation Procedure - CIP).

Mục tiêu và câu hỏi nghiên cứu của luận án được cụ thể hóa qua ba trục chính:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập công thức hàm dạng CIP tổng quát cho các phần tử 2D (tam giác 3 nút CT3, tứ giác 4 nút CQ4) và 3D (tứ diện 4 nút CTH4, lục diện 8 nút CHH8) nhằm đảm bảo trường ứng suất – biến dạng liên tục, trơn mịn mà vẫn bảo toàn thuộc tính Kronecker-delta và không làm tăng số bậc tự do tổng thể?
  2. Câu hỏi nghiên cứu 2 (RQ2): Phương pháp phần tử hữu hạn mở rộng nội suy kép (XCFEM/XCQ4) phát huy hiệu quả ra sao trong việc tính toán hệ số cường độ ứng suất tĩnh (SIFs), động (DSIFs) và mô phỏng quỹ đạo lan truyền vết nứt trong vật liệu biến tính chức năng (Functionally Graded Materials - FGM)?
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để giải quyết triệt để bài toán phi tuyến hình học 2D/3D và hiện tượng khóa thể tích đối với vật liệu gần như không nén được thông qua việc kết hợp CIP với các mô hình tích phân số mới thay thế tích phân Gauss truyền thống?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết H1: Hàm dạng CIP thông qua toán tử gradient trung bình nút có trọng số hình học sẽ làm mượt trường gradient và tăng tốc độ hội tụ năng lượng biến dạng vượt trội so với phần tử Lagrange cùng bậc.
  • Giả thuyết H2: Sự kết hợp giữa hàm làm giàu Heaviside, hàm dốc (ramp function) và phần tử CQ4 sẽ loại bỏ hiện tượng kỳ dị giả tạo, nắm bắt chính xác hiện tượng khúc xạ đường nứt trong môi trường FGM bất đồng nhất.
  • Giả thuyết H3: Mô hình tích phân điểm giữa phần tử 3D (3D-EM) và mặt giữa phần tử (EF) sẽ giảm thiểu số điểm tích phân từ 8 điểm Gauss xuống còn 1 điểm hoặc các điểm mặt, loại bỏ hiện tượng khóa thể tích và tiết kiệm đáng kể thời gian xử lý phi tuyến lớn.

Về quy mô và phạm vi, luận án khảo sát toàn diện từ các kết cấu mẫu 2D/3D đàn hồi tuyến tính, vật liệu composite FGM với hàm biến thiên cơ tính dạng mũ và tuyến tính (như hợp kim TiB/Ti), bài toán dao động riêng của hình trụ rỗng lệch tâm, đến các cấu hình phá hủy phức tạp (mẫu uốn 3 điểm TPB, uốn 4 điểm FPB, tấm chữ nhật có lỗ tròn và vết nứt nghiêng) và các bài toán phi tuyến hình học lớn (màng Cook 3D, kết cấu nhà cao tầng chịu tải xô ngang).

Literature Review và Positioning

Bối cảnh lý thuyết của cơ học tính toán đương đại ghi nhận ba dòng chảy chính trong nỗ lực nâng cao hiệu năng xấp xỉ số:

  1. Dòng nghiên cứu phương pháp đẳng hình học (Isogeometric Analysis - IGA): Khởi xướng bởi Hughes et al. (2005), sử dụng hàm B-spline và NURBS để đồng nhất hóa hình học CAD và giải tích phân tích. Mặc dù kiểm soát hoàn hảo độ trơn $C^p$ giữa các phần tử, IGA gặp rào cản lớn khi áp đặt điều kiện biên thiết yếu do thiếu thuộc tính Kronecker-delta, buộc phải sử dụng phương pháp nhân tử Lagrange hoặc hàm phạt phức tạp.
  2. Dòng nghiên cứu phương pháp không lưới (Meshfree/Meshless methods): Tiêu biểu là Element-Free Galerkin (Belytschko et al., 1994) và Reproducing Kernel Particle Method (Liu et al., 1995). Ưu điểm loại bỏ hoàn toàn mắt lưới nhưng lại phụ thuộc nặng nề vào các tham số bán kính miền hỗ trợ (support domain) mang tính kinh nghiệm, chi phí tính toán tích phân nền rất cao và cũng không thỏa mãn Kronecker-delta.
  3. Dòng nghiên cứu phần tử hữu hạn làm trơn (Smoothed FEM - S-FEM): Phát triển mạnh bởi Liu et al. (2007-2015) với các biến thể CS-FEM (cell-based), NS-FEM (node-based), ES-FEM (edge-based) và FS-FEM (face-based). Dù cải thiện độ chính xác gradient, S-FEM can thiệp vào trường biến dạng thông qua các ô làm mịn mở rộng (smoothing cells), làm tăng độ phức tạp của thuật toán xây dựng ma trận độ cứng.
       ┌─────────────────────────────────────────────────────────────┐
       │                Phương pháp số truyền thống (FEM)            │
       │           - C0 liên tục, gradient gián đoạn tại nút          │
       │           - Khóa thể tích, suy giảm bậc khi lưới méo        │
       └──────────────────────────────┬──────────────────────────────┘
                                      │
              ┌───────────────────────┴───────────────────────┐
              ▼                                               ▼
┌───────────────────────────┐                   ┌───────────────────────────┐
│ Phương pháp thay thế lưới │                   │ Cải tiến trên nền tảng    │
│ - IGA (Hughes, 2005)      │                   │ Phần tử hữu hạn (FEM)     │
│ - Meshfree (EFG, RKPM)    │                   │ - S-FEM (Liu et al., 2007)│
│ *Nhược điểm: Mất tính     │                   │ *Nhược điểm: Tạo thêm ô   │
│ Kronecker-delta, khó biên │                   │ làm mịn trung gian        │
└───────────────────────────┘                   └─────────────┬─────────────┘
                                                              │
                                                              ▼
                                                ┌───────────────────────────┐
                                                │ Phương pháp CFEM (Luận án)│
                                                │ - Thủ tục CIP toán học    │
                                                │ - Giữ nguyên bậc tự do DOF│
                                                │ - Bảo toàn Kronecker-delta│
                                                │ - Gradient liên tục, trơn │
                                                └───────────────────────────┘

Trong tương quan đó, thủ tục nội suy kép CIP được Zheng et al. (2010) công bố đầu tiên cho phần tử tam giác 2D (TFEM), sau đó được mở rộng phân tích nứt đàn hồi bởi Peng et al. (2017). Nhóm nghiên cứu của PGS.TS Bùi Quốc Tính (2014, 2015, 2016, 2017) đã mở rộng thành công cho phần tử tứ giác 4 nút CQ4, dao động vật liệu áp điện và cơ học phá hủy tĩnh/động.

Luận án này định vị tại mắt xích đột phá then chốt: Hệ thống hóa công thức toán học CIP thành cấu trúc tường minh tổng quát cho không gian 3 chiều (CTH4, CHH8), tích hợp toàn diện vào cơ học nứt FGM bất đẳng hướng dưới tải trọng động lực học cao tần, và tiên phong giải quyết bài toán phi tuyến hình học 2D/3D bằng sơ đồ tích phân số mới "Element Mid-point" (3D-EM) và "Element Mid-face" (EF).

So sánh với các nghiên cứu quốc tế tiêu chuẩn:

  • Đối chiếu với bài toán thanh biến tính FGM chịu kéo của Kim & Paulino (2002), luận án chứng minh phần tử tuyến tính CQ4 đạt độ chính xác tương đương phần tử bậc cao 8 nút Q8 nhưng tối ưu hóa vượt bậc về tài nguyên tính toán.
  • So sánh với mô hình XFEM cổ điển của Moës et al. (1999) trong cơ học phá hủy, việc tích hợp CIP vào phần tử XCQ4 triệt tiêu hoàn toàn sự nhảy vọt ứng suất cục bộ quanh đỉnh nứt, đem lại đường cong DSIFs ổn định theo thời gian mà không cần làm mịn lưới cực tiểu.

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

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

Luận án đóng góp trực tiếp vào lý thuyết đàn hồi vi phân và giải tích xấp xỉ phần tử hữu hạn thông qua việc thiết lập giải pháp toán học hoàn chỉnh cho trường chuyển vị nội suy kép: $$u(\mathbf{x}) = \sum_{I=1}^{n} R_I(\mathbf{x}) a_I = \mathbf{R} \mathbf{a}$$

Trong đó, vectơ hàm dạng CIP $R_I(\mathbf{x})$ được cấu thành từ toán tử hàm dạng Lagrange cổ điển $N_I$ kết hợp các đạo hàm trung bình nút $\bar{N}{I,x}, \bar{N}{I,y}, \bar{N}{I,z}$ cùng các hàm phụ trợ $\phi_i, \phi{ix}, \phi_{iy}, \phi_{iz}$: $$R_I(\mathbf{x}) = \sum_{i=1}^{n_e} \left( \phi_i(\mathbf{x}) N_I^{(e)} + \phi_{ix}(\mathbf{x}) \bar{N}{I,x}^{(e)} + \phi{iy}(\mathbf{x}) \bar{N}{I,y}^{(e)} + \phi{iz}(\mathbf{x}) \bar{N}_{I,z}^{(e)} \right)$$

Đạo hàm trung bình nút được tính toán thông qua trọng số thể tích/diện tích của tập hợp các phần tử lân cận chia sẻ chung nút $I$ ($S_I$): $$\bar{N}{I,w}^{(e)} = \sum{e \in S_I} w_e \left. \frac{\partial N_I}{\partial w} \right|e, \quad w_e = \frac{A_e}{\sum{j \in S_I} A_j}$$

Chứng minh toán học trong luận án đã khẳng định tính thỏa mãn tuyệt đối của hai bổ đề quan trọng:

  1. Bổ đề thỏa mãn tính Kronecker-delta: $\phi_i(\mathbf{x}p) = \delta{ip}$ tại mọi nút hình học cục bộ $p \in {i, j, k, m}$.
  2. Bổ đề triệt tiêu đạo hàm nút: $\phi_{i,w}(\mathbf{x}_p) = 0$ đối với mọi hướng tọa độ không gian $w \in {x, y, z}$.

Nhờ cấu trúc toán học này, luận án chứng minh rằng: "Thủ tục CIP can thiệp trực tiếp vào hàm dạng của phần tử hữu hạn truyền thống bằng các kỹ thuật liên quan đến đạo hàm trung bình của nút phần tử... không làm tăng số DOF nhưng vẫn đem lại hiệu quả cao trong khi FEM truyền thống có bậc thấp không mang lại được trừ khi dùng phần tử bậc cao."

                       ┌───────────────────────────────────────┐
                       │  Nút hình học I & Miền hỗ trợ S_I     │
                       │  - Tập hợp các phần tử con lân cận    │
                       └───────────────────┬───────────────────┘
                                           │
                    ┌──────────────────────┴──────────────────────┐
                    ▼                                             ▼
     ┌─────────────────────────────┐               ┌─────────────────────────────┐
     │ Đạo hàm trung bình nút      │               │ Hàm phụ trợ CIP (Phi_i)     │
     │ N_bar_{I,w} = sum(w_e * dN) │               │ - Thỏa mãn Kronecker-delta  │
     │ w_e = A_e / sum(A_j)        │               │ - Triệt tiêu đạo hàm nút    │
     └──────────────┬──────────────┘               └──────────────┬──────────────┘
                    │                                             │
                    └──────────────────────┬──────────────────────┘
                                           │
                                           ▼
                       ┌───────────────────────────────────────┐
                       │      Hàm dạng CIP bậc cao R_I(x)      │
                       │ - Tăng bậc trơn C1 nội bộ phần tử     │
                       │ - DOF không đổi = 2 (2D) hoặc 3 (3D)  │
                       │ - Khử hoàn toàn gián đoạn gradient    │
                       └───────────────────────────────────────┘

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

Khung phân tích của luận án là sự hợp nhất liên ngành chặt chẽ giữa ba trường phái lý thuyết:

  • Lý thuyết cơ học rạn nứt đàn hồi tuyến tính bất đồng nhất (Inhomogeneous LEFM): Mô tả vật liệu FGM có cơ tính phụ thuộc tọa độ $E(\mathbf{x}) = E_0 e^{\beta y}$, $\nu(\mathbf{x}) = \nu_0 e^{\kappa y}$.
  • Kỹ thuật làm giàu trường chuyển vị (XFEM Enrichment): Tích hợp hàm gián đoạn Heaviside cho thân vết nứt và 4 hàm tiệm cận đỉnh nứt không đồng nhất kết hợp hàm dốc (ramp function) để loại bỏ miền pha trộn không tương thích.
  • Lý thuyết phi tuyến hình học biến dạng hữu hạn (Finite Strain Formulation): Ứng dụng tensor biến dạng Green-Lagrange và tensor ứng suất Piola-Kirchhoff loại 2, kết hợp kỹ thuật tích phân số 3D-EM và EF nhằm triệt tiêu hoàn toàn hiện tượng khóa thể tích khi $\nu \to 0.5$.

Điều kiện biên và giới hạn khung phân tích được xác định rõ: Vật liệu đàn hồi tuyến tính tuân theo định luật Hooke tổng quát, hình học biến dạng lớn nhưng biến dạng vi phân nằm trong miền đàn hồi, bài toán nứt phẳng 2D và bài toán kết cấu khối 3D liên tụ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ủ nghiêm ngặt quan điểm nhận thức luận thực chứng tính toán (computational positivism), xây dựng quy trình nghiên cứu mô phỏng tất định kết hợp thẩm định đối chéo (cross-validation) đa tầng.

Cấp độ phân tích Đối tượng kết cấu Phương pháp phần tử Sơ đồ tích phân / Đánh giá
Tĩnh học 2D/3D Dầm công son chịu uốn, Dầm chữ T CT3, CQ4, CTH4, CHH8 Tích phân Gauss chuẩn; Đối chiếu nghiệm giải tích Timoshenko
Cơ tính FGM Thanh phức tạp TiB/Ti chịu kéo CQ4 So sánh chuẩn phần tử bậc cao Q8 (Kim & Paulino)
Động lực học nứt Tấm FGM có lỗ tròn nứt nghiêng (CCTT) XCQ4 Tích phân tương tác động $J$, tiêu chuẩn góc nứt MTS
Lan truyền nứt Dầm FGM uốn 3 điểm (TPB), 4 điểm (FPB) XCQ4 Kỹ thuật hàm dốc (ramp function) đối chiếu thực nghiệm
Phi tuyến hình học Màng Cook 3D, Nhà cao tầng xô ngang CHH8 Sơ đồ tích phân mới 3D-EM và EF; Arc-length algorithm

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thiết kế khép kín nhằm bảo đảm tính lặp lại (reproducibility) và kiểm chứng sai số tối đa:

  1. Rời rạc hóa miền tính toán: Tạo lập cấu trúc lưới đa cấp từ thô đến cực mịn (ví dụ dầm công son: từ $12\times3$ đến $192\times48$; khối 3D chữ T: 480, 960, 2240 và 6720 phần tử).
  2. Xây dựng ma trận độ cứng cải tiến: Tự động xác định các tập lân cận $S_I$ cho từng nút, tính toán ma trận Jacobi chuyển đổi $\mathbf{J}$ qua đạo hàm hàm dạng CIP: $$\mathbf{J} = \begin{bmatrix} \frac{\partial x}{\partial \xi} & \frac{\partial y}{\partial \xi} & \frac{\partial z}{\partial \xi} \ \frac{\partial x}{\partial \eta} & \frac{\partial y}{\partial \eta} & \frac{\partial z}{\partial \eta} \ \frac{\partial x}{\partial \zeta} & \frac{\partial y}{\partial \zeta} & \frac{\partial z}{\partial \zeta} \end{bmatrix}$$
  3. Định lượng phá hủy nứt gãy: Sử dụng tích phân tương tác miền $J$ (Interaction Integral) kết hợp hàm trọng số trơn $q(\mathbf{x})$ để phân tách chính xác hệ số cường độ ứng suất Mode I ($K_I$) và Mode II ($K_{II}$). Quỹ đạo phát triển vết nứt được điều khiển bởi tiêu chuẩn ứng suất tiếp tuyến cực đại (Maximum Tangential Stress - MTS): $$\theta_c = 2 \arctan \left( \frac{K_I - \sqrt{K_I^2 + 8 K_{II}^2}}{4 K_{II}} \right)$$
  4. Thuật toán phi tuyến: Sử dụng phương pháp Newton-Raphson kết hợp kỹ thuật kiểm soát độ dài cung (Arc-length method) để bám bắt chính xác các điểm tới hạn (limit points) và hiện tượng nhảy xuyên biến dạng (snap-through instability).

Data và phân tích

Toàn bộ thuật toán được tác giả lập trình độc lập 100% trên nền tảng MATLAB. Dữ liệu số được đối chiếu chuẩn hóa với phần mềm thương mại cao cấp ABAQUS sử dụng lưới phần tử 20 nút bậc cao (HH20) với hơn 10.800 phần tử ($>130.000$ bậc tự do).

Sai số năng lượng biến dạng toàn cục được định lượng chính xác qua chuẩn $L_2$: $$\text{Error}(E) = \sqrt{\frac{1}{2} \int_{\Omega} (\boldsymbol{\sigma} - \boldsymbol{\sigma}{\text{exact}}) : (\boldsymbol{\varepsilon} - \boldsymbol{\varepsilon}{\text{exact}}) , d\Omega}$$

Các chỉ số thống kê về thời gian tính toán, sai số tương đối (%) và tốc độ hội tụ (convergence rate) được biểu diễn trên hệ trục logarit kép ($\log-\log$) để minh chứng tính ổn định tiệm cận của thuật toán.

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

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

Thứ nhất: Tốc độ hội tụ và độ chính xác gradient ứng suất vượt trội. Trong bài toán dầm công son 2D chịu uốn, dữ liệu thực nghiệm số chỉ rõ: Để đạt cùng ngưỡng sai số năng lượng $0.3%$, phần tử CIP tứ giác (CQ4) chỉ cần $2.68,\text{s}$ trên lưới $24\times6$, trong khi phần tử FEM tiêu chuẩn (Q4) phải cần tới $5.58,\text{s}$ trên lưới mịn $96\times24$. Phương pháp đề xuất tiết kiệm hơn $52%$ chi phí thời gian tính toán thực tế.

Cấu hình lưới Số bậc tự do (DOFs) Sai số năng lượng Q4 (%) Thời gian Q4 (s) Sai số năng lượng CQ4 (%) Thời gian CQ4 (s)
$12\times3$ 104 5.21 0.072 1.14 1.004
$24\times6$ 350 2.68 0.215 0.29 2.680
$48\times12$ 1.274 1.35 0.980 0.08 8.420
$96\times24$ 4.850 0.31 5.580 0.02 28.150

Thứ hai: Khả năng mô phỏng liên tục trường ứng suất FGM và tái phân bố vị trí phá hủy. Đối với thanh kết cấu phức tạp FGM TiB/Ti ($E_0 = 375,\text{GPa}, \nu_0 = 0.14$ tại pha TiB và $E_1 = 107,\text{GPa}, \nu_1 = 0.34$ tại pha Ti), dữ liệu luận án chứng minh sai số ứng suất pháp $\sigma_{xx}$ tại điểm kỳ dị $B'$ của phần tử CQ4 chỉ đạt $0.07%$ so với chuẩn phần tử bậc cao Q8, trong khi phần tử Q4 tiêu chuẩn mắc sai số lên tới $1.63%$.

Đặc biệt, kết quả chỉ ra quy luật cơ học quan trọng: Sự phân cấp vật liệu FGM giúp tái phân bố trường ứng suất, làm giảm giá trị ứng suất đỉnh và dịch chuyển vị trí tập trung ứng suất nguy hiểm từ điểm $A'$ (vật liệu đồng nhất) ra biên ngoài điểm $B'$, tạo điều kiện lý tưởng cho công tác kiểm định và bảo trì kết cấu thực tế.

Vật liệu đồng nhất:     [Ứng suất đỉnh tại A' (Bên trong kết cấu)]  --> Nguy cơ phá hủy ngầm khó phát hiện
                                        │
                                        ▼ (Chuyển đổi sang vật liệu FGM)
Vật liệu biến tính FGM: [Ứng suất đỉnh dịch chuyển ra B' (Biên ngoài)] --> Dễ dàng kiểm tra và bảo trì

Thứ ba: Độ chính xác dao động riêng 3D trên kết cấu phức tạp. Trong phân tích dao động tự do của hình trụ 3D có lỗ khoét tròn lệch tâm, phần tử CHH8 thể hiện ưu thế vượt trội khi mô phỏng sự tách đôi tần số dao động do tính bất đối xứng hình học.

"Sai số lớn nhất quan sát được ở tần số tự nhiên khác không thứ tư trong trường hợp độ lệch tâm 80%, chỉ là 0.52% đối với 3072 phần tử CHH8. Giá trị tương ứng thu được của các phần tử HH8 sử dụng cùng một lưới là gần 1.76%."

Thứ tư: Độ tin cậy trong mô phỏng lan truyền nứt động và khúc xạ đường nứt. Phần tử XCQ4 kết hợp hàm dốc đã mô phỏng hoàn hảo quỹ đạo nứt gãy trong các mẫu dầm FGM chịu uốn 3 điểm (TPB) và 4 điểm (FPB). Kết quả đường nứt dự đoán trùng khớp hoàn toàn với dữ liệu thực nghiệm uốn phá hủy và nghiệm số từ phương pháp phần tử hữu hạn mở rộng rời rạc (XFEM-DET), chứng minh khả năng tự điều chỉnh hướng lan truyền khi vết nứt đi qua dải biến thiên mô đun đàn hồi.

Thứ năm: Khắc phục triệt để hiện tượng khóa thể tích trong phân tích phi tuyến. Sơ đồ tích phân 3D-EM và EF kết hợp phần tử CHH8 giải quyết thành công bài toán màng Cook 3D và tấm chịu uốn với vật liệu gần như không nén được ($\nu = 0.4999$). Trong khi tích phân Gauss tiêu chuẩn bị khóa cứng hoàn toàn và dự báo sai lệch độ võng đỉnh, sơ đồ tích phân mới cho kết quả chuyển vị tiệm cận chính xác nghiệm giải tích mà không phát sinh dao động số dạng hourglass.

Implications đa chiều

  • Về mặt học thuật lý thuyết: Xác lập nền tảng toán học vững chắc cho lớp phần tử hữu hạn có độ trơn gradient nội tại cao mà không phụ thuộc vào bậc đa thức Lagrange, mở đường cho việc tái cấu trúc các thư viện phần tử hữu hạn kinh điển.
  • Về phương pháp luận tính toán: Cung cấp sơ đồ tích phân 3D-EM và EF với số điểm tích phân tối thiểu, giảm tải tới $60-70%$ chi phí tính toán ma trận tiếp tuyến trong các bài toán phi tuyến lớn.
  • Về ứng dụng công nghiệp: Ứng dụng trực tiếp vào việc thiết kế tối ưu các cấu kiện hàng không vũ trụ (vỏ tên lửa, cánh tuabin FGM chịu nhiệt), tấm chắn bức xạ trong lò phản ứng hạt nhân và vật liệu sinh học y tế (xương nhân tạo gradient cơ tính).
  • Về chuyển giao công nghệ: Toàn bộ thuật toán mã nguồn mở trên MATLAB tạo tiền đề để các viện nghiên cứu và doanh nghiệp công nghiệp trong nước phát triển các phần mềm kỹ thuật hỗ trợ máy tính (CAE) tự chủ, giảm thiểu phụ thuộc vào các bộ phần mềm thương mại đắt đỏ của nước ngoài.

Limitations và Future Research

Luận án thừa nhận một cách khoa học 4 giới hạn nghiên cứu chính:

  1. Mô hình vật liệu: Mới tập trung vào trạng thái đàn hồi tuyến tính vi phân tuân theo định luật Hooke; chưa mở rộng sang các trường ứng xử phi tuyến vật liệu phức tạp như tính dẻo hữu hạn (large-strain elastoplasticity), tính nhớt đàn hồi hoặc hư hại vi mô liên tục (continuum damage mechanics).
  2. Phạm vi cơ học phá hủy: Các bài toán mô phỏng nứt động và lan truyền nứt mới được hoàn thiện trên miền phẳng 2D; chưa giải quyết trọn vẹn bài toán mặt nứt không gian 3D phức tạp với các hiện tượng xoắn nứt (crack twisting) và phân nhánh nứt đa trục.
  3. Chi phí thiết lập ma trận ban đầu: Do miền hỗ trợ nút $S_I$ trong CIP mở rộng hơn FEM truyền thống, thời gian tính toán hàm dạng và ghép nối ma trận phần tử ở bước tiền xử lý cao hơn FEM tiêu chuẩn trên cùng một cấp độ lưới.
  4. Tương tác đa trường vật lý: Chưa tích hợp đồng thời bài toán liên kết đa trường vật lý phức tạp (Multiphysics) như tương tác cơ nhiệt - điện từ - dòng chảy chất lỏng trong môi trường FGM nứt nẻ.

Chương trình nghiên cứu 10 năm tiếp theo được đề xuất:

  • Hướng 1: Mở rộng công thức XCFEM cho bài toán mặt nứt không gian 3D tổng quát trong vật liệu composite dị hướng.
  • Hướng 2: Phát triển phần tử CIP phi tuyến đàn - dẻo và từ biến nhiệt độ cao phục vụ công nghiệp luyện kim và nhiệt điện.
  • Hướng 3: Tối ưu hóa thuật toán song song hóa quy mô lớn trên nền tảng GPU (NVIDIA CUDA) cho sơ đồ tích phân 3D-EM.
  • Hướng 4: Tích hợp CFEM với mạng nơ-ron vật lý (Physics-Informed Neural Networks - PINNs) nhằm xây dựng mô hình đại diện thời gian thực (real-time digital twins) cho các kết cấu công trình chịu tải trọng bất thường.

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

  • Tác động học thuật: Luận án tạo ra bước tiến lớn trong chuyên ngành Cơ kỹ thuật tại Việt Nam, đóng góp nhiều công trình công bố trên các tạp chí ISI/Scopus Q1 uy tín như Acta Mechanica, Composite Structures. Tiềm năng trích dẫn quốc tế ước tính cao nhờ giải quyết đúng các điểm nghẽn của FEM truyền thống.
  • Đột phá công nghiệp: Cung cấp công cụ mô phỏng chuẩn xác độ bền mỏi và lan truyền nứt cho ngành công nghiệp chế tạo vật liệu tiên tiến FGM, tối ưu hóa quá trình sản xuất giảm thiểu khuyết tật tập trung ứng suất.
  • Chính sách và an toàn công trình: Cung cấp cơ sở khoa học định lượng phục vụ việc biên soạn các tiêu chuẩn thiết kế kết cấu cao tầng chịu tải trọng xô ngang bất thường (gió bão, động đất) và quy chuẩn kiểm định an toàn đập thủy điện, công trình ngầm.

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

  • Nghiên cứu sinh và Giới học thuật: Thừa hưởng hệ thống công thức giải tích chuẩn xác, tường minh về hàm dạng CIP 2D/3D và phương pháp XCFEM làm nền tảng phát triển các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư R&D Hàng không và Cơ khí chính xác: Sử dụng công cụ mô phỏng để phân tích ứng suất mịn và trường nhiệt trên các chi tiết FGM mà không cần đầu tư hệ thống máy chủ tính toán siêu mạnh để chia lưới cực mịn.
  • Chuyên gia thiết kế kết cấu công trình: Tiếp cận giải pháp phân tích phi tuyến hình học 3D ổn định, loại bỏ hoàn toàn rủi ro khóa số trong các kết cấu bê tông cốt thép hoặc vật liệu polyme gần như không nén được.

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? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc xây dựng thành công công thức toán học hàm dạng tổng quát CIP cho không gian 3 chiều ($R_I(\mathbf{x})$ cho CTH4 và CHH8). Công thức này mở rộng trực tiếp lý thuyết xấp xỉ phần tử hữu hạn Lagrange cổ điển bằng cách tích hợp toán tử gradient trung bình nút có trọng số hình học. Điểm đột phá nằm ở chỗ nó tạo ra trường xấp xỉ trơn bậc cao giữa các phần tử nhưng vẫn bảo toàn tuyệt đối thuộc tính Kronecker-delta ($\phi_i(\mathbf{x}p) = \delta{ip}$) và không làm tăng bất kỳ bậc tự do (DOF) nào của hệ phương trình tổng thể.

2. Đột phá phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó? Trả lời:

  • So với phương pháp S-FEM của Liu et al. (2007), CIP can thiệp trực tiếp vào cấu trúc giải tích của hàm dạng $R_I(\mathbf{x})$ thay vì tạo thêm các ô làm mịn trung gian (smoothing cells), giúp giảm thiểu độ phức tạp hình học và giữ nguyên cấu trúc ma trận thưa chuẩn của FEM.
  • So với phần tử bậc cao Q8 trong nghiên cứu FGM của Kim & Paulino (2002), phần tử tuyến tính CQ4 của luận án đạt độ chính xác tương đương (sai số $0.07%$ so với $1.63%$ của Q4) nhưng tiết kiệm hơn $50%$ số nút trên mỗi phần tử, giảm thiểu đáng kể băng thông ma trận độ cứng.

3. Phát hiện số bất ngờ nhất được hỗ trợ bởi dữ liệu thực nghiệm là gì? Trả lời: Phát hiện về sự dịch chuyển vị trí ứng suất nguy hiểm trong kết cấu FGM: Khi chuyển từ vật liệu đồng nhất sang FGM biến thiên dạng mũ (thanh TiB/Ti), vị trí ứng suất pháp cực đại $\sigma_{xx}$ chuyển đổi hoàn toàn từ góc lõm bên trong ($A'$) ra mép biên tự do bên ngoài ($B'$), đồng thời độ lớn ứng suất đỉnh giảm rõ rệt. Dữ liệu này chứng minh cơ chế tự tái phân phối ứng suất của FGM, mang lại lợi thế an toàn kết cấu cực lớn trong kỹ thuật thực tế.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Trả lời: Luận án cung cấp đầy đủ và chi tiết toàn bộ công thức giải tích của các hàm phụ trợ ($\phi_i, \phi_{ix}, \phi_{iy}, \phi_{iz}$), bảng trọng số và tọa độ điểm tích phân cho mô hình 3D-EM và EF, cùng thuật toán lặp phi tuyến Newton-Raphson/Arc-length. Toàn bộ mã nguồn được chuẩn hóa trên MATLAB cho phép cộng đồng nghiên cứu tái lập chính xác $100%$ các bảng số liệu và đồ thị hội tụ.

5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Chương trình 10 năm định hình lộ trình: Giai đoạn 1 (1-3 năm): Mở rộng XCFEM mô phỏng nứt mỏi 3D không gian; Giai đoạn 2 (3-6 năm): Xây dựng mô hình vật liệu đàn - dẻo phi tuyến và từ biến nhiệt độ cao; Giai đoạn 3 (6-10 năm): Đóng gói thư viện phần tử CIP tối ưu hóa song song trên GPU CUDA và tích hợp vào các phần mềm CAE mã nguồn mở quốc tế.

Kết luận

  1. Luận án đã thiết lập thành công cơ sở toán học tường minh cho Phương pháp Phần tử hữu hạn nội suy kép (CFEM) trên toàn bộ các họ phần tử 2D (CT3, CQ4) và 3D (CTH4, CHH8), duy trì trọn vẹn thuộc tính Kronecker-delta mà không làm tăng số bậc tự do của bài toán.
  2. Phát triển phương pháp phần tử hữu hạn mở rộng nội suy kép (XCFEM/XCQ4) kết hợp hàm dốc, giải quyết triệt để bài toán tập trung ứng suất, dao động động lực học và mô phỏng chính xác đường lan truyền vết nứt trong vật liệu biến tính chức năng FGM.
  3. Đề xuất thành công hai sơ đồ tích phân số mới 3D-EM và EF, khắc phục hoàn toàn hiện tượng khóa thể tích trong bài toán phi tuyến hình học lớn đối với vật liệu gần như không nén được ($\nu \to 0.5$) với chi phí thời gian tính toán tối ưu.
  4. Minh chứng tính ưu việt vượt bậc về độ chính xác và tốc độ hội tụ qua hệ thống kiểm chứng số đa dạng: Giảm $52%$ thời gian tính toán trên dầm 2D để đạt cùng ngưỡng sai số $0.3%$; giảm sai số ứng suất FGM xuống $0.07%$; đạt sai số tần số dao động riêng 3D chỉ $0.52%$.
  5. Mở ra ba hướng nghiên cứu học thuật mới đầy tiềm năng: Cơ học phá hủy 3D trong môi trường bất đồng nhất, bài toán tương tác đa trường vật lý phi tuyến và công nghệ tính toán song song hiệu năng cao cho phần tử nội suy kép.
  6. Đóng góp di sản khoa học thực tiễn vững chắc, nâng cao năng lực tự chủ công nghệ mô phỏng số của ngành Cơ kỹ thuật Việt Nam trên trường quốc tế.