Tổng quan về luận án

Sự phát triển của khoa học vật liệu và kỹ thuật cơ học kết cấu hiện đại đòi hỏi các giải pháp đột phá nhằm thay thế vật liệu truyền thống trong môi trường làm việc khắc nghiệt. Vật liệu có cơ tính biến thiên (Functionally Graded Materials – FGM), được các nhà khoa học Nhật Bản khởi xướng từ năm 1984, đại diện cho bước tiến vượt bậc của composite tiên tiến nhờ sự chuyển tiếp trơn và liên tục giữa các pha thành phần gốm và kim loại. Cấu trúc biến thiên này triệt tiêu hoàn toàn hiện tượng tập trung ứng suất và bong tách lớp vốn là nhược điểm chí tử của composite phân lớp kinh điển khi vận hành dưới gradient nhiệt độ cực hạn—lên tới 2400 K cùng mức chênh lệch 1600 K qua chiều dày kết cấu mỏng dưới 10 mm.

Khoảng trống nghiên cứu (research gap) trọng tâm nằm ở tính chất phức tạp của kết cấu vỏ hở thoải hai độ cong (doubly curved shallow shell). Trong khi phần lớn y văn thế giới tập trung vào tấm phẳng hoặc vỏ kín đối xứng trục (vỏ trụ, vỏ nón), các phân tích động lực học cho vỏ thoải hai độ cong chịu tác động đồng thời của tải trọng cơ học và trường nhiệt độ còn rất hạn chế. Đặc biệt, đa số các công bố trước đây tiếp cận bài toán động nhiệt theo "Loại 1" (chỉ xét sự biến đổi cơ tính theo nhiệt độ) hoặc đơn giản hóa ứng suất nhiệt thành hằng số độc lập với điều kiện biên. Luận án khắc phục triệt để khoảng trống này bằng việc triển khai cách tiếp cận "Loại 2": nhiệt độ vừa làm suy giảm mô đun đàn hồi, vừa gây ra trạng thái biến dạng và trường ứng suất tĩnh ban đầu; kết cấu sau đó thực hiện dao động quanh vị trí cân bằng mới với ứng suất trước phụ thuộc trực tiếp vào điều kiện biên.

Nghiên cứu tập trung giải quyết ba câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Phân bố nhiệt độ dừng (tăng đều, tuyến tính, phi tuyến) theo chiều dày ảnh hưởng như thế nào đến độ võng và trạng thái ứng suất tĩnh của vỏ thoải FGM hai độ cong?
  2. RQ2: Ứng suất trước do nhiệt độ và tải trọng tĩnh tác động như thế nào đến phổ tần số dao động riêng và đáp ứng chuyển vị thời gian khi giải bài toán dao động cưỡng bức có cản?
  3. RQ3: Làm thế nào để xây dựng một mô hình số hữu hiệu vượt qua giới hạn của nghiệm giải tích nhằm phân tích vỏ FGM có hình dạng bề mặt tùy ý $z = f(x, y)$ và điều kiện biên phức tạp?

Hai giả thuyết nghiên cứu được thiết lập:

  • H1: Trường ứng suất trước nhiệt $N_{xx}^0, N_{yy}^0, N_{xy}^0$ đóng vai trò như một ma trận độ cứng hình học $[K_{ini}]$, gây ra hiệu ứng mềm hóa kết cấu (structural softening) làm suy giảm đáng kể tần số dao động riêng cơ bản.
  • H2: Phần tử vỏ 3D suy biến 8 nút với 5 bậc tự do tại mỗi nút trên mặt tham chiếu có khả năng loại bỏ hiện tượng khóa cắt (shear locking), phản ánh chính xác đáp ứng động lực học của vỏ FGM so với nghiệm giải tích giải bằng phương pháp Galerkin.

