Tổng quan về luận án

Trong cơ học tính toán hiện đại, phương pháp Phần tử Hữu hạn (Finite Element Method - FEM) khởi xướng từ nghiên cứu nền tảng của Courant (1941) và được hệ thống hóa toán học bởi Strang và Fix (1973) đóng vai trò trụ cột trong việc giải các phương trình đạo hàm riêng mô tả cơ học vật rắn biến dạng. Dẫu vậy, FEM truyền thống vẫn tồn tại các rào cản nội tại nghiêm trọng: trường đạo hàm (ứng suất, biến dạng) bị gián đoạn tại các biên phần tử và nút liên kết; hiện tượng khóa thể tích (volumetric locking) và khóa cắt (shear locking) làm suy giảm độ chính xác; cũng như chi phí tính toán tăng vọt khi phải tinh chỉnh lưới cục bộ (remeshing) trong các bài toán lan truyền vết nứt. Nhằm khắc phục những giới hạn này, luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 9520101.01) của nghiên cứu sinh Nguyễn Đình Dư, 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 tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội (2023), đã tiên phong mở rộng và hoàn thiện Phương pháp Phần tử Hữu hạn Nội suy Kép (Consecutive-Interpolation Finite Element Method - CFEM).

Research gap cốt lõi mà luận án giải quyết là sự thiếu vắng một công thức hàm dạng tổng quát có thể áp dụng nhất quán cho các phần tử 2D và 3D bậc thấp mà không làm tăng bậc tự do (Degrees of Freedom - DOF), đồng thời thiếu một khung tính toán hiệu năng cao kết hợp giữa kỹ thuật làm giàu mở rộng (XFEM) với sơ đồ tích phân số thay thế để mô phỏng chính xác cơ học phá hủy tĩnh/động trong vật liệu biến đổi chức năng (Functionally Graded Materials - FGM) và phân tích phi tuyến hình học.

Nghiên cứu tập trung vào 4 câu hỏi và giả thuyết khoa học trọng tâm:

  1. RQ1 & H1: Liệu thủ tục nội suy kép (CIP) có thể xây dựng thành công trên công thức thống nhất cho phần tử 2D (tam giác CT3, tứ giác CQ4) và 3D (tứ diện CTH4, lục diện CHH8) nhằm triệt tiêu sự gián đoạn ứng suất nút mà vẫn giữ nguyên thuộc tính Kronecker-delta mà không phát sinh thêm DOF ẩn?
  2. RQ2 & H2: Sự kết hợp giữa CIP và kỹ thuật làm giàu mở rộng (phần tử XCQ4) kết hợp hàm dốc (ramp function) có vượt trội hơn XFEM truyền thống (XQ4) trong việc xác định hệ số cường độ ứng suất tĩnh/động ($K_I, K_{II}$, DSIFs) và quỹ đạo phát triển vết nứt trên vật liệu dị tính FGM hay không?
  3. RQ3 & H3: Trong phân tích phi tuyến hình học 2D/3D (biến dạng hữu hạn, chuyển vị lớn), việc tích hợp phần tử CIP với các mô hình tích phân thay thế (Element Mid-point - EM, Element Mid-face - EF) có đảm bảo độ chính xác ngang bằng tích phân Gauss đầy đủ nhưng giảm thiểu đáng kể thời gian CPU và kháng hiện tượng khóa phần tử hay không?
  4. RQ4 & H4: Ứng xử cơ học của FGM khi chịu tải cơ học phức tạp có thể giảm thiểu đỉnh ứng suất và chuyển dịch vị trí phá hủy nguy hiểm ra khỏi vùng chịu lực tập trung hay không?

Khung lý thuyết của công trình dựa trên Lý thuyết đàn hồi tuyến tính (Goodier, 1951), Cơ học phá hủy đàn hồi tuyến tính (LEFM - Griffith-Irwin), Cơ học môi trường liên tục phi tuyến (Belytschko et al., 2014), và Lý thuyết làm giàu phân vùng đơn vị (Partition of Unity - Melenk & Babuška, 1996).

Đóng góp đột phá của luận án được định lượng cụ thể: Phần tử CQ4 đạt cùng độ chính xác sai số năng lượng 0.3% chỉ mất 2.68 giây (lưới $24\times6$), trong khi phần tử Q4 truyền thống tiêu tốn 5.58 giây (lưới $96\times24$) – giảm hơn 51.9% thời gian xử lý; phần tử CHH8 đạt sai số năng lượng 0.53% với 960 phần tử (141 giây), so với 2240 phần tử HH8 (151 giây); và mô hình tích phân 3D-EM (9 điểm), EF (6 điểm) giúp tiết kiệm từ 60% đến 77.8% số điểm tích phân so với Gauss 27 điểm ($3\times3\times3$), triệt tiêu nguy cơ phân kỳ trong thuật toán phi tuyến Newton-Raphson.

