Tổng quan về luận án

Nghiên cứu của tác giả Trần Đại Hào (2024) thuộc chuyên ngành Cơ kỹ thuật (mã số: 9520101), Trường Đại học Xây dựng Hà Nội, dưới sự hướng dẫn khoa học của TS. Đặng Xuân Hùng và GS.TS. Trần Minh Tú, mang tên "Phân tích tĩnh, dao động riêng và đánh giá độ tin cậy của tấm composite lớp", đại diện cho một bước tiến tiên phong trong việc giải quyết đồng thời bài toán cơ học kết cấu tất định và bài toán đánh giá độ an toàn ngẫu nhiên của các kết cấu composite tiên tiến.

Trong bối cảnh khoa học kỹ thuật đương đại, vật liệu composite cốt sợi đồng phương truyền thống (như sợi carbon TR50S, sợi thủy tinh E-glass, sợi aramid Kevlar 129) và đặc biệt là vật liệu nano composite có cơ tính biến thiên gia cường bằng ống nano carbon đơn vách (Functionally Graded Carbon Nanotube Reinforced Composite - FG-CNTRC) đang trở thành tâm điểm của ngành hàng không vũ trụ, công trình biển, giao thông vận tải và năng lượng tái tạo. Tuy nhiên, khoảng trống nghiên cứu (research gap) cốt lõi vẫn tồn tại sâu sắc trên hai phương diện:

  1. Về mô hình hóa cơ học: Lý thuyết biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) tuy có cấu trúc toán học gọn gàng và hiệu quả tính toán cao hơn lý thuyết biến dạng cắt bậc cao (Higher-order Shear Deformation Theory - HSDT) hay lý thuyết nhiều lớp (Layerwise Theory), nhưng lại mắc phải nhược điểm cố hữu là giả thiết ứng suất cắt ngang không đổi theo chiều dày, đòi hỏi hệ số hiệu chỉnh cắt $K_s$ phức tạp và không thỏa mãn điều kiện biên ứng suất cắt bằng không tại bề mặt tự do. Mặt khác, các phương pháp giải tích cổ điển như phương pháp Navier hay phương pháp Levy chỉ áp dụng được cho các điều kiện biên lý tưởng rất hạn chế (như bốn biên tựa khớp SSSS hoặc hai biên đối diện tựa khớp).
  2. Về độ tin cậy và an toàn kết cấu: Các phương pháp thiết kế truyền thống chủ yếu dựa trên mô hình tất định (deterministic model), xem các tham số đầu vào như mô đun đàn hồi, kích thước hình học, góc phương sợi và tải trọng là những hằng số cố định. Đúng như tác giả đã khẳng định trong văn bản luận án: "trong thực tế, mặc dù các kết cấu đã được thiết kế theo các phương pháp tất định nhưng vẫn có xác suất không an toàn do chưa kể đến xác suất phá hỏng các cấu kiện, sự an toàn tổng thể của kết cấu công trình, các yếu tố không chắc chắn của các tham số thiết kế".

Để giải quyết triệt để các khoảng trống nêu trên, luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một mô hình bán giải tích vừa đảm bảo tính linh hoạt với mọi tổ hợp điều kiện biên phức tạp (SSSS, CCCC, CSCS, CFCF, SFSF), vừa khắc phục được nhược điểm ứng suất cắt ngang của lý thuyết FSDT cho tấm composite nhiều lớp và tấm FG-CNTRC?
  • Giả thuyết 1 (H1): Việc kết hợp phương pháp năng lượng Pb2-Ritz với hệ phương trình cân bằng vi phân của lý thuyết đàn hồi 3D sẽ tái lập chính xác trường ứng suất cắt ngang dạng parabol liên tục qua chiều dày mà không cần nâng bậc trường chuyển vị lên HSDT.
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố ống nano carbon (UD, FG-V, FG-O, FG-X), tỷ phần thể tích $V^*_{CNT}$, góc phương sợi và số lớp ảnh hưởng định lượng như thế nào đến độ võng và tần số dao động riêng của tấm composite lớp?
  • Giả thuyết 2 (H2): Dạng phân bố FG-X tập trung CNT ở hai bề mặt ngoài sẽ tối ưu hóa độ cứng chống uốn và nâng cao tần số dao động riêng cơ bản vượt trội so với phân bố đều (UD) hoặc FG-O.
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ biến động ngẫu nhiên của từng nhóm tham số vật liệu, hình học và tải trọng đóng góp vào xác suất phá hủy $P_f$ của tấm composite lớp theo tiêu chuẩn bền Tsai-Wu như thế nào?
  • Giả thuyết 3 (H3): Độ tin cậy của kết cấu nhạy cảm nhất với sự không chắc chắn của tải trọng tác dụng và mô đun đàn hồi, trong khi góc phương sợi tạo ra hiệu ứng phi tuyến mạnh lên xác suất không an toàn.

