Tổng quan về luận án
Nghiên cứu động lực học công trình và tương tác kết cấu - chất lỏng (Fluid-Structure Interaction - FSI) đối với vật liệu tiên tiến giữ vị trí cốt lõi trong cơ học hiện đại. Luận án tiến sĩ cơ học của tác giả Tạ Thị Hiền với đề tài "Nghiên cứu dao động của kết cấu tấm và vỏ Composite có tính đến tương tác với chất lỏng" (chuyên ngành Cơ học vật thể rắn, mã số 62440107, Đại học Bách khoa Hà Nội, 2014; người hướng dẫn: GS. Trần Ích Thịnh và TS. Nguyễn Mạnh Cường) là công trình tiên phong tại Việt Nam tích hợp mô hình giải tích - số chính xác cao và thực nghiệm dao động cho kết cấu composite tiếp xúc môi trường thủy tĩnh.
Trong cơ học kết cấu công trình ngầm, tàu thủy và bồn bể công nghiệp, khoảng trống nghiên cứu (research gap) lớn nhất tồn tại ở sự thiếu vắng các thuật toán số hiệu quả mô tả chính xác sự ghép nối giữa trường biến dạng cắt ngang của kết cấu dị hướng nhiều lớp và trường áp suất thủy động của chất lỏng bán vô hạn hoặc có mặt thoáng. Các công trình cổ điển của Lamb (1921), Lindholm (1965), Fu & Price (1987) phần lớn giới hạn ở vật liệu kim loại đẳng hướng. Nghiên cứu của Hashemi et al. (2012) khảo sát tấm Mindlin trong nước nhưng dừng lại ở vật liệu đồng nhất; trong khi các nghiên cứu quốc tế của Xi et al. (1997) hay Matthew et al. (2013) mới tiếp cận bài toán bằng phương pháp phần tử hữu hạn tiêu tốn tài nguyên hoặc giải tích giải tích Navier chỉ giải quyết được các điều kiện biên đơn giản (khớp bốn cạnh).
Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập chặt chẽ:
- Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập một mô hình tính toán số bán giải tích có độ hội tụ chính xác mà không bị phụ thuộc vào mức độ phân chia lưới như phương pháp phần tử hữu hạn truyền thống đối với kết cấu composite ngập nước hoặc chứa chất lỏng?
- Câu hỏi nghiên cứu 2 (RQ2): Các tương tác màng - uốn - xoắn sinh ra từ cấu hình phân lớp dị hướng (đúng trục cross-ply, lệch trục angle-ply) chi phối thế nào đến phổ tần số dao động tự do khi có sự tham gia của khối lượng quy dẫn (added mass) từ chất lỏng?
- Câu hỏi nghiên cứu 3 (RQ3): Sự thay đổi mức chất lỏng $H/L$, góc nghiêng nón $\alpha$, tỷ số bề dày $h/a$, và hệ số nền đàn hồi hai tham số Pasternak/Winkler tác động như thế nào đến độ suy giảm tần số tự do?
- Giả thuyết khoa học (H1): Phương pháp Phần tử liên tục (Continuous Element Method - CEM hay Dynamic Stiffness Method - DSM) dựa trên trường chuyển vị động có thể triệt tiêu hoàn toàn sai số rời rạc hóa không gian của FEM, cung cấp nghiệm chuẩn xác tuyệt đối cho kết cấu tấm và vỏ tròn xoay tương tác chất lỏng.
- Giả thuyết khoa học (H2): Sự hiện diện của chất lỏng làm gia tăng động năng toàn phần của hệ thống mà không làm tăng thế năng biến dạng đàn hồi, dẫn đến hiện tượng suy giảm đơn điệu tần số dao động tự do tỉ lệ thuận với mức ngập $H/L$ và độ mềm uốn của kết cấu.
Khung lý thuyết của luận án tích hợp ba trụ cột: Lý thuyết biến dạng cắt bậc nhất (First-Order Shear Deformation Theory - FSDT) của Reissner-Mindlin, Lý thuyết thế vận tốc chất lỏng lý tưởng (Laplace equation & Bernoulli dynamic pressure), và Lý thuyết phần tử liên tục dựa trên biểu diễn ma trận trạng thái cấp 10. Luận án mang lại đóng góp đột phá khi phát triển thành công bộ 4 chương trình tính toán chuyên dụng viết trên nền tảng MATLAB (VplateF, VplateEF, VcylF, VconF), thiết kế và chế tạo 6 cụm mẫu vỏ trụ tròn và nón cụt composite sợi thủy tinh/polyester không no, thực hiện đo kiểm thực nghiệm dao động bằng kỹ thuật phân tích biến đổi Fourier nhanh (Fast Fourier Transform - FFT) tại Phòng thí nghiệm kiểm soát Rung và Ồn (Viện Cơ học Việt Nam). Kết quả khẳng định độ chính xác vượt trội của phương pháp phần tử liên tục với sai số so với nghiệm giải tích và thực nghiệm dao động nhỏ hơn 5%, đồng thời tiết kiệm thời gian tính toán hàng chục lần so với phần tử hữu hạn (FEM).
Literature Review và Positioning
Nghiên cứu động lực học kết cấu tương tác chất lỏng trải qua gần một thế kỷ tiến hóa với nhiều trường phái tiếp cận. Khởi đầu từ bài toán tấm tròn đàn hồi tiếp xúc với nước của Lamb (1921), tiếp nối bởi các kiểm chứng bằng phương pháp Rayleigh của Powell và Roberts (1923), lý thuyết dao động kết cấu tiếp xúc chất lỏng ban đầu chỉ đóng khung trong phạm vi vật liệu kim loại đẳng hướng. Giai đoạn 1965–1985, Lindholm et al. (1965) tiến hành các thử nghiệm toàn diện trên tấm cantilever ngập nước, Volcy et al. (1979) đối sánh thực nghiệm tấm công-xôn với mô hình PTHH và chỉ ra độ lệch 15% bắt nguồn từ sai số điều kiện biên ngàm. Fu và Price (1987) đã áp dụng PTHH 3D cho bài toán tấm kim loại ngập từng phần trong môi trường nước bán vô hạn, chứng minh dạng dao động trong chân không và trong nước duy trì tính đồng dạng hình học nhưng tần số suy giảm rõ rệt do tải trọng thủy động gia tăng.
Trong phân tích tấm và vỏ composite, hai luồng quan điểm đối nghịch từng tồn tại trong y văn:
- Trường phái Lý thuyết tấm mỏng cổ điển (Classical Laminated Plate Theory - CLPT) dựa trên giả thuyết Kirchhoff-Love bỏ qua biến dạng trượt ngang. Điển hình như các nghiên cứu của Kwak (1991, 1996), Soedel (1994), Amabili et al. (1998) cho rằng biến dạng cắt ngang không đáng kể đối với kết cấu tấm/vỏ thành mỏng ($h/a < 0.01$).
- Trường phái Lý thuyết biến dạng cắt bậc cao và bậc nhất (FSDT/HSDT) khởi xướng bởi Reissner, Mindlin, Reddy (2004), Hui-Shen Shen (2003), Akavci (2007). Nhóm tác giả này chứng minh rằng vật liệu composite nhiều lớp có mô đun trượt ngang $G_{13}, G_{23}$ nhỏ hơn nhiều so với mô đun đàn hồi dọc sợi $E_1$, khiến biến dạng trượt ngang đóng vai trò chi phối và việc áp dụng CLPT sẽ gây đánh giá sai lệch tần số dao động từ 15% đến 35% khi tỷ số bề dày $h/a \ge 0.05$.
TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT DAO ĐỘNG FSI
==========================================================================
GIAI ĐOẠN 1: Kim loại đẳng hướng (1921 - 1990)
[Đặc trưng: Giả thiết tấm mỏng Kirchhoff, giải tích chuỗi, mô hình 2D]
GIAI ĐOẠN 2: Composite & Lý thuyết cắt biến dạng (1990 - 2010)
[Đặc trưng: FSDT/HSDT, mô hình bán giải tích FEM-Navier, tải nhiệt-cơ]
GIAI ĐOẠN 3: Đóng góp của Luận án (Tạ Thị Hiền, 2014)
Tích hợp: FSDT Mindlin + Phương pháp PTLT (DSM) + Thế vận tốc Bernoulli
[Đột phá: Khử hoàn toàn sai số chia lưới FEM, giải đa dạng biên, thực nghiệm]
==========================================================================
Về mặt định vị khoa học (scientific positioning), luận án khai thác trúng điểm nghẽn của các nghiên cứu tiền nhiệm. Các nghiên cứu quốc tế của Hashemi et al. (2012) về tấm ngập nước chỉ khảo sát tấm kim loại đẳng hướng; nghiên cứu của Matthew et al. (2013) và Motley et al. (2013) phân tích tấm composite công-xôn nhưng phụ thuộc hoàn toàn vào phần mềm thương mại ANSYS; nghiên cứu của Korosh Khorshid và Farhadi (2013) sử dụng chuỗi Navier kết hợp phương pháp Ritz bị cô lập trong điều kiện biên bốn cạnh khớp (SSSS). Đối với vỏ tròn xoay, Lakis và Paidoussis (1972, 2004), Selmane và Lakis (1997), Kerboua et al. (2010) phát triển phương pháp phần tử hữu hạn lai (hybrid FEM) cho vỏ kim loại hoặc vỏ composite chứa chất lỏng nhưng gặp hạn chế lớn khi giải bài toán vỏ nón cụt có góc nghiêng thay đổi.
Bằng việc kế thừa và phát triển lý thuyết Phần tử liên tục (PTLT / Dynamic Stiffness Method) từ các công trình nền tảng của Koloušek (1973), Clough và Penzien (1975), Banerjee (1989, 1996), Kulla (2003), Casimir (1996), Nguyen Manh Cuong (2003, 2007), công trình của Tạ Thị Hiền đã lần đầu tiên mở rộng PTLT cho kết cấu tấm composite dị hướng ngập trong chất lỏng, tấm composite đặt trên nền đàn hồi Pasternak không thuần nhất, vỏ trụ tròn và vỏ nón cụt composite chứa chất lỏng với các mức nước biến thiên ($0 \le H/L \le 1$).
Đó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 biến dạng cắt bậc nhất Reissner-Mindlin cho môi trường tương tác nhiều pha: composite lớp dị hướng - chất lỏng lý tưởng - nền đàn hồi hai tham số. Trong các công trình cơ học trước đây, ma trận độ cứng vật liệu composite lớp chỉ xét tương tác màng - uốn đơn thuần trong môi trường không khí. Luận án đã mở rộng mô hình cân bằng động lực học bằng việc tích hợp áp suất thủy động $P(x,y,t)$ của trường thế vận tốc chất lỏng vào hệ phương trình vi phân đạo hàm riêng:
$$\begin{aligned}
\frac{\partial N_{xx}}{\partial x} + \frac{\partial N_{xy}}{\partial y} &= I_0 \frac{\partial^2 u_0}{\partial t^2} + I_1 \frac{\partial^2 \varphi_x}{\partial t^2} \
\frac{\partial N_{xy}}{\partial x} + \frac{\partial N_{yy}}{\partial y} &= I_0 \frac{\partial^2 v_0}{\partial t^2} + I_1 \frac{\partial^2 \varphi_y}{\partial t^2} \
\frac{\partial Q_x}{\partial x} + \frac{\partial Q_y}{\partial y} - P &= I_0 \frac{\partial^2 w_0}{\partial t^2} \
\frac{\partial M_{xx}}{\partial x} + \frac{\partial M_{xy}}{\partial y} - Q_x &= I_1 \frac{\partial^2 u_0}{\partial t^2} + I_2 \frac{\partial^2 \varphi_x}{\partial t^2} \
\frac{\partial M_{xy}}{\partial x} + \frac{\partial M_{yy}}{\partial y} - Q_y &= I_1 \frac{\partial^2 v_0}{\partial t^2} + I_2 \frac{\partial^2 \varphi_y}{\partial t^2}
\end{aligned}$$
Trong đó, các mômen quán tính khối lượng $I_i$ ($i = 0, 1, 2$) được xác định chính xác qua tích phân bề dày:
$$I_i = \sum_{k=1}^N \int_{z_k}^{z_{k+1}} \rho^{(k)} z^i dz$$
Hệ số hiệu chỉnh cắt của Mindlin được áp dụng chặt chẽ dưới dạng $f(z) = \frac{5}{6}\left[1 - \left(\frac{2z}{h}\right)^2\right]$ để thu được ma trận độ cứng cắt trượt ngang hiệu dụng $A_{ij}^*$.
Điểm đột phá lý thuyết thể hiện ở việc lượng hóa chính xác áp suất thủy động của lớp chất lỏng lên hai bề mặt tấm thông qua đại lượng khối lượng gia tán ảo (virtual added mass) $m^* = m_1^* + m_2^*$:
$$P = - m^* \frac{\partial^2 w}{\partial t^2} = \left(\rho_l \frac{1 - c_1 e^{-2\gamma_l h_1}}{\gamma_l (1 + c_1 e^{-2\gamma_l h_1})} + \rho_l \frac{1 + e^{-2\gamma_l h_2}}{\gamma_l (1 - e^{-2\gamma_l h_2})}\right) \omega^2 w$$
Mô hình này tạo nên sự chuyển dịch giác hải (paradigm shift) từ các phép gần đúng số truyền thống sang phương pháp giải tích chính xác hóa trường chuyển vị động, giải quyết triệt để sự ghép nối phức tạp giữa các mode dao động màng - uốn - xoắn của tấm angle-ply và cross-ply.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp thành công ba lý thuyết cơ học nền tảng:
- Lý thuyết cơ học vật liệu composite nhiều lớp dị hướng (Laminated Composite Mechanics) của Jones (1999) và Reddy (2004).
- Lý thuyết thủy động lực học trường thế Bernoulli - Laplace cho chất lỏng lý tưởng không nén được, không nhớt và chuyển động phi xoáy.
- Lý thuyết nền đàn hồi hai tham số Pasternak - Winkler mô tả tương tác tiếp xúc nền biến dạng cắt ngang.
Điều kiện biên áp dụng trong khung phân tích được tổng quát hóa cho mọi tổ hợp biên tự do (Free - TD), khớp (Simply supported - K), và ngàm (Clamped - N). Đối với tấm chữ nhật, nghiệm được khai triển qua chuỗi Lévy biến tính, chuyển bài toán vi phân đạo hàm riêng hai chiều thành hệ phương trình vi phân thường cấp 10 đối với biến không gian $x$:
$$\frac{d\mathbf{y}_m(x)}{dx} = \mathbf{A}_m(\omega) \mathbf{y}_m(x)$$
Ma trận truyền $\mathbf{T}_m(\omega) = e^{\mathbf{A}_m a}$ được tính toán giải tích thông qua hàm mũ ma trận, cho phép thiết lập ma trận độ cứng động lực $\mathbf{K}(\omega)$ liên hệ trực tiếp giữa véctơ lực tổng quát và véctơ chuyển vị tại các biên:
$$\mathbf{K}(\omega) \mathbf{U}_m = \mathbf{F}_m$$
Biên giới hạn lý thuyết (boundary conditions & scope) được quy định tường minh: kết cấu làm việc trong miền đàn hồi tuyến tính, biến dạng bé, trực hướng từng lớp, chất lỏng tiếp xúc đồng bộ (không xuất hiện hiện tượng tách dòng hoặc sủi bọt khí - cavitation).
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ủ triệt để thế giới quan thực chứng (positivism) kết hợp chặt chẽ giữa suy diễn toán học hình thức (deductive mathematical modeling), mô phỏng số chính xác cao (advanced scientific computing) và thực nghiệm kiểm chứng độc lập (empirical validation).
SƠ ĐỒ THIẾT KẾ NGHIÊN CỨU
Thiết kế nghiên cứu đa tầng (multi-level research design) bao gồm 4 cấu trúc đối tượng khảo sát:
- Tầng 1: Tấm phẳng composite ngập hoàn toàn hoặc một phần trong nước, tấm trên nền Pasternak/Winkler.
- Tầng 2: Vỏ trụ tròn xoay composite lớp chứa chất lỏng có mặt thoáng.
- Tầng 3: Vỏ nón cụt composite lớp chịu điều kiện biên ngàm - tự do, khớp - khớp chứa các mức chất lỏng $H/L$ từ $0%$ đến $100%$.
- Tầng 4: Thực nghiệm chế tạo mẫu vật liệu composite sợi thủy tinh/polyester và đo phổ đáp ứng dao động.
Quy trình nghiên cứu rigorous
Quy trình chế tạo mẫu và thu thập dữ liệu thực nghiệm được chuẩn hóa nghiêm ngặt:
- Chế tạo mẫu: Hợp tác thực hiện cùng Viện Nghiên cứu và Chế tạo Tàu thủy (UNINSHIP) - Đại học Thủy sản Nha Trang. Vật liệu sử dụng là vải sợi thủy tinh dệt kết hợp nhựa polyester không no, gia công bằng phương pháp đắp tay định hình qua khuôn trụ và khuôn nón kim loại chính xác cao.
- Phân nhóm mẫu:
- 03 mẫu vỏ trụ tròn composite: Ký hiệu C1 ($L=600\text{ mm}$, $R=100\text{ mm}$, $h=2.2\text{ mm}$), C2 ($L=800\text{ mm}$, $R=100\text{ mm}$, $h=2.2\text{ mm}$), C3 ($L=1000\text{ mm}$, $R=100\text{ mm}$, $h=2.2\text{ mm}$).
- 03 mẫu vỏ nón cụt composite: Ký hiệu N1 ($L=600\text{ mm}$, $R_1=80\text{ mm}$, $R_2=140\text{ mm}$, $h=2.2\text{ mm}$), N2 ($L=800\text{ mm}$, $R_1=80\text{ mm}$, $R_2=160\text{ mm}$, $h=2.2\text{ mm}$), N3 ($L=1000\text{ mm}$, $R_1=80\text{ mm}$, $R_2=180\text{ mm}$, $h=2.2\text{ mm}$).
- Quy trình đo đạc: Mẫu được liên kết ngàm một đầu vào bệ thép khối lượng lớn thông qua vành đai kẹp bu-lông cường độ cao, đầu còn lại để tự do. Kích thích dao động bằng búa gõ xung lực (impact hammer) có gắn cảm biến lực tại các vị trí kinh tuyến xác định. Tín hiệu gia tốc dao động được thu nhận qua đầu đo áp điện gắn tại thành vỏ, truyền qua bộ chuyển đổi tín hiệu sang máy tính và xử lý phổ trong miền tần số bằng thuật toán biến đổi Fourier nhanh (Fast Fourier Transform - FFT).
Data và phân tích
Quá trình phân tích số và kiểm chứng độ vững (robustness checks) được thực hiện với các thông số vật liệu chi tiết:
- Tấm kim loại kiểm chuẩn: Thép đàn hồi $E = 207\text{ GPa}$, $\rho = 7850\text{ kg/m}^3$, $\nu = 0.3$; ngâm trong bể nước $\rho_n = 1000\text{ kg/m}^3$ kích thước $5\text{ m} \times 5\text{ m} \times 5\text{ m}$, tỷ số chiều dày $h/a = 0.1$.
- Vật liệu Composite thực nghiệm: Composite sợi thủy tinh/Polyester có $E_1 = 14.5\text{ GPa}$, $E_2 = 5.2\text{ GPa}$, $G_{12} = 2.1\text{ GPa}$, $\nu_{12} = 0.28$, $\rho = 1750\text{ kg/m}^3$.
- Vật liệu Composite tính toán số: Hợp chất Graphite/Epoxy cao cấp có $E_1 = 138\text{ GPa}$, $E_2 = 8.96\text{ GPa}$, $G_{12} = 7.1\text{ GPa}$, $G_{13} = 7.1\text{ GPa}$, $G_{23} = 2.84\text{ GPa}$, $\nu_{12} = 0.3$, $\rho = 1580\text{ kg/m}^3$.
| Trường hợp khảo sát |
Cấu hình / Thông số |
Tần số PTLT (Hz / Không thứ nguyên) |
Đối chứng Giải tích / FEM |
Sai số tương đối (%) |
| Tấm thép ngập nước (h1/a=0.5) |
Mindlin, Khớp 4 cạnh |
$\bar{\omega} = 36.425$ |
Hashemi (2012): $36.421$ |
$+0.011%$ |
| Tấm thép trong không khí |
SSSS, $h/a=0.1$ |
$\bar{\omega} = 48.338$ |
Hashemi (2012): $48.338$ |
$0.000%$ |
| Tấm Composite [0/90/90/0] |
SSSS, Graphite/Epoxy |
$\bar{\omega} = 12.148$ |
Reddy (2004): $12.147$ |
$+0.008%$ |
| Vỏ trụ C1 (Khô, H/L=0) |
Ngàm - Tự do, Thực nghiệm |
$f = 118.5\text{ Hz}$ |
Thực nghiệm FFT: $122.0\text{ Hz}$ |
$-2.87%$ |
| Vỏ trụ C1 (Đầy nước, H/L=1) |
Ngàm - Tự do, Thực nghiệm |
$f = 46.2\text{ Hz}$ |
Thực nghiệm FFT: $48.5\text{ Hz}$ |
$-4.74%$ |
| Vỏ nón N1 (Khô, H/L=0) |
$\alpha=6^\circ$, Ngàm - Tự do |
$f = 134.8\text{ Hz}$ |
Thực nghiệm FFT: $138.0\text{ Hz}$ |
$-2.32%$ |
| Vỏ nón N1 (Đầy nước, H/L=1) |
$\alpha=6^\circ$, Ngàm - Tự do |
$f = 52.6\text{ Hz}$ |
Thực nghiệm FFT: $55.0\text{ Hz}$ |
$-4.36%$ |
So sánh về hiệu năng tính toán trên cấu hình máy tính tiêu chuẩn: Khi phân tích tấm composite lệch trục $[0^\circ/45^\circ/0^\circ/45^\circ]$ biên TD-K-TD-K, phương pháp phần tử liên tục (PTLT) chỉ mất 0.82 giây với 01 phần tử duy nhất, trong khi phương pháp phần tử hữu hạn (PTHH) trong ANSYS với lưới chia hội tụ $40 \times 40$ phần tử tiêu tốn 36.45 giây mà vẫn tồn tại sai số tương đối 1.24%.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Quy luật suy giảm tần số phi tuyến theo mức ngập chất lỏng: Khi tỷ số mức nước $H/L$ tăng từ $0$ (vỏ khô) lên $1.0$ (chứa đầy nước), tần số dao động riêng cơ bản của cả vỏ trụ tròn và vỏ nón cụt composite sụt giảm mạnh mẽ từ $58%$ đến $62%$. Tốc độ suy giảm diễn ra nhanh nhất trong giai đoạn $0 \le H/L \le 0.6$, sau đó tiệm cận giá trị bão hòa khi $H/L \to 1.0$.
Tần số riêng f (Hz)
0.0 0.2 0.4 0.6 0.8 1.0
- Hiệu ứng gia cường của nền đàn hồi hai tham số Pasternak: Nền đàn hồi Pasternak làm tăng đáng kể tần số dao động tự do không thứ nguyên của tấm composite. Hệ số cắt của nền ($r_2$) có mức độ ảnh hưởng vượt trội so với hệ số nén Winkler thuần túy ($r_1$), giúp bù đắp lên tới $40%$ lượng suy giảm độ cứng do hiệu ứng tải trọng khối lượng nước gia tán gây ra.
- Đặc trưng cấu hình phân lớp và góc nghiêng vỏ nón: Với vỏ nón cụt composite, khi tăng góc nghiêng $\alpha$ từ $15^\circ$ lên $60^\circ$, tần số dao động riêng ở các mode sóng vòng $n \ge 2$ tăng mạnh do độ cứng hình học tiếp tuyến tăng lên. Cấu hình phân lớp đối xứng đúng trục $[0^\circ/90^\circ/90^\circ/0^\circ]$ thể hiện khả năng kháng uốn vượt trội và ít bị suy giảm tần số trong môi trường chất lỏng hơn cấu hình lệch trục phản xứng $[45^\circ/-45^\circ/45^\circ/-45^\circ]$.
- Sự đồng dạng dạng mode dao động: Mặc dù tần số dao động riêng bị suy giảm nghiêm trọng do tương tác với chất lỏng, hình dạng hình học của các mode dao động cơ bản ($m=1, n=1, 2, 3$) của kết cấu ngập trong chất lỏng hoàn toàn đồng dạng (isomorphic) với mode dao động trong không khí, minh chứng cho tính bảo toàn phân bố năng lượng biến dạng của hệ đàn hồi.
Implications đa chiều
- Về lý thuyết: Đặt nền móng vững chắc cho trường phái tính toán bán giải tích trong cơ học FSI, mở rộng biên giới của lý thuyết Dynamic Stiffness Method cho các kết cấu vỏ tròn xoay phức tạp.
- Về phương pháp luận: Cung cấp một giải pháp thay thế hoàn hảo cho phương pháp phần tử hữu hạn truyền thống, xóa bỏ định kiến về việc phải chia nhỏ lưới đối với các bài toán dao động bậc cao.
- Về ứng dụng thực tiễn: Cung cấp bộ công cụ tính toán chính xác và tức thời cho các kỹ sư đóng tàu, thiết kế bồn bể composite chứa hóa chất (như các sản phẩm thực tế của UNINSHIP và Crown), ngăn ngừa triệt để hiện tượng cộng hưởng nguy hiểm gây phá hủy công trình khi vận hành trong môi trường chất lỏng.
Limitations và Future Research
Mặc dù đạt được những bước tiến đột phá, công trình vẫn tồn tại một số giới hạn nghiên cứu cần được nhìn nhận khách quan:
- Giả thiết chất lỏng: Mô hình giả định chất lỏng là lý tưởng (không nhớt, không nén được, không xoáy). Giả thiết này chưa phản ánh được hiện tượng tắt dần dao động do ma sát nhớt (viscous damping) và hiện tượng sóng va (sloshing) tự do phi tuyến ở bề mặt thoáng khi kích động mạnh.
- Miền làm việc kết cấu: Mới giới hạn ở biến dạng bé tuyến tính hình học, chưa xét tới các hiệu ứng phi tuyến độ võng lớn von Kármán hoặc ảnh hưởng của vết nứt, phân lớp (delamination) cục bộ trong quá trình chế tạo thực nghiệm.
- Điều kiện biên thực nghiệm: Mối liên kết ngàm trong thực nghiệm tại xưởng vẫn tồn tại độ mềm đàn hồi nhất định tại chân đế, dẫn đến độ lệch nhỏ ($2.3% - 4.8%$) so với điều kiện ngàm lý thuyết tuyệt đối trong tính toán số.
Chương trình nghiên cứu tương lai (Future research directions):
- Phát triển phần tử liên tục cho bài toán tương tác chất lỏng phi tuyến có tính đến độ nhớt của môi trường (phương trình Navier-Stokes ghép nối).
- Tích hợp phân tích dao động ngẫu nhiên và kiểm soát rung động chủ động sử dụng các lớp vật liệu áp điện thông minh (Smart Piezoelectric Layers).
- Nghiên cứu dao động của kết cấu vỏ composite chức năng biến đổi cơ tính (Functionally Graded Materials - FGM) chứa chất lỏng dưới tác động đồng thời của trường nhiệt - cơ.
- Mở rộng bài toán xác định đáp ứng động lực học của kết cấu vỏ chịu tải trọng sóng va đập thủy động lực học ngẫu nhiên.
Tác động và ảnh hưởng
Luận án của tác giả Tạ Thị Hiền đã tạo dựng những giá trị khoa học và ứng dụng thực tiễn sâu rộ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í khoa học uy tín chuyên ngành Cơ học và các kỷ yếu hội nghị cơ học toàn quốc, đóng góp trực tiếp vào hệ thống giáo trình cao học và nghiên cứu sinh tại Đại học Bách khoa Hà Nội và Viện Cơ học Việt Nam.
- Công nghiệp & Quốc phòng: Thuật toán và phần mềm được ứng dụng trực tiếp trong tính toán thiết kế vỏ tàu composite hai thân, bàn đẩy tàu cánh ngầm, vòm che máy bay quân sự, các bể nuôi trồng thủy sản trên biển của UNINSHIP và hệ thống bồn chứa hóa chất chuyên dụng của Công ty Crown.
- Lợi ích xã hội & Môi trường: Tối ưu hóa trọng lượng và kết cấu composite giúp tăng tuổi thọ công trình chứa chất lỏng lên trên 30 năm, chống ăn mòn hóa chất và nước biển tuyệt đối, tiết kiệm hàng triệu USD chi phí bảo dưỡng định kỳ và bảo vệ môi trường sinh thái biển.
Đối tượng hưởng lợi
- Nghiên cứu sinh & Học viên Cao học ngành Cơ học/Xây dựng/Giao thông: Tiếp cận một phương pháp số bán giải tích chuẩn mực (PTLT/DSM), làm chủ kỹ thuật lập trình ma trận truyền và phương pháp xử lý toán học cho hệ phương trình vi phân dao động ghép nối.
- Các Giáo sư & Nhà nghiên cứu cao cấp: Sở hữu tài liệu tham khảo giá trị để phát triển các lý thuyết vỏ bậc cao (HSDT) và mở rộng sang các dòng vật liệu mới như FGM, CNTs-reinforced composites tiếp xúc chất lỏng.
- Kỹ sư R&D công nghiệp Tàu thủy & Hàng không Vũ trụ: Ứng dụng ngay 4 bộ chương trình MATLAB (
VplateF, VplateEF, VcylF, VconF) vào quy trình mô phỏng thiết kế kết cấu vỏ chứa nhiên liệu lỏng và thân tàu composite.
- Các nhà hoạch định chính sách & Tiêu chuẩn hóa công trình: Có cơ sở khoa học định lượng vững chắc để xây dựng các quy chuẩn, tiêu chuẩn kỹ thuật quốc gia về an toàn động lực học cho các công trình composite ngập nước.
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 độc đáo nhất là việc mở rộng thành công Phương pháp Độ cứng động (Dynamic Stiffness Method / PTLT) trên nền tảng Lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin cho bài toán ghép nối đa trường: kết cấu composite dị hướng lớp - trường áp suất thủy động Bernoulli/Laplace - nền đàn hồi Pasternak không thuần nhất. Công trình đã tổng quát hóa ma trận hệ số động lực học $\mathbf{A}_m$ kích thước $10 \times 10$, tích hợp đại lượng khối lượng gia tán ảo $m^*$ trực tiếp vào phương trình vi phân trạng thái chuyển vị - nội lực, thiết lập nghiệm giải tích chuẩn mực cho hàm thế vận tốc có mặt thoáng và đáy cứng.
2. Sự cách tân trong phương pháp nghiên cứu khi so sánh với các công trình quốc tế tiền nhiệm?
So với phương pháp phần tử hữu hạn truyền thống (Fu & Price, 1987; Matthew et al., 2013) vốn đòi hỏi chia hàng nghìn phần tử và dễ mắc sai số cắt số học, phương pháp PTLT của luận án biểu diễn trường chuyển vị bằng nghiệm động lực học chính xác, chỉ cần duy nhất 01 phần tử cho toàn bộ tấm hoặc đoạn vỏ mà vẫn triệt tiêu hoàn toàn sai số rời rạc hóa. So với phương pháp giải tích Navier cổ điển (Korosh Khorshid & Farhadi, 2013) chỉ giải được biên khớp 4 cạnh (SSSS), phương pháp của luận án xử lý chính xác tuyệt đối mọi điều kiện biên phức tạp (Ngàm - Khớp - Tự do) thông qua ma trận truyền $\mathbf{T}_m(\omega) = e^{\mathbf{A}_m a}$.
3. Phát hiện số liệu gây bất ngờ nhất trong nghiên cứu thực nghiệm và tính toán?
Phát hiện bất ngờ nhất là hiện tượng sụt giảm đột ngột của tần số dao động riêng ở các mức nước thấp ($0 \le H/L \le 0.4$): tần số riêng giảm tới hơn $45%$ tổng lượng suy giảm của toàn dải ngập, trong khi khoảng tăng mức nước từ $0.8 \le H/L \le 1.0$ chỉ làm tần số giảm thêm dưới $5%$. Điều này chứng minh rằng phần chất lỏng ở đáy vỏ đóng góp khối lượng gia tán động lực học lớn nhất do biên độ dao động uốn tại khu vực này đạt cực đại đối với liên kết ngàm - tự do.
(Vùng biến thiên cực đại) (Vùng bão hòa)
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Luận án cung cấp quy trình tái lập vô cùng chi tiết: đầy đủ các hệ thức giải tích của ma trận $\mathbf{A}_m$, quy tắc chuyển đổi biến trạng thái, cấu trúc giải thuật trên MATLAB, thông số hình học và cơ lý tính của mẫu thí nghiệm sợi thủy tinh/polyester, sơ đồ bố trí cảm biến gia tốc và thông số búa gõ kích động FFT. Mọi nhà nghiên cứu độc lập đều có thể tái lập chính xác các đường cong đáp ứng tần số và kết quả đo kiểm.
5. Chương trình nghị sự nghiên cứu 10 năm được vạch ra như thế nào?
Lộ trình phát triển học thuật bao gồm: (1) Hoàn thiện phần tử liên tục phi tuyến hình học cho kết cấu vỏ mỏng chịu tải trọng động lực học ngẫu nhiên; (2) Mở rộng thuật toán PTLT cho các dạng vỏ phức tạp như vỏ composite hyperbolic paraboloid và vỏ nón cụt ghép tầng; (3) Tích hợp trường đa vật lý: Nhiệt - Điện - Từ - Thủy động lực học cho vật liệu FGM và vật liệu thông minh; (4) Chuyển giao các module tính toán thành phần mềm thương mại độc lập phục vụ công nghiệp đóng tàu và công trình biển.
Kết luận
Luận án tiến sĩ của tác giả Tạ Thị Hiền đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn cốt lõi:
- Thiết lập thành công mô hình toán học giải tích - số chính xác cao dựa trên lý thuyết biến dạng cắt bậc nhất FSDT và phương pháp Phần tử liên tục (PTLT) để giải quyết trọn vẹn bài toán dao động tự do của tấm composite ngập nước, tấm trên nền đàn hồi và vỏ tròn xoay composite chứa chất lỏng.
- Xây dựng hoàn chỉnh hệ thống 4 chương trình tính toán số chuyên dụng trên môi trường MATLAB (
VplateF, VplateEF, VcylF, VconF), mang lại hiệu năng vượt bậc về tốc độ tính toán và độ chính xác so với các phần mềm thương mại PTHH (ANSYS).
- Khảo sát và lượng hóa toàn diện các quy luật cơ học: xác định rõ ảnh hưởng của mức nước $H/L$, tỷ số hình học $h/a, R/h$, góc nghiêng nón $\alpha$, đặc tính dị hướng $E_1/E_2$, cấu hình phân lớp $[0^\circ/90^\circ]$, và hai hệ số nền Pasternak $r_1, r_2$ đến phổ tần số dao động riêng.
- Thực hiện thành công nghiên cứu thực nghiệm chuyên sâu: thiết kế, chế tạo 6 cụm mẫu vỏ composite và tiến hành đo phổ rung động FFT trong môi trường chất lỏng biến thiên, cung cấp tập dữ liệu thực nghiệm chuẩn xác khẳng định độ tin cậy của thuật toán lý thuyết.
- Mở ra 3 hướng nghiên cứu khoa học mới: Động lực học công trình FSI phi tuyến cho vật liệu composite phân lớp chức năng FGM; Phân tích dao động tắt dần của kết cấu thông minh tiếp xúc chất lỏng nhớt; Thiết kế tối ưu hóa độ bền dao động cho các công trình công nghiệp và quốc phòng hiện đại.
