Tổng quan về luận án

Nghiên cứu của Tiến sĩ Vũ Văn Thẩm với đề tài "Phân tích tĩnh và động tấm, vỏ thoải hai độ cong composite nano carbon – áp điện" (Chuyên ngành: Cơ kỹ thuật, Mã số: 9520101, Trường Đại học Xây dựng Hà Nội, 2021; Người hướng dẫn: PGS. TS. Trần Hữu Quốc và GS. TS. Trần Minh Tú) đại diện cho một bước tiến tiên phong trong lĩnh vực cơ học tính toán vật liệu composite thông minh thế hệ mới tại Việt Nam và khu vực. Trong bối cảnh kỹ thuật hàng không vũ trụ, cơ điện tử chính xác và công trình phòng thủ hiện đại đòi hỏi các kết cấu siêu nhẹ, độ cứng cao cùng khả năng tự cảm biến và điều khiển chủ động dao động, vật liệu composite có cơ tính biến thiên được gia cường bởi ống nano carbon tích hợp các lớp áp điện (Piezoelectric Functionally Graded Carbon Nanotube-Reinforced Composite – PFG-CNTRC) đang nổi lên như một giải pháp vật liệu đột phá.

                    ┌──────────────────────────────────────────────────────────┐
                    │               Lớp Áp Điện Trên (Actuator)                │
                    ├──────────────────────────────────────────────────────────┤
                    │    Lõi FG-CNTRC (Phân bố CNT: UD, FG-V, FG-Λ, FG-X, FG-O)│
                    │      (Polyme đẳng hướng PMMA + Cốt sợi Nano Carbon)      │
                    ├──────────────────────────────────────────────────────────┤
                    │               Lớp Áp Điện Dưới (Sensor)                  │
                    └──────────────────────────────────────────────────────────┘

Khoảng trống học thuật (Research Gap) cốt lõi được luận án chỉ ra: Mặc dù các công trình quốc tế như Zhu et al. (2012) hay Shen & Zhang (2012) đã phân tích ứng xử của tấm FG-CNTRC phẳng đơn thuần, các công bố về kết cấu vỏ thoải hai độ cong (doubly-curved shallow shells) PFG-CNTRC đa trường cơ - điện vẫn còn vô cùng khiêm tốn. Về mặt lý thuyết mô hình hóa:

  1. Lý thuyết đàn hồi 3 chiều (3D Elasticity) tuy chính xác tuyệt đối nhưng biểu thức toán học cồng kềnh, bất khả thi khi giải các bài toán vỏ nhiều lớp có điều kiện biên phức tạp.
  2. Lý thuyết biến dạng cắt bậc nhất (FSDT) bắt buộc phải sử dụng hệ số hiệu chỉnh cắt ($k = 5/6$) vốn phụ thuộc phi tuyến vào cấu trúc lớp, hình học và điều kiện biên.
  3. Lý thuyết biến dạng cắt bậc ba (TSDT) của Reddy (1984) sử dụng 5 ẩn số chuyển vị, khiến hệ phương trình cân bằng và ma trận phần tử hữu hạn mở rộng cồng kềnh.

Để giải quyết triệt để các hạn chế trên, luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và 2 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để cải tiến lý thuyết biến dạng cắt bậc cao 4 ẩn chuyển vị (HSDST-4) cho kết cấu vỏ thoải hai độ cong đa lớp PFG-CNTRC nhằm triệt tiêu hoàn toàn ứng suất cắt ngang tại hai mặt tự do mà không cần sử dụng hệ số hiệu chỉnh cắt?
  • RQ2: Mức độ ảnh hưởng định lượng của các quy luật phân bố ống nano carbon (UD, FG-V, FG-$\Lambda$, FG-X, FG-O) và tỷ phần thể tích $V^*_{CNT}$ đến độ võng, trường ứng suất và phổ tần số dao động riêng của vỏ là bao nhiêu?
  • RQ3: Cơ chế điều khiển phản hồi cơ - điện thông qua điện áp đặt $V_t, V_b$ và trạng thái mạch (mạch kín/mạch hở) kiểm soát độ cứng động học và đáp ứng quá độ của kết cấu dưới các dạng tải trọng va đập theo thời gian như thế nào?
  • H1: Mô hình phân bố cốt sợi nano carbon dạng chữ X (FG-X) tập trung CNT mật độ cao tại hai biên ngoài sẽ tối ưu hóa mô-men quán tính hiệu dụng, gia tăng độ cứng uốn và nâng cao tần số dao động riêng vượt trội so với phân bố đều (UD) và phân bố dạng tâm (FG-O).
  • H2: Việc ghép nối trường điện thế áp đặt thông qua lớp áp điện đóng vai trò như một bộ kích hoạt điều khiển chủ động, cho phép giảm chấn và triệt tiêu biến dạng uốn cơ học tĩnh lên tới trên 30%.

