Tổng quan về luận án

Trong bức tranh tổng thể của 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) đóng vai trò nền tảng kể từ công trình khởi xướng của A. Courant (1943) và sự phát triển bùng nổ của J. H. Argyris, R. W. Clough, O. C. Zienkiewicz, P. G. Ciarlet, G. Strang và G. Fix (1973). Mặc dù thống trị các bộ giải thương mại như ABAQUS, ANSYS, COMSOL hay NASTRAN, FEM truyền thống vẫn tồn tại các khiếm khuyết nội tại nghiêm trọng. Như luận án đã trích dẫn: "Các trường đạo hàm thu được từ FEM, như trường ứng suất và biến dạng trong cơ vật rắn, thì không liên tục về mặt vật lý tại các nút của phần tử", đồng thời phương pháp này rất nhạy cảm với hiện tượng méo lưới, khóa thể tích (volumetric locking) khi hệ số Poisson tiến sát giới hạn $\nu \to 0.5$, và khóa cắt (shear locking) trong phân tích kết cấu mỏng. 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, đã thực hiện nghiên cứu tiên phong: "Nghiên cứu ứng dụng phương pháp CFEM trong cơ học vật rắn biến dạng".

Nghiên cứu giải quyết trực tiếp khoảng trống học thuật (research gap) về việc nâng cao bậc trơn và độ chính xác của trường biến dạng – ứng suất mà không làm gia tăng số bậc tự do (Degrees of Freedom - DOF), khắc phục triệt để nhược điểm của các phương pháp làm mịn (S-FEM) hay phương pháp đẳng hình học (IGA) vốn đòi hỏi kỹ thuật áp đặt điều kiện biên phức tạp. Đề tài tập trung vào 4 câu hỏi nghiên cứu và giả thuyết cốt lõi:

  1. RQ1 & H1: Liệu thủ tục nội suy liên tiếp (Consecutive Interpolation Procedure - CIP) có thể xây dựng lại hàm dạng Lagrange để đạt được tính liên tục $C^0-C^1$ và giữ nguyên tính chất Kronecker-delta trên các phần tử 2D (tam giác CT3, tứ giác CQ4) và 3D (tứ diện CTH4, lục diện CHH8)?
  2. RQ2 & H2: Sự kết hợp giữa CIP và kỹ thuật làm giàu trường chuyển vị (XCFEM/XCQ4) với hàm dốc (ramp function) có thể triệt tiêu hoàn toàn sự gián đoạn ứng suất và nâng cao độ chính xác của hệ số cường độ ứng suất tĩnh (SIF) và động (DSIF) trong vật liệu biến tính chức năng (Functionally Graded Materials - FGM)?
  3. RQ3 & H3: Phương pháp CFEM có giải quyết được bài toán phi tuyến hình học lớn và hiện tượng mất ổn định snap-through mà không bị khóa thể tích khi vật liệu gần như không nén được?
  4. RQ4 & H4: Các mô hình tích phân thay thế (Element Mid-point - EM, Element Mid-face - EF) có thể giảm thiểu số điểm tích phân và chi phí tính toán so với tích phân Gauss truyền thống trong không gian 3D?

Khung lý thuyết của luận án được xây dựng dựa trên cơ học môi trường liên tục phi tuyến, lý thuyết cơ học phá hủy đàn hồi tuyến tính (LEFM) cho vật liệu dị tính FGM và giải thuật tích phân số tối ưu. Phạm vi khảo sát trải rộng từ các bài toán chuẩn 2D/3D (dầm công son, dầm chữ T, khung chịu tải phân bố) đến các kết cấu FGM phức tạp (hợp kim TiB/Ti với mô đun đàn hồi biến thiên $E_0 = 375\text{ GPa}, \nu_0 = 0.14$ đến $E_1 = 107\text{ GPa}, \nu_1 = 0.34$), phân tích dao động riêng hình trụ rỗng lệch tâm ($e_c = 0%, 40%, 80%$) và mô phỏng lan truyền vết nứt uốn 3 điểm, 4 điểm.

Literature Review và Positioning

