Tổng quan luận án

Luận án tiến sĩ kỹ thuật "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 độ" do nghiên cứu sinh Trần Ngọc Cảnh thực hiện 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 Bộ môn Cơ học vật rắn, Khoa Cơ khí, Học viện Kỹ thuật Quân sự phối hợp cùng Trường Sĩ quan Kỹ thuật Quân sự. Công trình thuộc chuyên ngành Cơ kỹ thuật / Cơ học vật rắn.

Vật liệu composite phân lớp và kết cấu composite sandwich liên kết chốt (shear connector) được ứng dụng rộng rãi trong các công trình kỹ thuật dân dụng và quân sự (vỏ máy bay, tàu vũ trụ, xuồng cao tốc, vỏ bảo vệ chống nổ, mái che công trình kiến trúc). Trong điều kiện làm việc thực tế, các kết cấu này thường xuyên chịu tác động đồng thời của tải trọng động lực học cường độ cao từ sóng xung kích (do các vụ nổ gây ra) và trường nhiệt độ môi trường. Khi các lớp vật liệu có cơ tính khác biệt lớn ghép nối với nhau, hiện tượng trượt và bong tách giữa các lớp dễ dàng xuất hiện. Giải pháp sử dụng chốt liên kết cho phép các lớp trượt đàn hồi cục bộ nhằm giải tỏa ứng suất và ngăn ngừa phá hủy liên kết. Tuy nhiên, phần lớn các công trình đã công bố chỉ dừng lại ở bài toán dầm hoặc tấm phẳng chịu tải đơn lẻ; các nghiên cứu về mô hình toán học và phương pháp số tính toán cho vỏ thoải composite lớp bám dính tuyệt đối cũng như vỏ sandwich có liên kết chốt chịu đồng thời tải trọng sóng xung kích và nhiệt độ còn rất hạn chế.

Mục tiêu và các nội dung nghiên cứu của luận án được xác định cụ thể như sau:

  • Mục tiêu nghiên cứu: Góp phần hoàn thiện phương pháp tính toán lý thuyết các kết cấu composite khác nhau dưới tác dụng của các dạng tải trọng khác nhau, đáp ứng nhu cầu ứng dụng kết cấu trong thực tiễn.
  • Các nội dung nghiên cứu chính (được đánh số theo văn bản):
    1. Xây dựng mô hình bài toán, thiết lập các quan hệ ứng xử cơ học và thành lập phương trình dao động cho vỏ thoải composite lớp và vỏ thoải composite sandwich có liên kết chốt chịu tác dụng đồng thời tải trọng sóng xung kích và nhiệt độ.
    2. Xây dựng thuật toán phần tử hữu hạn (PTHH) và lập chương trình tính, để phân tích động lực học cho kết cấu vỏ thoải composite lớp và vỏ thoải composite sandwich.
    3. Khảo sát ảnh hưởng của một số yếu tố nhiệt độ, hình học, vật liệu,... đến phản ứng động của kết cấu chịu đồng thời tải trọng sóng xung kích và nhiệt độ. Trên cơ sở đó đề xuất các phương án hợp lý trong thiết kế, chế tạo, khai thác, sửa chữa đối với kết cấu loại này trong thực tiễn và trong kỹ thuật.

Đối tượng nghiên cứu của luận án là kết cấu vỏ thoải composite lớp (LCS) và vỏ thoải composite sandwich (SLCS) có liên kết chốt, ở hai dạng hình học cơ bản là panel vỏ trụ thoải và panel vỏ thoải hai độ cong (panel vỏ cầu). Phạm vi nghiên cứu giới hạn ở việc xác định đáp ứng động lực học (chuyển vị, góc xoay, trượt tương đối giữa các lớp và ứng suất) của kết cấu vỏ thoải làm việc trong miền đàn hồi tuyến tính, chịu tải trọng sóng xung kích do nổ hóa học trong không khí và trường nhiệt độ dừng phân bố đều.

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

