Tổng quan về luận án

Luận án tiến sĩ cơ học với tiêu đề "Phân tích ổn định và động lực phi tuyến của vỏ thoải sandwich cơ tính biến thiên chịu tác dụng của tải cơ và nhiệt" (Chuyên ngành: Cơ học vật rắn, Mã số: 62440107) do nghiên cứu sinh Đặng Thùy Đông thực hiện dưới sự hướng dẫn khoa học của PGS. Vũ Đỗ Long và cố GS. Đào Văn Dũng tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội. Công trình đại diện cho bước tiến phong trong lĩnh vực cơ học tính toán vật liệu và kết cấu composite tiên tiến tại Việt Nam, tập trung giải quyết bài toán ứng xử cơ học phi tuyến phức tạp của các kết cấu vỏ thoải sandwich có cơ tính biến thiên (Functionally Graded Material - FGM) gia cường gân FGM đa hướng và vỏ gấp nếp đặt trên nền đàn hồi chịu tác động đồng thời của trường nhiệt độ và tải trọng cơ học biến thiên theo thời gian.

Trong cơ học kết cấu hiện đại, các cấu trúc sandwich truyền thống thường đối mặt với nguy cơ phá hủy nguy hiểm do hiện tượng tập trung ứng suất cục bộ tại bề mặt tiếp giáp giữa các lớp vật liệu có cơ tính bất liên tục. Luận án trích dẫn trực tiếp thực trạng khoa học: "nhược điểm lớn nhất của loại kết cấu này là mỗi lớp có tính chất vật liệu khác nhau, dẫn đến tập trung ứng suất tại các vị trí tiếp giáp giữa các lớp, gây hiện tượng bong tách, phá hủy kết cấu. Để khắc phục vấn đề này, kết cấu vật liệu sandwich cơ tính biến thiên đã được đề xuất." Mặc dù các nghiên cứu quốc tế đã tiếp cận kết cấu FGM tấm và vỏ, phần lớn các công bố trước đó chỉ giới hạn ở vỏ trơn thuần nhất, vỏ FGM không gân hoặc chỉ gia cường bằng gân thuần nhất trực giao theo lý thuyết vỏ cổ điển (CST), chưa giải quyết toàn diện bài toán vỏ sandwich FGM đa lớp kết hợp hệ thống gân xiên FGM chịu nhiệt, chỏm cầu thoải biến dạng đối xứng trục theo biến dạng trượt giải bằng phương pháp chuyển vị, và kết cấu vỏ trống/vỏ trụ FGM gấp nếp (corrugated) có lõi gấp nếp đặt trên nền đàn hồi Winkler-Pasternak.

Mục tiêu khoa học của luận án được cụ thể hóa thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố cơ tính biến thiên mở rộng (Sigmoid mở rộng và Lũy thừa mở rộng) cùng sự phụ thuộc nhiệt độ $Pr(T, z)$ ảnh hưởng như thế nào đến ứng xử ổn định tĩnh phi tuyến và tải tới hạn vồng của vỏ thoải hai độ cong sandwich FGM?
  2. Câu hỏi nghiên cứu 2 (RQ2): Sự kết hợp giữa lý thuyết biến dạng trượt bậc cao (HSDT), bậc nhất (FSDT) với mô hình gân FGM xiên ($\theta$) và gân trực giao san đều (Lekhnitskiy) phản ánh chính xác đáp ứng dao động phi tuyến và tải tới hạn động (Budiansky-Roth) ra sao dưới tác dụng của tải cơ - nhiệt kết hợp?
  3. Câu hỏi nghiên cứu 3 (RQ3): Quy luật ứng xử động lực phi tuyến, hiện tượng phách điều hòa và quỹ đạo pha của chỏm cầu thoải sandwich FGM biến dạng đối xứng trục thay đổi như thế nào khi giải trực tiếp bằng phương pháp chuyển vị?
  4. Câu hỏi nghiên cứu 4 (RQ4): Cấu trúc vỏ trống và vỏ trụ sandwich FGM gấp nếp hình sin/hình thang kết hợp nền đàn hồi hai tham số Pasternak gia tăng độ cứng chống mất ổn định tĩnh và động đến mức độ định lượng nào so với kết cấu vỏ trơn tương đương?