Lịch sử phát triển các phương pháp số giải phương trình vi phân đạo hàm riêng trong cơ học vật rắn chứng kiến hai trường phái chính nhằm khắc phục hạn chế của FEM cổ điển:

  • Trường phái phát triển phương pháp số hoàn toàn mới: Tiêu biểu là phương pháp Đẳng hình học (IGA) khởi xướng bởi Hughes và cộng sự (2005), sử dụng hàm B-splines hoặc NURBS. Mặc dù IGA kiểm soát hoàn hảo bậc trơn giữa các phần tử, phương pháp này gặp khó khăn nghiêm trọng trong việc áp đặt điều kiện biên Dirichlet do thiếu tính chất Kronecker-delta, buộc phải sử dụng phương pháp nhân tử Lagrange hoặc hàm phạt rất tốn kém. Bên cạnh đó, họ phương pháp không lưới (Meshfree methods) như Element-Free Galerkin (EFG của Belytschko, 1994), Reproducing Kernel Particle Method (RKPM của Liu, 1995) và Radial Point Interpolation Method (RPIM của Wang, 2002) loại bỏ hoàn toàn khái niệm lưới nhưng lại phụ thuộc nặng nề vào các tham số mở rộng (dilation parameters) do người dùng tự chọn mà không có tiêu chuẩn xác định nghiệm tối ưu. Phương pháp Phần tử Biên (BEM) của Brebbia (1978) lại giới hạn ở các bài toán có nghiệm cơ bản (fundamental solutions).
  • Trường phái cải tiến trực tiếp từ nền tảng FEM: Điển hình là Phương pháp Phần tử Hữu hạn Làm mịn (Smoothed Finite Element Method - S-FEM) do G. R. Liu và cộng sự phát triển từ năm 2007 (bao gồm cell-based CS-FEM, node-based NS-FEM, edge-based ES-FEM, và face-based FS-FEM). Tuy nhiên, S-FEM can thiệp vào trường biến dạng thông qua các miền làm mịn mở rộng chứ không thay đổi hàm dạng ban đầu, dẫn đến việc xử lý cục bộ phức tạp. Phương pháp Phần tử Hữu hạn Mở rộng (XFEM) của Belytschko & Black (1999) và Moës et al. (1999) sử dụng phân hoạch đơn vị (Partition of Unity) để mô hình hóa vết nứt nhưng vẫn chịu ảnh hưởng bởi sự gián đoạn ứng suất tại các phần tử chuẩn.

Luận án định vị phương pháp Phần tử Hữu hạn Nội suy Kép (CFEM) là bước đột phá kết hợp ưu điểm của cả hai trường phái. Khởi nguồn từ công trình TFEM của C. Zheng (2010) và mở rộng cho bài toán nứt 2D bởi X. Peng (2017), kỹ thuật CIP được nhóm nghiên cứu của PGS. TS Bùi Quốc Tính và NCS Nguyễn Đình Dư chính thức chuẩn hóa và mở rộng cho phần tử tứ giác 4 nút (CQ4) năm 2014, phân tích dao động điện - cơ năm 2016, và vết nứt động năm 2017.

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

  1. So với nghiên cứu của Kim Jeong-Ho và Paulino G. H. (2002) về kết cấu FGM sử dụng phần tử tứ giác bậc cao 8 nút (Q8), phương pháp CQ4 của luận án đạt độ chính xác tương đương về trường ứng suất $\sigma_{xx}$ nhưng chỉ sử dụng phần tử 4 nút với chi phí xây dựng lưới thấp hơn đáng kể.
  2. So với mô hình phần tử mở rộng XFEM truyền thống của Moës et al. (1999), mô hình XCQ4 tích hợp hàm dốc trong luận án loại bỏ hiện tượng dao động số giả tạo tại lân cận đỉnh nứt và cho đường lan truyền vết nứt khớp chính xác với thực nghiệm uốn của dầm FGM.

Đó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 xấp xỉ số trong cơ học môi trường liên tục thông qua việc tổng quát hóa thủ tục nội suy liên tiếp CIP cho không gian 2D và 3D.

