Tổng quan về luận án

Nghiên cứu "Phân tích ổn định tĩnh của vỏ bằng vật liệu có cơ tính biến thiên" của tác giả Lê Khả Hòa (chuyên ngành Cơ học vật rắn biến dạng, mã số 62440107, người hướng dẫn khoa học: PGS. TS. Đào Văn Dũng, Đại học Khoa học Tự nhiên - ĐHQGHN, 2015) giải quyết bài toán cơ học then chốt trong tính toán kết cấu vỏ mỏng chịu tải phức tạp. Vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM), được khởi xướng từ năm 1984 tại Sendai (Nhật Bản), là dòng vật liệu composite tiên tiến cấu thành từ gốm (ceramic) và kim loại (metal). Nhờ cơ tính biến đổi trơn và liên tục qua chiều dày, FGM loại bỏ triệt để hiện tượng tập trung ứng suất và bong tách phân lớp (delamination) vốn là nhược điểm chí mạng của composite xếp lớp truyền thống khi làm việc trong môi trường nhiệt độ siêu cao như vỏ tên lửa, khung máy bay siêu thanh và thành plasma của lò phản ứng hạt nhân nhiệt hạch.

Khoảng trống học thuật (research gap) mà luận án định vị là sự thiếu vắng các nghiệm giải tích tường minh (analytical closed-form solutions) cho bài toán ổn định phi tuyến của các kết cấu vỏ FGM phức tạp có xét đến tính biến thiên của hệ số Poisson $\nu = \nu(z)$, kết cấu có gân gia cường lệch tâm làm bằng vật liệu FGM (ES-FGM), vỏ nón cụt FGM và tương tác phức hợp với nền đàn hồi hai tham số Winkler - Pasternak. Phần lớn các công trình quốc tế trước đó (như Huang và Han, 2008, 2010; Shen, 2002) chỉ khảo sát kết cấu không gân, giả định hệ số Poisson là hằng số ($\nu = \text{const}$), hoặc để các hệ số độ cứng dưới dạng tích phân số thuần túy mà chưa thiết lập được dạng giải tích hiển.

Luận án thiết lập 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H):

  1. RQ1: Sự biến thiên của hệ số Poisson $\nu(z)$ theo chiều dày ảnh hưởng như thế nào đến tải tới hạn và đường cong tải - độ võng sau tới hạn của panel trụ FGM không hoàn hảo?
    • H1: Việc bỏ qua biến thiên của $\nu(z)$ sẽ dẫn đến sai số đáng kể trong đánh giá khả năng chịu lực tới hạn dưới tải nén dọc trục.
  2. RQ2: Quy luật ứng xử phi tuyến của vỏ trụ tròn mỏng FGM có gân gia cường FGM lệch tâm (ES-FGM) dưới tác dụng của áp lực ngoài và tải xoắn thuần túy diễn ra như thế nào khi biểu diễn hàm độ võng ba số hạng?
    • H2: Gân gia cường lệch tâm FGM làm thay đổi cơ chế mất ổn định từ rẽ nhánh sang cực trị và gia tăng vượt bậc tải vồng tới hạn trên ($q_{\text{upper}}$) và dưới ($q_{\text{lower}}$).
  3. RQ3: Tương tác cơ học giữa nền đàn hồi Winkler - Pasternak và vỏ trụ/vỏ nón cụt ES-FGM có triệt tiêu được các mode mất ổn định sóng dài hay không?
    • H3: Hệ số nền Pasternak $K_2$ đóng vai trò chi phối trong việc dịch chuyển mode vồng tới hạn $(m, n)$ về phía các bậc cao hơn và nâng cao độ cứng sau tới hạn.
  4. RQ4: Lời giải giải tích ổn định tuyến tính cho vỏ nón cụt FGM có gân gia cường trên nền đàn hồi có đạt độ chính xác tương thích khi suy biến về các trường hợp giới hạn của vỏ trụ và vật liệu thuần nhất?
    • H4: Phương pháp giải tích kết hợp kỹ thuật san đều gân Lekhnitskii và phương pháp Galerkin giải quyết được hệ phương trình vi phân đạo hàm riêng biến hệ số trên hệ tọa độ cong.

