Luận án TS Cơ học: Phân tích ổn định tĩnh vỏ vật liệu có cơ tính biến thiên của Lê Khả Hòa

Tổng quan về luận án

Nghiên cứu về ứng xử cơ học của các kết cấu tiên tiến làm việc trong điều kiện khắc nghiệt là một trong những trọng tâm hàng đầu của Cơ học vật rắn biến dạng hiện đại. Luận án tiến sĩ của tác giả Lê Khả Hòa với tiê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" (chuyên ngành Cơ học vật rắn, mã số 62440107, thực hiện tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội dưới sự hướng dẫn của PGS. Đào Văn Dũng) giải quyết trọn vẹn bài toán xác định giới hạn chịu tải và ứng xử sau tới hạn của các dạng vỏ mỏng tiên tiến.

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 bởi các nhà khoa học tại Viện Sendai (Nhật Bản), đánh dấu bước đột phá khi kết hợp giữa gốm (ceramic) chịu nhiệt, chống ăn mòn và kim loại (metal) có độ bền dẻo cao. Như văn bản luận án đã trích dẫn: "Vật liệu cơ tính biến thiên thường được tạo nên từ hai loại vật liệu thành phần là gốm (ceramic) và kim loại (metal) với tỷ phần thể tích của mỗi thành phần được lựa chọn một cách hợp lý, biến đổi trơn và liên tục theo bề dầy của kết cấu". Đặc tính biến thiên trơn loại bỏ hoàn toàn sự tập trung ứng suất tại mặt tiếp xúc liên lớp – nhược điểm cố hữu của composite xếp lớp truyền thống.

Khoảng trống học thuật (research gap) cốt lõi mà luận án định vị bao gồm:

1. Hầu hết các nghiên cứu trước đây (như Shen, 2002; Huang & Han, 2008; Nguyễn Đình Đức & Hoàng Văn Tùng, 2011) chỉ xét hệ số Poisson là hằng số ($\nu = \text{const}$) và khảo sát kết cấu vỏ FGM không gân trơn, bỏ qua tính biến thiên của hệ số Poisson theo chiều dày $\nu = \nu(z)$ cũng như chưa thiết lập được hệ thức giải tích tường minh cho ma trận độ cứng $A_{ij}$.
2. Các công trình phân tích phi tuyến vỏ trụ và vỏ nón FGM có gân gia cường lệch tâm (Eccentrically Stiffened FGM - ES-FGM) đặt trên nền đàn hồi hai tham số Winkler-Pasternak dưới tải trọng cơ học phức tạp (nén dọc trục, áp lực ngoài, mômen xoắn thuần túy) bằng phương pháp giải tích với hàm độ võng 3 số hạng còn rất hạn chế.

Để giải quyết khoảng trống này, luận án thiết lập 4 câu hỏi nghiên cứu và giả thuyết khoa học:

• RQ1/H1: Sự phụ thuộc tọa độ của hệ số Poisson $\nu(z)$ theo chiều dày có làm thay đổi đáng kể ma trận độ cứng liên kết dãn-uốn và làm dịch chuyển giá trị tải tới hạn trên/dưới của panel trụ và vỏ trụ không hoàn hảo?
• RQ2/H2: Việc thiết lập các biểu thức giải tích hiển cho hệ số độ cứng $A_{ij}$ cho phép triệt tiêu sai số tích phân số và cung cấp nghiệm giải tích đóng chính xác cho phương trình vi phân phi tuyến Donnell-von Kármán?
• RQ3/H3: Gân gia cường FGM lệch tâm (ES-FGM) kết hợp kỹ thuật san đều gân Lekhnitsky làm tăng cường khả năng chịu tải tới hạn xoắn, nén và áp lực ngoài vượt trội so với gân thuần nhất?
• RQ4/H4: Bằng phương pháp biến phân Galerkin, bài toán ổn định tuyến tính của vỏ nón cụt FGM có gân với hệ phương trình vi phân đạo hàm riêng hệ số biến thiên có thể tìm được nghiệm hiển giải tích chính xác?

Khung lý thuyết của luận án tích hợp lý thuyết vỏ mỏng phi tuyến hình học Donnell-von Kármán, kỹ thuật san đều tác dụng gân Lekhnitsky, quy tắc hỗn hợp lũy thừa (P-FGM) và hàm mũ (E-FGM), cùng mô hình tương tác kết cấu - nền đàn hồi hai tham số Winkler - Pasternak. Phạm vi nghiên cứu bao quát các dạng hình học: panel trụ, vỏ trụ tròn mỏng và vỏ nón cụt FGM chịu tải nén trục ($p$), áp lực đều mặt ngoài ($q$), và xoắn thuần túy ($\tau$), với kích thước hình học khảo sát rộng ($R/h = 100 \div 200$, $b/h = 30 \div 75$, $a/b = 0.5 \div 1.5$, góc côn $\alpha = 15^\circ \div 60^\circ$).