Khung lý thuyết của luận án được xây dựng dựa trên sự liên kết giữa Nguyên lý biến phân Hamilton, Hệ phương trình điện từ trường Maxwell, Lý thuyết vỏ thoải hai độ cong bậc cao 4 ẩn chuyển vị cải tiến (HSDST-4) và Quy tắc pha trộn mở rộng (Extended Rule of Mixtures - ERM) với các hệ số hiệu dụng CNT ($\eta_1, \eta_2, \eta_3$). Phạm vi nghiên cứu bao quát 4 dạng hình học kinh điển: Tấm phẳng (PLA, $R_x=R_y=\infty$), Vỏ trụ thoải (CYL, $R_x=R, R_y=\infty$), Vỏ cầu thoải (SPH, $R_x=R_y$), và Vỏ paraboloid hyperbol/yên ngựa (HPR, $R_x=-R_y$) với dải thông số tỷ lệ hình học $a/h = 10 \div 50$, tỷ số độ cong $R_x/a = 5 \div 20$, tỷ lệ thể tích CNT $V^*_{CNT} = 0.11 \div 0.28$.

Literature Review và Positioning

Tổng quan tài liệu của luận án hệ thống hóa 3 dòng nghiên cứu (research streams) chính yếu trong cơ học vật liệu và kết cấu tiên tiến:

                            ┌─────────────────────────────────────────┐
                            │    VẬT LIỆU COMPOSITE NANO VÀ ÁP ĐIỆN   │
                            └────────────────────┬────────────────────┘
                                                 │
         ┌───────────────────────────────────────┼───────────────────────────────────────┐
         ▼                                       ▼                                       ▼
┌─────────────────────────────────┐   ┌─────────────────────────────────┐   ┌─────────────────────────────────┐
│     Dòng 1: Vật liệu FG-CNTRC   │   │  Dòng 2: Kết cấu Thông minh     │   │   Dòng 3: Mô hình Cơ học Tấm Vỏ │
├─────────────────────────────────┤   ├─────────────────────────────────┤   ├─────────────────────────────────┤
│• Shen (2009): Khái niệm FG-CNTRC│   │• Curie (1880), Voigt (1894)     │   │• Kirchhoff-Love (CST): Mỏng     │
│• Zhu et al. (2012): FSDT tấm    │   │• Mindlin (1972): Sóng áp điện   │   │• Reissner-Mindlin (FSDT): k=5/6 │
│• Alibeigloo (2014): Đàn hồi 3D  │   │• Tzou & Gadre (1989): Shell FEM │   │• Reddy (1984): TSDT 5 ẩn        │
│• Ke et al. (2012): Timoshenko   │   │• Kiani (2016): Tấm PFG-CNTRC    │   │• Shimpi & Patel (2006): HSDT-4  │
└─────────────────────────────────┘   └─────────────────────────────────┘   └─────────────────────────────────┘
  1. Dòng nghiên cứu vật liệu FG-CNTRC: Khởi xướng bởi Shen (2009), ý tưởng phân bố gradient nồng độ CNT dọc theo chiều dày ma trận polyme đẳng hướng đã mở ra hướng đi đột phá nhằm khắc phục nhược điểm kết tụ hạt nano. Zhu et al. (2012) sử dụng FSDT khảo sát uốn và dao động tự do của tấm FG-CNTRC. Alibeigloo & Liew (2013, 2014) thiết lập lời giải giải tích 3D cho tấm và vỏ trụ FG-CNTRC chịu tải cơ - nhiệt. Ke et al. (2012) và Lei et al. (2015, 2016) khảo sát ổn định và dao động phi tuyến của dầm/vỏ trụ bằng phương pháp không lưới kp-Ritz.
  2. Dòng nghiên cứu kết cấu thông minh tích hợp áp điện: Bắt nguồn từ khám phá của Jacques & Pierre Curie (1880) và lý thuyết áp điện tuyến tính của Voigt (1894), Mindlin (1972) đã đặt nền móng cho cơ học tấm áp điện. Tiếp đó, Tzou & Gadre (1989), Lee (1990), Kioua & Mirza (2000) phát triển mô hình phần tử hữu hạn cho vỏ mỏng áp điện. Gần đây, Kiani (2016) và Selim et al. (2016, 2017) đã mở rộng bài toán sang tấm FG-CNTRC gắn lớp áp điện sử dụng phương pháp Ritz và IMLS-Ritz.
  3. Dòng phát triển lý thuyết biến dạng cắt tấm và vỏ: Chuyển biến từ lý thuyết cổ điển Kirchhoff-Love (CST) sang lý thuyết biến dạng cắt bậc nhất (FSDT) của Reissner-Mindlin, và lý thuyết biến dạng cắt bậc ba (TSDT) của Reddy (1984). Đặc biệt, Shimpi & Patel (2006) đề xuất lý thuyết 4 ẩn chuyển vị cho tấm đẳng hướng, sau đó được mở rộng bởi Thai & Kim (2013), Daouadji et al. (2012).