Tương ứng với các câu hỏi trên, hệ thống 4 giả thuyết khoa học ($H_1 - H_4$) được xây dựng, định vị rằng việc tích hợp gân gia cường FGM xiên và cấu trúc lõi gấp nếp sẽ loại bỏ triệt để điểm kỳ dị tập trung ứng suất, nâng cao vượt bậc tải tới hạn động $P_{crd}$ và dịch chuyển tần số dao động tự do cơ bản $\omega_0$ ra khỏi miền cộng hưởng nguy hiểm dưới các vận tốc gia tải xung biến thiên ($c_p, c_q$).

Phạm vi nghiên cứu bao quát các dạng hình học phức tạp thỏa mãn điều kiện vỏ thoải của Vlasov: "Mặt giữa của vỏ thoải có độ cong Gauss rất nhỏ nên hình học của mặt này vốn là hình học Reimamn có thể thay thế bằng hình học Euclide trên mặt phẳng" với điều kiện biên hình học nghiêm ngặt $f_n \le \min(a, b)/5$ và $h \le \min(R_x, R_y)/5$. Các mô hình khảo sát bao gồm: vỏ thoải hai độ cong ($R_x, R_y$), panel cầu ($R_x=R_y$), panel trụ ($R_x=\infty$), vỏ hình yên ngựa hyperbolic ($R_x R_y < 0$), chỏm cầu thoải biến dạng đối xứng trục ($a, R, h$) và vỏ trống tròn xoay gấp nếp có gân gia cường đặt trên nền đàn hồi hai tham số ($K_1, K_2$).

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận sự phát triển mạnh mẽ của cơ học kết cấu FGM từ các công trình nền tảng của Shen (2009-2014) về đáp ứng vồng và sau vồng của panel trụ/cầu FGM đơn lớp sử dụng lý thuyết HSDT kết hợp kỹ thuật nhiễu hai bước (two-step perturbation technique). Song song đó, Tornabene và Viola (2007) phát triển phương pháp cầu phương vi phân tổng quát (GDQ) dựa trên FSDT cho vỏ thoải hai độ cong FGM, trong khi Matsunaga (2008) và Alijani et al. (2011) khảo sát dao động cưỡng bức phi tuyến của vỏ FGM bằng lý thuyết Donnell và HSDT. Về kết cấu sandwich FGM, Alibeigoo (2011), Alibeigoo và Liew (2014) thiết lập lời giải giải tích dựa trên lý thuyết đàn hồi 3D phân tích dao động tự do của panel trụ sandwich FGM, còn Pandey và Pradyumna (2015) phát triển mô hình phần tử hữu hạn (FEM) đánh giá ảnh hưởng của các mô hình đồng nhất hóa Voigt và Mori-Tanaka.

Trong trường phái nghiên cứu kết cấu gia cường gân, Đào Huy Bích và cộng sự (2013-2015), Nguyễn Đình Đức và cộng sự (2014-2017) đã tiên phong áp dụng kỹ thuật san đều tác dụng gân của Lekhnitskiy kết hợp tính phi tuyến hình học von Kármán cho tấm và vỏ FGM. Tuy nhiên, các nghiên cứu này chủ yếu dừng lại ở gân thuần nhất trực giao, khảo sát trên nền lý thuyết vỏ mỏng cổ điển CST hoặc gân FGM chưa xét đến biến dạng nhiệt đồng thời trong cả vỏ và gân dầm trượt.

Trường phái / Tác giả Mô hình kết cấu & Vật liệu Khung lý thuyết & Phương pháp Khoảng trống học thuật (Research Gap)
Shen et al. (2009–2014) Panel trụ / cầu FGM đơn lớp HSDT, Kỹ thuật nhiễu 2 bước, Galerkin Chưa xét kết cấu sandwich đa lớp, không có gân gia cường FGM xiên
Alibeigoo & Liew (2014) Panel trụ sandwich FGM Lý thuyết đàn hồi 3D giải tích, chuỗi Fourier Bỏ qua hệ thống gân gia cường và cấu trúc gấp nếp trên nền đàn hồi
Pandey & Pradyumna (2015) Tấm và vỏ thoải sandwich FGM FEM, FSDT/HSDT, mô hình Mori-Tanaka Chưa khảo sát ổn định động theo Budiansky-Roth dưới tải nhiệt phụ thuộc
Đào Huy Bích et al. (2013–2017) Vỏ thoải FGM, chỏm cầu FGM CST, hàm ứng suất Airy, Runge-Kutta 4 Chủ yếu dùng gân thuần nhất trực giao, chỏm cầu chưa giải theo chuyển vị FSDT
Đặng Thùy Đông (Luận án) Vỏ thoải sandwich FGM, chỏm cầu, vỏ trống gấp nếp HSDT & FSDT, Hàm ứng suất + Chuyển vị, Lekhnitskiy, Homogenization Xia Khắc phục toàn diện: Gân FGM xiên chịu nhiệt, giải chuyển vị FSDT, vỏ gấp nếp có nền Pasternak

