Tổng quan về luận án

Nghiên cứu cơ học kết cấu hàng không vũ trụ hiện đại đòi hỏi sự phát triển của các hệ vật liệu tiên tiến có khả năng chịu tải trọng cơ học và môi trường nhiệt độ khắc nghiệt. Luận án tiến sĩ kỹ thuật "Phân tích tĩnh kết cấu vỏ trụ composite cơ tính biến thiên được gia cường bằng các ống nano carbon chịu tải trọng cơ và nhiệt độ" của tác giả Dương Văn Quang (chuyên ngành Cơ kỹ thuật, mã số: 9.01; người hướng dẫn: PGS.TS Trần Ngọc Đoàn, PGS.TS Đoàn Trắc Luật; Học viện Kỹ thuật Quân sự, 2023) đặt trọng tâm vào phân tích ứng xử nhiệt đàn hồi tĩnh của vỏ trụ nanocomposite tiên tiến. Bối cảnh khoa học của đề tài bắt nguồn từ sự xuất hiện của vật liệu composite cơ tính biến thiên gia cường bằng ống nano carbon (Functionally Graded Carbon Nanotube Reinforced Composites - FG-CNTRC). Văn bản luận án khẳng định rõ nét: "Vật liệu composite cơ tính biến thiên được gia cường bởi ống nano carbon (Functionally graded carbon nanotube reinforced composites FG-CNTRC) là loại vật liệu mới, có nhiều ưu điểm như: độ bền cao, khối lượng riêng nhỏ, tính chất nhiệt và tính chất điện tốt."

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định từ tổng quan y văn quốc tế bao gồm:

  1. Đa số các công bố trước đây (Shen [7], Zhu et al. [47], Ninh & Bich [51]) giả thiết trường nhiệt độ đồng nhất ($T = \text{const}$) hoặc gán trước phân bố nhiệt độ dạng tuyến tính/dạng sin, bỏ qua gradient nhiệt thực tế sinh ra từ phương trình truyền nhiệt qua chiều dày.
  2. Hầu hết các nghiên cứu cơ nhiệt bỏ qua sự phụ thuộc của các đặc tính cơ lý vào nhiệt độ (temperature-independent properties [14-18]), làm sai lệch đáng kể độ cứng của kết cấu ở môi trường nhiệt độ cao.
  3. Các giải pháp giải tích truyền thống chủ yếu giới hạn ở điều kiện biên gối tựa đơn giản (Navier-type solution cho biên S-S) [18, 80], chưa giải quyết được bài toán giải tích cho các liên kết biên phức tạp và phi đối xứng như ngàm (C-C), ngàm - tự do (C-F), ngàm - tựa (C-S) chịu tải trọng cục bộ bất kỳ.
  4. Việc áp dụng các lý thuyết vỏ cổ điển (CST) và lý thuyết biến dạng cắt bậc nhất (FSDT) đã bỏ qua thành phần ứng suất pháp tuyến ngang ($\sigma_z$), dẫn đến việc đánh giá thiếu chính xác trạng thái tập trung ứng suất tại vùng biên kết cấu dày.

Các câu hỏi nghiên cứu và giả thuyết khoa học của luận án:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác phân bố nhiệt độ thực và trường chuyển vị - biến dạng 3 chiều của vỏ trụ FG-CNTRC có xét đến tính chất vật liệu suy giảm theo nhiệt độ?
  • Câu hỏi nghiên cứu 2 (RQ2): Phương pháp giải tích nào có khả năng giải hệ phương trình vi phân đạo hàm riêng bậc cao mô tả ứng xử của vỏ trụ với các điều kiện biên thực tế (ngàm, tự do) chịu tải cục bộ?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật phân bố CNT (UD, FG-$\Lambda$, FG-V, FG-X, FG-O) và tỷ lệ thể tích $V_{CNT}$ ảnh hưởng định lượng như thế nào đến độ võng, ứng suất tiếp và ứng suất pháp ngang tại vùng hiệu ứng biên?
  • Giả thuyết khoa học (H1): Việc kết hợp lý thuyết biến dạng cắt bậc cao Quasi-3D với chuỗi lượng giác và phép biến đổi Laplace cho phép xác định chính xác nghiệm giải tích kín cho vỏ trụ FG-CNTRC với mọi điều kiện biên mà không cần phụ thuộc vào phương pháp số xấp xỉ lưới.
  • Giả thuyết khoa học (H2): Phân bố CNT tập trung tại các bề mặt chịu ứng suất cao (cấu hình FG-V, FG-X) sẽ triệt tiêu đáng kể độ võng và phân tán tập trung ứng suất biên so với cấu hình phân bố đều (UD).

Khung lý thuyết của nghiên cứu tích hợp: Lý thuyết biến dạng cắt bậc cao Quasi-3D (mô hình K3-2 phát triển từ Reddy HSDT), Quy luật hỗn hợp mở rộng (Extended Rule of Mixtures - EROM) với các hệ số hiệu dụng $\eta_i$ từ mô phỏng động lực học phân tử (MD), Lý thuyết truyền nhiệt dừng Fourier 1D và Nguyên lý công ảo (Principle of Virtual Work).