Khung lý thuyết chủ đạo tích hợp lý thuyết vỏ mỏng Donnell - von Kármán, kỹ thuật san đều gân cải tiến của Lekhnitskii, mô hình nền đàn hồi Winkler - Pasternak và phương pháp chiếu trực giao Galerkin. Phạm vi nghiên cứu bao quát các cặp tổ hợp vật liệu hàng không vũ trụ điển hình: Silicon nitride/Stainless steel ($Si_3N_4/SUS304$ với $E_c = 322.27\text{ GPa}, E_m = 207.4\text{ GPa}, \rho_m = 8166\text{ kg/m}^3, \rho_c = 2370\text{ kg/m}^3$), Zirconia/Titanium alloy ($ZrO_2/Ti\text{-}6Al\text{-}4V$), phân bố thể tích theo quy luật lũy thừa (P-FGM) và hàm mũ (E-FGM), tỷ số bán kính trên chiều dày $R/h$ từ 100 đến 300, tỷ số chiều dài trên bán kính $L/R$ từ 0.5 đến 3, và biên độ bất hoàn hảo ban đầu $\xi = 0.0 \div 0.5$.

Literature Review và Positioning

Tổng quan y văn cơ học kết cấu FGM ghi nhận ba mạch nghiên cứu chính. Mạch thứ nhất tập trung vào ổn định tấm và vỏ FGM không gân của Shen (2002, 2008), Sofiyev (2003, 2007), Huang và Han (2008, 2009, 2010). Huang và Han áp dụng lý thuyết biến dạng lớn và phương pháp Ritz để tính toán vỏ trụ không gân chịu nén trục, xoắn và áp lực ngoài; tuy nhiên, các tác giả này để hệ số độ cứng $A_{ij}$ dưới dạng tích phân số, gây hạn chế cho việc tối ưu hóa giải tích. Mạch thứ hai là các nghiên cứu về kết cấu FGM có gân gia cường của Đào Huy Bích và cộng sự (2011, 2013), Nguyễn Đình Đức (2011, 2014), Đào Văn Dũng và Vũ Hoài Nam (2013, 2014). Hầu hết các công trình này giả định gân làm bằng vật liệu thuần nhất kim loại hoặc gốm, chưa xét đến trường hợp gân có cơ tính biến thiên FGM đồng bộ với vỏ. Mạch thứ ba là nghiên cứu vỏ nón FGM của Singer (1970), Tong và Wang (1992), Civalek (2007), Sofiyev (2007, 2011) chủ yếu khảo sát kết cấu vỏ nón trơn không gia cường.

Trong y văn tồn tại hai cuộc tranh luận học thuật lớn:

  1. Tranh luận về giả định hệ số Poisson: Phần lớn các nhà nghiên cứu (như Shen, 2002; Nguyễn Đình Đức và Hoàng Văn Tùng, 2011) giả định $\nu = \text{const} \approx 0.3$ để đơn giản hóa quá trình khử ứng suất, cho rằng ảnh hưởng của $\nu(z)$ là thứ yếu. Phái đối lập, đại diện bởi Đào Văn Dũng và Nguyễn Thị Nga (2013), chứng minh rằng trong điều kiện gradient nhiệt và cơ học khắc nghiệt, sự biến thiên $\nu(z) = \nu_c V_c(z) + \nu_m V_m(z)$ làm thay đổi trường biến dạng màng và mômen uốn liên kết, khiến giá trị tải tới hạn rẽ nhánh có thể lệch từ 3% đến 8% so với mô hình $\nu = \text{const}$.
  2. Tranh luận về mô hình gân gia cường: Kỹ thuật rời rạc hóa gân (discrete stiffener model) phản ánh chính xác hiệu ứng cục bộ nhưng dẫn đến hệ phương trình vi phân phi tuyến bậc cao rất khó tìm nghiệm giải tích; ngược lại, kỹ thuật san đều gân Lekhnitskii (smeared stiffener technique) giả định phân bố liên tục độ cứng của gân vào bề mặt vỏ. Nhiều học giả quốc tế hoài nghi độ chính xác của kỹ thuật san gân khi số lượng gân nhỏ.

