Tổng quan luận án

Kết cấu dạng tấm và vỏ là những thành phần chịu lực cơ bản, xuất hiện phổ biến trong nhiều lĩnh vực kỹ thuật như xây dựng công trình dân dụng (sàn, mái vòm, silo, bể chứa), cơ khí chế tạo (bình áp lực, đường ống, tháp giải nhiệt), công nghiệp đóng tàu (thân vỏ tàu ngầm, vách ngăn) và kỹ thuật hàng không vũ trụ (thân, cánh máy bay, thùng nhiên liệu). Dưới tác dụng của các dạng tải trọng tĩnh và động, trạng thái ứng xử cơ học của tấm và vỏ thể hiện tính chất phức tạp do sự kết hợp giữa biến dạng màng, uốn và trượt cắt ngang. Khi hình học kết cấu trở nên phức tạp hoặc sử dụng các hệ vật liệu mới, các phương pháp giải tích giải nghiệm chính xác gặp nhiều rào cản toán học, đòi hỏi phải áp dụng các phương pháp số, trong đó phương pháp phần tử hữu hạn (Finite Element Method - FEM) đóng vai trò trung tâm.

Tuy nhiên, việc áp dụng FEM truyền thống trong phân tích kết cấu tấm/vỏ vẫn tồn tại các hạn chế kỹ thuật:

  1. Hiện tượng khóa cắt (shear locking) xuất hiện khi chiều dày kết cấu rất nhỏ so với kích thước mặt phẳng;
  2. Hiện tượng khóa màng (membrane locking) phát sinh khi sử dụng lưới phần tử thô hoặc méo mó trong các bài toán uốn chiếm ưu thế;
  3. Hiện tượng dao động giả mạo hay đồng hồ cát (hourglass modes) khi sử dụng các kỹ thuật tích phân giảm thiếu kiểm soát;
  4. Sự gián đoạn của trường biến dạng và ứng suất qua ranh giới giữa các phần tử liền kề;
  5. Chi phí tính toán tăng cao và sự bất ổn định số khi thực hiện phân tích phi tuyến hình học.

Từ những hạn chế trên, luận án tập trung thực hiện các mục tiêu và nhiệm vụ nghiên cứu cụ thể:

  1. Phát triển và thiết lập công thức toán học cho các kỹ thuật phần tử hữu hạn mới nhằm phân tích ứng xử tĩnh, dao động tự do, ổn định và phi tuyến của kết cấu dạng tấm và vỏ;
  2. Xây dựng một tập hợp gồm 4 phần tử tứ giác 4 nút bậc thấp (SQ4H, SQ4T, SQ4C, SQ4P) có cấu trúc đơn giản, giảm thiểu tối đa ảnh hưởng của các hiện tượng khóa cơ học (khóa cắt, khóa màng, khóa thể tích);
  3. Ứng dụng các phần tử đề xuất vào việc phân tích nhiều mô hình kết cấu với các điều kiện biên, hình học và vật liệu phức tạp (vật liệu đẳng hướng, composite nhiều lớp, vật liệu phân lớp chức năng FGM và vật liệu xốp phân lớp chức năng gia cường tiểu cầu graphene FGP-GPLs);
  4. Lập trình tính toán số, kiểm chứng độ hội tụ, độ chính xác và đánh giá sai số so sánh với các nghiệm giải tích và các phần tử hữu hạn tham chiếu.

Đối tượng và phạm vi nghiên cứu của luận án:

  • Đối tượng nghiên cứu: Kết cấu dạng tấm (tấm phẳng hình vuông, hình bình hành, hình thang, tấm gấp các góc $60^\circ$, $90^\circ$, tấm có sườn gia cường) và kết cấu dạng vỏ (vỏ trụ, vỏ cầu, vỏ bán cầu khoét lỗ $18^\circ$, vỏ nửa yên ngựa).
  • Phạm vi nghiên cứu: Phân tích uốn tĩnh tuyến tính, dao động tự do, mất ổn định đàn hồi (chịu nén đơn trục và đa trục) và phân tích uốn phi tuyến hình học trong miền đàn hồi tuyến tính của vật liệu.

