Tổng quan về luận án

Dao động tần số thấp (<25 Hz, đặc biệt là <5 Hz) là nguyên nhân cốt lõi dẫn đến các tổn thương cơ xương khớp, thoái hóa cột sống ở người vận hành và gây suy giảm nghiêm trọng độ bền mỏi của thiết bị cơ khí chính xác. Các hệ thống cách ly thụ động truyền thống—gồm một phần tử đàn hồi tuyến tính kết hợp giảm chấn—luôn đối mặt với mâu thuẫn cơ bản: để mở rộng vùng cách ly về miền tần số cực thấp, độ cứng của hệ bắt buộc phải giảm; tuy nhiên, độ cứng thấp lại làm suy giảm nghiêm trọng khả năng chịu tải tĩnh và gây ra biến dạng tĩnh ($H_0$) không thể chấp nhận. Nhằm giải quyết triệt để rào cản này, công trình tiến sĩ chuyên ngành Kỹ thuật Cơ khí (mã ngành 9520103) của nghiên cứu sinh Võ Ngọc Yến Phương tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh đã tiên phong phát triển: "Hệ thống cách ly dao động có độ cứng gần bằng không" (Quasi-Zero Stiffness Adaptive Vibration Isolation Model - QSAVIM).

+-----------------------------------------------------------------------------+
|                          HỆ THỐNG QSAVIM TOÀN DIỆN                          |
+-------------------------------------+---------------------------------------+
| Cơ cấu đỡ tải (LBM - Độ cứng dương) | Cơ cấu hiệu chỉnh (SCM - Độ cứng âm)  |
| - Góc nêm alpha = 37 độ             | - Cam bán nguyệt (R=60mm, r=20mm)     |
| - Lò xo khí cao su / Xylanh khí     | - Khử độ cứng động học tại u=0 (DSEP) |
+-------------------------------------+---------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------+
|    TỔNG HỢP: Độ cứng tĩnh cao (đỡ tải lớn) + Độ cứng động học tiệm cận 0    |
|    Dải cách ly mở rộng: triệt tiêu dao động kích động từ omega > 18.5 rad/s |
+-----------------------------------------------------------------------------+

Nghiên cứu tập trung vào khoảng trống học thuật thực chứng (empirical research gap): các cơ cấu Quasi-Zero Stiffness (QZS) kim loại hoặc nam châm vĩnh cửu trước đây (như Carrella et al., 2007; Liu et al., 2013) chỉ đạt hiệu năng tối ưu tại một mức tải trọng tĩnh cố định; khi tải trọng thay đổi, độ cứng động học tại vị trí cân bằng tĩnh mong muốn (Desirable Static Equilibrium Position - DSEP) bị lệch pha, làm mất đặc tính QZS. Luận án đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • RQ1: Làm thế nào để mô hình hóa chính xác tương tác phi tuyến đa vật lý (khí nén, đàn nhớt phân số, ma sát trượt) của phần tử phi kim loại trong cơ cấu QZS?
  • RQ2: Quy luật điều chỉnh tỷ số áp suất $\mu = P_{wh2}/P_{wh1}$ và thể tích buồng phụ $V_{ac}$ tác động ra sao đến độ ổn định, hiện tượng rẽ nhánh (bifurcation) và hiện tượng nhảy tần số (frequency jump)?
  • RQ3: Mô hình QZS sử dụng phần tử khí nén có thể mở rộng băng thông cách ly ở tần số thấp hơn bao nhiêu so với mô hình truyền thống tương đương (ETVIM)?
  • H1: Việc kết hợp cơ cấu đỡ tải (LBM) độ cứng dương với cơ cấu hiệu chỉnh (SCM) độ cứng âm thông qua lò xo khí/xylanh khí sẽ triệt tiêu độ cứng động học tại DSEP ($K_{DSEP} \approx 0$) mà không làm tăng biến dạng tĩnh.
  • H2: Tích hợp mô hình ma sát Berg và đàn nhớt phân số Kelvin-Voigt cho phép dự báo chính xác hiện tượng trễ cơ học với sai số cực tiểu thông qua thuật toán tối ưu hóa di truyền (GA).
  • H3: Hệ thống QSAVIM mở rộng dải tần số cách ly hiệu dụng về miền dưới 20 rad/s ($<3.2$ Hz), giảm tần số cộng hưởng ít nhất 50% so với hệ thống ETVIM.