Khung lý thuyết tổng quát tích hợp lý thuyết biến dạng cắt FSDT, nguyên lý biến phân cực tiểu thế năng toàn phần Pb2-Ritz, lý thuyết đồng nhất hóa vật liệu vi mô (Eshelby-Mori-Tanaka và Rule of Mixture mở rộng), cùng phương pháp mô phỏng ngẫu nhiên Monte Carlo kết hợp thuyết bền phá hủy Tsai-Wu. Phạm vi nghiên cứu bao trùm các tấm composite chữ nhật với tỷ số kích thước $a/h \in [10, 50]$, tỷ số cạnh $b/a \in [1, 2]$, tỷ phần thể tích CNT từ $V^*_{CNT} = 0.11$ đến $0.17$, với quy mô mẫu mô phỏng Monte Carlo đạt $10^5$ đến $10^6$ lần lặp, đem lại độ chuẩn xác và tin cậy cao.

Literature Review và Positioning

Tổng quan y văn trong luận án hệ thống hóa các luồng nghiên cứu chủ đạo trong cơ học kết cấu composite qua hơn một thế kỷ:

  1. Lý thuyết tấm đơn lớp tương đương (ESL): Bắt đầu từ lý thuyết tấm cổ điển CPT của Kirchhoff (1850) và Love (1888), được áp dụng rộng rãi bởi Noor (1973) [97], Reddy và Khdeir (1989) [113]. Tuy nhiên, CPT bỏ qua hoàn toàn biến dạng trượt ngang nên chỉ chính xác với tấm rất mỏng. Nhằm khắc phục hạn chế này, Reissner (1945) và Mindlin (1951) [87] phát triển FSDT, được kế thừa bởi Whitney (1973) [153], Liew et al. (1995) [68], Zhu et al. (2012) [178], và Lei et al. (2013) [61]. Để loại bỏ sự phụ thuộc vào hệ số hiệu chỉnh cắt $K_s$, Reddy (1984) và Phan & Reddy (1985) [102] đề xuất lý thuyết biến dạng cắt bậc ba 5 ẩn (HSDT-5), tiếp tục được Shen & Zhang (2010) [130], Shen & Zhu (2012) [129], Selim et al. (2016) [124] ứng dụng cho tấm nano composite.
  2. Vật liệu nano composite FG-CNTRC: Sau phát hiện lịch sử của Sumio Iijima (1991, 1993) [48] về ống nano carbon đơn vách (SWCNT) và đa vách (MWCNT) với mô đun đàn hồi vượt ngưỡng 1 TPa và độ bền kéo đạt ~100 GPa, Shen (2009) lần đầu tiên đề xuất khái niệm FG-CNTRC. Hàng loạt nghiên cứu của Kiani (2016) [89], Zhang et al. (2014) [167], Wang & Shen (2012) [151] đã phân tích ứng xử phi tuyến nhiệt và dao động của loại vật liệu này.

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  • Tranh luận 1: ESL vs. Layerwise Theory: Các tác giả Robbins & Reddy (1993) [120] và Malekzadeh & Heydarpour (2014) [82] chỉ ra rằng ESL không mô tả chính xác ứng suất tại bề mặt tiếp giáp giữa các lớp có độ cứng lệch nhau quá lớn (như tấm sandwich). Ngược lại, những người ủng hộ ESL lập luận rằng Layerwise làm bùng nổ số ẩn số theo số lớp, khiến việc phân tích kết cấu thực tế trở nên bất khả thi về mặt tài nguyên tính toán.
  • Tranh luận 2: Phương pháp số (FEM/Meshless) vs. Phương pháp bán giải tích (Ritz variants): Trong khi phương pháp phần tử hữu hạn (FEM) linh hoạt về hình học nhưng đòi hỏi chia lưới phức tạp và dễ gặp hiện tượng "khóa cắt" (shear locking), các phương pháp bán giải tích như DQM-Ritz (Eftekhari & Jafari 2012 [30]), Chebyshev-Ritz (Kiani et al. 2013 [53]), Gram-Schmidt Ritz (Hu et al. 2004 [46]) cho độ chính xác cao nhưng gặp khó khăn khi mô tả biên tự do (Free edge).

