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 có cơ tính biến thiên (Functionally Graded Material - FGM) là loại vật liệu composite tiên tiến được hình thành từ hai hoặc nhiều thành phần vật liệu (chủ yếu là gốm và kim loại), trong đó tỉ lệ thể tích của các pha thành phần biến đổi liên tục và trơn từ bề mặt này sang bề mặt kia theo chiều dày kết cấu. Nhờ tính chất kháng nhiệt của gốm kết hợp với độ dẻo dai của kim loại, FGM khắc phục được hiện tượng tập trung ứng suất và tách lớp thường gặp ở composite phân lớp truyền thống khi làm việc trong môi trường nhiệt độ cao hoặc chịu truyền nhiệt lớn như trong công nghiệp hàng không vũ trụ, tên lửa, tàu thủy, lò phản ứng hạt nhân và kết cấu cơ nhiệt.

Kết cấu vỏ hai độ cong được ứng dụng phổ biến trong thực tế công kỹ thuật. Tuy nhiên, các nghiên cứu chuyên sâu về ổn định tĩnh và đáp ứng động lực phi tuyến của vỏ thoải hai độ cong FGM vẫn còn hạn chế, đặc biệt trong các trường hợp kết cấu có xét đến tính không hoàn hảo hình dáng ban đầu, sự phụ thuộc của cơ tính vật liệu vào nhiệt độ, có bố trí gân gia cường hoặc tích hợp các lớp áp điện khi đặt trên nền đàn hồi chịu tác dụng của các trường tải trọng cơ học, nhiệt và cơ - nhiệt đồng thời. Luận án của tác giả Trần Quốc Quân (2018) với đề tài "Phân tích ổn định tĩnh và đáp ứng động lực phi tuyến của vỏ thoải hai độ cong FGM trên nền đàn hồi" được thực hiện nhằm giải quyết các khoảng trống khoa học này.

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

Mục tiêu tổng quát của luận án là thiết lập mô hình toán học giải tích để giải quyết bài toán ổn định tĩnh và động lực phi tuyến của vỏ thoải hai độ cong FGM và vỏ thoải hai độ cong FGM áp điện trên nền đàn hồi dưới tác dụng của các điều kiện tải cơ và nhiệt; xác định tải trọng tới hạn, tần số dao động tự do, mối liên hệ biên độ – tần số tuyến tính và phi tuyến, đường cong tải trọng – biên độ độ võng và đường cong biên độ độ võng – thời gian; đồng thời đánh giá định lượng ảnh hưởng của các tham số vật liệu, hình học, gân gia cường, nền đàn hồi và nhiệt độ lên đáp ứng của kết cấu.

Các nhiệm vụ nghiên cứu cụ thể được phân định theo văn bản gồm bốn nội dung:

  1. Phân tích ổn định tĩnh và đáp ứng động lực phi tuyến của vỏ thoải hai độ cong FGM không gân gia cường trên nền đàn hồi với các cạnh tựa bản lề chịu tải cơ học và tải cơ – nhiệt kết hợp bằng lý thuyết vỏ cổ điển, có xét đến tính không hoàn hảo hình dáng ban đầu và sự phụ thuộc nhiệt độ của tính chất vật liệu.
  2. Phân tích ổn định tĩnh và đáp ứng động lực phi tuyến của vỏ thoải hai độ cong FGM không hoàn hảo có gân gia cường với các cạnh tựa bản lề chịu các tải cơ, nhiệt và cơ – nhiệt kết hợp bằng lý thuyết vỏ cổ điển, có xét đến tính không hoàn hảo hình dáng và sự phụ thuộc nhiệt độ.
  3. Phân tích ổn định tĩnh và động lực phi tuyến của vỏ thoải hai độ cong FGM không hoàn hảo có gân gia cường chịu tác động của tải cơ học, tải nhiệt và tải cơ – nhiệt kết hợp bằng lý thuyết vỏ biến dạng trượt bậc ba.
  4. Phân tích đáp ứng động lực phi tuyến của vỏ thoải hai độ cong FGM không hoàn hảo không gân gia cường được tích hợp hai lớp áp điện ở mặt trên và mặt dưới sử dụng lý thuyết vỏ biến dạng trượt bậc ba trong trường hợp chịu tác dụng tổng hợp của tải cơ, nhiệt và điện.

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

  • Đối tượng nghiên cứu: Vỏ thoải hai độ cong làm bằng vật liệu composite FGM dạng hàm lũy thừa (P-FGM), có và không có gân gia cường lệch tâm (gân thuần nhất, đẳng hướng), vỏ không hoàn hảo hình dáng ban đầu, tựa trên nền đàn hồi hai tham số (Winkler – Pasternak), bao gồm cả kết cấu vỏ mỏng và vỏ dày, có xét tới trường hợp tích hợp lớp áp điện.
  • Phạm vi nghiên cứu: Ứng xử cơ học phi tuyến trong miền đàn hồi (không xét phá hủy dẻo), ổn định tĩnh (mất ổn định kiểu cực trị và rẽ nhánh) và đáp ứng động lực phi tuyến dưới tác dụng của tải nén dọc trục, áp lực ngoài phân bố đều, tải nhiệt đồng đều và tải kết hợp cơ – nhiệt – điện; điều kiện biên tựa bản lề (tựa tự do và tựa cố định trong mặt phẳng).

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

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 cơ học kết cấu FGM theo bốn nhóm chủ đề chính:

Các nghiên cứu về kết cấu FGM không có gân gia cường

  • Nhóm nghiên cứu quốc tế:
    • Shen và cộng sự (2003–2012): Shen & Wang (2012) phân tích uốn và sau vồng của panel trụ FGM trên nền đàn hồi; Shen (2009) nghiên cứu nhiệt sau vồng của tấm FGM; Shen et al. (2010) nghiên cứu sau vồng vỏ trụ FGM trên nền Pasternak chịu áp lực trong môi trường nhiệt độ; Huang & Shen (2004) công bố kết quả dao động phi tuyến tấm FGM truyền nhiệt; Yang & Shen (2003) nghiên cứu mất ổn định động của panel trụ FGM dưới tải dọc trục tĩnh và tuần hoàn dựa trên lý thuyết biến dạng trượt bậc cao của Reddy với cơ tính phụ thuộc nhiệt độ.
    • Sofiyev và cộng sự (2003–2015): Tập trung vào phân tích ổn định phi tuyến của vỏ trụ tròn, vỏ nón và vỏ nón cụt FGM chịu tải dọc trục, áp lực ngoài, áp suất thủy tĩnh trên nền Winkler – Pasternak dựa trên lý thuyết Donnell, Von Kármán và phương pháp Galerkin (Sofiyev 2003, 2007, 2010; Sofiyev & Kuruoglu 2007, 2011, 2015).
    • Huang, Han và cộng sự (2006–2011): Huang & Han (2008, 2009, 2010), Huang et al. (2011), Zhang et al. (2006) nghiên cứu vồng, sau vồng và ổn định động của vỏ trụ FGM đàn hồi và đàn dẻo chịu tải nén dọc trục, tải uốn thuần túy và áp lực ngoài bằng phương pháp năng lượng Ritz và lý thuyết Donnell.
    • Alijani, Amabili và cộng sự (2011–2014): Alijani & Amabili (2013, 2014), Alijani et al. (2011) khảo sát dao động cưỡng bức, dao động phi tuyến hình học và mất ổn định động của tấm chữ nhật và vỏ hai độ cong FGM chịu kích động điều hòa và biến đổi nhiệt độ.
  • Nhóm nghiên cứu trong nước:
    • Nguyễn Đình Đức và cộng sự (2012–2016): Khảo sát đáp ứng phi tuyến của vỏ cầu thoải, vỏ trụ tròn và tấm dày S-FGM, P-FGM, E-FGM trên nền Pasternak chịu tải cơ – nhiệt bằng lý thuyết cổ điển và biến dạng trượt bậc cao của Reddy (Nguyễn Đình Đức et al. 2013, 2015; Nguyễn Đình Đức & Phạm Hồng Công 2013).
    • Đào Huy Bích, Đào Văn Dũng và cộng sự (2011–2013): Đào Huy Bích et al. (2011, 2013) nghiên cứu đáp ứng tĩnh, động và dao động phi tuyến của vỏ cầu thoải, vỏ trụ tròn và mảnh vỏ trống FGM trên nền đàn hồi có xét ảnh hưởng nhiệt độ theo lý thuyết Donnell cải tiến và nghiệm ba số hạng.
    • Nguyễn Xuân Hùng và cộng sự (2013–2014); Nguyễn Thời Trung và cộng sự (2013–2014): Phát triển các phương pháp phần tử hữu hạn đẳng hình học (IGA) và phần tử hữu hạn trơn (CS-FEM, ES-FEM) để phân tích tĩnh, dao động tự do và mất ổn định tấm FGM.

Các nghiên cứu về kết cấu FGM có gân gia cường

  • Nghiên cứu quốc tế: Satouri et al. (2014) khảo sát vỏ trụ FGM gia cường gân dọc; Golmakani (2013) nghiên cứu tấm tròn FGM gân vòng chuyển vị lớn; Farahani et al. (2012) phân tích giải tích vỏ trụ FGM gia cường gân ngoài; Sun et al. (2014) khảo sát vỏ trụ FGM gia cường gân lưới chịu tải nén – nhiệt theo lý thuyết Reddy; Ardestani et al. (2014) và Taczała et al. (2013) nghiên cứu tấm FGM có gân đồng tâm/lệch tâm bằng FSDT và FEM.
  • Nghiên cứu trong nước: Đào Huy Bích et al. (2013, 2014) phân tích ổn định động panel trụ và vỏ trụ FGM có gân theo lý thuyết cổ điển và kỹ thuật san đều gân Lekhnitskii; Đào Văn Dũng & Lê Khả Hòa (2015) đề xuất nghiệm ba số hạng cho vỏ trụ FGM có gân chịu xoắn; Đào Văn Dũng & Đỗ Quang Chấn (2016) khảo sát vỏ nón cụt FGM gân trực giao theo FSDT; Nguyễn Đình Đức và cộng sự (2014–2016) nghiên cứu tấm, vỏ trụ ellip và vỏ nón cụt FGM có gân gia cường lệch tâm trong môi trường nhiệt độ, có kể đến biến dạng nhiệt trong gân.