Phạm vi nghiên cứu bao quát các bài toán đàn hồi tuyến tính tĩnh/động, phân tích dao động tự do, cơ học nứt FGM 2D và phi tuyến hình học 2D/3D trên các mẫu kết cấu tiêu chuẩn quốc tế (dầm công-xôn, dầm Cook, tấm kéo/uốn 3 điểm TPB/4 điểm, vỏ trụ dày, khung snap-through, kết cấu nhà cao tầng 540 phần tử). Toàn bộ thuật toán được tác giả tự lập trình độc lập trên môi trường mã nguồn MATLAB kết hợp đối chuẩn kiểm chứng nghiêm ngặt với phần mềm thương mại ABAQUS.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp số giải bài toán cơ học vật rắn chứng kiến ba dòng tiếp cận chính. Dòng thứ nhất là Phương pháp Phần tử Hữu hạn mở rộng (XFEM) khởi xướng bởi Belytschko và Black (1999), Moës et al. (1999), phát triển tiếp bởi Fries và Belytschko (2010), cho phép mô phỏng vết nứt độc lập với cấu trúc lưới bằng cách bổ sung hàm Heaviside và hàm tiệm cận đỉnh nứt $r^{1/2}\sin(\theta/2)$. Tuy nhiên, XFEM vẫn kế thừa hàm dạng Lagrange tuyến tính của FEM gốc nên trường ứng suất quanh đỉnh nứt vẫn bị thô và phụ thuộc nặng vào mật độ lưới. Dòng thứ hai là Phân tích Đẳng hình học (Isogeometric Analysis - IGA) do Hughes et al. (2005), Cottrell et al. (2009) đề xuất, sử dụng hàm B-splines/NURBS để đạt độ liên tục bậc cao $C^p$. Mặc dù vậy, IGA gặp trở ngại lớn trong việc gán điều kiện biên do thiếu thuộc tính Kronecker-delta ($\phi_I(x_J) \neq \delta_{IJ}$), đòi hỏi các kỹ thuật phức tạp như nhân tử Lagrange hay hàm phạt (Nguyen-Thanh et al., 2017). Dòng thứ ba là các Phương pháp Không lưới (Meshfree Methods) như EFG (Belytschko et al., 1994), RKPM (Liu et al., 1995), RPIM (Wang & Liu, 2002), Bui et al. (2018). Mặc dù không cần lưới, meshfree lại phụ thuộc chặt vào bán kính hỗ trợ (support domain) do người dùng tự chọn mà không có tiêu chuẩn toán học tối ưu, đồng thời chi phí tích phân số rất cao.

Về hướng cải tiến trực tiếp trên nền tảng FEM, Phương pháp Phần tử Hữu hạn Làm mịn (Smoothed FEM - S-FEM) do Liu et al. (2016) hệ thống hóa qua các biến thể cell-based (CS-FEM), node-based (NS-FEM), edge-based (ES-FEM) và face-based (FS-FEM) đã làm mịn trường biến dạng qua các miền đệm. Tuy nhiên, S-FEM đòi hỏi các điểm nút mở rộng và làm phức tạp hóa thuật toán dựng ma trận độ cứng. Khắc phục nhược điểm này, thủ tục CIP nguyên bản được Zheng et al. (2010) đề xuất cho phần tử tam giác TFEM, và sau đó được PGS.TS Bùi Quốc Tính cùng nhóm cộng sự phát triển cho phần tử tứ giác CQ4 (Bui et al., 2014) và XCQ4 (Bui et al., 2015).

                      [PHƯƠNG PHÁP SỐ TRONG CƠ HỌC VẬT RẮN]
                                       |
       +-------------------------------+-------------------------------+
       |                               |                               |
  [FEM & XFEM]                 [IGA & MESHFREE]                 [CFEM / CIP]
  - Gián đoạn ứng suất         - Mất Kronecker-delta           - Liên tục ứng suất
  - Bị khóa thể tích/cắt       - Khó gán điều kiện biên        - Giữ nguyên số DOF
  - Remeshing phức tạp         - Chi phí tích phân cao         - Thuộc tính Kronecker-delta
       |                               |                               |
       +-------------------------------+-------------------------------+
                                       |
                       [ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN]
              - Công thức toán học thống nhất 2D (CT3, CQ4) & 3D (CTH4, CHH8)
              - Phần tử XCQ4 làm giàu hàm dốc cho phá hủy động FGM
              - Sơ đồ tích phân 3D-EM (9 pts) & EF (6 pts) cho phi tuyến hình học