Luận án đã tổng quan hệ thống các hướng nghiên cứu trong nước và quốc tế liên quan đến cơ học kết cấu tấm, vỏ composite:

  1. Hướng nghiên cứu tấm và vỏ composite phân lớp bám dính tuyệt đối trên thế giới:

    • Phương pháp giải tích và nửa giải tích cho vỏ mỏng và dày: Wu và Vison (1971) xét dao động tự do tấm composite có biến dạng trượt ngang; Calcote, Librescu phân tích tĩnh và động tấm vỏ dị hướng; Franklin khảo sát nghiệm màng vỏ composite tròn xoay; Rizzo và Vicano thiết lập bài toán 3 chiều cho ống dị hướng; Pagano nghiên cứu bài toán biến dạng trượt ngang cho dầm, tấm vuông và ống composite; Whitney áp dụng lý thuyết vỏ Ambartsumyan cho vỏ trụ đa lớp; Reuter dùng phương trình Donnell phân tích vỏ trụ có cốt lệch trục; Bert, Martin, Kunukkasseril nghiên cứu vỏ mỏng dị hướng, vỏ nón và vỏ trụ tròn xoay; Hildebrand, Ambartsumyan, Kalnins, Almed phát triển lý thuyết vỏ trực hướng có xét biến dạng trượt ngang. Vinson, Mehta, Renton, Summer, Walts, Zukas nghiên cứu các hàm xấp xỉ phân bố ứng suất trượt qua chiều dày.
    • Phương pháp số và phần tử hữu hạn: Toorani và Lakis dùng phần tử lai phân tích dao động phi tuyến vỏ trụ; Jansen công bố chuỗi công trình sử dụng lý thuyết vỏ phi tuyến Donnell kết hợp phương pháp phân tích đơn giản, mở rộng và khai triển nhiễu loạn để khảo sát dao động phi tuyến, mất ổn định động và hiện tượng flutter; Rahman và cộng sự kết hợp PTHH với phương pháp nhiễu loạn; Kubenko, Kovalchuk, Mallon, Bespalova, Urusova nghiên cứu lý thuyết và thực nghiệm ổn định động vỏ composite trực hướng và 3 lớp cốt sợi thủy tinh; Shariyat dùng công thức tọa độ tổng thể chồng chất kép cho vỏ composite đàn nhớt; Lugovoi, Dogan, Vaicaitis, Mohammadi, Sedaghati nghiên cứu dao động phi tuyến vỏ sandwich 3 lớp lõi mềm và đàn nhớt; Dong, Wang, Kozlov nghiên cứu lan truyền sóng và giảm chấn trong vỏ gắn lớp áp điện dưới tác dụng nhiệt. Cinefra và Valvano dùng phần tử 9 nút MITC chống khóa cắt; Francesco ứng dụng công thức hợp nhất Carrera (CUF) và phương pháp cầu phương vi phân (GDQ); Sushree và cs phân tích tĩnh và dao động vỏ thoải composite kết hợp so sánh phần mềm Ansys; Sandipan và cs phát triển phần tử 8 nút 7 bậc tự do cho vỏ hai độ cong; Wei Li và Haijun Shen dùng lý thuyết Layerwise cho vỏ composite áp điện; Emarti và Singha khảo sát panel trụ composite chịu tải nổ parabol theo tiêu chí phá hủy Tsai-Wu.
    • Bài toán chịu tác động nhiệt độ và nhiệt ẩm: Pipes và cs dùng biến đổi Laplace tính ứng suất nhiệt; Parhi, Bhattacharya, Sinha nghiên cứu tách lớp trong tấm composite carbon/epoxy; Krishna, Patel, Ganapathi dùng mô hình bậc cao; Garber, Mcdonaugh, Kliger, Daugherty, Vinson khảo sát composite pyrolic/graphite; Biswal và cs nghiên cứu thực nghiệm và mô phỏng vỏ composite hai độ cong Glass/Epoxy trong môi trường nhiệt ẩm; Rahmani và cs khảo sát vỏ trụ sandwich lõi mềm; Li và Qiao nghiên cứu vỏ trụ mỏng trên nền đàn hồi Pasternak; Chetan và cs nghiên cứu vỏ composite 2 độ cong chịu tải cơ nhiệt bằng phần tử 9 nút.

  2. Hướng nghiên cứu kết cấu tấm và vỏ composite tại Việt Nam:

    • Đào Huy Bích, Vũ Khắc Bảy ứng dụng thuần nhất hóa trong ổn định phi tuyến vỏ trụ; Phạm Thị Toan, Đào Huy Bích dùng lý thuyết mô đun hiệu quả tính vỏ lượn sóng và tấm vỏ có cấu trúc tựa tuần hoàn; Hoàng Xuân Lượng, Phạm Tiến Đạt, Nguyễn Minh Tuấn nghiên cứu bài toán tĩnh và động cho tấm vỏ mỏng và ống dày bằng FSDT và PTHH; Trần Ích Thịnh và Ngô Như Khoa nghiên cứu thực nghiệm môi trường hóa học, bài toán dao động tấm composite chịu nhiệt ẩm dùng mô hình biến dạng bậc 3; Nguyễn Hoa Thịnh khảo sát vỏ trụ trong môi trường không đồng nhất chịu tải xung; Tạ Thị Hiền nghiên cứu dao động tấm vỏ tương tác chất lỏng bằng phương pháp độ cứng động; Nguyễn Văn Hưng và cs phân tích dao động vỏ trụ thoải composite chịu sóng xung kích trong nước; Trần Minh Tú và cs tính vỏ 2 độ cong chịu uốn trong môi trường nhiệt độ bằng PTHH; Nguyễn Mạnh Cường và cs phân tích vỏ trụ và vỏ nón đối xứng trên nền Winkler bằng phần tử liên tục có xét gia tốc Coriolis; Nguyễn Đình Đức và cs giải tích dao động vỏ composite hệ số Poisson âm chịu tải nổ bằng phương pháp Galerkin và Runge-Kutta bậc 4; Trịnh Anh Tuấn và cs nghiên cứu dao động vỏ 2 độ cong có gân gia cường.

  3. Hướng nghiên cứu kết cấu nhiều lớp liên kết chốt hoặc lớp đàn hồi:

    • Về dầm nhiều lớp liên kết chốt: Đã có một số công trình công bố áp dụng phương pháp giải tích và phương pháp số (tài liệu [26], [32]-[35], [38], [41], [61], [63], [65], [70], [87], [90], [106]-[110], [120], [135]).
    • Về tấm và vỏ: Lartham và Toledano nghiên cứu uốn tấm hai lớp có trượt tuyến tính; Peng và cs đề xuất mô hình 3D kết hợp 1D cho tấm có lớp xen giữa; Richard và cs thực nghiệm kết cấu sandwich composite thép - bê tông - thép liên kết chốt chữ J chịu tải nổ.