Khung lý thuyết của nghiên cứu được xây dựng trên sự tích hợp của Nhiệt động lực học khí nén biến thiên, Cơ học vật liệu đàn nhớt phân số, Lý thuyết phân tích phi tuyến (Normal Form, Multi-Scale) và Lý thuyết động lực học hỗn loạn (Poincaré Section, Attractor Basin). Phạm vi thực nghiệm được giới hạn nghiêm ngặt trong dải tần số kích động 5–10 Hz (tương đương 32–63 rad/s) với áp suất làm việc cơ sở $P_{wh1} = 1.85 - 2.5$ bar và góc nghiêng nêm $\alpha = 37^\circ$. Kết quả định lượng đột phá cho thấy QSAVIM dời tần số bắt đầu cách ly từ 37.7 rad/s (ở hệ ETVIM) xuống 18.5 rad/s, đồng thời kéo giảm tần số cộng hưởng từ 27 rad/s xuống 12.5 rad/s.

Literature Review và Positioning

Nghiên cứu về giảm chấn và cách ly dao động tần số thấp trải qua nhiều giai đoạn tiến hóa. Dòng nghiên cứu khởi thủy tập trung vào việc giảm độ cứng tuyến tính của dầm đàn hồi Euler (Virgin et al., 2008; Chin et al., 2010; Winterflood et al., 2002). Mặc dù các kết cấu này hạ được tần số dao động riêng $\omega_n = \sqrt{k/m}$, nhược điểm chí mạng là độ võng tĩnh tăng vọt theo tỷ lệ nghịch với độ cứng ($H_0 = mg/k$), làm mất ổn định kết cấu chịu tải.

TIẾN TRÌNH PHÁT TRIỂN CỦA CÁC CÔNG NGHỆ CÁCH LY DAO ĐỘNG
=============================================================================
[1. Hệ tuyến tính cổ điển] -> Độ võng tĩnh lớn, không cách ly được tần số < 25 Hz
           |
           v
[2. QZS Cơ học / Từ tính]  -> Carrella (2007), Liu (2013), Ye (2020)
                              Hạn chế: Tải trọng cố định, không tùy biến thích nghi
           |
           v
[3. QSAVIM Khí nén Đột phá]-> Võ Ngọc Yến Phương (2022)
                              Tự do điều áp, triệt tiêu độ cứng động học thích nghi
=============================================================================

Để khắc phục nghịch lý trên, các hệ thống có độ cứng tĩnh cao – độ cứng động học thấp (High-Static-Low-Dynamic-Stiffness - HSLDS hay QZS) đã ra đời. Carrella et al. (2007) thiết lập nền tảng giải tích cho hệ QZS ba lò xo cơ khí; tiếp đó, Dong et al. (2012) và Wang et al. (2014) tích hợp các phần tử nam châm vĩnh cửu nhằm tạo lực đẩy âm. Ye et al. (2020) giới thiệu cơ cấu cam-con lăn tịnh tiến-quay, trong khi Yan et al. (2019) mở rộng cơ cấu Steward QZS 6 bậc tự do.