Luận án này định vị chính xác vào khoảng trống khoa học: mở rộng CIP thành một công thức toán học thống nhất cho toàn bộ hệ phần tử 2D/3D (CT3, CQ4, CTH4, CHH8), kết hợp cơ chế làm giàu đỉnh nứt FGM mới không dùng hàm tiệm cận phức tạp, và tích hợp sơ đồ tích phân rút gọn ổn định cho phân tích phi tuyến hình học.

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

  • So với nghiên cứu của Kim và Paulino (2002, 2004) sử dụng phần tử FGM đẳng tham số bậc cao Q8 với tích phân tương tác chuẩn: Phần tử CQ4 của luận án đạt độ chính xác tương đương Q8 (sai số ứng suất $\sigma_{xx}$ tại điểm B' trên thanh TiB/Ti chỉ 0.07% so với 1.25% của Q4) nhưng chỉ cần 4 nút, tiết kiệm 50% số DOF trên mỗi phần tử.
  • So với công trình của Song và Paulino (2006) và Chafi & Boulenouar (2019) về phá hủy động trong FGM: Luận án phát triển phần tử XCQ4 kết hợp hàm làm giàu dạng dốc (ramp function) thay thế 4 hàm tiệm cận $r^{1/2}$, giúp giảm 4 DOF phụ trên mỗi nút đỉnh nứt nhưng vẫn thu được giá trị DSIF bám sát giải tích và dữ liệu thực nghiệm của Parameswaran & Shukla (2001).
  • So với các nghiên cứu phi tuyến hình học của Bordas et al. (2017) và phần tử hiệu chỉnh biến dạng QT10MS (2018), NICE-H8 (2009): Mô hình tích phân thay thế 3D-EM và EF của luận án triệt tiêu chi phí dư thừa của tích phân Gauss 27 điểm mà không gây ra hiện tượng mất ổn định hay phân kỳ như Gauss 1 điểm.

Đó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à hoàn thiện các lý thuyết cơ học tính toán nền tảng thông qua các luận điểm khoa học có bằng chứng xác thực:

Thứ nhất, luận án đóng góp vào Lý thuyết nội suy hàm dạng phần tử hữu hạn bằng việc thiết lập công thức tổng quát 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{x}) \mathbf{a}$$ trong đó hàm dạng CIP $R_I(\mathbf{x})$ được hợp thành từ hàm dạng Lagrange tiêu chuẩn $N_I$ kết hợp với các đạo hàm trung bình nút $\bar{N}{I,x}, \bar{N}{I,y}, \bar{N}{I,z}$ thông qua các hàm bổ trợ phi tuyến $\Phi_I, \Phi{Ix}, \Phi_{Iy}, \Phi_{Iz}$: $$R_I(\mathbf{x}) = \Phi_I(\mathbf{x}) N_I + \Phi_{Ix}(\mathbf{x}) \bar{N}{I,x} + \Phi{Iy}(\mathbf{x}) \bar{N}{I,y} + \Phi{Iz}(\mathbf{x}) \bar{N}{I,z}$$ Đạo hàm trung bình nút được tính toán dựa trên tập hợp các phần tử lân cận $S_I$ có chung nút $I$ với trọng số thể tích/diện tích $w_e = A_e / \sum{j \in S_I} A_j$: $$\bar{N}{I,w} = \sum{e \in S_I} w_e N_{I,w}^{|e|}$$ Luận án đã chứng minh chặt chẽ về mặt giải tích toán học rằng hàm bổ trợ $\Phi$ thỏa mãn hoàn toàn điều kiện Kronecker-delta: $\Phi_I(\mathbf{x}J) = \delta{IJ}$ và triệt tiêu đạo hàm tại nút $\Phi_{I,w}(\mathbf{x}_J) = 0$, đảm bảo tính tương thích và liên tục của trường gradient giữa các biên phần tử.

Thứ hai, luận án làm phong phú Cơ học phá hủy đàn hồi tuyến tính FGM bằng cách tích hợp gradient vật liệu không gian $E(\mathbf{x}) = E_0 e^{\beta x}$ và $\rho(\mathbf{x}) = \rho_0 e^{\beta x}$ vào phương trình làm giàu trường chuyển vị: $$\mathbf{u}(\mathbf{x}) = \sum_{I \in \mathcal{N}} R_I(\mathbf{x}) \mathbf{u}I + \sum{J \in \mathcal{N}{cut}} R_J(\mathbf{x}) H(\mathbf{x}) \mathbf{a}J + \sum{K \in \mathcal{N}{tip}} R_K(\mathbf{x}) \left[ \sum_{\alpha=1}^{4} F_\alpha(r, \theta) \mathbf{b}_K^\alpha \right]$$ đồng thời chứng minh tính khả thi của kỹ thuật làm giàu bằng hàm dốc tuyến tính độc lập, giảm thiểu hiện tượng ma trận độ cứng bị suy biến (ill-conditioning).

Thứ ba, luận án tạo bước tiến trong Lý thuyết cơ học phi tuyến hình học khi chứng minh rằng việc tăng bậc liên tục của trường đạo hàm thông qua CIP giúp phần tử bậc thấp (CQ4, CHH8) đạt được độ linh hoạt chống lại hiện tượng khóa biến dạng mà trước đây bắt buộc phải dùng phần tử bậc hai (Q8, HH20).

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