Cốt lõi toán học của CFEM là hàm dạng mở rộng được biểu diễn dưới dạng: $$\mathbf{R}I(\mathbf{x}) = \phi_I(\mathbf{x})\mathbf{N}I + \phi{Ix}(\mathbf{x})\overline{\mathbf{N}}{I,x} + \phi_{Iy}(\mathbf{x})\overline{\mathbf{N}}{I,y} + \phi{Iz}(\mathbf{x})\overline{\mathbf{N}}{I,z}$$ trong đó $\mathbf{N}I$ là hàm dạng Lagrange truyền thống; $\overline{\mathbf{N}}{I,x}, \overline{\mathbf{N}}{I,y}, \overline{\mathbf{N}}{I,z}$ là các đạo hàm trung bình theo trọng số diện tích/thể tích của các phần tử lân cận chia sẻ nút $I$: $$\overline{\mathbf{N}}{I,x} = \sum_{e \in S_I} w_e \mathbf{N}{I,x}^{(e)}, \quad w_e = \frac{\Delta_e}{\sum{\bar{e} \in S_I} \Delta_{\bar{e}}}$$ với $S_I$ là tập hợp các phần tử liên kết tại nút $I$ và $\Delta_e$ là diện tích (2D) hoặc thể tích (3D) của phần tử thứ $e$.

Luận án đã chứng minh một cách tường minh rằng các hàm phụ trợ $\phi_I, \phi_{Ix}, \phi_{Iy}, \phi_{Iz}$ thỏa mãn điều kiện triệt tiêu đạo hàm tại nút và bảo toàn tuyệt đối thuộc tính Kronecker-delta: $$\phi_I(\mathbf{x}p) = \delta{Ip}, \quad \phi_{I,w}(\mathbf{x}_p) = 0 \quad (\forall w \in {x, y, z})$$ Đóng góp này mở ra một bước chuyển dịch hệ hình (paradigm shift): làm trơn trường gradient ứng suất - biến dạng mà không cần bổ sung nút phụ, không làm tăng kích thước ma trận độ cứng tổng thể và áp đặt điều kiện biên trực tiếp như FEM cổ điển.

+-----------------------------------------------------------------------------------+
|                        KHUNG PHÂN TÍCH TỔNG HỢP CFEM                              |
+-----------------------------------------------------------------------------------+
                                         |
     +-----------------------------------+-----------------------------------+
     |                                   |                                   |
     v                                   v                                   v
+-----------------------+   +-----------------------+   +-----------------------+
|  THỦ TỤC NỘI SUY CIP  |   |    CƠ HỌC PHÁ HỦY     |   |  PHI TUYẾN HÌNH HỌC   |
|      (2D & 3D)        |   |       (X-CFEM)        |   |  & TÍCH PHÂN THAY THẾ |
+-----------------------+   +-----------------------+   +-----------------------+
| - Phần tử CT3, CQ4    |   | - Hàm làm giàu Heaviside| | - Biến dạng lớn       |
| - Phần tử CTH4, CHH8  |   | - Hàm dốc (Ramp func) |   |   Green-Lagrange      |
| - Miền hỗ trợ S_I     |   | - Tích phân J tương tác | | - Thuật toán Arc-Length|
| - Giữ Kronecker-delta |   | - Tiêu chuẩn nứt MTS  |   | - Tích phân EM & EF   |
| - DOF không đổi       |   | - Vật liệu FGM dị tính|   | - Trị khóa thể tích   |
+-----------------------+   +-----------------------+   +-----------------------+
     |                                   |                                   |
     +-----------------------------------+-----------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
