Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu dao động của vỏ composite tròn xoay chứa chất lỏng" giải quyết bài toán động lực học của kết cấu vỏ tròn xoay liên hợp bằng vật liệu composite lớp có tiếp xúc với môi trường chất lỏng.

Tính cấp thiết và khoảng trống nghiên cứu

Các kết cấu vỏ tròn xoay bằng composite lớp cốt sợi/nền nhựa hữu cơ (vỏ trụ bậc, vỏ nón ghép trụ, vỏ nón–trụ–nón) được ứng dụng rộng rãi trong các ngành công nghiệp hàng không vũ trụ, đóng tàu, bồn bể hóa chất và xây dựng. Khi các kết cấu này chứa chất lỏng, đặc tính dao động của hệ thống thay đổi đáng kể so với khi làm việc trong không khí do tương tác thủy động lực học. Đồng thời, tính dị hướng cao của vật liệu composite gây ra các tương tác cơ học màng–uốn–xoắn phức tạp. Mặc dù bài toán dao động của kết cấu kim loại đẳng hướng tiếp xúc với chất lỏng đã được nghiên cứu từ nhiều thập kỷ trước, các công trình khảo sát định lượng dao động của vỏ composite tròn xoay có mặt cắt ngang thay đổi chứa chất lỏng vẫn còn hạn chế. Do đó, việc xây dựng mô hình lý thuyết phù hợp và phương pháp tính toán số có độ tin cậy cao để xác định tần số dao động tự do của các hệ vỏ này là yêu cầu kỹ thuật và khoa học cấp thiết.

Mục tiêu nghiên cứu

Luận án xác định 3 mục tiêu nghiên cứu cụ thể:

  1. Tìm được lời giải số có độ tin cậy bằng phương pháp phần tử liên tục (PTLT) và thực nghiệm về tần số dao động tự do của các kết cấu vỏ composite tròn xoay có mặt cắt ngang thay đổi (vỏ trụ bậc composite, vỏ nón–trụ composite, vỏ nón–trụ–nón composite, vỏ nón–nón–nón composite) chứa các mức chất lỏng khác nhau.
  2. Xác định định lượng ảnh hưởng của mức chất lỏng, khối lượng riêng của vật liệu vỏ và chất lỏng, kích thước hình học kết cấu vỏ, điều kiện biên đến tần số và dạng dao động của kết cấu composite, hỗ trợ tính toán và thiết kế kết cấu vỏ tròn xoay composite lớp chứa chất lỏng.
  3. Thực hiện quy trình thiết kế, chế tạo mẫu và đo đạc thực nghiệm tần số dao động riêng của một số mẫu vỏ tròn xoay dạng nón–trụ, nón–trụ–nón composite chứa và không chứa chất lỏng trong phòng thí nghiệm.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu:
    1. Vỏ trụ bậc composite chứa và không chứa chất lỏng.
    2. Vỏ nón cụt–trụ (gọi tắt là nón–trụ) composite chứa và không chứa chất lỏng.
    3. Vỏ nón cụt–trụ–nón cụt (gọi tắt là nón–trụ–nón) composite chứa và không chứa chất lỏng.
    4. Vỏ nón cụt–nón cụt–nón cụt (gọi tắt là nón–nón–nón) composite chứa và không chứa chất lỏng.
    5. Các mẫu thí nghiệm vỏ nón–trụ, nón–trụ–nón composite sợi thủy tinh/nền polyester không no, chứa nước và không chứa nước.
  • Phạm vi nghiên cứu:
    • Khảo sát dao động tự do của các kết cấu vỏ composite đặt đứng, chứa và không chứa chất lỏng.
    • Kết cấu làm việc trong giới hạn đàn hồi tuyến tính, biến dạng bé, vật liệu composite trực hướng đúng trục.
    • Sử dụng lý thuyết biến dạng cắt bậc nhất của Reissner–Mindlin (FSDT).
    • Môi trường chất lỏng được giả thiết là chất lỏng lý tưởng (không nén được, không nhớt, không có chuyển động xoáy). Vỏ và chất lỏng luôn dao động cùng nhau, áp suất chất lỏng tác dụng theo phương pháp tuyến với bề mặt tiếp xúc. Biểu thức áp suất chất lỏng cho các mức chất lỏng được xây dựng dựa trên công thức giải tích của trường hợp vỏ chứa đầy và phụ thuộc vào chuyển vị pháp tuyến $w$.

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

Các hướng nghiên cứu trước đây

