Tổng quan về luận án

Sự phát triển mạnh mẽ của khoa học vật liệu và kỹ thuật tính toán đã thúc đẩy sự ra đời của các kết cấu thông minh thế hệ mới, đặc biệt là các tấm composite đa tầng, vật liệu cơ tính biến thiên (FGM) có độ xốp và được gia cường bởi các tấm nano graphene (GPLs) tích hợp các lớp áp điện (Piezoelectric materials). Tuy nhiên, việc mô phỏng ứng xử tĩnh tuyến tính, phi tuyến hình học, dao động tự do và điều khiển chủ động các kết cấu này đang đối mặt với rào cản nghiêm trọng từ phương pháp phần tử hữu hạn truyền thống (FEM). Trong quy trình thiết kế và phân tích cơ học truyền thống (CAD/CAE), giai đoạn tạo lưới (meshing) và xử lý hình học thô chiếm tới khoảng 80% tổng thời gian phân tích kỹ thuật ("meshing and clean-up steps, which account approximately 80% of the total analysis time of a problem"). Quá trình xấp xỉ lưới trong FEM làm suy giảm độ chính xác hình học thực tế, đặc biệt đối với các kết cấu có biên cong phức tạp như tấm tròn, tấm elip, tấm vành khăn hay tấm có lỗ khoét hình học phức tạp.

+-------------------------------------------------------------------------------+
|                    KHUNG PHÂN TÍCH TỔNG QUAN CỦA LUẬN ÁN                      |
+-------------------------------------------------------------------------------+
|  1. MÔ HÌNH HÌNH HỌC & TÍNH TOÁN:                                             |
|     CAD / NURBS Functions  ==>  Bézier Extraction  ==>  Bernstein Elements    |
|                                                                               |
|  2. MÔ HÌNH CƠ HỌC TẤM:                                                       |
|     Lý thuyết UHSDT (Thư giãn ứng suất biên, không cần hệ số cắt SCFs)        |
|     Thế điện phi tuyến Maxwell: Phi(z) = Tuyến tính + Hàm điều hòa Cosine     |
|                                                                               |
|  3. ĐỐI TƯỢNG VẬT LIỆU CAO CẤP:                                               |
|     - Laminated Composite & Piezoelectric Smart Plates                        |
|     - PFGP-GPLs (Tấm xốp FGM gia cường 1% wt GPLs + Tích hợp áp điện)         |
|     - FGPMP (Tấm vật liệu áp điện FGM xốp PZT-4/PZT-5H)                       |
|                                                                               |
|  4. PHÂN TÍCH & ĐIỀU KHIỂN CHỦ ĐỘNG:                                          |
|     Tĩnh phi tuyến (von Kármán) + Dao động tự do + Quá độ (Step/Blast Loads)  |
|     Vòng điều khiển kín: Độ lợi dịch chuyển Gd & Độ lợi vận tốc Gv            |
+-------------------------------------------------------------------------------+

Nghiên cứu của NCS. Nguyễn Thị Bích Liễu, dưới sự hướng dẫn khoa học của GS. Nguyễn Xuân Hùng và PGS. Đặng Thiện Ngôn, đã giải quyết khoảng cách học thuật (research gap) cốt lõi: thiết lập khung tính toán Isogeometric Analysis (IGA) dựa trên toán tử trích xuất Bézier (Bézier extraction) kết hợp lý thuyết biến dạng cắt bậc cao không ràng buộc tổng quát (Generalized Unconstrained Higher-Order Shear Deformation Theory - UHSDT) để phân tích cơ học và điều khiển chủ động đáp ứng tĩnh/động lực học tuyến tính và phi tuyến cho các kết cấu tấm composite thông minh, tấm PFGP-GPLs và tấm FGPMP.

Mục tiêu khoa học và hệ thống giả thuyết nghiên cứu được xác lập đồng bộ:

  1. Câu hỏi nghiên cứu 1 ($RQ_1$): Làm thế nào để loại bỏ sai số chuyển đổi hình học CAD-CAE và tích hợp liền mạch hàm cơ sở NURBS trơn bậc cao $C^{p-1}$ vào cấu trúc mã nguồn phần tử hữu hạn tiêu chuẩn?
  2. Câu hỏi nghiên cứu 2 ($RQ_2$): Cơ chế tương tác giữa độ xốp ($e_0 \in [0, 0.4]$), quy luật phân tán hạt nano graphene (Patterns A, B, C) và trường điện thế ngoài ảnh hưởng như thế nào đến độ cứng uốn và ứng xử động lực học phi tuyến của kết cấu?
  3. Câu hỏi nghiên cứu 3 ($RQ_3$): Thuật toán điều khiển phản hồi vận tốc ($G_v$) và dịch chuyển ($G_d$) đạt hiệu quả tối ưu ra sao trong việc dập tắt dao động cưỡng bức dưới tải trọng động ngẫu nhiên và tải trọng sóng nổ (blast loading)?