Tổng quan tài liệu và vị trí của luận án

Trong lịch sử phát triển phương pháp phần tử hữu hạn cho kết cấu tấm và vỏ, mô hình phần tử đẳng tham số liên tục $C^0$ dựa trên lý thuyết Reissner-Mindlin được Irons và cộng sự cùng Zienkiewicz và cộng sự thiết lập vào những năm 1970-1980 đã cho phép nội suy độc lập trường chuyển vị và góc xoay. Dẫu vậy, mô hình này nhanh chóng bộc lộ nhược điểm khóa cắt khi chiều dày tấm tiến về 0. Việc áp dụng tích phân giảm nhằm triệt tiêu năng lượng biến dạng cắt thừa lại gây ra các mode năng lượng zero giả mạo (hiện tượng đồng hồ cát).

Để khắc phục hiện tượng khóa cắt, hàng loạt kỹ thuật đã được các nhà khoa học quốc tế đề xuất:

  • Nhóm phương pháp nội suy hỗn hợp các thành phần ten-xơ biến dạng tự nhiên MITC (Mixed Interpolation Tensorial Components) được phát triển bởi Bathe và cộng sự với các phần tử tứ giác và tam giác như MITC4, MITC8, MITC9, MITC16;
  • Phương pháp khoảng hở cắt rời rạc DSG (Discrete Shear Gap) do Bletzinger và cộng sự đề xuất, sau đó được phát triển kết hợp trơn biến dạng thành CS-DSG, ES-DSG, NS-DSG bởi Nguyen-Xuan, Nguyen-Thoi và cộng sự;
  • Dòng phần tử tấm Mindlin rời rạc MIN3, MIN4 phát triển bởi Tessler và cộng sự;
  • Phần tử tứ giác 4 nút bậc thấp MISQ20 và MISQ24 do Nguyen-Van và cộng sự cải tiến từ nền tảng MITC4 kết hợp kỹ thuật làm trơn biến dạng uốn và màng trên các miền con.

Đối với hiện tượng khóa màng trong phân tích kết cấu vỏ bằng phần tử phẳng, Lee và cộng sự đã đề xuất việc chia miền tứ giác thành các miền con tam giác để tính toán biến dạng màng và chuyển về các điểm buộc biên, tạo nên phần tử MITC4+. Ở hướng tiếp cận hàm nội suy bậc cao, phần tử phổ tấm PSE (Plate Spectral Element) sử dụng các điểm cầu phương Gauss-Legendre-Lobatto được nghiên cứu bởi Zrahia và cộng sự cùng Sprague và cộng sự. Các mô hình dựa trên đa thức Chebyshev trực giao cũng được khảo sát bởi Liu, He và Dang-Trung.

Song song đó, các phương pháp phần tử hữu hạn trơn SFEM (Smoothed Finite Element Methods) phát triển bởi Liu, Nguyen-Xuan và Nguyen-Thoi (bao gồm làm trơn trên nút NS-FEM, làm trơn trên cạnh ES-FEM, và làm trơn trên miền con CS-FEM) đã chứng minh ưu thế về độ hội tụ và giảm độ cứng giả tạo của mô hình số. Kỹ thuật nội suy kép (Twice Interpolation Strategy - TIS) đề xuất bởi Bui, Wu và Zheng đã cải thiện tính liên tục của trường ứng suất qua biên phần tử trong cơ học chất rắn phẳng. Phương pháp phân tích đẳng hình học IGA (IsoGeometric Analysis) dựa trên hàm B-spline/NURBS của Hughes, Bazilevs, Gómez, Thai-Hoang và Nguyen-Thanh liên kết trực tiếp hình học CAD với phân tích số, dù còn trở ngại về việc xử lý tích phân cục bộ và hàm dạng toàn cục.