Khung lý thuyết của luận án vận dụng Lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin (FSDT) với hệ số hiệu chỉnh cắt $k_c = 5/6$, mô hình hỗn hợp Voigt theo quy luật hàm lũy thừa P-FGM, mô hình phụ thuộc nhiệt độ Touloukian, nguyên lý thế năng toàn phần cực tiểu và nguyên lý Hamilton. Phạm vi nghiên cứu bao quát vỏ thoải FGM hai độ cong có hình chiếu bằng chữ nhật kích thước $a \times b$, chiều dày $h$, bán kính cong $R_x, R_y$, làm việc dưới các điều kiện biên tựa khớp (SSSS), khớp - ngàm (SCSC), và ngàm cứng (CCCC).

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận ba dòng nghiên cứu chính trong phân tích kết cấu FGM chịu nhiệt:

Trường phái nghiên cứu Tác giả tiêu biểu & Năm Phương pháp & Mô hình lý thuyết Phạm vi & Giới hạn chính
Lý thuyết biến dạng cắt & Giải tích giải tích Reddy & Cheng (1998, 2001); Vel & Batra (2004); Zenkour (2006, 2007); Tounsi et al. (2013); Swaminathan & Sangeetha (2017) Asymptotic method, chuỗi lũy thừa, SSDT, R-TSDT, lý thuyết 4 ẩn và HSDT 12 ẩn Chủ yếu khảo sát tấm phẳng FGM, tấm sandwich; giải quyết bài toán tĩnh cơ - nhiệt dừng
Động lực học & Ứng suất nhiệt vỏ FGM Haddadpour et al. (2007); Jooybar et al. (2016); Malekzadeh & Heydarpour (2012); Kiani et al. (2012); Yang & Shen (2003) Vi phân cầu phương (DQM), lý thuyết vỏ Love/Sanders, biến đổi Laplace-Fourier Tập trung vỏ kín (vỏ trụ, vỏ nón, panel trụ); tiếp cận nhiệt phần lớn theo Loại 1
Trường phái Cơ học tính toán tại Việt Nam Hoàng Văn Tùng (2010); Nguyễn Đình Đức và cs. (2013–2017); Đào Huy Bích và cs. (2011, 2013); Đào Văn Dũng và cs. (2013, 2015) Phương pháp giải tích và bán giải tích, giải phương trình vi phân phi tuyến Tập trung vào bài toán mất ổn định tĩnh phi tuyến, sau mất ổn định và vỏ FGM có gân gia cường

Tranh luận học thuật cốt lõi tồn tại giữa hai quan điểm xử lý trường nhiệt độ trong bài toán động:

  • Quan điểm Loại 1 (Malekzadeh & Beni, 2010; Jooybar et al., 2016): Giả định nhiệt độ chỉ làm suy thoái các hằng số đàn hồi $E(z, T), \alpha(z, T), \kappa(z, T)$. Theo đó, phương trình dao động thuần túy không chứa ma trận ứng suất trước, dẫn đến việc đánh giá quá cao độ cứng động lực học của kết cấu.
  • Quan điểm Loại 2 (Yang & Shen, 2003; Kim, 2005; Pradyumna & Bandyopadhyay, 2008): Khẳng định môi trường nhiệt gây ra trường lực màng ban đầu $N_{xx}^0, N_{yy}^0$. Tuy nhiên, các tác giả này hầu hết áp dụng công thức ứng suất nhiệt cục bộ đẳng hướng $\sigma_{xx}^T = -Q_{11}\alpha(z,T)\Delta T$, bỏ qua ảnh hưởng của các điều kiện kiềm chế biên.

Luận án định vị tại điểm giao thoa giữa cơ học giải tích chuẩn tắc và phương pháp phần tử hữu hạn mở rộng. So sánh với nghiên cứu của Kiani và cộng sự (2012) vốn chỉ giới hạn ở vỏ hai độ cong tựa khớp giải bằng biến đổi tích phân, luận án mở rộng nghiệm giải tích cho nhiều dạng điều kiện biên khác nhau (SSSS, SCSC, CCCC). Đồng thời, so với mô hình giải tích phi tuyến của Đào Huy Bích và Hoàng Văn Tùng (2011), nghiên cứu này thiết lập quy trình tính toán số bằng phần tử vỏ 3D suy biến, cho phép phân tích các kết cấu vỏ thoải không có bán kính cong cố định.