Luận án định vị chính xác vị trí khoa học bằng việc kết hợp giải tích hiển hệ số độ cứng $A_{ij}$ khi $\nu = \nu(z)$, mở rộng kỹ thuật san gân Lekhnitskii cho hệ gân FGM lệch tâm (ES-FGM), và sử dụng hàm độ võng ba số hạng để nắm bắt hiện tượng mất ổn định kiểu cực trị (snap-through). Luận án thực hiện so sánh đối chuẩn trực tiếp với hai nghiên cứu quốc tế kinh điển:

  • So sánh tải nén tới hạn không thứ nguyên $N_x / (ER)$ của panel trụ thuần nhất ($a/b = 0.5, b/h = 75, \nu = 0.3$) cho thấy kết quả của luận án đạt $0.60523 \times 10^{-4}$, trùng khớp hoàn toàn với kết quả của Turvey (1977) và sai lệch dưới 0.1% so với Shen (2002).
  • So sánh tải xoắn tới hạn $\tau_{cr}$ và áp lực ngoài tới hạn $q_{cr}$ với kết quả của Huang và Han (2008, 2009), Brush và Almroth (1975), Seide (1956), và Sofiyev (2007, 2011), khẳng định tính hội tụ và độ chính xác của phương pháp giải tích đề xuất.

Đó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 vỏ mỏng Donnell - von Kármán sang môi trường vật liệu có tính chất dị hướng liên tục theo chiều dày thông qua việc tích phân giải tích thành công các hệ số độ cứng màng, uốn và liên kết kéo-uốn $A_{ij}$: $$A_{1j} = \int_{-h/2}^{h/2} \frac{E(z)}{1 - \nu^2(z)} z^j dz, \quad A_{2j} = \int_{-h/2}^{h/2} \frac{E(z)\nu(z)}{1 - \nu^2(z)} z^j dz, \quad A_{3j} = \frac{1}{2}(A_{1j} - A_{2j})$$ với $j = 0, 1, 2$. Việc tìm ra biểu thức giải tích hiển cho các tích phân này giải quyết triệt để rào cản tính toán số lặp của các mô hình trước đây.

Mô hình lý thuyết thiết lập quan hệ phi tuyến giữa hàm ứng suất Airy $\Phi(x, y)$ và hàm độ võng $w(x, y)$, bao hàm đại lượng bất hoàn hảo hình học ban đầu $w^*(x, y) = \xi h \sin(m\pi x/a)\sin(n\pi y/b)$. Mô hình biểu diễn thành công sự dịch chuyển cơ chế mất ổn định: đối với vỏ hoàn hảo ($\xi = 0$), kết cấu mất ổn định theo kiểu rẽ nhánh (bifurcation buckling); đối với vỏ không hoàn hảo ($\xi > 0$), kết cấu mất ổn định theo kiểu cực trị (limit point buckling) với sự xuất hiện của tải tới hạn trên $q_{\text{upper}}$ và tải tới hạn dưới $q_{\text{lower}}$ được xác định từ điều kiện triệt tiêu đạo hàm: $$\frac{dp}{dw} = 0 \quad \text{hoặc} \quad \frac{dq}{dw} = 0$$

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

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc lý thuyết:

  1. Lý thuyết biến dạng lớn phi tuyến von Kármán: Mô tả chính xác các thành phần biến dạng mặt giữa $\epsilon_{x0}, \epsilon_{y0}, \gamma_{xy0}$ có chứa tích các số hạng gradient độ võng bậc hai $w_{,x}^2, w_{,y}^2, w_{,x}w_{,y}$ và hiệu ứng tương tác với độ không hoàn hảo $w^*$.
  2. Lý thuyết gân gia cường lệch tâm FGM (ES-FGM): Thiết lập phương trình cơ bản cho hệ gân dọc (stringers) và gân vòng (rings) có tỉ phần thể tích gốm $V_c(z) = (2z/h_s)^{k_2}$ biến thiên đồng bộ, giải quyết sự bất đối xứng hình học và cơ tính thông qua khoảng cách lệch tâm từ tâm gân đến mặt giữa vỏ ($z_s, z_r$).
  3. Mô hình tương tác nền hai tham số Pasternak: Nền đàn hồi được đặc trưng bởi hệ số Winkler $K_1$ (chống lún thẳng đứng) và hệ số Pasternak $K_2$ (chống trượt cắt), tạo nên phản lực nền $q_e = K_1 w - K_2 \nabla^2 w$.

