Tổng quan luận án

Dao động cơ học là hiện tượng xuất hiện phổ biến trong vận hành máy móc công nghiệp, thiết bị cơ khí cũng như các công trình kỹ thuật như nhà cao tầng, tháp cầu dây văng, trạm điều khiển không lưu và hệ thống giàn khoan biển. Các rung động quá mức dưới tác dụng của tải trọng ngoại lực, sóng hoặc gió thường gây suy giảm độ chính xác gia công, làm giảm tuổi thọ làm việc và có nguy cơ dẫn tới phá hủy kết cấu. Trong các giải pháp kiểm soát dao động thụ động, bộ giảm chấn động lực (Dynamic Vibration Absorber - DVA) được ứng dụng rộng rãi nhờ ưu điểm kết cấu đơn giản, vận hành độc lập và không đòi hỏi nguồn năng lượng bên ngoài. Tuy nhiên, việc áp dụng một bộ giảm chấn đơn lẻ (Single DVA) bộc lộ các hạn chế: kích thước và khối lượng bộ giảm chấn quá lớn đối với các kết cấu quy mô lớn, gây khó khăn cho việc chế tạo lò xo, phần tử cản cũng như công tác vận chuyển, lắp đặt; đồng thời, nếu bộ giảm chấn đơn lẻ gặp sự cố kỹ thuật thì toàn bộ khả năng kháng chấn của hệ thống sẽ bị mất. Do đó, hướng nghiên cứu kết hợp nhiều bộ giảm chấn động lực (DVAs), bao gồm hệ giảm chấn khối lượng – cản nhớt (Tuned Mass Dampers - TMDs) và hệ giảm chấn chất lỏng (Tuned Liquid Dampers - TLDs), mang lại khả năng phân tán kích thước, dễ chế tạo hàng loạt và duy trì độ tin cậy vận hành.

Khoảng trống nghiên cứu được luận án xác định bao gồm:

  • Thiếu các biểu thức giải tích tường minh xác định tham số tối ưu cho hệ nhiều bộ giảm chấn động lực đơn tần số (Single Dynamic Vibration Absorber - SDVA) lắp trên hệ chính một bậc tự do cho cả trường hợp không cản và có cản.
  • Chưa có nghiên cứu ứng dụng phương pháp thiết kế thực nghiệm số (phương pháp Taguchi) để giải quyết bài toán tối ưu hóa đa tham số cho hệ nhiều bộ giảm chấn động lực đa tần số (Multiple Dynamic Vibration Absorber - MDVA) lắp trên hệ chính có cản.
  • Thiếu mô hình động lực học tổng quát cho kết cấu dầm Euler - Bernoulli có kể đến cả cản trong và cản ngoài lắp nhiều bộ giảm chấn động lực, cũng như chưa có phân tích định lượng về việc tối ưu hóa vị trí lắp đặt các bộ giảm chấn dựa trên dạng dao động riêng của dầm.

Mục đích nghiên cứu chủ yếu của luận án là kết hợp nhiều bộ giảm chấn động lực để điều khiển giảm dao động cho hệ chính có cản, với nhiệm vụ trọng tâm là thiết kế tối ưu các tham số của các bộ giảm chấn động lực sao cho dao động của hệ chính đạt cực tiểu trong miền tần số cộng hưởng.

Nội dung nghiên cứu của luận án được triển khai qua 3 nhiệm vụ chính:

  1. Thiết lập hệ phương trình vi phân dao động và hàm đáp ứng tần số của hệ kết cấu và các bộ giảm chấn động lực.
  2. Tìm các tham số tối ưu của các bộ giảm chấn động lực để điều khiển giảm dao động có hại cho hệ ở tần số cộng hưởng và trong miền tần số cộng hưởng.
  3. Xác nhận các kết quả và hiệu quả điều khiển giảm dao động của các bộ giảm chấn thông qua các mô phỏng số về đáp ứng tần số và đáp ứng thời gian của hệ; so sánh với các kết quả đã biết để kiểm chứng độ tin cậy.

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

  • Đối tượng nghiên cứu: Hệ dao động có cản được mô hình hóa thành kết cấu một bậc tự do và kết cấu dạng dầm Euler - Bernoulli chịu kích động điều hòa có lắp nhiều bộ giảm chấn động lực (TMDs và TLDs).
  • Phạm vi nghiên cứu: Phân tích dao động tuyến tính cưỡng bức trong miền tần số xung quanh tần số cộng hưởng; quy đổi tương đương hệ nhiều bộ giảm chấn chất lỏng (TLDs) về hệ giảm chấn khối lượng (TMDs); khảo sát kết cấu dầm đàn hồi với hai điều kiện biên điển hình là dầm hai đầu bản lề và dầm một đầu ngàm một đầu tự do (dầm công xôn).

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