Luận án đã hệ thống hóa các hướng nghiên cứu trong và ngoài nước liên quan đến bài toán dao động của kết cấu vỏ tròn xoay:

  • Dao động của vỏ trụ kim loại và composite: Brenneman và Yang (1992) kết hợp ma trận độ cứng động lực và phần tử hữu hạn (PTHH) cho vỏ kim loại chứa chất lỏng; Han và Liu (1994) áp dụng lý thuyết vỏ mỏng Flügge và ma trận truyền; Amabili (1996, 2002) dùng lý thuyết vỏ cổ điển và phương pháp Ritz; Cho cùng cộng sự (2002), Kim (2002), Shekari (2009, 2010), Amiri cùng cộng sự (2012) dùng mô hình PTHH kết hợp phần tử biên khảo sát vỏ có vách ngăn và nắp bình. Đối với vỏ composite, Lakis và Laveau (1991, 1997), Lakis (2000), Toorani và Lakis (2001) nghiên cứu giải tích vỏ trụ composite dị hướng chịu chất lỏng chảy; Kochupillai, Ganesan và Padmanabhan (2002), Kadoli và Ganesan (2003) xây dựng thuật toán PTHH bán giải tích; Okazaki (2002, 2007) và Yu (2007) tính toán PTHH và thực nghiệm cho vỏ trụ composite chứa chất lỏng.
  • Dao động của vỏ trụ bậc: Warburton và Al-Najafi (1969) dùng PTHH cho vỏ hai bậc; Chang (1979) áp dụng phương pháp Rayleigh–Ritz; Zhou và Yang (1995) dùng phương pháp hàm truyền phân đoạn; Zhang và Xiang (2007) sử dụng kỹ thuật không gian trạng thái với lý thuyết Flügge; Qu (2013) dùng phương pháp phân miền; Ganesan và Sivadas (1990) khảo sát giải tích vỏ composite trực hướng có chiều dày thay đổi.
  • Dao động của vỏ nón và hệ vỏ tròn xoay ghép nối: El Damatty cùng cộng sự (2001, 2002, 2005) nghiên cứu lý thuyết và thực nghiệm vỏ nón, vỏ nón–trụ kim loại chứa lỏng; Senthil Kumar và Ganesan (2008), Ansary cùng cộng sự (2010, 2011) tối ưu hóa và khảo sát vỏ nón tương tác chất lỏng chảy; Lakis và cộng sự (1986, 1992, 2010) kết hợp PTHH và lý thuyết cổ điển cho vỏ nón dị hướng. Với các hệ ghép nối, Galletly và Mistry (1974) dùng sai phân hữu hạn; Irie cùng cộng sự (1984), Liang và Chen (2006) dùng ma trận truyền; Tavakoli và Singh (1989) dùng không gian trạng thái; Cheng và Nicolas (1992), Lee (2002), Qu cùng cộng sự (2013), Ma cùng cộng sự (2014) dùng phương pháp biến phân và chuỗi Fourier–Ritz; Efraim và Eisenberger (2006) xây dựng ma trận độ cứng động lực cho vỏ nón–trụ đẳng hướng; Patel cùng cộng sự (2000, 2014), Shang (2001), Kamat (2001), Kouchakzadeh và Shakouri (2014) phân tích vỏ nón–trụ, nón–nón composite bằng PTHH và lý thuyết Donnell.
  • Nghiên cứu thực nghiệm: Koval’chuk và Lakiza (1995) đo dao động vỏ trụ composite sợi thủy tinh khô; Hosokawa (2000), Wu (2010), Yu (2007) tiến hành đo đạc trên bình chứa; Tạ Thị Hiền (2014) thực nghiệm trên vỏ trụ và vỏ côn composite đơn lẻ chứa nước.
  • Nghiên cứu phương pháp độ cứng động lực và PTLT: Phương pháp được phát triển bởi Clough (1975), Hallauer (1985), Kulla (1985, 2003), Casimir (1996), Banerjee (1989, 2005). Tại Việt Nam, Đỗ Văn Hiến, Nguyễn Xuân Hùng, Nguyễn Tiến Khiêm, Trần Văn Liên nghiên cứu trên dầm; Trần Ích Thịnh, Nguyễn Mạnh Cường (2013, 2016) và Tạ Thị Hiền (2014) thiết lập ma trận độ cứng động lực cho tấm, vỏ trụ và vỏ nón composite chứa chất lỏng.

Khoảng trống nghiên cứu được luận án giải quyết

