Tổng quan về luận án

Nghiên cứu về độ bền và ổn định của kết cấu thành mỏng hoạt động trong điều kiện khắc nghiệt là trọng tâm của cơ học vật rắn biến dạng hiện đại. Luận án tiến sĩ cơ học của tác giả Lê Khả Hòa (2015) với đề tài "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" (chuyên ngành Cơ học vật rắn, mã số 62440107, do PGS.TS Đào Văn Dũng hướng dẫn tại Đại học Quốc gia Hà Nội) đại diện cho bước tiến giải tích trong việc dự báo ứng xử cơ học phi tuyến của kết cấu vỏ mỏng chế tạo từ vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM).

Trong bối cảnh vật liệu FGM được ứng dụng ngày càng rộng rãi trong công nghiệp hàng không vũ trụ (bộ phận tên lửa, khung thân máy bay siêu thanh), năng lượng hạt nhân (tường plasma lò phản ứng nhiệt hạch), việc loại bỏ sự tập trung ứng suất tiếp xúc và hiện tượng bong tách lớp (delamination) đặc trưng của composite lớp truyền thống là ưu điểm nổi bật: "Do các tính chất cơ lý biến đổi trơn và liên tục từ mặt này đến mặt kia nên các kết cấu FGM tránh được sự tập trung ứng suất trên bề mặt tiếp xúc giữa các lớp, tránh được sự bong tách và rạn nứt trong kết cấu." Tuy nhiên, khoảng trống nghiên cứu (research gap) lớn của tài liệu khoa học quốc tế trước thời điểm này là phần lớn các công bố giải tích chỉ giới hạn ở các dạng kết cấu đơn giản (thanh, tấm phẳng, vỏ thoải) không có gân gia cường và thường đơn giản hóa coi hệ số Poisson là hằng số ($\nu = \text{const}$). Các kết cấu vỏ phức tạp hơn như panel trụ, vỏ trụ tròn chịu xoắn và áp lực ngoài có gân FGM gia cường lệch tâm (ES-FGM), hoặc vỏ nón cụt FGM đặt trên nền đàn hồi Winkler - Pasternak hầu như chưa có lời giải giải tích tường minh hoàn chỉnh.

Luận án giải quyết 5 câu hỏi và giả thuyết khoa học cốt lõi:

  1. Câu hỏi nghiên cứu 1 (Q1): Sự phụ thuộc theo chiều dày của hệ số Poisson $\nu = \nu(z)$ ảnh hưởng như thế nào đến ma trận độ cứng giải tích và đường cong sau tới hạn của panel và vỏ trụ FGM? $\rightarrow$ Giả thuyết 1 (H1): Việc bỏ qua tính biến thiên của $\nu(z)$ dẫn đến sai lệch đáng kể trong việc xác định tải vồng cận trên ($r_{0\text{upper}}$) và cận dưới ($r_{0\text{lower}}$).
  2. Câu hỏi nghiên cứu 2 (Q2): Làm thế nào để thiết lập hệ phương trình vi phân phi tuyến giải tích cho panel trụ FGM không hoàn hảo với các điều kiện biên hỗn hợp (tựa đơn 4 cạnh vs. hai cạnh cong tựa đơn, hai cạnh thẳng ngàm trượt)? $\rightarrow$ Giả thuyết 2 (H2): Bằng việc sử dụng hàm ứng suất Airy và nghiệm hàm độ võng thích hợp, phương pháp Galerkin cho phép tìm được nghiệm giải tích tường minh xác định mối quan hệ tải - độ võng.
  3. Câu hỏi nghiên cứu 3 (Q3): Tác động tương hỗ giữa hệ gân gia cường lệch tâm FGM (ES-FGM) và vỏ FGM dưới tải áp lực ngoài và tải xoắn thuần túy diễn ra theo quy luật nào? $\rightarrow$ Giả thuyết 3 (H3): Gân FGM lệch tâm kết hợp với phương pháp san đều tác dụng gân Lekhnitskii sẽ triệt tiêu hiệu ứng tập trung ứng suất nhiệt - cơ và nâng cao rõ rệt tải tới hạn xoắn $\tau_{\text{cr}}$ và áp lực ngoài $q_{\text{cr}}$.
  4. Câu hỏi nghiên cứu 4 (Q4): Nền đàn hồi hai tham số Winkler - Pasternak ($K_1, K_2$) đóng vai trò định lượng thế nào trong việc kìm hãm sự mất ổn định hình học của vỏ trụ và vỏ nón FGM có gân? $\rightarrow$ Giả thuyết 4 (H4): Tham số lớp cắt Pasternak $K_2$ có hiệu quả gia cường độ cứng chống vồng vượt trội so với tham số nhún Winkler $K_1$.
  5. Câu hỏi nghiên cứu 5 (Q5): Có thể xây dựng một hệ thống công thức giải tích hiển cho bài toán ổn định tuyến tính của vỏ nón cụt FGM có gân với hệ số vi phân phụ thuộc biến không gian hay không? $\rightarrow$ Giả thuyết 5 (H5): Phép biến đổi tọa độ kết hợp phương pháp Galerkin có thể giải quyết triệt để hệ ba phương trình đạo hàm riêng có hệ số là hàm biến thiên theo tọa độ trục.

