Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Vật liệu composite cơ tính biến thiên được gia cường bằng ống nano carbon (Functionally Graded Carbon Nanotube Reinforced Composites - FG-CNTRC) là loại vật liệu composite mới với các đặc tính như độ bền cao, khối lượng riêng nhỏ, cùng tính chất cơ nhiệt và dẫn điện tốt. Trong lĩnh vực kỹ thuật hàng không vũ trụ, vỏ trụ là dạng kết cấu cơ bản trên các thiết bị bay, thường xuyên chịu tác động đồng thời của tải trọng cơ học và trường nhiệt độ.

Tổng quan tình hình nghiên cứu cho thấy phần lớn các công trình công bố về ứng xử nhiệt đàn hồi của kết cấu FG-CNTRC giả thiết tính chất vật liệu không phụ thuộc vào nhiệt độ (temperature-independent properties). Một số ít nghiên cứu có xét đến tính chất phụ thuộc nhiệt độ (temperature-dependent properties) nhưng lại chưa đánh giá ảnh hưởng của các điều kiện biên khác nhau đến bài toán nhiệt đàn hồi, hoặc chỉ khảo sát trường nhiệt độ đồng nhất. Bên cạnh đó, các phương pháp giải tích truyền thống thường giới hạn ở điều kiện biên tựa đơn giản (S-S). Hơn nữa, đa số các lý thuyết tấm vỏ bậc cao thông dụng thường bỏ qua ứng suất pháp tuyến ngang theo chiều dày vỏ, dẫn đến hạn chế khi phân tích trạng thái ứng suất - biến dạng tại các khu vực chuyển tiếp, vùng biên và vùng chịu tải tập trung - vốn là những vị trí dễ phát sinh phá hủy trong thực tế. Do đó, việc nghiên cứu phân tích tĩnh kết cấu vỏ trụ FG-CNTRC chịu tải trọng cơ và nhiệt độ bằng lý thuyết biến dạng cắt bậc cao có kể đến ứng suất pháp tuyến ngang và áp dụng phương pháp giải tích cho các điều kiện biên khác nhau là vấn đề khoa học có tính cấp thiết.

Mục tiêu nghiên cứu

Luận án xác định hai mục tiêu nghiên cứu cụ thể:

  1. Xây dựng mô hình toán học, phương pháp giải tích và chương trình tính toán tin cậy để phân tích tĩnh vỏ trụ bằng vật liệu FG-CNTRC chịu tác dụng đồng thời của tải trọng cơ học và nhiệt độ, góp phần bổ sung và hoàn thiện mô hình tính toán ứng xử cơ học của kết cấu FG-CNTRC.
  2. Khảo sát, đánh giá định lượng và định tính đáp ứng tĩnh của vỏ trụ FG-CNTRC dưới các điều kiện làm việc khác nhau (điều kiện biên, quy luật phân bố CNT, tỷ lệ thể tích CNT, các thông số hình học và dạng tải trọng), từ đó đưa ra các khuyến cáo khoa học phục vụ tính toán, thiết kế và khai thác kết cấu trong thực tiễn.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Kết cấu vỏ trụ tròn xoay làm bằng vật liệu nanocomposite cơ tính biến thiên gia cường bằng ống nano carbon (FG-CNTRC), chịu tác dụng của tải trọng cơ học và tải trọng nhiệt độ.
  • Phạm vi nghiên cứu:
    • Kết cấu vỏ trụ có các điều kiện biên khác nhau (ngàm hai đầu C-C, một đầu ngàm một đầu tự do C-F, một đầu ngàm một đầu tựa C-S, tựa hai đầu S-S) và các thông số hình học khác nhau.
    • Vật liệu cốt gia cường là ống nano carbon đơn vách (SWCNT) dạng thẳng, phân bố đồng phương theo trục dọc của vỏ, có tỷ lệ thể tích biến thiên theo chiều dày vỏ theo 5 quy luật (UD, FG-Ʌ, FG-V, FG-O, FG-X).
    • Trạng thái biến dạng của vật liệu nằm trong giới hạn đàn hồi tuyến tính.
    • Mô hình cơ học sử dụng lý thuyết biến dạng cắt bậc cao có xét đến ứng suất pháp tuyến ngang (mô hình quasi-3D ký hiệu K3-2) kết hợp với các đặc trưng cơ lý của vật liệu phụ thuộc vào nhiệt độ.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước

