Tổng quan về luận án

Dao động cơ học là hiện tượng phổ biến và thường mang tính phá hủy đối với các hệ kết cấu kỹ thuật như công trình cao tầng, tháp cầu dây văng, giàn khoan biển, cũng như các thiết bị máy móc chính xác trong công nghiệp và giao thông vận tải. Hiện tượng cộng hưởng xảy ra khi tần số kích động từ môi trường (gió, động đất, sóng biển, tải trọng động) trùng với tần số dao động riêng của kết cấu, gây ra biến dạng uốn lớn, phá hủy mỏi và làm suy giảm nghiêm trọng tuổi thọ làm việc. Để giải quyết triệt để vấn đề này, luận án tiến sĩ cơ học của tác giả Vũ Đức Phúc (chuyên ngành Cơ học, mã số 9440109) với đề tài "Điều khiển tối ưu dao động bằng kết hợp nhiều bộ giảm chấn động lực", dưới sự hướng dẫn khoa học của GS.TSKH Nguyễn Văn Khang và PGS.TS Nguyễn Phong Điền tại Trường Đại học Bách khoa Hà Nội (2020), đã phát triển hệ thống các giải pháp điều khiển thụ động tiên tiến thông qua việc tích hợp đa bộ giảm chấn động lực (Dynamic Vibration Absorbers - DVAs).

Tính tiên phong của luận án thể hiện ở việc chuyển dịch từ mô hình giảm chấn đơn lẻ (Single Dynamic Vibration Absorber - SDVA / STMD) với nguy cơ mất ổn định cao khi lệch pha tần số (detuning) sang mô hình tổ hợp nhiều bộ giảm chấn động lực đơn tần số và đa tần số (Multiple Dynamic Vibration Absorbers - MDVA / MTMD / MTLD) lắp đặt trên cả hệ tập trung một bậc tự do có cản và hệ liên tục dầm Euler - Bernoulli.

Khoảng trống nghiên cứu (research gap) mà luận án giải quyết bao gồm:

  1. Thiếu vắng các giải pháp giải tích và bán giải tích tường minh để tối ưu hóa đồng thời các tham số của nhiều bộ DVA đơn tần số lắp trên hệ chính có kể đến lực cản nhớt ($c_s \neq 0$).
  2. Sự hạn chế của các phương pháp tối ưu hóa kinh nghiệm cổ điển (như phương pháp của Fujino và đồng nghiệp) khi thiết kế bộ giảm chấn chất lỏng đa tần số (MTLD) và hệ giảm chấn khối lượng đa tần số (MTMD), đòi hỏi một thuật toán tối ưu hóa mạnh mẽ có khả năng xử lý phiếm hàm mục tiêu đa chiều.
  3. Thiếu khung lý thuyết kết hợp mô hình hóa liên tục bằng phương pháp biến phân Ritz - Galerkin và tối ưu hóa phân bố không gian cho dầm Euler - Bernoulli có cản chịu kích động điều hòa.

Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập đồng bộ:

  • RQ1: Làm thế nào để mở rộng lý thuyết kinh điển của Den Hartog để xác định tham số tối ưu ($\alpha_{opt}, \xi_{opt}$) cho hệ nhiều bộ giảm chấn đơn tần số gắn trên hệ có cản?
    • Hypothesis H1: Áp dụng phương pháp tuyến tính hóa tương đương kết hợp tiêu chuẩn bình phương tối thiểu (BPTT) và tiêu chuẩn đối ngẫu sẽ quy đổi chính xác hệ có cản về hệ không cản tương đương, triệt tiêu biên độ dao động tại tần số cộng hưởng.
  • RQ2: Thuật toán nào cho phép tối ưu hóa các tham số phân tán tần số của hệ MTMD/MTLD nhằm mở rộng dải tần dập tắt dao động và tăng tính bền vững (robustness)?
    • Hypothesis H2: Ứng dụng phương pháp thiết kế thực nghiệm Taguchi với mảng trực giao và tỷ số tín hiệu trên nhiễu (SNR) kết hợp hàm mục tiêu có trọng số sẽ tìm được nghiệm tối ưu toàn cục nhanh hơn các thuật giải tiến hóa truyền thống.
  • RQ3: Vị trí lắp đặt và phân bố tham số của các TMD ảnh hưởng như thế nào đến sự dập tắt các dạng dao động riêng của dầm liên tục Euler - Bernoulli?
    • Hypothesis H3: Việc định vị TMD tại các điểm bụng sóng của dạng dao động riêng thứ nhất thông qua khai triển chuỗi hàm trực giao Ritz - Galerkin sẽ mang lại hiệu quả giảm chấn vượt trội trên 90%.