Khung lý thuyết của nghiên cứu dựa trên lý thuyết vỏ mỏng phi tuyến Donnell - von Kármán, kỹ thuật san đều gân Lekhnitskii cải tiến cho vật liệu FGM, và mô hình tương tác môi trường đàn hồi hai tham số Pasternak. Phạm vi khảo sát trải rộng trên nhiều hệ thông số hình học và vật liệu: tỉ số bán kính trên chiều dày $R/h \in [100, 200]$, tỉ số kích thước mặt cắt $a/b \in [0.605, 1.5]$, $b/h \in [30, 75]$, góc nửa đỉnh nón $\alpha \in [10^\circ, 60^\circ]$, chỉ số tỉ phần thể tích $k \in [0, 5]$, và biên độ không hoàn hảo hình học ban đầu $\xi \in [0, 0.5]$ trên các tổ hợp vật liệu $\text{Si}_3\text{N}_4/\text{SUS304}$ ($E_c = 322.4 \text{ GPa}, E_m = 207.788 \text{ GPa}, \nu_c = 0.24, \nu_m = 0.3177$) và $\text{ZrO}_2/\text{Ti-6Al-4V}$.

Literature Review và Positioning

Lịch sử nghiên cứu ứng xử ổn định của kết cấu thành mỏng FGM được định hình bởi nhiều trường phái học thuật quốc tế và trong nước. Trường phái sử dụng lý thuyết biến dạng trượt bậc cao (HSDT) và phương pháp tham số bé được dẫn dắt bởi Shen (2002, 2007) đã phân tích trạng thái sau tới hạn của vỏ trụ FGM chịu nén trục, áp lực ngoài và xoắn trong môi trường nhiệt độ. Song song đó, Huang và Han (2008, 2009, 2010) áp dụng lý thuyết biến dạng lớn phi tuyến kết hợp phương pháp biến phân Ritz để khảo sát vỏ trụ FGM không gân. Tuy nhiên, một hạn chế lớn trong các công trình của Huang và Han là các hệ số độ cứng $A_{ij}$ vẫn phải duy trì dưới dạng tích phân số, gây khó khăn cho việc phân tích độ nhạy giải tích trực tiếp. Trong lĩnh vực vỏ nón FGM, Sofiyev (2005, 2010) và Sofiyev & Schnack (2004) đã triển khai phương pháp Galerkin dựa trên phương trình Donnell cải tiến để xác định tải tới hạn tuyến tính của vỏ nón cụt FGM không gân chịu nén và tải xoắn động.

Tại Việt Nam, nhóm nghiên cứu của Đào Huy Bích và cộng sự (2011, 2013) đã tiên phong phát triển phương pháp giải tích phân tích ổn định tĩnh và động của tấm và vỏ thoải FGM có gân gia cường lệch tâm bằng kỹ thuật san gân. Tiếp đó, Nguyễn Đình Đức và Hoàng Văn Tùng (2011) đã nghiên cứu ổn định phi tuyến của panel trụ FGM không hoàn hảo chịu nén dọc trục nhưng chỉ giới hạn trong trường hợp hệ số Poisson không đổi ($\nu = \text{const}$) và điều kiện biên bốn cạnh tựa đơn. Đào Văn Dũng và các cộng sự (2012, 2013, 2014) đã mở rộng nghiên cứu ổn định tĩnh của tấm và panel trụ FGM có gân lệch tâm trên nền đàn hồi.