Luận án đã hệ thống hóa các công trình nghiên cứu trong và ngoài nước liên quan đến vật liệu FG-CNTRC và ứng xử cơ học của kết cấu tấm vỏ:

  1. Nghiên cứu về đặc tính của ống nano carbon (CNT) và vật liệu composite gia cường CNT:
    • Iijima (1991) phát hiện ra CNT. Tiếp đó, Overney và cộng sự (1993), Wong và cộng sự (1997), Salvetat và cộng sự (1999), Lu (1997) xác định các đặc trưng cơ học vượt trội của CNT với mô-đun Young trên $1\text{ TPa}$ và độ giãn dài tới $12%$.
    • Các công trình thực nghiệm của Shaffer và Windle (1999), Qian và cộng sự (2000), Safadi và cộng sự (2002), Cadek và cộng sự (2004), Ruan và cộng sự (2006), Pötschke và cộng sự (2002) chứng minh việc bổ sung CNT vào nền polymer giúp cải thiện độ bền kéo, mô-đun đàn hồi và độ dẫn nhiệt.
    • Shen (2009) lần đầu đề xuất mô hình vật liệu nanocomposite có cơ tính biến thiên (FG-CNTRC) với các dạng phân bố tỷ lệ thể tích CNT theo chiều dày, được phát triển trong các công trình tiếp theo (2011, 2012, 2014).