Hệ thống giả thuyết bao gồm: ($H_1$) Kỹ thuật trích xuất Bézier cho phép bảo toàn 100% hình học CAD gốc đồng thời duy trì tính liên tục liên phần tử $C^{p-1}$, nâng cao tốc độ hội tụ vượt trội so với FEM truyền thống; ($H_2$) Việc bổ sung $1\text{ wt}%$ tấm nano GPLs phân bố tập trung tại bề mặt ngoài (Pattern A) có thể bù đắp hoàn toàn sự suy giảm độ cứng do hệ số rỗng $e_0 = 0.4$ gây ra trong tấm xốp FGM.

Phạm vi nghiên cứu bao gồm việc phân tích tĩnh, dao động tự do và động lực học quá độ phi tuyến trên 4 hệ vật liệu phức tạp, khảo sát đa dạng điều kiện biên (CCCC, SSSS, CFFF, CCFF) và cấu hình hình học từ tấm vuông, chữ nhật đến các biên phi kinh điển (tấm elip, tấm tròn, vành khăn, tấm có lỗ khoét phức tạp), tạo lập bộ dữ liệu chuẩn mực cho cơ học tính toán vật liệu tiên tiến.


Literature Review và Positioning

Lịch sử phát triển của hình học tính toán ghi nhận các mốc đột phá từ đa thức Bernstein của Pierre Bézier (1966), B-splines của Riesenfeld (1972) đến Non-Uniform Rational B-splines (NURBS) của Versprille (1975), trở thành tiêu chuẩn công nghiệp trong CAD. Năm 2005, Hughes, Cottrell và Bazilevs đã khởi xướng phương pháp Phân tích đẳng hình học (Isogeometric Analysis - IGA), đề xuất nguyên lý mang tính cách mạng: sử dụng trực tiếp các hàm cơ sở mô tả hình học trong CAD (B-splines/NURBS) làm hàm xấp xỉ trường chuyển vị trong phân tích cơ học ("instead of converting one system to another... one should substitute one system for the other so that the conversion is no longer needed").

Trong lĩnh vực cơ học kết cấu tấm, các dòng lý thuyết chính bao gồm:

  1. Lý thuyết tấm cổ điển (CPT/CLPT): Bỏ qua biến dạng trượt ngang, chỉ phù hợp cho tấm rất mỏng ($a/h \ge 100$).
  2. Lý thuyết biến dạng cắt bậc nhất (FSDT): Giả thiết góc quay không đổi theo chiều dày, đòi hỏi hệ số hiệu chỉnh cắt (Shear Correction Factors - SCFs) phụ thuộc phức tạp vào hình học và điều kiện biên.
  3. Lý thuyết biến dạng cắt bậc cao (HSDT) và lý thuyết từng lớp (Layer-wise Theory - LWT): Mô hình hóa ứng xử phi tuyến của ứng suất trượt ngang qua chiều dày.

Cuộc tranh luận học thuật kéo dài giữa hai trường phái: LWT cho kết quả ứng suất cục bộ chính xác cao nhưng làm bùng nổ số bậc tự do (DOFs) theo số lớp xếp; ngược lại, mô hình lớp đơn tương đương (Equivalent Single Layer - ESL) dựa trên HSDT giảm thiểu biến số nhưng thường đòi hỏi đạo hàm bậc cao $C^1$, gây khó khăn nghiêm trọng cho phần tử hữu hạn $C^0$ tiêu chuẩn.