Khoảng trống nghiên cứu được luận án lựa chọn giải quyết là: Thiết lập mô hình toán học và xây dựng thuật toán phần tử hữu hạn hoàn chỉnh để tính toán đáp ứng động lực học của vỏ thoải composite lớp bám dính tuyệt đối và vỏ thoải composite sandwich có liên kết chốt (kể đến chuyển vị trượt đàn hồi tương đối giữa các lớp tại chốt) dưới tác dụng đồng thời của tải trọng sóng xung kích trong không khí và trường nhiệt độ dừng.

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

Luận án xây dựng mô hình dựa trên các cơ sở lý thuyết và phương pháp sau:

  1. Mô hình hình học và trường chuyển vị:

    • Sử dụng lý thuyết vỏ thoải hai độ cong với kích thước mặt bằng $a \times b$, chiều dày $h$, bán kính cong chính $R_x, R_y$ và bán kính xoắn $R_{xy}$. Các trường hợp riêng gồm: panel trụ ($R_1 \neq 0, R_2 = \infty$), panel cầu ($R_1 = R_2$), vỏ có độ cong Gauss dương ($R_1 R_2 > 0$) và vỏ yên ngựa có độ cong Gauss âm ($R_1 R_2 < 0$).
    • Áp dụng lý thuyết biến dạng cắt bậc nhất (FSDT) theo giả thiết Reissner - Mindlin ($\sigma_z = 0$), trường chuyển vị tại điểm bất kỳ có tọa độ $(x, y, z)$ ở thời điểm $t$ có dạng: $$u(x, y, z, t) = u_0(x, y, 0, t) + z \phi_x(x, y, 0, t)$$ $$v(x, y, z, t) = v_0(x, y, 0, t) + z \phi_y(x, y, 0, t)$$ $$w(x, y, z, t) = w_0(x, y, 0, t)$$ trong đó $u_0, v_0, w_0$ là chuyển vị tại mặt trung bình; $\phi_x, \phi_y$ là góc xoay của tiết diện ngang quanh các trục $y$ và $x$.

  2. Quan hệ ứng suất - biến dạng và nội lực:

    • Vật liệu mỗi lớp là composite cốt sợi đồng phương, tuân theo định luật Hooke cho môi trường trực hướng phẳng với ma trận độ cứng màng $[Q_m]$ và ma trận độ cứng cắt $[Q_s]$.
    • Biến dạng nhiệt của lớp thứ $i$ được xác định qua các hệ số giãn nở nhiệt chính $\alpha_1, \alpha_2$, góc đặt sợi $\theta$ và độ biến thiên nhiệt độ $\Delta T = T - T_{ref}$.
    • Các thành phần nội lực phần tử gồm véc tơ lực màng $N$, mô men uốn - xoắn $M$ và lực cắt $Q$ được liên hệ với biến dạng thông qua các ma trận độ cứng $A$ (màng), $B$ (tương tác màng - uốn), $D$ (uốn - xoắn), $H$ (cắt) cùng các véc tơ lực nhiệt tương đương $N^T, M^T$.

  3. Mô hình tải trọng sóng xung kích (SXK) do nổ hóa học:

    • Áp suất đỉnh $\Delta P$ của sóng tới trong không khí xác định theo công thức thực nghiệm của Sa-đốp-sky theo khoảng cách $R$ và đương lượng nổ $C$: $$\Delta P = 0{,}84 \frac{\sqrt[3]{C}}{R} + 2{,}7 \frac{\sqrt[3]{C^2}}{R^2} + 7{,}0 \frac{C}{R^3} \quad (\text{kG/cm}^2)$$
    • Áp suất phản xạ cực đại $P_{px}$ tác dụng lên bề mặt kết cấu khi góc tới $\alpha < 40^\circ$: $$P_{px} = 2 \Delta P + \frac{6 \Delta P^2}{\Delta P + 7 P_0} \cos^2 \alpha$$
    • Luận án sử dụng mô hình tải trọng sóng xung kích tuyến tính hóa dạng tam giác cân bằng xung lượng với thời gian tác dụng hiệu dụng $\tau_2 = \frac{2 I^+}{P_{px}}$ để tính toán động lực học.

  4. Mô hình phần tử hữu hạn và liên kết chốt:

    • Vỏ composite lớp thông thường (LCS): Xây dựng phần tử vỏ thoải nhiều lớp dựa trên nguyên lý Hamilton và phương trình Lagrange.
    • Vỏ composite sandwich có chốt liên kết (SLCS): Kết cấu gồm 3 lớp phân tách (lớp trên $t$, lớp giữa $c$, lớp dưới $b$). Mô hình phần tử 8 nút 3 lớp được thiết lập, tại mỗi nút có các bậc tự do chuyển vị riêng cho từng lớp. Sự trượt đàn hồi tại mặt tiếp xúc giữa các lớp được kiểm soát bởi độ cứng chống cắt của chốt $k_{sc}$ ($k_s$), đóng góp vào thế năng biến dạng của chốt $U_e^{sc}$.
    • Phương trình vi phân dao động tổng thể có dạng: $$[M] {\ddot{q}} + [C] {\dot{q}} + [K] {q} = {P_{SXK}(t)} + {F^T}$$