Khoảng trống học thuật mà luận án lựa chọn giải quyết là: Tích hợp các kỹ thuật trơn biến dạng trên miền con (CS), chiến lược nội suy kép (TIS), kỹ thuật tổ hợp biến dạng khử khóa kép (màng và cắt), và chuỗi đa thức trực giao Chebyshev vào việc xây dựng các phần tử tứ giác 4 nút bậc thấp. Hướng đi này duy trì tính đơn giản trong thiết lập ma trận phần tử hữu hạn, tương thích với các chương trình tính toán thương mại hiện hành, đồng thời giải quyết triệt để các dạng khóa số học và cải thiện độ chính xác phân tích ứng suất cho cả kết cấu tấm phẳng, tấm gấp, tấm có sườn và vỏ cong mỏng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án xây dựng trên nền tảng cơ học kết cấu tấm/vỏ và cơ học vật liệu tiên tiến:

  • Lý thuyết biến dạng kết cấu: Luận án áp dụng Lý thuyết biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) và Lý thuyết biến dạng cắt bậc cao (Higher-order Shear Deformation Theory - HSDT). Nhằm khắc phục yêu cầu đạo hàm bậc cao đòi hỏi tính liên tục $C^1$, kỹ thuật cải biên dạng $C^0$-HSDT được sử dụng để đưa trường chuyển vị về tính liên tục $C^0$, phù hợp cho việc thiết lập hàm dạng 4 nút.
  • Mô hình vật liệu:
    1. Vật liệu đẳng hướng: Định nghĩa thông qua mô đun đàn hồi $E$, hệ số Poisson $\mu$ (giá trị thực tế từ 0,0 đến 0,5), mô đun trượt $G = E / [2(1 + \mu)]$, mật độ khối lượng $\rho$ và hệ số giãn nở nhiệt.
    2. Vật liệu composite nhiều lớp: Mô tả quan hệ ứng suất - biến dạng qua ma trận độ cứng $Q_{ij}$ trong hệ trục chính $(1, 2)$ và ma trận chuyển đổi $Q'_{ij}$ khi trục sợi hợp với trục tọa độ quy chiếu $(x, y)$ một góc $\theta$.
    3. Vật liệu phân lớp chức năng (P-FGM): Kết cấu hỗn hợp gốm (ceramic) và kim loại (metal) biến thiên theo chiều dày $z \in [-h/2, h/2]$ theo quy luật hàm lũy thừa: $$V_c(z) = \left(\frac{2z + h}{2h}\right)^n, \quad V_m(z) = 1 - V_c(z)$$ trong đó $n \ge 0$ là chỉ số tỷ phần thể tích.
    4. Vật liệu xốp phân lớp chức năng gia cường tấm tiểu cầu graphene (FGP-GPLs): Kết hợp ba dạng phân bố độ xốp (P-S: đối xứng, P-A: bất đối xứng, P-U: phân bố đều) và ba dạng phân bố graphene (GPL-S, GPL-A, GPL-U). Mối quan hệ giữa độ xốp $e_0 = 1 - E'_2/E'1$, hàm lượng graphene $W{F(GPL)}$, mô đun Young $E(z)$, mô đun cắt $G(z)$ và mật độ $\rho(z)$ được xác định thông qua mô hình vi cơ học Halpin-Tsai và quy tắc hỗn hợp mở rộng.