|     ỨNG DỤNG ĐA QUY MÔ: TĨNH, ĐỘNG, LAN TRUYỀN VẾT NỨT, DAO ĐỘNG RIÊNG, KHUNG VỎ  |
+-----------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liên ngành giữa 3 lý thuyết cơ bản:

  1. Lý thuyết đàn hồi dị tính FGM: Mô tả sự biến thiên liên tục của các thông số vật liệu $E(\mathbf{x}), \nu(\mathbf{x}), \rho(\mathbf{x})$ theo không gian: $$E(y) = E_0 e^{\beta_E y}, \quad \beta_E = \frac{1}{W} \ln\left(\frac{E_1}{E_0}\right)$$
  2. Lý thuyết phá hủy đàn hồi tương tác (Interaction Integral): Sử dụng tích phân $J$ tương tác trên miền hỗ trợ mở rộng của phần tử XCQ4 để tách biệt chính xác hệ số cường độ ứng suất Mode I ($K_I$) và Mode II ($K_{II}$), kết hợp tiêu chuẩn ứng suất tiếp lớn nhất (Maximum Tangential Stress - MTS) để xác định góc lệch nứt $\theta_0$: $$\theta_0 = 2 \arctan \left( \frac{K_I - \sqrt{K_I^2 + 8K_{II}^2}}{4K_{II}} \right)$$
  3. Lý thuyết tích phân thay thế tối ưu: Đề xuất kỹ thuật Element Mid-point (EM) và Element Mid-face (EF) cho phần tử 2D và 3D (3D-EM). Thay vì sử dụng lưới tích phân Gauss $2\times2\times2$ (8 điểm) trong phần tử lục diện CHH8, mô hình EM/EF giảm thiểu số điểm tích phân xuống các điểm giữa mặt và tâm phần tử, triệt tiêu hoàn toàn hiện tượng khóa số học (locking phenomena) trong vật liệu gần như không nén được ($\nu \to 0.5$).

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 thực chứng định lượng (positivism/quantitative paradigm), kết hợp giữa giải tích toán học giải tích cấu trúc phần tử và mô phỏng số học (numerical experiments). Toàn bộ hệ thống thuật toán và cấu trúc phần tử được phát triển độc lập và thực thi trên nền tảng phần mềm MATLAB, được kiểm chứng chéo (benchmarking) với các nghiệm giải tích chính xác và phần mềm thương mại tiêu chuẩn công nghiệp ABAQUS (sử dụng lưới siêu mịn phần tử bậc cao 20 nút HH20).

+-----------------------------------------------------------------------------------+
|                        QUY TRÌNH THỰC THI THUẬT TOÁN                              |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [Rời rạc hóa miền tính toán] ---> [Xác định tập lân cận S_I và trọng số w_e]      |
|                                                  |                                |
|                                                  v                                |
|  [Áp dụng tích phân EM/EF] <------- [Tính đạo hàm trung bình và hàm phụ trợ phi_I] |
|             |                                                                     |
|             v                                                                     |
|  [Thiết lập ma trận K, M] --------> [Giải hệ phương trình / Thuật toán Arc-Length]|
|                                                  |                                |
|                                                  v                                |
|  [Hậu xử lý trường ứng suất liên tục] <-- [Trích xuất SIF / Khảo sát snap-through]|
|                                                                                   |
+-----------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình mô phỏng số tuân thủ nghiêm ngặt các tiêu chuẩn kiểm định độ tin cậy và tính hợp thức:

  • Kiểm định tính hội tụ lưới (Mesh Convergence Analysis): Thực hiện phân tích trên đa cấp độ lưới (từ lưới thô 480 phần tử đến lưới siêu mịn 6720 phần tử đối với 3D; từ 52 nút đến 2425 nút đối với 2D).
  • Đánh giá sai số năng lượng biến dạng toàn cục (Strain Energy Error Norm): $$\text{Error}(E) = \frac{1}{2} \int_\Omega (\boldsymbol{\sigma} - \boldsymbol{\sigma}{\text{ref}})^T (\boldsymbol{\epsilon} - \boldsymbol{\epsilon}{\text{ref}}) d\Omega$$
  • Kiểm định độ nhạy của bước thời gian và mật độ lưới: Trong phân tích nứt động, khảo sát phản ứng động học dưới tải trọng hàm bước (step loading) và tải trọng điều hòa dạng sin (sinusoidal loading) với 4 bước thời gian khác nhau để loại trừ hiện tượng phân tán sóng số học (numerical wave dispersion).

Data và phân tích

Các bộ dữ liệu mô phỏng đại diện phản ánh sự ưu việt vượt trội của CFEM:

Dữ liệu so sánh hiệu quả tính toán bài toán dầm công son 2D (Sai số năng lượng và thời gian thực thi):

  • Lưới $12\times3$ (52 nút, 104 DOF): CQ4 đạt sai số $2.004%$ trong thời gian $1.33\text{ s}$, trong khi Q4 có sai số $6.85%$ thời gian $0.072\text{ s}$.
  • Lưới $24\times6$ (175 nút): CQ4 đạt sai số $0.300%$ chỉ với $2.68\text{ s}$.
  • Lưới $96\times24$ (2425 nút): Q4 phải mất tới $5.58\text{ s}$ và số lượng nút gấp gần 14 lần mới đạt được mức sai số tương đương ($0.38%$).