Vật liệu Khối lượng riêng $\rho$ ($\text{g/cm}^3$) Mô-đun Young $E$ ($\text{GPa}$) Độ bền kéo $\sigma_m$ ($\text{GPa}$) Độ giãn dài tới hạn $\varepsilon$ ($%$) Độ dẫn nhiệt $k$ ($\text{W/mK}$) Độ dẫn điện $\sigma$ ($\text{S/m}$)
SWCNT 1,6 1200 150 12 > 3000 $10^6 - 10^7$
Sợi carbon M60JB 1,7 588 3,8 0,6 140 - 800 $10^4 - 8,5\times 10^6$
Sợi thủy tinh kiểu S 2,46 85 4,6 5,4 0,13 $10^{-14}$
Thép 17-7 PH RH950 7,005 200 1,5 4 80 -

  1. Nghiên cứu ứng xử tấm vỏ FG-CNTRC chịu tải trọng cơ và nhiệt:

    • Shen, Xiang, Wang và Zhu (2010–2014) khảo sát uốn phi tuyến, dao động và mất ổn định nhiệt của tấm và vỏ trụ FG-CNTRC trên nền đàn hồi.
    • Các tác giả trong nước: Lê Thị Như Trang (2017) khảo sát ổn định nhiệt đàn hồi của tấm và panel FG-CNTRC; Vũ Thanh Long (2018) nghiên cứu tấm và panel sandwich FG-CNTRC bằng lý thuyết FSDT; Nguyễn Văn Thành (2020) khảo sát đáp ứng động và phi tuyến bằng HSDT của Reddy; Phạm Thanh Hiếu (2021) phân tích ổn định vỏ trụ và vỏ trống FG-CNTRC theo CST và FSDT; Vũ Văn Thẩm (2018) sử dụng lý thuyết tấm bốn ẩn cải tiến (HSDT4) cho kết cấu tấm vỏ có gắn lớp áp điện; Chung và cộng sự khảo sát tấm FG-CNTRC có gân chịu tải nổ; Dat và cộng sự nghiên cứu dao động tấm sandwich CNTRC.
    • Về tải nhiệt: Alibeigloo sử dụng lý thuyết đàn hồi 3 chiều cho tấm và panel tựa đơn; Pourasghar và Chen (2016), Moradi-Dastjerdi và Payganeh khảo sát vỏ trụ có xét gradient nhiệt và cơ tính phụ thuộc nhiệt độ cho biên tựa đơn.

  2. Lý thuyết kết cấu tấm vỏ:

    • Lý thuyết cổ điển (CST/CPT): Đã được áp dụng bởi Ninh và Bích, Văn Tùng và Trang, Phạm Thanh Hiếu, Đỗ Đức và cộng sự, nhưng bỏ qua biến dạng trượt ngang.
    • Lý thuyết biến dạng cắt bậc nhất (FSDT): Áp dụng bởi Liew, Mehrabadi, Malekzadeh, Lei, Ansari, Long và Tùng, có tính đến trượt ngang nhưng đòi hỏi hệ số hiệu chỉnh cắt.
    • Lý thuyết biến dạng cắt bậc cao (HSDT): Phùng Văn và cộng sự, Natarajan, Karami, Soni, Shen phát triển nhằm phản ánh chính xác phân bố ứng suất cắt.
    • Lý thuyết quasi-3D có kể đến ứng suất pháp ngang ($\sigma_z$): Được Đoàn Trắc Luật và cộng sự phát triển cho vỏ composite lớp và vỏ FGM.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

  • Đa số nghiên cứu giải tích chỉ giải quyết điều kiện biên tựa đơn giản (S-S), thiếu phương pháp giải tích tổng quát cho các điều kiện biên ngàm (C) và tự do (F).
  • Hầu hết các bài toán nhiệt đàn hồi giả thiết trường nhiệt độ đồng nhất hoặc dạng hàm tuyến tính cho trước; rất ít nghiên cứu giải bài toán truyền nhiệt một chiều theo chiều dày để xác định chính xác phân bố nhiệt độ phi tuyến và đồng thời xét tính chất cơ lý phụ thuộc nhiệt độ.
  • Các lý thuyết HSDT phổ biến bỏ qua thành phần ứng suất pháp tuyến ngang ($\sigma_z = 0$), không phản ánh đầy đủ trạng thái ứng suất tại vùng lân cận liên kết biên của vỏ dày.
  • Thiếu các nghiên cứu đánh giá định lượng và chi tiết hiện tượng tập trung ứng suất (hiệu ứng biên) của vỏ trụ FG-CNTRC dưới tác dụng của tải trọng cục bộ và tải trọng cơ nhiệt kết hợp.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và mô hình toán học

  1. Trường chuyển vị và mô hình quasi-3D (K3-2): Trường chuyển vị ba chiều $(u, v, w)$ trong hệ tọa độ cong trực giao $O\xi\theta z$ được khai triển theo đa thức bậc cao theo tọa độ chiều dày $z$: $$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}$$ Trong đó chuyển vị pháp tuyến $w$ là hàm bậc hai theo $z$, cho phép xác định thành phần biến dạng pháp ngang $\varepsilon_z \neq 0$ và ứng suất pháp ngang $\sigma_z \neq 0$ (khác với mô hình K3-0 có $w = w_0$).

  2. Định luật Hooke cho vật liệu trực hướng chịu nhiệt: Mối quan hệ giữa ứng suất và biến dạng có xét đến nhiệt độ được biểu diễn qua ma trận độ cứng $C_{ij}(T, z)$: $$\begin{bmatrix} \sigma_\xi \ \sigma_\theta \ \sigma_z \ \sigma_{\theta z} \ \sigma_{\xi z} \ \sigma_{\xi\theta} \end{bmatrix} = \begin{bmatrix} C_{11} & C_{12} & C_{13} & 0 & 0 & 0 \ C_{21} & C_{22} & C_{23} & 0 & 0 & 0 \ C_{31} & C_{32} & C_{33} & 0 & 0 & 0 \ 0 & 0 & 0 & C_{44} & 0 & 0 \ 0 & 0 & 0 & 0 & C_{55} & 0 \ 0 & 0 & 0 & 0 & 0 & C_{66} \end{bmatrix} \begin{bmatrix} \varepsilon_\xi - \alpha_{11}\Delta T \ \varepsilon_\theta - \alpha_{22}\Delta T \ \varepsilon_z - \alpha_{33}\Delta T \ \gamma_{\theta z} \ \gamma_{\xi z} \ \gamma_{\xi\theta} \end{bmatrix}$$ với $\Delta T = T(\xi, \theta, z) - T_{ref}$.

  3. Mô hình đồng nhất hóa vật liệu (Extended Rule of Mixture): Các đặc trưng cơ lý tương đương phụ thuộc nhiệt độ được xác định theo quy luật hỗn hợp mở rộng: $$E_{11}(T, z) = \eta_1 V_{CNT} E_{11}^{CNT}(T) + V_m E_m(T)$$ $$\frac{\eta_2}{E_{22}(T, z)} = \frac{V_{CNT}}{E_{22}^{CNT}(T)} + \frac{V_m}{E_m(T)}$$ $$\frac{\eta_3}{G_{12}(T, z)} = \frac{V_{CNT}}{G_{12}^{CNT}(T)} + \frac{V_m}{G_m(T)}$$ $$\nu_{12}(z) = V_{CNT}\nu_{12}^{CNT} + V_m\nu_m, \quad \rho(z) = V_{CNT}\rho^{CNT} + V_m\rho^m$$ $$\alpha_{11}(T, z) = \frac{V_{CNT}E_{11}^{CNT}\alpha_{11}^{CNT} + V_m E_m\alpha^m}{E_{11}(T, z)}$$ $$\alpha_{22}(T, z) = (1 + \nu_{12}^{CNT})V_{CNT}\alpha_{22}^{CNT} + (1 + \nu_m)V_m\alpha^m - \nu_{12}(z)\alpha_{11}(T, z)$$ trong đó $\eta_1, \eta_2, \eta_3$ là các hệ số hiệu dụng của CNT thu được từ đối chiếu với mô phỏng động lực học phân tử (MD). Tỷ lệ thể tích thỏa mãn $V_{CNT}(z) + V_m(z) = 1$ theo 5 dạng phân bố: UD, FG-Ʌ, FG-V, FG-O, FG-X.

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

  • Phương pháp tiếp cận: Nghiên cứu lý thuyết kết hợp giải tích và tính toán số.
  • Thiết lập phương trình cân bằng: Ứng dụng Nguyên lý công ảo ($\delta U - \delta A_s - \delta A_{e1} - \delta A_{e2} = 0$) để thiết lập hệ phương trình đạo hàm riêng (PDE) 11 ẩn chuyển vị $(u_0, u_1, u_2, u_3, v_0, v_1, v_2, v_3, w_0, w_1, w_2)$ cùng hệ thống các điều kiện biên tương ứng.
  • Kỹ thuật giải:
    1. Sử dụng chuỗi lượng giác đơn theo hướng vòng $\theta$ để biến đổi hệ phương trình đạo hàm riêng (PDE) thành hệ phương trình vi phân thường (ODE) theo biến tọa độ đường sinh $\xi$.
    2. Áp dụng phép biến đổi Laplace theo $\xi$ để chuyển hệ ODE sang không gian ảnh đại số, giải hệ đại số để tìm các hằng số tích phân từ điều kiện biên, sau đó thực hiện biến đổi Laplace ngược để thu được nghiệm giải tích chính xác.
    3. Để tính toán chính xác các thành phần ứng suất cắt $\sigma_{\xi z}, \sigma_{\theta z}$, luận án sử dụng các phương trình cân bằng vi phân của lý thuyết đàn hồi 3D để tích phân theo chiều dày, đảm bảo thỏa mãn điều kiện biên ứng suất tại mặt trong và mặt ngoài của vỏ.