Đó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 biến dạng cắt bậc nhất Reissner-Mindlin cho bài toán tương tác nhiệt - cơ động lực học của vỏ thoải FGM hai độ cong thông qua việc tích hợp ma trận ứng suất trước phụ thuộc biên. Mô hình trường chuyển vị được phân rã thành:

$$u^p(x, y, z, t) = u^t(x, y, z) + u^d(x, y, z, t)$$

Trong đó $u^t = {u_0^t, v_0^t, w_0^t, \phi_x^t, \phi_y^t}^T$ là chuyển vị tĩnh do tải trọng cơ học phân bố $q_z^t$ và tải trọng nhiệt $T(z)$ gây ra; $u^d = {u_0^d, v_0^d, w_0^d, \phi_x^d, \phi_y^d}^T$ là chuyển vị động gia tăng dao động quanh vị trí cân bằng tĩnh mới.

Bằng việc áp dụng nguyên lý Hamilton cho hệ thống có ứng suất trước:

$$\int_{t_1}^{t_2} \left( \delta T^d - \delta U^d - \delta U_0 + \delta W^d \right) dt = 0$$

Với thế năng do ứng suất ban đầu gây ra:

$$U_0 = \frac{1}{2} \int_A \int_{-h/2}^{h/2} \left[ \sigma_{xx}^0 \left(\frac{\partial w_0^d}{\partial x}\right)^2 + \sigma_{yy}^0 \left(\frac{\partial w_0^d}{\partial y}\right)^2 + 2\sigma_{xy}^0 \left(\frac{\partial w_0^d}{\partial x}\right)\left(\frac{\partial w_0^d}{\partial y}\right) \right] dz , dx , dy$$

Nghiên cứu đã chứng minh rằng trường ứng suất tĩnh $(\sigma_{xx}^0, \sigma_{yy}^0, \sigma_{xy}^0)$ thu được từ việc giải hệ 5 phương trình cân bằng tĩnh sẽ trực tiếp tham gia cấu thành ma trận độ cứng hình học $[K_{ini}]$. Đây là bằng chứng toán học xác lập bước chuyển từ tiếp cận nhiệt Loại 1 sang Loại 2 một cách hoàn chỉnh.

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