Trong bức tranh tổng quan đó, các mâu thuẫn học thuật (academic debates) nổi bật bao gồm:

  • Mâu thuẫn giữa tính chính xác 3D và hiệu quả tính toán 2D: Lời giải đàn hồi 3 chiều của Alibeigloo (2014) mang tính chuẩn mực nhưng giới hạn nghiêm ngặt ở các dạng hình học đối xứng và biên tựa đơn giản; ngược lại, các mô hình 2D tương đương (ESL) dễ triển khai nhưng thường hy sinh độ chính xác phân bố ứng suất cắt qua chiều dày.
  • Mâu thuẫn về số lượng biến trường trong mô hình cắt bậc cao: Các mô hình HSDT truyền thống của Reddy (1984) hay Loja et al. (2001) duy trì 5 biến độc lập ($u_0, v_0, w, \theta_x, \theta_y$), dẫn đến ma trận độ cứng lớn và chi phí tính toán cao khi giải bài toán vỏ phi tuyến hoặc động lực học phức tạp.

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

  • So với Zhu et al. (2012): Công trình của Zhu chỉ áp dụng FSDT cho tấm phẳng và phụ thuộc vào hệ số cắt $k=5/6$. Luận án của Vũ Văn Thẩm mở rộng cho vỏ thoải hai độ cong đa dạng (SPH, CYL, HPR) bằng HSDST-4, tự động thỏa mãn ứng suất tiếp triệt tiêu tại biên mà không cần hệ số hiệu chỉnh.
  • So với Kiani (2016) và Selim et al. (2016): Kiani sử dụng phương pháp Ritz với đa thức Chebyshev chỉ cho tấm PFG-CNTRC phẳng; nghiên cứu của tác giả thiết lập trọn vẹn cả lời giải giải tích Navier lẫn mô hình phần tử hữu hạn 4 nút 32 bậc tự do cho hệ vỏ cong hai chiều chịu tải trọng va đập phức tạp.

Đó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 lý thuyết biến dạng cắt bậc cao bốn ẩn chuyển vị (HSDST-4) từ kết cấu tấm phẳng sang hình học vỏ thoải hai độ cong có tương tác cơ - điện. Trường chuyển vị được thiết lập thông qua việc phân tách độ võng toàn phần $w(x,y,z,t)$ thành thành phần uốn $w_b$ và thành phần cắt $w_s$:

$$w(x, y, z, t) = w_b(x, y, t) + w_s(x, y, t)$$

Trường chuyển vị tổng quát $(u, v, w)$ tại điểm bất kỳ có tọa độ $(x,y,z)$ được định nghĩa chuẩn xác:

$$u(x, y, z, t) = u_0(x, y, t) - z\left(\frac{u_0}{R_x} + \frac{\partial w_b}{\partial x}\right) - f(z)\frac{\partial w_s}{\partial x}$$

$$v(x, y, z, t) = v_0(x, y, t) - z\left(\frac{v_0}{R_y} + \frac{\partial w_b}{\partial y}\right) - f(z)\frac{\partial w_s}{\partial y}$$

Hàm phân bố ứng suất cắt $f(z)$ được luận án cải tiến chọn lọc dưới dạng:

$$f(z) = z\left[-\frac{1}{8} + \frac{3}{2}\left(\frac{z}{h_{tam}}\right)^2\right]$$

với $g(z) = 1 - f'(z) = 1 - \left[-\frac{1}{8} + \frac{9}{2}\left(\frac{z}{h_{tam}}\right)^2\right]$, bảo đảm điều kiện biên ứng suất cắt $\tau_{xz}, \tau_{yz}$ triệt tiêu tuyệt đối tại bề mặt trên và dưới của kết cấu ($z = \pm h_{tam}/2$).

       Mặt trên (z = +h/2) ──────────────────────────  τ_xz = 0, τ_yz = 0
                                 ▲
                                 │ Phân bố parabol bậc 2 qua chiều dày:
                                 │ g(z) = 1 - f'(z)
                                 ▼
       Mặt dưới (z = -h/2) ──────────────────────────  τ_xz = 0, τ_yz = 0

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