Tổng quan lịch sử nghiên cứu bộ giảm chấn động lực ghi nhận ý tưởng đầu tiên của P. Watt (1883) về giảm rung lắc tàu thủy và sáng chế DVA không cản của H. Frahm (1909). Cơ sở lý thuyết giải tích cho DVA có cản nhớt lắp trên hệ chính không cản được J. P. Den Hartog (1928) xây dựng dựa trên phương pháp hai điểm cố định, sau đó được H. Hahnkamm (1932) và J. E. Brock (1946) hoàn thiện chứng minh toán học. H. Yamaguchi (1988) tiếp cận bài toán bằng tiêu chuẩn cực đại độ tắt lôga; T. Asami và O. Nishihara (2002) tìm ra nghiệm giải tích đại số chính xác.

Đối với hệ chính có cản lắp DVA có cản:

  • T. Ioi và K. Ikeda (1978) sử dụng phép lặp Newton xác định tỷ lệ tần số và tỷ số cản tối ưu thông qua hàm mục tiêu cân bằng hai đỉnh cộng hưởng.
  • S. E. Randall và đồng nghiệp (1981) đề xuất tiêu chuẩn min-max tìm cực tiểu trên bề mặt đáp ứng số.
  • T. Igusa và A. Der Kiureghian (1985), Y. Fujino và M. Abe (1993) áp dụng phương pháp nhiễu loạn để thu được các biểu thức giải tích xấp xỉ.
  • T. Asami và đồng nghiệp (2002) đưa ra lời giải giải tích dạng đóng theo tiêu chuẩn tối ưu $H_2$ và $H_\infty$.

Đối với hệ nhiều bộ giảm chấn động lực (DVAs/MTMD):

  • T. Igusa và K. Xu (1992, 1994) đề xuất phân bố tần số của các bộ DVAs xung quanh tần số riêng của kết cấu và chứng minh hệ DVAs giảm độ nhạy cảm với sai lệch điều chỉnh tham số.
  • Seung-Yong Oka và đồng nghiệp (2006) sử dụng kỹ thuật nội suy đường cong phi tuyến.
  • Lei Zuo và S. A. Nayfeh (2004, 2005) áp dụng phương pháp Gradient và điều khiển không gian trạng thái; Nam Hoang và P. Warnitchai (2007) sử dụng thuật toán phi tuyến trên cơ sở gradient và hàm Lyapunov.
  • M. Mohebbi và đồng nghiệp (2013) sử dụng thuật giải di truyền (GAs) tối ưu hệ 15 bộ DVA cho nhà cao tầng; Sung-Yong Kim và Cheol-Ho Lee (2018) thiết lập bản đồ đường bao giới hạn không gian hoạt động của các DVAs.
  • Hệ giảm chấn chất lỏng (TLDs/MTLD) dựa trên lý thuyết sóng nước nông được nghiên cứu bởi V. J. Modi, F. Welt, A. Kareem, Y. Fujino và đồng nghiệp từ cuối thế kỷ 20.

Đối với kết cấu dầm Euler - Bernoulli lắp DVAs:

  • J. C. Snowdon (1975) khảo sát dầm công xôn không cản lắp 2 bộ DVA; R. G. Jacquot (1978) dùng phương pháp giả định mode; M. Noori và E. Farshidianfar (2014) giải bài toán theo tiêu chuẩn $H_\infty$ và $H_2$.
  • H. N. Ozguven và B. Candir (1986) tối ưu 2 DVA lắp trên dầm có cản theo tiêu chuẩn min-max; M. Latas (2015, 2017) khảo sát hệ giảm chấn phân bố liên tục.