Khung phân tích của luận án là sự hợp nhất giữa ba trụ cột lý thuyết:

  1. Mô hình đồng nhất hóa vi cơ học Voigt: Xác định tính chất hiệu dụng $P_{eff}(z, T) = P_m(T) + [P_c(T) - P_m(T)] V_c(z)$ với hàm tỷ lệ thể tích $V_c(z) = (z/h + 1/2)^p$.
  2. Lý thuyết đàn nhiệt phụ thuộc nhiệt độ Touloukian: $P(T) = P_0 (1 + P_{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$ cho từng thành phần gốm và kim loại.
  3. Mô hình truyền nhiệt Fourier 1D dừng: Giải phương trình vi phân $\frac{d}{dz}\left[\kappa(z)\frac{dT}{dz}\right] = 0$ để trích xuất hàm phân bố nhiệt phi tuyến thực tế $T(z)$.

Điều kiện biên nghiên cứu được định nghĩa chặt chẽ qua các phương trình toán học giải tích:

  • SSSS (Bốn biên tựa khớp): $v_0 = w_0 = \phi_y = N_{xx} = M_{xx} = 0$ tại $x = 0, a$; $u_0 = w_0 = \phi_x = N_{yy} = M_{yy} = 0$ tại $y = 0, b$.
  • SCSC (Khớp - Ngàm đối xứng): Tựa khớp tại $x = 0, a$; ngàm cứng ($u_0 = v_0 = w_0 = \phi_x = \phi_y = 0$) tại $y = 0, b$.
  • CCCC (Bốn biên ngàm cứng): Triệt tiêu toàn bộ 5 thành phần chuyển vị trên toàn bộ đường bao.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng (Positivism) và duy thực phê phán (Critical Realism), kết hợp chặt chẽ giữa suy diễn giải tích toán học chính xác và mô phỏng số phần tử hữu hạn. Thiết kế nghiên cứu tuân thủ mô hình song hành đa phương pháp (Dual-Track Methodology):

Quy trình nghiên cứu rigorous

Quy trình giải tích áp dụng phương pháp trực giao Galerkin:

  1. Biểu diễn 5 ẩn số chuyển vị dưới dạng chuỗi lượng giác kép: $$u_0(x, y) = \sum_{m=1}^{\infty}\sum_{n=1}^{\infty} U_{mn} X_m'(x) Y_n(y); \quad w_0(x, y) = \sum_{m=1}^{\infty}\sum_{n=1}^{\infty} W_{mn} X_m(x) Y_n(y)$$
  2. Các hàm dạng không gian $X_m(x), Y_n(y)$ được lựa chọn chính xác theo từng điều kiện biên (hàm $\sin, \cos$ hoặc dạng lũy thừa kết hợp).
  3. Chiếu hệ phương trình vi phân cân bằng lên hệ hàm cơ sở thông qua tích phân diện tích: $$\int_0^b \int_0^a \text{Eq}_k \cdot \Psi_k(x, y) , dx , dy = 0 \quad (k = 1..5)$$
  4. Thu được hệ phương trình đại số tuyến tính $[K]{\Delta_{mn}^t} = {F_{ch}^t} + {F_{nd}^t}$.

Đối với phân tích động lực học cưỡng bức, ma trận cản Rayleigh $[C] = a_1 [M] + a_2 ([K] + [K_{ini}])$ được thiết lập từ hai tỷ số cản $\xi_1, \xi_2$ ứng với hai tần số riêng thấp nhất $\omega_1, \omega_2$:

$$a_1 = \frac{2\omega_1\omega_2(\xi_1\omega_2 - \xi_2\omega_1)}{\omega_2^2 - \omega_1^2}; \quad a_2 = \frac{2(\xi_2\omega_2 - \xi_1\omega_1)}{\omega_2^2 - \omega_1^2}$$

Hệ phương trình vi phân chuyển động cấp 2 được giải từng bước thời gian bằng thuật toán Newmark-$\beta$ với tham số gia tốc trung bình ($\gamma = 0.5, \beta = 0.25$), đảm bảo tính ổn định vô điều kiện.

Xây dựng mô hình PTHH sử dụng phần tử vỏ 3D suy biến (3D-Degenerated shell element) 8 nút:

  • Rút gọn từ phần tử khối 3D về mặt phẳng trung bình, giảm từ 24 bậc tự do xuống còn 5 bậc tự do/nút (3 chuyển vị tổng thể $u, v, w$ và 2 góc xoay vector chỉ phương $\mathbf{V}_1, \mathbf{V}_2$).
  • Thiết lập 4 hệ tọa độ chuyển đổi: Hệ tọa độ tổng thể $xyz$, hệ tọa độ địa phương $x'y'z'$, hệ tọa độ nút $\mathbf{V}_1\mathbf{V}_2\mathbf{V}_3$ và hệ tọa độ tự nhiên $\xi, \eta, \zeta \in [-1, 1]$.
  • Hàm dạng đẳng tham số bậc 2 Serendipity $N_i(\xi, \eta)$ được áp dụng với tích phân số Gauss $2 \times 2$ trong mặt phẳng và quy tắc Simpson theo phương chiều dày $\zeta$.
         η

Data và phân tích

Phân tích số sử dụng hệ vật liệu kim loại - gốm điển hình:

  • Ti-6Al-4V (Kim loại): $E_m(T) = 122.56 \times 10^9 (1 - 4.586 \times 10^{-4} T) \text{ Pa}$, $\alpha_m(T) = 7.5788 \times 10^{-6} (1 + 6.638 \times 10^{-4} T) \text{ K}^{-1}$, $\rho_m = 4429 \text{ kg/m}^3$, $\kappa_m = 1.0 \text{ W/m}\cdot\text{K}$.
  • Zirconia / $\text{ZrO}_2$ (Gốm): $E_c(T) = 244.27 \times 10^9 (1 - 3.270 \times 10^{-4} T) \text{ Pa}$, $\alpha_c(T) = 12.766 \times 10^{-6} (1 - 1.491 \times 10^{-4} T) \text{ K}^{-1}$, $\rho_c = 3000 \text{ kg/m}^3$, $\kappa_c = 2.09 \text{ W/m}\cdot\text{K}$.
  • Silicon Nitride / $\text{Si}_3\text{N}_4$ & SUS304: Sử dụng bộ thông số phụ thuộc nhiệt độ theo thực nghiệm của Touloukian.

Tác giả đã lập trình 03 bộ chương trình tính toán độc lập trên môi trường MATLAB:

  1. ShellpanelStatic(GT): Giải bài toán tĩnh, tính toán độ võng $w_0^t(x, y)$ và các thành phần ứng suất $\sigma_{xx}^t, \sigma_{yy}^t, \sigma_{xy}^t, \sigma_{xz}^t, \sigma_{yz}^t$.
  2. ShellpanelVibration(GT): Thiết lập ma trận $[K_{ini}]$, giải bài toán trị riêng suy rộng xác định tần số riêng $\omega_{mn}$ và dạng mode tương ứng.
  3. ShellpanelForcedvibration(GT): Tích phân Newmark-$\beta$ tính toán đáp ứng quá độ và đáp ứng dừng dưới kích thích điều hòa $F_0 \sin(\Omega t)$ hoặc tải trọng xung.

Quy trình kiểm chứng (robustness check) thực hiện qua kiểm tra độ hội tụ lưới phần tử ($4 \times 4, 8 \times 8, 16 \times 16, 20 \times 20$) và đối chiếu chéo giữa nghiệm giải tích và PTHH, ghi nhận sai số tương đối luôn nhỏ hơn $1.2%$.

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

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

  1. Hiệu ứng mềm hóa kết cấu do ứng suất nhiệt trước: Khi độ chênh nhiệt $\Delta T = T_c - T_m$ tăng từ $0\text{ K}$ lên $500\text{ K}$, tần số dao động riêng cơ bản $\omega_1$ của vỏ giảm từ $22%$ đến $38.5%$ tùy thuộc vào điều kiện kiềm chế biên. Điều này chứng minh rằng việc bỏ qua ma trận $[K_{ini}]$ trong cách tiếp cận Loại 1 dẫn đến sự sai lệch nguy hiểm trong đánh giá an toàn kết cấu.
  2. Quy luật phân bố ứng suất nhiệt phức hợp: Do sự bất đối xứng của cơ tính theo chiều dày, ứng suất màng nhiệt không triệt tiêu ngay cả khi không có tải trọng cơ học. Bề mặt giàu gốm chịu ứng suất nén nhiệt rất cao trong khi bề mặt giàu kim loại chịu ứng suất kéo. Biên ngàm (CCCC) làm tăng ứng suất nén nhiệt lên gấp $2.4$ lần so với biên tựa khớp (SSSS).
  3. Ưu thế của quy luật truyền nhiệt phi tuyến: Phân tích chỉ ra rằng quy luật truyền nhiệt tuyến tính hoặc tăng nhiệt đều đánh giá quá cao độ võng tĩnh (sai lệch $12 - 18%$). Nghiệm phương trình truyền nhiệt Fourier phi tuyến phản ánh chính xác hiệu ứng cách nhiệt của lớp gốm bề mặt, làm giảm ứng suất nhiệt cực đại $14.6%$.
  4. Tác động của chỉ số tỷ lệ thể tích $p$ và độ cong: Khi $p$ tăng từ $0$ (hoàn toàn gốm) đến $\infty$ (hoàn toàn kim loại), độ võng tĩnh tăng dần do mô đun đàn hồi giảm. Độ cong hai chiều ($R_x, R_y$) tạo hiệu ứng vòm (arch action), giúp vỏ cầu ($R_x = R_y$) có độ cứng uốn cao hơn đáng kể và độ võng giảm tới $41.2%$ so với tấm phẳng FGM có cùng kích thước.
  5. Hiện tượng cộng hưởng động lực học phụ thuộc nhiệt: Dưới tải trọng cưỡng bức điều hòa, tỷ số tần số $\Omega/\omega_1 = 1.0$ gây ra hiện tượng cộng hưởng với biên độ chuyển vị tăng vọt. Khi nhiệt độ môi trường tăng, do tần số riêng $\omega_1$ bị suy giảm, đỉnh cộng hưởng bị dịch chuyển về phía vùng tần số kích thích thấp hơn.

Implications đa chiều

  • Tiến bộ lý thuyết: Đặt nền móng hoàn chỉnh cho lý thuyết dao động kết cấu composite FGM có xét đến lịch sử chịu tải nhiệt - cơ tĩnh trước đó.
  • Đổi mới phương pháp luận: Khẳng định tính hiệu quả của phần tử 3D suy biến trong việc giải quyết các bài toán vỏ thoải có phương trình mặt cong $z = f(x, y)$ bất kỳ mà phương pháp giải tích không thể tiếp cận.
  • Ứng dụng công nghiệp: Cung cấp bộ công cụ tính toán và bảng tra cứu số liệu chuẩn xác cho kỹ sư hàng không vũ trụ khi thiết kế vỏ bảo vệ mũi tên lửa, ống xả tuabin khí và tấm chắn nhiệt lò phản ứng hạt nhân.
  • Khuyến nghị thiết kế: Khuyến cáo bắt buộc phải tính toán kiểm tra tần số dao động riêng trong điều kiện làm việc ở nhiệt độ cao thực tế để tránh hiện tượng cộng hưởng ngoài dự kiến do tần số riêng bị suy giảm.

Limitations và Future Research

Nghiên cứu thừa nhận một số giới hạn khoa học mang tính tiền đề:

  1. Giả thiết biến dạng nhỏ và quan hệ ứng suất - biến dạng tuyến tính đàn hồi, chưa bao hàm hiệu ứng phi tuyến hình học von Kármán lớn hoặc phi tuyến vật liệu (chảy dẻo nhiệt).
  2. Hệ số Poisson $\nu$ được coi là hằng số theo chiều dày ($\nu = \text{const} \approx 0.3$) để giảm bớt độ phức tạp giải tích, mặc dù sai số thực tế dưới $2%$.
  3. Giới hạn trong trường nhiệt độ tĩnh dừng 1D $T(z)$, chưa mở rộng cho bài toán truyền nhiệt không dừng phụ thuộc thời gian ($\partial T/\partial t \ne 0$) hoặc bài toán sốc nhiệt ba chiều.
  4. Chưa tích hợp ảnh hưởng của độ xốp (porosity) phát sinh trong quá trình thiêu kết vật liệu FGM thực tế.

Chương trình nghiên cứu tương lai định hướng 4 trọng tâm:

  • Phát triển mô hình phi tuyến hình học (Large Deflection) kết hợp lý thuyết biến dạng cắt bậc cao (HSDT/SSDT) cho vỏ thoải FGM.
  • Khảo sát đáp ứng động lực học của kết cấu vỏ FGM xốp (porous FGM) gia cường ống nano carbon (CNTRC) trong môi trường nhiệt ẩm (hygrothermal).
  • Nghiên cứu bài toán tương tác nhiệt - cơ - điện chủ động bằng cách dán các lớp áp điện (Piezoelectric FGM) để điều khiển dao động.
  • Tối ưu hóa đa mục tiêu quy luật phân bố hàm tỷ lệ thể tích $V_c(z)$ nhằm giảm thiểu trọng lượng và tối đa hóa tần số dao động riêng.

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

  • Tác động học thuật: Luận án đóng góp 07 công trình khoa học công bố trên các tạp chí và kỷ yếu hội nghị chuyên ngành uy tín trong nước; dự kiến đóng góp nguồn tài liệu tham khảo cho các nghiên cứu tiếp theo về cơ học kết cấu FGM tại Việt Nam và khu vực.
  • Chuyển dịch công nghệ: Mở ra khả năng làm chủ công nghệ tính toán thiết kế các cấu kiện chịu nhiệt độ cao trong ngành công nghiệp quốc phòng, sản xuất động cơ phản lực và phương tiện bay siêu thanh.
  • Lợi ích kinh tế - kỹ thuật: Tối ưu hóa quá trình chọn lựa vật liệu gốm/kim loại ($p, h, R_x, R_y$), giúp kéo dài tuổi thọ kết cấu thêm $25 - 30%$ và giảm thiểu nguy cơ phá hủy do mỏi nhiệt.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận khung lý thuyết giải tích mẫu mực và bộ chương trình MATLAB hoàn chỉnh phục vụ nghiên cứu phát triển.
  • Kỹ sư R&D Hàng không - Vũ trụ: Ứng dụng thuật toán và mô hình phần tử 3D suy biến vào các phần mềm mô phỏng kỹ thuật chuyên dụng (CAE/FEA).
  • Ngành Công nghiệp Năng lượng & Hạt nhân: Sử dụng cơ sở dữ liệu về ứng suất nhiệt để thiết kế vỏ lò phản ứng và đường ống dẫn khí nóng an toàn tuyệt đố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ì?

Đóng góp độc đáo nhất là việc thiết lập thành công mô hình toán học giải tích cho bài toán động lực học vỏ thoải FGM hai độ cong theo cách tiếp cận Loại 2. Khác với các nghiên cứu truyền thống xem ứng suất nhiệt là cục bộ và bất biến, luận án xác định trường ứng suất nhiệt trước $(\sigma_{xx}^0, \sigma_{yy}^0, \sigma_{xy}^0)$ trực tiếp từ nghiệm của hệ 5 phương trình vi phân cân bằng tĩnh có kể đến điều kiện biên. Trường ứng suất này được tích hợp chính xác vào thế năng ban đầu $U_0$, hình thành nên ma trận độ cứng hình học $[K_{ini}]$ trong phương trình chuyển động Hamilton.

2. Sự đổi mới trong phương pháp nghiên cứu so với các công trình quốc tế tiền nhiệm?

So với Kiani và cs. (2012) vốn dùng biến đổi tích phân chỉ áp dụng được cho biên khớp (SSSS), luận án đã dùng phương pháp Galerkin mở rộng với hệ hàm lượng giác chuyên biệt giải quyết được biên hỗn hợp phức tạp (SCSC, CCCC). So với mô hình PTHH tấm phẳng của Bhandari & Purohit (2014), luận án ứng dụng phần tử vỏ 3D suy biến 8 nút (Ahmad et al.), tích hợp vector chỉ phương theo chiều dày $\mathbf{V}_3$, loại bỏ hoàn toàn hiện tượng khóa cắt mà không cần tăng số lượng nút tính toán quá mức như phần tử khối 3D (Apalak, 2010 với 2500 lớp).

3. Phát hiện thực nghiệm/tính toán nào gây bất ngờ nhất?

Phát hiện bất ngờ nhất là hiện tượng suy giảm đột ngột của tần số dao động riêng dưới điều kiện biên ngàm cứng (CCCC) khi nhiệt độ tăng. Thông thường, biên ngàm làm tăng độ cứng kết cấu so với biên khớp. Tuy nhiên trong môi trường nhiệt, biên ngàm kiềm chế hoàn toàn sự dãn nở nhiệt, sinh ra ứng suất nén nội tại khổng lồ, khiến ma trận $[K_{ini}]$ mang giá trị âm vượt trội, làm triệt tiêu độ cứng uốn tổng thể nhanh hơn nhiều so với biên tựa khớp.

4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án công khai toàn bộ hệ thức giải tích vi phân (từ quan hệ biến dạng - chuyển vị, phương trình cân bằng tĩnh, hệ phương trình chuyển động động học), định nghĩa rõ ràng các tích phân độ cứng $A_{ij}, B_{ij}, D_{ij}$, lực nhiệt $N_{ij}^{nd}, M_{ij}^{nd}$, cùng lưu đồ thuật toán 4 bước chi tiết của 03 module MATLAB (ShellpanelStatic, ShellpanelVibration, ShellpanelForcedvibration).

5. Chương trình nghị sự nghiên cứu 10 năm được phác thảo ra sao?

Lộ trình 10 năm gồm 3 giai đoạn:

  • Giai đoạn 1 (1-3 năm): Mở rộng nghiệm cho bài toán phi tuyến hình học (Large Deflection) và ổn định động lực học phi tuyến.
  • Giai đoạn 2 (4-6 năm): Mô hình hóa vỏ FGM có lỗ xốp, composite gia cường hạt nano trong điều kiện nhiệt không dừng và sốc nhiệt động lực.
  • Giai đoạn 3 (7-10 năm): Tích hợp vật liệu thông minh (Smart FGM/Piezoelectric), xây dựng phần mềm phân tích kết cấu nhiệt - cơ tối ưu hóa phục vụ công nghiệp chế tạo thiết bị bay không người lái (UAV) và tên lửa đẩy.

Kết luận

Luận án "Phân tích tĩnh và động kết cấu vỏ thoải FGM hai độ cong trong môi trường nhiệt" đã đạt được các đóng góp then chốt:

  1. Thiết lập thành công hệ phương trình chủ đạo và nghiệm giải tích giải bằng phương pháp Galerkin cho bài toán tĩnh, dao động riêng và dao động cưỡng bức của vỏ thoải FGM hai độ cong trong môi trường nhiệt độ với các điều kiện biên SSSS, SCSC, CCCC.
  2. Hiện thực hóa trọn vẹn cách tiếp cận động lực học nhiệt Loại 2, đưa ma trận ứng suất trước do nhiệt độ phụ thuộc điều kiện biên vào phương trình chuyển động của vỏ.
  3. Xây dựng mô hình và thuật toán phần tử hữu hạn sử dụng phần tử vỏ 3D suy biến 8 nút, mở rộng khả năng tính toán cho các kết cấu vỏ có dạng bề mặt toán học $z = f(x, y)$ bất kỳ.
  4. Lập trình thành công bộ 03 chương trình tính toán chuyên dụng trên nền tảng MATLAB (ShellpanelStatic, ShellpanelVibration, ShellpanelForcedvibration), khảo sát toàn diện ảnh hưởng của các tham số vật liệu $p$, kích thước hình học $a/b, a/h, R_x, R_y$, quy luật truyền nhiệt và tỷ số cản Rayleigh.
  5. Chứng minh bản chất của hiệu ứng mềm hóa kết cấu do nhiệt độ và xác lập quy luật suy giảm tần số dao động riêng dưới tác động của trường nhiệt độ phi tuyến.
  6. Mở ra 3 nhánh nghiên cứu chuyên sâu về vỏ composite chức năng thông minh, phi tuyến hình học và điều khiển chủ động, đóng góp thiết thực cho sự phát triển của cơ học tính toán và kỹ thuật hàng không vũ trụ hiện đại.