Khung lý thuyết của luận án được xây dựng dựa trên sự kết hợp giữa Cơ học giải tích (phương trình vi phân Lagrange loại 2), Cơ học vật rắn biến dạng (lý thuyết dầm Euler - Bernoulli), Lý thuyết tối ưu hóa động lực học và Phương pháp thiết kế tham số Taguchi. Phạm vi nghiên cứu bao quát từ mô hình toán học giải tích, tính toán mô phỏng số trên MATLAB, đến kiểm nghiệm thực nghiệm trên hệ thống bàn rung điện động TV50350/LS-120 với dầm đồng đỏ có chiều dài $L = 0{,}639\text{ m}$, khối lượng $m_b = 1{,}28\text{ kg}$, đạt hiệu quả triệt tiêu biên độ dao động lên đến $96{,}84%$.

Literature Review và Positioning

Nghiên cứu về bộ giảm chấn động lực bắt đầu từ phát minh của Watt (1883) nhằm chống rung lắc cho tàu thủy và được Frahm (1909) nhận bằng sáng chế cho mô hình bộ giảm chấn không cản gồm khối lượng $m_a$ và lò xo $k_a$. Trong mô hình này, khi tần số ngoại lực $\Omega$ trùng với tần số riêng $\omega_a = \sqrt{k_a / m_a}$, dao động của hệ chính bị triệt tiêu hoàn toàn, tức $x_s = 0$. Tuy nhiên, hạn chế cốt tử của hệ DVA không cản là xuất hiện hai đỉnh cộng hưởng mới ở hai bên tần số $\omega_a$, khiến hệ thống trở nên cực kỳ nguy hiểm nếu tần số kích động thay đổi.

Bước đột phá lý thuyết xuất hiện khi Den Hartog (1928) bổ sung phần tử cản nhớt $c_a$ và phát triển "Lý thuyết hai điểm cố định" (Fixed Points Theory). Ông chứng minh rằng mọi đường cong đáp ứng tần số của hệ chính không cản khi thay đổi hệ số cản $\xi_a$ đều đi qua hai điểm cố định $P$ và $Q$. Bằng cách cân bằng biên độ đáp ứng tại hai điểm này, Den Hartog đã thiết lập công thức giải tích kinh điển xác định tỷ lệ tần số tối ưu $\alpha_{opt}$ và tỷ số cản tối ưu $\xi_{opt}$: $$\alpha_{opt} = \frac{1}{1 + \mu_a}, \quad \xi_{opt} = \sqrt{\frac{3\mu_a}{8(1 + \mu_a)}}$$ trong đó $\mu_a = m_a / m_s$ là tỷ lệ khối lượng.

Biên độ |H(Ω)|

Trong các giai đoạn tiếp theo, các dòng nghiên cứu quốc tế tập trung vào việc khắc phục hai hạn chế lớn của lý thuyết Den Hartog:

  1. Hệ chính có cản ($c_s \neq 0$): Brock (1946), Warburton (1982), Asami và Nishihara (2003) nhận thấy khi hệ chính có cản, hai điểm $P$ và $Q$ không còn cố định tuyệt đối. Các tác giả phải sử dụng các phương pháp tối ưu xấp xỉ hoặc phương pháp số phức tạp.
  2. Hệ nhiều bộ giảm chấn (MTMD / MTLD): Xu, Kareem (1993), Sun, Kareem (1995), và đặc biệt là nhóm nghiên cứu của Fujino và đồng nghiệp (1992) đã tiên phong phân tích hệ nhiều bộ giảm chấn chất lỏng (TLD/MTLD) bằng phương pháp thực nghiệm và hồi quy bán giải tích. Tuy nhiên, các công thức của Fujino chủ yếu dựa trên phân bố tần số đều và giả định kích động băng rộng, thiếu tính linh hoạt khi áp dụng cho kết cấu thực có cản chịu kích động điều hòa đơn tần hoặc đa tần.

Tại Việt Nam, các công trình của GS.TSKH Nguyễn Văn Khang và cộng sự đã nghiên cứu sâu về phương pháp giải tích và phương pháp tuyến tính hóa tương đương cho các hệ cơ học phi tuyến. Luận án của Vũ Đức Phúc đã định vị chính xác điểm giao thoa giữa các trường phái nghiên cứu: kế thừa tính tường minh giải tích của Den Hartog, tích hợp phương pháp tuyến tính hóa tương đương để giải quyết bài toán hệ có cản, đồng thời thay thế phương pháp giải tích cổ điển bằng thuật toán tối ưu Taguchi cho hệ đa bậc tự do và kết cấu liên tục.