Tuy nhiên, tồn tại một cuộc tranh luận học thuật lớn giữa hai trường phái:

  1. Trường phái QZS kim loại/từ tính cố định: Ưu tiên kết cấu gọn nhẹ, không cần nguồn cấp năng lượng nhưng hoàn toàn bất lực khi khối lượng tải cách ly $M$ biến động, dẫn đến hiện tượng trôi điểm cân bằng khỏi DSEP và làm hệ thống mất hoàn toàn đặc tính QZS.
  2. Trường phái QZS thích ứng khí nén (Pneumatic Adaptive QZS): Được khởi xướng và phát triển mạnh mẽ bởi các nghiên cứu chuyên sâu (Vo et al., 2020; Le et al., 2021). Nhóm tác giả khẳng định việc ứng dụng lò xo khí cao su (Rubber Air Spring - RAS) và xylanh khí nén (Pneumatic Cylinder - PC) kết hợp buồng tích áp phụ ($V_{ac}$) mang lại khả năng biến thiên độ cứng linh hoạt thông qua điều khiển áp suất khí nén cấp vào hệ thống.

So sánh với hai nghiên cứu quốc tế điển hình:

  • So với nghiên cứu của Liu et al. (2013) về bộ cách ly 5 lò xo có rãnh ngang: Mô hình của Liu et al. đạt QZS thuần túy cơ học nhưng độ cứng âm phụ thuộc hoàn toàn vào kích thước hình học cố định của lò xo kim loại, không thể tinh chỉnh khi tải trọng thay đổi $\pm 5%$. Ngược lại, mô hình QSAVIM của tác giả cho phép tái thiết lập trạng thái QZS tại DSEP bằng cách vi chỉnh tỷ số áp suất $\mu$.
  • So với nghiên cứu của Ye et al. (2020) về cơ cấu cam-con lăn: Cơ cấu của Ye et al. chịu tổn thất ma sát trượt lớn và không tích hợp được tính chất tiêu tán năng lượng nội tại của vật liệu. Luận án của Võ Ngọc Yến Phương tích hợp thành công mô hình ma sát Berg và tính đàn nhớt phân số của cao su, biến tổn hao nội tại thành cơ chế dập tắt dao động tự nhiên mà không cần bổ sung các cụm giảm chấn phụ trợ phức tạp.

Đó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 sâu sắc lý thuyết dao động phi tuyến và cơ học khí dung thông qua việc thiết lập phương trình cấu tạo động học tổng quát của hệ cách ly dao động thích ứng đa vật lý:

$$\mathbf{P1:}\quad F_s(x) = F_{air}(x) + F_{fri}(x) + F_{vie}(x)$$

Trong đó, lực đàn hồi tổng hợp $F_s$ được lượng hóa chính xác bằng nguyên lý cộng dồn ba trường lực:

  1. Trường lực nhiệt động đàn hồi đa biến $F_{air}(x) = A_e(P - P_{atm})$, với diện tích làm việc hiệu dụng $A_e(x)$ và thể tích hiệu dụng $V_e(x)$ được mô hình hóa dưới dạng hàm đa thức phi tuyến bậc cao: $$A_e(x) = (0.3209x^4 + \dots) \times 10^{-3}, \quad V_e(x) = (0.17 + \dots) \times 10^{-3}$$
  2. Trường lực ma sát vật liệu cao su theo mô hình phi tuyến Berg: $$F_{fri} = F_{ref} + \frac{x - x_{ref}}{x_2(1 - a) + |x - x_{ref}|}(F_{fri\max} - F_{ref})$$
  3. Trường lực đàn nhớt suy rộng theo toán tử vi tích phân phân số Kelvin-Voigt: $$F_{vie}(t) = K_e x(t) + b D_t^c x(t) \quad (0 < c < 1)$$

Nghiên cứu xác lập bước chuyển đổi hệ hình (paradigm shift) từ các mô hình QZS tĩnh sang mô hình QZS thích nghi nhiệt-động-nhớt, chứng minh rằng sự bất đối xứng của đường cong độ cứng quanh DSEP có thể được điều tiết thông qua tối ưu hóa thể tích buồng tích áp phụ vô hướng $\hat{V}{ac1}$ và $\hat{V}{ac2}$.

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