Tình hình nghiên cứu trong nước:

  • Đề tài cấp nhà nước KC05-30 (2005) do Nguyễn Đông Anh chủ trì đã nghiên cứu thiết kế, chế tạo TMD ứng dụng cho cầu và phát triển phần mềm phân tích kết cấu.
  • Nguyễn Đông Anh và Lã Đức Việt (2006) tổng kết các phương pháp tính toán tối ưu DVA trong chuyên khảo chuyên ngành.
  • Nguyễn Huy Thế (2010) khảo sát giải tích lý thuyết điểm cố định; Lã Đức Việt và Phan Thị Trà My (2013) tiếp cận theo hướng dòng năng lượng cực đại.
  • N. D. Anh và đồng nghiệp (2013, 2014, 2015, 2017) phát triển phương pháp tuyến tính hóa tương đương cho hệ có cản dựa trên tiêu chuẩn bình phương tối thiểu, tiêu chuẩn đối ngẫu và tiêu chuẩn đối ngẫu có trọng số.
  • Nguyễn Chỉ Sáng và đồng nghiệp (2012, 2014) nghiên cứu DVA cho hệ nhiều bậc tự do; Phạm Nhân Hòa và Chu Quốc Thắng (2009) dùng thuật toán Newmark phân tích cản chất lỏng nhớt; Nguyễn Anh Tuấn (2017) ứng dụng TMD giảm rung bảo vệ cơ thể người lao động.
  • Nguyễn Đức Thị Thu Định và đồng nghiệp (2010, 2012, 2016) nghiên cứu mô hình quy đổi TLD thành TMD áp dụng cho tháp cầu Bãi Cháy dưới dạng hệ một bậc tự do.

Vị trí của luận án: Kế thừa các nền tảng lý thuyết trên, luận án giải quyết 4 vấn đề còn khuyết:

  1. Mở rộng công thức giải tích Den Hartog cho hệ nhiều bộ DVA đơn tần số trên hệ chính không cản và áp dụng phương pháp tuyến tính hóa tương đương (tiêu chuẩn bình phương tối thiểu, tiêu chuẩn đối ngẫu) cho hệ chính có cản.
  2. Ứng dụng phương pháp thiết kế thực nghiệm Taguchi để tối ưu hóa đồng thời nhiều tham số cho hệ DVAs đa tần số và hệ TLDs quy đổi.
  3. Thiết lập mô hình động lực học hoàn chỉnh cho dầm Euler - Bernoulli có kể đến đồng thời cản trong và cản ngoài lắp nhiều bộ DVA.
  4. Tối ưu hóa tham số và vị trí lắp đặt các bộ DVAs trên dầm đàn hồi dựa trên đặc trưng hàm riêng dao động.

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

Luận án phối hợp đồng bộ giữa các cơ sở lý thuyết cơ học giải tích, phương pháp giải tích tương đương, phương pháp quy hoạch thực nghiệm số Taguchi và mô phỏng số học:

Khung lý thuyết / Phương pháp Nội dung và công cụ toán học áp dụng trong luận án
Phương trình Lagrange loại 2 Thiết lập hệ phương trình vi phân chuyển động liên kết dạng ma trận $[M]{\ddot{x}} + [D]{\dot{x}} + [K]{x} = {f(t)}$ cho hệ dao động rời rạc và hệ dầm liên tục rời rạc hóa.
Phương pháp biên độ phức Chuyển đổi hệ phương trình vi phân sang miền tần số với kích động phức $F_0 e^{i\Omega t}$ để suy ra hàm đáp ứng tần số phức $\hat{x}_s(\Omega)$ và hàm độ khuyếch đại $H(\Omega)$.
Lý thuyết hai điểm cố định mở rộng Thiết lập điều kiện cân bằng biên độ tại hai điểm cố định $S, T$ độc lập với tỷ số cản nhằm giải tích tường minh tỷ lệ tần số tối ưu $\gamma_{opt}$ và tỷ số cản tối ưu $\xi_{opt}$ cho hệ nhiều bộ DVA đơn tần số.
Phương pháp tuyến tính hóa tương đương Thay thế hệ chính có cản bằng hệ chính không cản tương đương dựa trên tiêu chuẩn bình phương sai lệch tối thiểu (BPTT) và tiêu chuẩn đối ngẫu.
Lý thuyết sóng nước nông & Quy đổi TLD - TMD Quy đổi mô hình dao động sóng chất lỏng trong thùng hở chữ nhật về hệ khối lượng - lò xo - cản nhớt tương đương (TMD) để tích hợp vào bài toán dao động kết cấu.
Phương pháp Ritz - Galerkin Khai triển độ võng dầm liên tục $w(x,t) = \sum_{k=1}^m X_k(x) q_k(t)$ theo hệ hàm riêng trực giao $X_k(x)$ để rời rạc hóa phương trình đạo hàm riêng của dầm Euler - Bernoulli có cản.
Phương pháp thiết kế thực nghiệm Taguchi Sử dụng các bảng mảng trực giao $L_{27}(3^{13})$ và $L_{18}(2^1 \times 3^7)$, đánh giá tỷ số tín hiệu trên nhiễu $SNR$ (đặc trưng chất lượng "càng nhỏ càng tốt" - smaller-the-better), phân tích trung bình $ANOM$, phương sai $ANOVA$ và thuật toán thu hẹp mức tham số để tìm nghiệm tối ưu đa biến.
Mô phỏng số và phân tích đáp ứng Lập trình tính toán trên phần mềm Matlab, khảo sát đáp ứng tần số $H(\Omega)$ và đáp ứng thời gian $x(t)$ dưới kích động điều hòa cưỡng bức.