Tiêu chí so sánh Nghiên cứu của Den Hartog (1928) Nghiên cứu của Fujino và CS (1992) Luận án của Vũ Đức Phúc (2020)
Đối tượng kết cấu Hệ 1 bậc tự do không cản ($c_s = 0$) Hệ kết cấu công trình chịu tải ngẫu nhiên Hệ 1-DOF có cản ($c_s \neq 0$) & Dầm Euler - Bernoulli
Loại giảm chấn 01 bộ DVA/TMD đơn lẻ Nhiều bộ TLDs (MTLD) Tổ hợp $N$ bộ DVAs/TMDs/TLDs đơn & đa tần số
Phương pháp tối ưu Giải tích hai điểm cố định $P, Q$ Thực nghiệm & hàm truyền xấp xỉ Tuyến tính hóa tương đương, BPTT & Taguchi SNR
Khả năng triệt tiêu Dải tần rất hẹp, nhạy cảm detuning Dải tần rộng nhưng kích thước thùng lớn Dải tần rộng, tối ưu vị trí không gian, giảm chấn $>96%$

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mang lại bốn đóng góp lý thuyết mang tính đột phá cho chuyên ngành Cơ học ứng dụng và Điều khiển dao động công trình:

  1. Mở rộng lý thuyết hai điểm cố định cho hệ $N$ bộ DVA đơn tần số: Luận án đã chứng minh một cách chặt chẽ rằng khi lắp $N$ bộ giảm chấn động lực có thông số hoàn toàn giống nhau ($m_j = m_a/N, k_j = k_a/N, c_j = c_a/N$) lên hệ chính không cản, đáp ứng tần số của hệ hoàn toàn đồng nhất với hệ lắp 1 bộ DVA có tổng khối lượng $m_a$. Tỷ lệ tần số tối ưu và tỷ số cản tối ưu của mỗi bộ DVA thành phần tuân theo công thức: $$\alpha_{opt, j} = \frac{1}{1 + \mu_a}, \quad \xi_{opt, j} = \sqrt{\frac{3\mu_a}{8(1 + \mu_a)}}$$ với $\mu_a = \sum_{j=1}^N m_j / m_s$.

  2. Thiết lập cơ sở giải tích quy đổi hệ có cản bằng phương pháp Tuyến tính hóa tương đương: Đối với hệ chính có cản ($c_s \neq 0$), luận án áp dụng nguyên lý công khả dĩ và phương pháp tuyến tính hóa tương đương, đưa ra hai tiêu chuẩn thiết kế tối ưu mới:

    • Tiêu chuẩn Bình phương tối thiểu (BPTT): Cực tiểu hóa sai lệch bình phương giữa năng lượng hao tán thực tế và năng lượng hao tán của hệ không cản quy đổi.
    • Tiêu chuẩn Đối ngẫu (Dual Criterion): Cân bằng đồng thời thế năng đàn hồi và hàm hao tán Rayleigh trong chu kỳ dao động bình ổn.

  3. Mô hình toán học tương đương giữa TLD và TMD: Luận án hoàn thiện khung lý thuyết quy đổi thủy động lực học sóng phi tuyến trong thùng chất lỏng (Tuned Liquid Damper - TLD) về mô hình dao động khối lượng - lò xo - cản nhớt (TMD) tương đương: $$m_{TMD} = m \frac{8\tanh(\pi h_0 / l)}{\pi^3 h_0 / l}, \quad k_{TMD} = \omega_{TLD}^2 m_{TMD}, \quad c_{TMD} = 2 \xi_{TLD} \omega_{TLD} m_{TMD}$$ trong đó tần số dao động tự nhiên của sóng chất lỏng bậc nhất là $\omega_{TLD} = \sqrt{\frac{\pi g}{l} \tanh\left(\frac{\pi h_0}{l}\right)}$, và tỷ số cản $\xi_{TLD} = \frac{1}{2h_0}\sqrt{\frac{\nu}{2\omega_{TLD}}}\left(1 + \frac{2h_0}{b} + S\right)$.

  4. Khung lý thuyết không gian - thời gian cho kết cấu dầm liên tục: Thiết lập thành công hệ phương trình vi phân chuyển động ma trận hóa cho dầm Euler - Bernoulli chịu liên kết biên tùy ý (hai đầu bản lề, một đầu ngàm một đầu tự do) có gắn $N$ bộ TMD tại các tọa độ không gian $\eta_j$. Khung phân tích tách biệt tường minh ma trận khối lượng $[\mathbf{M}]$, ma trận cản $[\mathbf{D}]$ và ma trận độ cứng $[\mathbf{K}]$ thông qua cơ sở hàm riêng trực giao.

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

Luận án tích hợp ba lý thuyết nền tảng để xây dựng khung phân tích động lực học kết cấu hoàn chỉnh:

Quy trình phân tích được thiết kế qua các bước logic:

  1. Mô tả hình học và động học: Biểu diễn độ võng của dầm $w(x,t) = \sum_{k=1}^m X_k(x) q_k(t)$, với $X_k(x)$ là các dạng dao động riêng thỏa mãn điều kiện biên hình học.
  2. Thiết lập hàm Lagrange và hàm hao tán Rayleigh: $$T = \frac{1}{2} \int_0^L \rho A \dot{w}^2 dx + \frac{1}{2} \sum_{j=1}^N m_j \dot{u}j^2$$ $$\Pi = \frac{1}{2} \int_0^L EI \left(\frac{\partial^2 w}{\partial x^2}\right)^2 dx + \frac{1}{2} \sum{j=1}^N k_j \left[u_j - w(\eta_j, t)\right]^2$$ $$\Phi = \frac{1}{2} \int_0^L c^{(e)} \dot{w}^2 dx + \frac{1}{2} \int_0^L c^{(i)} I \left(\frac{\partial^2 \dot{w}}{\partial x^2}\right)^2 dx + \frac{1}{2} \sum_{j=1}^N d_j \left[\dot{u}_j - \dot{w}(\eta_j, t)\right]^2$$
  3. Biến đổi Laplace/Fourier: Chuyển hệ phương trình vi phân cấp hai sang miền tần số phức để thu được ma trận hàm đáp ứng tần số $[\mathbf{H}(\Omega)] = [-\Omega^2 \mathbf{M} + i\Omega \mathbf{D} + \mathbf{K}]^{-1}$.
  4. Thiết lập hàm mục tiêu Taguchi: Luận án đề xuất hàm mục tiêu tổng quát có trọng số $f$: $$f = w_1 \max_{\Omega} |H(\Omega)| + w_2 \int_{\Omega_1}^{\Omega_2} |H(\Omega)| d\Omega$$ cho phép cân bằng giữa việc triệt tiêu đỉnh cực đại và làm phẳng toàn bộ vùng tần số lân cận cộng hưởng.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận thực chứng (Positivism) và chủ nghĩa hiện thực phản biện (Critical Realism), kết hợp chặt chẽ giữa suy diễn giải tích toán học, thuật toán tìm kiếm tham số tối ưu và kiểm chứng thực nghiệm vật lý độc lập. Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 (Giải tích thuần túy): Thiết lập nghiệm bình ổn dạng sóng điều hòa $x_s(t) = \hat{x}_s e^{i\Omega t}$ và đạo hàm các công thức biên độ phức.
  • Tầng 2 (Thuật toán Taguchi): Xây dựng quy trình tối ưu hóa số thích nghi bằng mảng trực giao $L_{27}(3^{13})$ cho hệ 5 bộ DVA trên hệ 1 bậc tự do và mảng trực giao $L_{18}(2^1 \times 3^7)$ cho dầm liên tục.
  • Tầng 3 (Mô phỏng động lực học): Lập trình trên môi trường MATLAB/Simulink khảo sát đáp ứng miền thời gian bằng thuật toán Runge - Kutta bậc 4-5 và đáp ứng miền tần số.
  • Tầng 4 (Thực nghiệm vật lý): Kiểm chứng mô hình dầm côngxôn gắn 2 bộ TMD trên bàn rung điều khiển tự động.

Quy trình nghiên cứu rigorous

Quy trình tối ưu hóa tham số dựa trên phương pháp Taguchi được chuẩn hóa qua các giai đoạn nghiêm ngặt:

  1. Khởi tạo mảng trực giao: Xác định các thông số điều khiển gồm độ cứng $k_j$ và hệ số cản $d_j$ cho từng bộ DVA ($j = 1, \dots, N$). Mỗi thông số được chia làm 3 phân mức giá trị (Mức 1, Mức 2, Mức 3).
  2. Tính toán tỷ số tín hiệu trên nhiễu (SNR): Áp dụng đặc trưng chất lượng "Càng nhỏ càng tốt" (Smaller-the-Better): $$\text{SNR}i = -10 \log{10} \left( \frac{1}{n} \sum_{k=1}^n y_{ik}^2 \right) = -10 \log_{10} (H_i^2)$$ trong đó $H_i$ là giá trị hàm mục tiêu tại thí nghiệm thứ $i$.
  3. Phân tích giá trị trung bình (ANOM): Xác định mức tối ưu cho từng tham số điều khiển bằng cách lấy giá trị cực đại của $(\text{SNR}_m)_p$ ứng với phân mức $p$ của tham số $m$.
  4. Thuật toán co cụm không gian tham số: Sau mỗi vòng lặp, nếu chưa thỏa mãn điều kiện hội tụ sai số phương sai $\sigma^2 \le \varepsilon$, thuật toán tự động thu hẹp dải biến thiên xung quanh điểm tối ưu cục bộ và tạo bảng phân mức mới cho vòng lặp tiếp theo.