Khung phân tích của luận án tích hợp liền mạch ba trụ cột lý thuyết:

  • Định luật thứ nhất Nhiệt động lực học: $dE_{ch} = dE_{in} - dE_{out} - dE_{ae} - \delta Q$, giải bài toán biến thiên năng lượng khí giãn nở đa biến ($n = 1.4$).
  • Mô hình ma sát lai Stribeck-Coulomb-Viscous: $F_{sf} = F_c + (F_{st} - F_c)e^{-(v_r/v_s)^{n_s}} + \sigma v_r$, giải bài toán ma sát tiếp xúc giữa cần piston và nòng xylanh.
  • Phương pháp giải tích xấp xỉ bậc cao (Normal Form và Multi-Scale): Cho phép giải phương trình vi phân phi tuyến bậc 5 dạng:

$$\ddot{u} + \epsilon \mu_d \dot{u} + \alpha_1 u + \alpha_2 u^2 + \alpha_3 u^3 + \alpha_4 u^4 + \alpha_5 u^5 = f_0 \cos(\Omega t)$$

MA TRẬN ĐỘ CỨNG ĐỘNG HỌC KHÔNG THỨ NGUYÊN TẠI VỊ TRÍ CÂN BẰNG (DSEP)
=============================================================================
K_DSEP(mu) = A * (sqrt(1 - H_o^2) - 1) + D + mu * [B * (sqrt(1 - H_o^2) - 1) - C] + 1
Điều kiện đạt Quasi-Zero Stiffness: Cho K_DSEP(mu) = 0
==> Xác định chính xác tỷ số áp suất tối ưu mu = P_wh2 / P_wh1
=============================================================================

Điều kiện biên và phạm vi ứng dụng của mô hình được kiểm soát chặt chẽ: góc nêm $\alpha \in [30^\circ, 45^\circ]$, bán kính con lăn $r \in [10, 60]$ mm, bán kính cam $R = 60$ mm, biến dạng tĩnh $H_0 = 25.6$ mm, và áp suất vận hành tối đa không vượt quá $P_{\max} = 6$ bar.

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

+-----------------------------------------------------------------------------+
|                     QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU                        |
+-------------------------------------+---------------------------------------+
| PHÂN TÍCH LÝ THUYẾT & MÔ PHỎNG      | THỰC NGHIỆM ĐO ĐẠC & NHẬN DẠNG        |
| - Nhiệt động lực học khí nén        | - Card thu thập dữ liệu NI-PCI 6221   |
| - Đạo hàm phân số Kelvin-Voigt      | - Cảm biến lực Loadcell & LVDT        |
| - Phân tích Normal Form & RK4       | - Bàn rung thủy lực Excitable Table   |
| - Mô hình ảo AMEsim xác định ma sát | - Thuật toán Di truyền (GA) tối ưu J  |
+-------------------------------------+---------------------------------------+
                                      |
                                      v
+-----------------------------------------------------------------------------+
| KIỂM CHỨNG & SO SÁNH: QSAVIM (Đề xuất) vs. ETVIM (Truyền thống tương đương) |
+-----------------------------------------------------------------------------+

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (positivism) với phương pháp tam giác hóa toàn diện (methodological and data triangulation). Hệ thống được phân rã và kiểm chứng theo 3 cấp độ:

  1. Cấp độ phần tử (Component level): Đánh giá độc lập đặc tính đàn hồi của lò xo khí cao su và xylanh khí nén có buồng phụ.
  2. Cấp độ cơ cấu (Mechanism level): Phân tích độc lập cơ cấu nâng đỡ mang tải (Load Bearing Mechanism - LBM) và cơ cấu hiệu chỉnh độ cứng âm (Stiffness Corrected Mechanism - SCM).
  3. Cấp độ toàn hệ thống (System level): Khảo sát tích hợp mô hình QSAVIM hoàn chỉnh trên bàn rung kích động thủy lực dưới các tín hiệu điều hòa quét tần (swept-sine) và tín hiệu kích động đa tần phức hợp.

Quy trình nghiên cứu rigorous