Luận án định vị chính xác khoảng trống bằng cách sử dụng phương pháp Pb2-Ritz với chuỗi đa thức hai chiều đầy đủ kết hợp hàm biên toán học. Định vị này tạo ra sự vượt trội rõ rệt khi so sánh với ít nhất hai nghiên cứu quốc tế:

  • So sánh với nghiên cứu của Ramtekkar (2002) [111] về tải trọng phá hủy lớp đầu tiên ($q_{th}$) theo thuyết bền Tsai-Wu: mô hình của luận án đạt độ khớp nối với sai số dưới 1.18% nhưng giảm thiểu thời gian tính toán nhờ tránh được việc phân chia phần tử 3D cồng kềnh.
  • So sánh với công bố của Liew et al. (1995) [68] và Shen & Zhang (2010) [130]: mô hình Pb2-Ritz trong luận án xử lý mượt mà các tổ hợp biên hỗn hợp không đối xứng (như CSCS, CFCF, SFSF) mà các phương pháp Navier hay Levy thông thường hoàn toàn bất lực.

Đó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 Reissner-Mindlin thông qua việc thiết lập thuật toán khôi phục ứng suất cắt ngang dựa trên hệ phương trình cân bằng vi phân của lý thuyết đàn hồi ba chiều: $$\frac{\partial \sigma_{xx}}{\partial x} + \frac{\partial \tau_{xy}}{\partial y} + \frac{\partial \tau_{xz}}{\partial z} = 0$$ $$\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_{yy}}{\partial y} + \frac{\partial \tau_{yz}}{\partial z} = 0$$

Bằng cách tích phân ứng suất pháp $\sigma_{xx}, \sigma_{yy}$ và ứng suất tiếp $\tau_{xy}$ thu được từ lời giải Pb2-Ritz theo tọa độ chiều dày $z$, trường ứng suất cắt ngang $\tau_{xz}(z)$ và $\tau_{yz}(z)$ được tái tạo lại dưới dạng hàm bậc hai (parabol), tự động thỏa mãn điều kiện biên bằng không tại hai mặt tự do trên ($z = h/2$) và dưới ($z = -h/2$) của tấm: $$\tau_{xz}(z) = -\int_{-h/2}^{z} \left(\frac{\partial \sigma_{xx}}{\partial x} + \frac{\partial \tau_{xy}}{\partial y}\right) dz$$ $$\tau_{yz}(z) = -\int_{-h/2}^{z} \left(\frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_{yy}}{\partial y}\right) dz$$

Đóng góp này giải quyết mâu thuẫn tồn tại hàng thập kỷ giữa tính đơn giản của FSDT và độ chính xác của ứng suất cắt phân lớp mà không phải trả giá bằng việc tăng số biến chuyển vị lên 5 hoặc 12 ẩn như trong HSDT.

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 ba trụ cột lý thuyết:

  1. Lý thuyết đồng nhất hóa mở rộng: Sử dụng mô hình quy tắc hỗn hợp mở rộng bổ sung hệ số ảnh hưởng CNT ($\eta_j$) phụ thuộc vào tỷ phần thể tích $V^*{CNT}$ để tính toán chính xác các mô đun $E{11}, E_{22}, G_{12}$: $$E_{11} = \eta_1 V_{CNT} E_{11}^{CNT} + V_m E_m$$ $$\frac{\eta_2}{E_{22}} = \frac{V_{CNT}}{E_{22}^{CNT}} + \frac{V_m}{E_m}, \quad \frac{\eta_3}{G_{12}} = \frac{V_{CNT}}{G_{12}^{CNT}} + \frac{V_m}{G_m}$$
  2. Phương pháp năng lượng biến phân Pb2-Ritz: Trường chuyển vị $u_0, v_0, w_0, \psi_x, \psi_y$ được biểu diễn bằng tích của hàm biên $\phi(x,y)$ và chuỗi đa thức trực giao 2 chiều đầy đủ bậc $M \times N$: $$u_0(x,y) = \phi_u(x,y) \sum_{m=1}^M \sum_{n=1}^N c_{mn}^u x^{m-1} y^{n-1}$$ Hàm biên $\phi(x,y) = \prod_{i=1}^4 [\Gamma_i(x,y)]^{k_i}$ tự động thỏa mãn các điều kiện biên hình học nghiêm ngặt với $k_i = 1$ cho biên ngàm (C) hoặc khớp (S), và $k_i = 0$ cho biên tự do (F).
  3. Mô hình ngẫu nhiên hóa Monte Carlo phi tuyến: Tích hợp điều kiện an toàn kết cấu thông qua hàm trạng thái giới hạn $g(\mathbf{X}) = R(\mathbf{X}) - S(\mathbf{X})$ với tiêu chuẩn phá hủy tương tác đa trục Tsai-Wu: $$F_1 \sigma_1 + F_2 \sigma_2 + F_{11} \sigma_1^2 + F_{22} \sigma_2^2 + F_{66} \tau_{12}^2 + 2 F_{12} \sigma_1 \sigma_2 = 1$$ trong đó các hệ số $F_i, F_{ij}$ được xác định từ các giới hạn bền kéo/nén dọc ($X_t, X_c$), ngang ($Y_t, Y_c$) và cắt phẳng ($S_{12}$).

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

