Tổng quan về luận án

Nghiên cứu về cơ học kết cấu tiên tiến trong thế kỷ 21 chứng kiến sự chuyển dịch mang tính bước ngoặt từ các vật liệu đẳng hướng đồng nhất truyền thống (kim loại, hợp kim) và composite nhiều lớp (laminated composites) sang vật liệu có cơ tính biến đổi – Functionally Graded Materials (FGM). Kể từ khi được đề xuất lần đầu tiên vào năm 1984 bởi nhóm các nhà khoa học tại Viện Nghiên cứu Sendai (Nhật Bản), FGM đã chứng minh tính ưu việt vượt trội trong việc loại trừ hiện tượng tập trung ứng suất nhiệt tại bề mặt tiếp xúc và ngăn ngừa sự phân tách lớp (delamination) nhờ sự biến thiên thành phần trơn và liên tục trong không gian. Tuy nhiên, việc ứng dụng FGM vào các kết cấu dạng vỏ phức tạp chịu môi trường khắc nghiệt đòi hỏi những công cụ giải tích và mô hình toán học chính xác cao để dự đoán các hiện tượng mất ổn định phi tuyến. Luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 62520101) của nghiên cứu sinh Vũ Thị Thùy Anh với đề tài "Phân tích ổn định phi tuyến của vỏ cầu làm bằng vật liệu composite FGM", dưới sự hướng dẫn khoa học của GS.TSKH Nguyễn Đình Đức tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội (2017), là công trình tiên phong giải quyết trọn vẹn bài toán ổn định tĩnh phi tuyến, xác định tải vồng tới hạn và ứng xử sau mất ổn định của các dạng kết cấu vỏ cầu FGM dưới tác động đồng thời của trường lực cơ học và nhiệt độ cao.

Khoảng trống học thuật (research gap) mà luận án định vị bắt nguồn từ thực trạng: phần lớn các nghiên cứu trước đây chủ yếu tập trung vào các kết cấu FGM hình học cơ bản như dầm, tấm phẳng hoặc vỏ trụ tròn, vỏ nón cụt; trong khi các công trình về vỏ cầu thoải (shallow spherical shells), đặc biệt là các dạng hình học đặc biệt như vỏ cầu nhẫn (spherical toroid shells) và mảnh cầu nhẫn (spherical panel shells) làm từ FGM phân bố theo quy luật hàm lũy thừa (P-FGM) và hàm Sigmoid (S-FGM), có xét đến gân gia cường (stiffeners) và đặt trên nền đàn hồi hai tham số (Winkler – Pasternak) chịu tải cơ nhiệt kết hợp, hầu như còn rất khan hiếm hoặc chưa có lời giải giải tích tường minh.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật biến thiên thành phần vật liệu (P-FGM, S-FGM) và tính chất vật liệu phụ thuộc nhiệt độ tác động như thế nào đến độ võng phi tuyến, tải tới hạn vồng trên ($q_u$) và vồng dưới ($q_l$) của kết cấu vỏ cầu?
  2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế tương tác giữa nền đàn hồi hai tham số Winkler – Pasternak ($k_1, k_2$), tính không hoàn hảo hình học ban đầu ($\mu$), và các điều kiện biên (ngàm tựa tự do - FM so với ngàm tựa cố định - IM) ảnh hưởng như thế nào đến hiện tượng hóp (snap-through) của vỏ cầu thoải?
  3. Câu hỏi nghiên cứu 3 (RQ3): Mô hình hóa gân gia cường trực giao theo lý thuyết san đều tác dụng gân của Lekhnitsky nâng cao khả năng mang tải sau vồng của vỏ cầu nhẫn và mảnh cầu nhẫn FGM trong môi trường nhiệt độ ra sao?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết H1: Việc tích hợp cấu trúc phân bố Sigmoid (S-FGM) và giảm chỉ số tỷ lệ thể tích $k$ sẽ làm tăng độ cứng chống uốn hiệu dụng, từ đó nâng cao tải trọng tới hạn và giảm thiểu độ khắc nghiệt của ứng xử mất ổn định kiểu cực trị.
  • Giả thuyết H2: Tham số tương tác cắt của nền Pasternak ($k_2$) có đóng góp áp đảo so với tham số nén Winkler ($k_1$) trong việc kiềm chế biến dạng phi tuyến và nâng cao ổn định sau vồng.
  • Giả thuyết H3: Kỹ thuật san đều gân Lekhnitsky kết hợp với phương pháp hàm ứng suất và Bubnov – Galerkin cho phép thiết lập các biểu thức giải tích đóng (closed-form solutions) chính xác cao để dự đoán trạng thái cân bằng sau mất ổn định của vỏ cầu FGM.