Quy trình nhận dạng tham số vật lý và kiểm chuẩn thực nghiệm được tiến hành nghiêm ngặt:

  • Phần cứng thu thập dữ liệu: Sử dụng card xử lý tín hiệu tốc độ cao National Instruments NI-Card 6221 kết nối máy tính chủ; cảm biến chuyển vị tuyến tính LVDT; cảm biến gia tốc áp điện trở; loadcell đo lực nén kéo động; đồng hồ áp suất số và cụm van điều áp tỷ lệ khí nén điện tử.
  • Nhận dạng tham số mô hình đàn nhớt phân số: Thiết lập thí nghiệm kích động cưỡng bức ở các dải tần $f \in {0.5, 3, 4, 5, 6, 7}$ Hz với biên độ cố định 8 mm. Tích hợp thuật toán Di truyền (Genetic Algorithm - GA) để tối thiểu hóa hàm mục tiêu sai số bình phương trung bình:

$$J = \frac{1}{m} \sum_{i=1}^{m} \left[ \left( \frac{F_{vie1i}}{x_{i0}} - (K_{vie} + b_i \cos(c\pi/2)) \right)^2 + \left( \frac{F_{vie2i}}{x_{i0}} - b_i \sin(c\pi/2) \right)^2 \right]$$

Xác định chính xác bộ tham số tối ưu: $K_{vie} = 1.05 \times 10^4$ N/m, $b = 850$ N$\cdot\text{s}^c$/m, $c = 0.65$.

Data và phân tích

  • Mô phỏng nguyên mẫu ảo (Virtual Prototyping): Xây dựng mô hình số động lực học chất lưu và ma sát trong phần mềm LMS Imagine.Lab AMEsim nhằm nhận dạng chính xác đường đặc tính ma sát Stribeck của xylanh khí mà không bị nhiễu bởi rung động cơ học ngoại vi.
  • Giải tích phi tuyến và kiểm chứng số: Tích hợp thuật toán Runge-Kutta bậc 4 (RK4) với bước thời gian phân chia cực mịn $\Delta t = T/100$ (với $T = 2\pi/\Omega$).
  • Mặt cắt Poincaré (Poincaré Section) & Bản đồ Lưu vực Hấp dẫn (Attractor-Basin Phase Portrait): Khảo sát tính ổn định nghiệm trạng thái dừng, xác định tọa độ các điểm cân bằng hút (attractor points) và các quỹ đạo pha phân nhánh trong không gian trạng thái $[u, \dot{u}]$.

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