+-------------------------------------------------------------------------------+
|             BẢNG ĐỐI CHIẾU TIẾP CẬN VỚI CÁC CÔNG TRÌNH QUỐC TẾ                |
+-------------------------------------------------------------------------------+
| Nghiên cứu & Tác giả | Phương pháp / Mô hình        | Giới hạn chính          |
+-------------------------------------------------------------------------------+
| Mitchell & Reddy     | Giải tích Navier + TSDT      | Chỉ giải được tấm chữ   |
| (1995)               | cho composite áp điện        | nhật biên tựa đơn SSSS  |
+-------------------------------------------------------------------------------+
| Hwang & Park         | Phần tử DKQ + FSDT           | Xấp xỉ C0, mất độ trơn  |
| (2000)               | Điều khiển vận tốc âm        | hình học ở biên cong    |
+-------------------------------------------------------------------------------+
| Kitipornchai et al.  | Phương pháp Ritz +           | Giới hạn cho kết cấu dầm|
| (2017)               | Timoshenko cho FG xốp        | 1D, chưa xét tấm 2D/3D  |
+-------------------------------------------------------------------------------+
| Barati et al.        | Giải tích 4 biến             | Không xử lý được hình   |
| (2017)               | Tấm FGPMP áp điện xốp        | học có lỗ khoét phức tạp|
+-------------------------------------------------------------------------------+
| Luận án NCS.         | IGA Bézier Extraction +      | Khắc phục toàn bộ các   |
| Nguyễn Thị Bích Liễu | UHSDT + Tích hợp áp điện     | hạn chế trên: Đa hình   |
| (2019)               | điều khiển chủ động          | học, phi tuyến, C(p-1)  |
+-------------------------------------------------------------------------------+

Luận án định vị tiên phong bằng cách giải quyết triệt để khoảng trống nghiên cứu: kết hợp IGA trên nền tảng trích xuất Bézier với lý thuyết UHSDT, mở rộng phân tích phi tuyến hình học von Kármán và áp dụng thuật toán điều khiển kín cho kết cấu tấm xốp FGM gia cường hạt nano GPLs tích hợp áp điện với hình học phức tạp đa vùng (multi-patch).


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

                      +-----------------------------+
                      |   LÝ THUYẾT ĐÓNG GÓP UHSDT   |
                      +--------------+--------------+
                                     |
             +-----------------------+-----------------------+
             |                                               |
             v                                               v
+---------------------------+                   +---------------------------+
| Thư giãn ứng suất trượt   |                   | Đề xuất 2 hàm phân bố f(z)|
| sigma_xz = sigma_yz = 0   |                   | tối ưu, loại bỏ hoàn toàn |
| tại mặt trên/dưới tấm     |                   | hệ số hiệu chỉnh cắt SCFs |
+---------------------------+                   +---------------------------+

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

Luận án phát triển khung lý thuyết biến dạng cắt bậc cao không ràng buộc tổng quát (Generalized UHSDT). Khác với lý thuyết TSDT truyền thống của Reddy vốn áp đặt điều kiện ứng suất cắt bằng không tại hai mặt tự do $\sigma_{xz}(\pm h/2) = \sigma_{yz}(\pm h/2) = 0$, lý thuyết UHSDT tổng quát hóa trường chuyển vị thông qua hàm phân bố $f(z)$ tùy chọn:

$$u(x,y,z,t) = u_0(x,y,t) - z \frac{\partial w_0}{\partial x} + f(z)\theta_x(x,y,t)$$

$$v(x,y,z,t) = v_0(x,y,t) - z \frac{\partial w_0}{\partial y} + f(z)\theta_y(x,y,t)$$

$$w(x,y,z,t) = w_0(x,y,t)$$

Đóng góp lý thuyết then chốt nằm ở việc đề xuất 2 dạng hàm phân bố trượt $f(z)$ mới có tính chất toán học vượt trội, loại bỏ nhu cầu sử dụng hệ số hiệu chỉnh cắt (SCFs) mà vẫn phản ánh chính xác phân bố phi tuyến thực tế của ứng suất trượt ngang. Đồng thời, mô hình thế điện áp điện phi tuyến $\Phi(z,t)$ được xây dựng bằng sự kết hợp giữa hàm tuyến tính và hàm điều hòa cosine:

$$\Phi(z,t) = -\left(\frac{2z}{h}\right) V_0(t) + \cos\left(\frac{\pi z}{h}\right) \tilde{\phi}(x,y,t)$$

Công thức này thỏa thuận hoàn toàn với phương trình Maxwell tĩnh điện trong điều kiện quasi-static mà không làm tăng biến số tính toán.

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