Khung phân tích của luận án tích hợp chặt chẽ 3 lý thuyết nền tảng:

  1. Lý thuyết vật liệu micro-nano mechanics (Shen, 2009): Xác định các mô-đun đàn hồi hiệu dụng $E_{11}(z), E_{22}(z), G_{12}(z), \nu_{12}$ của từng lớp nano composite theo quy tắc pha trộn mở rộng với các thông số hiệu quả $\eta_1, \eta_2, \eta_3$.
  2. Lý thuyết vỏ thoải hai độ cong Donnell-Mushtari-Vlasov: Kết hợp các thành phần độ cong $1/R_x, 1/R_y$ vào ten-xơ biến dạng màng và mô-men uốn.
  3. Lý thuyết điện từ trường áp điện (Piezoelectricity Theory): Mô hình hóa cảm ứng điện dịch $D_i$ và điện trường $E_i = -\partial \Phi / \partial x_i$, tích hợp phương trình Maxwell $\text{div} \mathbf{D} = 0$.

Điều kiện biên và giới hạn ứng dụng của mô hình: Áp dụng chuẩn xác cho các vỏ thoải có độ cong nông ($a/R_x \le 0.2, b/R_y \le 0.2$), giả định liên kết giữa các lớp composite và lớp áp điện là dính bám lý tưởng tuyệt đối (perfect bonding), bỏ qua biến dạng pháp tuyến qua chiều dày ($\epsilon_{zz} = 0$).

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học thực chứng (Positivism) và hệ hình diễn dịch tính toán cơ học (Computational Deductive Paradigm). Bản chất của nghiên cứu là tiếp cận mô hình hóa liên tục (Continuum Mechanics), thiết lập các phương trình vi phân đạo hàm riêng chi phối và giải quyết thông qua hai nhánh phương pháp song song: Phương pháp giải tích chính xác (Analytical Navier Method) và Phương pháp số (Finite Element Method - FEM).

                        ┌────────────────────────────────────────────────────────┐
                        │                 NGUYÊN LÝ BIẾN PHÂN HAMILTON           │
                        │                    δ ∫ (T - U + W) dt = 0              │
                        └───────────────────────────┬────────────────────────────┘
                                                    │
                 ┌──────────────────────────────────┴──────────────────────────────────┐
                 ▼                                                                     ▼
┌──────────────────────────────────────────────┐     ┌─────────────────────────────────────────────────┐
│     TIẾP CẬN GIẢI TÍCH (Chương 2)            │     │       TIẾP CẬN PHẦN TỬ HỮU HẠN (Chương 3)       │
├──────────────────────────────────────────────┤     ├─────────────────────────────────────────────────┤
│• Chuỗi kép Fourier (Nghiệm Navier)           │     │• Phần tử chữ nhật 4 nút, 8 DOF/nút (C1)         │
│• Điều kiện biên tựa khớp 4 cạnh (SSSS)       │     │• Đa dạng điều kiện biên (CCCC, SCSC, CFFF...)   │
│• Giải bài toán trị riêng tĩnh và dao động    │     │• Tích hợp tải va đập và mạch điều khiển hồi tiếp│
└──────────────────────────────────────────────┘     └─────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

  1. Thiết lập phương trình chủ đạo: Từ nguyên lý biến phân Hamilton $\int_{t_1}^{t_2} (\delta T - \delta U + \delta W) dt = 0$, hệ 6 phương trình cân bằng liên hợp cơ - điện phi tuyến được dẫn xuất đầy đủ.
  2. Lời giải giải tích Navier (Chương 2): Áp dụng cho bài toán vỏ 4 cạnh tựa khớp (SSSS). Các hàm chuyển vị và điện thế được khai triển chuỗi lượng giác: $$u_0 = U_{mn} \cos\left(\frac{m\pi x}{a}\right)\sin\left(\frac{n\pi y}{b}\right)e^{i\omega t}$$ $$v_0 = V_{mn} \sin\left(\frac{m\pi x}{a}\right)\cos\left(\frac{n\pi y}{b}\right)e^{i\omega t}$$ $$w_b = W_{bmn} \sin\left(\frac{m\pi x}{a}\right)\sin\left(\frac{n\pi y}{b}\right)e^{i\omega t}$$ $$w_s = W_{smn} \sin\left(\frac{m\pi x}{a}\right)\sin\left(\frac{n\pi y}{b}\right)e^{i\omega t}$$ Hệ phương trình quy về bài toán đại số tuyến tính $[S]{\Delta} - \omega^2 [M]{\Delta} = {0}$.
  3. Xây dựng mô hình Phần tử hữu hạn 4 nút (Chương 3):
    • Sử dụng phần tử vỏ chữ nhật 4 nút với 8 bậc tự do tại mỗi nút: $$\mathbf{d}i = \left[ u{0i}, v_{0i}, w_{bi}, \frac{\partial w_{bi}}{\partial x}, \frac{\partial w_{bi}}{\partial y}, w_{si}, \frac{\partial w_{si}}{\partial x}, \frac{\partial w_{si}}{\partial y} \right]^T$$
    • Tổng cộng 32 bậc tự do cho một phần tử. Các trường $u_0, v_0$ được nội suy bằng hàm dạng Lagrange bậc nhất ($C^0$), trong khi các trường $w_b, w_s$ được nội suy bằng hàm Hermite bậc ba ($C^1$), thỏa mãn sự liên tục của độ võng và độ dốc tại biên phần tử.
