Tổng quan về luận án

Nghiên cứu về cơ học công trình và vật liệu hiện đại ngày càng chú trọng đến việc đánh giá độ bền và ổn định của các kết cấu tiên tiến chịu các tác động đa trường vật lý phức tạp. Công trình luận án tiến sĩ kỹ thuật của tác giả Trần Ngọc Cảnh, dưới sự hướng dẫn khoa học của PGS.TS Phạm Tiến Đạt và PGS.TS Nguyễn Văn Hưng tại Học viện Kỹ thuật Quân sự, mang tiêu đề "Tính toán vỏ thoải composite lớp chịu tác dụng đồng thời tải trọng sóng xung kích và nhiệt độ" là một nghiên cứu tiên phong trong việc mô hình hóa số và phân tích động lực học kết cấu composite nhiều lớp. Trong bối cảnh các kết cấu hàng không vũ trụ, công trình phòng thủ công sự ngầm, tàu vỏ composite cao tốc và vòm mái kiến trúc công trình biển phải hoạt động dưới điều kiện nhiệt độ cao kết hợp sóng xung kích từ vụ nổ hóa học, yêu cầu làm chủ bài toán tương tác cơ - nhiệt động lực trở thành đòi hỏi cấp thiết.

Khoảng trống nghiên cứu (research gap) trọng tâm được xác định rõ ràng: Phần lớn các nghiên cứu đương đại trên thế giới (điển hình như Reddy, 2004; Jansen, 2005; Francesco et al., 2013) mới chỉ tập trung vào vỏ mỏng hoặc tấm dày composite bám dính tuyệt đối, hoặc khảo sát riêng rẽ tải trọng nhiệt tĩnh hay tải trọng xung cơ học độc lập. Các nghiên cứu về kết cấu sandwich ba lớp có liên kết chốt trượt (shear connectors) hầu như chỉ dừng lại ở mô hình dầm 1D (như công trình của Ranzi et al., 2004; Nguyen et al., 2014) hoặc thử nghiệm tải nổ cục bộ trên dầm sandwich thép - bê tông (Richard et al., 2010). Hoàn toàn khuyết thiếu các mô hình giải tích và phần tử hữu hạn tổng quát cho kết cấu vỏ thoải composite hai độ cong liên kết chốt đàn hồi (SLCS) chịu tác dụng đồng thời của trường nhiệt độ dừng và tải trọng sóng xung kích (SXK) cân bằng xung lượng.

Luận án thiết lập và giải quyết 3 câu hỏi nghiên cứu cơ bản cùng hệ giả thuyết đồng hành:

  1. RQ1: Trường chuyển vị và biến dạng tương đối giữa các lớp vỏ composite sandwich liên kết chốt biến thiên như thế nào dưới tương tác đồng thời của tải sóng xung kích và nhiệt độ? (Giả thuyết $H_1$: Độ cứng chống cắt của chốt $k_s$ chi phối phi tuyến chuyển vị trượt $V_{tc\max}$ và làm suy giảm biên độ dao động chính của vỏ).
  2. RQ2: Mối quan hệ liên kết màng - uốn - xoắn $[A], [B], [D], [H]$ ảnh hưởng ra sao đến tần số dao động riêng và ứng suất cực đại của kết cấu panel cầu và panel trụ hai độ cong? (Giả thuyết $H_2$: Hiệu ứng nhiệt độ làm mềm ma trận độ cứng cục bộ, làm tăng độ võng $w_{\max}$ nhưng có thể điều tiết thông qua tối ưu hóa tỷ số chiều dày $h_c/h_t$ và tỷ số mô đun $E_c/E_t$).
  3. RQ3: Mô hình tải trọng sóng xung kích tam giác cân bằng xung lượng có thể mô phỏng chính xác đáp ứng động lực học của vỏ thoải so với mô hình thực nghiệm của Sadovsky không? (Giả thuyết $H_3$: Mô hình xung lượng tương đương với thời gian duy trì hiệu dụng $\tau_2 = 2I/P_{px}$ phản ánh hội tụ và bảo toàn năng lượng dao động của hệ).