Nội dung chính theo từng chương

Chương 1. TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU

Chương 1 trình bày tổng quan cấu trúc, đặc tính của CNT đơn vách (SWCNT) và đa vách (MWCNT), các cơ chế tương tác phân tử giữa cốt CNT và nền polymer. Chương này phân tích chi tiết 5 mô hình phân bố tỷ lệ thể tích CNT theo chiều dày (UD, FG-Ʌ, FG-V, FG-O, FG-X) do Shen đề xuất. Đồng thời, tác giả tổng hợp tình hình nghiên cứu kết cấu tấm vỏ FG-CNTRC trong và ngoài nước theo ba nhóm vấn đề chính: ảnh hưởng của tải trọng nhiệt và tính chất phụ thuộc nhiệt độ; các lý thuyết tấm vỏ từ cổ điển (CST), biến dạng cắt bậc nhất (FSDT), biến dạng cắt bậc cao (HSDT) đến lý thuyết đàn hồi 3D; và các phương pháp giải bài toán tấm vỏ. Trên cơ sở chỉ ra các hạn chế của các công bố trước đó, chương 1 xác định rõ mục tiêu, đối tượng, phạm vi và 3 nội dung nghiên cứu trọng tâm của đề tài luận án.

Chương 2. XÂY DỰNG MÔ HÌNH TÍNH TOÁN VÀ PHƯƠNG PHÁP GIẢI