Data và phân tích

Thiết lập thực nghiệm được triển khai tại phòng thí nghiệm chuyên sâu với hệ thống phần cứng và phần mềm đo kiểm chuẩn quốc tế:

  • Cơ hệ mẫu: Thanh dầm bằng đồng đỏ tiết diện chữ nhật: Chiều dài $L = 0{,}639\text{ m}$, chiều rộng $b = 0{,}04\text{ m}$, chiều dày $h = 0{,}005\text{ m}$, khối lượng $m_b = 1{,}28\text{ kg}$, mômen quán tính tiết diện $I = 4{,}1667 \times 10^{-10}\text{ m}^4$, môđun đàn hồi xác định bằng thực nghiệm $E = 77{,}402 \times 10^9\text{ N/m}^2$, hệ số cản trong $c^{(i)} = 0{,}002\text{ Ns/m}$. Tần số dao động riêng thứ nhất đo được là $f_1 = 5{,}5\text{ Hz}$.
  • Hệ thống tạo rung: Bàn rung điện động TIRA TV50350/LS-120 (lực kích rung cực đại 18 N, dải tần 2 - 4000 Hz, hành trình dịch chuyển $\pm 2\text{ cm}$, tải trọng mẫu 3 kg) tích hợp hệ thống bù tải khí nén.
  • Khối điều khiển & khuếch đại: Bộ điều khiển số đa kênh VR-8500 kết nối Ethernet, dải động 100 dB, bộ biến đổi 24-bit A/D và D/A, điều khiển qua phần mềm View Vibration; Bộ khuếch đại công suất tương tự BAA 120 TIRA công suất định mức 120 W.
  • Hệ thu thập và xử lý tín hiệu: Cảm biến gia tốc áp điện siêu nhạy của hãng Brüel & Kjær, hệ thống thu nhận dữ liệu đa kênh và phần mềm phân tích phổ Pulse LabShop phiên bản 10.5 với độ phân giải lên đến 4 tỷ điểm đo/kênh.

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

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

  1. Hiệu năng dập tắt dao động vượt trội của hệ 2 TMD lắp trên dầm: Thực nghiệm đo gia tốc tại vị trí $x = 47{,}5\text{ cm}$ khi dầm chịu kích động điều hòa tại ngàm ở tần số cộng hưởng $f = 5{,}5\text{ Hz}$ với các biên độ kích động khác nhau ($y_0 = 0{,}5\text{ mm}; 1{,}0\text{ mm}; 1{,}5\text{ mm}$) đã chứng minh hiệu quả giảm chấn kinh ngạc của hệ 2 TMD (TMD1 tại $35\text{ cm}$, TMD2 tại $63{,}9\text{ cm}$, $m_j = 0{,}09\text{ kg}$, $k_j = 100\text{ N/m}$, $d_j = 0{,}00101\text{ Ns/m}$):
    • Khi $y_0 = 0{,}5\text{ mm}$: Gia tốc giảm từ $8{,}77\text{ m/s}^2$ (không TMD) xuống $0{,}33\text{ m/s}^2$ (lý thuyết) và $0{,}307\text{ m/s}^2$ (thực nghiệm), đạt hiệu quả giảm chấn $96{,}24%$ (lý thuyết) và $96{,}84%$ (thực nghiệm), sai số giữa lý thuyết và thực nghiệm chỉ là $0{,}63%$.
    • Khi $y_0 = 1{,}0\text{ mm}$: Gia tốc giảm từ $17{,}53\text{ m/s}^2$ xuống $0{,}66\text{ m/s}^2$ (lý thuyết) và $0{,}634\text{ m/s}^2$ (thực nghiệm), hiệu quả đạt $96{,}24%$ (lý thuyết) và $95{,}99%$ (thực nghiệm), sai số $0{,}26%$.
    • Khi $y_0 = 1{,}5\text{ mm}$: Gia tốc giảm từ $26{,}30\text{ m/s}^2$ xuống $0{,}99\text{ m/s}^2$ (lý thuyết) và $0{,}973\text{ m/s}^2$ (thực nghiệm), hiệu quả đạt $96{,}24%$ (lý thuyết) và $96{,}12%$ (thực nghiệm), sai số $0{,}12%$.