Nhóm vật liệu Quy luật phân bố cơ tính theo chiều dày Đặc trưng tính toán chính
Đẳng hướng Đồng nhất trên toàn tiết diện $E, G, \mu, \rho$ là hằng số
Composite lớp Từng lớp đồng phương hoặc lệch trục $\theta$ Ma trận biến đổi độ cứng $Q'_{ij}$
P-FGM Hàm lũy thừa thể tích gốm/kim loại qua $z$ Chỉ số gradient $n$
FGP-GPLs Kết hợp hàm phân bố xốp và tỷ phần GPLs Thông số $e_0, W_{F(GPL)}$, quy luật P-S/P-A/P-U
  • Phương pháp nghiên cứu: Luận án sử dụng phương pháp nghiên cứu định lượng hoàn toàn dựa trên mô phỏng số và phân tích cơ học giải tích - số:
    1. Thiết lập phương trình cân bằng biến phân và ma trận độ cứng phần tử $\mathbf{K}$, ma trận khối lượng $\mathbf{M}$, ma trận độ cứng hình học $\mathbf{K}_g$;
    2. Xây dựng các thuật toán phần tử hữu hạn cải tiến trong môi trường tính toán số;
    3. Lập trình giải các bài toán uốn tĩnh $(\mathbf{K}\mathbf{q} = \mathbf{F})$, bài toán dao động tự do $((\mathbf{K} - \omega^2\mathbf{M})\mathbf{q} = \mathbf{0})$, bài toán ổn định đàn hồi $((\mathbf{K} - \lambda\mathbf{K}_g)\mathbf{q} = \mathbf{0})$ và thuật toán lặp phi tuyến hình học (quan hệ phi tuyến độ võng - tải trọng);
    4. Nguồn dữ liệu kiểm chứng: Sử dụng nghiệm giải tích chính xác (Navier, Levy), dữ liệu từ các bài báo khoa học quốc tế uy tín (ISI) và kết quả mô phỏng bằng các phần tử đối chứng chuẩn (MITC4, SQUAD4, RDKQ, HBQ8, KUMBA).

Nội dung chính theo từng chương

Chương 1. Khái quát chung về kết cấu tấm/vỏ và phương pháp phần tử hữu hạn

Chương này tổng kết định nghĩa hình học và cơ học của kết cấu tấm/vỏ (phân loại 1D, 2D, 3D dựa trên độ dày thành mỏng và mặt trung hòa). Tác giả tổng hợp các dạng phần tử cơ bản: phần tử màng, phần tử tấm uốn thuần túy, và phần tử vỏ phẳng (kết hợp trạng thái màng và uốn). Chương 1 phân tích chi tiết cơ chế gây khóa cắt, khóa màng và các phương pháp xử lý trong lịch sử tính toán kết cấu; đồng thời xác định rõ bố cục 8 chương và 4 đóng góp cốt lõi của luận án.

Chương 2. Cơ sở lý thuyết kết cấu tấm/vỏ và công thức phần tử hữu hạn

Chương 2 thiết lập hệ thống phương trình chủ đạo cho các mô hình vật liệu (đẳng hướng, composite trực hướng nhiều lớp, P-FGM, FGP-GPLs). Chương này trình bày chi tiết động học biến dạng theo lý thuyết FSDT và HSDT, nguyên lý thế năng toàn phần dừng, và kỹ thuật chuyển đổi biến dạng phục vụ rời rạc hóa phần tử hữu hạn. Công thức ma trận phần tử vỏ phẳng 4 nút được xây dựng tổng quát cho các bài toán phân tích uốn tĩnh, tần số dao động tự do riêng $\omega$ và tải trọng tới hạn mất ổn định $P_{cr}$.

Chương 3. Phần tử SQ4H

Chương 3 trình bày việc phát triển phần tử tứ giác 4 nút SQ4H dựa trên kỹ thuật làm trơn biến dạng trên miền con (CS) kết hợp lý thuyết biến dạng cắt bậc ba cải biên $C^0$-HSDT.

  • Đặc điểm kỹ thuật: Miền phần tử được chia thành các miền con ($n_c = 1, 2, 4$) để tính toán biến dạng trơn, loại bỏ đòi hỏi đạo hàm bậc cao $C^1$ và khử hiện tượng khóa cắt.
  • Kết quả khảo sát số: Phần tử SQ4H được kiểm chứng qua bài toán phân tích uốn phi tuyến hình học cho tấm vuông, tấm hình bình hành và tấm hình thang ngàm nhiều lớp chịu tải phân bố đều; phân tích uốn phi tuyến tấm gấp góc $60^\circ$ và $90^\circ$ ngàm các biên. Kết quả cho thấy SQ4H cải thiện độ chính xác, giảm thiểu hiện tượng mất ổn định số trong quá trình lặp phi tuyến khi tỷ số cạnh/chiều dày thay đổi.