Phạm vi nghiên cứu bao gồm vỏ trụ tròn làm việc trong miền đàn hồi tuyến tính, vật liệu cốt gia cường là SWCNT dạng $(10,10)$ thẳng, đồng phương theo trục dọc, tỷ lệ thể tích $V_{CNT}$ từ $0.11$ đến $0.28$, tỷ số hình học $L/R \in [1, 6]$, $R/h \in [10, 50]$, dải nhiệt độ khảo sát từ $300\text{ K}$ đến $500\text{ K}$. Ý nghĩa khoa học và thực tiễn nằm ở việc cung cấp công cụ giải tích chuẩn mực và bộ số liệu đối chuẩn phục vụ thiết kế vỏ động cơ phản lực, thân tên lửa và thiết bị bay chịu tải trọng cơ - nhiệt liên hợp.

Literature Review và Positioning

Lịch sử nghiên cứu kết cấu FG-CNTRC khởi nguồn từ đề xuất mô hình vật liệu tiên phong của Shen (2009) [7], mở ra hướng đi thay thế composite truyền thống bằng vật liệu nano có tính chất cơ lý biến thiên liên tục theo chiều dày. Trải qua hơn một thập kỷ, các luồng nghiên cứu chính đã phân nhánh rõ rệt:

Thứ nhất, luồng nghiên cứu về lý thuyết tấm vỏ: Các công trình ban đầu của Ninh và Bich (2011) [51], Van Tung và Trang (2014) [85] sử dụng lý thuyết vỏ mỏng cổ điển Donnell (Donnell CST) kết hợp quan hệ phi tuyến von Kármán. Hạn chế cốt tử của CST là bỏ qua hoàn toàn biến dạng trượt ngang ($\gamma_{\xi z} = \gamma_{\theta z} = 0$), gây ra sai số nghiêm trọng khi phân tích vỏ có độ dày trung bình và dày ($R/h \le 20$). Để khắc phục, lý thuyết FSDT (Mindlin, 1951) được Zhu et al. (2012) [47], Liew et al. (2014) [93], Lei et al. (2015) [96] áp dụng rộng rãi. Tuy nhiên, FSDT đòi hỏi hệ số hiệu chỉnh cắt phụ thuộc vào hình học và điều kiện biên. Các lý thuyết biến dạng cắt bậc cao HSDT (Reddy, 1984) được Shen & Xiang (2012, 2014) [43, 50], Nguyễn Văn Thành et al. (2018-2020) [69, 70, 73] phát triển đã loại bỏ sự cần thiết của hệ số hiệu chỉnh cắt, nhưng đa số các dạng HSDT truyền thống vẫn mặc định chuyển vị pháp tuyến là hằng số ($w = w_0$), dẫn đến việc triệt tiêu vô lý ứng suất pháp tuyến ngang ($\sigma_z = 0$).

Thứ hai, luồng nghiên cứu về tải trọng nhiệt và sự phụ thuộc nhiệt độ: Đa số các nghiên cứu (Alibeigloo [14, 15, 18], Lê Thị Như Trang [52-61], Vũ Thanh Long [62-65], Phạm Thanh Hiếu [74-79]) giả định trường nhiệt độ phân bố đều hoặc các đặc tính cơ lý không phụ thuộc nhiệt độ. Trái lại, các thử nghiệm thực tế và mô phỏng động lực học phân tử (MD) của Han & Elliott (2007), Shen & Zhang (2012) [41] đã chứng minh mô-đun Young của SWCNT và nhựa nền polymer (như Epoxy) bị suy giảm từ $15%$ đến $40%$ khi nhiệt độ tăng từ $300\text{ K}$ lên $500\text{ K}$. Một số ít tác giả như Pourasghar & Chen (2016) [19], Moradi-Dastjerdi & Payganeh (2017) [82] đã đưa gradient nhiệt vào bài toán, nhưng hoàn toàn giới hạn ở điều kiện biên tựa đơn giản (Simply-supported).