Nguồn dữ liệu tính toán trong luận án được chuẩn hóa từ các thông số kết cấu cơ học tiêu chuẩn: hệ một bậc tự do có khối lượng $m_s = 250\text{ kg}$, độ cứng $k_s = 1.500.000\text{ N/m}$, lực kích động $F_0 = 250\text{ N}$, tần số riêng $\omega_s = 77,46\text{ rad/s}$; kết cấu dầm Euler - Bernoulli có thông số chiều dài, mặt cắt, mô đun đàn hồi, hệ số cản trong $c^{(i)}$ và cản ngoài $c^{(e)}$ xác định.

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

Chương 1: Tổng quan nghiên cứu điều khiển dao động bằng bộ giảm chấn động lực

Chương 1 hệ thống hóa các nguyên lý cơ bản của DVA từ mô hình không cản ban đầu: $$m_s \ddot{x}s + (k_s + k_a) x_s - k_a x_a = F_0 \sin(\Omega t)$$ $$m_a \ddot{x}a - k_a x_s + k_a x_a = 0$$ Tại tần số kích động $\Omega$, biên độ dao động hệ chính triệt tiêu hoàn toàn khi thỏa mãn điều kiện $k_a = m_a \Omega^2$. Từ đó, chương phân tích sự phát triển sang mô hình có cản nhớt của Den Hartog với các tham số tối ưu cho hệ không cản lắp 1 DVA: $$\gamma{opt} = \frac{1}{1 + \mu}, \quad \xi{opt} = \sqrt{\frac{3\mu}{8(1 + \mu)}}$$ Chương tổng thuật 131 tài liệu tham khảo trong và ngoài nước, phân định rõ các ưu nhược điểm của các phương pháp số, phương pháp giải tích, phương pháp năng lượng và chỉ ra 4 nhiệm vụ khoa học trọng tâm của luận án.

Chương 2: Điều khiển tối ưu dao động của hệ chính có cản bằng nhiều bộ giảm chấn động lực đơn tần số dựa trên phương pháp của Den-Hartog

Chương 2 thiết lập mô hình toán học cho hệ một bậc tự do lắp $n_a$ bộ giảm chấn động lực đơn tần số (SDVA). Sử dụng phương pháp Lagrange loại 2, hệ phương trình vi phân chuyển động được viết dưới dạng ma trận cỡ $(n_a + 1) \times (n_a + 1)$: $$[M]{\ddot{x}} + [D]{\dot{x}} + [K]{x} = {f(t)}$$

Bằng phương pháp biên độ phức, hàm đáp ứng tần số của hệ được thiết lập: $$H(\Omega) = \frac{1}{\sqrt{a^2 + b^2}}$$ với các hệ số $a, b$ phụ thuộc vào tỷ số tần số $\lambda = \Omega/\omega_s$, tỷ lệ tần số $\gamma = \omega_a/\omega_s$, tỷ số cản $\xi_a$ và tỷ lệ khối lượng $\mu_a = m_a/m_s$.