Khung lý thuyết của công trình tích hợp giữa Lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin (FSDT), Lý thuyết vỏ thoải hai độ cong Donnell - Mushtari - Vlasov, Nguyên lý năng lượng biến phân Hamilton và Mô hình sóng xung kích nổ hóa học Sadovsky. Phạm vi nghiên cứu khảo sát các kết cấu panel trụ ($R_1 = R_x, R_2 = \infty$), panel cầu ($R_1 = R_2 = R$), với các tỷ số hình học $a/h \in [10, 100]$, $R/a \in [2, 10, \infty]$, tỷ số mô đun đàn hồi $E_c/E_t \in [0.1, 10]$, và hệ số độ cứng chốt $k_s/E_t \in [10^{-4}, 10^2]$ dưới các điều kiện biên ngàm bốn cạnh (C-C-C-C), tựa đơn bốn cạnh (S-S-S-S), và ngàm - tự do (C-F-C-F). Đóng góp mang tính đột phá của luận án nằm ở việc phát triển hoàn chỉnh thuật toán phần tử hữu hạn 8 nút đẳng tham số cho kết cấu vỏ composite sandwich có xét đến chuyển vị trượt trễ đàn hồi giữa các lớp, cung cấp công cụ tính toán chính xác có độ tin cậy được kiểm chứng với các kết quả quốc tế công bố trên Journal of Composite Structures và Computers & Structures.

Literature Review và Positioning

Cơ học kết cấu composite nhiều lớp đã phát triển qua nhiều thập kỷ với hai nhánh tiếp cận chính: Lý thuyết đơn lớp tương đương (Equivalent Single Layer - ESL) và Lý thuyết phân lớp (Layerwise Theory).

  • Nhánh ESL, khởi xướng từ lý thuyết tấm cổ điển (CLPT) theo giả thuyết Kirchhoff-Love, đã được mở rộng bởi Reissner (1945), Mindlin (1951) với lý thuyết biến dạng cắt bậc nhất (FSDT), và sau đó là lý thuyết biến dạng cắt bậc cao (HSDT) của Reddy (2004). Mặc dù ESL giảm thiểu đáng kể số lượng ẩn số chuyển vị, lý thuyết này bộc lộ hạn chế lớn khi không thể mô tả sự gián đoạn của đạo hàm chuyển vị (hiệu ứng zig-zag) và biến dạng cắt qua chiều dày tại mặt phân giới các lớp có cơ tính chênh lệch cao.
  • Nhánh Layerwise, được Pagano (1969, 1970) đề xuất qua nghiệm giải tích 3D đàn hồi và sau đó là Reddy (2004), Carrera (2002) với Khung công thức hợp nhất (Carrera Unified Formulation - CUF), cho phép mô tả chính xác trường ứng suất trượt ngang nhưng đòi hỏi chi phí tính toán cực lớn, đặc biệt khi chuyển sang các bài toán động lực học phi tuyến và vỏ cong hai chiều.
                      TIẾN TRÌNH LÝ THUYẾT CƠ HỌC COMPOSITE
                      
  Lý thuyết Cổ điển (CLPT)         Lý thuyết Đơn lớp tương đương (ESL/FSDT)
  (Kirchhoff-Love)                (Mindlin, Reissner, Reddy)
  Lý thuyết Phân lớp (Layerwise)    Lý thuyết Zig-Zag cải tiến
  (Pagano 1970, Reddy 2004)        (Carrera 2002, Cinefra 2014)
      KHUNG MÔ HÌNH HÓA ĐỘT PHÁ CỦA LUẬN ÁN (Trần Ngọc Cảnh, 2021)
      - Vỏ thoải composite sandwich 2 độ cong (SLCS)
      - Bổ sung bậc tự do trượt đàn hồi tại mặt tiếp xúc (k_s, k_sc)
      - Tích hợp đồng thời trường nhiệt dừng & Tải sóng nổ Sadovsky

Trong dòng nghiên cứu về dao động và đáp ứng động lực học của kết cấu vỏ, các công trình của Jansen (2005, 2007) sử dụng lý thuyết vỏ phi tuyến Donnell kết hợp phương pháp nhiễu loạn và xấp xỉ Galerkin để phân tích mất ổn định động của vỏ trụ tròn. Toorani & Lakis (2006) phát triển phần tử lai xét đến biến dạng trượt và quán tính quay. Về đáp ứng cơ - nhiệt, Biswal et al. (2016) khảo sát dao động vỏ thoải trong môi trường nhiệt ẩm sử dụng phần tử 8 nút; Chetan et al. (2017) phân tích ứng xử tĩnh của vỏ hai độ cong chịu tải nhiệt bằng FSDT và phần tử 9 nút MITC. Dưới tác dụng của tải trọng nổ, Emarti & Singha (2012) khảo sát panel trụ composite chịu tải nổ áp dụng tiêu chí phá hủy Tsai-Wu, trong khi Nguyễn Đình Đức và cộng sự (2017) giải tích dao động cưỡng bức vỏ composite có hệ số Poisson âm dưới tải sóng xung kích.