Tranh luận học thuật sâu sắc tồn tại giữa hai quan điểm: (1) Quan điểm sử dụng lý thuyết biến dạng trượt bậc nhất (FSDT) với hệ số hiệu chỉnh cắt cố định $K_s = 5/6$ (như Zenkour, 2006; Bodaghi & Shakeri, 2012) mang lại tính đơn giản trong giải tích nhưng phụ thuộc vào hệ số hiệu chỉnh biên; và (2) Quan điểm bắt buộc sử dụng lý thuyết biến dạng trượt bậc cao (HSDT của Reddy, 2004) để triệt tiêu hoàn toàn ứng suất trượt ngang trên các mặt tự do của vỏ dày trung bình mà không cần hệ số $K_s$. Luận án định vị chính xác vị thế học thuật bằng cách xây dựng song song hai hệ phương trình chủ đạo theo cả FSDT và HSDT, chứng minh sự tương thích và chỉ rõ biên độ sai lệch của FSDT khi chiều dày tương đối $h/\min(a,b)$ gia tăng.

So sánh trực tiếp với hai nghiên cứu quốc tế chuẩn tắc:

  • So sánh lực nén dọc trục tới hạn không thứ nguyên $\bar{P} = P_{cr} b^2 / (\pi^2 D)$ của tấm/panel trụ với Zenkour (2006): Kết quả phân tích của luận án theo HSDT cho sai số tuyệt đối dưới $0.42%$, khẳng định độ tin cậy tuyệt đối của thuật toán giải tích.
  • So sánh tham số tần số dao động cơ bản không thứ nguyên $\bar{\omega} = \omega h \sqrt{\rho_0/E_0} (a/h)^2$ của chỏm cầu thoải với công bố kinh điển của Varadan & Pandalai (1979) và Sathyamoorthy (1978): Sai số tương đối dao động từ $0.15%$ đến $0.68%$, hoàn toàn vượt trội về độ chính xác khi mô tả tính phi tuyến biên dạng cong.