Hiệu quả giảm chấn (%)

  1. Khảo sát tại điểm đo $x = 31\text{ cm}$ với cấu hình TMD phân bố khác: Khi lắp 2 TMD tại tọa độ $40\text{ cm}$ và $63{,}9\text{ cm}$ ($k_j = 143{,}9\text{ N/m}$), hiệu quả điều khiển dao động vẫn duy trì ở mức xuất sắc:

    • Khi $y_0 = 1{,}0\text{ mm}$: Hiệu quả giảm chấn lý thuyết đạt $90{,}51%$, thực nghiệm đạt $90{,}17%$ (sai số $0{,}37%$).
    • Khi $y_0 = 1{,}5\text{ mm}$: Hiệu quả giảm chấn lý thuyết đạt $90{,}50%$, thực nghiệm đạt $91{,}33%$ (sai số $0{,}91%$).

  2. Vượt trội của thuật toán Taguchi so với phương pháp kinh nghiệm của Fujino: Khi tối ưu hóa hệ 5 bộ DVA đa tần số (MTMD) trên hệ chính 1 bậc tự do, thuật toán Taguchi tìm được bộ thông số phân tán tần số tối ưu giúp đường cong đáp ứng tần số $|H(\Omega)|$ hạ thấp đồng đều trên toàn bộ dải kích động $\Omega \in [60, 90]\text{ rad/s}$. So với phương pháp thực nghiệm của Fujino, phương pháp Taguchi triệt tiêu hoàn toàn hiện tượng nhô cao của hai đỉnh cộng hưởng phụ, giảm biên độ cực đại $|H_{max}|$ thêm $14{,}8%$ với cùng tổng tỷ lệ khối lượng $\mu_a = 5%$.

  3. Tác động của hàm mục tiêu có trọng số ($f$): Việc đưa vào thành phần tích phân năng lượng $\int |H(\Omega)| d\Omega$ trong hàm mục tiêu $f$ ngăn chặn triệt để hiện tượng "tối ưu hóa cục bộ" (over-tuning) tại đúng một điểm tần số mà làm tăng đột biến dao động ở các tần số lân cận.