Literature Review và Positioning

Tổng quan y văn trong luận án tổng hợp ba dòng nghiên cứu chính trong cơ học kết cấu FGM:

Dòng nghiên cứu 1: Ổn định phi tuyến của tấm và vỏ FGM không gân Khởi đầu từ các nghiên cứu kinh điển về vỏ đẳng hướng của Brush & Almroth (1975) và Seide (1956, 1962), lý thuyết ổn định vỏ phát triển mạnh cho vật liệu dị tính. Shen (2002, 2009) cùng Shen & Noda (2005) sử dụng phương pháp tham số bé và lý thuyết biến dạng trượt bậc cao (HSDT) để khảo sát sau tới hạn của vỏ FGM nhưng phụ thuộc vào giải pháp số. Huang & Han (2008, 2009, 2010) áp dụng lý thuyết biến dạng lớn và phương pháp Ritz để phân tích vỏ trụ FGM chịu nén, áp lực và xoắn, tuy nhiên các hệ số độ cứng $A_{ij}$ vẫn duy trì ở dạng tích phân số, chưa giải quyết được biểu thức giải tích hiển. Nguyễn Đình Đức & Hoàng Văn Tùng (2011) nghiên cứu giải tích panel trụ FGM tựa đơn nhưng đóng khung trong giả thiết hệ số Poisson không đổi ($\nu = \text{const}$).

Dòng nghiên cứu 2: Kết cấu FGM có gân gia cường và nền đàn hồi Najafizadeh và cộng sự (2009) khảo sát ổn định tuyến tính vỏ trụ gia cường gân FGM theo chuyển vị. Đào Huy Bích cùng các cộng sự (2011, 2013) tiên phong tiếp cận giải tích phi tuyến cho tấm và vỏ thoải FGM có gân gia cường lệch tâm sử dụng lý thuyết vỏ cổ điển kết hợp kỹ thuật san gân. Sofiyev (2007, 2010) và Bagherizadeh et al. (2011) phát triển các mô hình vỏ FGM trên nền Pasternak. Đào Văn Dũng và Vũ Hoài Nam (2013, 2014) khảo sát ổn định động của vỏ trụ gia cường. Tuy nhiên, việc áp dụng hàm độ võng ba số hạng để bắt trọn trạng thái sau vồng phi tuyến của vỏ trụ ES-FGM chịu xoắn và áp lực ngoài đồng thời vẫn chưa được hoàn thiện tường minh.

Dòng nghiên cứu 3: Ổn định và dao động vỏ nón FGM Singer (1966), Tong & Wang (1992), Sofiyev (2005, 2009, 2011) và Malekzadeh & Heydarpour (2012) phân tích vỏ nón FGM không gân. Đối với vỏ nón có gân, Mecitoglu (1996) và Rao & Reddy (2008) phân tích dao động và thiết kế tối ưu trên vỏ đẳng hướng. Luận án chỉ ra khoảng trống: việc phân tích ổn định giải tích vỏ nón cụt FGM có gân FGM lệch tâm đặt trên nền Pasternak hoàn toàn chưa có lời giải hiển trong y văn quốc tế.

                  ┌─────────────────────────────────────────────────────────┐
                  │                 TIẾN TRÌNH NGHIÊN CỨU                   │
                  └────────────────────────────┬────────────────────────────┘
                                               │
             ┌─────────────────────────────────┴─────────────────────────────────┐
             ▼                                                                   ▼
┌─────────────────────────┐                                         ┌─────────────────────────┐
│     DÒNG TRUYỀN THỐNG   │                                         │    ĐỘT PHÁ CỦA LUẬN ÁN  │
│  (Shen; Huang & Han;    │                                         │      (Lê Khả Hòa)       │
│    Nguyễn Đình Đức)     │                                         │                         │
├─────────────────────────┤                                         ├─────────────────────────┤
│ • ν = const             │ ─── Tranh luận 1: Biến thiên ν(z) ────► │ • ν = ν(z) tường minh   │
│ • Tích phân số A_ij     │ ─── Tranh luận 2: Hệ số giải tích ────► │ • Biểu thức hiển A_ij   │
│ • Vỏ trơn không gân     │ ─── Tranh luận 3: Gân lệch tâm ES ────► │ • ES-FGM + Smeared Rib  │
│ • Võng 1 số hạng        │ ─── Tranh luận 4: Hàm độ võng ────────► │ • Độ võng 3 số hạng     │
│ • Vỏ trụ đơn giản       │ ─── Tranh luận 5: Hình học phức tạp ──► │ • Vỏ nón cụt trên nền   │
└─────────────────────────┘                                         └─────────────────────────┘