Chương 2 thiết lập bài toán tổng quát cho vỏ trụ FG-CNTRC có chiều dài $L$, bán kính trung bình $R$, chiều dày $h$ trong hệ tọa độ $O\xi\theta z$. Nội dung chương gồm:

  • Xác định hệ thức phụ thuộc nhiệt độ của vật liệu composite cốt SWCNT (10,10) và nền polymer theo Quy luật hỗn hợp mở rộng.
  • Xây dựng quan hệ hình học biến dạng - chuyển vị và quan hệ vật lý ứng suất - biến dạng theo lý thuyết biến dạng cắt bậc cao quasi-3D (mô hình K3-2 gồm 11 thành phần chuyển vị).
  • Áp dụng nguyên lý công ảo để dẫn xuất tường minh hệ 11 phương trình vi phân cân bằng cơ nhiệt và hệ điều kiện biên tại hai đầu vỏ ($\xi = 0$ và $\xi = L/R$).
  • Xây dựng quy trình giải tích hai bước: khai triển các hàm chuyển vị và tải trọng theo chuỗi Fourier lượng giác hướng vòng $\theta$ để khử đạo hàm theo $\theta$, tiếp đó áp dụng phép biến đổi Laplace theo tọa độ trục $\xi$ để tìm biểu thức nghiệm giải tích dạng tường minh, cho phép áp dụng đồng thời cho nhiều tổ hợp điều kiện biên (C-C, C-F, C-S, S-S) chịu tải trọng cục bộ phân bố bất kỳ.
  • Xây dựng thuật toán và chương trình tính toán số trên máy tính.

Chương 3. NGHIÊN CỨU KHẢO SÁT VỎ TRỤ FG-CNTRC CHỊU TẢI TRỌNG CƠ

Chương 3 tập trung khảo sát đáp ứng tĩnh của vỏ trụ FG-CNTRC chịu tải trọng cơ học hướng kính:

  • Kiểm chứng độ tin cậy: So sánh kết quả chuyển vị và ứng suất tính toán với các công bố của tác giả khác cho bài toán vỏ FGM và vỏ FG-CNTRC, khẳng định độ chính xác của mô hình và chương trình tính.
  • Khảo sát ảnh hưởng của điều kiện biên: So sánh 4 loại liên kết biên (C-C, C-F, C-S, S-S). Kết quả cho thấy chuyển vị và ứng suất phụ thuộc chặt chẽ vào độ cứng liên kết biên; biên C-C cho chuyển vị bé nhất, biên C-F cho chuyển vị lớn nhất.
  • Hiệu ứng biên và tập trung ứng suất: Khảo sát chi tiết vùng lân cận biên ngàm ($\xi = 0$ và $\xi = L/R$). Kết quả chỉ ra rằng ứng suất pháp dọc trục $\sigma_\xi$ và ứng suất tiếp $\sigma_{\xi z}$ tăng vọt cục bộ tại sát mép ngàm và suy giảm nhanh khi đi vào vùng giữa nhịp vỏ.
  • Ảnh hưởng của thông số vật liệu:
    • Kiểu phân bố CNT: Dạng FG-V và FG-X làm giảm đáng kể ứng suất tại bề mặt chịu tải so với dạng UD và FG-Ʌ.
    • Tỷ lệ thể tích $V^_{CNT}$: Khi tăng $V^_{CNT}$ (khảo sát từ $0,12$ đến $0,28$), độ võng của vỏ giảm rõ rệt và độ cứng kết cấu tăng lên.
  • Ảnh hưởng của tải trọng: Đánh giá mức độ tập trung tải trọng, vị trí đặt tải trọng cục bộ và dạng hàm phân bố tải trọng (đều, tam giác, hình sin) đến phân bố trường ứng suất của vỏ.

Chương 4. NGHIÊN CỨU KHẢO SÁT VỎ TRỤ FG-CNTRC CHỊU TẢI TRỌNG CƠ NHIỆT