Tuy nhiên, tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai quan điểm mô hình hóa liên kết phân lớp:

  1. Trường phái bám dính lý tưởng (Perfect Bonding): Giả định vector chuyển vị liên tục tuyệt đối tại mặt phân giới ($u_t = u_c = u_b$ tại $z = z_i$), được áp dụng trong hầu hết các mô hình của Reddy (2004) và Biswal et al. (2016). Trường phái này cho rằng hiệu ứng trượt giữa các lớp là không đáng kể đối với kết cấu composite mỏng dán keo đồng nhất.
  2. Trường phái liên kết linh hoạt/trượt đàn hồi (Flexible/Shear Slip Bonding): Được bảo vệ bởi Lartham & Toledano (1988), Richard et al. (2010), Ranzi et al. (2004), lập luận rằng khi kết cấu ghép nối các lớp vật liệu có độ cứng chênh lệch cực lớn (như composite sandwich có lớp lõi mềm hoặc liên kết thép - bê tông sử dụng chốt neo), sự trượt đàn hồi tương đối $\Delta u = u_k - u_{k-1} \neq 0$ xuất hiện tất yếu và chi phối toàn bộ cơ chế hấp thụ năng lượng va đập. Bỏ qua độ mềm của chốt liên kết sẽ dẫn đến đánh giá sai lệch đáng kể (lên đến 25-40%) về độ võng động và ứng suất cắt cục bộ.

Luận án của tác giả Trần Ngọc Cảnh đã định vị chính xác vào điểm giao thoa chưa được khai phá này: Kế thừa tính ưu việt của lý thuyết FSDT trong việc giảm thiểu bậc tự do tính toán, nhưng bứt phá khỏi giả định bám dính tuyệt đối bằng cách thiết lập trường chuyển vị trượt độc lập cho từng phân lớp thông qua ma trận độ cứng liên kết chốt $k_s$.

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

  • So với mô hình của Emarti & Singha (2012) trên Composite Structures (chỉ khảo sát tải nổ trên panel trụ bám dính tuyệt đối không xét trường nhiệt), mô hình của luận án mở rộng cho vỏ hai độ cong bất kỳ chịu tải cơ - nhiệt đồng thời và có xét cơ cấu chốt trượt.
  • So với nghiên cứu của Chetan et al. (2017) (chỉ giải bài toán tĩnh cơ nhiệt trên vỏ liên tục), luận án đã giải quyết trọn vẹn bài toán động lực học quá độ (transient dynamic response) trong miền thời gian thực dưới xung áp suất nổ phi tuyến.

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

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

Công trình luận án mang lại những bước tiến lý thuyết quan trọng cho ngành cơ học vật rắn và tính toán kết cấu composite:

                  CƠ CHẾ BIẾN DẠNG TRƯỢT VỎ SANDWICH CÓ CHỐT

  1. Mở rộng lý thuyết biến dạng cắt bậc nhất FSDT cho hệ vỏ composite sandwich liên kết chốt trượt (SLCS): Luận án đã tổng quát hóa trường chuyển vị Reissner - Mindlin cho kết cấu 3 phân lớp (lớp mặt trên $t$, lớp lõi $c$, lớp mặt dưới $b$). Vector chuyển vị tại một điểm bất kỳ trong lớp thứ $k$ ($k=t, c, b$) được định nghĩa: $$u_k(x,y,z,t) = u_{0k}(x,y,0,t) + z \psi_{xk}(x,y,0,t)$$ $$v_k(x,y,z,t) = v_{0k}(x,y,0,t) + z \psi_{yk}(x,y,0,t)$$ $$w_k(x,y,z,t) = w_0(x,y,0,t)$$ Trong đó, chuyển vị pháp tuyến $w_0$ được giả thiết liên tục qua chiều dày, nhưng các chuyển vị tiếp tuyến $u_{0k}, v_{0k}$ và các góc xoay mặt cắt ngang $\psi_{xk}, \psi_{yk}$ của từng lớp là độc lập, cho phép mô phỏng chính xác hiện tượng trượt trượt tương đối tại mặt phân giới: $$\Delta u_{tc} = u_t(z=-h_t/2) - u_c(z=+h_c/2) \neq 0$$

  2. Thiết lập biểu thức giải tích tường minh cho thế năng biến dạng đàn hồi của hệ thống chốt liên kết: Thế năng biến dạng của chốt neo chống cắt trong một phần tử ($e_{sc}$) được xây dựng trực tiếp thông qua hệ số độ cứng cắt $k_s$ và $k_{sc}$: $$\Pi_{sc}^e = \frac{1}{2} \iint_{A^e} \left[ k_s^{(tc)} (\Delta u_{tc}^2 + \Delta v_{tc}^2) + k_s^{(cb)} (\Delta u_{cb}^2 + \Delta v_{cb}^2) \right] dx dy$$ Đóng góp này cho phép tích hợp trực tiếp cơ chế truyền lực cắt gián đoạn vào phiếm hàm Lagrange toàn cục của hệ thống kết cấu vỏ.

  3. Mô hình hóa tương tác liên kết cơ - nhiệt trong hệ phương trình vi phân dao động đa biến: Xây dựng hệ thức liên hệ biến dạng màng và mô men uốn tổng quát có xét đến hiệu ứng giãn nở nhiệt: $$\mathbf{N} = \mathbf{A}\boldsymbol{\varepsilon}^0_m + \mathbf{B}\boldsymbol{\kappa} - \mathbf{N}^T$$ $$\mathbf{M} = \mathbf{B}\boldsymbol{\varepsilon}^0_m + \mathbf{D}\boldsymbol{\kappa} - \mathbf{M}^T$$ $$\mathbf{Q} = \mathbf{H}\boldsymbol{\gamma}^s$$ Trong đó, vector lực màng nhiệt $\mathbf{N}^T$ và mô men nhiệt $\mathbf{M}^T$ được tích phân chính xác qua từng lớp theo hệ số giãn nở nhiệt định hướng $[\alpha_{x}, \alpha_{y}, \alpha_{xy}]^T_i$.

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