Khung phân tích tích hợp nhất quán ba hệ thống lý thuyết nền tảng:

  1. Lý thuyết môi trường liên tục cơ - điện ghép nối (Electro-mechanical continuum mechanics): Tích hợp hiệu ứng áp điện trực tiếp và nghịch đảo thông qua hệ phương trình cấu tạo nhiệt - cơ - điện.
  2. Cơ học vi mô vật liệu xốp FGM gia cường nano (Micromechanics of Porous FGM-GPLs): Sử dụng mô hình Halpin-Tsai cải tiến kết hợp định luật hỗn hợp mở rộng để xác định mô đun đàn hồi $E(z)$, khối lượng riêng $\rho(z)$ và hệ số Poisson $\nu(z)$ thay đổi liên tục qua chiều dày theo 3 dạng phân bố lỗ rỗng và 3 quy luật phân tán GPLs.
  3. Toán tử trích xuất Bézier (Bézier extraction operator): Phân rã toán tử chiếu toàn cục của NURBS thành tích đại số giữa ma trận trích xuất $\mathbf{C}^e$ và hàm cơ sở đa thức Bernstein cục bộ:

$$\mathbf{N}^e(\xi) = \mathbf{C}^e \mathbf{B}(\xi)$$

Khung phân tích này phá vỡ rào cản tính toán đơn miền, cho phép thiết lập điều kiện liên tục $C^{p-1}$ liên phần tử và liên kết trơn tru các miền tham số đa mảng (multi-patch) cho các kết cấu phức tạp.


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

+-------------------------------------------------------------------------------+
|                       QUY TRÌNH TÍNH TOÁN IGA - BÉZIER                        |
+-------------------------------------------------------------------------------+
|  1. THIẾT KẾ HÌNH HỌC CAD (NURBS Domain)                                      |
|     Xác lập Knot Vector: Xi = {xi_1, xi_2, ..., xi_{n+p+1}}                   |
|     Tập điểm kiểm soát P_{i,j} & Trọng số w_{i,j}                             |
|                                                                               |
|  2. TRÍCH XUẤT BÉZIER (Bézier Extraction)                                     |
|     Knot Insertion liên tiếp ==> N^e(xi) = C^e * B(xi)                        |
|     Cấu trúc dữ liệu tương thích hoàn toàn Element-by-Element của FEM         |
|                                                                               |
|  3. TĂNG CƯỜNG HÀM CƠ SỞ (k-refinement)                                       |
|     Nâng bậc hàm cơ sở p (p-refinement) TRƯỚC ==> Chèn nút h (h-refinement)  |
|     Đạt độ trơn tối đa C^{p-1} với số bậc tự do tối ưu nhất                   |
|                                                                               |
|  4. TÍCH PHÂN SỐ GAUSS QUADRATURE & GIẢI HỆ PHƯƠNG TRÌNH                      |
|     Thiết lập ma trận: M u'' + (C + C_act) u' + (K + K_act) u = F_ext        |
|     Vòng lặp Newmark / Newton-Raphson cho phân tích phi tuyến tĩnh/động       |
+-------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Nghiên cứu vận hành trên hệ hình nhận thức luận thực chứng (Positivist Paradigm), dựa trên cơ sở mô hình hóa giải tích số học chính xác và xác thực chéo (benchmarking). Điểm cốt lõi về mặt phương pháp luận là khai thác kỹ thuật $k$-refinement độc nhất vô nhị của IGA. Khác với FEM (chỉ có $h$-refinement chia nhỏ lưới hoặc $p$-refinement tăng bậc đa thức tạo liên kết $C^0$), $k$-refinement trong IGA thực hiện nâng bậc đa thức $p$ trước, sau đó mới tiến hành chèn nút $h$. Kết quả mang lại không gian xấp xỉ có bậc cao và độ trơn liên phần tử cực đại $C^{p-1}$, giảm triệt để sai số tán xạ sóng (dispersion error) và khóa trượt (shear locking).

Quy trình nghiên cứu rigorous

Quy trình tính toán được thực thi theo các giai đoạn nghiêm ngặt:

  1. Thiết lập miền tham số: Khởi tạo vector nút mở $\Xi = {\xi_1, \xi_2, ..., \xi_{n+p+1}}$ và lưới điểm kiểm soát (control net) $\mathbf{P}{i,j}$ cùng trọng số $w{i,j}$.
  2. Chuyển đổi Bézier: Thực hiện phép chèn nút lặp lại cho đến khi mỗi nút trong đạt bội số bằng bậc $p$. Phép phân rã Bézier cô lập từng phần tử nhưng bảo toàn ma trận liên kết trích xuất $\mathbf{C}^e$, biến đổi hàm cơ sở NURBS toàn cục thành các hàm Bernstein cục bộ.
  3. Thiết lập dạng biến phân yếu (Weak form): Áp dụng nguyên lý Hamilton mở rộng cho hệ liên hợp cơ - điện:

$$\int_{t_1}^{t_2} \left( \delta \mathcal{T} - \delta \mathcal{U} + \delta \mathcal{W}_{ext} \right) dt = 0$$

  1. Tích phân số học: Sử dụng phép cầu nghiệm Gauss (Gaussian quadrature) chuẩn hóa trên từng phần tử Bézier cha $\hat{\Omega}^e = [-1, 1]^2$.

Data và phân tích

Hệ phương trình động lực học tổng quát sau rời rạc hóa có dạng:

$$\mathbf{M}\ddot{\mathbf{u}} + \mathbf{C}\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} - \mathbf{K}{u\phi}\boldsymbol{\Phi} = \mathbf{F}{m}$$

$$\mathbf{K}{\phi u}\mathbf{u} + \mathbf{K}{\phi\phi}\boldsymbol{\Phi} = \mathbf{Q}_{e}$$

Trong đó ma trận độ cứng $\mathbf{K}$, ma trận khối lượng $\mathbf{M}$, ma trận ghép nối áp điện $\mathbf{K}{u\phi}$ và ma trận điện dung $\mathbf{K}{\phi\phi}$ được tích hợp chính xác. Thuật toán điều khiển chủ động đưa vào các hệ số khuếch đại hồi tiếp vận tốc $G_v$ và dịch chuyển $G_d$:

$$\mathbf{V}_a(t) = -G_d \mathbf{w}(t) - G_v \dot{\mathbf{w}}(t)$$

Tạo ra ma trận giảm chấn chủ động $\mathbf{C}{active} = G_v \mathbf{K}{u\phi} \mathbf{K}{\phi\phi}^{-1} \mathbf{K}{\phi u}$ và ma trận độ cứng chủ động $\mathbf{K}{active} = G_d \mathbf{K}{u\phi} \mathbf{K}{\phi\phi}^{-1} \mathbf{K}{\phi u}$. Thuật toán tích phân Newmark kết hợp lặp Newton-Raphson được thiết lập để xử lý các bài toán phi tuyến hình học von Kármán.


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

+-------------------------------------------------------------------------------+
|              TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CỦA NGHIÊN CỨU                   |
+-------------------------------------------------------------------------------+
| 1. HIỆU ỨNG GIA CƯỜNG GPLs (PATTERN A TỐI ƯU):                                |
|    Chỉ 1% wt GPLs phân tán tập trung tại bề mặt ngoài giúp khôi phục hoàn toàn|
|    độ cứng bị suy giảm bởi 40% độ xốp (e0 = 0.4), vượt trội hơn Pattern B & C.|
|                                                                               |
| 2. ĐIỀU KHIỂN CHỦ ĐỘNG TRIỆT ĐỂ (Gd & Gv):                                    |
|    - Độ lợi Gd làm giảm trực tiếp độ võng uốn tĩnh tuyến tính và phi tuyến.   |
|    - Độ lợi Gv dập tắt dao động quá độ cưỡng bức cực nhanh dưới tải sóng nổ.  |
|                                                                               |
| 3. ĐIỀU BIẾN TẦN SỐ BẰNG ĐIỆN THẾ NGOÀI V0:                                  |
|    Điện áp dương (+V0) làm mềm kết cấu (giảm tần số), điện áp âm (-V0) làm   |
|    tăng độ cứng tương đương; cho phép 'tune' tần số cộng hưởng tức thời.      |
|                                                                               |
| 4. ĐỘ CHÍNH XÁC VƯỢT TRỘI TRÊN BIÊN CONG & LỖ KHOÉT PHỨC TẠP:                |
|    IGA Bézier hội tụ nhanh hơn 40-60% so với FEM bậc 1, giải quyết chính xác  |
|    tấm vành khăn, tấm elip, tấm có lỗ khoét phức tạp mà không bị méo lưới.   |
+-------------------------------------------------------------------------------+

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

  • Quy luật phân tán GPLs tối ưu (Pattern A): Kết quả mô phỏng số chỉ ra rằng sự phân tán của các tấm nano graphene theo mẫu A (tập trung mật độ cao nhất ở hai bề mặt trên và dưới của tấm) mang lại hiệu quả gia cường vượt bậc. Cụ thể, đối với tấm cantilever PFGP-GPLs chịu tải trọng phân bố đều, khi đưa hàm lượng $1\text{ wt}%$ GPLs vào nền xốp có hệ số rỗng $e_0 = 0.2$, độ võng đỉnh giảm tới hơn $65%$ so với trường hợp không có GPLs. Cấu trúc vật liệu xốp gia cường GPLs kết hợp ưu thế giảm trọng lượng bản thân mà không làm suy giảm độ cứng tổng thể.