Tranh luận học thuật cốt lõi:

• Tranh luận 1 (Hệ số Poisson): Quan điểm truyền thống cho rằng ảnh hưởng của $\nu(z)$ là thứ yếu so với môđun Young $E(z)$. Luận án chứng minh khi kết cấu chịu biến dạng phi tuyến lớn và có độ không hoàn hảo hình học, sự biến thiên của $\nu(z)$ làm thay đổi trường ứng suất nội tại và độ cứng chống uốn, dẫn đến sai lệch đáng kể ở đường cong sau tới hạn.
• Tranh luận 2 (Xấp xỉ hàm độ võng trong tải xoắn): Nghiệm một số hạng (Sofiyev & Schnack, 2004) không mô tả được trạng thái mất ổn định cục bộ và tính bất đối xứng khi chịu xoắn thuần túy; luận án khẳng định bắt buộc phải sử dụng hàm độ võng dạng ba số hạng để xác định chính xác tải vồng cận dưới ($q_{lower}, \tau_{lower}$).

So sánh định lượng với các nghiên cứu quốc tế:

• So sánh tải tới hạn không thứ nguyên $N_x / (ER)$ cho panel trụ thuần nhất chịu nén trục với Turvey (1977) và Shen (2002): tại $a/b = 0.60523$, $b/h = 75$, mode vồng $(m,n) = (2,2)$, kết quả của Turvey đạt $0.3842 \times 10^{-4}$, Shen đạt $0.3842 \times 10^{-4}$, trong khi luận án đạt $0.3841 \times 10^{-4}$ (độ sai lệch tuyệt đối $< 0.03%$).
• So sánh tải nén tới hạn $P_{cr}$ (kN) của panel trụ $Si_3N_4/SUS304$ hoàn hảo bốn cạnh tựa bản lề với Shen (2002): tại $k = 0$, Shen đạt $26.31\text{ kN}$ – Luận án đạt $26.30\text{ kN}$; tại $k = 1$, Shen đạt $20.30\text{ kN}$ – Luận án đạt $20.29\text{ kN}$; tại $k = 5$, Shen đạt $18.15\text{ kN}$ – Luận án đạt $18.14\text{ kN}$.

Đó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 phi tuyến vỏ mỏng Donnell-von Kármán cho môi trường vật chất dị tính liên tục qua chiều dày, mang lại các đột phá lý thuyết:

1. Thiết lập biểu thức hiển giải tích cho ma trận độ cứng $A_{ij}$: Thay vì sử dụng tích phân số phức tạp qua chiều dày, luận án tích phân giải tích trực tiếp quy luật hàm lũy thừa và hàm mũ với hệ số Poisson $\nu(z)$, tạo ra bộ công thức đóng chính xác trong Phụ lục A. Như 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".
2. Khái quát hóa kỹ thuật san đều gân Lekhnitsky cho gân FGM lệch tâm (ES-FGM): Tích hợp hiệu ứng lệch tâm của hệ thống gân dọc (stiffeners) và gân vành (rings) làm bằng FGM, phản ánh chính xác sự dịch chuyển trục trung hòa và tương tác tương hỗ giữa gân và vỏ.
3. Phân biệt rạch ròi hai cơ chế mất ổn định: Phân tích giải tích đồng thời mất ổn định rẽ nhánh (Loại 1 - Bifurcation buckling) cho bài toán tuyến tính và mất ổn định cực trị (Loại 2 - Limit point / Snap-through buckling) cho bài toán phi tuyến thông qua việc tìm các điểm cực trị $dp/dw = 0$, xác định tải vồng cận trên ($q_{upper}$) và cận dưới ($q_{lower}$).

                             ┌──────────────────────────────────────┐
                             │       KHUNG PHÂN TÍCH ĐỘC ĐÁO        │
                             └──────────────────┬───────────────────┘
                                                │
         ┌──────────────────────────────────────┼──────────────────────────────────────┐
         ▼                                      ▼                                      ▼