Khung phân tích của luận án kết hợp 3 trụ cột lý thuyết: Lý thuyết cơ học môi trường liên tục dị hướng, Lý thuyết động lực học kết cấu mỏng hai độ cong, và Phương pháp biến phân năng lượng.

Thành phần khung phân tích Cơ sở lý thuyết Bản chất đóng góp mới
Trường động học (Kinematics) FSDT + Độc lập chuyển vị tiếp tuyến từng phân lớp Mô tả chính xác chuyển vị trượt trôi tại mặt tiếp xúc mà không cần tăng bậc đa thức biến dạng qua chiều dày.
Mô hình tải nổ (Blast Dynamic) Quy luật biến thiên áp suất phản xạ $p_{px}(t) = P_{px}(1 - t/\tau_2)$ Cân bằng xung lượng thực tế từ công thức thực nghiệm Sadovsky, loại bỏ sai số bước nhảy vô hạn.
Trường nhiệt (Thermal Field) Nhiệt độ dừng đồng nhất $\Delta T = T - T_{ref}$ Xác định sự suy giảm độ cứng tổng thể và ứng suất dư màng trước khi chịu tải nổ động.

Biên điều kiện biên và giới hạn ứng dụng (Boundary conditions & Validity domain): Mô hình được xác lập chặt chẽ trong giới hạn biến dạng đàn hồi tuyến tính của vật liệu composite cốt sợi polymer (GFRP/CFRP), bỏ qua biến dạng co ngót thể tích theo phương pháp tuyến ($\varepsilon_z = 0$), và áp dụng cho vỏ thoải có tỷ số độ võng hình học thỏa mãn điều kiện vỏ nông ($w/a \ll 1, a/R \le 0.5$).

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ lập trường nhận thức luận thực chứng (Positivism) và cơ học tính toán tất định (Deterministic Computational Mechanics). Thiết kế nghiên cứu đa cấp độ (Multi-level computational design) được xây dựng từ quy mô lớp vật liệu đơn cơ sở (micro/macro-mechanics level) đến cấp độ phân tử phần tử hữu hạn và cấp độ toàn thể kết cấu kết hợp đa trường vật lý.

                   QUY TRÌNH PHÂN TÍCH PHẦN TỬ HỮU HẠN CƠ - NHIỆT

Quy trình nghiên cứu rigorous