Khung lý thuyết tổng thể của luận án tích hợp chặt chẽ: Lý thuyết vỏ thoải phi tuyến Von Kármán – Donnell, Lý thuyết biến dạng trượt bậc nhất của vỏ (FSDT), Lý thuyết san đều tác dụng gân Lekhnitsky, Mô hình nền đàn hồi Winkler – Pasternak và Phương trình truyền nhiệt Fourier 1D trạng thái dừng qua chiều dày vỏ cầu. Phạm vi nghiên cứu bao quát các kết cấu vỏ cầu thoải đối xứng trục, vỏ cầu nhẫn và mảnh cầu nhẫn với dải thông số hình học rộng ($R/h = 70 \div 90$, $r_0/R = 0.94$, góc mở $\theta$), hệ vật liệu gốm – kim loại điển hình ($Si_3N_4/SUS304$, $Al/Al_2O_3$) với các hằng số nhiệt cơ phụ thuộc nhiệt độ.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về cơ học kết cấu FGM cho thấy sự phát triển mạnh mẽ qua ba dòng nghiên cứu chính:

  1. Dòng nghiên cứu ổn định vỏ FGM dạng trụ và nón: Các công trình tiêu biểu của Shen (2009, 2012) áp dụng phương pháp giải tích nhiễu loạn tham số bé và đa thức Chebyshev để khảo sát ứng xử sau mất ổn định của vỏ trụ tròn FGM chịu tải cơ - nhiệt; nhóm nghiên cứu của Sofiyev và cộng sự (2007, 2010, 2011) tập trung vào mất ổn định tĩnh và động của vỏ nón cụt mỏng FGM trên nền đàn hồi Pasternak chịu áp lực ngoài và tải nén dọc trục biến thiên; Huang và Han (2008, 2010) khảo sát mất ổn định phi tuyến vỏ trụ FGM bằng phương pháp hàm ứng suất.
  2. Dòng nghiên cứu ổn định vỏ cầu và tấm tròn FGM: Shahsiah và Eslami (2006) áp dụng lý thuyết Donnell – Mushtari – Vlasov để phân tích mất ổn định nhiệt tuyến tính của vỏ cầu thoải FGM; Ganapathi (2007) khảo sát mất ổn định nhiệt độ võng nhỏ; Chao và Lin (2006) cùng với Dumir và cộng sự (1991, 1999) nghiên cứu ứng xử vồng của nắp cầu cạn (spherical caps) và tấm vành khuyên đẳng hướng/trực hướng; Chih-Ping Wu và Yi-Hwa Tsai (2007) sử dụng phương pháp vi phân vuông phân (Differential Quadrature - DQ) cho vỏ cầu nhẫn FGM.
  3. Dòng nghiên cứu cơ học FGM tại Việt Nam: Nhóm nghiên cứu tiên phong của GS.TSKH Nguyễn Đình Đức và cộng sự (từ 2007, 2008, 2014, 2016) đã công bố hệ thống các công trình nền tảng về ổn định tĩnh, động lực học và dao động phi tuyến của tấm, vỏ thoải hai độ cong, vỏ trụ FGM và S-FGM có gân gia cường trong môi trường nhiệt độ; GS. Đào Huy Bích, PGS. Đào Văn Dũng và TS. Hoàng Văn Tùng (2010, 2011, 2013) phát triển các bài toán ổn định phi tuyến tĩnh và động của vỏ cầu thoải FGM chịu tải cơ và áp lực ngoài kết hợp nhiệt độ.
                                KHÔNG GIAN HỌC THUẬT VỀ KẾT CẤU FGM