Đối với hệ chính không cản ($\xi_s = 0$), luận án mở rộng thành công lý thuyết hai điểm cố định $S, T$, chứng minh tung độ tại hai điểm này không phụ thuộc vào hệ số cản và tìm ra hệ công thức giải tích tường minh: $$\gamma_{opt} = \frac{1}{1 + n_a \mu_a}$$ $$\xi_{opt} = \sqrt{\frac{3 n_a \mu_a}{8(1 + n_a \mu_a)}}$$ Độ cứng và hệ số cản tối ưu của mỗi bộ giảm chấn trong hệ $n_a$ bộ được xác định: $$k_a = \frac{k_s \mu_a}{(1 + n_a \mu_a)^2}, \quad c_a = 2 m_s \omega_s \mu_a \sqrt{\frac{3 n_a \mu_a}{8(1 + n_a \mu_a)^3}}$$ Biên độ đáp ứng cực đại tại hai đỉnh cộng hưởng cân bằng đạt giá trị: $$H_S = H_T = \sqrt{\frac{n_a \mu_a + 2}{n_a \mu_a}}$$

Kết quả khảo sát số học đối với hệ chính ($m_s = 250\text{ kg}$, $k_s = 1.500.000\text{ N/m}$, $\omega_s = 77,46\text{ rad/s}$) được tổng hợp trong các bảng số liệu sau:

Cấu hình khảo sát Tổng tỷ lệ khối lượng $m_a$ mỗi bộ (kg) $k_a$ mỗi bộ (N/m) $c_a$ mỗi bộ (Ns/m) Biên độ $x_s$ (mm) $H_{max}$
1 DVA $\mu = 1%$ 2,50 14704,44 23,3657 2,06 14,18
2 DVA $\sum \mu = 1%$ 1,25 7352,22 11,6829 2,06 14,18
5 DVA $\sum \mu = 1%$ 0,50 2940,88 4,6731 2,06 14,18

Số liệu trên chứng minh tính tương đương hoàn toàn giữa 1 bộ DVA lớn và hệ $n_a$ bộ DVA nhỏ khi có cùng tổng khối lượng, tổng độ cứng và tổng hệ số cản.

Khi đánh giá độ tin cậy vận hành của hệ trong trường hợp mỗi bộ có $\mu_a = 1%$ ($m_a = 2,5\text{ kg}$):

  • Lắp 5 DVA: $x_s = 0,41\text{ mm}$, $H_{max} = 6,41$.
  • Khi hỏng giảm còn 2 DVA: $x_s = 1,03\text{ mm}$, $H_{max} = 10,05$.
  • Khi hỏng giảm còn 1 DVA: $x_s = 2,06\text{ mm}$, $H_{max} = 14,18$. Đường cong đáp ứng vẫn giữ tính đối xứng tại hai đỉnh, khẳng định khi một số bộ bị hỏng thì các bộ còn lại vẫn duy trì hiệu quả kháng chấn.

Đối với hệ chính có cản ($\xi_s \neq 0$), luận án áp dụng phương pháp tuyến tính hóa tương đương, chuyển đổi hệ có cản về hệ không cản tương đương có tần số $\omega_e^2 = \omega_s^2 + \gamma$ thông qua phiếm hàm sai lệch: $$S = \int_D (2\xi_s \omega_s \dot{x}_s - \gamma x_s)^2 dt \to \min$$ Từ tiêu chuẩn bình phương tối thiểu (BPTT) và tiêu chuẩn đối ngẫu, các công thức tối ưu cho hệ có cản được thiết lập đầy đủ.

Chương 3: Điều khiển tối ưu dao động của hệ chính có cản bằng nhiều bộ giảm chấn động lực đa tần số dựa trên phương pháp Taguchi

Chương 3 giải quyết bài toán tối ưu cho hệ nhiều bộ giảm chấn động lực đa tần số (MDVA/MTMD) và hệ giảm chấn chất lỏng (TLDs/MTLD):

  1. Quy đổi tương đương TLD thành TMD: Dựa trên lý thuyết sóng nước nông trong thùng chữ nhật có chiều dài $l$, chiều rộng $b$, chiều cao mực nước $h_0$, khối lượng chất lỏng $m$, tần số cơ bản của khối chất lỏng dao động được xác định: $$f_{TLD} = \frac{1}{2\pi}\sqrt{\frac{\pi g}{l} \tanh\left(\frac{\pi h_0}{l}\right)}$$ Khối lượng chất lỏng tham gia dao động và hệ số cản nhớt tương đương được quy đổi chính xác sang mô hình TMD tương đương.
  2. Thuật toán tối ưu hóa đa thông số Taguchi: Thiết kế mảng trực giao $L_{27}(3^{13})$ để tối ưu đồng thời 10 tham số điều khiển ($k_1, c_1, k_2, c_2, \dots, k_5, c_5$) của hệ 5 bộ giảm chấn.
  • Tỷ số tín hiệu trên nhiễu được tính theo tiêu chuẩn "càng nhỏ càng tốt": $$SNR = -10 \log_{10}\left(\frac{1}{n}\sum_{i=1}^n y_i^2\right)$$
  • Hàm mục tiêu khảo sát gồm hàm đáp ứng biên độ tần số cực đại và hàm mục tiêu có trọng số $f$.
  • Tiến hành phân tích giá trị trung bình ($ANOM$) và phân tích phương sai ($ANOVA$) sau các bước lặp để xác định mức ảnh hưởng của từng tham số và tự động thu hẹp miền phân mức mới cho đến khi hội tụ.
  1. So sánh kết quả:
  • Hệ 5 DVA đa tần số cho vùng tần số hiệu quả rộng hơn và giảm biên độ cộng hưởng đồng đều hơn so với hệ đơn tần số.
  • Kết quả tối ưu theo thuật toán Taguchi đạt độ chính xác cao khi so sánh đối chứng với phương pháp giải tích thực nghiệm của Fujino và cộng sự.