Node 4 (0,b) ─────────────────────────── Node 3 (a,b)
     │                                        │
     │      Phần tử chữ nhật 4 nút (C1)       │
     │        8 bậc tự do / 1 nút             │
     │        Tổng: 32 DOFs / Element         │
     │                                        │
Node 1 (0,0) ─────────────────────────── Node 2 (a,0)
  1. Tích phân số và kỹ thuật xử lý: Tích phân Gauss-Legendre $2\times 2$ được sử dụng để tính toán các ma trận độ cứng cơ học $[K_{uu}]$, ma trận tương tác cơ - điện $[K_{u\phi}]$, ma trận điện dung $[K_{\phi\phi}]$ và ma trận khối lượng $[M]$. Khử biến điện thế ở cấp độ phần tử thông qua kỹ thuật ngưng tụ tĩnh: $$[K_{piezo}] = [K_{u\phi}][K_{\phi\phi}]^{-1}[K_{\phi u}]$$

Data và phân tích

Toàn bộ thuật toán giải tích và phần tử hữu hạn được lập trình trực tiếp trên môi trường MATLAB.

  • Thông số vật liệu kiểm chứng:
    • Pha nền Polyme (PMMA): $E_m = 2.1\text{ GPa}$, $\nu_m = 0.34$, $\rho_m = 1150\text{ kg/m}^3$.
    • Ống nano carbon SWCNT $(10,10)$: $E_{11}^{CNT} = 5.6466\text{ TPa}$, $E_{22}^{CNT} = 7.08\text{ TPa}$, $G_{12}^{CNT} = 1.9445\text{ TPa}$, $\nu_{12}^{CNT} = 0.175$, $\rho_{CNT} = 1400\text{ kg/m}^3$.
    • Các hệ số hiệu dụng $\eta_1, \eta_2, \eta_3$ ứng với các tỷ phần thể tích $V^*_{CNT}$:
      • $V^*_{CNT} = 0.11 \implies \eta_1 = 0.149, \eta_2 = 0.934, \eta_3 = 0.934$
      • $V^*_{CNT} = 0.14 \implies \eta_1 = 0.150, \eta_2 = 0.941, \eta_3 = 0.941$
      • $V^*_{CNT} = 0.17 \implies \eta_1 = 0.149, \eta_2 = 1.381, \eta_3 = 1.381$
      • $V^*_{CNT} = 0.28 \implies \eta_1 = 0.285, \eta_2 = 1.109, \eta_3 = 1.109$
    • Lớp áp điện PZT-4 / PZT-5A: $E_p = 63\text{ GPa}$, $\nu_p = 0.3$, $\rho_p = 7600\text{ kg/m}^3$, $e_{31} = e_{32} = -4.1\text{ C/m}^2$, $p_{11} = p_{22} = 1.53\times 10^{-8}\text{ F/m}$, $p_{33} = 1.50\times 10^{-8}\text{ F/m}$.
  • Kiểm tra độ hội tụ: Lưới phần tử hữu hạn $20\times 20$ (441 nút, 3528 ẩn số) đạt độ hội tụ đơn điệu với sai số tương đối so với lời giải giải tích Navier dưới $0.12%$.

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

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

        Tần số dao động riêng (Hz) theo kiểu phân bố CNT (V*CNT = 28%)
        ┌─────────────────────────────────────────────────────────────┐
   FG-X │ ████████████████████████████████████████████ 100% (Cao nhất)│
     UD │ ██████████████████████████████ 78%                          │
   FG-V │ ██████████████████████████ 68%                              │
   FG-Λ │ ██████████████████████████ 68%                              │
   FG-O │ ██████████████████ 52% (Thấp nhất)                          │
        └─────────────────────────────────────────────────────────────┘
  1. Quy luật phân bố FG-X sở hữu hiệu năng cơ học vượt trội: Với cùng một tỷ phần thể tích $V^*{CNT} = 0.28$, kết cấu có kiểu phân bố CNT dạng FG-X làm tăng tần số dao động riêng thứ nhất lên 38.5% – 48.2% so với kiểu FG-O, và cao hơn 16.4% – 22.1% so với phân bố đều UD. Đồng thời, độ võng uốn tĩnh chính giữa kết cấu của kiểu FG-X giảm hơn 42% so với FG-O. Nguyên nhân lý thuyết: FG-X tập trung mật độ ống nano carbon cao nhất tại hai bề mặt ngoài, nơi có ứng suất uốn $\sigma{xx}, \sigma_{yy}$ cực đại, tối ưu hóa mô-men kháng uốn tổng thể.
  2. Hiệu ứng phi tuyến của tỷ lệ thể tích $V^*_{CNT}$: Khi tăng $V^*_{CNT}$ từ $11%$ lên $28%$, độ cứng uốn tăng vọt khiến độ võng không thứ nguyên $\bar{w}$ của vỏ cầu thoải giảm từ $0.0482$ xuống $0.0195$ (giảm hơn 59.5%), trong khi tần số dao động riêng tăng tương ứng trên 62%.
  3. Ảnh hưởng quyết định của dạng hình học và độ cong vỏ: Khảo sát 4 dạng kết cấu cho thấy độ cứng giảm dần theo thứ tự: Vỏ cầu thoải (SPH) > Vỏ trụ thoải (CYL) > Tấm phẳng (PLA) > Vỏ yên ngựa (HPR). Khi tỷ số bán kính cong $R_x/a$ giảm từ 20 xuống 5 (vỏ có độ cong lớn hơn), tần số riêng của vỏ cầu tăng 28.4% do xuất hiện hiệu ứng màng (membrane effect) chuyển hóa năng lượng uốn thành năng lượng kéo/nén trong mặt phẳng.
  4. Cơ chế điều khiển chủ động bằng điện thế áp đặt ($V_t, V_b$):
    • Khi đặt điện áp kích thích $V_t > 0$ lên lớp áp điện trên, độ võng uốn của vỏ giảm tuyến tính. Tại mức điện thế $V_t = 200\text{V} \div 400\text{V}$, biến dạng uốn do tải trọng cơ học gây ra có thể được triệt tiêu hoàn toàn về mức 0.
    • Ứng suất kéo tiếp tuyến tại mặt đáy giảm 28% – 35%, chuyển trạng thái chịu lực của kết cấu sang nén đồng đều an toàn.
  5. Hiệu ứng độ cứng áp điện cảm ứng (Piezoelectric Stiffening Effect): Tần số dao động riêng của vỏ PFG-CNTRC ở trạng thái mạch hở (open circuit) luôn cao hơn trạng thái mạch kín (closed circuit) từ 2.5% đến 6.8%. Bản chất là do mạch hở không cho phép giải phóng điện tích, sinh ra trường phản điện trường làm tăng ma trận độ cứng tương đương của kết cấu.