Điều kiện biên (Boundary Conditions) được xây dựng chặt chẽ cho hai trường hợp kỹ thuật:

  • Bốn cạnh tựa đơn (Simply Supported - SSSS): $$w = 0, \quad M_x = 0, \quad N_{xy} = 0, \quad N_x = -r_0 h \quad \text{tại } x = 0, a$$ $$w = 0, \quad M_y = 0, \quad N_{xy} = 0, \quad N_y = -p_0 h \quad \text{tại } y = 0, b$$
  • Hai cạnh cong tựa đơn, hai cạnh thẳng ngàm trượt (Clamped-Sliding - CSCS): $$w = 0, \quad w_{,y} = 0, \quad N_y = 0, \quad N_{xy} = 0 \quad \text{tại } y = 0, b$$

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

Thiết kế nghiên cứu

Nghiên cứu vận hành trên nền tảng bản thể luận hiện thực và nhận thức luận thực chứng diễn dịch (Positivist-Deductive Paradigm) của cơ học môi trường liên tục. Bài toán biên phi tuyến của vật lý toán được mô hình hóa từ các định luật bảo toàn động lượng và tương thích biến dạng:

[Mô hình vật liệu FGM: E(z), ν(z)] + [Hình học vỏ Donnell - von Kármán]
                           │
                           ▼
          [Thiết lập quan hệ Ứng suất - Biến dạng (Hooke)]
                           │
                           ▼
[Hệ phương trình vi phân đạo hàm riêng phi tuyến: w(x,y) & Φ(x,y)]
                           │
                           ▼
 [Áp dụng Kỹ thuật san gân Lekhnitskii & Nền Winkler - Pasternak]
                           │
                           ▼
  [Phép chiếu trực giao Galerkin với hàm độ võng ba số hạng]
                           │
                           ▼
[Hệ phương trình đại số phi tuyến hiển: r0(W), p(W), q(W), τ(W)]
                           │
                           ▼
[Xác định tải tới hạn (Bifurcation / Snap-through) & Vẽ đường cong sau tới hạn]

Quy trình nghiên cứu rigorous

Quy trình giải tích được chuẩn hóa qua 4 giai đoạn nghiêm ngặt:

  1. Thiết lập phương trình tương thích và cân bằng: Đưa vào hàm ứng suất Airy $\Phi(x, y)$ thỏa mãn tự động hai phương trình cân bằng màng $N_{x,x} + N_{xy,y} = 0$ và $N_{xy,x} + N_{y,y} = 0$, biến đổi hệ phương trình cân bằng vỏ mỏng về dạng hai phương trình chủ đạo ghép phi tuyến: $$\nabla_A^4 \Phi + \nabla_B^4 w - \frac{1}{R} w_{,xx} = w_{,xy}^2 - w_{,xx} w_{,yy} + 2 w_{,xy} w_{,xy}^* - w_{,xx} w_{,yy}^* - w_{,yy} w_{,xx}^$$ $$L_1(\Phi) + L_2(w) + \frac{1}{R}\Phi_{,xx} + \Phi_{,yy}(w_{,xx} + w_{,xx}^) - 2\Phi_{,xy}(w_{,xy} + w_{,xy}^) + \Phi_{,xx}(w_{,yy} + w_{,yy}^) + q_0 - q_e = 0$$
  2. Lựa chọn nghiệm dạng giải tích: Hàm độ võng $w(x, y)$ và độ bất hoàn hảo $w^*(x, y)$ được biểu diễn qua chuỗi lượng giác nhiều số hạng thỏa mãn chính xác các điều kiện biên động học và tĩnh học. Đối với bài toán xoắn và áp lực ngoài của vỏ trụ ES-FGM, hàm độ võng ba số hạng được cấu tạo đặc biệt: $$w(x, y) = f_0 + f_1 \sin\frac{m\pi x}{L} \sin\frac{n y}{R} + f_2 \sin^2\frac{m\pi x}{L} \sin^2\frac{n y}{R}$$
  3. Phép chiếu trực giao Galerkin: Thực hiện tích phân giải tích Galerkin trên toàn bộ miền quy chiếu của vỏ $\int_0^a \int_0^b \Omega(x, y) \sin(m\pi x/a)\sin(n\pi y/b) dx dy = 0$ nhằm chuyển đổi hệ phương trình vi phân đạo hàm riêng phi tuyến sang hệ phương trình đại số phi tuyến theo các tham số biên độ độ võng không thứ nguyên $W = f/h$.
  4. Triệt tiêu nghiệm và giải tích cực trị: Khử hoàn toàn các tham số hàm ứng suất để thu được một phương trình đại số duy nhất liên hệ trực tiếp giữa tải trọng tác dụng ($r_0, p, q, \tau$) và độ võng cực đại $W_{max}/h$.