Thông số / Thành phần Ký hiệu Ý nghĩa cơ học
Kích thước hình học vỏ $a, b, h$ Chiều rộng, chiều dài và tổng chiều dày của vỏ
Bán kính độ cong $R_x, R_y, R_{xy}$ Bán kính cong chính theo trục $x, y$ và bán kính xoắn
Chiều dày từng lớp sandwich $h_t, h_c, h_b$ Chiều dày lớp trên ($t$), lớp giữa ($c$) và lớp dưới ($b$)
Mô đun đàn hồi từng lớp $E_t, E_c, E_b$ Mô đun đàn hồi kéo/nén của các lớp thành phần
Độ cứng chống cắt của chốt $k_s$ ($k_{sc}$) Khả năng chống trượt đàn hồi giữa các lớp liên kết
Tải trọng sóng xung kích $p_{px}(t), P_{px}, \tau_2$ Áp suất tức thời, áp suất đỉnh phản xạ và thời gian hiệu dụng
Tải trọng nhiệt độ $\Delta T, {N^T}, {M^T}$ Biến thiên nhiệt độ, véc tơ lực màng nhiệt và mô men nhiệt

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 lý thuyết về vật liệu composite phân lớp, các mô hình tính toán cơ học vĩ mô và vi mô, lý thuyết đơn lớp tương đương (ESL), lý thuyết lớp (Layerwise) và lý thuyết dích dắc (zig-zag). Tác giả phân tích các phương pháp tính toán vỏ mỏng và vỏ dày composite, điều kiện Codaxi-Gauss đối với vỏ có độ vồng lớn và tính ưu việt của lý thuyết vỏ thoải hai độ cong. Chương này tổng kết cơ chế lan truyền sóng xung kích nổ hóa học trong không khí, công thức thực nghiệm Sa-đốp-sky và quy luật biến thiên áp suất sóng phản xạ. Đồng thời, tác giả điểm lại chi tiết các công trình nghiên cứu tấm, vỏ composite và kết cấu liên kết chốt trong và ngoài nước, từ đó xác định rõ 3 nội dung trọng tâm cần giải quyết của luận án.