Khung phân tích của luận án tích hợp liên hoàn ba trụ cột: (1) Kỹ thuật nội suy kép CIP giải quyết tính gián đoạn trường ứng suất; (2) Phương pháp Tập mức (Level Set Method - LSM) và Tích phân tương tác $J$ giải quyết bài toán kỳ dị đỉnh nứt; (3) Sơ đồ tích phân trung điểm thay thế giải quyết độ phức tạp tính toán phi tuyến.

                      +-----------------------------------------+
                      |   KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN   |
                      +-----------------------------------------+
                                           |
         +---------------------------------+---------------------------------+
         |                                 |                                 |
         v                                 v                                 v
+------------------+             +-------------------+             +------------------+
|  THỦ TỤC NỘI SUY |             | LÀM GIÀU MỞ RỘNG  |             | TÍCH PHÂN SỐ MỚI |
|    KÉP (CIP)     |             |    (X-CFEM/LSM)   |             | (3D-EM & EF)     |
| CT3, CQ4, CTH4,  |             | XCQ4 kết hợp      |             | 9 điểm (EM) &    |
| CHH8             |             | Ramp Function &   |             | 6 điểm (EF) thay |
| Liên tục gradient|             | Dynamic J-integral|             | thế 27 điểm Gauss|
+------------------+             +-------------------+             +------------------+
         |                                 |                                 |
         +---------------------------------+---------------------------------+
                                           |
                                           v
                       +---------------------------------------+
                       |    HIỆU NĂNG TÍNH TOÁN VƯỢT TRỘI      |
                       | - Giữ nguyên số DOF như FEM gốc       |
                       | - Độ chính xác tương đương bậc cao   |
                       | - Tiết kiệm 50-77% chi phí CPU        |
                       +---------------------------------------+

Boundary conditions được xác lập rõ ràng: Vật liệu làm việc trong miền đàn hồi tuyến tính đẳng hướng hoặc FGM liên tục theo quy luật hàm mũ/tuyến tính; biến dạng hình học lớn tuân theo mô tả Total Lagrangian hoặc Co-rotational; trạng thái ứng suất phẳng/biến dạng phẳng cho 2D và không gian 3 chiều thuần túy cho 3D.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Triết học Thực chứng (Positivism) kết hợp Chủ nghĩa Hiện thực Cực đoan (Critical Realism) trong khoa học kỹ thuật tính toán: kiểm chứng các quy luật vật lý khách quan thông qua mô hình toán học giải tích, mô phỏng số học định lượng và đối chuẩn với dữ liệu thực nghiệm tiêu chuẩn.

Thiết kế nghiên cứu đa cấp độ (Multi-level computational design):

  • Cấp độ phần tử (Local element level): Kiểm tra tính hội tụ năng lượng biến dạng, thuộc tính Kronecker-delta, đạo hàm triệt tiêu tại nút và khả năng biến dạng không khóa.
  • Cấp độ kết cấu 2D (Structural 2D level): Phân tích tĩnh dầm công-xôn chịu cắt parabol, thanh FGM chịu kéo phức tạp, tấm FGM nứt tâm (CCTT) chịu tải xung bậc thang và tải điều hòa sin, dầm nứt uốn 3 điểm (TPB) và uốn 4 điểm.
  • Cấp độ kết cấu 3D (Structural 3D level): Dầm công-xôn mặt cắt chữ T, hình trụ rỗng có lỗ khoan lệch tâm (dao động tự do), dầm Cook 3D gần như không nén được ($\nu = 0.4999$), và kết cấu khung nhà cao tầng 3D chịu tải xô ngang $q = 1600 \text{ kN/m}^2$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số hóa bao gồm các giao thức kiểm chuẩn đa tầng (Cross-validation protocols):

+-----------------------------------------------------------------------------------+
|                        QUY TRÌNH NGHIÊN CỨU & KIỂM CHUẨN                          |
+-----------------------------------------------------------------------------------+
| 1. XÂY DỰNG CÔNG THỨC TOÁN HỌC (CIP 2D/3D, Ma trận Jacobian, Hàm hỗ trợ Phi)      |
|                                         |                                         |
| 2. THIẾT LẬP MÔ HÌNH SỐ (MATLAB Engine: Dựng ma trận độ cứng K, Vector lực F)     |
|                                         |                                         |
| 3. XỬ LÝ PHI TUYẾN & ĐỘNG LỰC (Thuật toán Newton-Raphson, Arc-Length, Newmark-beta)|
|                                         |                                         |
| 4. KIỂM ĐỊNH TAM GIÁC (Triangulation):                                            |
|    - Lời giải giải tích chính xác (Goodier, Timoshenko)                           |
|    - Mô phỏng phần tử bậc cao thương mại (ABAQUS HH20, 10800 phần tử siêu mịn)    |
|    - Dữ liệu thực nghiệm nứt quốc tế (Parameswaran & Shukla, 2001; Tilbrook, 2009)|
+-----------------------------------------------------------------------------------+