Data và phân tích

Luận án xây dựng thuật toán và chương trình tính toán số hóa trên môi trường phần mềm giải tích biểu thức (Symbolic Computational Software - Maple và Matlab). Dữ liệu đầu vào sử dụng các thông số thực nghiệm chuẩn của vật liệu $Si_3N_4/SUS304$ và $ZrO_2/Ti\text{-}6Al\text{-}4V$:

Đại lượng cơ lý Kim loại ($SUS304$) Gốm ($Si_3N_4$) Hợp kim ($Ti\text{-}6Al\text{-}4V$) Gốm Zirconia ($ZrO_2$)
Môđun Young $E$ (GPa) 207.4 322.27 105.7 168.06
Hệ số Poisson $\nu$ 0.3177 0.2400 0.2981 0.2300
Khối lượng riêng $\rho$ (kg/m$^3$) 8166 2370 4429 5600
Hệ số giãn nở nhiệt $\alpha$ ($10^{-6}/\text{K}$) 15.32 5.87 9.03 10.45

Kiểm tra tính vững (robustness checks) được thực hiện thông qua việc khảo sát sự biến thiên của bước lưới mode vồng $(m, n)$ từ $(1, 1)$ đến $(2, 20)$, quét toàn bộ không gian tham số để tìm giá trị tải tới hạn cực tiểu tuyệt đối $p_{\text{cr}} = \min_{(m,n)} p(m, n)$.

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

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

  1. Hiệu ứng phi tuyến của hệ số Poisson $\nu(z)$: Kết quả số tại Bảng 2.7 chỉ ra rằng khi tính đến $\nu = \nu(z)$, tải tới hạn $p_{cr}$ của vỏ trụ FGM ($L/R = 1, R/h = 100, k = 1$) đạt 138.45 MPa, trong khi mô hình $\nu = \text{const} = 0.3$ cho giá trị 132.18 MPa (chênh lệch 4.75%). Khi chỉ số $k$ tăng lên $k = 5$, sai số đạt mức 6.2%, chứng minh giả định hệ số Poisson hằng số gây đánh giá thấp khả năng chịu tải thực tế.
  2. Hiệu ứng gia cường vượt bậc của gân ES-FGM: Đối với vỏ trụ chịu áp lực ngoài (Bảng 3.3 và Hình 3.9), việc bố trí hệ gân FGM ($n_s = 20, n_r = 20$) làm tăng tải tới hạn $q_{cr}$ từ 2.34 lần đến 4.12 lần so với vỏ FGM không gân cùng chiều dày tương đương. Gân đặt ở mặt trong (internal stiffeners) cho khả năng chịu nén dọc trục và áp lực ngoài cao hơn từ 12% đến 18% so với gân đặt ở mặt ngoài (external stiffeners).
  3. Hiện tượng chuyển đổi độ cứng sau tới hạn: Đồ thị quan hệ tải - độ võng (Hình 2.2 và 2.3) phát hiện hiện tượng đảo chiều cơ tính: Ở giai đoạn biến dạng bé ($W/h < 0.5$), độ không hoàn hảo $\xi$ làm suy giảm nghiêm trọng độ cứng và tải tới hạn của kết cấu. Tuy nhiên, khi độ võng đủ lớn ($W/h > 1.2$), biến dạng màng chiếm ưu thế làm xuất hiện hiệu ứng hóa cứng hình học (geometric strain hardening), khiến đường cong tải - độ võng của kết cấu có $\xi = 0.5$ vươn lên cao hơn kết cấu hoàn hảo $\xi = 0$.
  4. Quy luật chi phối của nền Pasternak: Nền đàn hồi Pasternak ($K_2$) có hiệu quả ổn định hóa cao gấp 3.5 đến 5.2 lần so với nền Winkler ($K_1$). Dưới tác dụng của nền đàn hồi hai tham số, mode vồng nguy hiểm nhất của vỏ nón FGM bị dịch chuyển từ mode sóng chu vi thấp ($n = 2, 3$) sang mode sóng chu vi cao ($n = 6, 8$), đồng thời triệt tiêu hoàn toàn vùng mất ổn định nhảy bước (snap-through) nguy hiểm.
  5. Ảnh hưởng phi tuyến của góc côn $\alpha$ ở vỏ nón FGM: Lực nén tới hạn $P_{cr}$ và áp lực ngoài $q_{cr}$ của vỏ nón FGM có gân giảm phi tuyến theo quy luật hàm $\cos^2\alpha$ khi góc bán đỉnh $\alpha$ tăng từ $0^\circ$ (vỏ trụ) lên $60^\circ$ (Hình 4.16 và 4.17).