Dữ liệu phân tích dầm 3D tiết diện chữ T ($E = 210\text{ GPa}, \nu = 0.3$ chịu lực phân bố $10\text{ N/mm}^2$):

  • Phần tử CHH8 với lưới 960 phần tử tiêu tốn $141\text{ s}$ để đạt sai số năng lượng biến dạng chỉ $0.53%$.
  • Phần tử chuẩn HH8 với lưới 2240 phần tử tiêu tốn $151\text{ s}$ nhưng vẫn chịu mức sai số lớn hơn ($0.75%$), và ở mức lưới 480 phần tử sai số của HH8 lên tới $2.72%$.

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

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

  1. Khả năng làm mịn tự nhiên trường ứng suất không gián đoạn: Như luận án đã nhấn mạnh: "Cách tiếp cận với thủ tục CIP 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". Trên mô hình thanh FGM TiB/Ti, sai số ứng suất chính $\sigma_{xx}$ tại điểm B' thu được từ phần tử CQ4 chỉ là $0.07%$ so với nghiệm chuẩn Q8 của Kim & Paulino (2002), trong khi phần tử chuẩn Q4 chịu sai số lên tới $1.15%$.
  2. Hiện tượng tái phân bố ứng suất cực đại trong kết cấu FGM: Kết quả mô phỏng số chỉ ra rằng, đối với vật liệu đồng nhất, ứng suất tập trung cực đại xuất hiện tại đáy rãnh góc nhọn (điểm A'), trong khi đối với kết cấu biến tính FGM, điểm cực đại dịch chuyển ra biên ngoài (điểm B'). Phát hiện này cung cấp cơ sở cơ học then chốt cho thiết kế chống phá hủy và bảo trì công trình kỹ thuật.
  3. Độ chính xác vượt bậc trong phân tích dao động riêng hình trụ rỗng lệch tâm: Khảo sát 4 mode dao động đầu tiên có tần số khác không trên kết cấu hình trụ khoét rỗng với độ lệch tâm $e_c = 80%$, sai số lớn nhất của phần tử CHH8 (lưới 3072 phần tử) so với nghiệm giải tích của Leissa et al. chỉ là $0.52%$, vượt trội hoàn toàn so với sai số xấp xỉ $1.1%$ của phần tử HH8 cùng mật độ lưới.
  4. Mô phỏng chính xác quỹ đạo lan truyền vết nứt động và bài toán uốn 3 điểm/4 điểm: Mô hình phần tử mở rộng XCQ4 kết hợp hàm dốc tái hiện hoàn hảo đường nứt thực nghiệm trên mẫu FGM mà không xảy ra hiện tượng phân nhánh giả tạo hay phụ thuộc hướng lưới.
  5. Triệt tiêu hiện tượng khóa thể tích trong bài toán phi tuyến hình học Cook's Membrane: Khi hệ số Poisson tiến tới giới hạn không nén được ($\nu = 0.4999$), phương pháp tích phân thay thế 3D-EM kết hợp CHH8 duy trì đường cong tải trọng - chuyển vị chuẩn xác và hội tụ ổn định qua điểm tới hạn snap-through bằng giải thuật Arc-length, trong khi tích phân Gauss truyền thống bị khóa cứng hoàn toàn.

Implications đa chiều

  • Ý nghĩa học thuật: Thiết lập nền tảng toán học tổng quát cho việc chuyển đổi mọi mã nguồn FEM hiện có sang hệ phần tử nội suy kép CIP mà không cần tái cấu trúc cơ sở dữ liệu nút.
  • Đổi mới phương pháp luận: Cung cấp mô hình tích phân thay thế 3D-EM và EF giúp tiết kiệm từ $30%$ đến $50%$ chi phí tích phân số trong các bài toán biến dạng lớn và tiếp xúc phi tuyến.
  • Ứng dụng thực tiễn công nghiệp: Cho phép các kỹ sư hàng không vũ trụ và năng lượng hạt nhân mô phỏng chính xác ứng xử nhiệt - cơ và độ bền mỏi của các tấm vỏ composite/FGM chịu tải trọng khắc nghiệt.
  • Khuyến nghị chính sách công nghệ: Định hướng phát triển các gói phần mềm tính toán kỹ thuật nội địa tự chủ (CAE software), giảm thiểu sự phụ thuộc vào các giấy phép thương mại nước ngoài đắt đỏ.