Từ tổng quan tài liệu, luận án xác định 3 khoảng trống khoa học chưa được công bố trước thời điểm nghiên cứu:

  1. Chưa có công trình nào nghiên cứu dao động tự do của kết cấu vỏ trụ bậc composite chứa chất lỏng bằng phương pháp PTLT.
  2. Chưa có nghiên cứu nào ứng dụng phương pháp PTLT để khảo sát hệ vỏ tròn xoay liên hợp phức tạp gồm vỏ nón–trụ, nón–trụ–nón, nón–nón–nón composite chứa các mức chất lỏng khác nhau.
  3. Chưa có công bố thực nghiệm nào đo đạc tần số dao động tự do của kết cấu vỏ liên hợp nón–trụ và nón–trụ–nón composite chứa nước.
Hướng nghiên cứu Các công trình tiêu biểu trước đây Vị trí và đóng góp của luận án
Vỏ trụ bậc composite chứa lỏng Zhang & Xiang (2007 - vỏ kim loại khô), Qu (2013 - vỏ kim loại), Ganesan & Sivadas (1990 - giải tích vỏ khô) Ứng dụng phương pháp PTLT thiết lập ma trận độ cứng động lực cho vỏ trụ bậc composite có chứa chất lỏng
Vỏ nón–trụ, nón–trụ–nón composite chứa lỏng Patel (2000 - PTHH vỏ khô), Efraim & Eisenberger (2006 - ma trận độ cứng động vỏ kim loại), Sweedan (2009 - PTHH đơn giản hóa) Xây dựng thuật toán PTLT ghép nối đa phần tử nón–trụ, nón–trụ–nón, nón–nón–nón composite chứa nhiều mức chất lỏng
Nghiên cứu thực nghiệm hệ vỏ liên hợp composite Tạ Thị Hiền (2014 - thực nghiệm vỏ trụ, vỏ nón đơn lẻ chứa nước) Thiết kế, chế tạo mẫu và đo thực nghiệm tần số riêng của vỏ nón–trụ (CT-14, CT-27) và nón–trụ–nón (CTC-9) chứa nước

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

Khung lý thuyết cơ học và thủy động lực học

  • Lý thuyết vỏ Reissner–Mindlin: Trường chuyển vị của vỏ xét đến ảnh hưởng của biến dạng cắt bậc nhất được mô tả qua công thức: $$u(x, \theta, z, t) = u_0(x, \theta, t) + z \varphi_x(x, \theta, t)$$ $$v(x, \theta, z, t) = v_0(x, \theta, t) + z \varphi_\theta(x, \theta, t)$$ $$w(x, \theta, z, t) = w_0(x, \theta, t)$$ trong đó $u_0, v_0, w_0$ là chuyển vị tại mặt trung bình; $\varphi_x, \varphi_\theta$ là các góc xoay pháp tuyến; hệ số hiệu chỉnh cắt được chọn là $f = 5/6$.
  • Quan hệ ứng suất – biến dạng và nội lực của vật liệu composite trực hướng: $$\begin{Bmatrix} {N} \ {M} \ {Q} \end{Bmatrix} = \begin{bmatrix} [A] & [B] & 0 \ [B] & [D] & 0 \ 0 & 0 & [F] \end{bmatrix} \begin{Bmatrix} {\varepsilon} \ {\kappa} \ {\gamma} \end{Bmatrix}$$ với $[A]$ là ma trận độ cứng màng, $[B]$ là ma trận tương tác màng–uốn–xoắn, $[D]$ là ma trận độ cứng uốn, và $[F]$ là ma trận độ cứng cắt.
  • Mô hình chất lỏng: Dòng chảy không nhớt, không nén được và không xoáy được mô tả qua hàm thế vận tốc $\Phi(r, \theta, x, t) = \phi(r) f(x) \cos(m\theta) e^{i\omega t}$ thỏa mãn phương trình Laplace $\nabla^2 \Phi = 0$. Nghiệm của $\phi(r)$ được biểu diễn qua các hàm Bessel biến đổi loại 1 và loại 2 ($I_m, K_m$). Áp suất thủy động tác dụng lên mặt tiếp xúc vỏ được xác định qua phương trình Bernoulli: $P = -\rho_f \frac{\partial \Phi}{\partial t} \Big|_{r=R}$, chuyển hóa thành số hạng khối lượng chất lỏng tham gia dao động $m^*$ bổ sung vào phương trình chuyển động của thành vỏ.

Phương pháp Phần tử liên tục (PTLT)