So sánh tải tới hạn nén trục Pcr (MN) của vỏ nón FGM theo cấu hình gân và nền:
[Vỏ FGM trơn, không nền]          : █ 1.45 MN
[Vỏ FGM + Nền Pasternak]          : ████ 3.82 MN
[Vỏ ES-FGM + Gân trong]           : ████████ 7.15 MN
[Vỏ ES-FGM + Gân trong + Nền Pas] : ████████████████████ 14.62 MN

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng thành công hệ tiên đề Donnell - von Kármán cho môi trường FGM có hệ số Poisson biến thiên và kết cấu vỏ nón cụt có gân gia cường, hoàn thiện lý thuyết tấm và vỏ dị hướng chức năng.
  • Về mặt phương pháp luận: Cung cấp quy trình giải tích chuẩn hóa kết hợp kỹ thuật san gân Lekhnitskii với phương pháp Galerkin ba số hạng, có thể chuyển giao trực tiếp để phân tích các dạng kết cấu phức tạp khác như vỏ cầu thoải FGM, vỏ hypebolit hay tấm lượn sóng FGM.
  • Về mặt ứng dụng thực tiễn: Cho phép các kỹ sư thiết kế tối ưu hóa hình học gân (chiều cao $h_s, h_r$, chiều rộng $b_s, b_r$), lựa chọn gradient thành phần gốm $k$ tối ưu nhằm giảm trọng lượng kết cấu hàng không vũ trụ tới 25-35% mà vẫn đảm bảo an toàn ổn định tĩnh tuyệt đối.
  • Về mặt chính sách và quy chuẩn: Đóng góp cơ sở toán học thực nghiệm để xây dựng các tiêu chuẩn tính toán kết cấu chịu nhiệt độ cao trong ngành công nghiệp quốc phòng và năng lượng nguyên tử quốc gia.

Limitations và Future Research

Luận án duy trì tính khách quan khoa học khi xác định 4 giới hạn nội tại:

  1. Lý thuyết vỏ mỏng cổ điển: Việc áp dụng lý thuyết Donnell - Kirchhoff bỏ qua biến dạng trượt cắt ngang ($\gamma_{xz} = \gamma_{yz} = 0$). Giả định này chỉ bảo đảm độ chính xác cao cho vỏ mỏng ($R/h \ge 50$), nhưng sẽ gây sai số vượt quá 10% khi áp dụng cho vỏ dày ($R/h < 20$) do ảnh hưởng của hiệu ứng trượt cắt bậc cao.
  2. Mô hình hóa gân bằng phương pháp san đều: Kỹ thuật Lekhnitskii phân bố độ cứng gân liên tục trên bề mặt, do đó không thể nắm bắt được hiện tượng mất ổn định cục bộ của bản gân (rib local buckling) hoặc hiện tượng tập trung ứng suất tại chân gân khi khoảng cách giữa các gân quá lớn.
  3. Phạm vi tải trọng tĩnh và đàn hồi tuyến tính: Nghiên cứu chỉ giới hạn ở bài toán tĩnh và vật liệu làm việc trong miền đàn hồi tuyến tính của Hooke, chưa xét đến biến dạng dẻo, hiện tượng rão (creep) ở nhiệt độ cao và đáp ứng động lực học chịu tải xung kích va đập.
  4. Điều kiện biên lý tưởng hóa: Các điều kiện biên tựa đơn (SSSS) và ngàm trượt (CSCS) được giải bằng chuỗi lượng giác giải tích; các điều kiện biên phức tạp thực tế như đàn hồi phi tuyến cục bộ chưa được tích hợp hoàn toàn.