Chương 4 mở rộng khảo sát bài toán cơ nhiệt kết hợp:

  • Xác định phân bố nhiệt độ: Thiết lập và giải phương trình truyền nhiệt dừng một chiều theo chiều dày vỏ $\frac{d}{dz}\left(k_{33}(z)\frac{dT}{dz}\right) = 0$ với điều kiện biên nhiệt độ mặt trong $T_{in}$ và mặt ngoài $T_{out}$, thu được phân bố nhiệt độ phi tuyến chính xác qua chiều dày.
  • Kiểm chứng mô hình: Đối chiếu kết quả phân bố nhiệt độ, chuyển vị $w$ và ứng suất $\sigma_\theta$ với kết quả của các công trình công bố quốc tế về vỏ FGM và vỏ FG-CNTRC gia cường CNT lượn sóng trong môi trường nhiệt.
  • Khảo sát đáp ứng tĩnh cơ nhiệt:
    • Khảo sát các trường nhiệt độ khác nhau ($T_{in} = 300\text{ K}, T_{out} = 400\text{ K}$ và ngược lại) kết hợp tải áp suất mặt trong/mặt ngoài.
    • Ảnh hưởng của kiểu phân bố CNT: Phân bố dạng FG-Ʌ và FG-V thể hiện sự phân hóa ứng suất nhiệt rất rõ giữa mặt trong và mặt ngoài do độ dẫn nhiệt $k_{33}(z)$ và hệ số giãn nở nhiệt biến thiên bất đối xứng theo chiều dày.
    • Ảnh hưởng của tỷ lệ thể tích $V^*_{CNT}$: Tăng tỷ lệ thể tích CNT giúp hạn chế chuyển vị nhiệt và giảm đỉnh ứng suất kéo do nhiệt gây ra.
    • Đánh giá sự tương tác giữa tải trọng nhiệt và tải trọng cơ: Tải trọng nhiệt đóng góp thành phần ứng suất nén/kéo ban đầu làm thay đổi đáng kể tổng trạng thái ứng suất tĩnh của vỏ trụ.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận và khoa học

  1. Thiết lập thành công mô hình toán học giải tích cho bài toán phân tích tĩnh vỏ trụ FG-CNTRC chịu tải trọng cơ nhiệt dựa trên lý thuyết biến dạng cắt bậc cao kiểu quasi-3D có kể đến ứng suất pháp tuyến ngang ($\sigma_z \neq 0$), kết hợp đồng thời với gradient nhiệt độ thực tế xác định từ phương trình truyền nhiệt và tính chất cơ lý của vật liệu phụ thuộc vào nhiệt độ.
  2. Xây dựng giải pháp giải tích kết hợp chuỗi Fourier lượng giác và phép biến đổi Laplace, cho phép giải hệ phương trình vi phân đạo hàm riêng bậc cao của vỏ trụ FG-CNTRC với các điều kiện biên phức tạp khác nhau (C-C, C-F, C-S, S-S) chịu tải trọng cục bộ phân bố không đều.
  3. Ứng dụng các phương trình cân bằng vi phân của lý thuyết đàn hồi 3 chiều để xác định chính xác các thành phần ứng suất cắt ngang ($\sigma_{\xi z}, \sigma_{\theta z}$), đảm bảo thỏa mãn điều kiện biên không ứng suất tiếp tại các bề mặt tự do của vỏ.

Đóng góp mới về mặt thực tiễn và kiến nghị ứng dụng

  1. Đưa ra bức tranh định lượng chi tiết về hiện tượng tập trung ứng suất và hiệu ứng biên tại vùng chuyển tiếp liên kết của vỏ trụ FG-CNTRC, giúp chỉ ra vị trí nguy hiểm nhất về mặt độ bền của kết cấu.
  2. Đánh giá định lượng hiệu quả gia cường của 5 quy luật phân bố CNT (UD, FG-Ʌ, FG-V, FG-O, FG-X) và các mức tỷ lệ thể tích $V^*_{CNT}$, chứng minh rằng cấu hình phân bố CNT tập trung nhiều ở bề mặt chịu tải (dạng FG-V hoặc FG-X) mang lại khả năng chịu lực tối ưu và giảm ứng suất cục bộ tốt nhất.
  3. Cung cấp bộ chương trình tính toán số có độ chính xác cao làm công cụ tính toán, kiểm tra bền và hỗ trợ tối ưu hóa thiết kế kết cấu dạng vỏ trụ chịu tải cơ nhiệt làm việc trong môi trường hàng không vũ trụ và công nghiệp chịu nhiệt.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài

  • Nghiên cứu giới hạn trong phạm vi ứng xử đàn hồi tuyến tính của vật liệu, chưa xét đến tính chất phi tuyến vật liệu hoặc phi tuyến hình học khi kết cấu chịu tải trọng lớn gây biến dạng lớn.
  • Mô hình giả định các ống nano carbon (CNT) thẳng tuyệt đối và định hướng hoàn toàn đồng phương theo trục dọc của vỏ trụ, chưa xét đến yếu tố uốn lượn (waviness) ngẫu nhiên hoặc sự kết tụ (agglomeration) của CNT trong nền polymer.
  • Luận án tập trung vào phương pháp nghiên cứu lý thuyết và tính toán số, chưa thực hiện các khảo sát thực nghiệm chế tạo và đo đạc trực tiếp trên mẫu vỏ trụ FG-CNTRC thực tế.