Chương 4. Phần tử SQ4T

Chương 4 thiết lập phần tử tứ giác 4 nút SQ4T dựa trên chiến lược nội suy kép (Twice Interpolation Strategy - TIS).

  • Đặc điểm kỹ thuật: Hàm dạng bậc cao được xây dựng thông qua việc xét ảnh hưởng của nhóm $n_{sp}$ nút lân cận nằm trong miền ảnh hưởng của phần tử, kết hợp giá trị nút và gradient trung bình nút.
  • Kết quả khảo sát số:
    • Phân tích uốn tĩnh: Kiểm chứng trên kết cấu màng Cook, tấm vuông composite nhiều lớp chịu tải đều và tải hình sin, vỏ trụ chịu lực và tấm hình bình hành FGM. Trường ứng suất $\sigma$ và độ võng chuẩn hóa $w^*$ hội tụ mượt mà, khắc phục hiện tượng gián đoạn ứng suất qua biên phần tử so với phần tử MITC4.
    • Phân tích dao động và ổn định: Đánh giá tần số chuẩn hóa cho tấm vuông $[0^\circ/90^\circ/90^\circ/0^\circ]$ ($a/h = 5, E_1/E_2 = 40$), tính toán tải trọng tới hạn nén đơn trục cho tấm đa lớp và vỏ trụ 5 lớp $[0^\circ/90^\circ/0^\circ/90^\circ/0^\circ]$.
    • Phân tích phi tuyến: Xây dựng đường cong tải trọng - độ võng phi tuyến cho tấm FGM hình bình hành với điều kiện biên ngàm và tựa đơn.

Chương 5. Phần tử SQ4C

Chương 5 phát triển phần tử tứ giác 4 nút SQ4C ứng dụng kỹ thuật tổ hợp biến dạng (Combined strain strategy) tích hợp làm trơn miền con:

  • Đặc điểm kỹ thuật: Áp dụng 4 điểm buộc biên để tính toán biến dạng cắt ngang (khử khóa cắt), đồng thời thiết lập 4 điểm buộc màng kết hợp nút ảo số 5 ở tâm phần tử (khử khóa màng), sau đó áp dụng phép trơn biến dạng trên các miền con tứ giác.
  • Kết quả khảo sát số:
    • Kiểm tra tính uốn tĩnh trên vỏ trụ có màng cứng hai đầu với lưới chia chữ nhật đều và lưới méo, vỏ cầu chịu tải đều, và tấm vuông tựa đơn có gia cường sườn ngang.
    • Phân tích dao động tự do xác định 6 dạng dao động đầu tiên của tấm composite dày ($a/h = 5$) và vỏ cầu ngàm.
    • Phân tích ổn định tấm vuông tựa đơn có $1$ hoặc $n_{st}$ sườn ngang cách đều chịu nén đơn trục; khảo sát sự phụ thuộc của tải trọng tới hạn chuẩn hóa theo tỷ số $a/h$, tỷ số mô đun $E_1/E_2$, và độ cứng sườn.

Chương 6. Phần tử SQ4P

Chương 6 xây dựng phần tử tứ giác 4 nút SQ4P dựa trên chuỗi đa thức trực giao Chebyshev bậc cao ($p_1, p_2$ từ bậc 3 đến bậc 6):

  • Đặc điểm kỹ thuật: Các hàm dạng 1D và 2D xây dựng từ đa thức Chebyshev cho phép tăng bậc nội suy linh hoạt mà không cần tăng thêm nút hình học ngoài 4 nút góc ban đầu.
  • Kết quả khảo sát số:
    • Khảo sát độ võng chuẩn hóa $w^$ và ứng suất $\sigma_x^$ của tấm vuông tựa đơn/ngàm khi tỷ số $a/h$ biến đổi từ 10 đến 100.000 (khẳng định không bị khóa cắt ở dải tấm siêu mỏng).
    • Đánh giá chuyển vị điểm A trên bài toán vỏ bán cầu khoét lỗ $18^\circ$ và kết cấu nửa vỏ yên ngựa.
    • Phân tích toàn diện ứng xử uốn tĩnh, tần số mode 1 và lực tới hạn mất ổn định đơn trục/đa trục của tấm và vỏ làm bằng vật liệu FGP-GPLs với các tổ hợp phân bố độ xốp và tỷ phần graphene (P-S/GPL-S, P-A/GPL-A, P-U/GPL-U) khi thay đổi hệ số xốp $e_0$ và $W_{F(GPL)}$.