Trong bức tranh tổng quan đó tồn tại hai luồng quan điểm và tranh luận khoa học lớn:

  • Tranh luận về ảnh hưởng của hệ số Poisson $\nu(z)$: Phần lớn các nghiên cứu quốc tế bỏ qua sự thay đổi của hệ số Poisson theo chiều dày để đơn giản hóa hệ phương trình vi phân, viện dẫn rằng ảnh hưởng của $\nu$ là thứ yếu so với mô đun Young $E(z)$. Ngược lại, quan điểm của luận án chứng minh rằng việc xem $\nu(z)$ biến đổi theo hàm mũ hoặc lũy thừa làm xuất hiện các số hạng liên kết phức tạp trong ma trận độ cứng $A_{ij}$, tạo ra sự thay đổi định lượng rõ ràng lên giá trị tải tới hạn rẽ nhánh và điểm cực trị sau vồng.
  • Tranh luận về phương pháp xử lý hệ gân gia cường: Việc mô hình hóa gân FGM gia cường lệch tâm theo lý thuyết phần tử hữu hạn rời rạc đòi hỏi chi phí tính toán rất lớn và khó biểu diễn bản chất vật lý tổng quát. Luận án đã định vị đóng góp bằng cách kế thừa và nâng tầm kỹ thuật san đều gân Lekhnitskii, mở rộng từ gân thuần nhất sang gân FGM lệch tâm có cơ tính biến thiên đồng bộ với vỏ.

So sánh với các nghiên cứu kinh điển quốc tế:

  1. So sánh với nghiên cứu của Turvey (1977) và Shen (2002) cho panel trụ thuần nhất chịu nén trục ($a/b=0.605, b/h=75$), kết quả không thứ nguyên $N_x/(ER)$ của luận án đạt $0.60523 \times 10^{-4}$ tại mode vồng $(m,n)=(2,2)$, hoàn toàn trùng khớp với dữ liệu của Turvey và Shen.
  2. So sánh với lời giải giải tích của Brush và Almroth (1975) cho vỏ trụ tròn thuần nhất chịu áp lực ngoài, sai số tương đối của mô hình giải tích trong luận án là dưới $0.15%$.
  3. So sánh với kết quả của Seide (1962) và Sofiyev (2010) đối với vỏ nón cụt thuần nhất chịu nén dọc trục, độ lệch tải tới hạn tới hạn của luận án hoàn toàn nằm trong giới hạn hội tụ giải tích cho phép.

Đó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 khung lý thuyết vỏ mỏng phi tuyến Donnell - von Kármán bằng việc thiết lập tường minh ma trận độ cứng phụ thuộc biến không gian $z$ khi cả $E(z)$ và $\nu(z)$ cùng biến thiên theo luật lũy thừa và luật hàm mũ: $$E(z) = E_c e^{\frac{1}{2}\ln\left(\frac{E_c}{E_m}\right)\left(1+\frac{2z}{h}\right)}, \quad \nu(z) = \nu_c e^{\frac{1}{2}\ln\left(\frac{\nu_c}{\nu_m}\right)\left(1+\frac{2z}{h}\right)}$$

Đóng góp lý thuyết then chốt được tác giả khẳng định: "Các biểu thức hiển này giúp cho việc lập trình tính toán trở nên dễ dàng hơn việc sử dụng phương pháp số để tính các hệ số độ cứng $A_{ij}$ từ công thức (2.7) mà các tài liệu khác vẫn làm. Đây cũng là một phần đóng góp mới của luận án."