Chương trình nghiên cứu 5-10 năm tới bao gồm:

  • Mở rộng phân tích ổn định phi tuyến sử dụng Lý thuyết biến dạng trượt bậc cao (HSDT) và Lý thuyết biến dạng trượt bậc nhất (FSDT) cho vỏ dày FGM.
  • Phát triển bài toán ổn định động nhiệt điện từ (magneto-electro-thermo-mechanical dynamic stability) có xét đến hiệu ứng áp điện của vật liệu FGM thông minh.
  • Phân tích hiện tượng mất ổn định đàn dẻo (elasto-plastic buckling) và phá hủy mỏi của kết cấu FGM có lỗ khoét và khuyết tật tế vi (porosities).
  • Tối ưu hóa topo phân bố gân FGM tự do bằng phương pháp giải tích kết hợp trí tuệ nhân tạo (AI/Machine Learning).

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

  • Ảnh hưởng học thuật: Luận án và các công bố liên quan trên các tạp chí chuyên ngành uy tín (như Vietnam Journal of Mechanics, các tạp chí cơ học quốc tế ISI/Scopus) tạo tiền đề cho hàng chục trích dẫn khoa học trong các nghiên cứu tiếp theo về cơ học kết cấu FGM tại Việt Nam và quốc tế.
  • Chuyển giao công nghiệp: Mô hình giải tích giúp rút ngắn thời gian tính toán mô phỏng từ hàng giờ trên các phần mềm phần tử hữu hạn (ANSYS, ABAQUS) xuống còn vài giây thông qua các công thức hiển, phục vụ đắc lực cho khâu thiết kế sơ bộ (conceptual design) trong công nghiệp hàng không, tên lửa và chế tạo bồn bể áp lực composite.
  • Ý nghĩa xã hội và an ninh năng lượng: Nâng cao độ tin cậy và hệ số an toàn cho các kết cấu vỏ lò phản ứng hạt nhân và bồn chứa khí hóa lỏng dung tích lớn, ngăn ngừa nguy cơ sự cố thảm họa do mất ổn định oằn vỏ gây ra.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp giải tích chặt chẽ, nắm bắt kỹ thuật thiết lập hàm ứng suất Airy và giải hệ phương trình vi phân phi tuyến biến hệ số bằng phương pháp Galerkin.
  • Các nhà khoa học cơ học tính toán: Sử dụng các công thức giải tích hiển $A_{ij}, r_0(W), q_{cr}, \tau_{cr}$ làm nghiệm chuẩn (benchmark solutions) tin cậy để kiểm chứng tính chính xác của các phần mềm và thuật toán phần tử hữu hạn mới.
  • Kỹ sư R&D hàng không vũ trụ và quốc phòng: Ứng dụng trực tiếp công thức tính toán độ dày tối ưu $h$, tỷ số gân $h_s/b_s$, và chỉ số phân bố vật liệu $k$ để thiết kế vỏ nhẹ chịu áp lực và tải xoắn cực hạn.
  • Cơ quan quản lý an toàn công trình: Có thêm cơ sở khoa học định lượng để ban hành quy chuẩn kiểm định kết cấu bồn bể, đường ống chịu áp lực làm từ vật liệu mới.

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 thiết lập thành công hệ nghiệm giải tích hiển cho hệ số độ cứng $A_{ij}$ trong điều kiện hệ số Poisson biến thiên liên tục theo chiều dày $\nu = \nu(z)$, kết hợp mở rộng lý thuyết vỏ Donnell - von Kármán cho kết cấu vỏ trụ và vỏ nón có gân gia cường FGM lệch tâm (ES-FGM). Nghiên cứu đã giải quyết được hạn chế tồn tại nhiều năm trong các công trình của Huang & Han (2008) và Shen (2002) vốn phải dựa vào tích phân số hoặc giả định $\nu = \text{const}$.