Hướng nghiên cứu tiếp theo

  • Mở rộng mô hình quasi-3D để phân tích bài toán ổn định (mất ổn định tĩnh và sau mất ổn định) cũng như bài toán dao động tự do, dao động cưỡng bức của vỏ trụ FG-CNTRC có xét đến ứng suất pháp tuyến ngang.
  • Nghiên cứu kết cấu vỏ FG-CNTRC có kể đến sự phân tán không hoàn hảo của CNT (hiện tượng kết tụ, hướng phân bố ngẫu nhiên 2D/3D, hình học CNT lượn sóng).
  • Nghiên cứu bài toán tối ưu hóa cấu trúc phân bố CNT và bài toán điều khiển dao động chủ động khi kết hợp thêm các lớp vật liệu thông minh (áp điện, từ giảo).
  • Tiến hành các nghiên cứu thực nghiệm chế tạo và kiểm chứng cơ tính thực tế của kết cấu vỏ composite nano biến thiên.

Giá trị tham khảo

Đối tượng thụ hưởng và ứng dụng

  • Nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học vật rắn, Kỹ thuật Hàng không vũ trụ: Tài liệu tham khảo hữu ích về phương pháp giải tích giải hệ phương trình vi phân đạo hàm riêng bậc cao bằng biến đổi Laplace và khai triển chuỗi Fourier.
  • Giảng viên và nhà nghiên cứu: Cung cấp hệ thống tài liệu tổng quan toàn diện về vật liệu FG-CNTRC, lý thuyết biến dạng cắt quasi-3D và các mô hình cơ nhiệt mở rộng.
  • Kỹ sư thiết kế kết cấu thiết bị bay, ống áp lực và vỏ động cơ: Cung cấp cơ sở dữ liệu số và các quy luật phân bố ứng suất, định hướng lựa chọn kiểu phân bố và tỷ lệ thể tích CNT phù hợp nhằm tăng khả năng chịu lực tại các vùng liên kết biên nguy hiểm.

Các nội dung có giá trị tham khảo cao nhất

  • Chương 2: Toàn bộ quy trình toán học thiết lập hệ 11 phương trình cân bằng từ nguyên lý công ảo và kỹ thuật áp dụng phép biến đổi Laplace tìm nghiệm giải tích cho các điều kiện biên khác nhau.
  • Chương 3 và Chương 4: Hệ thống các bảng số liệu kiểm chứng, các đồ thị phân bố ứng suất theo chiều dày và dọc theo đường sinh vỏ, đặc biệt là các quy luật định lượng về hiệu ứng biên và tương tác cơ nhiệt.

Câu hỏi thường gặp

1. Tại sao luận án cần sử dụng lý thuyết biến dạng cắt bậc cao quasi-3D (mô hình K3-2) thay vì lý thuyết cổ điển hay FSDT?

Lý thuyết vỏ cổ điển (CST) hoàn toàn bỏ qua biến dạng trượt ngang, trong khi lý thuyết biến dạng cắt bậc nhất (FSDT) đòi hỏi hệ số hiệu chỉnh cắt và giả thiết chuyển vị pháp tuyến là hằng số theo chiều dày ($\varepsilon_z = 0, \sigma_z = 0$). Đối với vỏ có chiều dày từ trung bình đến dày hoặc làm từ vật liệu composite nanocomposite có cơ tính biến thiên, ứng suất pháp ngang $\sigma_z$ đóng vai trò quan trọng. Mô hình quasi-3D K3-2 khai triển chuyển vị $w$ bậc hai theo $z$, cho phép xác định chính xác trạng thái ứng suất 3 chiều và mô tả đầy đủ hiện tượng tập trung ứng suất tại vùng biên mà các lý thuyết thông thường bỏ qua.

2. Vật liệu FG-CNTRC trong luận án gồm những thành phần nào và có bao nhiêu quy luật phân bố CNT?

Vật liệu FG-CNTRC trong luận án được cấu tạo từ nền polymer (như epoxy) gia cường bằng các ống nano carbon đơn vách SWCNT (10,10) phân bố theo trục dọc vỏ. Luận án khảo sát 5 quy luật phân bố tỷ lệ thể tích CNT theo chiều dày vỏ ($z \in [-h/2, h/2]$):

  • UD: Phân bố đều ($V_{CNT} = V^*_{CNT}$).
  • FG-Ʌ: Phân bố giảm tuyến tính từ mặt trong ra mặt ngoài ($V_{CNT} = 2(1/2 - z/h)V^*_{CNT}$).
  • FG-V: Phân bố tăng tuyến tính từ mặt trong ra mặt ngoài ($V_{CNT} = 2(1/2 + z/h)V^*_{CNT}$).
  • FG-O: Phân bố tập trung ở giữa và giảm dần về hai bề mặt ($V_{CNT} = 2(1 - 2|z|/h)V^*_{CNT}$).
  • FG-X: Phân bố tập trung nhiều ở hai bề mặt và ít nhất ở mặt giữa ($V_{CNT} = 4(|z|/h)V^*_{CNT}$).