Cấu hình dầm & Biên độ kích động Gia tốc không TMD ($m/s^2$) Gia tốc có 2 TMD lý thuyết ($m/s^2$) Gia tốc có 2 TMD thực nghiệm ($m/s^2$) Hiệu quả thực nghiệm (%) Sai số lý thuyết - thực nghiệm (%)
$x = 47{,}5\text{ cm}, y_0 = 0{,}5\text{ mm}$ 9,73 0,33 0,307 96,84% 0,63%
$x = 47{,}5\text{ cm}, y_0 = 1{,}0\text{ mm}$ 15,90 0,66 0,634 95,99% 0,26%
$x = 47{,}5\text{ cm}, y_0 = 1{,}5\text{ mm}$ 25,10 0,99 0,973 96,12% 0,12%
$x = 31{,}0\text{ cm}, y_0 = 1{,}0\text{ mm}$ 8,07 0,827 0,793 90,17% 0,37%
$x = 31{,}0\text{ cm}, y_0 = 1{,}5\text{ mm}$ 13,50 1,241 1,170 91,33% 0,91%

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính khả thi của việc giải bài toán tối ưu động lực học đa chiều bằng các phương pháp thiết kế tham số thống kê (Taguchi/SNR) thay vì phụ thuộc vào các công thức giải tích đóng vốn chỉ áp dụng được cho mô hình lý tưởng hóa cực kỳ đơn giản.
  • Về phương pháp luận: Cung cấp quy trình khép kín từ mô hình hóa toán học liên tục (Ritz - Galerkin), tối ưu hóa tham số (Taguchi), mô phỏng số (MATLAB) đến đo kiểm thực nghiệm (Pulse LabShop), có thể áp dụng trực tiếp cho các dạng kết cấu phức tạp hơn như tấm, vỏ, khung không gian.
  • Về ứng dụng thực tiễn: Cho phép ngành xây dựng và cơ khí chế tạo thay thế các khối cầu giảm chấn TMD khổng lồ đơn lẻ (vốn chiếm nhiều diện tích hữu dụng trên đỉnh tháp và tạo áp lực cục bộ lớn lên khung chịu lực) bằng hàng loạt các bộ TMD nhỏ hoặc các két nước sinh hoạt TLD phân tán tại nhiều tầng, dễ lắp đặt, bảo trì và chế tạo hàng loạt.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án cũng chỉ rõ các giới hạn mang tính học thuật:

  1. Giả định tuyến tính: Mô hình dầm Euler - Bernoulli và các bộ TMD đều giả định biến dạng nhỏ và lực cản tuyến tính nhớt. Hiện tượng phi tuyến hình học khi dầm dao động với biên độ lớn hoặc lực cản phi tuyến của dòng chất lỏng chảy rối trong TLD chưa được mô hình hóa toàn diện.
  2. Bỏ qua cản khí động ngoài trong mô hình dầm gốc: Trong phân tích lý thuyết dầm không lắp TMD, việc bỏ qua cản ngoài $c^{(e)}$ dẫn đến sai số lý thuyết so với thực nghiệm ban đầu xấp xỉ $10{,}94%$ (ở $y_0 = 0{,}5\text{ mm}$), dù sau đó khi lắp TMD sai số này giảm xuống dưới $1%$.
  3. Kích động điều hòa một chiều: Luận án tập trung chủ yếu vào kích động điều hòa đơn tần số và dải hẹp, chưa khảo sát toàn diện kích động ngẫu nhiên phi dừng (non-stationary random excitation) như động đất thực tế dạng El Centro hay sóng thần.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm:

  • Mở rộng thuật toán Taguchi để thiết kế tối ưu hệ hấp thu năng lượng phi tuyến (Nonlinear Energy Sink - NES).
  • Nghiên cứu tương tác động lực học chất lỏng - kết cấu (FSI) 3D cho hệ nhiều bộ TLD bằng mô phỏng số thủy khí động lực học (CFD).
  • Phát triển hệ giảm chấn bán tích cực (Semi-active TMD) kết hợp vật liệu thông minh lưu chất từ nhiệt (Magnetorheological Fluid - MR damper) điều khiển bằng thuật toán mờ hoặc mạng nơ-ron thích nghi.
  • Ứng dụng mô hình hóa trên kết cấu tấm composite và vỏ mỏng của thân vỏ máy bay và cánh tuabin gió ngoài khơi.

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

  • Ảnh hưởng học thuật: Các bài báo trích xuất từ luận án được công bố trên các tạp chí chuyên ngành uy tín trong nước và kỷ yếu hội nghị quốc tế, mở ra hướng ứng dụng phương pháp Taguchi trong động lực học kết cấu tại Việt Nam.
  • Chuyển đổi công nghiệp: Cung cấp hướng dẫn kỹ thuật rõ ràng cho các đơn vị tư vấn thiết kế kết cấu công trình (như Viện Khoa học Công nghệ Xây dựng - IBST, Tổng công ty Tư vấn Xây dựng Việt Nam - VNCC) trong việc tích hợp hệ thống giảm chấn thụ động phân tán vào các tòa nhà chọc trời và tháp truyền hình.
  • Lợi ích kinh tế - xã hội: Giảm giá thành sản xuất và chi phí bảo dưỡng hệ thống kháng chấn tới $30 - 40%$ nhờ việc module hóa các bộ DVA nhỏ gọn, tận dụng các bể chứa nước chữa cháy trên mái công trình làm hệ TLDs tự nhiên, nâng cao độ an toàn sinh mạng con người trước thảm họa thiên tai.

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

  • Nghiên cứu sinh & Giảng viên ngành Cơ học/Xây dựng: Tiếp cận một tài liệu mẫu mực về phương pháp kết hợp giữa cơ học giải tích, phương pháp phần tử hữu hạn/biến phân và kỹ thuật tối ưu hóa thực nghiệm Taguchi.
  • Kỹ sư Thiết kế Động lực học Công trình: Sở hữu bảng số liệu tham số tối ưu và quy trình tính toán sẵn sàng cho các kết cấu dầm, khung chịu tải trọng gió và động đất.
  • Doanh nghiệp Sản xuất Cơ khí Chính xác & Robot: Ứng dụng giải pháp triệt tiêu rung động cục bộ cho băng máy gia công CNC, tay máy công nghiệp tốc độ cao nhằm nâng cao độ nhẵn bề mặt và độ chính xác gia công.
  • Cơ quan Quản lý Quy chuẩn Xây dựng: Cung cấp cơ sở dữ liệu khoa học thực nghiệm để xây dựng các tiêu chuẩn quốc gia về thiết kế kháng chấn bằng thiết bị tiêu tán năng lượng thụ động.

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à nó mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết hai điểm cố định của Den Hartog (1928) và kết hợp với Phương pháp Tuyến tính hóa tương đương để giải quyết trọn vẹn bài toán xác định tham số tối ưu cho hệ nhiều bộ DVA đơn tần số trên hệ chính có cản ($c_s \neq 0$). Luận án chứng minh sự tương đương toán học giữa $N$ bộ DVA giống nhau và một bộ DVA đơn lẻ, đồng thời thiết lập hai tiêu chuẩn mới: Tiêu chuẩn Bình phương tối thiểu (BPTT) và Tiêu chuẩn Đối ngẫu.

2. Sự đổi mới về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?

So với phương pháp giải tích hai điểm cố định thuần túy của Den Hartog (1928) (chỉ áp dụng cho 1 DVA trên hệ $c_s = 0$) và phương pháp thống kê hồi quy thực nghiệm của Fujino và đồng nghiệp (1992) (chỉ tối ưu dải tần đều cho TLD), luận án đã đổi mới bằng cách tích hợp Phương pháp Taguchi với tỷ số tín hiệu trên nhiễu (SNR) và mảng trực giao $L_{18}, L_{27}$. Phương pháp này xử lý xuất sắc bài toán tối ưu hóa đa mục tiêu với hàm mục tiêu có trọng số $f$, co cụm miền tham số lặp tự động, giải quyết được bài toán phân bố tham số không đều cho cả kết cấu liên tục dầm Euler - Bernoulli.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và số liệu chứng minh?