Khung lý thuyết thiết lập 3 mệnh đề cơ học cốt lõi:

  • Mệnh đề 1 (Về tính bất đối xứng vật liệu): Sự biến thiên của mô đun đàn hồi và hệ số Poisson theo chiều dày phá vỡ tính đối xứng của trường ứng suất qua mặt giữa, làm xuất hiện các thành phần liên kết uốn - kéo nén $A_{11}, A_{21}, A_{31}$ khác không, dẫn đến hiện tượng uốn khớp nối ngay cả khi kết cấu chỉ chịu lực nén thuần túy trong mặt phẳng.
  • Mệnh đề 2 (Về điểm rẽ nhánh và điểm cực trị): Vỏ FGM hoàn hảo chịu mất ổn định theo kiểu rẽ nhánh (bifurcation buckling) với tải vồng trên xác định từ trạng thái cân bằng lân cận tuyến tính, trong khi vỏ FGM không hoàn hảo hình học ($\xi > 0$) mất ổn định theo kiểu điểm giới hạn cực trị (snap-through limit point buckling) với điều kiện giới hạn triệt tiêu đạo hàm bậc nhất $dp/dw = 0$.
  • Mệnh đề 3 (Về tính tương đương cơ học của gân FGM): Hệ gân FGM bố trí lệch tâm (ES-FGM) có thể được quy đổi giải tích thành các lớp vật liệu bổ sung có độ cứng tương đương thông qua tích phân mặt cắt mở rộng Lekhnitskii, bảo toàn quan hệ ứng suất - biến dạng phi tuyến toàn cục.

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

Khung phân tích của luận án tích hợp đồng bộ ba lý thuyết nền tảng:

  1. Lý thuyết vỏ mỏng Donnell với quan hệ phi tuyến hình học von Kármán mô tả trường chuyển vị lớn và độ võng ban đầu $w^*(x,y)$.
  2. Kỹ thuật đồng nhất hóa giải tích san đều gân gia cường Lekhnitskii tích hợp cơ tính FGM.
  3. Mô hình tương tác cơ học môi trường đàn hồi hai tham số Winkler - Pasternak ($q_e = K_1 w - K_2 \nabla^2 w$).

Hệ phương trình chủ đạo phi tuyến được dẫn xuất thành công về hệ hai phương trình tương thích biến dạng và cân bằng phi tuyến liên kết giữa độ võng $w(x,y)$ và hàm ứng suất Airy $\Phi(x,y)$: $$\nabla^4 \Phi = C_1 \nabla^4 w + C_2 \left( w_{,xy}^2 - w_{,xx}w_{,yy} - \frac{1}{R}w_{,xx} - 2w_{,xy}w_{,xy}^* + w_{,xx}w_{,yy}^* + w_{,yy}w_{,xx}^* \right)$$ $$\frac{1}{R}\Phi_{,xx} - C_3 \nabla^4 \Phi + C_4 \nabla^4 w + L(w, \Phi) + q_0 = 0$$

Điều kiện biên được giải quyết cho hai trường hợp phức tạp: (i) Bốn cạnh tựa đơn bản lề ($w=0, M_x=0, N_{xy}=0, N_x=-r_0 h$); (ii) Hai cạnh cong tựa đơn và hai cạnh thẳng ngàm trượt ($w=0, w_{,y}=0, N_y=0, N_{xy}=0$), kết hợp hàm độ võng ba số hạng đảm bảo hội tụ chính xác hình thái sau vồng.

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

Thiết kế nghiên cứu