Đó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 và phát triển đột phá 3 hệ thống lý thuyết nền tảng trong cơ học vật rắn biến dạng:

  1. Lý thuyết vật liệu cơ tính biến thiên đa lớp (S-FGM & P-FGM Generalized): Thiết lập mô hình giải tích 4 biến thể cấu trúc sandwich FGM gồm Mô hình 1A (FGM – Kim loại – FGM), Mô hình 1B (FGM – Gốm – FGM) theo quy luật Sigmoid mở rộng; và Mô hình 2A (Gốm – FGM – Kim loại), Mô hình 2B (Kim loại – FGM – Gốm) theo quy luật Lũy thừa mở rộng. Cụ thể, hàm cơ tính hiệu dụng phân bố theo chiều dày $z \in [-h/2, h/2]$ được chuẩn hóa: $$Pr^{sh}(z, T) = \begin{cases} Pr_i + Pr_{ji}\left(\frac{2z+h}{2h_t}\right)^{k_t}, & -h/2 \le z \le -h/2+h_t \ Pr_j, & -h/2+h_t \le z \le h/2-h_b \ Pr_i + Pr_{ji}\left(\frac{2z-h}{2h_b}\right)^{k_b}, & h/2-h_b \le z \le h/2 \end{cases}$$ với $Pr \in {E, \rho, \alpha}$ và các tham số phụ thuộc nhiệt độ phi tuyến $Pr_i(T) = P_0(1 + P_{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$.
  2. Lý thuyết dầm trượt san đều Lekhnitskiy cải biên cho gân FGM xiên đa hướng: Luận án mở rộng kỹ thuật san đều gân kinh điển của Lekhnitskiy (1968) sang hệ thống gân làm bằng vật liệu FGM định hướng nghiêng một góc $\theta$ bất kỳ ($0^\circ \le \theta \le 90^\circ$). Mô hình tính đến biến dạng trượt ngang của gân và sự phân bố nhiệt độ đồng thời trong kết cấu dầm - vỏ, loại bỏ hoàn toàn giả thiết gân thuần nhất đẳng nhiệt trước đây.
  3. Mở rộng lý thuyết phi tuyến hình học von Kármán cho vỏ thoải có độ không hoàn hảo ban đầu: Đưa vào biểu thức độ võng ban đầu $w^*(x,y)$ vào toàn bộ phương trình tương thích biến dạng và phương trình cân bằng/chuyển động động học theo HSDT, cho phép mô phỏng chính xác hiện tượng biến dạng phân nhánh (bifurcation) và hiện tượng vồng sụp (snap-through).

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

Khung phân tích của luận án là sự tích hợp ma trận hữu cơ giữa lý thuyết biến dạng trượt bậc cao HSDT của Reddy, lý thuyết vỏ thoải Donnell-Vlasov, phương pháp hàm ứng suất Airy và phương pháp chuyển vị trực tiếp. Điểm độc đáo then chốt thể hiện ở quy trình chuyển đổi hệ phương trình vi phân đạo hàm riêng phi tuyến bậc cao sang hệ phương trình vi phân thường phi tuyến theo thời gian $t$ bằng kỹ thuật biến phân Galerkin:

[Mô hình hình học & Vật liệu Sandwich FGM + Gân FGM + Nền Pasternak]
[Hệ thức Động học HSDT / FSDT + Phi tuyến von Kármán + Sai lệch w*]
[Định luật Hooke mở rộng chịu nhiệt E(z,T), α(z,T) + San đều Lekhnitskiy]
[Thiết lập Phương trình Tương thích Biến dạng & Cân bằng Động lượng]
[Phương pháp Hàm ứng suất Airy (f) / Phương pháp Chuyển vị (u, v, w)]
[Trực giao hóa Galerkin trên miền biên tựa đơn/ngàm (FM, IM)]
[Hệ PT Vi phân Thường Phi tuyến Động lực học (Duffing-type Equations)]
[Phân tích Ổn định Tĩnh]        [Tích phân Runge-Kutta Bậc 4]
                                [Tiêu chuẩn Budiansky-Roth]
                                (Xác định Tải tới hạn động P_crd)

Điều kiện biên nghiên cứu bao gồm 3 cấu hình thực tế: Bốn cạnh tựa đơn tự do di trượt (FM), Bốn cạnh tựa đơn cố định không di trượt (IM), và Biên hỗn hợp (hai cạnh tự do, hai cạnh cố định), thiết lập các giới hạn biên toán học tường minh cho bài toán.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng duy lý (positivism) kết hợp mô hình hóa giải tích cơ học giải tích nghiêm ngặt (rigorous analytical modeling) và mô phỏng số kiểm chứng (numerical simulation). Thiết kế đa cấp độ (multi-level mechanical design) bao quát từ cấp độ vi mô (tỷ phần thể tích pha gốm/kim loại $V_c, V_m$, gradient phân bố hàm mũ/sigmoid), cấp độ trung mô (tương tác tiếp xúc liên tục hoàn hảo giữa các lớp vỏ và dầm gân gia cường, đồng nhất hóa độ cứng tương đương nếp gấp của Xia et al.), đến cấp độ vĩ mô (ứng xử tĩnh - động của toàn bộ kết cấu vỏ đặt trên nền đàn hồi hai hệ số $K_1 = 10^6 \text{ N/m}^3, K_2 = 10^5 \text{ N/m}$).

Quy trình nghiên cứu rigorous

Quy trình giải quyết bài toán cơ học phi tuyến trong luận án tuân thủ chuỗi thao tác đại số và giải tích toán học chuẩn xác:

  1. Thiết lập trường chuyển vị và biến dạng: Theo HSDT, trường chuyển vị tại điểm $(x,y,z)$ được định nghĩa chuẩn xác: $$u(x,y,z,t) = u_0(x,y,t) + z\phi_x(x,y,t) - c_1 z^3\left(\phi_x + \frac{\partial w}{\partial x}\right)$$ $$v(x,y,z,t) = v_0(x,y,t) + z\phi_y(x,y,t) - c_1 z^3\left(\phi_y + \frac{\partial w}{\partial y}\right)$$ $$w(x,y,z,t) = w_0(x,y,t)$$ với hệ số hình học $c_1 = 4/(3h^2)$, đảm bảo điều kiện triệt tiêu ứng suất cắt tiếp tuyến tại mặt trên và mặt dưới ($z = \pm h/2$).
  2. Khử biến dạng màng và thiết lập hệ phương trình đạo hàm riêng: Áp dụng phương pháp hàm ứng suất $f(x,y,t)$ với $N_x = \partial^2 f/\partial y^2, N_y = \partial^2 f/\partial x^2, N_{xy} = -\partial^2 f/\partial x \partial y$, thỏa mãn tự động hai phương trình cân bằng nội lực trong mặt phẳng.
  3. Trực giao hóa biến phân Galerkin: Sử dụng hàm dạng dao động thỏa mãn trực tiếp các điều kiện biên: $$w(x,y,t) = W(t) \sin(\lambda_m x) \sin(\lambda_n y), \quad w^(x,y) = W^ \sin(\lambda_m x) \sin(\lambda_n y)$$ với $\lambda_m = m\pi/a, \lambda_n = n\pi/b$. Chiếu hệ phương trình liên tục lên không gian hàm thử để thu được phương trình vi phân chuyển động phi tuyến bậc hai dạng Duffing mở rộng đối với biên độ độ võng không thứ nguyên $\xi(t) = W(t)/h$: $$m_{eff} \frac{d^2\xi}{dt^2} + 2\epsilon m_{eff} \frac{d\xi}{dt} + k_1 \xi + k_2 \xi^2 + k_3 \xi^3 = F(t, T)$$
  4. Tiêu chuẩn ổn định động Budiansky-Roth: Tải tới hạn động ($P_{crd}$ hoặc $q_{crd}$) được xác định tại điểm mà tốc độ biến thiên độ võng cực đại theo tải trọng đạt bước nhảy đột biến (inflection point of the dynamic response curve): $$\frac{d(\xi_{max})}{dP} \to \infty \quad \text{hoặc} \quad \frac{d^2(\xi_{max})}{dP^2} = 0$$

Data và phân tích

Dữ liệu tính toán số sử dụng hệ vật liệu gốm nhiệt độ cao Silicon Nitride ($Si_3N_4$) và hợp kim Thép không gỉ (SUS304) hoặc Hợp kim Titan (Ti-6Al-4V). Các hệ số phụ thuộc nhiệt độ chuẩn hóa từ thực nghiệm nhiệt cơ:

  • $Si_3N_4$: $E = 348.43 \times 10^9 \text{ Pa}, \nu = 0.24, \rho = 2370 \text{ kg/m}^3, \alpha = 5.87 \times 10^{-6} \text{ K}^{-1}$
  • SUS304: $E = 201.04 \times 10^9 \text{ Pa}, \nu = 0.3, \rho = 8166 \text{ kg/m}^3, \alpha = 12.33 \times 10^{-6} \text{ K}^{-1}$

Các thuật toán tích phân số được lập trình toàn diện trên môi trường tính toán biểu thức giải tích và số học tự động (MATLAB / Maple), thực hiện tích phân Runge-Kutta bậc 4 với bước thời gian hội tụ $\Delta t = 10^{-5} \text{ s}$, đảm bảo tính ổn định số tuyệt đối (numerical stability) và triệt tiêu sai số làm tròn trong toàn bộ miền thời gian khảo sát $t \in [0, 0.1 \text{ s}]$.

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

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

  1. Hiệu ứng gia cường của hệ thống gân FGM xiên ($\theta = 45^\circ$): Luận án chứng minh định lượng rằng việc bố trí gân xiên FGM $45^\circ$ kết hợp gân trực giao (hệ gân lưới) làm tăng tải tới hạn nén dọc trục $P_{cr}$ của panel trụ sandwich FGM lên đến $42.8% - 67.3%$ so với kết cấu vỏ không gân, đồng thời vượt trội hơn $18.5%$ so với việc chỉ bố trí gân trực giao truyền thống có cùng thể tích vật liệu gia cường.
  2. Quy luật phân bố Sigmoid (Mô hình 1B) tối ưu hóa khả năng chịu nhiệt: Kết cấu vỏ sandwich FGM Mô hình 1B (lõi gốm thuần nhất, hai lớp phủ FGM) thể hiện độ ổn định nhiệt cơ vượt trội nhất. Khi nhiệt độ tăng từ $T_0 = 300\text{ K}$ lên $T_0 + \Delta T = 500\text{ K}$, tải tới hạn tĩnh của Mô hình 1B chỉ suy giảm $14.2%$, trong khi Mô hình 2A (lớp phủ gốm/kim loại thuần nhất, lõi FGM) bị suy giảm đến $38.6%$ do hiện tượng mềm hóa nhiệt cục bộ của lớp kim loại mặt dưới.
  3. Hiện tượng phách điều hòa (Harmonic Beating) và chuyển pha phi tuyến trong chỏm cầu thoải: Lần đầu tiên phát hiện hiện tượng phách điều hòa rõ nét trong đáp ứng dao động phi tuyến của chỏm cầu thoải sandwich FGM dưới tác dụng của tải điều hòa $q(t) = Q_0 \sin(\Omega t)$ khi tần số kích thích $\Omega \approx 2500 \text{ rad/s}$ tiệm cận tần số dao động tự do $\omega_0$. Đường cong quỹ đạo pha (phase plane trajectory $\dot{W} - W$) chuyển đổi từ quỹ đạo ellipse tuần hoàn khép kín sang các vòng xoắn hỗn loạn mở rộng khi biên độ tải vượt ngưỡng áp lực tới hạn tĩnh $Q_0 > 10^6 \text{ N/m}^2$.
  4. Hiệu năng vượt trội của kết cấu vỏ trống gấp nếp hình thang và lõi gấp nếp: Kết cấu vỏ trống sandwich FGM lõi gấp nếp làm tăng khả năng chịu tải nén dọc trục và áp lực ngoài tới hạn lên $2.3 - 3.1$ lần so với vỏ trống sandwich lõi đặc có cùng khối lượng. Nền đàn hồi Pasternak ($K_1 = 10^6 \text{ N/m}^3, K_2 = 10^5 \text{ N/m}$) đóng vai trò nâng đỡ đặc biệt quan trọng, làm tăng tải tới hạn động $P_{crd}$ thêm $31.5%$ và triệt tiêu nhanh dao động quá độ nhờ lớp cản nhớt.
  5. Ảnh hưởng phi tuyến của tốc độ gia tải cơ học ($c_p, c_q$): Dưới tốc độ gia tải nhanh ($c_p = 3 \text{ GPa/s}$, $c_q = 10^6 \text{ N/m}^2\text{s}$), tải tới hạn động $P_{crd}$ luôn lớn hơn tải tới hạn tĩnh $P_{cr}^{upper}$ từ $15.8%$ đến $48.2%$ do hiệu ứng quán tính cản trở sự phát triển tức thời của các mode mất ổn định dạng sóng.

Implications đa chiều

  • Về mặt lý thuyết: Thiết lập hệ thống công thức giải tích dạng hiển (closed-form solutions) cho tần số dao động tự do, quan hệ tải trọng - độ võng tĩnh, và phương trình vi phân chuyển động cho 4 lớp kết cấu vỏ thoải phức tạp, làm phong phú kho tàng lý thuyết vỏ phi tuyến thế giới.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình kết hợp kỹ thuật san đều dầm trượt Lekhnitskiy với lý thuyết biến dạng trượt bậc cao HSDT và tiếp cận giải tích hàm chuyển vị trực tiếp cho các bài toán vỏ thoải đối xứng trục và vỏ tròn xoay tròn biến dạng phức tạp.
  • Về mặt thực tiễn kỹ thuật: Cung cấp bộ công cụ tính toán, biểu đồ thiết kế và tiêu chuẩn kiểm tra an toàn chịu tải xung cơ - nhiệt cực hạn cho các kỹ sư kết cấu hàng không vũ trụ (vỏ tên lửa, thân tàu con thoi), công trình ngầm chịu áp lực lớn và bồn bể chứa khí hóa lỏng nhiệt độ cao.
  • Về mặt chính sách và công nghệ vật liệu: Định hướng chiến lược phát triển công nghệ sản xuất vật liệu composite phân lớp chức năng FGM tại Việt Nam, phục vụ trực tiếp cho ngành công nghiệp quốc phòng và giao thông vận tải hiện đại.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật xuất sắc, luận án vẫn tồn tại một số giới hạn biên khoa học:

  1. Giả thiết liên kết hoàn hảo giữa các lớp (Perfect Bonding Interface): Luận án giả thiết các lớp FGM và lõi liên kết tuyệt đối lý tưởng, chưa xét đến hiện tượng phá hủy do trượt dính cục bộ (interlaminar delamination) hoặc nứt tế vi dưới chu kỳ tải nhiệt lặp (thermal fatigue).
  2. Bỏ qua biến dạng xoắn thuần túy của gân gia cường: Áp dụng mô hình dầm biến dạng trượt cho gân FGM nhưng bỏ qua thành phần độ cứng xoắn tự do St. Venant của gân mảnh.
  3. Hệ số Poisson giả thiết là hằng số: Giả thiết $\nu = \text{const}$ dọc theo chiều dày để đơn giản hóa quá trình lấy tích phân giải tích, mặc dù trong thực tế hệ số Poisson có sự biến thiên nhỏ giữa pha kim loại và gốm.

Chương trình nghiên cứu phát triển 5-10 năm tới (Future Research Agenda):

  • Phát triển mô hình phần tử hữu hạn phi đẳng hình học (Isogeometric Analysis - IGA) kết hợp cơ học phá hủy liên tục (Phase-field fracture) để mô phỏng sự bong tách lớp và lan truyền vết nứt nhiệt trong vỏ sandwich FGM.
  • Mở rộng bài toán sang vỏ thoải sandwich FGM có xét đến vật liệu áp điện (piezoelectric) hoặc vật liệu từ - điện - nhiệt (magneto-electro-elastic) phục vụ kết cấu thông minh (smart structures) có khả năng tự kiểm soát dao động chủ động.
  • Nghiên cứu thực nghiệm kiểm chứng (Experimental validation) chế tạo mẫu vỏ sandwich FGM bằng công nghệ in 3D laser kim loại - gốm (Additive Manufacturing) và thử nghiệm bàn rung động lực học.

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

Công trình luận án của tác giả Đặng Thùy Đông đã tạo nên dấu ấn học thuật sâu sắc với chuỗi 12 công trình khoa học công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (Q1/Q2) và các tạp chí chuyên ngành hàng đầu trong nước (Tạp chí Cơ học, Vietnam Journal of Mechanics).

  • Tác động học thuật (Academic Impact): Hệ thống bài báo của tác giả được trích dẫn rộng rãi bởi cộng đồng nghiên cứu cơ học tính toán quốc tế tại Mỹ, Châu Âu, Trung Quốc, Iran và Hàn Quốc, đóng góp quan trọng vào chỉ số ảnh hưởng học thuật của Đại học Quốc gia Hà Nội.
  • Chuyển giao công nghệ và R&D công nghiệp: Cung cấp cơ sở khoa học chính xác cho Viện Công nghệ Vũ trụ, Tập đoàn Công nghiệp Than - Khoáng sản, các đơn vị thiết kế kết cấu vỏ tàu biển composite và vỏ bảo vệ chịu nhiệt siêu nhẹ trong kỹ thuật hàng không.
  • Tác động đào tạo sau đại học: Đóng vai trò là tài liệu tham khảo mẫu mực về phương pháp giải tích phi tuyến và kỹ thuật xử lý bài toán cơ nhiệt cho các nghiên cứu sinh, học viên cao học chuyên ngành Cơ học vật rắn, Cơ kỹ thuật và Kỹ thuật Xây dựng/Công trình Giao thông.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ (Doctoral Researchers): Tiếp cận khung phương pháp luận giải tích mẫu mực, nắm vững kỹ thuật giải hệ phương trình đạo hàm riêng phi tuyến bằng phương pháp Galerkin, tích phân Runge-Kutta và tiêu chuẩn ổn định Budiansky-Roth.
  • Giáo sư & Chuyên gia Cơ học tính toán (Senior Academics): Sở hữu nguồn dữ liệu số chuẩn xác (benchmark datasets) để so sánh, kiểm chứng các thuật toán phần tử hữu hạn (FEM), phần tử biên (BEM) hoặc các phương pháp không lưới (meshfree methods) mới.
  • Kỹ sư trưởng R&D Công nghiệp Quốc phòng & Hàng không: Ứng dụng trực tiếp các công thức tính toán tải vồng tới hạn và tần số dao động để tối ưu hóa hình học kết cấu vỏ tên lửa, tấm chắn nhiệt turbine khí và thân tàu ngầm.
  • Cơ quan Quản lý Tiêu chuẩn & Hoạch định Công nghệ: Sử dụng dữ liệu độ tin cậy để xây dựng quy chuẩn tính toán thiết kế kết cấu composite phi truyền thống tại Việt Nam.

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? Trả lời: Đóng góp độc đáo nhất là việc thiết lập thành công mô hình giải tích vỏ thoải hai độ cong sandwich FGM đa biến thể gia cường bằng hệ thống gân FGM xiên ($\theta$) chịu nhiệt, giải quyết trọn vẹn theo Lý thuyết biến dạng trượt bậc cao (HSDT). Luận án đã mở rộng lý thuyết san đều gân của Lekhnitskiy từ dạng gân thuần nhất đẳng nhiệt sang dầm FGM dị hướng chịu nhiệt có xét đến biến dạng trượt ngang.

2. Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đó? Trả lời: Luận án đã khắc phục sự phụ thuộc vào hệ số hiệu chỉnh cắt $K_s$ của FSDT (như trong nghiên cứu của Zenkour, 2006) bằng cách phát triển hệ phương trình chủ đạo giải tích theo HSDT đầy đủ; đồng thời là công trình đầu tiên tại Việt Nam giải quyết bài toán động lực chỏm cầu thoải sandwich FGM đối xứng trục bằng phương pháp chuyển vị trực tiếp kết hợp thuật toán Runge-Kutta và tiêu chuẩn Budiansky-Roth, thay thế cho cách tiếp cận hàm ứng suất truyền thống gặp nhiều bế tắc toán học khi xét đến biến dạng trượt.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì? Trả lời: Hiện tượng phách điều hòa (harmonic beating) xuất hiện rõ nét ở chỏm cầu thoải sandwich FGM tại các vùng tần số kích thích cận cộng hưởng nhưng không dẫn đến phá hủy vồng tức thì nếu biên độ tải dưới ngưỡng tới hạn; cùng với việc góc xiên gân $\theta = 45^\circ$ tạo ra hiệu ứng tăng cứng phi đối xứng cực đại, giúp tăng tải tới hạn động lên đến $67.3%$, vượt xa hiệu quả của hệ gân trực giao kinh điển.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Hoàn toàn đầy đủ và minh bạch. Luận án cung cấp chi tiết toàn bộ các ma trận hệ số độ cứng ($A_{ij}, B_{ij}, D_{ij}, E_{ij}, F_{ij}, G_{ij}, L_{ij}$), hệ số khối lượng quán tính ($I_1 \dots I_7$), bảng thông số vật liệu phụ thuộc nhiệt độ ($P_0 \dots P_3$), các hàm dạng dao động trực giao và thuật toán Runge-Kutta, cho phép các nhà khoa học độc lập lập trình tái lập chính xác $100%$ các kết quả số.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Tập trung vào 3 trục chính: (1) Tích hợp trường đa vật lý Điện - Từ - Nhiệt - Cơ cho vỏ sandwich FGM thông minh; (2) Phát triển mô hình số vi cơ học đánh giá độ bền mỏi nhiệt và bong tách giao diện lớp; (3) Sản xuất bồi đắp in 3D mẫu thực nghiệm và chuẩn hóa quy trình thử nghiệm tải xung nhiệt siêu âm.

Kết luận

Luận án tiến sĩ của Đặng Thùy Đông là một công trình khoa học xuất sắc, mẫu mực và hoàn chỉnh về mặt học thuật trong chuyên ngành Cơ học vật rắn, đúc kết 6 đóng góp cốt lõi:

  1. Xây dựng thành công hệ phương trình chủ đạo phi tuyến phân tích ổn định tĩnh, dao động tự do và động lực học phi tuyến cho 4 cấu hình vỏ thoải sandwich FGM đa lớp chịu tải cơ - nhiệt theo HSDT và FSDT.
  2. Đề xuất và hoàn thiện mô hình gân FGM xiên đa hướng chịu nhiệt, giải quyết triệt để bài toán san đều độ cứng gân trượt trên nền đàn hồi Winkler-Pasternak.
  3. Khám phá toàn diện quy luật ứng xử động lực chỏm cầu thoải FGM bằng phương pháp chuyển vị trực tiếp, làm sáng tỏ hiện tượng phách điều hòa và cơ chế phân nhánh pha.
  4. Tiên phong mô hình hóa giải tích kết cấu vỏ trống và vỏ trụ FGM gấp nếp/lõi gấp nếp, chứng minh hiệu năng nâng cao tải trọng chịu lực từ 2 đến 3 lần.
  5. Xác định chính xác tải tới hạn động $P_{crd}$ thông qua tiêu chuẩn Budiansky-Roth hiện đại dưới các tốc độ gia tải xung biến thiên theo thời gian.
  6. Công bố quốc tế ấn tượng với 12 bài báo ISI/Scopus, khẳng định năng lực nghiên cứu đỉnh cao và mở ra các hướng phát triển bền vững cho ngành cơ học vật liệu tiên tiến tại Việt Nam trên trường quốc tế.