Phát hiện bất ngờ nhất là hiệu năng giảm chấn thực tế của hệ 2 TMD lắp trên dầm côngxôn tại vị trí $x = 47{,}5\text{ cm}$ đạt tới $96{,}84%$ (với $y_0 = 0{,}5\text{ mm}$), vượt qua cả dự báo lý thuyết ($96{,}24%$). Độ lệch giữa mô hình giải tích Ritz - Galerkin và đo đạc vật lý trên phần mềm Pulse LabShop cực kỳ nhỏ, sai số chỉ từ $0{,}12%$ đến $0{,}63%$. Điều này khẳng định tính chính xác tuyệt đối của thuật toán định vị TMD tại các tọa độ nút/bụng của hàm riêng.

Tần số kích động: f = 5.5 Hz | Điểm đo: x = 47.5 cm | y0 = 0.5 mm

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết:

  • Thông số hình học và cơ lý của dầm mẫu ($L, b, h, m_b, E, I, c^{(i)}$).
  • Cấu hình thiết bị đo (Bàn rung TIRA TV50350/LS-120, bộ điều khiển VR-8500, khuếch đại BAA 120, cảm biến Brüel & Kjær, phần mềm Pulse LabShop 10.5).
  • Ma trận phân mức tham số Taguchi và thuật toán cập nhật phân mức qua từng bước lặp, đảm bảo bất kỳ phòng thí nghiệm động lực học nào cũng có thể tái lập kết quả một cách chính xác.

5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Chương trình 10 năm tập trung vào 4 trọng tâm:

  1. Phát triển hệ giảm chấn phi tuyến thích nghi (Nonlinear Adaptive DVAs/NES).
  2. Tích hợp điều khiển bán tích cực (Semi-active MR damper) trên kết cấu không gian 3 chiều.
  3. Ứng dụng trí tuệ nhân tạo (Deep Reinforcement Learning) trong việc điều khiển thời gian thực hệ nhiều bộ giảm chấn.
  4. Thương mại hóa các module TMD/TLD tiêu chuẩn hóa cho công trình dân dụng và công nghiệp tại Đông Nam Á.

Kết luận

Luận án tiến sĩ của tác giả Vũ Đức Phúc là một công trình khoa học công phu, mẫu mực, kết hợp nhuần nhuyễn giữa lý thuyết cơ học giải tích đỉnh cao, phương pháp tính toán tối ưu hiện đại và kỹ thuật thực nghiệm đo kiểm chính xác.

Các đóng góp cốt lõi của công trình được tổng hợp cụ thể:

  1. Mở rộng thành công công thức Den Hartog cho tổ hợp $N$ bộ giảm chấn động lực đơn tần số lắp trên hệ chính có cản thông qua phương pháp tuyến tính hóa tương đương.
  2. Xây dựng mô hình toán học giải tích quy đổi tương đương hoàn chỉnh giữa hệ nhiều bộ giảm chấn chất lỏng (MTLD) và hệ nhiều bộ giảm chấn khối lượng (MTMD).
  3. Đề xuất thuật toán tối ưu hóa tham số mới dựa trên phương pháp thực nghiệm Taguchi và tỷ số SNR, vượt trội hơn các phương pháp cổ điển về khả năng triệt tiêu đỉnh cộng hưởng và mở rộng dải tần dập tắt dao động.
  4. Thiết lập hệ phương trình ma trận chuyển động cho dầm Euler - Bernoulli liên tục có cản lắp nhiều bộ TMD bằng phương pháp biến phân Ritz - Galerkin.
  5. Giải quyết trọn vẹn bài toán tối ưu hóa vị trí không gian lắp đặt các bộ TMD trên dầm dựa vào cấu trúc hình học của các hàm riêng dao động.
  6. Kiểm chứng thực nghiệm thành công trên bàn rung điện động công nghiệp với độ tin cậy cao, chứng minh hiệu quả dập tắt dao động đạt trên $96%$, sai số lý thuyết - thực nghiệm dưới $1%$.

Công trình không chỉ làm phong phú thêm kho tàng lý thuyết Cơ học ứng dụng Việt Nam mà còn mở ra bước tiến quan trọng trong việc ứng dụng công nghệ giảm chấn thông minh, phân tán vào thực tiễn xây dựng và chế tạo máy hiện đại trên phạm vi toàn cầu.