Chương 2: 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 ĐỘ

Chương 2 thiết lập mô hình toán học cho vỏ thoải composite lớp bám dính tuyệt đối, các lớp bố trí đối xứng qua mặt trung bình, vật liệu cốt sợi đồng phương. Trên cơ sở lý thuyết biến dạng cắt bậc nhất Mindlin và nguyên lý năng lượng, tác giả xây dựng công thức tính thế năng biến dạng màng, uốn, cắt, thế năng do nhiệt độ, động năng và công ngoại lực của sóng xung kích. Thuật toán phần tử hữu hạn được xây dựng để rời rạc hóa kết cấu thành hệ phương trình vi phân dao động ma trận. Tác giả lập chương trình tính toán số, kiểm tra độ tin cậy của chương trình thông qua so sánh tần số dao động riêng và đáp ứng động với các công bố của tác giả khác. Chương này tiến hành tính toán số, phân tích đáp ứng động (chuyển vị, ứng suất, góc xoay) cho panel trụ và panel cầu composite lớp dưới các mức nhiệt độ và đặc trưng tải trọng sóng xung kích khác nhau.

Chương 3: TÍNH TOÁN VỎ THOẢI COMPOSITE SANDWICH CÓ CHỐT LIÊN KẾT CHỊU TÁC DỤNG ĐỒNG THỜI TẢI TRỌNG SÓNG XUNG KÍCH VÀ NHIỆT ĐỘ

Chương 3 mở rộng mô hình cho kết cấu vỏ thoải composite sandwich 3 phân lớp (lớp trên $t$, lớp lõi $c$, lớp dưới $b$) được liên kết với nhau bằng hệ thống chốt đàn hồi có độ cứng chống cắt $k_s$. Trường chuyển vị được thiết lập riêng cho từng lớp, cho phép tính toán chuyển vị trượt tương đối tại mặt tiếp xúc $z = z_t$ và $z = z_b$. Biểu thức thế năng biến dạng được bổ sung phần thế năng trượt của chốt $U_e^{sc}$. Tác giả xây dựng phần tử hữu hạn 8 nút 3 lớp đặc thù, tích hợp đầy đủ các ma trận khối lượng, ma trận độ cứng kết cấu, độ cứng chốt và véc tơ tải trọng nút do sóng xung kích và nhiệt độ. Tác giả kiểm tra sự hội tụ lưới, kiểm tra độ tin cậy của chương trình và thực hiện tính toán số đối với panel trụ composite 3 lớp có liên kết chốt.

Chương 4: KHẢO SÁT MỘT SỐ YẾU TỐ ẢNH HƯỞNG ĐẾN ĐÁP ỨNG ĐỘNG VỎ COMPOSITE SANDWICH CÓ LIÊN KẾT CHỐT