Chương 4: Điều khiển tối ưu dao động uốn cưỡng bức của dầm Euler - Bernoulli có cản bằng nhiều bộ giảm chấn động lực

Chương 4 mở rộng bài toán sang hệ liên tục dầm Euler - Bernoulli chiều dài $L$ có kể đến đồng thời hệ số cản trong kết cấu $c^{(i)}$ và hệ số cản ngoài môi trường $c^{(e)}$: $$\rho A \frac{\partial^2 w(x,t)}{\partial t^2} + c^{(e)} \frac{\partial w(x,t)}{\partial t} + c^{(i)} I \frac{\partial^5 w(x,t)}{\partial x^4 \partial t} + E I \frac{\partial^4 w(x,t)}{\partial x^4} = p(x,t) - \sum_{j=1}^{n_a} F_j(t) \delta(x - x_j)$$ trong đó $F_j(t) = d_j [\dot{w}(x_j,t) - \dot{u}_j(t)] + k_j [w(x_j,t) - u_j(t)]$ là lực tương tác từ bộ giảm chấn thứ $j$ đặt tại vị trí $x_j$.

Áp dụng phương pháp Ritz - Galerkin, chuyển vị dầm được biểu diễn qua chuỗi hàm riêng: $$w(x,t) = \sum_{k=1}^m X_k(x) q_k(t)$$ Hệ phương trình vi phân chuyển động của dầm liên kết với các TMDs được đưa về hệ phương trình ma trận tuyến tính và thiết lập hàm đáp ứng tần số tổng quát.

Luận án ứng dụng phương pháp Taguchi với mảng trực giao $L_{18}(2^1 \times 3^7)$ để tối ưu hóa các thông số độ cứng $k_j$ và hệ số cản $d_j$ của 3 bộ TMD:

  • Dầm hai đầu bản lề: Hàm riêng dạng sin $X_k(x) = \sin(k\pi x / L)$. Kết quả phân tích cho thấy lắp các bộ TMD tập trung tại vị trí giữa dầm ($x = L/2$ - nơi có biên độ hàm riêng thứ nhất cực đại) mang lại hiệu quả giảm chấn uốn cao nhất.
  • Dầm một đầu ngàm một đầu tự do: Hàm riêng kết hợp hàm lượng giác và hàm hyperbol. Khảo sát vị trí lắp đặt chỉ ra rằng đặt các bộ giảm chấn tại đầu tự do ($x = L$) triệt tiêu biên độ dao động uốn tại tần số riêng thứ nhất hiệu quả nhất.

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

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

  1. Mở rộng thành công phương pháp hai điểm cố định của Den Hartog cho hệ nhiều bộ giảm chấn động lực đơn tần số (SDVA) lắp trên hệ chính một bậc tự do không cản; thiết lập các công thức giải tích tường minh xác định tỷ lệ tần số tối ưu $\gamma_{opt}$, tỷ số cản tối ưu $\xi_{opt}$, độ cứng tối ưu $k_a$ và hệ số cản tối ưu $c_a$.
  2. Ứng 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 sai lệch tối thiểu và tiêu chuẩn đối ngẫu để giải quyết bài toán thiết kế tối ưu hệ nhiều bộ DVA đơn tần số lắp trên hệ chính có cản.
  3. Đề xuất thuật toán tối ưu hóa đa tham số dựa trên phương pháp quy hoạch thực nghiệm Taguchi (mảng trực giao $L_{27}$ và $L_{18}$, tỷ số $SNR$, phân tích $ANOM/ANOVA$), cho phép giải quyết bài toán tối ưu nhiều biến phức tạp cho hệ DVAs đa tần số và hệ dầm liên tục có cản.
  4. Xây dựng mô hình động lực học tổng quát dầm Euler - Bernoulli có đồng thời cản trong và cản ngoài lắp nhiều bộ DVA bằng phương pháp Ritz - Galerkin; xác lập nguyên lý bố trí vị trí lắp đặt tối ưu của các bộ DVA dựa trên cực đại của dạng mode dao động riêng.