┌──────────────────────────────┐ ┌──────────────────────────────┐ ┌──────────────────────────────┐
│  MÔ HÌNH VẬT LIỆU & HÌNH HỌC │ │   CƠ SỞ CƠ HỌC & ĐỘNG HỌC    │ │     KỸ THUẬT BIẾN PHÂN       │
├──────────────────────────────┤ ├──────────────────────────────┤ ├──────────────────────────────┤
│ • P-FGM / E-FGM Gradation    │ │ • Donnell-von Kármán Shell   │ │ • Hàm ứng suất Airy phi tuyến│
│ • Poisson thay đổi: ν = ν(z) │ │ • Kỹ thuật san gân Lekhnitsky│ │ • Hàm độ võng 3 số hạng      │
│ • Độ lệch hình học: w*(x,y)  │ │ • Nền Winkler-Pasternak      │ │ • Trực giao Galerkin         │
└──────────────────────────────┘ └──────────────────────────────┘ └──────────────────────────────┘
                                                │
                                                ▼
                               ┌──────────────────────────────────┐
                               │   NGHIỆM GIẢI TÍCH HIỂN ĐÓNG     │
                               │  (p_cr, q_upper, q_lower, τ_cr)  │
                               └──────────────────────────────────┘

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

Khung phân tích của luận án kết hợp 3 trụ cột lý thuyết:

• Cơ học môi trường liên tục phi tuyến: Thiết lập phương trình tương thích biến dạng và phương trình cân bằng phi tuyến theo hàm ứng suất Airy $\phi(x,y)$ và hàm độ võng $w(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)$$

• Mô hình nền đàn hồi hai tham số: Tích hợp phản lực nền $q_e = K_1 w - K_2 \nabla^2 w$, trong đó $K_1$ là hệ số Winkler mô tả độ lún đàn hồi cục bộ, $K_2$ là hệ số Pasternak phản ánh tương tác cắt của lớp nền.
• Điều kiện biên đa dạng: Giải quyết triệt để hai dạng điều kiện biên phức tạp:
  • Bốn cạnh tựa đơn (bản lề): $w = 0, M_x = 0, N_x = r_0 h$ tại $x = 0, a$; $w = 0, M_y = 0, N_y = -p_0 h$ tại $y = 0, b$.
  • Hai cạnh cong tựa đơn, hai cạnh thẳng ngàm trượt: $w = 0, M_x = 0, N_x = r_0 h$ tại $x = 0, a$; $w = 0, \frac{\partial w}{\partial y} = 0, N_y = 0, N_{xy} = 0$ tại $y = 0, b$.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (Positivism) và chủ nghĩa duy lý diễn dịch (Deductive Rationalism) trong Cơ học giải tích. Phương pháp luận dựa trên suy diễn thuần túy từ phương trình vi phân đạo hàm riêng của cơ học môi trường liên tục đến nghiệm giải tích tường minh thông qua kỹ thuật xấp xỉ biến phân Galerkin, sau đó kiểm chứng chéo (cross-validation) với các bộ dữ liệu số chuẩn quốc tế.

┌──────────────────────────────────────────────────────────────────────────────────┐
│                           QUY TRÌNH PHÂN TÍCH GIẢI TÍCH                          │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ 1. Thiết lập trường biến dạng von Kármán & Định luật Hooke với E(z), ν(z)        │
│    ε_x = ε_x^0 - z k_x;  ε_y = ε_y^0 - z k_y;  γ_xy = γ_xy^0 - 2z k_xy           │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ 2. Tích phân năng lượng xác định nội lực & ma trận độ cứng A_ij giải tích hiển   │
│    (N_x, N_y, N_xy, M_x, M_y, M_xy) = ∫_{-h/2}^{h/2} (σ_x, σ_y, σ_xy)(1, z) dz    │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ 3. Thiết lập hệ phương trình chủ đạo Donnell-von Kármán theo hàm φ và độ võng w  │
│    - Phương trình tương thích biến dạng                                          │
│    - Phương trình cân bằng phi tuyến có xét gân ES-FGM và nền Pasternak          │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ 4. Lựa chọn dạng nghiệm thỏa mãn điều kiện biên & Trực giao hóa Galerkin         │
│    ∫∫ [Phương trình cân bằng] × δw dx dy = 0                                     │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
                                         ▼
┌──────────────────────────────────────────────────────────────────────────────────┐
│ 5. Khử hàm ứng suất φ ➔ Phương trình quan hệ Tải - Độ võng phi tuyến (p - W)     │
│    Xác định cực trị: dp/dW = 0 ➔ Tải tới hạn trên (p_upper) và dưới (p_lower)    │
└──────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

• Xử lý độ không hoàn hảo ban đầu: Hàm độ cong không hoàn hảo được đưa vào hệ thống động học theo dạng:

$$w^*(x,y) = \xi h \sin\left(\frac{m\pi x}{a}\right)\sin\left(\frac{n\pi y}{b}\right)$$

với hệ số $\xi \in [0, 0.5]$, mô tả độ lệch hình học chế tạo ban đầu.

• Hàm độ võng phi tuyến 3 số hạng: Để khảo sát bài toán tải xoắn thuần túy và áp lực ngoài, hàm độ võng được chọn có dạng:

$$w(x,y) = f_0 + f_1 \sin(\alpha x + \beta y) + f_2 \sin(\alpha x - \beta y)$$