Implications đa chiều

  • Về mặt lý thuyết & phương pháp: Khẳng định tính đúng đắn và ưu việt của lý thuyết HSDST-4 cải tiến; chứng minh việc giảm từ 5 biến xuống 4 biến chuyển vị hoàn toàn không làm suy giảm độ chính xác mà còn tiết kiệm 20% bộ nhớ tính toán và giảm thời gian giải ma trận trên máy tính.
  • Về ứng dụng kỹ thuật thực tiễn: Cung cấp công thức và biểu đồ thiết kế tối ưu cho kỹ sư kết cấu hàng không vũ trụ (vỏ thân máy bay phản lực, cánh lái tên lửa, vỏ vệ tinh) và kỹ thuật biển, nơi các tấm/vỏ composite chịu tải trọng va đập và yêu cầu triệt tiêu rung động tức thời thông qua mạch điều khiển hồi tiếp $G_d, G_v$.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu (limitations):

  1. Giả thiết tuyến tính hình học: Luận án giới hạn ở bài toán biến dạng bé và góc xoay nhỏ; chưa xét đến phi tuyến hình học Von Kármán trong các điều kiện tải trọng động lực học cực hạn gây mất ổn định lớn.
  2. Giả thiết liên kết lớp lý tưởng: Chưa xét đến hiện tượng trượt biên, nứt tế vi nền polymer hoặc hiện tượng tách lớp (delamination) giữa lớp áp điện và lớp lõi composite khi chịu tải chu kỳ dài hạn.
  3. Phạm vi nhiệt - ẩm: Chưa tích hợp ảnh hưởng đồng thời của trường nhiệt độ biến thiên phi tuyến và độ ẩm môi trường (Hygo-thermal effects) đến sự suy giảm cơ tính của mạng liên kết polymer-CNT.
  4. Giới hạn vỏ thoải: Mô hình chỉ áp dụng cho vỏ thoải ($R_x, R_y \ge 5a$); chưa mở rộng cho các kết cấu vỏ sâu (deep shells) hoặc vỏ có độ dày lớn ($a/h < 10$) đòi hỏi lý thuyết 3D phân lớp liên tục (Layer-wise Theory).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Phát triển mô hình HSDST-4 phi tuyến hình học đầy đủ kết hợp thuật toán lặp Newmark/Newton-Raphson khảo sát hiện tượng mất ổn định động lực học (dynamic buckling).
  • Hướng 2: Áp dụng phương pháp phần tử hữu hạn mở rộng (XFEM) hoặc mô hình trường pha (Phase-field model) để mô phỏng sự lan truyền vết nứt và tách lớp cơ - điện.
  • Hướng 3: Tối ưu hóa phân bố gradient ống nano carbon hai chiều và ba chiều (2D/3D FG-CNTRC) kết hợp công nghệ in 3D vật liệu composite chức năng.

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

  • Tác động học thuật: Các kết quả của luận án đã được công bố trên các tạp chí khoa học chuyên ngành Cơ học và Xây dựng uy tín. Dữ liệu số trong luận án cung cấp hệ thống nghiệm chuẩn (Benchmark solutions) có độ tin cậy cao để các nhà nghiên cứu trong nước và quốc tế đối chiếu, kiểm chứng các thuật toán số và mô hình phần tử mới.
  • Tác động công nghiệp và chuyển giao công nghệ: Kết quả nghiên cứu tạo tiền đề khoa học vững chắc cho các viện thiết kế hàng không vũ trụ, công nghiệp quốc phòng và giao thông vận tải trong việc lựa chọn cấu hình vật liệu thông minh (tỷ lệ $V^*_{CNT}$, kiểu FG-X, thông số mạch áp điện) nhằm tối ưu hóa tải trọng bản thân và tăng độ bền mỏi công trình.

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

                        ┌────────────────────────────────────────────────────────┐
                        │                 ĐỐI TƯỢNG HƯỞNG LỢI                    │
                        └───────────────────────────┬────────────────────────────┘
                                                    │
         ┌──────────────────────────────────────────┼──────────────────────────────────────────┐
         ▼                                          ▼                                          ▼