Các nghiên cứu về kết cấu FGM áp điện và vỏ hai độ cong

  • Kết cấu FGM áp điện: Shen (2005), Dai et al. (2007), Ebrahimi et al. (2008, 2009), Shariyat (2009), Hashemi et al. (2010), Xia & Shen (2008), Sheng & Wang (2009), Rafiee et al. (2013), Kiani et al. (2011) cùng các nghiên cứu trong nước của Trần Ích Thịnh et al. (2006, 2011), Nguyễn Đình Đức et al. (2014, 2015) đã thiết lập các mô hình tính toán dao động, điều khiển chủ động và đáp ứng động học của tấm, vỏ trụ FGM tích hợp các lớp áp điện chịu tải nhiệt – điện – cơ.
  • Vỏ hai độ cong FGM: Fadaee et al. (2013) đưa ra nghiệm Levy cho dao động tự do; Đào Văn Dũng et al. (2013) phân tích sau vồng vỏ thoải hai độ cong FGM có gân theo lý thuyết Reddy; Najafi et al. (2011) khảo sát va chạm tốc độ thấp; Sayyaadi et al. (2015) nghiên cứu vỏ hai độ cong FGM áp điện; Shen et al. (2014) nghiên cứu dao động biên độ lớn theo HSDT; Đào Huy Bích et al. (2013) khảo sát động lực phi tuyến vỏ hai độ cong FGM có gân theo lý thuyết vỏ cổ điển.

Vị trí và khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết bài toán phân tích ổn định tĩnh và đáp ứng động lực phi tuyến cho vỏ thoải hai độ cong FGM dưới dạng nghiệm giải tích tường minh, bao quát đồng thời các yếu tố: độ không hoàn hảo hình dáng ban đầu, sự phụ thuộc nhiệt độ của vật liệu (T-D), hệ gân gia cường chịu ứng suất nhiệt, lớp phủ áp điện, mô hình nền hai tham số (Winkler – Pasternak) dưới cả hai khung lý thuyết vỏ cổ điển (CST) và lý thuyết biến dạng trượt bậc ba (TSDT).

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

Mô hình vật liệu FGM và sự phụ thuộc nhiệt độ

Luận án nghiên cứu mô hình P-FGM cấu tạo từ gốm ($c$) và kim loại ($m$). Tỉ lệ thể tích thành phần thay đổi theo chiều dày $z \in [-h/2, h/2]$ theo hàm lũy thừa: $$V_m(z) = \left(\frac{2z + h}{2h}\right)^N, \quad V_c(z) = 1 - V_m(z)$$ trong đó $N \ge 0$ là hệ số tỷ lệ thể tích. Khi $N = 0$, kết cấu thuần kim loại; khi $N \to \infty$, kết cấu thuần gốm.

Các đặc tính cơ lý hiệu dụng $Pr_{eff}(z)$ (mô đun đàn hồi $E$, hệ số Poisson $\nu$, hệ số giãn nở nhiệt $\alpha$, mật độ khối lượng $\rho$) được xác định theo quy tắc hỗn hợp: $$Pr_{eff}(z) = Pr_c V_c(z) + Pr_m V_m(z) = (Pr_m - Pr_c)\left(\frac{2z + h}{2h}\right)^N + Pr_c$$