Tính hợp lệ (Validity) và độ tin cậy (Reliability):

  • Construct Validity: Được xác lập thông qua việc chứng minh ma trận Jacobian khả nghịch ($\det(\mathbf{J}) > 0$) trên toàn miền tự nhiên $(\xi, \eta, \zeta) \in [-1, 1]^3$ và sự triệt tiêu hoàn toàn của các đạo hàm hàm bổ trợ tại nút.
  • Internal Validity: Kiểm soát tuyệt đối sai số làm tròn (round-off error) và sai số rời rạc hóa (discretization error) bằng việc thực hiện phân tích hội tụ trên 4 đến 5 cấp độ lưới phân giải tăng dần ($h, h/2, h/4, h/8$).
  • External Validity & Triangulation: Tam giác hóa dữ liệu qua 3 nguồn độc lập: (1) Lời giải giải tích lý thuyết; (2) Kết quả số từ các phần tử bậc cao trong ABAQUS; (3) Dữ liệu thực nghiệm cơ học nứt quốc tế.

Data và phân tích

Dữ liệu nghiên cứu là tập hợp các đại lượng cơ học chuẩn hóa thu được từ hàng trăm kịch bản tính toán số:

  • Phần mềm phát triển: Tác giả tự xây dựng toàn bộ mã nguồn xử lý phần tử hữu hạn trên MATLAB (R2018b/R2020a), giao tiếp và đối chiếu dữ liệu chuẩn với ABAQUS Standard/Explicit.
  • Thuật toán giải phi tuyến: Newton-Raphson biến đổi kết hợp kỹ thuật kiểm soát độ dài cung (Arc-length method) để bắt chính xác đường cong tải trọng - chuyển vị vượt qua điểm tới hạn lật kết cấu (snap-through phenomenon).
  • Động lực học phi tuyến: Sử dụng sơ đồ tích phân thời gian ẩn Newmark-$\beta$ ($\gamma = 0.5, \beta = 0.25$) đảm bảo tính ổn định vô điều kiện cho bài toán transient dynamic crack.
  • Kiểm tra tính vững (Robustness checks): Đánh giá hành vi phần tử trên cả lưới có quy tắc (structured mesh) và lưới bất quy tắc/méo mó nghiêm trọng (distorted/unstructured mesh).

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

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

Luận án đã đạt được 5 phát hiện cốt lõi mang tính đột phá khoa học kỹ thuật:

+-----------------------------------------------------------------------------------------------+
|                            TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CHÍNH                               |
+-----------------------------------------------------------------------------------------------+
| 1. KHẮC PHỤC TRIỆT ĐỂ SỰ GIÁN ĐOẠN ỨNG SUẤT NÚT:                                              |
|    - CQ4, CHH8 cho trường ứng suất trơn mịn liên tục tuyệt đối, không cần thuật toán lọc mượt.|
|                                                                                               |
| 2. HIỆU QUẢ TÍNH TOÁN & HỘI TỤ VƯỢT BẬC:                                                      |
|    - CQ4 đạt sai số 0.3% chỉ mất 2.68s (lưới 24x6) so với Q4 mất 5.58s (lưới 96x24).          |
|    - CHH8 đạt sai số 0.53% tốn 141s (960 phần tử), HH8 tốn 151s (2240 phần tử).              |
|                                                                                               |
| 3. TỐI ƯU HÓA ỨNG SUẤT TRONG VẬT LIỆU FGM:                                                   |
|    - FGM (TiB/Ti) giúp giảm đỉnh ứng suất và dịch chuyển vị trí phá hủy từ A' sang B'.       |
|    - Sai số sigma_xx của CQ4 so với Q8 chỉ là 0.07% (trong khi Q4 sai số 1.25%).              |
|                                                                                               |
| 4. CHÍNH XÁC CAO TRONG TÍNH TOÁN CƠ HỌC NỨT ĐỘNG (DSIFs):                                     |
|    - Phần tử XCQ4 làm giàu hàm dốc nắm bắt hoàn hảo sóng ứng suất sigma_yy và dao động DSIF.   |
|                                                                                               |
| 5. ĐỘT PHÁ SƠ ĐỒ TÍCH PHÂN PHI TUYẾN 3D:                                                      |
|    - 3D-EM (9 điểm) và EF (6 điểm) tiết kiệm 66.7 - 77.8% số điểm tích phân so với Gauss (27).|
|    - Đảm bảo hội tụ tuyệt đối và triệt tiêu nguy cơ phân kỳ của Gauss rút gọn 1 điểm.         |
+-----------------------------------------------------------------------------------------------+

Finding 1: Xóa bỏ hoàn toàn sự gián đoạn trường ứng suất mà không cần hậu xử lý. Trích dẫn trực tiếp từ luận án: "Các trường ứng suất được tính toán bởi các phần tử tăng cường CIP là mượt mà, trong khi các trường bằng phương pháp FEM truyền thống là không liên tục về mặt vật lý tại nút." Kết quả phân tích dầm công-xôn chữ T và tấm FGM cho thấy đường đẳng ứng suất $\sigma_{xx}, \sigma_{zz}$ của phần tử CQ4 và CHH8 hoàn toàn liên tục qua các biên phần tử, loại bỏ hoàn toàn các bước nhảy ứng suất giả tạo của Q4 và HH8.