Thuật toán PTLT giải quyết bài toán qua 7 bước:

  1. Thiết lập các phương trình vi phân quan hệ lực – chuyển vị và phương trình chuyển động của phần tử vỏ tròn xoay chứa chất lỏng.
  2. Chọn véc tơ trạng thái 10 thành phần: ${y}m = {u_0, v_0, w_0, \varphi_x, \varphi\theta, N_x, N_{x\theta}, Q_x, M_x, M_{x\theta}}^T$.
  3. Tính đạo hàm của véc tơ trạng thái theo tọa độ đường sinh $x$.
  4. Đưa về hệ phương trình vi phân bậc nhất dạng ma trận: $\frac{d{y}_m}{dx} = [A_m] {y}_m$, trong đó $[A_m]$ là ma trận kích thước $10 \times 10$.
  5. Tích phân hệ phương trình vi phân bằng phương pháp vi phân ma trận truyền để xác định ma trận truyền $[T(x, \omega)] = \exp\left(\int_0^L [A_m] dx\right)$, liên hệ véc tơ trạng thái giữa hai đầu phân đoạn.
  6. Chuyển đổi ma trận truyền thành ma trận độ cứng động lực phần tử cỡ $10 \times 10$ dạng: ${F} = [K(\omega)] {U}$.
  7. Ghép nối ma trận độ cứng động lực các phần tử theo điều kiện liên tục, thu được ma trận độ cứng động lực tổng thể của hệ vỏ kết hợp.

Kỹ thuật xác định tần số riêng: Phương pháp đường cong đáp ứng

Thay vì giải trực tiếp phương trình phi tuyến $\det|[K(\omega)]| = 0$ (dễ gặp khó khăn do chứa các hàm siêu việt) hoặc thuật toán Wittrick–Williams (tính toán phức tạp và thời gian dài đối với kết cấu ghép nối), luận án sử dụng phương pháp đường cong đáp ứng (Response curve method):

  • Đặt một lực kích thích đơn vị dạng điều hòa tại một điểm biên tự do (ví dụ tại điểm $M$, hoành độ $x = L$): $Q_x(L) = \frac{1}{R} \cos(m\theta) e^{i\omega t}$.
  • Lập trình trên Matlab giải hệ phương trình động lực học với dải tần số kích thích $\omega$ biến thiên.
  • Vẽ đồ thị đáp ứng chuyển vị pháp tuyến $w$ theo tần số kích thích $\omega$. Hoành độ của các đỉnh nhọn trên đường cong đáp ứng (nơi xảy ra cộng hưởng) chính là các tần số dao động tự do riêng cần tìm.

Phương pháp nghiên cứu thực nghiệm

  • Chế tạo các mẫu vỏ thực nghiệm bằng vật liệu composite sợi thủy tinh/nền nhựa polyester không no với cấu hình đối xứng $[0^\circ/90^\circ/0^\circ/90^\circ]$.
  • Sử dụng đồ gá cố định điều kiện biên (ngàm một đầu), bố trí các cảm biến gia tốc trên thân vỏ và kích thích xung lực bằng búa gõ chuyên dụng.
  • Dữ liệu gia tốc dao động được thu nhận qua card chuyển đổi tín hiệu, xử lý bằng phần mềm chuyên dụng để chuyển đổi tín hiệu từ miền thời gian sang miền tần số qua phép biến đổi Fourier nhanh (FFT), từ đó xác định phổ tần số dao động riêng.

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

Chương 1: Tổng quan về vấn đề nghiên cứu

Chương 1 tổng hợp tình hình nghiên cứu lý thuyết, phương pháp số và thực nghiệm đối với dao động của vỏ tròn xoay đơn lẻ và ghép nối, có hoặc không tương tác với chất lỏng. Luận án phân tích chi tiết các kết quả nghiên cứu trên thế giới và tại Việt Nam (các nhóm nghiên cứu tại Đại học Bách khoa Hà Nội, Đại học Quốc gia Hà Nội, Đại học Xây dựng, Học viện Kỹ thuật Quân sự). Tác giả chỉ ra tính khả thi và ưu điểm của phương pháp ma trận độ cứng động lực (PTLT) trong việc giảm bậc tự do so với PTHH và xử lý điều kiện liên tục đơn giản hơn phương pháp giải tích. Từ việc xác định các khoảng trống khoa học, chương này thiết lập mục tiêu giải quyết 5 bài toán cụ thể tương ứng với các cấu trúc vỏ liên hợp và thực nghiệm đo đạc.

Chương 2: Nghiên cứu dao động tự do của kết cấu vỏ trụ bậc composite chứa và không chứa chất lỏng