Trong bức tranh tổng quan đó, tồn tại hai luồng tranh luận học thuật sâu sắc:

  • Tranh luận 1 (Bản chất mất ổn định): Tranh luận giữa trường phái phân tích tuyến tính điểm rẽ nhánh (bifurcation buckling theo quan niệm Euler - Poincaré loại 1) và trường phái phân tích phi tuyến mất ổn định kiểu cực trị (extremum/snap-through buckling theo quan niệm loại 2). Đối với kết cấu dạng vỏ mỏng và vỏ thoải, tiếp cận tuyến tính đánh giá quá cao khả năng chịu lực thực tế do bỏ qua biến dạng lớn và hình thái hóp đột ngột khi đi qua điểm cực trị.
  • Tranh luận 2 (Ảnh hưởng của hệ số Poisson biến thiên): Một số tác giả cho rằng hệ số Poisson $\nu(z)$ thay đổi theo chiều dày ảnh hưởng lớn đến trạng thái ứng suất, tuy nhiên các nghiên cứu thực chứng của Nguyễn Đình Đức và cộng sự (2014, 2016) đã chứng minh sai số khi xem $\nu = \text{const}$ là không đáng kể (dưới 1-2%), tạo cơ sở vững chắc cho việc đơn giản hóa giải tích mà vẫn bảo toàn độ chính xác cơ học.

Vị thế khoa học của luận án được xác lập thông qua việc giải quyết trọn vẹn khoảng trống: thiết lập hệ phương trình phi tuyến hoàn chỉnh và tìm ra lời giải giải tích tường minh cho vỏ cầu thoải P-FGM/S-FGM, vỏ cầu nhẫn và mảnh cầu nhẫn có gân gia cường, đặt trên nền đàn hồi hai tham số, chịu tải cơ - nhiệt kết hợp và có xét đến độ không hoàn hảo hình dáng ban đầu. Khi so sánh với các công trình quốc tế, luận án vượt lên trên nghiên cứu của Chao và Lin (2006) (vốn chỉ xét nắp cầu trực hướng không có biến thiên cơ tính và nền đàn hồi) và tổng quát hóa vượt bậc so với Chih-Ping Wu và Yi-Hwa Tsai (2007) (vốn chỉ áp dụng phương pháp số DQ cục bộ mà không đưa ra được phương trình trạng thái tải - độ võng tường minh).