bắt trọn dạng vồng xoắn bất đối xứng dọc theo đường sinh vỏ trụ.

• Tính toán ma trận độ cứng gân ES-FGM: Tích hợp các tham số hình học gân (chiều rộng $b_s, b_r$, chiều cao $h_s, h_r$, bước gân $d_s, d_r$) cùng sự biến thiên cơ tính của gân theo chiều dày gân.
• Kiểm chứng độ tin cậy (Triangulation & Validation): Luận án thực hiện 8 so sánh độc lập với các công trình quốc tế tiêu biểu:
  • Brush & Almroth (1975) cho vỏ trụ thuần nhất chịu áp lực ngoài.
  • Seide (1956) và Sofiyev (2005) cho vỏ nón cụt thuần nhất chịu nén trục ($P_{cr}$).
  • Huang & Han (2008, 2010) cho vỏ trụ FGM chịu áp lực ngoài và xoắn thuần túy.
  • Shen (2002) và Turvey (1977) cho panel trụ FGM chịu nén.

Data và phân tích

Hệ thống số liệu tính toán dựa trên các cặp vật liệu thực tế dùng trong công nghiệp hàng không vũ trụ và hạt nhân:

• Hệ vật liệu gốm kim loại Si3N4/SUS304:
  • Kim loại (SUS304): $E_m = 207.788\text{ GPa}$, $\nu_m = 0.3177$, $\rho_m = 8166\text{ kg/m}^3$, $\alpha_m = 15.32 \times 10^{-6}\text{ K}^{-1}$.
  • Gốm ($Si_3N_4$): $E_c = 322.271\text{ GPa}$, $\nu_c = 0.2400$, $\rho_c = 2370\text{ kg/m}^3$, $\alpha_c = 7.47 \times 10^{-6}\text{ K}^{-1}$.
• Hệ vật liệu ZrO2/Ti-6Al-4V:
  • Hợp kim Titan (Ti-6Al-4V): $E_m = 115.8\text{ GPa}$, $\nu_m = 0.321$, $\rho_m = 4429\text{ kg/m}^3$.
  • Gốm Zirconia ($ZrO_2$): $E_c = 168.1\text{ GPa}$, $\nu_c = 0.298$, $\rho_c = 5600\text{ kg/m}^3$.
• Hệ vật liệu Al2O3/Al:
  • Nhôm (Al): $E_m = 70\text{ GPa}$, $\nu_m = 0.33$, $\rho_m = 2707\text{ kg/m}^3$.
  • Gốm Alumina ($Al_2O_3$): $E_c = 380\text{ GPa}$, $\nu_c = 0.26$, $\rho_c = 3800\text{ kg/m}^3$.

Môi trường tính toán: Toàn bộ thuật toán giải tích và quét tham số số học được lập trình tự động hóa trên phần mềm đại số biểu thức MAPLE và MATLAB, đảm bảo độ chính xác số học tuyệt đối.

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

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

Luận án rút ra 5 phát hiện then chốt có bằng chứng định lượng rõ ràng:

1. Ảnh hưởng suy giảm của chỉ số thể tích $k$ đến khả năng chịu tải: Khi chỉ số $k$ tăng từ 0 ($100%$ gốm $Si_3N_4$) lên $k = 1$ và $k = 5$, hàm lượng kim loại trong kết cấu tăng lên làm giảm môđun đàn hồi tương đương, dẫn đến tải nén tới hạn $P_{cr}$ giảm mạnh từ $26.30\text{ kN}$ ($k=0$) xuống $20.29\text{ kN}$ ($k=1$, giảm $22.85%$) và đạt $18.14\text{ kN}$ ($k=5$, giảm $31.03%$).
2. Hiện tượng nghịch lý về độ không hoàn hảo ($\xi$) trong vùng sau tới hạn:
   • Ở vùng độ võng nhỏ ($W/h < 0.5$), độ không hoàn hảo hình học $\xi$ tăng từ $0$ lên $0.5$ làm suy giảm nghiêm trọng tải vồng ban đầu của panel trụ.
   • Ngược lại, khi độ võng đủ lớn ($W/h > 1.0$), các đường cong tải - độ võng sau tới hạn của panel có độ không hoàn hảo lớn ($\xi = 0.5$) lại nằm cao hơn panel hoàn hảo ($\xi = 0$), chứng minh kết cấu có khả năng tái phân phối ứng suất màng và duy trì độ cứng mang tải sau mất ổn định.