Chương 2 tập trung xây dựng cơ sở toán học và thuật toán PTLT cho kết cấu vỏ trụ bậc composite:

  • Xây dựng ma trận độ cứng động lực: Mô hình hóa vỏ trụ composite trực hướng gồm $N$ lớp chịu áp suất thủy động của chất lỏng. Thiết lập hệ 10 phương trình vi phân bậc nhất cho véc tơ trạng thái ${y}m$, xác định ma trận hệ số $[A_m]$ kích thước $10 \times 10$ phụ thuộc vào các hằng số độ cứng $A{ij}, B_{ij}, D_{ij}, F_{ij}$ và thông số khối lượng quy đổi $m^*$ của chất lỏng.
  • Thiết lập điều kiện liên tục giữa các bậc: Tại các mặt cắt chuyển tiếp giữa hai phân đoạn trụ $i$ và $j$, điều kiện liên tục yêu cầu cân bằng nội lực và tương thích chuyển vị trên 10 thành phần: $$N_{xi} = N_{xj}, \quad Q_{xi} = Q_{xj}, \quad N_{x\theta i} = N_{x\theta j}, \quad M_{xi} = M_{xj}, \quad M_{x\theta i} = M_{x\theta j}$$ $$u_i = u_j, \quad v_i = v_j, \quad w_i = w_j, \quad \varphi_{xi} = \varphi_{xj}, \quad \varphi_{\theta i} = \varphi_{\theta j}$$
  • Ghép nối ma trận tổng thể: Kết cấu vỏ trụ bậc được chia thành các phần tử liên tục tương ứng với từng đoạn chiều dày và trạng thái chứa lỏng/khô. Với vỏ trụ 4 bậc (5 phân đoạn phần tử), ma trận độ cứng động lực ghép nối đạt kích thước $30 \times 30$.
  • Kiểm tra độ tin cậy và khảo sát tham số:
    • So sánh tần số không thứ nguyên $\Omega = \omega R \sqrt{\rho(1-\nu^2)/E}$ của vỏ trụ kim loại hai bậc với các công bố của Warburton, Zhang & Xiang và kết quả PTHH (phần mềm ANSYS).
    • So sánh tần số dao động riêng (Hz) của vỏ trụ 4 bậc composite với các điều kiện biên ngàm–tự do và ngàm–ngàm.
    • Đánh giá ảnh hưởng của mức chất lỏng $H/L$, tỷ lệ hình học $L/R$, số mode sóng vòng $m$, số mode dọc $n$, và hai cấu hình xếp lớp $[0^\circ/90^\circ/0^\circ/90^\circ]_2$ và $[0^\circ/90^\circ/0^\circ/90^\circ]_s$. Kết quả chỉ ra khi mức chất lỏng tăng, tần số dao động riêng giảm mạnh; điều kiện biên ngàm–ngàm cho tần số cao hơn đáng kể so với biên ngàm–tự do.

Chương 3: Nghiên cứu dao động tự do của kết cấu vỏ tròn xoay nón–trụ, nón–trụ–nón và nón–nón–nón composite chứa và không chứa chất lỏng

Chương 3 mở rộng phương pháp PTLT cho các hệ vỏ tròn xoay liên hợp có đường sinh nghiêng:

  • Mô hình phần tử vỏ nón composite: Thiết lập phương trình vi phân trạng thái cho phần tử vỏ nón có bán kính thay đổi theo quy luật $R(x) = R_1 + x \sin\alpha$ ($\alpha$ là nửa góc nón). Ma trận hệ số của vỏ nón chứa các số hạng biến thiên theo tọa độ $x$, được giải bằng kỹ thuật ma trận truyền vi phân chia nhỏ phân đoạn.
  • Ghép nối hệ nón–trụ, nón–trụ–nón và nón–nón–nón: Xây dựng ma trận độ cứng động lực ghép nối dựa trên các phương trình điều kiện biên tại các mặt cắt giao tiếp giữa phần nón và phần trụ, hoặc giữa các phần nón có góc nghiêng khác nhau ($\alpha_1, \alpha_2, \alpha_3$).
  • Kiểm tra độ tin cậy: So sánh tần số dao động nhỏ nhất $\Omega_1 = \omega R_1 \sqrt{\rho h / A_{11}}$ và số mode vòng $m$ của vỏ nón–trụ kim loại và composite trực hướng với công bố của Patel (2000) và Kouchakzadeh (2014) ở các góc nón $\alpha = 0^\circ, 30^\circ, 60^\circ$.
  • Phân tích tham số:
    • Khảo sát ảnh hưởng của tính phân kỳ và hội tụ của vỏ nón: kết cấu nón phân kỳ có tần số dao động tự do khác biệt rõ rệt so với nón hội tụ.
    • Khảo sát ảnh hưởng của góc nón $\alpha_1, \alpha_2, \alpha_3$ lên vỏ nón–nón–nón: góc nón lớn làm giảm độ cứng uốn tổng thể của hệ vỏ, dẫn đến giảm tần số dao động riêng.
    • So sánh ảnh hưởng của hai môi trường chất lỏng có khối lượng riêng khác nhau (nước và rượu) cùng khối lượng riêng của vật liệu vỏ đến phổ tần số dao động. Môi trường chất lỏng có khối lượng riêng lớn hơn tạo hiệu ứng khối lượng gia tăng mạnh hơn, làm tần số giảm sâu hơn.