SO SÁNH ĐẶC TÍNH ĐỘNG LỰC HỌC GIỮA HAI HỆ THỐNG
=============================================================================
Đặc tính động học                 ETVIM (Truyền thống)    QSAVIM (Nghiên cứu)
-----------------------------------------------------------------------------
Tần số cộng hưởng riêng           27.0 rad/s (~4.3 Hz)    12.5 rad/s (~1.99 Hz)
Tần số bắt đầu cách ly (T < 1)   37.7 rad/s (~6.0 Hz)    18.5 rad/s (~2.94 Hz)
Hiện tượng rẽ nhánh khi lệch tải  Không có                Xuất hiện chu kỳ đôi
Khả năng tự điều chỉnh độ cứng   Không thể               Bằng áp suất khí nén
=============================================================================

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

  1. Mở rộng vượt bậc dải tần số cách ly hiệu dụng: Hệ thống QSAVIM sử dụng lò xo khí cao su với tỷ số áp suất $\mu = 1.8$ bắt đầu cách ly dao động tại tần số $\Omega > 18.5$ rad/s, vượt trội hoàn toàn so với mô hình truyền thống ETVIM vốn bất lực ở mọi tần số dưới $37.7$ rad/s. Tần số cộng hưởng đỉnh giảm mạnh từ $27$ rad/s xuống $12.5$ rad/s (mức giảm $53.7%$).
  2. Phát hiện và định vị chính xác hiện tượng nhảy tần số (Frequency Jump): Khi quét tần số thuận nghịch từ $0 \to 3$ rad/s và ngược lại, đường cong đáp ứng biên độ-tần số bẻ cong về phía bên phải (hiện tượng làm cứng - hardening behavior), làm phát sinh bước nhảy biên độ lên tại $\Omega = 1.7$ rad/s và bước nhảy xuống tại $\Omega = 2.2$ rad/s.
  3. Hiện tượng rẽ nhánh chu kỳ phụ thuộc độ lệch tải trọng tĩnh: Khi tải trọng cách ly lệch $\pm 5%$ so với tải trọng thiết kế tối ưu ($M_{opt} = 4.747$ kg, thay đổi thành $4.510$ kg hoặc $4.984$ kg), hệ thống chuyển pha từ dao động đơn chu kỳ (Period-1) sang dao động đa chu kỳ (Period-2 và Period-3) trong các dải tần 8–10.5 rad/s và 13–14.5 rad/s, làm tăng đáng kể biên độ dao động nếu không được điều áp bù trừ.
  4. Hiệu ứng bất đối xứng độ cứng của xylanh khí buồng tích áp: Khác với lò xo khí cao su có đường cong độ cứng đối xứng dạng parabol lõm hoàn hảo, việc tăng thể tích buồng phụ $\hat{V}{ac1}$ ở cơ cấu LBM làm giảm độ cứng tổng thể nhưng tạo ra độ bất đối xứng quanh u=0, làm dịch chuyển điểm có độ cứng cực tiểu từ $u = 0$ sang vị trí $u{ls} \approx 0.084$.

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung cơ chế mô tả toán học chính xác cho hệ dao động tự điều chỉnh độ cứng phi tuyến, kết nối giữa cơ học môi trường liên tục và lý thuyết điều khiển khí nén.
  • Về mặt phương pháp luận: Cung cấp quy trình thiết kế chuẩn hóa 6 bước (từ xác định tải trọng $M \to$ lựa chọn hình học $R, r, \alpha \to$ tính toán áp suất $P_{wh1} \to$ xác định tỷ số áp suất $\mu \to$ tối ưu hóa thể tích buồng phụ $V_{ac} \to$ mô phỏng kiểm chứng chuyển vị và gia tốc).
  • Về mặt ứng dụng thực tiễn: Cung cấp giải pháp chế tạo hoàn chỉnh có khả năng nội địa hóa 100% tại Việt Nam cho:
    1. Hệ thống treo cabin và ghế ngồi tài xế xe tải hạng nặng, xe công trình nhằm triệt tiêu rung chấn bảo vệ cột sống.
    2. Bệ cách ly chống rung cho kính hiển vi điện tử quét (SEM), máy đo biên dạng 3D quang học và thiết bị gia công siêu chính xác.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật và kỹ thuật còn tồn tại:

  • Biên độ cộng hưởng tại tần số thấp còn cao: Dù tần số cộng hưởng giảm sâu về 12.5 rad/s, hệ số khuếch đại dao động tại đỉnh cộng hưởng vẫn tồn tại do hệ số cản nhớt nội tại của vật liệu cao su ở chế độ thụ động chưa đủ lớn để dập tắt hoàn toàn năng lượng tích tụ.
  • Đặc tính ma sát trượt của xylanh khí: Sự hiện diện của ma sát tĩnh $F_{st}$ và ma sát Coulomb $F_c$ tại các phớt làm kín của piston gây ra hiện tượng giật cục (stick-slip) ở biên độ dịch chuyển vi mô ($<0.5$ mm).
  • Chế độ thích nghi bán chủ động: Quá trình điều chỉnh tỷ số áp suất $\mu$ trong nghiên cứu hiện tại được thực hiện thủ công thông qua van giảm áp, chưa tích hợp vòng điều khiển phản hồi kín tự động thời gian thực (Real-time Closed-loop Control).