3. Hiệu ứng gia cường vượt bậc của gân FGM lệch tâm (ES-FGM):
   • Đối với vỏ trụ chịu áp lực ngoài, việc bổ sung gân ES-FGM làm tăng tải tới hạn $q_{cr}$ từ $150%$ đến hơn $320%$ so với vỏ không gân có cùng khối lượng tương đương.
   • Số lượng gân và cách bố trí (gân đặt trong, gân đặt ngoài hoặc gân kết hợp trong-ngoài) tạo ra sự phân hóa rõ rệt: gân đặt ngoài mang lại hiệu quả kháng xoắn cao hơn gân đặt trong do gia tăng mômen quán tính đối với trục uốn mặt ngoài.
4. Ảnh hưởng quyết định của hệ số nền cắt Pasternak ($K_2$):
   • Sự hiện diện của nền đàn hồi làm tăng tải tới hạn trên $q_{upper}$ và tải tới hạn dưới $q_{lower}$.
   • Nền Pasternak ($K_2 > 0$) có tác động ổn định hóa kết cấu mạnh gấp $2.5 \div 4$ lần so với mô hình nền chỉ có thành phần Winkler ($K_1 > 0, K_2 = 0$), do lớp cắt Pasternak truyền ứng suất tiếp xúc ra toàn bộ diện tích vỏ nón và vỏ trụ.
5. Đặc trưng mất ổn định của vỏ nón cụt FGM: Tải nén tới hạn $P_{cr}$ và áp lực ngoài tới hạn $q_{cr}$ của vỏ nón cụt giảm mạnh khi góc côn $\alpha$ tăng từ $15^\circ$ lên $60^\circ$, đồng thời độ nhạy vồng phụ thuộc phi tuyến vào tỷ số bán kính trên chiều dày $R/h$.

Implications đa chiều

• Về mặt học thuật: Luận án cung cấp bộ nghiệm giải tích đối chuẩn (closed-form benchmark solutions) tin cậy cao, đóng vai trò chuẩn hóa cho việc kiểm tra các phần mềm phần tử hữu hạn (ANSYS, ABAQUS) khi mô phỏng kết cấu FGM phức tạp.
• Về mặt phương pháp luận: Quy trình kết hợp trực giao Galerkin với hàm ứng suất Airy trên hệ tọa độ cong của vỏ nón mở ra hướng tiếp cận chuẩn tắc cho các bài toán vỏ có độ cong thay đổi.
• Về mặt kỹ thuật ứng dụng: Cung cấp công thức tính nhanh trực tiếp cho kỹ sư thiết kế thân tên lửa đẩy, vỏ tàu ngầm lặn sâu và ống trao đổi nhiệt trong lò phản ứng hạt nhân chịu áp lực cao và chênh lệch nhiệt độ lớn.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu:

1. Giới hạn lý thuyết vỏ mỏng cổ điển (CST): Việc áp dụng lý thuyết Donnell-von Kármán bỏ qua biến dạng trượt ngang ($\gamma_{xz}, \gamma_{yz}$), do đó kết quả chỉ đảm bảo độ chính xác cao cho vỏ mỏng ($R/h \ge 50$), kém chính xác hơn đối với vỏ dày trung bình ($R/h < 20$).
2. Kỹ thuật san đều gân (Smeared Stiffener Technique): Kỹ thuật san đều gân giả thiết số lượng gân đủ lớn để xem như độ cứng phân bố đều, do đó chưa mô tả được hiện tượng mất ổn định cục bộ của tấm vỏ giữa hai gân kề nhau (local skin buckling) hoặc sự phá hủy do trượt bóc tách liên kết gân-vỏ.
3. Môi trường tải tĩnh đơn thuần: Chưa mở rộng phân tích bài toán ổn định động dưới tải trọng xung va đập hoặc tải nhiệt phân bố không dừng theo thời gian.

Định hướng phát triển 5-10 năm tới:

• Phát triển mô hình giải tích dựa trên Lý thuyết biến dạng trượt bậc cao (HSDT) và Lý thuyết biến dạng trượt giả tam giác để giải bài toán vỏ FGM dày.
• Mở rộng nghiên cứu vật liệu FGM có chứa lỗ rỗng (porous FGM) hoặc vật liệu composite phân lớp tăng cường ống nano carbon (CNTRC).
• Phân tích ổn định động phi tuyến và đáp ứng chuyển tiếp của vỏ nón FGM có gân dưới tác dụng của sóng xung kích nổ và trường nhiệt phi tuyến kết hợp.

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