Chương 7. Đánh giá sai số chung giữa các phần tử

Chương 7 thực hiện việc so sánh định lượng trực tiếp về độ chính xác, tốc độ hội tụ và thời gian tính toán giữa 4 phần tử đề xuất (SQ4H, SQ4T, SQ4C, SQ4P) cùng các phần tử chuẩn:

  • Bài toán 1: Tấm vuông đẳng hướng chịu tải phân bố đều với biên tựa đơn (SSSS) và ngàm (CCCC), khảo sát sai số độ võng khi thay đổi tỷ số $a/h = 10, 100, 1000, 10000$.
  • Bài toán 2: Phân tích 4 tần số dao động đầu tiên của tấm vuông đẳng hướng và so sánh thời gian thực thi tính toán (giây).
  • Bài toán 3: Độ võng điểm giữa của vỏ cầu đẳng hướng tựa đơn chịu tải phân bố đều.
Phần tử Kỹ thuật cốt lõi Phạm vi ứng dụng tối ưu Ưu điểm nổi bật
SQ4H Trơn miền con + $C^0$-HSDT Tấm phẳng, tấm gấp, phi tuyến hình học Triệt tiêu khóa cắt, tính ổn định số cao trong bài toán phi tuyến
SQ4T Nội suy kép (TIS) qua nút lân cận Tấm/vỏ composite, FGM, phi tuyến Trường ứng suất liên tục qua biên phần tử, giảm sai số cục bộ
SQ4C Điểm buộc màng/cắt + Trơn miền con Tấm/vỏ có sườn gia cường, vỏ cong Khắc phục đồng thời khóa cắt và khóa màng trên lưới chia méo
SQ4P Chuỗi đa thức Chebyshev ($p_1, p_2$) Tấm/vỏ FGP-GPLs, kết cấu siêu mỏng Tăng độ chính xác linh hoạt qua bậc đa thức $p$, không bị khóa khi $a/h$ lớn

Chương 8. Kết luận và hướng phát triển

Chương 8 tổng kết toàn bộ các kết quả đạt được của luận án, khẳng định tính khả thi của 4 phần tử tứ giác 4 nút mới, đồng thời vạch ra các hướng mở rộng nghiên cứu trong tương lai.

Kết quả và những đóng góp mới

Luận án đã đạt được các đóng góp mới về mặt khoa học tính toán và ứng dụng kỹ thuật như sau:

  1. Xây dựng thành công phần tử tứ giác 4 nút SQ4H: Ứng dụng kỹ thuật làm trơn biến dạng trên các miền con kết hợp với lý thuyết biến dạng cắt bậc cao dạng cải biên liên tục $C^0$-HSDT. Phần tử này nâng cao độ chính xác khi mô phỏng uốn phi tuyến của kết cấu tấm phẳng và tấm gấp, giảm thiểu hiện tượng bất ổn định số trong các bước lặp phi tuyến hình học.
  2. Thiết lập phần tử tứ giác 4 nút SQ4T: Vận dụng chiến lược nội suy kép (TIS) để mở rộng hàm nội suy bậc cao dựa trên trường chuyển vị và gradient trung bình của các nút lân cận. Phần tử SQ4T khắc phục sự gián đoạn của trường biến dạng và ứng suất tiếp giáp giữa các biên phần tử trong phân tích tuyến tính và phi tuyến kết cấu tấm/vỏ.
  3. Phát triển phần tử tứ giác 4 nút SQ4C: Đề xuất giải pháp tổ hợp biến dạng toàn diện gồm 4 điểm buộc biến dạng cắt, 4 điểm buộc biến dạng màng (kết hợp nút ảo tâm) và kỹ thuật trơn biến dạng miền con. Phần tử này loại bỏ hiện tượng khóa màng và khóa cắt khi phân tích vỏ cong và kết cấu tấm có sườn gia cường, duy trì độ hội tụ ổn định ngay cả khi lưới phần tử bị méo.
  4. Hình thành phần tử tứ giác 4 nút SQ4P: Kết hợp chuỗi đa thức trực giao Chebyshev vào việc xây dựng trường hàm dạng xấp xỉ cho phần tử tấm/vỏ. Mô hình cho phép điều chỉnh bậc đa thức $p$ linh hoạt, đạt độ chính xác cao trên dải tỷ số độ mảnh cực rộng ($a/h$ từ 10 đến 100.000) và ứng dụng hiệu quả cho vật liệu tiên tiến FGP-GPLs.