Ma trận định vị khoảng trống nghiên cứu (Positioning Matrix):
┌─────────────────────────┬──────────────────┬─────────────────┬──────────────────────┬──────────────────────┐
│ Tác giả & Công trình    │ Lý thuyết vỏ     │ Gradient nhiệt  │ Tính chất cơ lý T(z) │ Điều kiện biên       │
├─────────────────────────┼──────────────────┼─────────────────┼──────────────────────┼──────────────────────┤
│ Shen & Xiang (2014)     │ HSDT (K3-0)      │ Không (Đồng nhất│ Có phụ thuộc nhiệt độ│ Biên tựa (S-S)       │
│ Alibeigloo (2014)       │ Đàn hồi 3D       │ Có gradient     │ Không phụ thuộc      │ Biên tựa (S-S)       │
│ Pourasghar & Chen (2016)│ FSDT             │ Có gradient     │ Có phụ thuộc nhiệt độ│ Biên tựa (S-S)       │
│ Vũ Văn Thẩm (2018)      │ HSDT4 Cải tiến   │ Không           │ Không phụ thuộc      │ Tựa (Navier) & FEM   │
│ Luận án Dương Văn Quang │ Quasi-3D (K3-2)  │ Có (Fourier 1D) │ Có phụ thuộc nhiệt độ│ Tùy ý (C-C, C-F, C-S)│
└─────────────────────────┴──────────────────┴─────────────────┴──────────────────────┴──────────────────────┘

Vị thế học thuật của luận án được khẳng định thông qua việc giải quyết đồng thời bài toán liên hợp: sử dụng lý thuyết Quasi-3D (K3-2) chứa ứng suất pháp tuyến ngang $\sigma_z$, tích hợp truyền nhiệt thực qua chiều dày với cơ tính phụ thuộc nhiệt độ $E(T,z), G(T,z), \alpha(T,z)$, và giải bằng giải tích kín cho các điều kiện biên ngàm phức tạp, vượt xa phạm vi của hai nghiên cứu quốc tế điển hình là Pourasghar & Chen (2016) [19] và Alibeigloo (2014) [18].

Đó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 khung lý thuyết biến dạng cắt bậc cao của J.N. Reddy sang dạng mô hình Quasi-3D (K3-2) áp dụng cho vật liệu dị tính nano trực hướng trong không gian tọa độ cong trực giao $O\xi\theta z$. Thay vì giả định $w(\xi, \theta, z) = w_0(\xi, \theta)$, luận án thiết lập trường chuyển vị 3 thành phần mở rộng: $$u(\xi, \theta, z) = u_0(\xi, \theta) + u_1(\xi, \theta)z + u_2(\xi, \theta)\frac{z^2}{2} + u_3(\xi, \theta)\frac{z^3}{6}$$ $$v(\xi, \theta, z) = v_0(\xi, \theta) + v_1(\xi, \theta)z + v_2(\xi, \theta)\frac{z^2}{2} + v_3(\xi, \theta)\frac{z^3}{6}$$ $$w(\xi, \theta, z) = w_0(\xi, \theta) + w_1(\xi, \theta)z + w_2(\xi, \theta)\frac{z^2}{2}$$

Mô hình chuyển vị này đóng góp 3 bước đột phá lý thuyết:

  1. Phục hồi trạng thái ứng suất 3 chiều đầy đủ: Cho phép biến dạng pháp tuyến ngang $\varepsilon_z = \frac{\partial w}{\partial z} = w_1 + w_2 z \neq 0$, từ đó kích hoạt trực tiếp thành phần ứng suất pháp ngang $\sigma_z$ thông qua định luật Hooke trực hướng mở rộng: $$\sigma_z = C_{31}(T, z)\varepsilon_\xi + C_{32}(T, z)\varepsilon_\theta + C_{33}(T, z)\varepsilon_z - \beta_z(T, z)\Delta T$$
  2. Loại bỏ hiện tượng khóa cắt (shear locking) mà không cần hệ số hiệu chỉnh: Các hàm bậc cao $z^2, z^3$ đảm bảo sự biến thiên parabol của ứng suất tiếp trượt $\tau_{\xi z}, \tau_{\theta z}$ trên mặt cắt ngang.
  3. Cơ chế cân bằng nội tại vi mô: Kết hợp hệ phương trình vi phân đàn hồi 3D để hiệu chỉnh ứng suất tiếp, thỏa mãn tuyệt đối điều kiện biên ứng suất mặt trong và mặt ngoài: $\tau_{\xi z}|{z=\pm h/2} = 0, \tau{\theta z}|_{z=\pm h/2} = 0$.

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

Khung phân tích của luận án tích hợp liền mạch 3 mô hình lý thuyết cốt lõi:

  • Quy luật hỗn hợp mở rộng (EROM): Xác định các mô-đun đàn hồi $E_{11}(T, z), E_{22}(T, z), G_{12}(T, z)$ và hệ số giãn nở nhiệt $\alpha_{11}(T, z), \alpha_{22}(T, z)$ phụ thuộc tọa độ $z$ và nhiệt độ $T(z)$, được hiệu chuẩn qua 3 hệ số hiệu dụng $\eta_1, \eta_2, \eta_3$ trích xuất từ dữ liệu mô phỏng MD của SWCNT $(10,10)$: $$E_{11}(T, z) = \eta_1 V_{CNT}(z) E_{11}^{CNT}(T) + V_m(z) E_m(T)$$ $$\frac{\eta_2}{E_{22}(T, z)} = \frac{V_{CNT}(z)}{E_{22}^{CNT}(T)} + \frac{V_m(z)}{E_m(T)}$$
  • Mô hình truyền nhiệt dừng 1 chiều phi tuyến: Phân bố nhiệt độ $T(z)$ không áp đặt trước mà được giải trực tiếp từ phương trình vi phân truyền nhiệt: $$\frac{d}{dz}\left[ k_{22}(z, T) \frac{dT(z)}{dz} \right] = 0, \quad \text{với } T(-h/2) = T_{in}, \quad T(h/2) = T_{out}$$
  • Nguyên lý biến phân công ảo (Principle of Virtual Displacements): Thiết lập hệ phương trình cân bằng tĩnh học tổng quát thông qua biến phân năng lượng: $$\delta U - \delta A_s - \delta A_{e1} - \delta A_{e2} = 0$$

Điều kiện biên điều chỉnh được xác định hoàn toàn tường minh, bao gồm 11 thành phần nội lực tổng quát ($N_i, M_i, N^_i, M^_i, Q_i, S_i, Q^*_i$) liên kết với 11 biến chuyển vị suy rộng ($u_0, u_1, u_2, u_3, v_0, v_1, v_2, v_3, w_0, w_1, w_2$).

Sơ đồ khung phân tích liên hợp cơ - nhiệt của luận án:
┌────────────────────────────────────────────────────────────────────────┐
│                        THÔNG SỐ ĐẦU VÀO CƠ - NHIỆT                     │
│  - Hình học: R, L, h (R/h, L/R)       - Nhiệt độ: Tin, Tout, Tref      │
│  - Vật liệu: SWCNT (10,10) + Epoxy    - Tải trọng: q+(ξ, θ), q-(ξ, θ)  │
│  - Phân bố CNT: UD, FG-Λ, FG-V, FG-X, FG-O; Tỷ lệ thể tích VCNT        │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│             ĐỒNG NHẤT HÓA VẬT LIỆU VÀ PHÂN BỐ NHIỆT ĐỘ T(z)            │
│  - Truyền nhiệt dừng Fourier: d/dz [k22(z, T) dT/dz] = 0 ──► T(z)      │
│  - EROM phụ thuộc nhiệt độ: E11(T,z), E22(T,z), G12(T,z), αii(T,z)     │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│               MÔ HÌNH QUASI-3D (K3-2) & NGUYÊN LÝ CÔNG ẢO              │
│  - Trường chuyển vị 11 ẩn số: u(ξ,θ,z), v(ξ,θ,z), w(ξ,θ,z) (chứa z^2) │
│  - Quan hệ biến dạng - chuyển vị & Ứng suất Hooke trực hướng 3D         │
│  - Nguyên lý công ảo: δU - δAs - δAe1 - δAe2 = 0                       │
│    ──► Hệ 11 phương trình vi phân đạo hàm riêng (PDEs) cân bằng tĩnh    │
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│               PHƯƠNG PHÁP GIẢI TÍCH KÍN TOÀN PHẦN                      │
│  - Khai triển chuỗi Fourier theo phương vòng θ                         │
│    ──► Chuyển hệ PDEs thành hệ ODEs bậc cao theo ξ                     │
│  - Phép biến đổi tích phân Laplace L{·} theo tọa độ đường sinh ξ       │
│    ──► Giải đại số ma trận trong miền ảnh s & Biến đổi ngược Laplace   │
│  - Áp đặt điều kiện biên tùy ý: C-C, C-F, C-S, S-S tại ξ = 0 và ξ = L/R│
└───────────────────────────────────┬────────────────────────────────────┘
                                    │
                                    ▼
┌────────────────────────────────────────────────────────────────────────┐
│                       KẾT QUẢ ĐẦU RA VÀ ĐỐI CHUẨN                      │
│  - Trường chuyển vị: u(ξ,θ,z), v(ξ,θ,z), w(ξ,θ,z)                       │
│  - Trường ứng suất 3D: σξ, σθ, σz, τξθ, τξz, τθz                      │
│  - Phân tích hiệu ứng biên & Khuyến nghị thiết kế kết cấu hàng không    │
└────────────────────────────────────────────────────────────────────────┘

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 suy diễn chuẩn tắc (Deductive Positivism), tiếp cận bài toán cơ học vật rắn biến dạng bằng phương pháp giải tích toán học thuần túy (Analytical Approach). Quy trình nghiên cứu loại bỏ hoàn toàn các sai số làm trơn hình học, suy biến phần tử hay hội tụ số vốn là nhược điểm cố hữu của phương pháp phần tử hữu hạn (FEM) hoặc phương pháp lưới không cấu trúc khi xử lý gradient ứng suất lớn tại vùng biên.

Mô hình toán học kiểm soát đầy đủ 5 cấu hình phân bố tỷ lệ thể tích CNT qua chiều dày $z \in [-h/2, h/2]$:

  1. Dạng UD (Đều): $V_{CNT}(z) = V_{CNT}^*$
  2. Dạng FG-$\Lambda$ (Tập trung mặt trong): $V_{CNT}(z) = 2\left(\frac{1}{2} - \frac{z}{h}\right)V_{CNT}^*$
  3. Dạng FG-V (Tập trung mặt ngoài): $V_{CNT}(z) = 2\left(\frac{1}{2} + \frac{z}{h}\right)V_{CNT}^*$
  4. Dạng FG-O (Tập trung ở mặt trung hòa): $V_{CNT}(z) = 2\left(1 - \frac{2|z|}{h}\right)V_{CNT}^*$
  5. Dạng FG-X (Tập trung ở hai bề mặt biên): $V_{CNT}(z) = 4\frac{|z|}{h}V_{CNT}^*$

Trong đó $V_{CNT}^* = \frac{w_{CNT}}{w_{CNT} + (\rho_{CNT}/\rho_m) - (\rho_{CNT}/\rho_m)w_{CNT}}$ là phân số thể tích tương đương với tỷ lệ phần trăm khối lượng $w_{CNT}$.

Quy trình nghiên cứu rigorous

Quy trình giải quyết hệ phương trình vi phân chủ đạo bao gồm hai giai đoạn biến đổi giải tích chặt chẽ:

Giai đoạn 1: Phân tách biến bằng chuỗi lượng giác Fourier đơn. Chuyển vị và tải trọng cục bộ tuần hoàn theo hướng vòng $\theta$ được biểu diễn dưới dạng chuỗi trực giao Fourier: $$u_i(\xi, \theta) = \sum_{m=0}^{\infty} U_{im}(\xi)\cos(m\theta), \quad v_i(\xi, \theta) = \sum_{m=1}^{\infty} V_{im}(\xi)\sin(m\theta), \quad w_j(\xi, \theta) = \sum_{m=0}^{\infty} W_{jm}(\xi)\cos(m\theta)$$ $$(i = 0, 1, 2, 3; \quad j = 0, 1, 2)$$

Phép biến đổi này chuyển đổi hệ 11 phương trình vi phân đạo hàm riêng (PDEs) phức tạp thành hệ 11 phương trình vi phân thường (ODEs) tuyến tính bậc cao theo tọa độ không thứ nguyên $\xi = x/R$.

Giai đoạn 2: Giải tích hệ ODEs bằng phép biến đổi Laplace. Áp dụng toán tử biến đổi Laplace $\mathcal{L}{f(\xi)} = \bar{f}(s) = \int_0^\infty f(\xi)e^{-s\xi}d\xi$ lên hệ ODEs theo biến $\xi$, hệ vi phân được đưa về hệ phương trình đại số tuyến tính trong không gian ảnh Laplace: $$[A(s)]{\bar{X}_m(s)} = {\bar{P}_m(s)} + [B(s)]{X_m(0)}$$

Nghiệm chuyển vị tổng quát thu được bằng cách lấy biến đổi ngược Laplace ($\mathcal{L}^{-1}$): $${X_m(\xi)} = [K(\xi)]{X_m(0)} + {F_m(\xi)}$$ Các hằng số chưa biết tại gốc ${X_m(0)}$ (chuyển vị hoặc nội lực tại $\xi = 0$) được xác định chính xác thông qua việc giải hệ điều kiện biên tại hai đầu vỏ trụ ($\xi = 0$ và $\xi = L/R$):

  • Biên ngàm (Clamped - C): $u_0 = u_1 = u_2 = u_3 = v_0 = v_1 = v_2 = v_3 = w_0 = w_1 = w_2 = 0$
  • Biên tự do (Free - F): $N_\xi = M_\xi = N^_\xi = M^\xi = N{\xi\theta} = M_{\xi\theta} = N^_{\xi\theta} = M^{\xi\theta} = Q\xi = S_\xi = Q^*_\xi = 0$
  • Biên tựa bản lề (Simply supported - S): $v_0 = w_0 = N_\xi = M_\xi = N^_\xi = M^_\xi = 0$

Data và phân tích

Tính xác thực và độ tin cậy của thuật toán giải tích được chứng minh thông qua kiểm chuẩn 3 cấp độ (Triangulation of Validation):

  1. Kiểm chuẩn cơ học thuần túy: So sánh kết quả chuyển vị không thứ nguyên và ứng suất với công bố của Doan et al. [24] cho vỏ trụ FGM ngàm hai đầu (C-C) chịu áp suất trong $P_{in}$. Sai số tương đối nhỏ hơn $0.12%$.
  2. Kiểm chuẩn nhiệt cơ đối xứng: So sánh ứng suất nhiệt và độ võng với kết quả giải tích đàn hồi 3 chiều của Alibeigloo [18] và Moradi-Dastjerdi [81] cho panel trụ FG-CNTRC có biên tựa đơn (S-S). Sự trùng khớp đạt mức $99.8%$.
  3. Phân tích độ nhạy tham số: Chương trình tính toán được lập trình tự động hóa giải hàng trăm trường hợp thay đổi liên tục các tỷ số $L/R$ ($1$ đến $6$), $R/h$ ($10$ đến $50$), phân số thể tích $V_{CNT}^*$ ($0.11, 0.14, 0.17, 0.28$) và chênh lệch nhiệt độ $\Delta T = T_{out} - T_{in}$ ($0\text{ K}$ đến $200\text{ K}$).

Dữ liệu đặc tính vật liệu của SWCNT $(10,10)$ trích xuất trực tiếp từ bảng cơ sở dữ liệu của luận án:

  • Mật độ: $\rho_{CNT} = 1.4 \times 10^3\text{ kg/m}^3$, $\rho_m = 1.15 \times 10^3\text{ kg/m}^3$ (nhựa Epoxy).
  • Đặc tính cơ nhiệt phụ thuộc $T$: $E_{11}^{CNT}(T) = 6.18 - 2.86 \times 10^{-3}T + 4.228 \times 10^{-6}T^2 - 2.27 \times 10^{-9}T^3\text{ (TPa)}$.
  • Dẫn nhiệt: $k_{11}^{CNT} = 3000\text{ W/mK}$, $k_m = 0.21\text{ W/mK}$; Hệ số giãn nở nhiệt: $\alpha_{11}^{CNT} = -1.12 \times 10^{-6}\text{ K}^{-1}$.

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

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

So sánh đáp ứng cơ - nhiệt giữa các cấu hình phân bố CNT (Vỏ C-C, L/R = 2, R/h = 20, Tin = 300K, Tout = 400K):
┌────────────────┬──────────────────────────┬────────────────────────────┬─────────────────────────────┐
│ Cấu hình CNT   │ Độ võng cực đại w_max/h  │ Ứng suất vòng σθ_max (MPa) │ Ứng suất biên σξ_edge (MPa) │
├────────────────┼──────────────────────────┼────────────────────────────┼─────────────────────────────┤
│ UD (Đều)       │ 1.000 (Chuẩn đối sánh)   │ 1.000 (Chuẩn đối sánh)     │ 1.000 (Chuẩn đối sánh)      │
│ FG-O           │ +18.4% (Kém nhất)        │ +14.2%                     │ +22.7%                      │
│ FG-Λ           │ -12.6%                   │ -9.8%                      │ -15.4%                      │
│ FG-V           │ -28.5% (Rất tốt)         │ -21.3%                     │ -32.1%                      │
│ FG-X           │ -34.2% (Tối ưu nhất)     │ -26.7%                     │ -38.6%                      │
└────────────────┴──────────────────────────┴────────────────────────────┴─────────────────────────────┘
  1. Hiệu năng vượt trội của cấu hình phân bố FG-X và FG-V: Luận án chứng minh định lượng rằng cấu hình phân bố đối xứng FG-X (CNT tập trung tối đa ở hai bề mặt ngoài) và FG-V (CNT tập trung ở bề mặt ngoài) làm giảm độ võng cực đại tương ứng $34.2%$ và $28.5%$ so với cấu hình phân bố đều UD ở cùng tỷ lệ thể tích $V_{CNT}^* = 0.17$. Ngược lại, cấu hình FG-O làm tăng độ võng lên $18.4%$, biểu thị sự lãng phí tài nguyên nano khi đưa cốt gia cường vào vùng trung hòa chịu ứng suất uốn nhỏ.
  2. Bản chất của hiệu ứng biên (Edge Boundary Layer Effect): Dữ liệu giải tích chỉ ra rằng tại các liên kết ngàm (C-C, C-F), ứng suất dọc trục $\bar{\sigma}_\xi$ và ứng suất pháp ngang $\bar{\sigma}z$ tăng đột biến cục bộ trong một khoảng cách suy giảm biên hẹp $\xi_0 \approx \beta \sqrt{h/R}$. Tại vị trí sát mép ngàm $\xi = 0$, ứng suất dọc trục $\sigma\xi$ đạt giá trị cực đại gấp $3.2$ đến $4.8$ lần so với vùng giữa thân vỏ, khẳng định cơ chế phá hủy luôn tiềm ẩn tại các mép liên kết.
  3. Hiện tượng mềm hóa do nhiệt (Thermal Softening Degradation): Khi xét đến tính chất phụ thuộc nhiệt độ $E(T)$, độ võng hướng kính của vỏ tăng thêm từ $16.8%$ đến $31.4%$ so với mô hình giả định $E = \text{const}$. Sự suy giảm mô-đun đàn hồi của nền Epoxy ở nhiệt độ $T > 380\text{ K}$ làm chuyển dịch trục trung hòa biến dạng, làm trầm trọng thêm mức độ tập trung ứng suất trượt ngang $\tau_{\xi z}$.
  4. Ảnh hưởng phi đối xứng của gradient nhiệt độ: Gradient nhiệt độ thực truyền từ trong ra ngoài ($T_{in} > T_{out}$) tạo ra trường ứng suất nén nhiệt ở mặt trong và ứng suất kéo nhiệt ở mặt ngoài. Khi kết hợp với tải trọng áp suất trong, cấu hình FG-$\Lambda$ phát huy hiệu quả triệt tiêu ứng suất kéo cục bộ hiệu quả hơn cấu hình FG-V.

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính tất yếu của việc đưa thành phần $\sigma_z$ vào phân tích ứng xử tĩnh của vỏ composite dày, xác lập giới hạn ứng dụng của lý thuyết FSDT và bổ sung một nhánh mở rộng hoàn chỉnh cho cơ học vật liệu FGM nano.
  • Về phương pháp luận: Thiết lập khung thuật toán tích hợp Fourier - Laplace giải quyết triệt để bài toán biên hỗn hợp không đối xứng cho hệ vi phân liên hợp bậc cao, có thể chuyển giao trực tiếp để phân tích các dạng kết cấu khác như vỏ nón, vỏ cầu, tấm sandwich nhiều lớp.
  • Về ứng dụng thực tiễn: Cung cấp giải pháp tối ưu hóa hình học phân bố CNT: khuyến nghị ưu tiên sử dụng cấu hình FG-X cho các kết cấu chịu uốn - nén hỗn hợp và cấu hình FG-V/FG-$\Lambda$ cho các ống dẫn khí nóng, thân động cơ tên lửa chịu gradient nhiệt hướng tâm một chiều.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ rõ 4 giới hạn nghiên cứu nội tại:

  1. Giới hạn đàn hồi tuyến tính: Vật liệu được giả thiết làm việc hoàn toàn trong miền đàn hồi tuyến tính của định luật Hooke trực hướng, chưa tính đến biến dạng dẻo hoặc tổn thương tích lũy của nền polymer khi vượt quá nhiệt độ chuyển pha thủy tinh ($T_g$).
  2. Giả thiết liên kết pha lý tưởng: Quy luật hỗn hợp mở rộng (EROM) ngầm định liên kết giữa ống nano carbon và nền Epoxy là hoàn hảo tuyệt đối (perfect interfacial bonding), bỏ qua hiện tượng trượt tiếp xúc, tạo cụm nano (CNT agglomeration) hoặc vi nứt tế vi tại bề mặt tiếp xúc.
  3. Mô hình hình học lý tưởng: Phân tích thực hiện trên vỏ trụ tròn đều trơn nhẵn, chưa xét đến các bất toàn hình học ban đầu (initial geometric imperfections), lỗ khoét kỹ thuật hoặc các gân tăng cứng (stiffeners) thường gặp trong kết cấu hàng không thực tế.
  4. Truyền nhiệt dừng 1 chiều: Trường nhiệt độ chỉ được xét biến thiên qua chiều dày $T(z)$, chưa giải bài toán truyền nhiệt 2 chiều hoặc 3 chiều không dừng ($T(\xi, \theta, z, t)$) dưới tác dụng của xung nhiệt động học.

Chương trình nghiên cứu tương lai (5-10 năm tới) đề xuất 4 hướng phát triển:

  • Mở rộng mô hình Quasi-3D cho bài toán phân tích phi tuyến hình học lớn (geometrically non-linear postbuckling) và động lực học va chạm va đập tốc độ cao.
  • Tích hợp hiệu ứng tập hợp hạt nano (CNT agglomeration) dựa trên lý thuyết vi cơ học Mori-Tanaka cải tiến kết hợp mô phỏng đa tỷ lệ (multiscale modeling).
  • Phát triển phương pháp giải tích cho vỏ trụ composite thông minh tích hợp các lớp cảm biến/chấp hành áp điện (piezoelectric layers) và vật liệu từ-điện-nhiệt đàn hồi để điều khiển chủ động ứng suất biên.
  • Tiến hành thực nghiệm chế tạo mẫu vỏ FG-CNTRC bằng công nghệ in 3D nano và đo đạc ứng suất tập trung bằng phương pháp tương quan ảnh số (DIC) để kiểm chứng thực nghiệm.

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

  • Tác động học thuật: Luận án mở ra một chuẩn đối sánh giải tích (analytical benchmark solutions) độ chính xác cao cho cộng đồng cơ học tính toán quốc tế. Các công thức đóng và bảng số liệu của luận án là công cụ tiêu chuẩn để kiểm chứng độ chính xác và tốc độ hội tụ của các phần mềm thương mại (ABAQUS, ANSYS) khi phát triển các phần tử hữu hạn vỏ đặc biệt.
  • Chuyển đổi công nghiệp hàng không & vũ trụ: Cung cấp cơ sở khoa học định lượng cho các viện nghiên cứu thiết kế tên lửa, khí tài bay và ống phóng quân sự trong việc tối ưu hóa tỷ lệ pha tạp nano, giúp giảm $20-30%$ khối lượng kết cấu chịu nhiệt mà vẫn đảm bảo độ bền biên an toàn.
  • Lợi ích xã hội & quốc phòng: Đóng góp trực tiếp vào mục tiêu làm chủ công nghệ vật liệu tiên tiến tự chủ, nâng cao tuổi thọ và độ tin cậy của các hệ thống vũ khí công nghệ cao trong điều kiện hoạt động khắc nghiệt.

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

  • Nghiên cứu sinh & Giảng viên ngành Cơ học/Kỹ thuật Hàng không: Tiếp cận phương pháp toán giải tích chuẩn mực (chuỗi Fourier kết hợp biến đổi Laplace) để xử lý các hệ vi phân đạo hàm riêng bậc cao trong cơ học kết cấu.
  • Kỹ sư R&D vật liệu Composite: Nắm bắt quy luật thiết kế phân bố CNT theo chiều dày (FG-X, FG-V) nhằm định hướng công nghệ chế tạo màng vật liệu biến tính tối ưu khả năng chịu nhiệt.
  • Chuyên gia thiết kế kết cấu thiết bị bay: Nhận diện chính xác vị trí và biên độ tập trung ứng suất tại các khớp nối ngàm, từ đó đưa ra giải pháp gia cường kết cấu biên cục bộ hợp lý, ngăn ngừa hiện tượng phá hủy mỏi và tách lớp nhiệ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 độc đáo nhất là việc thiết lập thành công mô hình giải tích tĩnh Quasi-3D (mô hình K3-2) mở rộng từ lý thuyết biến dạng cắt bậc cao của Reddy kết hợp phương trình vi cơ học EROM phụ thuộc nhiệt độ. Khác biệt cốt tử nằm ở việc khai triển chuyển vị pháp tuyến chứa số hạng bậc hai theo chiều dày ($w = w_0 + w_1 z + w_2 z^2/2$), cho phép trực tiếp tính toán ứng suất pháp ngang $\sigma_z$ và ứng suất trượt 3D liên tục, điều mà các lý thuyết tấm vỏ cổ điển CST và FSDT hoàn toàn không thực hiện được.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm? Trả lời: So với các nghiên cứu giải tích của Shen & Xiang [43, 50] và Alibeigloo [18] (chỉ giải được cho biên tựa đơn S-S bằng nghiệm Navier kép), luận án kết hợp chuỗi Fourier đơn theo phương vòng với phép biến đổi Laplace toàn phần theo phương trục. Kỹ thuật này chuyển hệ PDEs thành hệ phương trình đại số trong không gian ảnh $s$, cho phép giải chính xác nghiệm kín cho các điều kiện biên phi đối xứng phức tạp (C-C, C-F, C-S) chịu tải trọng cục bộ phân bố bất kỳ.

3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì? Trả lời: Dữ liệu chỉ ra rằng việc tăng tỷ lệ thể tích $V_{CNT}$ không phải lúc nào cũng làm giảm ứng suất tại mọi vị trí. Tại khu vực sát mép ngàm $\xi \approx 0$, việc tăng $V_{CNT}$ ở cấu hình FG-V làm tăng độ cứng chống uốn cục bộ quá mức, dẫn đến mức độ tập trung ứng suất uốn dọc trục $\sigma_\xi$ tại biên tăng cao hơn $24%$ so với khi phân bố đều, tạo ra nguy cơ nứt vỡ giòn cục bộ nếu không có biện pháp vát mép chuyển tiếp.

4. Quy trình kiểm chuẩn và tái lập kết quả (Replication Protocol) có được đảm bảo? Trả lời: Hoàn toàn minh bạch. Toàn bộ các hệ số độ cứng $C_{ij}(T,z)$, các phương trình biến phân công ảo, ma trận chuyển đổi Laplace và bảng thông số vật liệu phụ thuộc nhiệt độ của SWCNT $(10,10)$ và Epoxy đều được cung cấp tường minh trong các chương 2, 3 và 4 của luận án, cho phép tái lập độc lập 100% kết quả trên bất kỳ nền tảng tính toán ký hiệu nào (Mathematica, MATLAB, Maple).

5. Chương trình nghị sự nghiên cứu 10 năm tới được định hình ra sao? Trả lời: Luận án vạch ra lộ trình mở rộng từ phân tích tĩnh sang phân tích ổn định động phi tuyến (Dynamic Postbuckling), phân tích mỏi nhiệt cơ (Thermal-Mechanical Fatigue) có xét đến hiệu ứng kích thước nano phi cục bộ (Nonlocal Strain Gradient Elasticity) và tiến tới tối ưu hóa đa mục tiêu cấu trúc liên kết nano FG-CNTRC trong điều kiện bay siêu âm ($M > 5$).

Kết luận

Công trình nghiên cứu của NCS Dương Văn Quang đã hoàn thành xuất sắc các mục tiêu đề ra với 5 đóng góp học thuật cốt lõi:

  1. Thiết lập mô hình toán học Quasi-3D K3-2 hoàn chỉnh cho vỏ trụ FG-CNTRC chịu tải trọng cơ - nhiệt liên hợp, tính toán đầy đủ thành phần ứng suất pháp tuyến ngang $\sigma_z$.
  2. Xây dựng phương pháp giải tích giải thuật kín dựa trên chuỗi Fourier và biến đổi Laplace, phá vỡ giới hạn truyền thống chỉ giải được biên tựa đơn giản.
  3. Mô hình hóa chính xác gradient nhiệt dừng 1D và sự suy giảm đặc tính cơ lý theo nhiệt độ dựa trên quy luật hỗn hợp mở rộng hiệu chuẩn bằng mô phỏng MD.
  4. Khám phá và định lượng quy luật phân bố ứng suất vùng hiệu ứng biên, xác lập cấu hình phân bố FG-X và FG-V là giải pháp tối ưu hóa kết cấu chịu uốn nhiệt.
  5. Cung cấp hệ thống dữ liệu đối chuẩn giải tích chuẩn xác cao, đóng góp thiết thực cho sự phát triển của cơ học tính toán và công nghệ hàng không vũ trụ tại Việt Nam.

Luận án khẳng định bước tiến vững chắc trong việc làm chủ lý thuyết cơ học kết cấu tiên tiến, mở ra các hướng nghiên cứu liên ngành đột phá về vật liệu nanocomposite thông minh trong tương lai.