Mặt trên:   =========================== (Mật độ GPLs cao nhất - Pattern A)
Lõi giữa:   . . . . [ Lỗ xốp e0 ] . . . (Mật độ GPLs thấp, xốp nhẹ)
Mặt dưới:   =========================== (Mật độ GPLs cao nhất - Pattern A)
==> Tăng tối đa mô-men kháng uốn I_xx với cùng một khối lượng vật liệu!
  • Hiệu năng dập tắt dao động của điều khiển phản hồi ($G_d, G_v$): Khảo sát phản ứng quá độ dưới các dạng tải trọng động (tải điều hòa, tải bước nhảy Step load, tải xung tam giác và tải sóng nổ Blast load) chứng minh rằng hệ số hồi tiếp vận tốc $G_v$ đóng vai trò quyết định trong việc dập tắt biên độ dao động. Khi tăng $G_v$, thời gian suy giảm dao động về trạng thái cân bằng giảm tới 75-80% mà không gây hiện tượng mất ổn định số học.
  • Hiệu ứng điều biến tần số dao động bằng điện trường ngoài ($V_0$): Đối với tấm FGPMP làm từ PZT-4/PZT-5H, điện áp phân cực $V_0$ tác động trực tiếp lên trường ứng suất ban đầu. Điện áp dương tạo ra lực nén trong mặt phẳng làm giảm tần số dao động riêng không thứ nguyên $\bar{\omega}$, trong khi điện áp âm tạo hiệu ứng căng giúp tăng tần số cơ bản, mở ra khả năng chủ động tránh miền tần số cộng hưởng nguy hiểm trong thời gian thực.
  • Sự vượt trội của IGA Bézier trên hình học phức tạp: Phân tích dao động tự do của tấm FGPMP hình vuông có lỗ khoét phức tạp (complicated hole/cutout) và tấm vành khăn (annular plate) sử dụng lưới phần tử Bézier bậc hai ($15\times 15$ phần tử, 336 đến 840 điểm kiểm soát) đạt độ chính xác tương đương nghiệm giải tích chuẩn mực, trong khi phần tử hữu hạn thương mại đòi hỏi hàng chục nghìn bậc tự do để khử sai số xấp xỉ biên cong.

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng thành công biên giới của cơ học tính toán vật liệu composite thông minh đa pha, cung cấp cơ sở dữ liệu số chuẩn mực (benchmark solutions) cho cộng đồng khoa học quốc tế.
  • Về mặt phương pháp luận: Khẳng định tính khả thi của việc nhúng trực tiếp toán tử Bézier extraction vào các phần mềm thương mại hiện nay (như ABAQUS, ANSYS) mà không cần viết lại cấu trúc giải ma trận phần tử.
  • Về mặt kỹ thuật thực tiễn: Tạo tiền đề thiết kế các cảm biến, cơ cấu chấp hành siêu nhẹ, cánh bay phản lực thông minh, cấu kiện hàng không vũ trụ và robot nano có khả năng tự cảm biến và tự triệt tiêu rung chấn.

Limitations và Future Research

Nghiên cứu duy trì tính khách quan học thuật bằng việc thừa nhận các giới hạn nội tại:

  1. Giả thiết liên kết lớp lý tưởng: Mô hình giả định độ dính bám hoàn hảo giữa các lớp composite, lớp áp điện và nền xốp FGM, chưa tính đến hiện tượng tách lớp (delamination) hoặc phá hủy ma trận polymer dưới tải trọng động chu kỳ dài.
  2. Xấp xỉ tĩnh điện cho trường thế: Áp dụng phương trình Maxwell trong điều kiện quasi-static, chưa bao hàm hoàn toàn hiệu ứng sóng điện từ lan truyền ở dải tần số siêu cao (GHz).
  3. Giới hạn kích thước micro/nano: Mô hình cơ học liên tục vĩ mô chưa tích hợp lý thuyết biến dạng gradient (strain gradient elasticity) hay ứng suất phi cục bộ (nonlocal elasticity) để mô tả các hiệu ứng phụ thuộc kích thước (size-dependent effects) của hạt nano GPLs ở thang đo dưới nanomet.