Nghiên cứu được định vị vững chắc trên hệ hình nhận thức luận thực chứng (positivist paradigm) và cơ học môi trường liên tục tất định (deterministic continuum mechanics). Thiết kế nghiên cứu đa tầng (multi-level analytical framework) liên kết chặt chẽ ba cấp độ cấu trúc:

  • Cấp độ vi mô/trung mô: Quy tắc hỗn hợp thể tích (Voigt rule of mixtures) và quy luật hàm mũ/lũy thừa mô tả sự biến thiên trơn của các pha gốm/kim loại.
  • Cấp độ phần tử cấu trúc: Tích phân liên tục tensor ứng suất và mômen qua chiều dày $h$ để xác lập ma trận độ cứng giải tích $A_{ij}, B_{ij}, D_{ij}$.
  • Cấp độ hệ thống kết cấu: Thiết lập bài toán giá trị biên phi tuyến cho panel trụ, vỏ trụ tròn và vỏ nón cụt có gân gia cường.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu phân tích giải tích tuân thủ các bước nghiêm ngặt:

  1. Mô hình hóa hình học và biến dạng: Biểu diễn trường biến dạng mặt giữa $\epsilon_x^0, \epsilon_y^0, \gamma_{xy}^0$ có kể đến độ cong ban đầu $1/R$ và độ lệch hình học $w^*(x,y) = \xi h \sin(m\pi x/a)\sin(n\pi y/b)$.
  2. Khử ẩn số chuyển vị qua hàm ứng suất Airy: Sử dụng định nghĩa $N_x = \Phi_{,yy}$, $N_y = \Phi_{,xx}$, $N_{xy} = -\Phi_{,xy}$ để thỏa mãn tự động hai phương trình cân bằng trong mặt phẳng.
  3. Chiếu Galerkin trực giao: Tác động toán tử Galerkin trên miền xác định không gian $\Omega$: $$\int_0^a \int_0^b \mathcal{R}(w, \Phi) \sin\left(\frac{m\pi x}{a}\right)\sin\left(\frac{n\pi y}{b}\right) dx dy = 0$$ chuyển đổi hệ phương trình vi phân phi tuyến đạo hàm riêng thành hệ phương trình đại số phi tuyến bậc cao đối với biên độ độ võng $W/h$.
  4. Trích xuất nghiệm giải tích tường minh: Khử ẩn số hàm ứng suất để thu được công thức hiển trực tiếp cho tải vồng cận trên $r_{0\text{upper}}$, tải vồng cận dưới $r_{0\text{lower}}$, áp lực ngoài tới hạn $q_{\text{cr}}$, và mômen xoắn tới hạn $\tau_{\text{cr}}$.

Data và phân tích

Dữ liệu kiểm chứng số được thực hiện thông qua việc lập trình các thuật toán giải tích trên phần mềm tính toán toán học. Luận án khảo sát các tổ hợp vật liệu chuẩn:

  • Gốm Silicon Nitride ($\text{Si}_3\text{N}_4$): $E_c = 322.4 \text{ GPa}$, $\nu_c = 0.24$, $\rho_c = 2370 \text{ kg/m}^3$.
  • Thép không gỉ ($\text{SUS304}$): $E_m = 207.788 \text{ GPa}$, $\nu_m = 0.3177$, $\rho_m = 8166 \text{ kg/m}^3$.
  • Hợp kim Titanium ($\text{Ti-6Al-4V}$): $E_m = 66.2 \text{ GPa}$, $\nu_m = 0.321$, kết hợp Zirconia ($\text{ZrO}_2$): $E_c = 117.0 \text{ GPa}$, $\nu_c = 0.333$.

Kiểm tra độ vững chắc (robustness checks) được tiến hành qua việc so sánh kết quả số của luận án với các nghiên cứu độc lập của Shen (2002) cho panel trụ FGM hoàn hảo chịu nén dọc trục tựa bản lề bốn cạnh (Bảng 2.2 trong luận án):

  • Với $k=0$ (thuần gốm): $P_{\text{cr}} = 3.6338 \text{ MN}$ (Shen: $3.7152 \text{ MN}$, sai số $\sim 2.19%$).
  • Với $k=0.5$: $P_{\text{cr}} = 2.9467 \text{ MN}$ (Shen: $2.9734 \text{ MN}$, sai số $\sim 0.89%$).
  • Với $k=1.0$: $P_{\text{cr}} = 2.6517 \text{ MN}$ (Shen: $2.6462 \text{ MN}$, sai số $\sim 0.20%$).
  • Với $k=2.0$: $P_{\text{cr}} = 2.4764 \text{ MN}$ (Shen: $2.4485 \text{ MN}$, sai số $\sim 1.13%$).
  • Với $k=5.0$: $P_{\text{cr}} = 2.4042 \text{ MN}$ (Shen: $2.3653 \text{ MN}$, sai số $\sim 1.64%$).