Quy trình giải thuật phần tử hữu hạn được thiết lập thông qua các bước toán học nghiêm ngặt:

  1. Rời rạc hóa miền tính toán: Sử dụng phần tử vỏ thoải 8 nút đẳng tham số dạng tấm/vỏ cong bậc hai (Serendipity). Mỗi nút phần tử mang các bậc tự do chuyển vị dài và góc xoay tương ứng với cấu trúc phân lớp. Đối với vỏ sandwich 3 lớp liên kết chốt, vector chuyển vị nút phần tử bao gồm 13 bậc tự do: $$\mathbf{q}j = \left[ u{0t}, v_{0t}, \psi_{xt}, \psi_{yt}, u_{0c}, v_{0c}, \psi_{xc}, \psi_{yc}, u_{0b}, v_{0b}, \psi_{xb}, \psi_{yb}, w_0 \right]^T_j$$
  2. Nội suy hàm dạng (Shape Functions): Chuyển vị trong phần tử được xấp xỉ qua các hàm dạng bậc hai $N_j(\xi, \eta)$: $$\mathbf{u}^e = \sum_{j=1}^8 N_j(\xi, \eta) \mathbf{q}_j^e$$
  3. Thành lập phương trình Lagrange loại 2: $$\frac{d}{dt}\left( \frac{\partial L_e}{\partial \mathbf{\dot{q}}e} \right) - \frac{\partial L_e}{\partial \mathbf{q}e} = 0$$ Với hàm Lagrange $L_e = T_e - \Pi_e$, trong đó $T_e$ là động năng phần tử có xét quán tính quay, $\Pi_e = \Pi{esh} + \Pi{esc} - A_e - A_T$ là tổng thế năng biến dạng của vỏ, thế năng chốt, công ngoại lực sóng nổ và công tải nhiệt.
  4. Hệ phương trình vi phân dao động động lực học toàn cục: $$\mathbf{M} \mathbf{\ddot{q}}(t) + \mathbf{C} \mathbf{\dot{q}}(t) + \mathbf{K} \mathbf{q}(t) = \mathbf{P}(t) + \mathbf{F}^T$$ Trong đó $\mathbf{M}$ là ma trận khối lượng nhất quán, $\mathbf{C} = \alpha_R \mathbf{M} + \beta_R \mathbf{K}$ là ma trận cản Rayleigh, $\mathbf{K}$ là ma trận độ cứng tổng thể (bao gồm độ cứng màng, uốn, cắt của vỏ và độ cứng cắt của chốt), $\mathbf{P}(t)$ là vector tải trọng sóng xung kích nút, và $\mathbf{F}^T$ là vector lực nút tương đương do biến dạng nhiệt.

Data và phân tích

Để thực thi tính toán số, tác giả đã lập trình toàn bộ thuật toán trên nền tảng tính toán khoa học chuyên dụng (MATLAB/Fortran code tự phát triển), xây dựng module giải phương trình vi phân dao động theo phương pháp tích phân từng bước Newmark-$\beta$ ($\beta = 1/4, \gamma = 1/2$ đảm bảo tính ổn định vô điều kiện) kết hợp kỹ thuật khử Gauss giải hệ đại số tuyến tính quy mô lớn.

Dữ liệu đặc trưng vật liệu composite kiểm chuẩn:

  • Lớp composite sợi đồng phương (Graphite/Epoxy hoặc Glass/Epoxy): $E_1 = 132.38 \text{ GPa}$, $E_2 = 10.756 \text{ GPa}$, $G_{12} = G_{13} = 5.6537 \text{ GPa}$, $G_{23} = 3.378 \text{ GPa}$, $\nu_{12} = 0.24$, khối lượng riêng $\rho = 1578 \text{ kg/m}^3$, hệ số giãn nở nhiệt $\alpha_1 = 1.2 \times 10^{-6} \text{ /K}$, $\alpha_2 = 24.0 \times 10^{-6} \text{ /K}$.
  • Lớp lõi sandwich đàn nhớt/xốp: $E_c = 0.28 \text{ GPa}$, $G_c = 0.11 \text{ GPa}$, $\rho_c = 100 \text{ kg/m}^3$.

Kiểm tra độ tin cậy và hội tụ lưới phần tử: Lưới phần tử hữu hạn được khảo sát từ $4 \times 4$, $8 \times 8$, $10 \times 10$ đến $16 \times 16$ phần tử trên toàn miền vỏ. Kết quả tần số dao động riêng không thứ nguyên $\bar{\omega} = \omega (a^2/h)\sqrt{\rho/E_2}$ của panel cầu composite 4 lớp $[0^\circ/90^\circ/90^\circ/0^\circ]$ đạt độ hội tụ hoàn hảo tại lưới $8 \times 8$ (sai số $< 0.05%$). So sánh trực tiếp với kết quả của Reddy (2004) và phần mềm ANSYS cho thấy độ lệch tuyệt đối không vượt quá $0.82%$, chứng minh độ chính xác tuyệt đối của chương trình tính tự lập trình.

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

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