3. Phương pháp giải tích của luận án xử lý các điều kiện biên ngàm (C) và tự do (F) như thế nào?

Để giải quyết bài toán với các điều kiện biên khác nhau, tác giả sử dụng chuỗi lượng giác đơn theo hướng vòng $\theta$ để chuyển hệ phương trình đạo hàm riêng (PDE) thành hệ phương trình vi phân thường (ODE) theo tọa độ dọc trục $\xi$. Sau đó, áp dụng phép biến đổi Laplace đối với biến $\xi$ để chuyển hệ ODE thành hệ phương trình đại số trong không gian Laplace. Các giá trị biên tại $\xi = 0$ và $\xi = L/R$ (như chuyển vị triệt tiêu ở biên ngàm C, hoặc nội lực triệt tiêu ở biên tự do F) được sử dụng để xác định chính xác các hằng số tích phân trước khi thực hiện biến đổi Laplace ngược về miền thực.

4. Tải trọng nhiệt độ trong luận án được mô hình hóa như thế nào?

Khác với các nghiên cứu giả định nhiệt độ là hằng số hoặc phân bố tuyến tính đơn giản, luận án giải phương trình truyền nhiệt dừng một chiều theo phương chiều dày vỏ: $$\frac{d}{dz}\left(k_{33}(z)\frac{dT}{dz}\right) = 0$$ với điều kiện nhiệt độ bề mặt trong $T_{in}$ và mặt ngoài $T_{out}$. Do hệ số dẫn nhiệt $k_{33}(z)$ thay đổi liên tục theo quy luật phân bố CNT, trường nhiệt độ $T(z)$ thu được là một hàm phi tuyến theo chiều dày. Đồng thời, toàn bộ các mô-đun đàn hồi $E_{11}, E_{22}, G_{12}$ và hệ số giãn nở nhiệt $\alpha_{11}, \alpha_{22}$ của vật liệu đều được cập nhật phụ thuộc vào giá trị nhiệt độ cục bộ tại từng điểm tính toán.

5. Kết quả khảo sát cho thấy kiểu phân bố CNT nào có khả năng chịu lực và giảm ứng suất tốt nhất?

Kết quả tính toán ở Chương 3 và Chương 4 chỉ ra rằng các kiểu phân bố CNT dạng FG-V (khi tải tác dụng ở mặt ngoài) hoặc FG-X (tập trung CNT ở cả hai bề mặt ngoài) giúp tăng cường độ cứng uốn cục bộ của vỏ trụ đáng kể nhất. Các cấu hình này làm giảm biên độ chuyển vị hướng kính $w$ và giảm nồng độ ứng suất tập trung tại vùng liên kết biên ngàm so với kiểu phân bố đều (UD) và kiểu phân bố tập trung ở giữa (FG-O).

Kết luận

Luận án tiến sĩ của tác giả Dương Văn Quang đã giải quyết trọn vẹn bài toán phân tích tĩnh kết cấu vỏ trụ nanocomposite cơ tính biến thiên gia cường ống nano carbon (FG-CNTRC) chịu tác dụng đồng thời của tải trọng cơ học và nhiệt độ. Bằng việc xây dựng thành công mô hình lý thuyết biến dạng cắt bậc cao quasi-3D có xét đến ứng suất pháp tuyến ngang và áp dụng phương pháp giải tích biến đổi Laplace kết hợp chuỗi Fourier, công trình đã làm sáng tỏ quy luật phân bố ứng suất và hiện tượng hiệu ứng biên dưới các điều kiện biên phức tạp (C-C, C-F, C-S, S-S). Các kết quả tính toán và phát hiện khoa học trong luận án là cơ sở dữ liệu tin cậy phục vụ công tác tính toán, thiết kế và tối ưu hóa kết cấu vỏ composite tiên tiến trong ngành kỹ thuật hàng không vũ trụ và các ngành công nghiệp kỹ thuật cao.