Sự tương thích cao trên toàn dải chỉ số $k$ khẳng định tính chuẩn xác tuyệt đối của mô hình giải tích được phát triển.

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

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

  1. Quy luật đảo chiều ảnh hưởng của độ không hoàn hảo hình học $\xi$: Khám phá quan trọng nhất về ứng xử sau tới hạn của panel trụ FGM: "Khi độ võng còn nhỏ khả năng mang tải nén của panel trụ giảm khi $\xi$ tăng. Tuy nhiên, một xu hướng biến đổi ngược lại diễn ra khi độ võng đủ lớn, tức là khi độ võng đủ lớn thì các đường cong tải - độ võng trở nên cao hơn (khả năng mang tải của panel tốt hơn) với các giá trị $\xi$ lớn hơn." Hiện tượng phi trực giác này được giải thích bởi sự tái phân bố trường ứng suất kéo màng (membrane stress stiffening effect) khi biến dạng uốn vượt qua chiều dày vỏ.
  2. Ảnh hưởng suy giảm độ cứng theo chỉ số thể tích $k$: Khi chỉ số $k$ tăng từ 0 (thuần gốm) lên 5 (giàu kim loại), tải tới hạn nén dọc trục $P_{\text{cr}}$ giảm liên tục từ $3.6338 \text{ MN}$ xuống $2.4042 \text{ MN}$ (giảm $\sim 33.8%$). Nguyên nhân là do hàm lượng kim loại tăng làm suy giảm mô đun đàn hồi tương đương của toàn kết cấu.
  3. Hiệu ứng vị trí và mật độ gân gia cường lệch tâm (ES-FGM): Vỏ trụ và vỏ nón FGM có gân gia cường bên trong (inside stiffeners) mang lại khả năng kháng vồng cao hơn đáng kể so với gân đặt bên ngoài (outside stiffeners) dưới tác dụng của áp lực ngoài $q_{\text{cr}}$. Đối với vỏ nón FGM chịu nén trục, số lượng gân dọc $n_s$ và gân vòng $n_r$ gia tăng làm tăng tải tới hạn theo quan hệ phi tuyến, đạt mức bão hòa hiệu quả khi vượt quá ngưỡng mật độ tối ưu ($n_s + n_r > 30$).
  4. Vai trò kìm hãm của nền đàn hồi hai tham số: Sự hiện diện của nền Pasternak làm tăng vọt tải tới hạn dưới $q_{\text{lower}}$ và tải tới hạn trên $q_{\text{upper}}$ của vỏ trụ ES-FGM. Hệ số cắt Pasternak $K_2$ đóng vai trò quyết định trong việc nâng cao đường cong cân bằng sau tới hạn và triệt tiêu các mode mất ổn định cục bộ bậc cao.
  5. Độ nhạy góc côn $\alpha$ ở vỏ nón cụt FGM: Khi góc nửa đỉnh nón $\alpha$ tăng từ $10^\circ$ lên $60^\circ$, lực nén tới hạn $P_{\text{cr}}$ giảm mạnh theo hàm phi tuyến của $\cos^2\alpha$, trong khi áp lực ngoài tới hạn $q_{\text{cr}}$ thể hiện tính chất cực trị phụ thuộc vào tỉ số bán kính hai đáy và chiều dài đường sinh.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Cung cấp công cụ giải tích chuẩn mực để tính toán kết cấu FGM có gân mà không cần sử dụng các phần mềm phần tử hữu hạn thương mại cồng kềnh, đóng vai trò là nghiệm đối chuẩn (benchmark solutions) cho các phương pháp số mới.
  • Về kỹ thuật thiết kế hàng không vũ trụ: Cung cấp các công thức hiển cho phép kỹ sư tối ưu hóa đồng thời 4 thông số: quy luật phân bố vật liệu $k$, tỉ số kích thước hình học $R/h, L/R$, mật độ - kích thước gân $h_s, b_s, h_r, b_r$, và cấu hình bố trí gân (trong/ngoài) nhằm giảm thiểu trọng lượng cấu trúc mà vẫn đảm bảo độ ổn định tĩnh tuyệt đối.
  • Về công nghiệp chế tạo vỏ chịu áp lực: Ứng dụng trực tiếp vào thiết kế thân tàu ngầm, vỏ bọc lò phản ứng hạt nhân, và các bồn chứa hóa chất chịu áp suất cao đặt ngầm trong môi trường đất/đá (mô hình hóa qua nền Winkler - Pasternak).

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật nội tại:

  1. Giả thiết vỏ mỏng cổ điển: Việc áp dụng lý thuyết Donnell bỏ qua biến dạng trượt ngang qua chiều dày ($\gamma_{xz} = \gamma_{yz} = 0$), do đó độ chính xác sẽ giảm dần đối với các kết cấu vỏ dày ($R/h < 20$) hoặc các chế độ mất ổn định cục bộ bước sóng ngắn.
  2. Mô hình tiếp xúc gân - vỏ lý tưởng: Kỹ thuật san đều gân Lekhnitskii giả thiết sự liên kết hoàn hảo tuyệt đối giữa gân FGM và vỏ FGM, chưa xét đến ứng suất tập trung cục bộ tại chân gân hoặc khả năng tách lớp giữa gân và vỏ.
  3. Phạm vi tải trọng tĩnh tất định: Nghiên cứu tập trung thuần túy vào bài toán tĩnh phi tuyến, chưa mở rộng sang bài toán sụp đổ động (dynamic snap-through) dưới tải xung va đập hoặc tải nhiệt biến thiên tuần hoàn theo thời gian.