Thiết kế nghiên cứu

Nghiên cứu đứng vững trên lập trường triết học thực chứng (Positivism) và chủ nghĩa hiện thực khoa học (Scientific Realism), khẳng định các quy luật cơ học của vật liệu và kết cấu là khách quan, có thể đo lường, mô hình hóa toán học và kiểm chứng thực nghiệm số độc lập.

Thiết kế nghiên cứu đa cấp độ (Multi-level design) được triển khai qua 4 nấc thang:

  • Cấp độ Nano/Micro: Mô phỏng tương tác giữa ống nano carbon SWCNT và vật liệu nền polymer Poly[(m-phenylenevinylene)-co-(2,5-dioctoxy-p-phenylenevinylene)] (PmPV).
  • Cấp độ Trung mô (Meso-scale): Từng đơn lớp composite trực hướng (orthotropic lamina) với góc xoay sợi $\theta_k$ và quy luật phân bố cơ tính theo chiều dày (UD, FG-V, FG-O, FG-X).
  • Cấp độ Vĩ mô (Macro-scale): Tấm composite nhiều lớp tổng thể chịu tải uốn và dao động tự do với ma trận độ cứng liên kết $[A, B, D, A_s]$.
  • Cấp độ Ngẫu nhiên (Stochastic level): Mô phỏng không gian tham số phân tán của tải trọng và vật liệu để đánh giá độ tin cậy.

Quy trình nghiên cứu rigorous

Quy trình giải bài toán tĩnh và dao động riêng bằng Pb2-Ritz:

  1. Thiết lập phiếm hàm năng lượng: Thế năng toàn phần $\Pi = U - V$ với $U = \frac{1}{2} \int_A (\boldsymbol{\epsilon}^0)^T \mathbf{N} + (\boldsymbol{\kappa})^T \mathbf{M} + (\boldsymbol{\gamma})^T \mathbf{Q} , dA$, công ngoại lực $V = \int_A q(x,y) w_0 , dA$.
  2. Cực tiểu hóa biến phân: Áp dụng điều kiện cân bằng Euler-Lagrange $\frac{\partial \Pi}{\partial c_i} = 0$, đưa về hệ phương trình đại số tuyến tính $[K]{C} = {F}$.
  3. Bài toán trị riêng dao động: Áp dụng nguyên lý Hamilton $\delta \int_{t_1}^{t_2} (K - U) dt = 0$ với động năng $K$, dẫn đến phương trình dao động tự do: $$([K] - \omega^2 [M]) {C} = 0$$ Nghiệm của định thức $|[K] - \omega^2 [M]| = 0$ cung cấp các tần số dao động riêng $\omega_i$ và dạng dao động (mode shapes) tương ứng.

Data và phân tích

Toàn bộ hệ thống thuật toán được hiện thực hóa bằng mã nguồn tự phát triển trên nền tảng MATLAB.

Bảng 1: Thông số cơ lý của vật liệu nền và vật liệu gia cường trong mô hình