2. Đột phá về phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đó? So với nghiên cứu của Najafizadeh và cộng sự (2009) chỉ phân tích tuyến tính vỏ trụ FGM có gân, luận án đã xây dựng lời giải phi tuyến tĩnh với hàm độ võng ba số hạng, cho phép mô tả chính xác điểm rẽ nhánh và điểm cực trị (snap-through). So với Sofiyev (2007, 2011) chỉ khảo sát vỏ nón FGM không gân bằng phương pháp một số hạng đơn giản, luận án đã tích hợp đồng thời kỹ thuật san đều gân Lekhnitskii cải tiến cho hệ gân FGM và tương tác nền đàn hồi hai tham số Pasternak.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm tính toán là gì? Phát hiện bất ngờ nhất là hiện tượng đảo chiều đáp ứng độ võng sau tới hạn (Hình 2.2): Ở độ võng bé, độ không hoàn hảo ban đầu làm giảm tải trọng chịu đựng; nhưng khi độ võng vượt ngưỡng $W/h > 1.2$, kết cấu có độ không hoàn hảo lớn lại có khả năng mang tải sau vồng cao hơn kết cấu hoàn hảo nhờ sự kích hoạt sớm của trường ứng suất màng kéo căng (membrane tensile stress field).

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Có. Toàn bộ các hệ thức cấu tạo vật liệu, các tích phân tường minh $A_{ij}$, hệ số đại số $C_1 \div C_4$, biểu thức giải tích của tải vồng cận trên $r_{0\text{upper}}$ (công thức 2.26) và tải vồng cận dưới $r_{0\text{lower}}$ (công thức 2.27), cùng hệ thống bảng số liệu cơ lý tính chi tiết tại Chương 2, 3, 4 và Phụ lục A đều được công khai hoàn chỉnh, cho phép tái lập 100% các kết quả số và đường cong đồ thị.

5. Chương trình nghiên cứu 10 năm được định hình như thế nào? Chương trình nghiên cứu hướng tới: (1) Nâng cấp mô hình từ vỏ mỏng Donnell lên lý thuyết biến dạng trượt bậc cao Reddy (HSDT); (2) Khảo sát mất ổn định nhiệt đàn dẻo và phá hủy liên kết trong điều kiện cháy nổ; (3) Tích hợp cấu trúc rỗng xốp micro/nano (porous FGM); (4) Tự động hóa thiết kế tối ưu hóa cấu trúc FGM thông qua thuật toán di truyền và mạng nơ-ron nhân tạo.

Kết luận

  1. Hoàn thiện mô hình giải tích hệ số Poisson biến thiên: Thiết lập thành công các biểu thức hiển cho các hệ số độ cứng $A_{ij}$ khi $\nu = \nu(z)$, chứng minh sai số lên tới 6.2% nếu giả định $\nu = \text{const}$ trong bài toán ổn định vỏ mỏng FGM.
  2. Tiên phong giải bài toán vỏ trụ FGM có gân FGM lệch tâm (ES-FGM): Xây dựng hệ phương trình chủ đạo phi tuyến và tìm ra biểu thức hiển xác định tải tới hạn nén dọc trục, áp lực ngoài và tải xoắn thuần túy với hàm độ võng ba số hạng.
  3. Giải quyết bài toán vỏ nón cụt FGM có gân trên nền đàn hồi: Thiết lập phương pháp giải tích tường minh giải hệ phương trình đạo hàm riêng biến hệ số trên hệ tọa độ cong, chứng minh hiệu quả triệt tiêu mất ổn định sóng dài của nền Pasternak.
  4. Khám phá quy luật cơ học sau tới hạn: Làm sáng tỏ cơ chế chuyển đổi từ mất ổn định rẽ nhánh sang mất ổn định cực trị (snap-through) dưới tác động đồng thời của độ không hoàn hảo hình học ban đầu $\xi$ và chỉ số tỉ phần thể tích $k$.
  5. Cung cấp hệ thống nghiệm chuẩn quốc tế: Đóng góp bộ dữ liệu nghiệm giải tích hiển có độ chính xác cao, được kiểm chứng trùng khớp với các công trình kinh điển của Turvey, Shen, Brush & Almroth, Seide, và Sofiyev, tạo nền tảng vững chắc cho sự phát triển của cơ học tính toán vật liệu composite tiên tiến tại Việt Nam.