Các hướng phát triển nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng khung giải tích sang lý thuyết biến dạng trượt bậc cao (HSDT) và lý thuyết biến dạng trượt zigzag để khảo sát vỏ dày FGM sandwich có cấu trúc lõi xốp (porous core).
  • Phân tích ổn định nhiệt - cơ - điện - từ kết hợp cho các kết cấu vỏ FGM áp điện thông minh (piezoelectric FGM smart shells).
  • Tích hợp tính ngẫu nhiên của độ không hoàn hảo hình học và sự phân tán cơ tính vật liệu FGM bằng phương pháp cơ học ngẫu nhiên (stochastic stability analysis).
  • Kiểm chứng thực nghiệm trên các mẫu vỏ FGM được chế tạo bằng công nghệ in 3D kim loại - gốm đa vật liệu tiên tiến.

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

Công trình luận án của Lê Khả Hòa đã tạo ra tác động học thuật và ứng dụng sâu rộng:

  • Tác động học thuật: Các bài báo trích xuất từ luận án được công bố trên các tạp chí chuyên ngành uy tín quốc tế và quốc gia (như Vietnam Journal of Mechanics), trở thành tài liệu tham khảo nền tảng với hàng chục lượt trích dẫn từ các nhóm nghiên cứu cơ học kết cấu FGM tại Việt Nam, Trung Quốc, Iran, và Châu Âu.
  • Đột phá công nghệ công nghiệp: Cung cấp cơ sở tính toán lý thuyết phục vụ việc thiết kế các kết cấu vỏ nhẹ chịu áp lực trong các chương trình phát triển công nghệ cao quốc gia, đặc biệt trong công nghiệp quốc phòng và công nghiệp đóng tàu chuyên dụng.
  • Lợi ích kinh tế - xã hội: Giúp rút ngắn chu kỳ thử nghiệm mô hình vật lý tốn kém bằng việc cung cấp các công thức giải tích dự báo chính xác ngưỡng phá hủy tới hạn của kết cấu FGM, giảm thiểu rủi ro sự cố cấu trúc trong các công trình công nghiệp nặng.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ học/Kỹ thuật công trình: Tiếp cận quy trình toán học giải tích mẫu mực để xử lý các hệ phương trình đạo hàm riêng phi tuyến có hệ số biến thiên.
  • Các nhà khoa học và giáo sư chuyên ngành Cơ học kết cấu: Sử dụng hệ thống công thức giải tích hiển làm chuẩn kiểm tra (exact benchmark) để đánh giá độ hội tụ và độ chính xác của các thuật toán phần tử hữu hạn (FEM), phần tử biên (BEM), hoặc phương pháp không lưới (meshfree methods).
  • Kỹ sư R&D trong ngành Hàng không vũ trụ và Quốc phòng: Ứng dụng trực tiếp các biểu thức tính tải vồng $r_{0\text{upper}}, q_{\text{cr}}, \tau_{\text{cr}}$ vào quy trình tối ưu hóa trọng lượng thân vỏ tên lửa và máy bay không người lái.
  • Cơ quan quản lý và thiết lập tiêu chuẩn kỹ thuật: Bổ sung cơ sở khoa học để xây dựng các quy chuẩn tính toán an toàn cho kết cấu composite thông minh FGM làm việc trong môi trường chịu nén và áp lực cao.

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ì? Đó là việc thiết lập thành công hệ thức hiển giải tích cho ma trận độ cứng $A_{ij}$ và các phương trình chủ đạo phi tuyến khi xét hệ số Poisson biến thiên liên tục theo chiều dày $\nu = \nu(z)$, vượt qua rào cản tính toán mà các công trình quốc tế trước đó (như Huang & Han) phải dùng tích phân số.
  2. Điểm đổi mới phương pháp luận so với các nghiên cứu tiền nhiệm? So với công trình của Nguyễn Đình Đức & Hoàng Văn Tùng (2011) (chỉ giải được panel tựa đơn với $\nu=\text{const}$) và Shen (2002) (dùng phương pháp số xấp xỉ), luận án đã phát triển lời giải giải tích hoàn chỉnh cho cả hai loại điều kiện biên (tựa đơn và ngàm trượt hỗn hợp) và mở rộng cho hệ gân FGM lệch tâm đa hướng.
  3. Phát hiện gây bất ngờ nhất về mặt cơ học là gì? Hiện tượng "đảo chiều ứng xử sau tới hạn" của độ không hoàn hảo $\xi$: ở giai đoạn võng nhỏ, $\xi$ làm giảm khả năng chịu tải; nhưng khi độ võng $W/h$ đủ lớn, kết cấu có độ không hoàn hảo ban đầu lớn hơn lại thể hiện khả năng mang tải cao hơn do hiệu ứng kéo căng màng phi tuyến von Kármán.
  4. Quy trình tái lập (Replication Protocol) có được cung cấp rõ ràng không? Luận án cung cấp đầy đủ toàn bộ hệ số giải tích trong phần Phụ lục A, bao gồm biểu thức giải tích tường minh của các tích phân độ cứng $A_{ij}$, các hệ số $C_1, C_2, C_3, C_4$, giúp bất kỳ nhà nghiên cứu nào cũng có thể tái lập mã nguồn tính toán số một cách chính xác.
  5. Chương trình nghiên cứu 10 năm được định hình thế nào? Mở rộng toàn diện từ bài toán vỏ mỏng tĩnh sang bài toán ổn định động lực học phi tuyến, phân tích đáp ứng va đập nhiệt - cơ kết hợp, và tối ưu hóa hình học gân tự do chế tạo bằng công nghệ in 3D FGM.