Vật liệu / Cấu phần Mô đun Young ($E$) Hệ số Poisson ($\nu$) Khối lượng riêng ($\rho$) Độ bền kéo/nén ($\sigma_u$)
Nền Polymer (PmPV) $E_m = 2.1\text{ GPa}$ $\nu_m = 0.34$ $\rho_m = 1150\text{ kg/m}^3$ $\sigma_m = 80\text{ MPa}$
SWCNT (10,10) Armchair $E_{11}^{CNT} = 5.646\text{ TPa}$
$E_{22}^{CNT} = 7.080\text{ TPa}$		 $\nu_{12}^{CNT} = 0.175$ $\rho_{CNT} = 1400\text{ kg/m}^3$ $\sigma_{CNT} \approx 100\text{ GPa}$
$G_{12}^{CNT} = 1.944\text{ TPa}$
Sợi Carbon TR50S $E = 235\text{ GPa}$ $\nu = 0.20$ $\rho = 1820\text{ kg/m}^3$ Cường độ kéo: $4.90\text{ GPa}$
Sợi Aramid Kevlar 129 $E = 100\text{ GPa}$ $\nu = 0.35$ $\rho = 1440\text{ kg/m}^3$ Cường độ kéo: $3.38\text{ GPa}$

Khảo sát sự hội tụ số (Convergence test) chứng minh rằng khi bậc đa thức của hàm Pb2-Ritz tăng từ $P = 4$ đến $P = 10$ ($M = N = 10$), kết quả độ võng không thứ nguyên $\bar{w}$ và ba tần số dao động đầu tiên $\bar{\omega}_1, \bar{\omega}_2, \bar{\omega}_3$ đạt mức ổn định tuyệt đối với sai số dao động nhỏ hơn $0.05%$.

Trong mô phỏng Monte Carlo, các biến ngẫu nhiên gồm:

  • $E_1, E_2, G_{12}, \nu_{12}$: phân phối chuẩn (Gaussian) hoặc phân phối log-normal với hệ số biến động $CoV \in [0.05, 0.20]$.
  • Góc phương sợi $\theta_k$: phân phối chuẩn quanh giá trị danh định ($0^\circ, 45^\circ, 90^\circ$) với độ lệch chuẩn $\sigma_\theta \in [1^\circ, 5^\circ]$.
  • Chiều dày lớp $h_k$: $CoV \in [0.02, 0.10]$.
  • Tải trọng ngoài $q_0$: phân phối chuẩn với $CoV_q \in [0.10, 0.30]$. Số lần mô phỏng $N_{MC} = 10^5 \div 10^6$ bảo đảm hệ số biến thiên của xác suất phá hủy $\text{CoV}(P_f) = \sqrt{\frac{1-P_f}{N_{MC} \cdot P_f}} < 5%$.

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

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

  1. Ưu thế tuyệt đối của quy luật phân bố FG-X: Trong 4 quy luật phân bố ống nano carbon qua chiều dày (UD, FG-V, FG-O, FG-X), tấm nano composite có cấu hình FG-X luôn đạt độ võng nhỏ nhất và tần số dao động riêng cơ bản cao nhất ở mọi điều kiện biên. Cụ thể, với tỷ số $a/h = 10$ và $V^*{CNT} = 0.17$, tấm cấu hình FG-X giảm độ võng cực đại tới 34.2% so với phân bố đều (UD) và 48.7% so with phân bố FG-O. Điều này được giải thích về mặt cơ học: FG-X tập trung mật độ sợi CNT tối đa tại hai bề mặt ngoài ($z = \pm h/2$), nơi ứng suất uốn pháp $\sigma{xx}, \sigma_{yy}$ đạt giá trị lớn nhất, tối ưu hóa mômen quán tính kháng uốn hiệu dụng của tấm.
  2. Hiệu ứng liên kết góc phương sợi và số lớp: Đối với tấm composite nhiều lớp $[0^\circ/90^\circ]n$ và $[45^\circ/-45^\circ]n$, khi tăng số lớp từ $n = 2$ lên $n = 8$ ở cùng tổng chiều dày $h$, độ võng không thứ nguyên giảm mạnh và tiệm cận về giá trị ổn định. Hiện tượng này chứng minh rằng hiệu ứng liên kết uốn-kéo ($B{ij}$) và uốn-xoắn ($D{16}, D_{26}$) bị triệt tiêu nhanh chóng khi số lớp tăng, làm tăng đáng kể độ cứng chống uốn tổng thể.
  3. Phục hồi chính xác ứng suất cắt ngang: Phân bố ứng suất tiếp $\bar{\tau}{xz}, \bar{\tau}{yz}$ tính toán từ phương trình cân bằng 3D thể hiện dạng parabol trơn tru, liên tục qua các mặt phân giới lớp và triệt tiêu chính xác tại các mặt biên ngoài. Ngược lại, ứng suất cắt FSDT thuần túy tạo ra các bước nhảy gián đoạn phi thực tế giữa các lớp liền kề.
  4. Độ nhạy của độ tin cậy kết cấu ($P_f$): Kết quả mô phỏng Monte Carlo chỉ ra rằng xác suất không an toàn $P_f$ nhạy cảm phi tuyến mạnh nhất với hệ số biến động của tải trọng ($CoV_q$) và mô đun đàn hồi $E_1$. Khi $CoV_q$ tăng từ 0.10 lên 0.25, xác suất phá hủy $P_f$ tăng vọt từ $10^{-4}$ lên $4.8 \times 10^{-2}$ (tăng gần 500 lần). Biến động góc phương sợi $\sigma_\theta > 3^\circ$ làm suy giảm độ tin cậy nghiêm trọng ở các cấu hình góc xiên $[45^\circ/-45^\circ]$ do làm lệch hướng trục chịu lực chính.
  5. Kiểm chứng chính xác cao: So sánh kết quả tải trọng phá hủy lớp đầu tiên $q_{th}$ (kPa) với công bố kinh điển của Ramtekkar [111] cho thấy sự trùng khớp tuyệt vời:

Bảng 2: So sánh tải trọng phá hủy lớp đầu tiên $q_{th}$ (kPa) theo thuyết bền Tsai-Wu

Cấu hình tấm Tỷ số $a/h$ Ramtekkar [111] Luận án (Trần Đại Hào) Độ lệch tương đối (%)
$[0^\circ/90^\circ]$ 10 185.42 187.15 +0.93%
$[0^\circ/90^\circ/0^\circ]$ 10 242.10 244.80 +1.11%
$[45^\circ/-45^\circ]$ 10 156.80 158.20 +0.89%

Implications đa chiều

  • Về mặt học thuật: Thiết lập một khung tính toán chuẩn hóa kết hợp phương pháp biến phân Pb2-Ritz và lý thuyết độ tin cậy ngẫu nhiên, mở ra công cụ phân tích hiệu năng cao cho các cấu kiện composite thành mỏng.
  • Về mặt kỹ thuật công nghệ: Cung cấp nguyên lý thiết kế tối ưu hóa chức năng (functionally graded design) cho các kỹ sư: ưu tiên bố trí vật liệu nano theo quy luật FG-X thay vì trộn đều đồng nhất (UD), giúp tiết kiệm tới 30-40% lượng hạt nano CNT đắt tiền mà vẫn đảm bảo khả năng chịu lực tương đương.
  • Về quy chuẩn thiết kế: Khuyến nghị các cơ quan soạn thảo tiêu chuẩn (như TCVN, ASTM, Eurocode) chuyển dịch từ hệ số an toàn tĩnh đơn lẻ sang phương pháp thiết kế dựa trên độ tin cậy (Reliability-Based Design Optimization - RBDO) với hệ số an toàn khuyến nghị $n \ge 1.45$ khi hệ số biến động tải trọng vượt quá $15%$.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật nội tại:

  1. Giả thiết vật liệu composite và nano composite làm việc hoàn toàn trong miền đàn hồi tuyến tính vi mô, chưa xét đến biến dạng dẻo cục bộ của nền polymer khi chịu tải trọng cực hạn.
  2. Giả thiết liên kết giữa các lớp composite là lý tưởng tuyệt đối, chưa mô hình hóa hiện tượng tách lớp (delamination) hoặc phá hủy bề mặt tiếp giáp sợi-nền (debonding).
  3. Chưa tính đến ảnh hưởng đồng thời của môi trường ẩm nhiệt phức tạp (hygrothermal coupling effects) và tải trọng va đập động lực học tốc độ cao.
  4. Chi phí tính toán của phương pháp Monte Carlo nguyên bản (Direct MCS) với $10^6$ mẫu còn tương đối lớn, chưa tích hợp các thuật toán giảm mẫu nâng cao như Latin Hypercube Sampling (LHS) hay kỹ thuật mặt đáp ứng (Response Surface Methodology).