Chương 4: Nghiên cứu thực nghiệm xác định tần số dao động tự do của kết cấu vỏ tròn xoay nón–trụ, nón–trụ–nón composite chứa chất lỏng

Chương 4 trình bày toàn bộ quá trình thực nghiệm kiểm chứng trong phòng thí nghiệm:

  • Chế tạo mẫu: Chế tạo các mẫu vỏ composite trực hướng cốt sợi thủy tinh/nền polyester không no theo quy trình quấn và đắp lớp thủ công có ép khuôn, cấu hình xếp lớp $[0^\circ/90^\circ/0^\circ/90^\circ]$. Mẫu thí nghiệm bao gồm: vỏ nón–trụ ký hiệu CT-27 (góc nón $27^\circ$), vỏ nón–trụ ký hiệu CT-14 (góc nón $14^\circ$) và vỏ nón–trụ–nón ký hiệu CTC-9.
  • Thiết bị và bố trí thí nghiệm: Mẫu được gắn trên bệ gá thép vững chắc với liên kết ngàm ở đáy. Thiết bị đo gồm búa gõ kích động lực tích hợp cảm biến lực, các đầu đo gia tốc áp điện siêu nhẹ gắn tại các vị trí xác định dọc theo đường sinh vỏ, thiết bị thu nhận tín hiệu đa kênh và máy tính cài đặt phần mềm phân tích phổ dao động.
  • Kết quả thực nghiệm và đối sánh: Luận án đo đạc tín hiệu dao động trong miền tần số ở các trạng thái: vỏ khô (không chứa nước), vỏ chứa nước một phần và vỏ chứa đầy nước. Các đường cong đáp ứng thu được từ thực nghiệm được đặt đối sánh trực tiếp với đường cong đáp ứng tính toán từ thuật toán PTLT trên Matlab.
Mẫu thí nghiệm Trạng thái chứa lỏng Tần số thực nghiệm (Hz) Tần số PTLT lý thuyết (Hz) Đánh giá sai lệch
Vỏ nón–trụ CT-14 Khô (0% nước) Xuất hiện các đỉnh phổ rõ nét Khớp các vị trí đỉnh cộng hưởng Sai số nhỏ giữa đỉnh lý thuyết và đỉnh thực nghiệm
Vỏ nón–trụ CT-14 Chứa đầy nước Tần số dịch chuyển về vùng thấp Đường cong đáp ứng dịch chuyển tương ứng Hiệu ứng khối lượng nước làm giảm tần số trùng khớp quy luật
Vỏ nón–trụ CT-27 Các mức nước khác nhau Tần số giảm dần theo mức nước Tần số giảm đơn điệu theo mức nước Xác nhận tính tương đồng về quy luật suy giảm tần số
Vỏ nón–trụ–nón CTC-9 Chứa nước và khô Đỉnh phổ đo đạc rõ ràng Phù hợp với đường cong đáp ứng PTLT Kiểm chứng độ tin cậy của thuật toán ghép nối đa phần tử

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

Đóng góp mới về mặt lý luận và khoa học

  1. Xây dựng thành công mô hình toán học và thuật toán dựa trên phương pháp phần tử liên tục (PTLT) để giải bài toán dao động tự do của các hệ vỏ tròn xoay composite lớp có mặt cắt ngang thay đổi (vỏ trụ bậc, vỏ nón–trụ, nón–trụ–nón, nón–nón–nón) chứa và không chứa chất lỏng.
  2. Thiết lập hệ phương trình vi phân trạng thái bậc nhất cỡ $10 \times 10$ cho phần tử vỏ nón và trụ composite lớp dựa trên lý thuyết biến dạng cắt bậc nhất Reissner–Mindlin, kết hợp tương tác áp suất thủy động của chất lỏng không nén được thông qua hàm thế vận tốc và hàm Bessel biến đổi.
  3. Ứng dụng thành công kỹ thuật đường cong đáp ứng kết hợp với thuật toán tính ma trận truyền vi phân mũ ma trận để xác định chính xác tần số dao động riêng của hệ ghép nối đa phần tử, khắc phục hiện tượng tràn số và khó khăn trong việc tìm nghiệm của các hàm siêu việt.