Finding 2: Tốc độ hội tụ và hiệu năng thời gian vượt trội của phần tử CIP. Trong bài toán dầm công-xôn 2D chịu cắt parabol, để đạt cùng mức sai số năng lượng 0.3%, phần tử CQ4 chỉ cần lưới $24\times6$ với thời gian chạy mã lệnh là 2.68 giây, trong khi phần tử Q4 truyền thống bắt buộc phải chia lưới mịn tới $96\times24$ và tiêu tốn 5.58 giây (tiết kiệm 51.9% thời gian tính toán). Trong bài toán dầm công-xôn 3D chữ T, phần tử CHH8 chỉ cần 960 phần tử để đạt sai số năng lượng 0.53% trong 141 giây, trong khi phần tử HH8 tiêu chuẩn phải dùng đến 2240 phần tử và mất 151 giây.

Finding 3: Khám phá cơ chế phân phối lại ứng suất và dịch chuyển vị trí phá hoại trên kết cấu FGM. Đối với thanh kết cấu FGM phức tạp cấu tạo từ hợp kim Titan tinh khiết và Titan Monoboride (Ti/TiB, $E_0=375 \text{ GPa}, \nu_0=0.14$ đến $E_1=107 \text{ GPa}, \nu_1=0.34$), luận án phát hiện: "Ứng suất lớn nhất trong trường hợp vật liệu FGM có giá trị nhỏ hơn ứng suất lớn nhất trong vật liệu đồng chất. Điều này là do vật liệu FGM có sự phân phối lại ứng suất. Mặt khác, vị trí ứng suất cực đại cũng có sự khác biệt, vật liệu FGM thì tại điểm B' trong khi vật liệu đồng nhất là tại điểm A'." Phát hiện này có ý nghĩa thực tiễn to lớn trong việc chủ động điều hướng điểm phá hủy ra khỏi các vị trí tập trung ứng suất nguy hiểm. Sai số ứng suất chính $\sigma_{xx}$ của phần tử CQ4 so với phần tử bậc cao Q8 chỉ là 0.07%, vượt trội hoàn toàn so với mức sai số 1.25% của Q4.

Finding 4: Độ chính xác vượt bậc trong phân tích dao động tự do và động lực học nứt FGM. Trong bài toán dao động riêng của hình trụ rỗng có lỗ khoan lệch tâm (độ lệch tâm $e_c = 80%$), sai số tần số dao động mode 4 không thứ nguyên ($\omega \sqrt{\rho/G}$) của phần tử CHH8 (lưới 3072 phần tử) chỉ là 0.52% so với lời giải giải tích của Mirzaei (2012), trong khi HH8 sai số lên tới 1.83%. Trong bài toán phá hủy động tấm FGM chịu tải xung, phần tử XCQ4 kết hợp hàm dốc đã mô phỏng chính xác sự biến thiên của hệ số cường độ ứng suất động DSIF ($K_I(t), K_{II}(t)$) theo thời gian, bám sát các kết quả đối chuẩn quốc tế của Song & Paulino (2006).

Finding 5: Đột phá hiệu năng của sơ đồ tích phân 3D-EM và EF trong phân tích phi tuyến. Trong bài toán kết cấu nhà cao tầng 3D chịu tải trọng xô ngang $q = 1600 \text{ kN/m}^2$ (10 bước tải) và dầm Cook biến dạng lớn, sơ đồ tích phân 3D-EM (9 điểm) và EF (6 điểm) cho đường cong tải trọng - chuyển vị tại đỉnh trùng khít với tích phân Gauss đầy đủ 27 điểm ($3\times3\times3$), nhưng giảm hơn 66.7% đến 77.8% số điểm tích phân, giảm thời gian tính toán ma trận độ cứng tiếp tuyến từ 60.5 giây xuống còn 16.6 giây.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp nền tảng toán học mở rộng hoàn chỉnh cho họ phần tử CIP, chứng minh tính khả thi của việc nâng cao bậc hàm dạng mà không làm tăng ma trận DOF tổng thể.
  • Về phương pháp luận: Thiết lập một quy trình chuẩn hóa kết hợp giữa làm giàu mở rộng (XFEM) với nội suy kép (CIP), có thể chuyển giao trực tiếp để giải quyết các bài toán ghép đa trường vật lý (Multiphysics) như nhiệt - cơ, điện - từ - đàn hồi.
  • Về ứng dụng thực tiễn: Cung cấp công cụ mô phỏng độ tin cậy cao cho các kỹ sư kết cấu hàng không vũ trụ, hạt nhân và công trình biển trong việc đánh giá tuổi thọ mỏi, cơ chế nứt gãy và độ ổn định phi tuyến của vật liệu composite FGM.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn khoa học cần tiếp tục hoàn thiện:

  1. Hạn chế về khóa thể tích trong vật liệu siêu đàn hồi phi tuyến: Đối với vật liệu gần như không nén được (hệ số Poisson $\nu \to 0.4999$, mô hình neo-Hookean), phần tử CHH8 dù vượt trội hơn phần tử tuyến tính HH8 nhưng vẫn xuất hiện sai số biến dạng nhất định khi so sánh với phần tử bậc hai HH20 hoặc các phần tử chuyên dụng tích hợp biến dạng giả định như QT10MS (2018) hay NICE-H8 (2009).
  2. Hiện tượng phân kỳ của tích phân Gauss rút gọn 1 điểm: Việc sử dụng tích phân Gauss 1 điểm trong phân tích phi tuyến hình học 3D dẫn đến hiện tượng ma trận độ cứng bị suy biến dạng zero-energy modes, gây phân kỳ thuật toán Newton-Raphson.
  3. Phạm vi cấu tạo vật liệu: Luận án giới hạn ở mô hình vật liệu đàn hồi tuyến tính và phi tuyến hình học, chưa tích hợp mô hình phi tuyến vật liệu (plasticity/viscoplasticity) hoặc cơ chế phá hủy dẻo kết hợp hư hại liên tục (damage mechanics).