Đó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à hoàn thiện hệ thống lý thuyết cơ học phi tuyến cho kết cấu vỏ FGM thông qua các đóng góp nền tảng:

  • Mở rộng lý thuyết phi tuyến Von Kármán – Donnell: Luận án đã tích hợp thành công các tham số biến thiên cơ tính P-FGM và S-FGM vào quan hệ biến dạng - chuyển vị phi tuyến dạng cực đối xứng và phi đối xứng cho vỏ cầu, vỏ cầu nhẫn và mảnh cầu nhẫn.

  • Hình thành hệ đại lượng độ cứng tích phân tổng quát: Thiết lập hệ thức giải tích cho các mô đun độ cứng pháp tuyến ($E_1$), mô đun liên kết uốn - kéo ($E_2$), mô đun độ cứng chống uốn ($E_3$) và độ cứng uốn quy đổi $D$:

    $$E_1 = \int_{-h/2}^{h/2} E(z)dz, \quad E_2 = \int_{-h/2}^{h/2} z E(z)dz, \quad E_3 = \int_{-h/2}^{h/2} z^2 E(z)dz, \quad D = \frac{E_1 E_3 - E_2^2}{E_1(1 - \nu^2)}$$

    Đặc biệt, luận án chứng minh rằng đối với vỏ cấu trúc Sigmoid đối xứng (S-FGM ba lớp gốm - kim loại - gốm hoặc ngược lại), thành phần ghép nối uốn - kéo triệt tiêu ($E_2 = 0$), giúp khử bỏ hiện tượng tập trung ứng suất lệch trục đặc trưng của P-FGM thông thường.

  • Xây dựng tiêu chuẩn ổn định kiểu cực trị giải tích: Thiết lập phương trình bậc ba liên hệ giữa áp lực ngoài $q$ và độ võng không thứ nguyên $W/h$. Bằng việc áp dụng điều kiện cực trị $dq/dW = 0 \Leftrightarrow A - 2BW + C W^2 = 0$, luận án đưa ra công thức giải tích xác định chính xác điểm vồng trên ($q_u$) và điểm vồng dưới ($q_l$):

    $$q_{u, l} = q(W_{1,2}), \quad W_{1,2} = \frac{B \mp \sqrt{B^2 - AC}}{C}$$

    Cường độ hóp (mức độ mất ổn định khắc nghiệt) được định lượng tường minh thông qua biểu thức:

    $$\Delta q = q_u - q_l = \frac{4(B^2 - AC)^{3/2}}{3C^2}$$

                                  ĐƯỜNG CONG TẢI - ĐỘ VÕNG & HIỆN TƯỢNG HÓP
    Áp lực ngoài q
        0                                                                 Độ võng W/h

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

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

  1. Lý thuyết cơ học vi mô FGM: Mô hình hóa hàm biến thiên tỷ phần thể tích $V_c(z) = ((2z+h)/2h)^k$ cho P-FGM và hàm Sigmoid hai miền cho S-FGM, kết hợp mô hình nhiệt đàn hồi phụ thuộc nhiệt độ phi tuyến theo phương trình Touloukian:

    $$P_j(T) = P_0(1 + P_{-1}T^{-1} + P_1 T + P_2 T^2 + P_3 T^3)$$

  2. Lý thuyết môi trường tương tác: Mô hình hóa tương tác đàn hồi giữa kết cấu vỏ và môi trường xung quanh thông qua tải trọng tương đương Winkler – Pasternak: $q_e = k_1 w - k_2 \nabla^2 w$.

  3. Lý thuyết kết cấu gia cường trực giao: Áp dụng giải pháp của Đào Huy Bích dựa trên kỹ thuật san đều tác dụng gân Lekhnitsky, quy đổi hệ gân bố trí theo phương kinh tuyến và vĩ tuyến thành độ cứng phân bố đều liên tục mà không làm phức tạp hóa miền tính toán hình học.

Điều kiện biên và phạm vi ứng dụng được xác định rõ: Khảo sát chi tiết hai trạng thái biên thực tế là ngàm trượt tựa tự do (FM: $N_{r0} = 0$) và ngàm cứng tựa cố định (IM: triệt tiêu chuyển vị $u = 0$ tại biên $r = r_0$ theo nghĩa tích phân trung bình).

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng (positivism) 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 sử dụng phương pháp bán giải tích kết hợp giải tích chính xác (analytical and semi-analytical framework), giải quyết hệ phương trình vi phân đạo hàm riêng phi tuyến bậc cao thông qua kỹ thuật hàm ứng suất Airy và phép chiếu trực giao Bubnov – Galerkin. Thiết kế mô hình đa tầng bao gồm:

  • Tầng vi mô/vật liệu: Quy tắc hỗn hợp (Voigt scheme) xác định mô đun đàn hồi $E(z)$, hệ số giãn nở nhiệt $\alpha(z)$, hệ số dẫn nhiệt $K(z)$.
  • Tầng kết cấu/hình học: Hệ tọa độ cầu $(r, \theta, z)$ với phép xấp xỉ vỏ thoải ($r = R \sin\varphi \approx R\varphi$, $\cos\varphi \approx 1$, $Rd\varphi = dr$).

Quy trình nghiên cứu rigorous

Quy trình thực thi toán học được kiểm soát nghiêm ngặt qua 4 bước:

  1. Thiết lập hệ phương trình cơ bản: Xuất phát từ hệ thức Von Kármán – Donnell, định luật Hooke cho vật liệu FGM chịu nhiệt và phương trình cân bằng nội lực, đưa vào hàm ứng suất Airy $F(r, \theta)$ thỏa mãn đồng nhất các phương trình cân bằng trong mặt phẳng:

    $$N_r = \frac{1}{r}\frac{\partial F}{\partial r} + \frac{1}{r^2}\frac{\partial^2 F}{\partial \theta^2}, \quad N_\theta = \frac{\partial^2 F}{\partial r^2}, \quad N_{r\theta}