Chương 4 tập trung phân tích số toàn diện các yếu tố ảnh hưởng đến tần số dao động riêng và đáp ứng động lực học của vỏ thoải composite sandwich liên kết chốt (đặc biệt là panel cầu sandwich điều kiện biên tựa đơn - ngàm - tựa đơn S-C-S) chịu sóng xung kích trong môi trường nhiệt độ:

  • Khảo sát ảnh hưởng của các tham số hình học: Tỉ số chiều rộng trên chiều dày ($a/h$), tỉ số chiều dày lớp lõi trên lớp mặt ($h_c/h_t$), tỉ số cạnh mặt bằng ($a/b$), và tỉ số bán kính cong trên kích thước mặt bằng ($R/a$).
  • Khảo sát ảnh hưởng của các tham số vật liệu và liên kết: Độ cứng chống cắt của chốt liên kết ($k_s/E_t$), tỉ số mô đun đàn hồi giữa các lớp ($E_c/E_t$ và $E_b/E_t$).
  • Các bảng số liệu và đồ thị trong chương phản ánh rõ quy luật biến thiên của chuyển vị võng cực đại $w_{max}$, chuyển vị trượt tương đối lớn nhất giữa các lớp $V^{tc}{max}$, góc xoay $\phi_y$ và ứng suất lớn nhất $\sigma{max}$.

Dưới đây là bảng tổng hợp danh mục các bảng số liệu tính toán và kiểm chứng được thiết lập trong luận án:

Bảng số Tên nội dung bảng số liệu trong luận án
Bảng 2.1 Tần số dao động riêng đầu tiên không thứ nguyên
Bảng 2.2 Tám tần số dao động riêng đầu tiên
Bảng 3.1 Tần số dao động riêng không thứ nguyên
Bảng 3.2 Ảnh hưởng độ cứng chốt và chiều dày đến vị trí $w_{max}$ (mm)
Bảng 3.3 Ảnh hưởng độ cứng chốt và chiều dày đến trượt $V^{tc}_{max}$ (mm)
Bảng 3.4 Ảnh hưởng độ cứng chốt và chiều dày đến $\phi_y$ (rad)
Bảng 4.1 – 4.4 Tần số dao động riêng đầu tiên không thứ nguyên theo các tham số
Bảng 4.5 – 4.6 Ảnh hưởng tỉ số $a/h$ và $R/a$ đến chuyển vị và ứng suất lớn nhất
Bảng 4.7 – 4.8 Ảnh hưởng tỉ số $h_c/h_t$ và $R/a$ đến chuyển vị và ứng suất lớn nhất
Bảng 4.9 – 4.10 Ảnh hưởng tỉ số $a/b$ và $R/a$ đến chuyển vị và ứng suất lớn nhất
Bảng 4.11 – 4.12 Ảnh hưởng tỉ số $k_s/E_t$ và $R/a$ đến chuyển vị và ứng suất lớn nhất
Bảng 4.13 – 4.14 Ảnh hưởng tỉ số $E_c/E_t$ và $R/a$ đến chuyển vị và ứng suất lớn nhất
Bảng 4.15 – 4.16 Ảnh hưởng tỉ số $E_b/E_t$ và $R/a$ đến chuyển vị và ứng suất lớn nhất

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