Về mặt ứng dụng thực tiễn, các công thức và quy trình thuật toán được thiết lập ở dạng phần tử 4 nút đơn giản, có khả năng tích hợp trực tiếp vào các module phần tử người dùng (UEL/UMAT) trong các phần mềm thương mại hoặc các phần mềm mã nguồn mở để tính toán kết cấu kỹ thuật.

Hạn chế và hướng nghiên cứu tiếp

Dựa trên phạm vi và nội dung được trình bày trong văn bản luận án, các hạn chế và hướng mở rộng nghiên cứu bao gồm:

  • Hạn chế: Các mô hình phần tử đề xuất trong luận án được phát triển dưới dạng phần tử vỏ phẳng 4 nút dựa trên giả thuyết chuyển vị 2D, chưa mở rộng sang phần tử khối không gian 3D; các phân tích phi tuyến mới dừng lại ở phi tuyến hình học trong miền đàn hồi, chưa xét đến tính chất phi tuyến vật liệu (dẻo, từ biến, phá hủy nứt) hoặc bài toán tương tác chất lưu - kết cấu (Fluid-Structure Interaction).
  • Hướng phát triển tiếp theo:
    1. Mở rộng các kỹ thuật phần tử SQ4H, SQ4T, SQ4C, SQ4P cho bài toán phân tích động lực học phi tuyến theo miền thời gian thực;
    2. Nghiên cứu ứng xử phá hủy, lan truyền vết nứt và cơ chế tách lớp trong kết cấu composite và vật liệu phân lớp chức năng;
    3. Áp dụng các phần tử đề xuất vào bài toán tối ưu hóa hình học, tối ưu hóa tỷ phần thể tích vật liệu và sườn gia cường cho kết cấu tấm/vỏ;
    4. Mở rộng tính toán kết cấu tấm/vỏ làm từ các loại vật liệu thông minh (vật liệu áp điện, vật liệu áp từ, vật liệu nhớ hình).

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu đối với các nhóm đối tượng:

  • Nghiên cứu sinh và học viên cao học chuyên ngành Cơ kỹ thuật, Cơ học tính toán, Xây dựng công trình và Cơ khí chế tạo: Cung cấp cơ sở toán học chi tiết về các kỹ thuật trơn biến dạng (SFEM), kỹ thuật nội suy kép (TIS), kỹ thuật khử khóa màng/cắt và ứng dụng đa thức Chebyshev trong phần tử hữu hạn.
  • Giảng viên và nhà nghiên cứu: Cung cấp hệ thống ví dụ số chuẩn hóa, các bảng số liệu so sánh độ võng, tần số dao động và tải trọng tới hạn trên nhiều dạng kết cấu (tấm gấp, vỏ trụ, vỏ cầu, kết cấu FGP-GPLs có sườn gia cường) dùng để đối chuẩn (benchmark) cho các nghiên cứu phát triển phần tử mới.
  • Kỹ sư tính toán kết cấu: Tham khảo phương pháp xây dựng phần tử vỏ phẳng bậc thấp không bị khóa để nâng cao độ chính xác khi phân tích ứng suất cục bộ trong các kết cấu tấm và vỏ mỏng phức tạp.