Limitations và Future Research

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

  1. Phạm vi phi tuyến vật liệu: Nghiên cứu hiện tại tập trung chủ yếu vào trạng thái đàn hồi tuyến tính (Hookean elasticity) và phi tuyến hình học, chưa mở rộng cho mô hình vật liệu dẻo (elasto-plasticity) hoặc đàn - nhớt - dẻo biến dạng lớn.
  2. Quy mô bài toán nứt 3D: Mô hình phá hủy mở rộng X-CFEM mới chỉ được hoàn thiện cho bài toán phẳng 2D và dầm FGM, chưa giải quyết trọn vẹn bài toán lan truyền mặt nứt không gian 3D phức tạp với sự tương tác nhiều vết nứt cùng lúc.
  3. Chi phí khởi tạo hàm dạng CIP cục bộ: Do miền hỗ trợ $S_I$ lớn hơn FEM truyền thống, thời gian tính toán hàm dạng ban đầu tại mỗi phần tử cao hơn FEM cùng bậc (mặc dù tổng thời gian hội tụ toàn cục giảm).

Chương trình nghiên cứu tương lai (5-10 năm):

  • Phát triển phần tử vỏ nội suy kép (CIP-Shell) dựa trên lý thuyết biến dạng trượt bậc cao (HSDT) cho kết cấu không gian phức tạp.
  • Kết hợp CFEM với thuật toán trí tuệ nhân tạo (Physics-Informed Neural Networks - PINNs) để tối ưu hóa hàm dạng tự thích nghi.
  • Mở rộng thuật toán giải bài toán tương tác đa trường vật lý: Cơ - Nhiệt - Điện - Từ trong vật liệu áp điện FGM và vật liệu biến tính nano.

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

Luận án khẳng định vị thế dẫn đầu trong chuyên ngành Cơ học tính toán tại Việt Nam với các tác động đa tầng:

  • Tác động học thuật: Các công trình công bố trực tiếp từ luận án trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Acta Mechanica, Composite Structures) đã khẳng định năng lực hội nhập khoa học quốc tế của nhóm nghiên cứu GS. TSKH Nguyễn Đình Đức và PGS. TS Bùi Quốc Tính.
  • Tái thiết quy trình R&D công nghiệp: Tạo tiền đề cho việc tích hợp mô-đun CFEM vào các hệ thống thiết kế hỗ trợ máy tính (CAD/CAE), tối ưu hóa thiết kế vật liệu FGM trong ngành công nghiệp quốc phòng, vỏ động cơ tên lửa và tuabin khí thế hệ mới.
  • Chuyển giao và đào tạo: Hình thành khung bài giảng chuyên đề cao học và nghiên cứu sinh về phương pháp số tiên tiến trong Cơ học vật rắn biến dạng tại Đại học Quốc gia Hà Nội.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận phương pháp luận toán học chặt chẽ về kỹ thuật nội suy hàm dạng và giải thuật xử lý kỳ dị vết nứt.
  • Các nhà khoa học cơ học tính toán: Sử dụng công thức toán học mở của CIP và mô hình tích phân 3D-EM làm công cụ phát triển các lớp phần tử hữu hạn mới.
  • Kỹ sư thiết kế kết cấu công trình và hàng không: Nắm bắt quy luật phân bố lại ứng suất của vật liệu FGM để thiết kế các chi tiết chịu lực tối ưu, giảm trọng lượng kết cấu nhưng tăng độ bền phá hủy.
  • Các nhà phát triển phần mềm mô phỏng (CAE Developers): Khai thác mã nguồn thuật toán tối ưu hóa thời gian tính toán và giải quyết triệt để lỗi khóa thể tích trong phân tích kết cấu cao tầng chịu tải ngang và vật liệu cao su/polymer.

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ì?
Đó là việc thiết lập và chứng minh thành công công thức toán học tổng quát cho hàm dạng nội suy liên tiếp CIP trong không gian 3 chiều (CTH4, CHH8). Công thức này nâng bậc trơn của trường chuyển vị và đạo hàm ứng suất lên mức $C^0-C^1$ liên tục giữa các biên phần tử, bảo toàn nguyên vẹn tính chất Kronecker-delta ($\phi_I(\mathbf{x}p) = \delta{Ip}$) và không làm gia tăng bất kỳ bậc tự do (DOF) nào trên ma trận độ cứng tổng thể.