Chương trình nghiên cứu tiếp nối trong giai đoạn 5–10 năm tới bao gồm:

  1. Phát triển thuật toán điều khiển phi tuyến thích nghi bền vững ($H_\infty$, Sliding Mode Control kết hợp mạng Neural nhân tạo) để điều khiển đóng mở van servo khí nén tốc độ cao thời gian thực.
  2. Tích hợp vật liệu lưu chất từ trường (Magnetorheological Fluid - MRF) vào buồng giảm chấn nhằm dập tắt tức thời biên độ đỉnh cộng hưởng.
  3. Mở rộng mô hình tính toán từ 1 bậc tự do (1-DOF) lên hệ không gian 6 bậc tự do (6-DOF) ứng dụng cho bệ phóng hàng không vũ trụ và bàn thử nghiệm mô phỏng động đất.

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

  • Tác động học thuật: Mở ra một hướng nghiên cứu liên ngành chuyên sâu giữa Kỹ thuật Cơ khí, Điều khiển Tự động và Khoa học Vật liệu phi kim loại; công bố các bài báo khoa học chất lượng cao trên các tạp chí quốc tế uy tín (như Mechanical Systems and Signal Processing - MSSP, Journal of Vibroengineering).
  • Chuyển giao công nghệ công nghiệp: Tạo tiền đề kỹ thuật vững chắc cho các nhà sản xuất ô tô, xe cơ giới và thiết bị công nghiệp trong nước (như THACO, VinFast, các viện nghiên cứu cơ khí) chế tạo các cụm giảm xóc khí nén thông minh nội địa giá thành thấp thay thế hàng nhập khẩu đắt đỏ.
  • Lợi ích xã hội và y sinh: Giảm thiểu nguy cơ mắc bệnh nghề nghiệp về cột sống và tổn thương thần kinh vận động cho hơn 1.5 triệu tài xế xe buýt, xe tải đường dài và công nhân vận hành máy kéo tại Việt Nam.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học/Cơ điện tử: Tiếp cận phương pháp luận giải tích phi tuyến tiên tiến (Normal Form, Multi-Scale, Poincaré Map) kết hợp tối ưu hóa thuật toán di truyền (GA) ứng dụng trong dao động kỹ thuật.
  • Kỹ sư R&D trong ngành Công nghiệp Ô tô & Cơ khí chính xác: Bản vẽ thiết kế nguyên lý, quy chuẩn lựa chọn kích thước cam-nêm và công thức tính toán áp suất khí nén để ứng dụng trực tiếp vào sản phẩm thương mại.
  • Các nhà hoạch định chính sách lao động và an toàn giao thông: Căn cứ khoa học định lượng về rung động tần số thấp (<5 Hz) phục vụ xây dựng và cập nhật Tiêu chuẩn Quốc gia (TCVN/QCVN) về giới hạn rung động cho phép trên phương tiện giao thông vận tải.

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 xây dựng thành công mô hình giải tích độ cứng động học phi tuyến tổng quát cho hệ QZS sử dụng phần tử đàn hồi phi kim loại, mở rộng đồng thời lý thuyết nhiệt động lực học khí nén trong buồng thể tích biến thiên và lý thuyết đàn nhớt phân số Kelvin-Voigt ($F_{vie} = K_e x + b D_t^c x$). Luận án chứng minh chính xác toán học rằng việc kết hợp cơ cấu cam-nêm với phần tử khí biến thiên cho phép triệt tiêu độ cứng động học ($K_{DSEP} = 0$) tại bất kỳ vị trí cân bằng tĩnh nào thông qua điều chỉnh tỷ số áp suất $\mu$.

2. Phương pháp nghiên cứu có tính đổi mới như thế nào so với các nghiên cứu trước?

So với nghiên cứu của Carrella et al. (2007) và Liu et al. (2013), nghiên cứu này đổi mới toàn diện ở phương pháp nhận dạng đa vật lý: kết hợp mô hình hóa nguyên mẫu ảo trên phần mềm chuyên dụng AMEsim để bóc tách ma sát Stribeck của xylanh khí, đồng thời ứng dụng thuật toán di truyền (GA) để khớp nối tham số mô hình Berg và Kelvin-Voigt từ dữ liệu đo thực nghiệm thực tế với card NI-PCI 6221.