Chương trình nghiên cứu 5-10 năm tới mở ra các hướng phát triển chiến lược:

  • Lồng ghép kỹ thuật biến dạng giả định (Assumed Strain / Enhanced Assumed Strain - EAS) vào công thức CIP của phần tử CQ4 và CHH8 để triệt tiêu hoàn toàn hiện tượng khóa thể tích khi $\nu \to 0.5$.
  • Mở rộng phương pháp X-CFEM sang bài toán lan truyền vết nứt 3D không gian phức tạp và bài toán mỏi nhiệt - cơ động (Thermo-mechanical dynamic fracture) trong môi trường nhiệt độ cao.
  • Phát triển module phần mềm thương mại hóa CIP như một plugin/add-on tích hợp trực tiếp vào các hệ thống CAE tiêu chuẩn như ABAQUS UEL (User Element) hoặc OpenFOAM.

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

Nghiên cứu tạo ra tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật (Academic Impact): Các kết quả nghiên cứu trong luận án đã được công bố trên các tạp chí quốc tế ISI/Scopus hàng đầu (Q1) thuộc danh mục uy tín của ngành Cơ học tính toán như Applied Mathematical Modelling (Elsevier, 2021), Theoretical and Applied Fracture Mechanics, và Finite Elements in Analysis and Design. Dự báo công trình sẽ thu hút hàng trăm trích dẫn khoa học trong các nghiên cứu về phần tử gradient nâng cao và cơ học phá hủy FGM.
  • Chuyển đổi công nghiệp (Industry Transformation): Giúp ngành công nghiệp chế tạo hàng không vũ trụ và năng lượng tối ưu hóa quy trình thiết kế cánh tuabin, vỏ tàu vũ trụ bằng vật liệu FGM, giảm thiểu hệ số an toàn dư thừa và tiết kiệm từ 15% đến 20% chi phí thử nghiệm phá hủy thực tế nhờ mô phỏng số chính xác cao.
  • Chính sách và đào tạo: Cung cấp bằng chứng khoa học và công cụ mô phỏng phụ vụ việc xây dựng tiêu chuẩn tính toán kết cấu FGM tại Việt Nam, đồng thời đóng góp giáo trình chuyên đề chuyên sâu cho bậc đào tạo thạc sĩ, tiến sĩ ngành Cơ học kỹ thuật và Kỹ thuật công trình.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học tính toán: Tiếp cận mã nguồn thuật toán, công thức toán học tường minh về hàm dạng CIP 2D/3D và phương pháp làm giàu đỉnh nứt để phát triển các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư R&D trong công nghiệp Hàng không, Ô tô, Xây dựng: Ứng dụng công cụ tính toán hiệu năng cao để phân tích độ bền mỏi, dự đoán chính xác đường lan truyền vết nứt và phân tích ổn định phi tuyến của kết cấu phức tạp mà không bị phụ thuộc vào bản quyền phần mềm đắt đỏ.
  • Các nhà hoạch định chính sách và Viện nghiên cứu: Sử dụng dữ liệu khoa học tin cậy để thẩm định an toàn các công trình trọng điểm quốc gia, kết cấu cao tầng chịu tải trọng gió bão, động đất.

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à đã mở rộng lý thuyết nào? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc thiết lập Công thức toán học thống nhất cho Thủ tục Nội suy Kép (CIP) áp dụng đồng thời cho hệ phần tử 2D (CT3, CQ4) và 3D (CTH4, CHH8). Công trình đã mở rộng Lý thuyết nội suy Lagrange trong FEM cổ điển (Strang & Fix, 1973) bằng cách nhúng các đạo hàm trung bình nút $\bar{N}_{I,w}$ vào hàm dạng thông qua các hàm bổ trợ $\Phi$ phi tuyến, tạo ra trường ứng suất liên tục $C^0$-gradient mà hoàn toàn không phát sinh thêm bậc tự do ẩn (DOFs) và bảo toàn trọn vẹn thuộc tính Kronecker-delta.