Đóng góp về mặt thực tiễn và giải pháp đề xuất

  1. Chứng minh tính tương đương động lực học giữa một bộ DVA kích thước lớn và nhiều bộ DVA kích thước nhỏ có cùng tổng tham số, mở ra giải pháp phân tán kết cấu kháng chấn, giảm tải trọng tập trung, thuận lợi cho chế tạo, vận chuyển và thi công lắp đặt trên công trình thực tế.
  2. Nâng cao độ an toàn vận hành: việc phân chia thành nhiều bộ DVA giúp hệ thống duy trì hiệu năng dập tắt dao động khi một hoặc một vài bộ giảm chấn gặp sự cố hư hỏng.
  3. Cung cấp quy trình tính toán quy đổi tương đương từ hệ giảm chấn chất lỏng (TLDs) sang hệ khối lượng (TMDs) phục vụ thiết kế kháng chấn cho công trình tháp, cầu dây văng và nhà cao tầng.

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

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

  • Các mô hình dao động kết cấu và phần tử cản nhớt mới dừng lại ở phạm vi lý thuyết tuyến tính; chưa xét tới tính phi tuyến hình học, phi tuyến vật liệu hoặc cản phi tuyến khi kết cấu chịu biến dạng lớn.
  • Dạng tải trọng kích động chủ yếu là lực kích động điều hòa xác định trong miền tần số cộng hưởng; chưa mở rộng khảo sát các dạng kích động ngẫu nhiên băng rộng, kích động động đất dạng ồn trắng lọc (Kanai - Tajimi) hoặc tải trọng gió bão phi dừng.
  • Vị trí lắp đặt bộ giảm chấn trên dầm được khảo sát theo phương pháp so sánh các điểm rời rạc định trước dựa trên mode riêng, chưa xây dựng bài toán tối ưu đồng thời cả tọa độ không gian liên tục $x_j$ và các tham số cơ học ($k_j, c_j$) trong cùng một hàm mục tiêu giải tích.

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

  • Mở rộng thuật toán Taguchi và phương pháp tuyến tính hóa tương đương cho hệ kết cấu nhiều bậc tự do phi tuyến, kết cấu tấm, vỏ đàn hồi.
  • Nghiên cứu đáp ứng động lực học của dầm và kết cấu dưới tác dụng của tải trọng ngẫu nhiên và tải trọng đoàn xe di động.
  • Nghiên cứu thực nghiệm mở rộng mô hình dầm liên tục lắp hệ giảm chấn chất lỏng đa tần số (MTLD) trong điều kiện phòng thí nghiệm quy mô lớn.

Giá trị tham khảo

Luận án là tài liệu tham khảo hữu ích cho các nhóm đối tượng chuyên môn:

  • Nghiên cứu sinh, học viên cao học và giảng viên ngành Cơ học kỹ thuật, Kỹ thuật cơ khí, Kỹ thuật xây dựng công trình, Cầu hầm: Tham khảo phương pháp mở rộng lý thuyết Den Hartog (Chương 2), kỹ thuật rời rạc hóa dầm Euler - Bernoulli có cản bằng phương pháp Ritz - Galerkin (Chương 4) và thuật toán quy hoạch thực nghiệm số Taguchi áp dụng trong cơ học dao động (Chương 3 và 4).
  • Kỹ sư thiết kế kết cấu công trình và máy công nghiệp: Ứng dụng trực tiếp hệ công thức giải tích tính toán nhanh thông số lò xo, cản nhớt cho hệ giảm chấn phân tán; tham khảo quy trình quy đổi thiết kế bể giảm chấn chất lỏng (TLD) cho các công trình cao tầng và tháp cầu.

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