Đóng góp mới về mặt thực tiễn

  1. Cung cấp bộ dữ liệu định lượng có độ tin cậy về tần số dao động tự do và dạng dao động của các kết cấu vỏ tròn xoay composite ghép nối, làm rõ mức độ suy giảm tần số khi thay đổi mực chất lỏng, góc nghiêng của vỏ nón, tỷ lệ kích thước $L/R$ và cấu hình vật liệu.
  2. Hoàn thành quy trình công nghệ chế tạo mẫu và thực nghiệm đo đạc dao động của các hệ vỏ liên hợp nón–trụ và nón–trụ–nón composite chứa nước. Bộ số liệu thực nghiệm đo đạc được đóng vai trò là dữ liệu kiểm chuẩn thực tế cho các phương pháp tính toán số.

Kiến nghị và giải pháp đề xuất

  • Trong tính toán và thiết kế các bồn chứa, tháp hóa chất, silo composite có mặt cắt ngang thay đổi tiếp xúc với chất lỏng, bắt buộc phải tính toán tần số dao động tự do có xét đến khối lượng tham gia của chất lỏng để tránh nguy cơ xảy ra hiện tượng cộng hưởng gây phá hủy kết cấu.
  • Cần chú ý điều chỉnh góc nón $\alpha$ và phân bố chiều dày các bậc vỏ trụ nhằm tối ưu hóa độ cứng uốn và dịch chuyển tần số dao động riêng của công trình ra khỏi dải tần số kích động của động đất, gió hoặc dòng chảy.

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

Hạn chế của luận án

Dựa trên phạm vi nghiên cứu được tác giả xác định trong văn bản, luận án còn một số hạn chế:

  • Các mô hình chỉ khảo sát kết cấu làm việc trong phạm vi đàn hồi tuyến tính và biến dạng bé; chưa xét đến ứng xử phi tuyến hình học hoặc phi tuyến vật liệu.
  • Vật liệu composite được giới hạn ở dạng trực hướng đúng trục; chưa khảo sát các trường hợp góc bện sợi bất kỳ hoặc vật liệu composite cơ tính biến thiên (FGM) có phân bố phức tạp.
  • Mô hình chất lỏng xem là chất lỏng lý tưởng (không nhớt, không nén được, không chuyển động xoáy); áp suất chất lỏng tác dụng vuông góc lên thành vỏ được tính toán theo công thức giải tích của trường hợp chứa đầy và bỏ qua ảnh hưởng của dao động sóng bề mặt tự do (sloshing) cũng như áp suất thủy tĩnh ban đầu.

Hướng nghiên cứu tiếp theo

  • Mở rộng phương pháp PTLT để nghiên cứu dao động của các kết cấu vỏ composite tròn xoay chịu tải trọng động lực học phức tạp hoặc dao động cưỡng bức.
  • Phát triển mô hình tương tác chất lỏng có xét đến tính nhớt của chất lỏng và hiệu ứng sóng mặt tự do.
  • Nghiên cứu dao động của vỏ composite tròn xoay có gân gia cường hoặc làm từ vật liệu composite tiên tiến đa pha.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong lĩnh vực Cơ học vật rắn biến dạng, Cơ học kết cấu và Kỹ thuật công trình:

  • Đối với nghiên cứu sinh và học viên cao học: Luận án cung cấp quy trình chi tiết thiết lập ma trận truyền vi phân, ma trận độ cứng động lực PTLT và phương pháp giải phương trình vi phân bằng hàm mũ ma trận trên Matlab.
  • Đối với giảng viên và nhà khoa học: Các bảng số liệu so sánh tần số dao động tự do không thứ nguyên và tần số Hertz giữa phương pháp PTLT, PTHH và nghiệm giải tích là nguồn dẫn liệu đối chuẩn tin cậy cho các nghiên cứu về kết cấu mỏng.
  • Đối với kỹ sư thiết kế công nghiệp: Quy trình thí nghiệm thực tế tại Chương 4 và các kết luận về ảnh hưởng của mức chất lỏng, góc nón và tỷ lệ hình học cung cấp cơ sở kỹ thuật để thiết kế tối ưu các tháp hấp thụ, bồn bể hóa chất, kết cấu vỏ tàu thuyền và vòm composite chịu tải trọng chất lỏng.

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

1. Phương pháp Phần tử liên tục (PTLT) trong luận án khác phương pháp Phần tử hữu hạn (PTHH) như thế nào?