Sự phụ thuộc nhiệt độ của các tính chất vật liệu thành phần được mô hình hóa theo hàm phi tuyến: $$Pr(T) = P_0 (1 + P_{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$$ với $T = T_0 + \Delta T$, $T_0 = 300\text{ K}$ (nhiệt độ phòng), $\Delta T$ là gia số nhiệt độ, và $P_0, P_{-1}, P_1, P_2, P_3$ là các hệ số thực nghiệm đặc trưng cho từng loại vật liệu.

Vật liệu Tính chất $P_0$ $P_{-1}\text{ (K)}$ $P_1\text{ (K}^{-1}\text{)}$ $P_2\text{ (K}^{-2}\text{)}$ $P_3\text{ (K}^{-3}\text{)}$
Gốm $\text{Si}_3\text{N}_4$ $E\text{ (Pa)}$ $348{,}43 \times 10^9$ $0$ $-3{,}070 \times 10^{-4}$ $2{,}160 \times 10^{-7}$ $-8{,}946 \times 10^{-11}$
$\alpha\text{ (K}^{-1}\text{)}$ $5{,}8723 \times 10^{-6}$ $0$ $9{,}095 \times 10^{-4}$ $0$ $0$
$\rho\text{ (kg/m}^3\text{)}$ $2370$ $0$ $0$ $0$ $0$
$\nu$ $0{,}24$ $0$ $0$ $0$ $0$
Kim loại SUS304 $E\text{ (Pa)}$ $201{,}04 \times 10^9$ $0$ $3{,}079 \times 10^{-4}$ $-6{,}534 \times 10^{-7}$ $0$
$\alpha\text{ (K}^{-1}\text{)}$ $12{,}330 \times 10^{-6}$ $0$ $8{,}086 \times 10^{-4}$ $0$ $0$
$\rho\text{ (kg/m}^3\text{)}$ $8166$ $0$ $0$ $0$ $0$
$\nu$ $0{,}3$ $0$ $0$ $0$ $0$

Mô hình nền đàn hồi và gân gia cường

  • Nền đàn hồi Pasternak: Phản lực tương tác giữa nền và đáy vỏ được biểu diễn: $$q_e = k_1 w - k_2 \nabla^2 w = k_1 w - k_2 \left(\frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2}\right)$$ trong đó $k_1$ là hệ số độ cứng nền Winkler và $k_2$ là hệ số lớp trượt Pasternak.
  • Gân gia cường: Sử dụng phương pháp san đều tác dụng gân Lekhnitskii đối với lý thuyết vỏ cổ điển và công thức tác dụng gân mở rộng có xét biến dạng nhiệt và ứng suất nhiệt trong gân đối với lý thuyết biến dạng trượt bậc cao. Gân được đặt lệch tâm, làm bằng vật liệu đồng nhất tương ứng với mặt gắn gân.

Mô hình lý thuyết vỏ và phương pháp giải

  1. Lý thuyết vỏ cổ điển (CST - Donnell): Dựa trên các giả thiết Kirchhoff - Love, quan hệ phi tuyến hình học Von Kármán, giả thiết Volmir đối với chuyển vị mặt phẳng ($\partial^2 u/\partial t^2 \approx 0, \partial^2 v/\partial t^2 \approx 0$), đưa hệ phương trình chuyển động và tương thích biến dạng về dạng phương trình vi phân phi tuyến theo hàm độ võng $w(x,y,t)$ và hàm ứng suất $f(x,y,t)$. Độ lệch hình dáng ban đầu được mô hình hóa qua hàm $w^*(x,y) = \mu h \sin \alpha_m x \sin \beta_n y$, với $\mu \in [0, 1]$ là hệ số biên độ không hoàn hảo.
  2. Lý thuyết vỏ biến dạng trượt bậc ba (TSDT - Reddy): Loại bỏ giả thiết pháp tuyến thẳng, xét đến sự biến thiên phi tuyến bậc ba của biến dạng trượt ngang qua chiều dày kết cấu, thỏa mãn điều kiện ứng suất trượt bằng không tại hai bề mặt ngoài mà không cần hệ số hiệu chỉnh trượt.
  3. Kỹ thuật tính toán:
    • Sử dụng dạng nghiệm xấp xỉ thỏa mãn điều kiện biên tựa bản lề: $$w(x, y, t) = W(t) \sin \alpha_m x \sin \beta_n y$$ trong đó $\alpha_m = m\pi/a, \beta_n = n\pi/b$, $W(t)$ là biên độ độ võng.
    • Áp dụng phương pháp Galerkin để biến đổi hệ phương trình vi phân đạo hàm riêng thành phương trình vi phân thường phi tuyến theo biến thời gian $W(t)$.
    • Giải bài toán ổn định tĩnh bằng việc triệt tiêu số hạng quán tính, xác định cực trị giải tích $dq/dW = 0$.
    • Giải bài toán động lực phi tuyến bằng phương pháp tích phân số Runge – Kutta bậc 4.
  4. Tiêu chuẩn xác định mất ổn định:
    • Ổn định tĩnh: Tiêu chuẩn tĩnh học dựa trên đường cong cân bằng lệch và các điểm cực trị ($q_{upper}, q_{lower}$).
    • Ổn định động: Tiêu chuẩn Budiansky – Roth xác định thời điểm tới hạn $t_{cr}$ tại vị trí mà đường cong biên độ độ võng – thời gian $W(t)$ đạt giá trị cực đại đầu tiên sau khi có sự tăng vọt độ dốc.

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

CHƯƠNG 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU

Chương 1 trình bày hệ thống cơ sở lý luận về vật liệu FGM, phương pháp chế tạo, phân loại (P-FGM, S-FGM, E-FGM), và phạm vi ứng dụng trong các ngành công nghệ cao. Tác giả phân tích hai dạng mất ổn định cơ bản trong cơ học kết cấu:

  • Mất ổn định loại I (rẽ nhánh): Xuất hiện điểm phân nhánh trên đường cong cân bằng, phát sinh trạng thái cân bằng mới khác trạng thái ban đầu.
  • Mất ổn định loại II (cực trị): Đường cong tải trọng – độ võng không phân nhánh, kết cấu võng liên tục từ khi đặt tải và mất ổn định khi đạt tải cực trị trên $q_{upper}$, sau đó chuyển sang trạng thái sau vồng với tải cực trị dưới $q_{lower}$.

Chương này cũng tổng quan toàn diện các công trình nghiên cứu trong nước và quốc tế về tấm, panel trụ, vỏ trụ, vỏ nón, vỏ nón cụt và vỏ hai độ cong FGM; phân tích cơ chế làm việc của vật liệu áp điện (biến đổi năng lượng điện – cơ); đồng thời xác định cụ thể 4 mục tiêu nghiên cứu của đề tài.

CHƯƠNG 2: PHÂN TÍCH ỔN ĐỊNH TĨNH VÀ ĐÁP ỨNG ĐỘNG LỰC PHI TUYẾN CỦA VỎ THOẢI HAI ĐỘ CONG FGM KHÔNG HOÀN HẢO KHÔNG GÂN GIA CƯỜNG SỬ DỤNG LÝ THUYẾT VỎ CỔ ĐIỂN

Chương 2 thiết lập hệ phương trình cơ bản cho vỏ thoải hai độ cong FGM không gân gia cường, kích thước $a \times b$, chiều dày $h$, bán kính cong $R_x, R_y$, tựa trên nền đàn hồi Winkler – Pasternak dựa trên lý thuyết vỏ cổ điển Donnell và giả thiết Volmir. Hệ phương trình được đặt theo ẩn độ võng $w$ và hàm ứng suất $f$: $$\begin{aligned} P_1 \nabla^4 f + P_2 \nabla^4 w - \left(\frac{1}{R_x}\frac{\partial^2 f}{\partial y^2} + \frac{1}{R_y}\frac{\partial^2 f}{\partial x^2}\right) - \left[\frac{\partial^2 f}{\partial y^2}\frac{\partial^2(w+w^)}{\partial x^2} - 2\frac{\partial^2 f}{\partial x\partial y}\frac{\partial^2(w+w^)}{\partial x\partial y} + \frac{\partial^2 f}{\partial x^2}\frac{\partial^2(w+w^)}{\partial y^2}\right] \ = q - k_1 w + k_2 \nabla^2 w - I_1 \frac{\partial^2 w}{\partial t^2} \end{aligned}$$ và phương trình tương thích biến dạng: $$P_3 \nabla^4 w + P_4 \nabla^4 f + \left(\frac{1}{R_x}\frac{\partial^2 w}{\partial y^2} + \frac{1}{R_y}\frac{\partial^2 w}{\partial x^2}\right) = \left(\frac{\partial^2 w}{\partial x\partial y}\right)^2 - \frac{\partial^2 w}{\partial x^2}\frac{\partial^2 w}{\partial y^2} + 2\frac{\partial^2 w}{\partial x\partial y}\frac{\partial^2 w^}{\partial x\partial y} - \frac{\partial^2 w}{\partial x^2}\frac{\partial^2 w^}{\partial y^2} - \frac{\partial^2 w^}{\partial x^2}\frac{\partial^2 w}{\partial y^2}$$

Sử dụng phương pháp Galerkin, tác giả dẫn ra phương trình quan hệ tải trọng – biên độ độ võng giải tích cho vỏ chịu áp lực ngoài $q$: $$q = b_{11} \bar{W} + b_{21} \bar{W}(\bar{W} + \mu) + b_{31} \bar{W}(\bar{W} + 2\mu) + b_{41} \bar{W}(\bar{W} + \mu)(\bar{W} + 2\mu)$$ với $\bar{W} = W/h$. Đối với vỏ hoàn hảo ($\mu = 0$), biểu thức tải tới hạn cực đại $q_{upper}$ được dẫn ra dưới dạng hiển.

Các phát hiện và kết quả số chính:

  • Khi hệ số tỷ lệ thể tích $N$ tăng, hàm lượng kim loại tăng làm độ cứng của vỏ giảm, dẫn đến độ võng tăng và tải trọng tới hạn tĩnh/động giảm.
  • Nền đàn hồi (đặc biệt là hệ số trượt Pasternak $k_2$) và độ cong của vỏ giúp tăng độ cứng tổng thể, nâng cao tải tới hạn và giảm biên độ dao động.
  • Tính không hoàn hảo hình dáng ban đầu ($\mu$) làm suy giảm đáng kể tải tới hạn trên và làm thay đổi dạng đường cong tải trọng – độ võng (biến mất điểm rẽ nhánh rõ rệt, chuyển hẳn sang đáp ứng uốn phi tuyến).
  • Gia số nhiệt độ $\Delta T$ và tính chất phụ thuộc nhiệt độ (T-D) làm mềm kết cấu, giảm tần số dao động tự nhiên và giảm tải tới hạn so với trường hợp độc lập nhiệt độ (T-ID).

CHƯƠNG 3: PHÂN TÍCH ỔN ĐỊNH TĨNH VÀ ĐỘNG LỰC PHI TUYẾN CỦA VỎ THOẢI HAI ĐỘ CONG FGM KHÔNG HOÀN HẢO CÓ GÂN GIA CƯỜNG SỬ DỤNG LÝ THUYẾT VỎ CỔ ĐIỂN

Chương 3 mở rộng bài toán cho vỏ thoải hai độ cong FGM được gia cường bởi hệ gân trực giao dọc và ngang đặt lệch tâm, sử dụng kỹ thuật san đều tác dụng gân của Lekhnitskii.

Các luận điểm và phát hiện chính:

  • Bố trí hệ gân gia cường giúp tăng cường độ cứng kết cấu, làm tăng đáng kể tải tới hạn tĩnh và động, đồng thời giảm mạnh biên độ độ võng trong đáp ứng thời gian mà chỉ làm tăng một lượng nhỏ khối lượng tổng thể.
  • Luận án khảo sát định lượng ảnh hưởng của các điều kiện biên trong mặt phẳng: trường hợp các cạnh tựa cố định ($u=0, v=0$ tại biên) cho giá trị tải tới hạn và tần số dao động tự nhiên cao hơn rõ rệt so với trường hợp các cạnh tựa tự do di trượt ($N_{x0}=0, N_{y0}=0$).
  • Khảo sát tải nhiệt thuần túy và tải cơ – nhiệt kết hợp: khi nhiệt độ môi trường tăng, hiện tượng mất ổn định nhiệt xảy ra do ứng suất nén nhiệt sinh ra trong vỏ và hệ gân; việc tăng kích thước gân làm tăng nhiệt độ tới hạn của kết cấu.
  • Tốc độ đặt tải động $s$ ($q = s \cdot t$) có ảnh hưởng lớn đến giá trị tải tới hạn động theo tiêu chuẩn Budiansky – Roth: tốc độ đặt tải càng lớn thì tải tới hạn động nhận được càng cao.

CHƯƠNG 4: PHÂN TÍCH ỔN ĐỊNH TĨNH VÀ ĐÁP ỨNG ĐỘNG LỰC PHI TUYẾN CỦA VỎ THOẢI FGM HAI ĐỘ CONG KHÔNG HOÀN HẢO SỬ DỤNG LÝ THUYẾT VỎ BIẾN DẠNG TRƯỢT BẬC CAO

Chương 4 nâng cấp mô hình phân tích sang lý thuyết biến dạng trượt bậc ba của Reddy (TSDT) để áp dụng cho cả kết cấu vỏ mỏng và vỏ dày, gồm hai nội dung chính:

  1. Vỏ thoải hai độ cong FGM có gân gia cường: Sử dụng mô hình tác dụng gân cải tiến có xét đến ứng suất nhiệt trong gân. Kết quả cho thấy lý thuyết cổ điển đánh giá cao hơn độ cứng của kết cấu so với lý thuyết bậc cao do bỏ qua biến dạng trượt ngang; sai số giữa hai lý thuyết càng tăng khi tỷ số chiều dày trên nhịp ($h/b$) tăng.
  2. Vỏ thoải hai độ cong FGM áp điện: Cấu trúc gồm lớp lõi FGM và hai lớp áp điện phủ ở mặt trên và mặt dưới. Điện trường ngoài đặt vào các lớp áp điện tạo ra hiệu ứng biến dạng ngược, sinh ra mô men và lực nội tại trong mặt phẳng.
    • Khi đặt hiệu điện thế âm ($V_a < 0$), kết cấu bị nén thêm làm giảm tần số dao động và giảm khả năng chịu tải.
    • Khi đặt hiệu điện thế dương ($V_a > 0$), lớp áp điện tạo lực kéo căng giúp tăng độ cứng hiệu dụng, tăng tần số dao động tự do và triệt tiêu bớt biên độ dao động cưỡng bức, chứng minh khả năng điều khiển chủ động đáp ứng động lực của kết cấu thông qua điện thế điều khiển.

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

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

  1. Thiết lập thành công hệ phương trình cơ bản và dẫn ra các công thức giải tích dạng tường minh xác định tải tới hạn tĩnh, tải tới hạn động, tần số dao động tự do và quan hệ phi tuyến tải trọng – độ võng, biên độ – tần số cho vỏ thoải hai độ cong FGM và vỏ FGM áp điện trên nền đàn hồi.
  2. Tích hợp đồng thời trong mô hình giải tích các yếu tố phi tuyến hình học Von Kármán, tính không hoàn hảo hình dáng ban đầu, cơ tính phụ thuộc nhiệt độ (T-D), hệ gân gia cường lệch tâm có xét ứng suất nhiệt, và tương tác nền đàn hồi hai tham số (Winkler – Pasternak).
  3. Đánh giá chính xác phạm vi áp dụng và sự khác biệt giữa lý thuyết vỏ cổ điển (CST) và lý thuyết biến dạng trượt bậc ba (TSDT) đối với kết cấu vỏ thoải FGM khi thay đổi độ dày và độ cong.

Đóng góp mới về mặt thực tiễn

  1. Cung cấp bộ công cụ tính toán giải tích thuận tiện, chính xác cho các kỹ sư kết cấu, giúp dự báo nhanh khả năng chịu tải và độ ổn định của vỏ thoải hai độ cong FGM mà không cần thiết lập các mô hình phần tử hữu hạn 3D phức tạp.
  2. Làm sáng tỏ quy luật ảnh hưởng của các tham số thiết kế: hệ số phân bố thể tích $N$, tỷ lệ kích thước hình học ($b/h, b/a, R_x/h$), mật độ gân gia cường, độ cứng nền đàn hồi ($k_1, k_2$), gia số nhiệt độ ($\Delta T$) và hiệu điện thế điều khiển áp điện ($V_a$).

Kiến nghị và giải pháp đề xuất

  • Trong tính toán thiết kế kết cấu FGM làm việc ở môi trường nhiệt độ cao, bắt buộc phải sử dụng mô hình cơ tính phụ thuộc nhiệt độ (T-D) để tránh đánh giá quá cao khả năng chịu tải của kết cấu.
  • Đối với vỏ có độ dày trung bình và vỏ dày ($b/h \le 20$), cần áp dụng lý thuyết biến dạng trượt bậc cao (TSDT) để đảm bảo độ chính xác.
  • Tận dụng hệ gân gia cường và lớp phủ áp điện điều khiển điện áp dương để tối ưu hóa khối lượng kết cấu, nâng cao tải trọng tới hạn và giảm dao động trong môi trường chịu tải phức tạp.

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

Hạn chế của luận án

  • Luận án chỉ giới hạn nghiên cứu quy luật phân bố vật liệu dạng hàm lũy thừa (P-FGM); chưa khảo sát các dạng phân bố hàm Sigmoid (S-FGM) hay hàm số mũ (E-FGM) cho vỏ hai độ cong trong cùng hệ thống công thức.
  • Giả thiết vật liệu hoàn toàn đàn hồi tuyến tính, chưa xét đến biến dạng dẻo hoặc hiện tượng rạn nứt, phá hủy vật liệu dưới tải trọng lớn.
  • Mô hình nền đàn hồi giả thiết tiếp xúc liên tục lý tưởng, chưa xét đến hiện tượng vỏ tách rời khỏi nền cục bộ.
  • Điều kiện biên chủ yếu tập trung vào liên kết tựa bản lề bốn cạnh (tựa tự do và tựa cố định), chưa mở rộng cho biên ngàm hoặc các biên hỗn hợp phức tạp.

Hướng nghiên cứu phát triển

  • Mở rộng phân tích ổn định tĩnh và động lực cho vỏ hai độ cong làm bằng vật liệu FGM có lỗ rỗng (porous FGM) hoặc vật liệu composite gia cường ống nano carbon (CNTRC).
  • Nghiên cứu bài toán tương tác chất lưu – kết cấu (FSI) và vỏ FGM chứa chất lỏng hoặc chịu tải trọng sóng nổ, va chạm tốc độ cao.
  • Phát triển mô hình tính toán ổn định phi tuyến có xét đến quá trình phát triển vết nứt, phân lớp cơ nhiệt và hiện tượng tiếp xúc phi tuyến với nền đàn hồi.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Kỹ thuật Cơ khí, Kỹ thuật Xây dựng và Công trình biển: Luận án là tài liệu tham khảo chi tiết về quy trình thiết lập phương trình vi phân phi tuyến giải tích, phương pháp giải Galerkin, phương pháp Runge – Kutta và kỹ thuật áp dụng tiêu chuẩn ổn định động Budiansky – Roth.
  • Đối với giảng viên và nhà nghiên cứu: Cung cấp hệ thống nghiệm giải tích chuẩn xác (benchmark analytical solutions) phục vụ việc kiểm chứng và so sánh độ tin cậy của các thuật toán số, phần mềm phần tử hữu hạn thương mại hoặc các phương pháp số mới (Meshfree, IGA, S-FEM).
  • Đối với kỹ sư thiết kế trong ngành hàng không vũ trụ, đóng tàu và công trình chịu nhiệt: Cung cấp cơ sở định lượng và các biểu đồ quan hệ tải trọng – biến dạng giúp tối ưu hóa việc lựa chọn tỷ lệ vật liệu gốm/kim loại, cấu hình gân gia cường và thông số kích động áp điện nhằm đảm bảo an toàn kết cấu.

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

1. Vật liệu FGM trong luận án được mô hình hóa theo quy luật nào và đặc trưng phụ thuộc nhiệt độ được biểu diễn ra sao?

Vật liệu trong luận án được mô hình hóa theo quy luật hàm lũy thừa (P-FGM), trong đó tỷ lệ thể tích kim loại $V_m(z) = [(2z+h)/(2h)]^N$ và gốm $V_c(z) = 1 - V_m(z)$ với $N$ là hệ số tỷ lệ thể tích ($0 \le N < \infty$). Các đặc tính cơ lý phụ thuộc nhiệt độ (T-D) của gốm ($\text{Si}_3\text{N}_4, \text{Al}_2\text{O}_3, \text{ZrO}2$) và kim loại ($\text{SUS304}, \text{Ti-6Al-4V}$) được biểu diễn qua hàm phi tuyến bậc ba theo nhiệt độ: $Pr(T) = P_0 (1 + P{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$ với các hệ số thực nghiệm $P_i$ xác định.

2. Luận án sử dụng những tiêu chuẩn nào để xác định tải tới hạn trong bài toán ổn định tĩnh và ổn định động?

  • Trong bài toán ổn định tĩnh, luận án sử dụng tiêu chuẩn tĩnh học (Static criterion), khảo sát trạng thái cân bằng lệch và tìm điểm cực trị của đường cong quan hệ tải trọng – độ võng thông qua điều kiện triệt tiêu đạo hàm $dq/dW = 0$ để xác định tải tới hạn trên $q_{upper}$ và tải tới hạn dưới $q_{lower}$.
  • Trong bài toán ổn định động, luận án áp dụng tiêu chuẩn ổn định động Budiansky – Roth, trong đó thời điểm tới hạn $t_{cr}$ được xác định tại điểm mà đường cong biên độ độ võng theo thời gian $W(t)$ đạt giá trị cực đại đầu tiên sau giai đoạn tăng vọt độ dốc, từ đó suy ra giá trị tải tới hạn động tương ứng.

3. Phương pháp giải tích nào được áp dụng để giải hệ phương trình vi phân phi tuyến của vỏ hai độ cong FGM?

Luận án sử dụng phương pháp đặt bài toán theo ứng suất với hai ẩn hàm là độ võng $w(x,y,t)$ và hàm ứng suất $f(x,y,t)$. Dạng nghiệm hàm lượng giác thỏa mãn điều kiện biên tựa bản lề được thay vào hệ phương trình đạo hàm riêng. Sau đó, phương pháp Bubnov – Galerkin được sử dụng để khử các biến không gian, đưa hệ phương trình về dạng phương trình vi phân thường phi tuyến theo biến thời gian. Để giải đáp ứng động lực theo thời gian, phương pháp tích phân số Runge – Kutta bậc 4 được áp dụng kết hợp.

4. Sự khác biệt cơ bản giữa việc phân tích vỏ hai độ cong FGM theo lý thuyết vỏ cổ điển (Chương 2, 3) và lý thuyết biến dạng trượt bậc cao (Chương 4) là gì?

Lý thuyết vỏ cổ điển (CST) dựa trên giả thiết Kirchhoff bỏ qua biến dạng trượt ngang ($\gamma_{xz} = \gamma_{yz} = 0$), do đó có xu hướng đánh giá độ cứng của vỏ cao hơn thực tế và chỉ phù hợp cho kết cấu vỏ mỏng. Lý thuyết biến dạng trượt bậc ba của Reddy (TSDT) ở Chương 4 xét đến sự phân bố phi tuyến parabol của ứng suất tiếp qua chiều dày, triệt tiêu ứng suất trượt tại hai bề mặt tự do mà không cần hệ số hiệu chỉnh trượt, cho phép mô tả chính xác ứng xử ổn định và dao động của cả vỏ mỏng và vỏ dày ($b/h \le 20$).

5. Hiệu ứng áp điện được tích hợp và có vai trò như thế nào trong phân tích động lực ở Chương 4?

Kết cấu được tích hợp hai lớp áp điện ở mặt trên và mặt dưới lớp lõi FGM. Khi đặt hiệu điện thế $V_a$ vào các lớp áp điện, hiệu ứng áp điện nghịch sinh ra trường lực và mô men trong mặt phẳng vỏ. Hiệu điện thế dương ($V_a > 0$) tạo ứng suất kéo giúp tăng độ cứng hiệu dụng của kết cấu, làm tăng tần số dao động tự do và giảm biên độ dao động cưỡng bức, đóng vai trò như một cơ chế điều khiển chủ động dao động và ổn định động lực của vỏ.

Kết luận

Luận án tiến sĩ của tác giả Trần Quốc Quân đã giải quyết bài toán phân tích ổn định tĩnh và đáp ứng động lực phi tuyến của vỏ thoải hai độ cong FGM không hoàn hảo trên nền đàn hồi chịu tải cơ, nhiệt và điện. Bằng việc kết hợp lý thuyết vỏ cổ điển Donnell và lý thuyết biến dạng trượt bậc ba Reddy với phương pháp giải tích Galerkin và Runge – Kutta, luận án đã dẫn ra hệ thống công thức nghiệm giải tích tường minh xác định tải tới hạn tĩnh, tải tới hạn động và các đường cong đáp ứng phi tuyến của kết cấu. Các kết quả nghiên cứu đã làm sáng tỏ ảnh hưởng tương hỗ của hệ số vật liệu, gân gia cường, nền đàn hồi, nhiệt độ và hiệu điện thế áp điện, cung cấp cơ sở khoa học tin cậy cho công tác tính toán và thiết kế kết cấu FGM trong thực tế kỹ thuật.