1. Công thức giải tích mở rộng Den Hartog cho hệ $n_a$ bộ DVA đơn tần số lắp trên hệ không cản có dạng như thế nào?
Tỷ lệ tần số tối ưu và tỷ số cản tối ưu của mỗi bộ giảm chấn được xác định qua biểu thức: $$\gamma_{opt} = \frac{1}{1 + n_a \mu_a}, \quad \xi_{opt} = \sqrt{\frac{3 n_a \mu_a}{8(1 + n_a \mu_a)}}$$ trong đó $\mu_a = m_a/m_s$ là tỷ lệ khối lượng của một bộ DVA so với hệ chính.

2. Việc thay thế 1 bộ DVA lớn bằng $n_a$ bộ DVA nhỏ có làm suy giảm hiệu quả dập tắt dao động không?
Không. Kết quả phân tích toán học và mô phỏng số tại Chương 2 chứng minh rằng nếu tổng khối lượng, tổng độ cứng và tổng hệ số cản của $n_a$ bộ DVA bằng đúng các thông số tương ứng của 1 bộ DVA đơn lẻ ($\sum m_j = m_a, \sum k_j = k_a, \sum c_j = c_a$) thì đáp ứng tần số $H(\Omega)$ và đáp ứng thời gian $x(t)$ của hai hệ thống là hoàn toàn đồng nhất ($H_{max} = 14,18$, $x_s = 2,06\text{ mm}$ tại $\mu = 1%$).

3. Khi một hoặc một số bộ DVA trong hệ đơn tần số bị hỏng hóc, hệ thống hoạt động ra sao?
Hệ thống vẫn duy trì khả năng kháng chấn ổn định. Khi mỗi bộ DVA có tỷ lệ khối lượng $\mu_a = 1%$, nếu số bộ làm việc giảm từ 5 bộ xuống 2 bộ và 1 bộ, biên độ dao động hệ chính tăng tương ứng từ $0,41\text{ mm}$ lên $1,03\text{ mm}$ và $2,06\text{ mm}$, đường cong đáp ứng tần số vẫn duy trì tính đối xứng cân bằng tại hai đỉnh cộng hưởng.

4. Phương pháp Taguchi được triển khai như thế nào để tối ưu hóa hệ 5 bộ DVA đa tần số?
Luận án sử dụng bảng mảng trực giao $L_{27}(3^{13})$ với 10 tham số điều khiển ($k_1, c_1, \dots, k_5, c_5$), mỗi tham số chia 3 mức giá trị. Thuật toán tính toán tỷ số $SNR$ theo tiêu chuẩn "càng nhỏ càng tốt" cho hàm đáp ứng biên độ tần số, thực hiện phân tích $ANOM$ và $ANOVA$ để đánh giá mức độ ảnh hưởng của từng tham số, sau đó tiến hành thu hẹp miền phân mức qua các vòng lặp để xác định bộ tham số tối ưu toàn cục.

5. Vị trí lắp đặt bộ giảm chấn nào mang lại hiệu quả kháng chấn cao nhất trên kết cấu dầm?
Hiệu quả kháng chấn của các bộ giảm chấn phụ thuộc trực tiếp vào giá trị của hàm riêng $X_k(x)$ tại tọa độ lắp đặt. Đối với dầm hai đầu bản lề, vị trí tối ưu nằm tại giữa dầm ($x = L/2$). Đối với dầm một đầu ngàm một đầu tự do (dầm công xôn), vị trí tối ưu nằm tại đầu tự do ($x = L$).

Kết luận

Luận án tiến sĩ của tác giả Vũ Đức Phúc đã giải quyết trọn vẹn bài toán điều khiển tối ưu dao động bằng cách kết hợp nhiều bộ giảm chấn động lực (SDVA, MDVA, TMDs, TLDs) cho hệ một bậc tự do và dầm Euler - Bernoulli có cản. Công trình đã đóng góp hệ thống công thức giải tích mở rộng Den Hartog tường minh cho hệ nhiều DVA đơn tần số, đồng thời xây dựng quy trình tối ưu hóa đa thông số tin cậy dựa trên phương pháp thiết kế thực nghiệm Taguchi. Các kết quả mô phỏng số đã kiểm chứng tính tương đương động lực học của việc phân tán giảm chấn và xác lập cơ sở khoa học cho việc ứng dụng các hệ thống giảm chấn đa phần tử an toàn, gọn nhẹ trong kỹ thuật cơ khí và xây dựng công trình.