+-------------------------------------------------------------------------------+
|                       ĐỊNH HƯỚNG NGHIÊN CỨU TƯƠNG LAI                         |
+-------------------------------------------------------------------------------+
| 1. Cơ học phá hủy & Tách lớp: Phát triển XIGA mô phỏng vết nứt giao diện.      |
| 2. Môi trường ghép nối đa trường: Phân tích Nhiệt - Cơ - Điện - Từ (MCET).   |
| 3. Tối ưu hóa Topo (Topology Optimization): Tối ưu hóa vi cấu trúc xốp.      |
| 4. Kiểm chứng thực nghiệm: Chế tạo mẫu in 3D FGM-GPLs & đo đạc laser Doppler. |
+-------------------------------------------------------------------------------+

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

Luận án tạo ra tác động sâu rộng trên cả phương diện học thuật và chuyển giao công nghệ công nghiệp:

                  +-----------------------------------+
                  | TÁC ĐỘNG TOÀN DIỆN CỦA NGHIÊN CỨU |
                  +-----------------+-----------------+
                                    |
     +------------------------------+------------------------------+
     |                              |                              |
     v                              v                              v
+--------------------+    +--------------------+    +--------------------+
| KHOA HỌC HÀN LÂM   |    | CÔNG NGHIỆP CAO    |    | CHÍNH SÁCH VÀ MÔI  |
| Dự báo 200+ trích  |    | Hàng không vũ trụ, |    | TRƯỜNG: Giảm phát  |
| dẫn; Nền tảng cho  |    | MEMS/NEMS, Vỏ bọc  |    | thải, tiết kiệm    |
| nghiên cứu IGA     |    | vệ tinh siêu nhẹ   |    | 20-30% nhiên liệu  |
+--------------------+    +--------------------+    +--------------------+
  • Ảnh hưởng học thuật: Các công bố quốc tế trích xuất từ luận án trên các tạp chí ISI hàng đầu (như Composite Structures, Computer Methods in Applied Mechanics and Engineering) định hình xu hướng nghiên cứu mới trong việc ứng dụng Bézier IGA cho vật liệu nanocomposite.
  • Chuyển đổi công nghiệp: Cung cấp công cụ mô phỏng chuẩn xác cho ngành công nghiệp hàng không vũ trụ và sản xuất thiết bị y sinh, cho phép thiết kế vỏ thân máy bay, tấm pin mặt trời vệ tinh siêu nhẹ giảm từ 20-30% trọng lượng kết cấu nhưng tăng gấp đôi tuổi thọ mỏi nhờ cơ chế tự giảm chấn áp điện chủ động.
  • Lợi ích kinh tế - môi trường: Giảm tiêu hao nhiên liệu hàng không, hạ thấp chi phí thử nghiệm phá hủy mẫu thực nghiệm đắt đỏ nhờ các giải pháp tính toán mô phỏng tin cậy.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung toán học hoàn chỉnh, phương pháp lập trình IGA trích xuất Bézier rõ ràng để phát triển các bài toán phân tích vỏ mỏng, chất lưu - kết cấu (FSI) và cơ học rạn nứt.
  • Kỹ sư R&D Hàng không & Quốc phòng: Sử dụng các bảng dữ liệu chuẩn (benchmark tables) và thuật toán điều khiển $G_d, G_v$ để tích hợp vào hệ thống điều khiển kết cấu cánh bay thông minh thích ứng.
  • Nhà phát triển phần mềm CAE: Thu nhận giải thuật trích xuất Bézier để tích hợp các module phân tích IGA trơn bậc cao trực tiếp vào lõi phần tử hữu hạn hiện hành mà không cần tái cấu trúc phần mềm.

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 Lý thuyết biến dạng cắt bậc cao không ràng buộc tổng quát (Generalized UHSDT) với các hàm phân bố trượt $f(z)$ mới, kết hợp trường thế điện phi tuyến $\Phi(z)$ dạng điều hòa cosine. Mô hình này triệt tiêu hoàn toàn sự phụ thuộc vào hệ số hiệu chỉnh cắt (SCFs), thư giãn điều kiện biên ứng suất cắt bằng không mà vẫn bảo đảm tính chính xác tuyệt đối của phân bố ứng suất trượt qua chiều dày tấm composite và FGM đa tầng.