Kết quả khảo sát tham số số trị quy mô lớn trong luận án đã phát hiện những quy luật cơ học đột phá:

             ĐÁP ỨNG ĐỘNG THEO ĐỘ CỨNG CHỐT k_s VÀ NHIỆT ĐỘ
             
  Độ võng w_max (mm)                       Trượt tương đối V_tc (mm)
  1. Hiệu ứng phi tuyến của độ cứng chốt liên kết ($k_s$) đến chuyển vị trượt và độ võng động: Khi độ cứng chống cắt của chốt tăng từ $k_s/E_t = 10^{-4}$ (liên kết trượt tự do) lên $k_s/E_t = 10^2$ (liên kết bám dính cứng hoàn toàn), độ võng cực đại chính giữa vỏ $w_{\max}$ giảm mạnh tới $38.4%$, trong khi độ trượt tương đối giữa các lớp $V_{tc\max}$ giảm tiệm cận về $0$. Kết quả này khẳng định chốt liên kết đóng vai trò truyền mô men và lực cắt chuyển tiếp, kiểm soát trực tiếp trạng thái biến dạng tổng thể của hệ.
  2. Quy luật tương tác bất lợi giữa nhiệt độ môi trường và tải nổ: Khi nhiệt độ tăng từ $\Delta T = 0\text{ K}$ lên $\Delta T = 100\text{ K}$, tần số dao động riêng bậc 1 của vỏ thoải giảm từ $12.6%$ đến $18.3%$ do hiệu ứng giãn nở nhiệt gây ra ứng suất nén màng cục bộ làm giảm độ cứng chống uốn tổng thể. Dưới tác dụng đồng thời của sóng nổ, ứng suất kéo - nén lớn nhất $\sigma_{\max}$ tăng vọt $24.7%$ so với trường hợp chỉ chịu tải nổ ở nhiệt độ phòng.
  3. Ảnh hưởng của tỷ số độ cong $R/a$ và hình học vỏ: Vỏ có độ cong Gauss dương (panel cầu $R_x = R_y$) có khả năng chịu tải nổ vượt trội hơn panel trụ ($R_y = \infty$) và tấm phẳng ($R \to \infty$). Khi tỷ số độ cong $R/a$ giảm từ $10$ xuống $2$ (vỏ vồng cao hơn), độ võng cực đại $w_{\max}$ giảm đến $52.1%$ nhờ hiệu ứng truyền lực dọc trục màng (membrane action) thay thế cho cơ chế chịu uốn thuần túy.
  4. Tác động của cấu trúc chiều dày lõi ($h_c/h_t$): Tăng tỷ số chiều dày lõi trên chiều dày lớp mặt $h_c/h_t$ từ $2$ lên $8$ làm tăng mô men quán tính uốn của vỏ sandwich, giúp giảm biên độ rung động $w_{\max}$ tới $31.2%$ mà khối lượng kết cấu chỉ tăng không đáng kể.
Trường hợp khảo sát $w_{\max}$ (mm) $\sigma_{\max}$ (MPa) $V_{tc\max}$ (mm) $\bar{\omega}_1$ (không thứ nguyên)
Panel cầu ($k_s/E_t = 10^{-4}, \Delta T = 0\text{K}$) $4.82$ $185.4$ $0.142$ $14.22$
Panel cầu ($k_s/E_t = 1.0, \Delta T = 0\text{K}$) $3.25$ $210.8$ $0.038$ $18.65$
Panel cầu ($k_s/E_t = 10^2, \Delta T = 0\text{K}$) $2.97$ $228.1$ $0.002$ $19.84$
Panel cầu ($k_s/E_t = 1.0, \Delta T = 100\text{K}$) $4.05$ $262.9$ $0.045$ $15.82$
Panel trụ ($k_s/E_t = 1.0, \Delta T = 0\text{K}$) $4.68$ $245.2$ $0.052$ $15.11$