• Tác động học thuật: Các kết quả của luận án đã được công bố trên các tạp chí chuyên ngành uy tín (Vietnam Journal of Mechanics - VAST, Tạp chí Khoa học ĐHQGHN), đóng góp nguồn trích dẫn quan trọng cho cộng đồng nghiên cứu Cơ học tính toán và Vật liệu composite trong và ngoài nước.
• Chuyển đổi công nghiệp:
  • Công nghiệp hàng không vũ trụ: Ứng dụng tối ưu hóa chiều dày và mật độ gân nón cụt làm đầu mũi tên lửa, vỏ khoang hồi quyển (re-entry capsules).
  • Công nghiệp quốc phòng & Hàng hải: Thiết kế vỏ chịu áp lực cao cho phương tiện lặn tự hành ngầm (AUV) và tàu ngầm vỏ composite FGM.
  • Công nghiệp năng lượng & Hạt nhân: Ứng dụng chế tạo vách ngăn lò phản ứng nhiệt hạch chịu thông lượng nhiệt siêu cao và áp lực chân không ngoài.
• Tiêu chuẩn thiết kế: Cung cấp cơ sở lý thuyết vững chắc phục vụ việc xây dựng các tiêu chuẩn tính toán độ bền và ổn định kết cấu vật liệu composite chức năng tại Việt Nam.

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

• Nghiên cứu sinh và Học viên cao học chuyên ngành Cơ học: Thừa hưởng khung giải tích mẫu mực, bộ phương trình Donnell-von Kármán cải tiến và mã nguồn thuật toán Galerkin tường minh để mở rộng sang các bài toán dao động và sóng.
• Các nhà khoa học đầu ngành: Tiếp cận bộ nghiệm chuẩn (exact benchmark solutions) để thẩm định các phương pháp số mới (Isogeometric Analysis - IGA, Meshfree methods, Finite Element Methods).
• Kỹ sư R&D hàng không vũ trụ và đóng tàu: Sở hữu công thức đóng cho phép đánh giá tức thời sức chịu tải giới hạn ($p_{cr}, q_{cr}, \tau_{cr}$) mà không cần tiêu tốn hàng tuần chạy mô phỏng phần tử hữu hạn phi tuyến phức tạp.
• Các nhà hoạch định chính sách khoa học & công nghệ: Có thêm luận cứ kỹ thuật để định hướng đầu tư phát triển công nghệ chế tạo vật liệu FGM và in 3D composite tiên tiến phục vụ công nghiệp quốc phòng tự chủ.

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 lý thuyết độc đáo nhất là việc thiết lập hệ thức giải tích đóng tường minh cho ma trận độ cứng $A_{ij}$ khi tích hợp đồng thời sự biến thiên liên tục của hệ số Poisson $\nu = \nu(z)$ và môđun Young $E(z)$ theo quy luật lũy thừa và hàm mũ, từ đó giải quyết trọn vẹn hệ phương trình vi phân phi tuyến Donnell-von Kármán cho panel trụ và vỏ trụ không hoàn hảo. Luận án đã mở rộng lý thuyết vỏ mỏng phi tuyến von Kármán từ vật liệu đẳng hướng đồng nhất sang môi trường FGM dị tính đa chiều có gân gia cường lệch tâm (ES-FGM).

2. Điểm mới về phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế là gì?

So với Huang & Han (2008, 2010) (phải tính $A_{ij}$ bằng tích phân số) và Shen (2002) (dùng phương pháp số xấp xỉ), luận án đã kết hợp kỹ thuật san gân Lekhnitsky mở rộng với phép trực giao Galerkin trên hàm độ võng 3 số hạng:

$$w(x,y) = f_0 + f_1 \sin(\alpha x + \beta y) + f_2 \sin(\alpha x - \beta y)$$

Phương pháp này cho phép rút gọn trực tiếp hệ phương trình đạo hàm riêng phi tuyến bậc cao thành phương trình đại số xác định quan hệ Tải - Độ võng, từ đó tìm được biểu thức giải tích của điểm cực trị $dp/dw = 0$ cho tải vồng cận dưới ($q_{lower}, \tau_{lower}$) mà không gặp phải các khó khăn về hội tụ số.

3. Phát hiện số học nào bất ngờ và mang tính phản trực giác nhất trong luận án?

Phát hiện phản trực giác nhất nằm ở ứng xử sau tới hạn của panel trụ FGM không hoàn hảo: Ở giai đoạn đầu chịu nén ($W/h < 0.5$), độ lệch hình học ban đầu $\xi$ làm suy giảm nghiêm trọng khả năng chịu tải của panel; tuy nhiên, khi biến dạng tiến sâu vào miền sau tới hạn ($W/h > 1.0$), panel có độ lệch lớn ($\xi = 0.5$) lại có đường cong tải - độ võng nằm cao hơn panel hoàn hảo ($\xi = 0$). Điều này chứng minh sự tương tác phi tuyến hình học trong vỏ mỏng FGM kích hoạt hiệu ứng kéo màng (tensile membrane action), giúp kết cấu tự tăng cường khả năng chịu lực sau khi đã bị vồng ban đầu.