2. Điểm cải tiến phương pháp luận vượt trội so với các nghiên cứu tiền nhiệm?

So với nghiên cứu giải tích Navier của Mitchell & Reddy (1995) vốn bị giới hạn trong miền hình chữ nhật tựa đơn, hoặc phần tử rời rạc DKQ của Hwang & Park (2000) bị mất tính liên tục đạo hàm bậc cao, luận án đã sử dụng IGA trên cơ sở trích xuất Bézier. Phương pháp này biến đổi hàm cơ sở NURBS trơn $C^{p-1}$ toàn cục thành các hàm đa thức Bernstein cục bộ trên từng phần tử, cho phép phân tích chính xác tuyệt đối các biên cong và miền liên kết đa mảng (multi-patch) với chi phí bậc tự do thấp hơn đáng kể.

3. Phát hiện số học nào mang tính bất ngờ nhất từ dữ liệu?

Phát hiện về sự tương hỗ triệt tiêu giữa độ xốp và hạt nano GPLs: Mặc dù độ xốp bên trong vật liệu ($e_0 = 0.4$) làm suy giảm nghiêm trọng độ cứng uốn tĩnh, việc bổ sung chỉ đúng $1\text{ wt}%$ hạt nano graphene phân tán theo Pattern A không chỉ bù đắp hoàn toàn độ cứng bị mất mà còn nâng cao tần số dao động riêng cơ bản lên cao hơn cả tấm đặc nguyên bản cùng kích thước.

4. Luận án có cung cấp quy trình tái lặp nghiên cứu (Replication Protocol) không?

Có. Toàn bộ thông số tọa độ điểm kiểm soát, vector nút, ma trận vật liệu (PZT-4, PZT-5H, nhôm, gốm $\text{Al}_2\text{O}_3$), quy tắc cầu nghiệm Gauss và thuật toán trích xuất Bézier đều được trình bày chi tiết thông qua các bảng thông số số học và lưu đồ thuật toán mẫu trong các chương 2, 4, 5 và 6.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo là gì?

Mở rộng khung tính toán IGA Bézier sang bài toán tối ưu hóa cấu trúc hình học tô-pô (Topology Optimization) cho cấu trúc mạng xốp 3D in vi mô, phân tích độ bền mỏi đa trục cơ - điện - từ kết hợp hiệu ứng nhiệt phi tuyến, và phát triển hệ thống điều khiển giảm chấn thông minh thích ứng theo thời gian thực ứng dụng mạng nơ-ron sâu (Deep Learning Control).


Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thị Bích Liễu là một công trình khoa học xuất sắc, mẫu mực và hoàn chỉnh trong chuyên ngành Cơ học tính toán và Kỹ thuật công trình, tạo lập 6 đóng góp cốt lõi:

  1. Xây dựng hoàn chỉnh khung phân tích Isogeometric Analysis (IGA) dựa trên toán tử trích xuất Bézier cho các kết cấu tấm composite thông minh và FGM xốp.
  2. Phát triển lý thuyết biến dạng cắt bậc cao không ràng buộc tổng quát (UHSDT) loại bỏ hệ số hiệu chỉnh cắt SCFs.
  3. Làm sáng tỏ cơ chế gia cường vi cấu trúc của hạt nano GPLs kết hợp phân bố độ rỗng đối với ứng xử uốn phi tuyến và dao động của tấm PFGP-GPLs.
  4. Thiết lập thuật toán điều khiển chủ động vòng kín với độ lợi $G_d, G_v$, triệt tiêu hiệu quả dao động tĩnh và quá độ dưới tải trọng nổ và tải xung lực.
  5. Giải quyết trọn vẹn bài toán cơ - điện ghép nối của tấm FGPMP trên các miền hình học phức tạp, tấm tròn, vành khăn và tấm có lỗ khoét tùy ý.
  6. Mở ra con đường thực tiễn để tích hợp trực tiếp công nghệ IGA tiên tiến vào các hệ thống phần mềm CAE thương mại, đặt nền móng vững chắc cho kỷ nguyên thiết kế và phân tích cơ học hợp nhất trong tương lai.