Các kết quả và đóng góp mới của luận án được rút ra từ nội dung nghiên cứu gồm:

  1. Về mặt lý thuyết và phương pháp luận:

    • Xây dựng thành công mô hình toán học và hệ phương trình vi phân dao động hoàn chỉnh cho vỏ thoải composite lớp bám dính tuyệt đối và vỏ thoải composite sandwich có liên kết chốt chịu đồng thời tải trọng sóng xung kích (theo mô hình xung lượng nổ hóa học) và trường nhiệt độ dừng.
    • Đề xuất mô hình phần tử hữu hạn vỏ thoải 3 lớp 8 nút đặc thù, tích hợp độ cứng chống cắt của chốt liên kết, cho phép mô tả chính xác hiện tượng trượt đàn hồi tương đối giữa các phân lớp mà không cần sử dụng các giả thiết đơn lớp tương đương cứng bức.

  2. Về mặt thực tiễn và ứng dụng tính toán:

    • Xây dựng bộ chương trình tính toán số tin cậy, giải quyết bài toán động lực học phi dừng phức tạp cho kết cấu vỏ composite.
    • Làm sáng tỏ quy luật ứng xử động lực học của kết cấu: Nhiệt độ môi trường làm giảm độ cứng tổng thể của vỏ, dẫn đến giảm tần số dao động riêng và làm tăng biên độ chuyển vị động; độ cứng cắt của chốt $k_s$ có vai trò quyết định đến khả năng cản trượt và phân bố lại ứng suất giữa các lớp mặt và lớp lõi; việc lựa chọn tỉ số hình học ($a/h, h_c/h_t, R/a$) và tỉ số mô đun ($E_c/E_t, E_b/E_t$) hợp lý giúp giảm thiểu nguy cơ phá hủy cục bộ do ứng suất cắt và ứng suất uốn.
    • Cung cấp cơ sở khoa học và số liệu kỹ thuật hỗ trợ các kỹ sư trong việc tính toán, thiết kế, lựa chọn loại chốt liên kết và tối ưu hóa kết cấu vỏ composite sandwich chịu tải nổ trong môi trường nhiệt.

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

  • Hạn chế của nghiên cứu:

    • Bài toán được giới hạn trong phạm vi biến dạng đàn hồi tuyến tính của vật liệu, chưa xét đến biến dạng dẻo hoặc tổn thương tích lũy của chốt và vật liệu composite.
    • Các đặc trưng cơ học của vật liệu ($E, G, \nu$) được giả thiết không thay đổi theo nhiệt độ; trường nhiệt độ khảo sát là trường dừng phân bố đều.
    • Tải trọng sóng xung kích được mô hình hóa theo dạng xung lượng tam giác tương đương trong không khí vô hạn, chưa xét đến các hiệu ứng tương tác phức tạp của phản xạ nhiều lần hay tải trọng nổ trong môi trường ngập nước/đất đá đối với vỏ sandwich có chốt.

  • Hướng phát triển tiếp theo:

    • Mở rộng bài toán tính toán động lực học có xét đến yếu tố phi tuyến hình học (chuyển vị lớn) và phi tuyến vật liệu.
    • Nghiên cứu bài toán với trường nhiệt độ không dừng (biến thiên theo thời gian và không gian qua chiều dày vỏ) và cơ tính vật liệu phụ thuộc nhiệt độ.
    • Khảo sát kết cấu vỏ composite có gân gia cường, vỏ làm bằng vật liệu áp điện hoặc vật liệu có cơ tính biến thiên (FGM) kết hợp liên kết chốt.
    • Thực hiện các thử nghiệm nổ thực nghiệm trên mô hình vỏ sandwich liên kết chốt để kiểm chứng và hoàn thiện thêm các hệ số đàn - dẻo của chốt liên kết.

Giá trị tham khảo