Trong phương pháp PTHH truyền thống, trường chuyển vị xấp xỉ trong phần tử là trường chuyển vị tĩnh, đòi hỏi phải chia lưới rất mịn đối với kết cấu phức tạp để giảm thiểu sai số. Trong khi đó, phương pháp PTLT (ma trận độ cứng động lực) sử dụng các hàm dạng là nghiệm chuyển vị động chính xác thu được từ hệ phương trình vi phân chuyển động trong miền tần số. Do đó, mỗi phần tử liên tục là một miền liên tục có vô số bậc tự do, giúp giảm đáng kể số lượng phần tử ghép nối và hạn chế việc chia lưới phức tạp mà vẫn đạt được độ chính xác cao.

2. Sự hiện diện của chất lỏng ảnh hưởng như thế nào đến tần số dao động tự do của vỏ composite?

Chất lỏng bên trong vỏ đóng vai trò như một thành phần khối lượng bổ sung (khối lượng thủy động gia tăng $m^*$) nhưng không đóng góp vào độ cứng của kết cấu vỏ. Do đó, khi mức chất lỏng tăng lên ($H/L$ tăng), tần số dao động tự do của kết cấu ở tất cả các mode dao động đều giảm mạnh so với trạng thái vỏ khô. Ngoài ra, chất lỏng có khối lượng riêng càng lớn (ví dụ nước so với rượu) thì mức độ suy giảm tần số dao động riêng càng rõ rệt.

3. Kỹ thuật đường cong đáp ứng xác định tần số riêng trong luận án được thực hiện ra sao?

Phương pháp đường cong đáp ứng mô phỏng lại nguyên lý đo dao động thực nghiệm. Tác giả đặt một lực kích thích điều hòa đơn vị vuông góc với mặt trung bình của vỏ tại một điểm biên tự do: $Q_x(L) = \frac{1}{R} \cos(m\theta) e^{i\omega t}$. Sau đó, giải phương trình động lực học tổng thể ${F} = [K(\omega)] {U}$ trên phần mềm Matlab với tần số $\omega$ quét qua dải khảo sát để thu nhận chuyển vị pháp tuyến $w$. Đồ thị quan hệ giữa $w$ và $\omega$ sẽ xuất hiện các đỉnh cộng hưởng đột biến, và hoành độ của các đỉnh này chính là các tần số dao động tự do của kết cấu.

4. Điều kiện liên tục giữa các phân đoạn vỏ trụ bậc hoặc nón–trụ được thiết lập gồm những thành phần nào?

Tại mặt cắt ghép nối giữa hai phân đoạn, luận án thiết lập điều kiện liên tục tương thích hoàn toàn trên 10 thành phần của véc tơ trạng thái, bao gồm 5 thành phần lực/mô men và 5 thành phần chuyển vị/góc xoay: lực màng pháp tuyến ($N_x$), lực cắt ($Q_x$), lực màng tiếp tuyến ($N_{x\theta}$), mô men uốn ($M_x$), mô men xoắn ($M_{x\theta}$), chuyển vị dọc trục ($u$), chuyển vị vòng ($v$), chuyển vị pháp tuyến ($w$), góc xoay quanh trục vòng ($\varphi_x$) và góc xoay quanh trục dọc ($\varphi_\theta$).

5. Các mẫu thí nghiệm trong Chương 4 có thông số cấu tạo và vật liệu như thế nào?

Các mẫu thí nghiệm được chế tạo từ vật liệu composite trực hướng cốt sợi thủy tinh kết hợp nhựa nền polyester không no, xếp lớp đối xứng theo cấu hình $[0^\circ/90^\circ/0^\circ/90^\circ]$. Luận án chế tạo và đo đạc trên 3 mẫu hình học cụ thể gồm: mẫu vỏ nón–trụ CT-27 (nửa góc nón $27^\circ$), mẫu vỏ nón–trụ CT-14 (nửa góc nón $14^\circ$) và mẫu vỏ nón–trụ–nón CTC-9, được khảo sát ở các trạng thái không chứa nước và chứa các mức nước khác nhau.

Kết luận

Luận án đã giải quyết bài toán dao động tự do của các kết cấu vỏ tròn xoay liên hợp composite lớp (vỏ trụ bậc, nón–trụ, nón–trụ–nón, nón–nón–nón) chứa và không chứa chất lỏng thông qua việc xây dựng thuật toán phần tử liên tục (PTLT) và lập trình tính toán trong môi trường Matlab. Kết quả nghiên cứu làm rõ các quy luật suy giảm tần số dưới tác động thủy động của chất lỏng, góc nghiêng nón, điều kiện biên và cấu hình vật liệu. Độ tin cậy của mô hình lý thuyết được kiểm chứng chặt chẽ thông qua việc đối sánh với các công bố quốc tế, mô hình phần tử hữu hạn và hệ số liệu thực nghiệm đo đạc trực tiếp trên các mẫu vỏ composite trong phòng thí nghiệm.