2. Đổi mới 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ế tiền nhiệm? Trả lời: So với XFEM truyền thống của Moës et al. (1999), luận án tạo ra phần tử XCQ4 kết hợp hàm dốc tuyến tính (ramp function) giúp giảm 4 DOF tại mỗi nút đỉnh nứt nhưng vẫn nắm bắt chính xác điểm kỳ dị $r^{-1/2}$; so với Phân tích phi tuyến của Bordas et al. (2017), luận án phát triển sơ đồ tích phân thay thế 3D-EM (9 điểm) và EF (6 điểm) thay thế cho tích phân Gauss 27 điểm, giúp giảm hơn 66% thời gian lập ma trận độ cứng tiếp tuyến mà không bị suy biến năng lượng.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì và giải thích cơ chế lý thuyết? Trả lời: Phát hiện bất ngờ nhất là hiện tượng phân phối lại ứng suất và dịch chuyển vị trí tập trung ứng suất cực đại trong kết cấu FGM (TiB/Ti). Trong khi vật liệu đồng nhất đạt đỉnh ứng suất tại điểm hình học góc lõm $A'$, vật liệu FGM dịch chuyển vị trí cực đại sang mép biên ngoài $B'$ đồng thời giảm biên độ ứng suất đỉnh. Cơ chế lý thuyết là do gradient mô đun đàn hồi $E(y)$ và hệ số Poisson $\nu(y)$ biến thiên không gian liên tục đã làm mềm cục bộ vùng chịu lực, triệt tiêu sự tích tụ năng lượng biến dạng tập trung.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) hay không? Trả lời: Luận án cung cấp đầy đủ và tường minh toàn bộ giao thức toán học: bảng tọa độ tự nhiên, trọng số của các điểm tích phân trong mô hình EM và EF (Bảng 3.4, 3.5, 3.6), công thức giải tích của các hàm bổ trợ $\Phi_I, \Phi_{Ix}, \Phi_{Iy}, \Phi_{Iz}$, các tham số vật liệu cơ sở ($E, \nu, \rho$), cùng quy trình thiết lập ma trận Jacobian chuyển đổi hệ tọa độ, cho phép các nhà nghiên cứu độc lập tái lập 100% thuật toán trên bất kỳ ngôn ngữ lập trình nào (MATLAB, Fortran, C++ hoặc Python).

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Lộ trình 10 năm tập trung vào 3 trục chính: (1) Tích hợp lý thuyết biến dạng giả định (Assumed Strain) vào phần tử CHH8 để kháng khóa thể tích tuyệt đối cho vật liệu siêu đàn hồi và mô hình cao su sinh học; (2) Mở rộng X-CFEM sang bài toán lan truyền vết nứt 3D động lực học kết hợp hiệu ứng nhiệt - áp điện (Thermo-electro-mechanical fracture); (3) Tối ưu hóa thuật toán tính toán song song trên GPU và đóng gói thành thư viện mã nguồn mở cho cộng đồng cơ học tính toán toàn cầu.

Kết luận

Luận án tiến sĩ của NCS Nguyễn Đình Dư đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp cốt lõi:

  1. Thiết lập thành công công thức toán học thống nhất cho kỹ thuật nội suy kép (CIP) trên toàn bộ các phần tử hữu hạn chuẩn 2D (CT3, CQ4) và 3D (CTH4, CHH8), giải quyết triệt để vấn đề gián đoạn trường ứng suất nút mà không làm tăng số bậc tự do của bài toán.
  2. Phát triển phần tử tứ giác mở rộng nội suy kép XCQ4, ứng dụng thành công cho cơ học phá hủy đàn hồi tuyến tính FGM, chứng minh tính ưu việt của kỹ thuật làm giàu bằng hàm dốc thay thế cho các hàm tiệm cận truyền thống phức tạp.
  3. Đề xuất và hiện thực hóa thành công các sơ đồ tích phân số mới 3D-EM (9 điểm) và EF (6 điểm), giảm từ 60% đến 77.8% số điểm tích phân so với sơ đồ Gauss 27 điểm, mang lại bước đột phá về hiệu năng tính toán trong phân tích phi tuyến hình học.
  4. Chứng minh tính vượt trội của phần tử CIP bậc thấp trong việc đạt độ chính xác tương đương phần tử bậc hai (Q8, HH20) khi phân tích dầm công-xôn, dao động tự do hình trụ lệch tâm và kết cấu nhà cao tầng 3D chịu tải xô ngang.
  5. Phát hiện quy luật phân phối lại ứng suất trong kết cấu FGM, cung cấp cơ sở khoa học vững chắc cho việc thiết kế giảm thiểu rủi ro tập trung ứng suất phá hủy trong vật liệu composite tiên tiến.
  6. Mở ra 3 hướng nghiên cứu mới về cơ học phi tuyến vật liệu kết hợp biến dạng giả định, cơ học nhiệt nứt FGM 3D, và phát triển công cụ phần mềm tính toán kỹ thuật hiệu năng cao mang thương hiệu Việt Nam trên trường quốc tế.