Implications đa chiều

  • Tiến bộ lý thuyết: Đặt nền móng giải tích - số hoàn chỉnh cho phân hệ kết cấu nhiều lớp có liên kết đàn hồi gián đoạn, mở rộng biên giới của lý thuyết FSDT sang miền bài toán động lực học phi bám dính.
  • Đổi mới phương pháp luận: Cung cấp thuật toán phần tử hữu hạn 8 nút hiệu quả cao, tiết kiệm hơn $60%$ chi phí thời gian tính toán so với mô hình 3D Solid Layerwise trong ANSYS/Abaqus mà vẫn phản ánh trung thực trường trượt tiếp tuyến.
  • Ứng dụng thực tiễn công kỹ thuật: Cung cấp biểu đồ và công thức tra cứu thực dụng cho kỹ sư kết cấu trong việc lựa chọn mật độ chốt neo, độ cứng chốt $k_s$, cấu hình xếp lớp $[0^\circ/90^\circ]_\text{sym}$ và tỷ lệ chiều dày lõi tối ưu khi thiết kế vòm bảo vệ công sự, vỏ xe bọc thép, kết cấu tàu chiến và panel bảo vệ công trình năng lượng hạt nhân.

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 chỉ ra các giới hạn nghiên cứu nội tại:

  1. Giả thiết cơ tính vật liệu không phụ thuộc vào nhiệt độ ($E(T) = \text{const}$) và trường nhiệt độ là trường dừng đồng nhất. Trong thực tế các vụ nổ nhiệt độ cao tức thời, trường nhiệt biến thiên phi tuyến theo chiều dày và thời gian ($T(z, t)$).
  2. Chưa xét đến hiệu ứng phi tuyến hình học Von Kármán (biến dạng lớn uốn lớn) và sự suy giảm độ cứng dẻo/tổn thương tích lũy của chốt liên kết theo thời gian chịu tải động.
  3. Nghiên cứu tập trung vào phương pháp tính toán lý thuyết và mô phỏng số, chưa có điều kiện tiến hành thử nghiệm nổ thực địa trong hầm nổ quân sự do chi phí và điều kiện an toàn đặc thù.

Chương trình nghiên cứu tiếp nối khuyến nghị 4 hướng mở rộng cụ thể:

  • Phát triển mô hình bài toán truyền nhiệt không dừng phi tuyến $T(x,y,z,t)$ kết hợp cơ tính biến thiên theo nhiệt độ $E(T), \alpha(T)$.
  • Tích hợp mô hình phá hủy động học liên tục (Continuum Damage Mechanics) và tiêu chí tách lớp phá hủy chốt neo dưới tải lặp xung kích.
  • Mở rộng thuật toán cho kết cấu vỏ composite vật liệu có cơ tính biến thiên (FGM) gia cường hạt nano hoặc cấu trúc mạng rỗng (auxetic metamaterials).
  • Tiến hành thực nghiệm đo biến dạng động bằng cảm biến sợi quang FBG và camera quang học tốc độ cao (DIC) trên mẫu vỏ sandwich chịu nổ khí nén (shock tube testing).

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

Công trình luận án mang lại giá trị tác động sâu rộng trên nhiều phương diện:

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

  • Nghiên cứu sinh & Học viên cao học chuyên ngành Cơ kỹ thuật / Cơ học vật rắn: Thừa hưởng hệ thống công thức giải tích chuẩn xác, tường minh ma trận phần tử hữu hạn và mã nguồn giải thuật động lực học đa trường vật lý.
  • Các nhà khoa học & Giảng viên đại học: Sử dụng làm tài liệu tham khảo chuyên sâu chuẩn mực cho các học phần Tính toán kết cấu composite tiên tiến, Động lực học công trình và Phương pháp phần tử hữu hạn nâng cao.
  • Kỹ sư R&D tại các Viện nghiên cứu Quân sự & Công nghiệp Quốc phòng: Nắm vững công cụ định lượng để tính toán độ bền kết cấu kháng chấn, kháng nổ chịu nhiệt mà không phụ thuộc hoàn toàn vào các phần mềm thương mại đắt tiền từ nước ngoài.
  • Các nhà hoạch định tiêu chuẩn thiết kế xây dựng công sự: Sở hữu cơ sở khoa học tin cậy để biên soạn các quy chuẩn kỹ thuật quốc gia về thiết kế kết cấu nhiều lớp chịu tác động nổ hóa học.

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 lý thuyết độc đáo nhất là việc thiết lập trường chuyển vị và hàm thế năng biến dạng giải tích cho liên kết chốt trượt đàn hồi ($k_s, k_{sc}$) trên kết cấu vỏ thoải hai độ cong. Luận án đã mở rộng Lý thuyết biến dạng cắt bậc nhất Reissner - Mindlin (FSDT) từ trạng thái bám dính đơn lớp tương đương (ESL) truyền thống sang mô hình 3 phân lớp có bậc tự do chuyển vị tiếp tuyến độc lập, giải quyết triệt để bài toán trượt gián đoạn giữa các lớp mà vẫn giữ nguyên tính kinh tế trong tài nguyên tính toán.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 02 nghiên cứu quốc tế đã công bố?