3. Phát hiện nào gây bất ngờ nhất về mặt động lực học phi tuyến và bằng chứng thực nghiệm là gì?

Phát hiện bất ngờ nhất là sự dịch chuyển của lưu vực hấp dẫn (attractor basin) và hiện tượng rẽ nhánh chu kỳ đôi khi khối lượng tải chỉ lệch $5%$ so với thiết kế tĩnh ($M = 4.984$ kg so với $M_{opt} = 4.747$ kg). Tại dải tần 1.7–2.2 rad/s, hệ thống tồn tại đồng thời hai trạng thái nghiệm: cộng hưởng và không cộng hưởng. Trạng thái dừng cuối cùng hoàn toàn bị chi phối bởi điều kiện ban đầu $(u_0, \dot{u}_0)$, được chứng minh qua sơ đồ Poincaré và mô phỏng số RK4.

4. Luận án có cung cấp quy trình nhân bản (replication protocol) chi tiết không?

Có. Luận án cung cấp quy trình thiết kế 6 bước chi tiết kèm theo công thức giải tích tường minh (phương trình 4.29 - 4.31 và 5.28 - 5.29), các thông số hình học chuẩn ($R = 60$ mm, $r = 20$ mm, $\alpha = 37^\circ$, $H_0 = 25.6$ mm) và sơ đồ mạch đo thực nghiệm phần cứng hoàn chỉnh kết nối qua card NI-6221.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình ra sao?

Lộ trình 10 năm gồm 3 giai đoạn: (1) 2023–2025: Phát triển bộ điều khiển thích nghi thời gian thực tích hợp van servo khí nén thông minh; (2) 2026–2028: Chế tạo và thử nghiệm hệ thống treo bán chủ động kết hợp giảm chấn từ trường MRF trên xe thực tế; (3) 2029–2033: Mở rộng sang hệ thống cách ly đa hướng 6 bậc tự do phục vụ công nghiệp bán dẫn và bệ thử hàng không vũ trụ.

Kết luận

Công trình nghiên cứu tiến sĩ của tác giả Võ Ngọc Yến Phương đã đạt được 5 đóng góp học thuật cốt lõi:

  1. Thiết lập thành công mô hình toán học giải tích toàn diện mô tả đặc tính lực đàn hồi đa thành phần (khí nén, ma sát Berg, đàn nhớt phân số Kelvin-Voigt) của lò xo khí cao su và xylanh khí nén có buồng phụ.
  2. Thiết kế và chế tạo thành công kết cấu cơ khí QSAVIM tích hợp cơ cấu nêm-cam-con lăn, cho phép đạt độ cứng động học tiệm cận không ($K_{DSEP} \approx 0$) tại vị trí cân bằng tĩnh mong muốn.
  3. Làm sáng tỏ bản chất động lực học phi tuyến phức tạp của hệ thống, bao gồm hiện tượng nhảy tần số (jump-up tại 1.7 rad/s, jump-down tại 2.2 rad/s), hiện tượng rẽ nhánh dao động khi lệch tải trọng tĩnh và sự phân bố lưu vực hấp dẫn theo điều kiện ban đầu.
  4. Chứng minh bằng thực nghiệm trên bàn rung thủy lực khả năng vượt trội của QSAVIM: hạ tần số cộng hưởng từ 27 rad/s xuống 12.5 rad/s và mở rộng dải tần cách ly bắt đầu từ 18.5 rad/s (so với ngưỡng 37.7 rad/s của hệ truyền thống ETVIM).
  5. Xây dựng quy trình thiết kế và chế tạo hoàn chỉnh có tính khả thi cao, mở đường cho việc sản xuất và thương mại hóa các hệ thống cách ly dao động thích nghi tần số thấp chất lượng cao tại Việt Nam.