2. Đổi mới phương pháp luận của CFEM so với hai nghiên cứu quốc tế điển hình là gì?
So với phương pháp phần tử làm mịn S-FEM của G. R. Liu (cần chia thêm các miền tế bào làm mịn smoothing domains) và phương pháp đẳng hình học IGA của Hughes (dùng hàm NURBS làm mất tính chất Kronecker-delta), CFEM can thiệp trực tiếp vào đạo hàm trung bình tại nút $\overline{\mathbf{N}}_{I,w}$, cho phép nhúng trực tiếp vào cấu trúc dữ liệu của các mã nguồn FEM cổ điển mà không cần kỹ thuật gán điều kiện biên phức tạp.

3. Phát hiện thực nghiệm số gây bất ngờ nhất được hỗ trợ bởi dữ liệu là gì?
Trong phân tích kéo tấm FGM (TiB/Ti) có khía rãnh, phần tử CQ4 chỉ ra rằng ứng suất cực đại $\sigma_{\max}$ không nằm tại đáy rãnh góc nhọn (điểm A') như vật liệu đồng nhất mà dịch chuyển hoàn toàn ra vị trí biên ngoài (điểm B') với sai số cực nhỏ chỉ $0.07%$ so với nghiệm bậc cao Q8, chứng minh cơ chế tự phân phối lại ứng suất ưu việt của cấu trúc vật liệu biến tính chức năng.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?
Toàn bộ hệ thống phương trình hàm dạng, tọa độ điểm tích phân, trọng số cho mô hình tích phân EM/EF/3D-EM, các tham số vật liệu và điều kiện biên của từng bài toán thử nghiệm (dầm công son, màng Cook, dầm chữ T, mẫu nứt TPB) đều được tường minh chi tiết dưới dạng công thức giải tích và thuật toán thực thi bằng mã lệnh MATLAB.

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?
Chương trình tập trung vào 3 mũi nhọn: (1) Mở rộng CFEM cho bài toán nứt không gian 3D tương tác đa pha; (2) Phát triển các phần tử tấm/vỏ phi tuyến hình học bậc cao chịu tải nhiệt - cơ khắc nghiệt; (3) Tích hợp học máy (Machine Learning) để tối ưu hóa lưới phần tử và nhận dạng vết nứt tự động thời gian thực.

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 5 đóng góp mang tính đột phá khoa học:

  1. Hoàn thiện lý thuyết nội suy kép CIP cho không gian 2D và 3D: Xây dựng thành công hệ phần tử CT3, CQ4, CTH4, CHH8 với trường đạo hàm liên tục trơn mịn và giữ nguyên tính chất Kronecker-delta.
  2. Đột phá phương pháp luận cơ học phá hủy với XCQ4: Tích hợp thành công kỹ thuật làm giàu trường chuyển vị hàm dốc trong vật liệu biến tính FGM, loại bỏ dao động số tại đỉnh nứt và nâng cao độ chính xác của các hệ số SIF/DSIF.
  3. Phát triển thành công mô hình tích phân thay thế 3D-EM và EF: Giảm thiểu đáng kể số điểm tích phân và thời gian tính toán, giải quyết triệt để hiện tượng khóa thể tích khi phân tích phi tuyến hình học vật liệu gần như không nén được ($\nu \to 0.5$).
  4. Chứng minh tính ưu việt vượt bậc về hiệu năng tính toán: Tiết kiệm hơn $50%$ thời gian tính toán để đạt cùng cấp độ chính xác $0.3%$ sai số so với FEM truyền thống trên cùng cấu hình bài toán.
  5. Mở ra hướng nghiên cứu mới về vật liệu tiên tiến: Đặt nền móng vững chắc cho việc mô phỏng số chính xác các bài toán dao động riêng, phân tích phi tuyến và lan truyền hư hại trong các kết cấu composite/FGM hiện đại phục vụ kỹ thuật hàng không vũ trụ và công trình thế kỷ.