Các hướng nghiên cứu tiếp nối trong tương lai:

  • Mở rộng phương pháp Pb2-Ritz sang phân tích phi tuyến hình học lớn (Large deflection theo giả thiết von Kármán) và phân tích mất ổn định động học (Dynamic buckling).
  • Tích hợp mô hình học máy (Machine Learning Surrogate Models như Kriging hoặc Deep Neural Networks) để dự đoán tức thời xác suất phá hủy với thời gian thực.
  • Nghiên cứu thực nghiệm chế tạo và kiểm chứng cơ tính tấm FG-CNTRC sử dụng công nghệ in 3D tiên tiến.
  • Đánh giá độ tin cậy mỏi (Fatigue reliability) dưới tác động của tải trọng tuần hoàn ngẫu nhiên.

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

  • Tác động học thuật: Các kết quả của luận án đóng góp vào kho tàng cơ học tính toán vật liệu tiên tiến tại Việt Nam và quốc tế, với tiềm năng trích dẫn cao trong các tạp chí chuyên ngành hàng đầu nhóm Q1/Q2 (như Composite Structures, Thin-Walled Structures, Computers & Structures).
  • Chuyển đổi công nghiệp: Định hình giải pháp ứng dụng vật liệu FG-CNTRC trong chế tạo cánh quạt turbine gió công suất lớn ngoài khơi, vỏ tàu cao tốc, thân vỏ máy bay không người lái (UAV) và các bình áp lực chịu khí nén hydrogen.
  • Lợi ích kinh tế - xã hội: Giảm trọng lượng kết cấu công trình từ 20-35% so với vật liệu kim loại hoặc composite thông thường, giúp tiết kiệm nhiên liệu trong vận tải, cắt giảm phát thải carbon và nâng cao tuổi thọ công trình hạ tầng.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận một khung mẫu mực về phương pháp năng lượng Pb2-Ritz, mã nguồn MATLAB chuẩn hóa và phương pháp luận kết hợp tất định - ngẫu nhiên.
  • Giảng viên & Chuyên gia cơ học kết cấu: Nguồn tư liệu chuyên sâu giảng dạy các học phần Cao học và Tiến sĩ về Lý thuyết tấm vỏ nâng cao, Cơ học composite và Độ tin cậy công trình.
  • Kỹ sư R&D trong công nghiệp: Bộ công cụ tra cứu định lượng về ảnh hưởng của góc dệt sợi, tỷ phần thể tích nano và hệ số an toàn thực tế khi thiết kế panel composite.
  • Cơ quan quản lý & Ban soạn thảo tiêu chuẩn: Cung cấp bằng chứng khoa học thực chứng để xây dựng các chỉ dẫn kỹ thuật tính toán kết cấu composite theo phương pháp trạng thái giới hạn.

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?

Đóng góp độc đáo nhất là việc hoàn thiện mô hình FSDT thông qua việc tích hợp hệ phương trình vi phân cân bằng lý thuyết đàn hồi 3D để khôi phục trường ứng suất cắt ngang parabol thỏa mãn hoàn hảo điều kiện biên không ứng suất tại hai bề mặt ngoài. Công trình đã mở rộng lý thuyết biến dạng cắt bậc nhất của Reissner-Mindlin, loại bỏ sự phụ thuộc tiêu cực vào hệ số hiệu chỉnh cắt $K_s$ mà không làm phát sinh thêm các bậc tự do phức tạp như trong lý thuyết HSDT của Reddy.

2. Sự đổi mới về phương pháp luận (Methodology innovation) khi so sánh với các công trình trước đây?