Câu hỏi thường gặp

1. Luận án phát triển những phần tử tứ giác 4 nút cụ thể nào?
Luận án phát triển 4 phần tử tứ giác 4 nút mới gồm: SQ4H (dựa trên kỹ thuật trơn miền con và lý thuyết biến dạng cắt bậc cao $C^0$-HSDT), SQ4T (dựa trên chiến lược nội suy kép TIS), SQ4C (dựa trên kỹ thuật tổ hợp biến dạng khử khóa màng và khóa cắt), và SQ4P (dựa trên chuỗi đa thức Chebyshev).

2. Điểm đặc thù trong công thức của phần tử SQ4T là gì?
Phần tử SQ4T sử dụng kỹ thuật nội suy kép (TIS), trong đó hàm dạng bậc cao được xây dựng có xét đến ảnh hưởng của nhóm $n_{sp}$ nút lân cận trong miền ảnh hưởng của phần tử, tính toán kết hợp cả giá trị chuyển vị tại nút và gradient trung bình tại nút. Điều này giúp trường biến dạng và ứng suất đạt tính liên tục cao qua ranh giới phần tử.

3. Phần tử SQ4C khử hiện tượng khóa màng và khóa cắt bằng giải pháp nào?
Phần tử SQ4C sử dụng kỹ thuật tổ hợp biến dạng: bố trí 4 điểm buộc biên để tính toán biến dạng cắt ngang (khử khóa cắt), đồng thời chia phần tử thành các miền tam giác con để xác định biến dạng màng qua 4 điểm buộc màng kết hợp nút ảo trung tâm số 5 (khử khóa màng), sau đó áp dụng phép làm trơn biến dạng trên các miền con.

4. Vì sao phần tử SQ4P có thể phân tích chính xác kết cấu từ dày đến siêu mỏng mà không bị khóa cắt?
SQ4P xây dựng hàm dạng dựa trên chuỗi đa thức Chebyshev trực giao với bậc $p_1, p_2$ có thể tùy biến ($p = 3 \div 6$). Tính chất xấp xỉ của đa thức Chebyshev giúp phần tử duy trì độ chính xác và hội tụ ổn định ngay cả khi tỷ số cạnh/chiều dày $a/h$ lên tới 100.000 mà không gặp hiện tượng khóa cắt.

5. Luận án đã khảo sát những mô hình vật liệu tiên tiến nào cho kết cấu tấm/vỏ?
Ngoài vật liệu đẳng hướng và composite nhiều lớp, luận án khảo sát vật liệu phân lớp chức năng P-FGM (gốm - kim loại biến thiên theo hàm lũy thừa) và vật liệu xốp phân lớp chức năng gia cường tiểu cầu graphene FGP-GPLs với 3 dạng phân bố độ xốp (P-S, P-A, P-U) kết hợp 3 dạng phân bố graphene (GPL-S, GPL-A, GPL-U).

Kết luận

Luận án của tác giả Tôn Thất Hoàng Lân đã nghiên cứu và phát triển thành công 4 phần tử phần tử hữu hạn tứ giác 4 nút mới (SQ4H, SQ4T, SQ4C, SQ4P) dùng cho phân tích kết cấu dạng tấm và vỏ. Bằng việc kết hợp các kỹ thuật trơn biến dạng miền con, nội suy kép, tổ hợp biến dạng khử khóa và đa thức Chebyshev, các phần tử đề xuất đã giải quyết hiệu quả các hiện tượng khóa cắt, khóa màng và gián đoạn ứng suất. Toàn bộ các mô hình số được kiểm chứng qua các bài toán uốn tĩnh, dao động tự do, ổn định đàn hồi và uốn phi tuyến trên nhiều hệ vật liệu từ truyền thống đến tiên tiến, đóng góp giải pháp tính toán hữu ích cho lĩnh vực Cơ kỹ thuật và Cơ học công trình.