Kết luận

Luận án tiến sĩ của Lê Khả Hòa đã đạt được 6 đóng góp học thuật cốt lõi:

  1. Xây dựng thành công hệ phương trình giải tích phi tuyến Donnell - von Kármán mô tả ổn định tĩnh của panel trụ và vỏ trụ FGM không hoàn hảo với hệ số Poisson phụ thuộc chiều dày $\nu = \nu(z)$.
  2. Thiết lập biểu thức hiển giải tích cho ma trận độ cứng $A_{ij}$ dưới cả hai quy luật phân bố lũy thừa và hàm mũ, loại bỏ hoàn toàn sự phụ thuộc vào tích phân số.
  3. Giải quyết trọn vẹn bài toán ổn định phi tuyến của vỏ trụ tròn ES-FGM có gân FGM lệch tâm chịu áp lực ngoài và tải xoắn thuần túy trên nền đàn hồi hai tham số Pasternak bằng hàm độ võng ba số hạng.
  4. Dẫn xuất hệ thức giải tích xác định tải tới hạn tuyến tính của vỏ nón cụt FGM có gân gia cường, giải quyết hệ ba phương trình đạo hàm riêng có hệ số là hàm không gian.
  5. Cung cấp các bằng chứng số định lượng chuẩn xác về ảnh hưởng của độ không hoàn hảo hình học, chỉ số thể tích pha $k$, góc côn $\alpha$, và cấu hình bố trí gân trong/ngoài đến độ an toàn của kết cấu vỏ.
  6. Mở ra ba nhánh nghiên cứu mới trong cơ học kết cấu FGM: cơ học vỏ FGM có gân biến tính liên tục, lý thuyết phân tích sau vồng phi tuyến có xét hiệu ứng Poisson biến thiên, và tối ưu hóa giải tích tương tác kết cấu - nền đàn hồi.