So với phương pháp Navier (Noor 1973 [97]) vốn chỉ giải được bài toán bốn biên tựa khớp SSSS, và phương pháp phần tử hữu hạn FEM (Zhu et al. 2012 [178]) vốn đòi hỏi chi phí chia lưới lớn và dễ gặp lỗi khóa cắt, phương pháp Pb2-Ritz sử dụng tích của hàm biên $\phi(x,y)$ và chuỗi đa thức 2 chiều hoàn chỉnh cho phép giải quyết chính xác mọi tổ hợp biên bất kỳ (SSSS, CCCC, CSCS, CFCF, SFSF) với tốc độ hội tụ cực nhanh ($M = N = 10$) và độ chính xác tương đương giải tích thuần túy.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là quy luật phân bố FG-X không chỉ vượt trội về độ cứng uốn mà còn có khả năng kháng biến động tải trọng ngẫu nhiên tốt hơn nhiều so với phân bố đều (UD). Dữ liệu cho thấy ở cùng mức độ biến động tải trọng $CoV_q = 0.20$, xác suất không an toàn $P_f$ của tấm FG-X thấp hơn cấu hình UD tới 3.2 lần, khẳng định việc tối ưu hóa cơ tính theo chiều dày mang lại hiệu quả gia cường kép cả về mặt tĩnh lực học lẫn độ tin cậy ngẫu nhiên.

4. Luận án có cung cấp giao thức tái lập kết quả nghiên cứu (Replication protocol) không?

Có. Luận án cung cấp toàn bộ hệ thống công thức giải tích tường minh của các ma trận độ cứng $[K]$, ma trận khối lượng $[M]$, véc-tơ tải trọng ${F}$, cùng hệ số quy đổi hằng số đàn hồi vi mô $\eta_j$. Đồng thời, phần phụ lục đính kèm toàn bộ 5 chương trình tính toán hoàn chỉnh viết trên nền MATLAB (từ phân tích tĩnh, dao động riêng đến mô phỏng Monte Carlo với thuyết Tsai-Wu), cho phép cộng đồng khoa học độc lập kiểm chứng và tái lập 100% các bảng biểu, đồ thị.

5. Lộ trình nghiên cứu phát triển 10 năm (10-year research agenda) được vạch ra như thế nào?

Lộ trình 10 năm gồm 3 giai đoạn chiến lược:

  • Giai đoạn 1 (1-3 năm): Phát triển mô hình từ tấm phẳng sang vỏ cong composite phức hợp (cylindrical, spherical, conical shells) chịu tải ẩm nhiệt và áp điện thông minh (piezoelectric layers).
  • Giai đoạn 2 (4-6 năm): Ứng dụng trí tuệ nhân tạo (AI-driven Surrogate Modeling) để tối ưu hóa đa mục tiêu kết cấu composite có độ tin cậy cao trong thời gian thực.
  • Giai đoạn 3 (7-10 năm): Thực nghiệm hóa toàn diện, xây dựng phòng lab kiểm định độ bền mỏi - phá hủy vi cơ học và chuyển giao công nghệ chế tạo tấm sandwich FG-CNTRC cho ngành công nghiệp hàng không vũ trụ và năng lượng tái tạo.

Kết luận

Luận án tiến sĩ của tác giả Trần Đại Hào đã 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. Xây dựng thành công mô hình tất định phân tích tĩnh và dao động riêng cho tấm composite cốt sợi truyền thống và nano composite FG-CNTRC bằng phương pháp Pb2-Ritz trên nền lý thuyết FSDT.
  2. Thiết lập thuật toán khôi phục ứng suất cắt ngang chính xác từ hệ phương trình cân bằng đàn hồi 3D, khắc phục triệt để nhược điểm của FSDT mà không làm tăng khối lượng ẩn số tính toán.
  3. Đánh giá toàn diện sự ảnh hưởng của 4 quy luật phân bố CNT, chứng minh quy luật FG-X đem lại hiệu năng kháng uốn và tần số dao động cao nhất.
  4. Tích hợp thành công mô hình cơ học tất định với thuật toán ngẫu nhiên Monte Carlo và thuyết bền phá hủy tương tác Tsai-Wu để lượng hóa xác suất không an toàn $P_f$.
  5. Thực hiện phân tích độ nhạy (Sensitivity analysis), xác định rõ tải trọng và mô đun đàn hồi là các tham số chi phối lớn nhất đến độ tin cậy của tấm composite.
  6. Cung cấp bộ mã nguồn tính toán MATLAB chuẩn xác và bộ dữ liệu đối chuẩn phong phú phục vụ công tác nghiên cứu, đào tạo và thiết kế kỹ thuật thực tiễn.

Công trình tạo ra bước chuyển đổi quan trọng từ tư duy thiết kế tất định truyền thống sang phương pháp tiếp cận cơ học ngẫu nhiên hiện đại, mở ra chân trời mới cho việc ứng dụng an toàn và hiệu quả các loại vật liệu nano composite tiên tiến trong tương lai.