Luận án là tài liệu tham khảo có giá trị chuyên môn cao đối với:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Cơ học vật rắn, Cơ kỹ thuật, Kỹ thuật xây dựng công trình ngầm/công sự, Kỹ thuật cơ khí và Hàng không - Vũ trụ: Cung cấp phương pháp thiết lập phương trình vi phân dao động vỏ thoải và kỹ thuật xây dựng phần tử hữu hạn nhiều lớp có xét trượt chốt liên kết (tham khảo chi tiết tại Chương 2 và Chương 3).
  • Các kỹ sư thiết kế kết cấu công trình quân sự và dân dụng: Tham khảo các biểu đồ đáp ứng động, số liệu khảo sát tham số hình học và độ cứng chốt liên kết tại Chương 4 để lựa chọn cấu hình vật liệu, tỷ lệ lớp sandwich và bố trí chốt liên kết tối ưu nhằm nâng cao khả năng kháng chấn và chịu tải trọng nổ của công trình.
  • Các nhà nghiên cứu phát triển phần mềm tính toán kết cấu: Tham khảo thuật toán phần tử hữu hạn và quy trình tích hợp tải trọng sóng xung kích kết hợp nhiệt độ để xây dựng các mô-đun phân tích chuyên dụng.

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

  1. Tại sao tác giả lại sử dụng lý thuyết vỏ thoải thay vì lý thuyết vỏ tổng quát hai độ cong?

    • Trả lời: Do kết cấu panel vỏ có độ vồng nhỏ được sử dụng rất phổ biến trong thực tế kỹ thuật. Lý thuyết vỏ thoải cho phép biểu diễn các quan hệ cơ học tương tự như kết cấu tấm nhưng vẫn giữ được độ cứng hình học do độ cong mang lại, giúp đơn giản hóa hệ phương trình vi phân dao động, tránh sự phức tạp của điều kiện Codaxi-Gauss và hệ phương trình cân bằng cồng kềnh mà vẫn đảm bảo độ chính xác tin cậy.

  2. Mô hình tải trọng sóng xung kích được tác giả thiết lập dựa trên cơ sở nào?

    • Trả lời: Tác giả sử dụng mô hình sóng xung kích do nổ hóa học trong không khí. Áp suất dư cực đại của sóng tới $\Delta P$ được tính theo công thức thực nghiệm Sa-đốp-sky, áp suất phản xạ cực đại $P_{px}$ được xác định với giả thiết góc tới $\alpha < 40^\circ$. Quy luật biến thiên áp suất theo thời gian được tuyến tính hóa thành biểu đồ tam giác cân bằng xung lượng với thời gian tác dụng hiệu dụng $\tau_2 = 2 I^+ / P_{px}$.

  3. Ưu điểm của giải pháp sử dụng chốt liên kết (shear connector) trong vỏ sandwich composite là gì?

    • Trả lời: Khi ghép nối các lớp composite có cơ tính khác biệt rõ rệt (như thép - bê tông - composite), hiện tượng tập trung ứng suất dễ gây bong tách lớp. Chốt liên kết cho phép các lớp có thể trượt đàn hồi tương đối với nhau mà kết cấu vẫn làm việc bình thường, khắc phục nguy cơ phá hủy liên kết dán, hỗ trợ lắp ghép linh hoạt và tăng khả năng tiêu tán năng lượng từ tải trọng nổ.

  4. Phần tử hữu hạn được sử dụng để mô hình hóa vỏ sandwich có chốt ở Chương 3 có đặc điểm gì?

    • Trả lời: Đó là mô hình phần tử vỏ thoải 8 nút gồm 3 phân lớp riêng biệt (lớp trên, lớp giữa, lớp dưới). Mỗi nút chứa các bậc tự do chuyển vị dài và góc xoay riêng cho từng lớp. Năng lượng biến dạng do trượt tương đối tại mặt tiếp xúc giữa các lớp thông qua chốt được tích hợp trực tiếp vào ma trận độ cứng của phần tử thông qua hệ số độ cứng cắt $k_s$.

  5. Nhiệt độ môi trường tác động như thế nào đến đáp ứng động của vỏ composite?

    • Trả lời: Nhiệt độ môi trường gây ra trường biến dạng nhiệt ban đầu và các lực màng nhiệt, mô men nhiệt tương đương trong kết cấu. Kết quả tính toán số chỉ ra rằng sự gia tăng nhiệt độ làm giảm tần số dao động riêng của vỏ và làm tăng biên độ dao động chuyển vị cũng như ứng suất động dưới tác dụng của sóng xung kích.

Kết luận

Luận án của nghiên cứu sinh Trần Ngọc Cảnh đã giải quyết trọn vẹn bài toán tính toán động lực học cho vỏ thoải composite lớp và vỏ thoải composite sandwich có liên kết chốt chịu tác dụng đồng thời của tải trọng sóng xung kích và nhiệt độ. Bằng việc kết hợp lý thuyết biến dạng cắt bậc nhất, mô hình tải nổ Sa-đốp-sky và phương pháp phần tử hữu hạn với phần tử 8 nút 3 lớp đặc thù, công trình đã mô tả chính xác hiện tượng trượt đàn hồi tại chốt và ứng xử cơ nhiệt của kết cấu. Các kết quả khảo sát số toàn diện trong luận án cung cấp cơ sở khoa học và định lượng quan trọng cho việc thiết kế, chế tạo và ứng dụng kết cấu composite sandwich chịu tải trọng đặc biệt trong thực tiễn.