Tải nén p
   ▲
   │                   /  Đường cong sau tới hạn của vỏ lệch hình học (ξ = 0.5)
   │                  /  [Hiệu ứng màng kéo tăng tải khi W/h > 1.0]
   │    Điểm cực trị /
   │     (p_upper)  / 
   │       ┌───┐   /  /
   │      /     \ /  /   Đường cong vỏ hoàn hảo (ξ = 0)
   │     /       X  /
   │    /       / \/ 
   │   /       /  /\
   │  /       /  /  \ 
   │ /       /  /    \____ Điểm tải cận dưới (p_lower)
   │/_______/  /
   └────────────────────────────────────────► Độ võng W/h
   0      0.5       1.0       1.5

4. Luận án có cung cấp quy trình tái lập kết quả (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết:

• Toàn bộ các hệ số giải tích của ma trận độ cứng $A_{1j}, A_{2j}, A_{3j}$ trong Phụ lục A.
• Biểu thức giải tích của các hệ số $C_1, C_2, C_3, C_4$ trong phương trình chủ đạo.
• Tọa độ mode vồng $(m, n)$ tương ứng với từng tỷ số hình học ($a/b, R/h, L/R$).
• Các thông số cơ tính vật lý của $Si_3N_4/SUS304$ và $ZrO_2/Ti-6Al-4V$ với đầy đủ bảng hệ số nhiệt độ, giúp bất kỳ nhà nghiên cứu nào cũng có thể lập trình tái lập chính xác $100%$ các bảng số liệu và đồ thị sau tới hạn.

5. Chương trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án được định hình ra sao?

Chương trình 10 năm tới bao gồm:

• Giai đoạn 1 (1-3 năm): Tích hợp Lý thuyết biến dạng trượt bậc cao (HSDT) vào mô hình vỏ ES-FGM nhằm bao quát dải kết cấu vỏ dày trung bình.
• Giai đoạn 2 (4-6 năm): Nghiên cứu ổn định nhiệt - cơ - điện từ ghép nối (Multiphysics Thermo-Electro-Mechanical Buckling) cho vỏ FGM áp điện và vật liệu có độ xốp biến thiên (porous FGM).
• Giai đoạn 3 (7-10 năm): Nghiên cứu tối ưu hóa topology đa mục tiêu cấu trúc gân FGM và chế tạo thực nghiệm kiểm chứng bằng công nghệ in 3D kim loại - gốm đa vật liệu.

Kết luận

Luận án tiến sĩ của tác giả Lê Khả Hòa là một công trình khoa học công phu, chuẩn mực và có giá trị học thuật xuất sắc trong lĩnh vực Cơ học vật rắn biến dạng, với các đóng góp cốt lõi:

1. Xây dựng thành công phương pháp giải tích phân tích ổn định phi tuyến của panel trụ và vỏ trụ tròn FGM không hoàn hảo, tính đến sự biến thiên liên tục của hệ số Poisson $\nu = \nu(z)$ và thiết lập các biểu thức đóng hiển cho ma trận độ cứng $A_{ij}$.
2. Phát triển lời giải giải tích phi tuyến cho vỏ trụ tròn FGM có gân gia cường lệch tâm (ES-FGM) trên nền đàn hồi hai tham số Winkler - Pasternak dưới tác dụng của áp lực ngoài và tải mômen xoắn thuần túy, sử dụng hàm độ võng ba số hạng.
3. Thiết lập hệ phương trình ổn định và tìm được nghiệm giải tích xác định tải tới hạn của vỏ nón cụt FGM có gân gia cường trên nền đàn hồi Pasternak – một bài toán phức tạp chưa từng có lời giải đóng trong y văn trước đó.
4. Chứng minh định lượng các quy luật cơ học sâu sắc về ảnh hưởng của chỉ số gradient $k$, độ lệch hình học $\xi$, tỷ số kích thước hình học ($R/h, L/R$), góc nghiêng côn $\alpha$ và các tham số nền đàn hồi đến khả năng chịu tải tới hạn và ứng xử sau tới hạn.
5. Mở ra 3 hướng nghiên cứu mới đầy tiềm năng: cơ học vỏ FGM cấu trúc xốp, phân tích ổn định đa trường nhiệt - cơ động lực học, và tối ưu hóa hình học gân gia cường cho các kết cấu hàng không vũ trụ hiện đại.
6. Thiết lập hệ thống kết quả đối chuẩn giải tích có độ tin cậy tuyệt đối, phục vụ trực tiếp cho việc thẩm định các thuật toán phần tử hữu hạn thương mại và nâng cao năng lực tự chủ tính toán thiết kế kết cấu tiên tiến trong công nghiệp quốc phòng và dân dụng.