So với nghiên cứu của Biswal et al. (2016) (chỉ giải bài toán dao động riêng nhiệt ẩm trên vỏ bám dính hoàn toàn) và Emarti & Singha (2012) (phân tích nổ trên panel trụ thuần cơ học), luận án đã:

  • Tích hợp đồng thời bài toán ghép đa trường: Trường nhiệt độ dừng $\Delta T$ gây ứng suất trước + Tải trọng nổ phản xạ thời gian thực $p_{px}(t)$.
  • Xây dựng phần tử vỏ 8 nút 13 bậc tự do/nút tích hợp độ cứng chống cắt của chốt neo liên kết, mô phỏng trực tiếp chuyển vị trượt trôi $V_{tc}(t)$ trong miền thời gian thực thông qua thuật toán Newmark-$\beta$.

3. Phát hiện bất ngờ nhất (surprising finding) rút ra từ dữ liệu tính toán là gì?

Phát hiện bất ngờ nhất là tồn tại một ngưỡng bão hòa của độ cứng chốt liên kết ($k_s/E_t \approx 10^1 \div 10^2$). Khi $k_s/E_t < 1.0$, độ võng và ứng suất trượt nhạy cảm cực cao với sự biến thiên độ cứng chốt (giảm $38%$ độ võng khi tăng $k_s$). Tuy nhiên, khi vượt qua ngưỡng $10^2$, việc tăng thêm độ cứng chốt hầu như không làm giảm thêm độ võng $w_{\max}$ ($< 1.2%$) nhưng lại làm tăng vọt tập trung ứng suất cắt cục bộ quanh vị trí liên kết, dễ gây phá hủy giòn cục bộ của lớp mặt composite.

                  NGƯỠNG BÃO HÒA ĐỘ CỨNG CHỐT LIÊN KẾT (k_s)
                  
  Độ võng w_max (mm)
                           10^2

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

Có. Luận án cung cấp đầy đủ:

  • Hệ thống công thức giải tích tường minh các ma trận $[A], [B], [D], [H]$ và vector tải trọng nhiệt $\mathbf{N}^T, \mathbf{M}^T$.
  • Biểu thức giải tích của các ma trận độ cứng phần tử vỏ $[\mathbf{K}_e^{sh}]$ và chốt $[\mathbf{K}_e^{sc}]$.
  • Bảng thông số vật liệu chuẩn xác, điều kiện biên toán học, quy luật tải trọng sóng xung kích Sadovsky và sơ đồ khối thuật toán chương trình tính chi tiết.

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

  • Giai đoạn 1 (1-3 năm): Mở rộng giải thuật cho bài toán phi tuyến hình học Von Kármán và phi tuyến vật liệu của chốt trượt chịu tải nổ lặp chu kỳ.
  • Giai đoạn 2 (4-6 năm): Tích hợp trường nhiệt độ không dừng $T(z,t)$ và tối ưu hóa đa mục tiêu cấu trúc composite sandwich FGM auxetic.
  • Giai đoạn 3 (7-10 năm): Chế tạo thực nghiệm, phát triển công nghệ cảm biến thông minh tự chẩn đoán tách lớp và thương mại hóa phần mềm phần tử hữu hạn chuyên dụng cho công trình quốc phòng tại Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Trần Ngọc Cảnh là một công trình khoa học mẫu mực, hoàn chỉnh và có tính đột phá cao trong lĩnh vực cơ học vật rắn và kỹ thuật kết cấu composite. Sáu đóng góp cụ thể cốt lõi bao gồm:

  1. Xây dựng thành công mô hình toán học giải tích cho vỏ thoải composite lớp và vỏ thoải composite sandwich 3 phân lớp có liên kết chốt trượt hai độ cong chịu tác dụng đồng thời của tải trọng sóng xung kích và nhiệt độ.
  2. Thiết lập hệ phương trình vi phân dao động cưỡng bức đa trường vật lý dựa trên nguyên lý biến phân Hamilton và lý thuyết biến dạng cắt bậc nhất FSDT.
  3. Phát triển giải thuật và chương trình tính phần tử hữu hạn 8 nút đẳng tham số với độ chính xác cao và tốc độ hội tụ vượt trội.
  4. Phát hiện quy luật tác động đa thông số: Tỷ số độ cong $R/a$, độ cứng chốt $k_s/E_t$, tỷ số chiều dày lõi $h_c/h_t$ và độ chênh nhiệt độ $\Delta T$ đến đáp ứng động lực học của hệ kết cấu.
  5. Xác lập ngưỡng bão hòa độ cứng chốt neo, cung cấp cơ sở khoa học cho thiết kế tối ưu hóa liên kết trong công trình phòng thủ và hàng không.
  6. Mở ra 3 nhánh nghiên cứu chuyên sâu mới về cơ học kết cấu composite phi bám dính, vật liệu kháng nổ thông minh và động lực học công trình đa trường vật lý tại Việt Nam.