┌─────────────────────────────────┐        ┌─────────────────────────────────┐        ┌─────────────────────────────────┐
│       Nghiên cứu sinh &         │        │       Kỹ sư R&D Công nghiệp     │        │     Chuyên gia Cơ học           │
│       Học viên Cao học          │        │       Hàng không - Vũ trụ       │        │     và Cơ sở Tiêu chuẩn         │
├─────────────────────────────────┤        ├─────────────────────────────────┤        ├─────────────────────────────────┤
│• Tài liệu tham khảo chuẩn mực   │        │• Hướng dẫn lựa chọn phân bố CNT │        │• Nền tảng xây dựng quy chuẩn    │
│• Mã nguồn MATLAB giải tích/FEM  │        │• Sơ đồ điều khiển áp điện       │        │  tính toán kết cấu thông minh   │
│• Phương pháp luận HSDST-4       │        │• Tối ưu hóa khối lượng kết cấu  │        │• Dữ liệu kiểm chứng chuẩn       │
└─────────────────────────────────┘        └─────────────────────────────────┘        └─────────────────────────────────┘
  1. Nghiên cứu sinh và Giảng viên ngành Cơ học/Xây dựng/Hàng hải: Tiếp cận phương pháp luận mô hình hóa kết cấu vỏ cong 4 ẩn chuyển vị cải tiến và toàn bộ thuật toán lập trình MATLAB mẫu được cung cấp chi tiết trong phụ lục luận án.
  2. Kỹ sư R&D vật liệu tiên tiến: Khai thác dữ liệu định lượng về hiệu năng của 5 kiểu phân bố CNT để thiết kế các chi tiết vỏ composite thông minh chịu tải rung động khắc nghiệt.
  3. Các cơ quan quản lý và Viện nghiên cứu Tiêu chuẩn: Làm cơ sở khoa học để từng bước biên soạn hệ thống tiêu chuẩn kỹ thuật thiết kế kết cấu composite nano và kết cấu tích hợp vật liệu áp điện tại Việt Nam.

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ì?

Đóng góp độc đáo nhất là việc phát triển thành công Lý thuyết vỏ thoải hai độ cong bậc cao 4 ẩn chuyển vị (HSDST-4) cải tiến cho vật liệu PFG-CNTRC. Bằng cách phân tách thông minh chuyển vị ngang $w = w_b + w_s$ và thiết lập hàm phân bố $f(z) = z[-1/8 + 3/2(z/h)^2]$, mô hình vừa đảm bảo ứng suất tiếp triệt tiêu tại hai mặt biên tự do, vừa loại bỏ hoàn toàn hệ số hiệu chỉnh cắt $k_s$, đồng thời giảm số ẩn chuyển vị từ 5 (như TSDT của Reddy) xuống còn 4, tối ưu hóa ma trận độ cứng tổng thể.

2. Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với mô hình FSDT của Zhu et al. (2012) vốn phải giả định hệ số cắt $k=5/6$ bất kể góc xếp lớp, và mô hình TSDT 5 ẩn của Reddy (1984) hoặc Selim et al. (2016), phương pháp phần tử hữu hạn 4 nút $C^1$ (32 DOFs) của luận án tích hợp trực tiếp chuyển dịch điện tích áp điện và kỹ thuật ngưng tụ ma trận tĩnh, cho phép phân tích đồng thời 4 dạng vỏ cong (SPH, CYL, HPR, PLA) với độ chính xác tương đương nghiệm đàn hồi 3D nhưng thời gian tính toán giảm hơn 20%.

3. Phát hiện số nào gây bất ngờ và có ý nghĩa kỹ thuật cao nhất?

Phát hiện về khả năng triệt tiêu ứng suất kéo chính và đảo chiều biến dạng uốn thông qua điện thế áp đặt $V_t$. Khi điện thế $V_t$ tăng từ 0V lên 200V–400V, nó sinh ra một mô-men uốn ngược chiều hoàn toàn bù trừ mô-men do tải trọng cơ học phân bố gây ra, đưa độ võng đỉnh vỏ cầu về 0 và làm giảm đỉnh ứng suất kéo nguy hiểm tới trên 30%, biến kết cấu thành một hệ thống tự thích ứng (Self-adapting smart structure).

4. Luận án có cung cấp đầy đủ giao thức để tái lập nghiên cứu (Replication Protocol)?

Có. Luận án trình bày chi tiết toàn bộ các ma trận cơ sở $[A], [B], [D], [B_s], [D_s], [H_s], [A_s]$, các công thức quy đổi ten-xơ áp điện, các hàm nội suy Hermite/Lagrange, và cung cấp toàn văn 5 chương trình tính toán hoàn chỉnh viết bằng ngôn ngữ MATLAB trong phần Phụ lục (từ Phụ lục B1 đến B5: phân tích tĩnh giải tích, dao động riêng giải tích, phân tích tĩnh PTHH, dao động riêng PTHH, và đáp ứng động lực học có điều khiển hồi tiếp).

5. Lộ trình phát triển học thuật 10 năm tiếp theo từ kết quả luận án?

Lộ trình 10 năm mở rộng bao gồm:

  • Giai đoạn 2021–2024: Tích hợp hiệu ứng phi tuyến hình học Von Kármán và tải trọng nhiệt - ẩm môi trường.
  • Giai đoạn 2025–2028: Xây dựng mô hình phá hủy cơ - điện, nứt và tách lớp liên hợp sử dụng phương pháp XFEM và cơ học vi mô phá hủy.
  • Giai đoạn 2029–2031: Tối ưu hóa cấu trúc topology đa mục tiêu cho vật liệu in 3D composite nano carbon tích hợp mạng lưới cảm biến áp điện thông minh không dây.

Kết luận

Luận án Tiến sĩ của tác giả Vũ Văn Thẩm đã 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. Hoàn thiện hệ thống lý thuyết biến dạng cắt bậc cao 4 ẩn chuyển vị cải tiến (HSDST-4) áp dụng cho vỏ thoải hai độ cong PFG-CNTRC.
  2. Dẫn xuất thành công hệ phương trình vi phân chuyển động liên hợp cơ - điện và thiết lập lời giải giải tích Navier chính xác cho vỏ 4 biên tựa khớp.
  3. Xây dựng phần tử hữu hạn chữ nhật 4 nút $C^1$ (32 bậc tự do) tin cậy, giải quyết linh hoạt mọi điều kiện biên phức tạp (CCCC, SCSC, CFFF).
  4. Khám phá định lượng quy luật vượt trội của phân bố FG-X, giúp tăng tần số dao động riêng lên đến 48.2% và giảm độ võng uốn 42% so với FG-O.
  5. Chứng minh cơ chế điều khiển chủ động của điện thế kích hoạt và trạng thái mạch điện trong việc triệt tiêu biến dạng và giảm chấn động lực học.
  6. Cung cấp bộ mã nguồn MATLAB chuẩn mực và hệ thống nghiệm số kiểm chứng (benchmark) có giá trị lâu dài cho cộng đồng khoa học cơ học tính toán.

Công trình không chỉ nâng cao vị thế học thuật của ngành Cơ học ứng dụng Việt Nam trên trường quốc tế mà còn mở ra những hướng ứng dụng đột phá trong thiết kế kết cấu công nghệ cao thế hệ mới.