Tổng quan về luận án

Nghiên cứu động lực học máy hiện đại đang đứng trước những thách thức kỹ thuật to lớn khi xu hướng thiết kế công nghiệp đòi hỏi tốc độ vận hành ngày càng cao cùng với yêu cầu tối ưu hóa, giảm nhẹ trọng lượng kết cấu. Điều này dẫn đến sự xuất hiện tất yếu của các hiện tượng dao động phi tuyến mạnh và kích động tham số phức tạp ở chế độ chuyển động bình ổn. Luận án tiến sĩ kỹ thuật "Tính toán dao động tuần hoàn và rẽ nhánh của một số mô hình dao động trong máy" do NCS. Hoàng Mạnh Cường thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Văn Khang (Chuyên ngành Cơ học kỹ thuật, Mã số: 62.01, Trường Đại học Bách khoa Hà Nội, 2011) là một công trình tiên phong giải quyết trọn vẹn bài toán xác định chế độ dao động tuần hoàn phi tuyến và phân tích rẽ nhánh trong các hệ thống máy thực tế.

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án tập trung giải quyết xuất phát từ giới hạn căn bản của các phương pháp giải tích truyền thống:

"Ngay từ nửa đầu thế kỷ 20, nhiều nhà khoa học đã quan tâm nghiên cứu các bài toán dao động và ổn định của các hệ phi tuyến yếu... Tuy vậy việc ứng dụng các kết quả nghiên cứu về động lực học phi tuyến của các nhà toán học vào các hệ vật lý và các hệ cơ học vẫn còn là vấn đề đang được nhiều người quan tâm nghiên cứu. Đặc biệt là việc áp dụng các phương pháp số, tính toán dao động tuần hoàn và rẽ nhánh của các mô hình dao động phi tuyến mạnh trong máy là vấn đề rất thời sự hiện nay."

Các câu hỏi nghiên cứu và giả thuyết khoa học được định hình cụ thể:

  • RQ1: Làm thế nào để xác định chính xác điều kiện đầu $x(0) = x(T)$ của nghiệm tuần hoàn chu kỳ $T$ trong các hệ vi phân tuyến tính hệ số tuần hoàn và phi tuyến mạnh mà không bị phân kỳ hoặc tích lũy sai số khi tích phân số dài hạn?
  • RQ2: Cơ chế mất ổn định và quy luật rẽ nhánh của nghiệm tuần hoàn trong các mô hình dao động cơ học máy biến đổi như thế nào khi các tham số kích động vượt qua giá trị tới hạn?
  • H1: Việc kết hợp phương pháp tích phân số Runge-Kutta bậc 4 với định lý Floquet cho phép tính toán trực tiếp ma trận đơn đạo $C = e^{BT}$, từ đó xác định chính xác phổ nhân tử đặc trưng $\lambda_k$ để đánh giá độ ổn định Liapunov.
  • H2: Ứng dụng phương pháp bắn đơn (single shooting method) kết hợp thuật toán lặp Newton-Raphson cho phép tìm nghiệm tuần hoàn ổn định và không ổn định trong các hệ cơ học phi tuyến mạnh, đồng thời nhận diện chính xác các điểm rẽ nhánh Hopf, rẽ nhánh nhân đôi chu kỳ (period-doubling) và rẽ nhánh phá vỡ tính đối xứng (symmetry-breaking).

Khung lý thuyết của luận án được xây dựng vững chắc trên nền tảng lý thuyết ổn định chuyển động của A.M. Liapunov (1892), lý thuyết Floquet (1883) cho hệ phương trình vi phân tuyến tính tuần hoàn, phương pháp trung bình hoá tiệm cận Bogoliubov-Mitropolsky (1961) và lý thuyết rẽ nhánh hiện đại của Poincaré-Andronov. Luận án khảo sát phạm vi 4 mô hình dao động máy điển hình: bộ truyền bánh răng một cấp có xét đến sai số bước răng và khe hở, cơ cấu cam trong bộ điều chỉnh máy rèn dập, dao động xoắn hệ trục máy tàu thuỷ 4 xi lanh DTSC50, và các hệ dao động tự kích/tuyến tính từng khúc (Van der Pol-Duffing, piecewise linear systems).

Literature Review và Positioning

Nghiên cứu dao động phi tuyến đã trải qua lịch sử phát triển sâu rộng. Dòng nghiên cứu giải tích cổ điển gắn liền với các công trình của Krylov & Bogoliubov (1937), Bogoliubov & Mitropolsky (1961), và Nayfeh & Mook (1979), tập trung vào các phương pháp tiệm cận, phương pháp tham số bé, phương pháp trung bình hoá và phương pháp nhiều quy mô thời gian (method of multiple scales). Tuy nhiên, các phương pháp này bộc lộ hạn chế nghiêm trọng khi hệ số kích động không nhỏ hoặc khi đặc tính phi tuyến mang tính gián đoạn (như va đập, khe hở ăn khớp bánh răng, ma sát khô).

   [Các phương pháp giải tích cổ điển]
   (Krylov-Bogoliubov, Nayfeh, Mitropolsky)
         │  (Hạn chế: Chỉ áp dụng cho hệ phi tuyến yếu ε << 1)
         ▼
┌─────────────────────────────────────────────────────────────┐
│ TRANH LUẬN HỌC THUẬT: GIẢI TÍCH vs. PHƯƠNG PHÁP SỐ          │
├──────────────────────────────┬──────────────────────────────┤
│ Tích phân thời gian thuần túy│ Phương pháp giải tích mở rộng│
│ - Tốn kém thời gian tính toán│ - Không xử lý được hệ gián   │
│ - Không bắt được nghiệm mất  │   đoạn (khe hở, va đập)      │
│   ổn định & điểm rẽ nhánh    │ - Giới hạn bậc tự do nhỏ     │
└──────────────┬───────────────┴──────────────┬───────────────┘
               │                              │
               └──────────────┬───────────────┘
                              ▼
        [ĐỊNH VỊ ĐỘT PHÁ CỦA LUẬN ÁN (HOÀNG MẠNH CƯỜNG, 2011)]
        - Tích hợp Phương pháp Bắn đơn (Single Shooting)
        - Lý thuyết ổn định Floquet & Ma trận đơn đạo C = e^(BT)
        - Thuật toán tích phân Runge-Kutta-Nyström tối ưu
        - Phân tích đa dạng rẽ nhánh (Hopf, Flip, Fold, Symmetry)

Sự phát triển của lý thuyết hệ động lực phi tuyến hiện đại (Guckenheimer & Holmes, 1983; Parker & Chua, 1989; Seydel, 1994) đã mở ra các hướng tiếp cận số học như phương pháp cân bằng điều hoà gia lượng (incremental harmonic balance), phương pháp ánh xạ Poincaré và phương pháp bắn (shooting method).

Luận án định vị vị trí học thuật thông qua việc so sánh đối chuẩn trực tiếp với hai nghiên cứu quốc tế điển hình:

  1. Nghiên cứu của Kahraman & Singh (1990, 1991): Khảo sát động lực học phi tuyến của cặp bánh răng có khe hở và độ cứng thay đổi theo thời gian bằng phương pháp cân bằng điều hoà và tích phân số. Luận án của Hoàng Mạnh Cường đã mở rộng khi không chỉ xác định quỹ đạo dao động mà còn xây dựng thuật toán bắn đơn tự động xác định điều kiện đầu chính xác, đồng thời theo vết nhân tử Floquet để vẽ biểu đồ rẽ nhánh chi tiết khi thay đổi tần số quay và biên độ kích động.
  2. Nghiên cứu của Theodossiades & Natsiavas (2000): Nghiên cứu hành vi rẽ nhánh và chuyển động hỗn độn trong hệ cơ học tuyến tính từng khúc. Luận án đã phát triển thuật toán số dựa trên phương pháp Runge-Kutta-Nyström cho phép khảo sát chính xác hiện tượng rẽ nhánh nhân đôi chu kỳ (period-doubling cascade) và rẽ nhánh phá vỡ tính đối xứng trong hệ tuyến tính từng khúc nhiều bậc tự do.

Đó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 và hoàn thiện khung lý thuyết ổn định Liapunov và Floquet đối với các hệ phương trình vi phân tham số có ma trận quán tính biến thiên theo thời gian $M(t)\ddot{q} + C(t)\dot{q} + K(t)q = h(t)$.

Hệ phương trình vi phân chuyển động tổng quát của cơ hệ nhiều vật được thiết lập thông qua phương trình Lagrange loại II: $$\frac{d}{dt}\left(\frac{\partial T}{\partial \dot{q}_i}\right) - \frac{\partial T}{\partial q_i} + \frac{\partial \Phi}{\partial \dot{q}_i} + \frac{\partial \Pi}{\partial q_i} = Q_i \quad (i = 1, \dots, f)$$

Đóng góp lý thuyết then chốt nằm ở việc hệ thống hoá các điều kiện hình học và phổ của nhân tử Floquet $\lambda_k$ (nghiệm của phương trình đặc trưng $\det[C - \lambda E] = 0$) trên mặt phẳng phức đối với vòng tròn đơn vị ${z \in \mathbb{C} : |z| = 1}$:

  • Mệnh đề 1 (Ổn định tiệm cận): Hệ ổn định tiệm cận khi và chỉ khi tất cả các nhân tử đặc trưng nằm hoàn toàn trong vòng tròn đơn vị, tức $|\lambda_k| < 1, \forall k = 1, \dots, n$.
  • Mệnh đề 2 (Rẽ nhánh nút-yên ngựa / Fold): Một nhân tử Floquet thực rời vòng tròn đơn vị tại điểm $+1$ ($\lambda = +1$).
  • Mệnh đề 3 (Rẽ nhánh nhân đôi chu kỳ / Flip): Một nhân tử Floquet thực rời vòng tròn đơn vị tại điểm $-1$ ($\lambda = -1$), làm xuất hiện chu trình giới hạn mới với chu kỳ $2T$.
  • Mệnh đề 4 (Rẽ nhánh Neimark-Sacker / Hopf loại 2): Một cặp nhân tử Floquet phức liên hợp rời vòng tròn đơn vị tại $\lambda_{1,2} = e^{\pm i\theta} \ (\theta \neq 0, \pi)$, dẫn đến sự hình thành một xuyến hai chiều (torus) trong không gian pha.

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

Khung phân tích của luận án kết hợp đa phương pháp:

  1. Dạng chuẩn Lagrange-Bogoliubov và phương pháp trung bình hoá: Áp dụng cho các hệ phi tuyến yếu dạng $\dot{y} = \varepsilon f(t, y)$, giải tích hóa quy luật rẽ nhánh Hopf và rẽ nhánh nút-yên ngựa của các hệ Duffing, Mathieu và Van der Pol-cưỡng bức.
  2. Phương pháp bắn đơn tìm nghiệm biên hai điểm: Đưa bài toán tìm nghiệm tuần hoàn $x(t)$ với chu kỳ $T$ về bài toán giải hệ phương trình phi tuyến đối với điều kiện đầu: $$H(x_0) = \phi(T, x_0) - x_0 = 0$$ với $\phi(T, x_0)$ là trạng thái của hệ tại thời điểm $t = T$ xuất phát từ $x(0) = x_0$.
  3. Điều kiện biên xác định (Boundary Conditions): Thiết lập ma trận Jacobi quá độ $J = \frac{\partial \phi(T, x_0)}{\partial x_0} - E$ phục vụ vòng lặp hiệu chỉnh bậc hai Newton-Raphson: $$x_0^{(k+1)} = x_0^{(k)} - \left[\frac{\partial \phi(T, x_0^{(k)})}{\partial x_0} - E\right]^{-1} \left(\phi(T, x_0^{(k)}) - x_0^{(k)}\right)$$

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 lý nghiên cứu thực chứng kỹ thuật (positivist / deductive paradigm), kết hợp giữa mô hình hoá giải tích cơ học giải tích với tính toán số học ma trận và phân tích phổ tín hiệu thực nghiệm.

Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 (Mô hình hoá vật lý - toán): Thiết lập phương trình vi phân chuyển động phi tuyến thông qua hình thức giải tích Lagrange.
  • Tầng 2 (Thuật toán số hoá ma trận): Tuyến tính hoá quanh quỹ đạo cơ bản, biểu diễn các hàm kích động dưới dạng chuỗi Fourier hữu hạn nhiều phần tử.
  • Tầng 3 (Giải thuật tích phân & Ổn định): Tích phân số Runge-Kutta-Nyström, lập ma trận đơn đạo $C$, thuật toán lặp Von-Mises tìm trị riêng cực đại $|\lambda_{\max}|$.
  • Tầng 4 (Thực nghiệm đối chứng): Đối chiếu phổ tần số tính toán với phổ gia tốc đo đạc trên bệ thử nghiệm động lực học.

Quy trình nghiên cứu rigorous

Quy trình xác định nghiệm tuần hoàn của hệ phương trình tham số tuyến tính $M(t)\ddot{x} + C(t)\dot{x} + K(t)x = d(t)$ được thực hiện chặt chẽ qua các bước:

  1. Đưa phương trình cấp 2 về hệ phương trình trạng thái cấp 1: $\dot{x} = P(t)x + f(t)$.
  2. Chia chu kỳ $T$ thành $m$ khoảng bằng nhau với bước tích phân $h = T/m$.
  3. Áp dụng sơ đồ Runge-Kutta bậc 4, thiết lập quan hệ chuyển tiếp trạng thái: $$x_i = A_{i-1}x_{i-1} + b_{i-1} \quad (i = 1, \dots, m)$$
  4. Tích lũy ma trận chuyển trạng thái trên toàn chu kỳ để nhận được: $$x_m = \left(\prod_{i=1}^m A_{m-i}\right) x_0 + c_m$$
  5. Do điều kiện tuần hoàn $x(T) = x_m = x_0$, giải hệ phương trình đại số tuyến tính: $$\left[E - \prod_{i=1}^m A_{m-i}\right] x_0 = c_m$$ để thu được chính xác véc-tơ điều kiện đầu $x_0$.
  6. Áp dụng thuật toán lặp Von-Mises để tìm giá trị riêng lớn nhất của ma trận đơn đạo: $$a^{(k+1)} = C a^{(k)}, \quad \lambda_{\max} \approx \frac{(a^{(k+1)})^T a^{(k)}}{(a^{(k)})^T a^{(k)}}$$
   ┌──────────────────────────────────────────────────────────┐
   │    HỆ PHƯƠNG TRÌNH TRẠNG THÁI: dx/dt = P(t)x + f(t)       │
   └────────────────────────────┬─────────────────────────────┘
                                │
                                ▼
   ┌──────────────────────────────────────────────────────────┐
   │ Rời rạc hóa [0, T] thành m khoảng chia: h = T/m          │
   │ Tích phân Runge-Kutta bậc 4: x_i = A_(i-1) x_(i-1) + b_  │
   └────────────────────────────┬─────────────────────────────┘
                                │
                                ▼
   ┌──────────────────────────────────────────────────────────┐
   │ THIẾT LẬP ĐIỀU KIỆN BIÊN TUẦN HOÀN: x(0) = x(T) = x_m    │
   │ Giải hệ đại số: [E - Π A_i] x_0 = c_m                    │
   └────────────────────────────┬─────────────────────────────┘
                                │
                                ▼
   ┌──────────────────────────────────────────────────────────┐
   │ XÁC ĐỊNH MA TRẬN ĐƠN ĐẠO C = Φ(T) VÀ NHÂN TỬ FLOQUET    │
   │ Thuật toán lặp Von-Mises: Xác định λ_max                 │
   └────────────────────────────┬─────────────────────────────┘
                                │
        ┌───────────────────────┴───────────────────────┐
        ▼                                               ▼
┌──────────────────────────────┐        ┌──────────────────────────────┐
│       |λ_max| < 1            │        │       |λ_max| ≥ 1            │
│  HỆ ỔN ĐỊNH TIỆM CẬN         │        │  HỆ MẤT ỔN ĐỊNH / RẼ NHÁNH   │
│  Tích phân Runge-Kutta-      │        │  Phân tích rẽ nhánh:         │
│  Nyström vẽ quỹ đạo xác lập  │        │  Fold (+1), Flip (-1), Hopf  │
└──────────────────────────────┘        └──────────────────────────────┘

Data và phân tích

Luận án triển khai tính toán trên tập thông số thực nghiệm kỹ thuật cụ thể:

Mô hình nghiên cứu Thông số hình học & Khối lượng Thông số đàn hồi & Cản Kích động & Tần số
1. Bộ truyền bánh răng 1 cấp $J_1 = 0.093\text{ kg}\cdot\text{m}^2$, $J_2 = 0.272\text{ kg}\cdot\text{m}^2$, $r_1 = 30.46\times 10^{-3}\text{ m}$, $z_1 = 14, z_2 = 39$ $k_0 = 8.1846\times 10^6\text{ N/m}$, 10 hài Fourier cho $k(t)$, Độ cản Lehr $D = 0.024$ $\omega_1 = 60\text{ rad/s}$, Sai số tĩnh $q_0 = 1.2\times 10^{-5}\text{ m}$, 6 hài sai số $e(t)$
2. Bộ điều chỉnh Cam máy dập $I_1 = 1.11\text{ kg}\cdot\text{m}^2$, $m_2 = 136\text{ kg}$ $k_1 = 7692\text{ Nm/rad}$, $k_2 = 10^6\text{ N/m}$, $c_1 = 18.5\text{ Nms}$, $c_2 = 500\text{ Ns/m}$ Tốc độ quay $n = 50\text{ rpm}$ ($\omega = 5.236\text{ rad/s}$), Hàm truyền $U(\omega t)$ gồm 12 hài
3. Hệ trục tàu thuỷ DTSC50 $J_{01} = 0.03825$, $J_{02} = 0.04022$, $J_{03} = 0.03992$, $J_{04} = 0.03831\text{ kg}\cdot\text{m}^2$ $k_1 = k_2 = k_3 = 1.2\times 10^6\text{ Nm/rad}$, $k_5 = 0.1166\times 10^6$, $k_6 = 0.2917\times 10^6$, $c_i = 1.484$ $n = 1700\text{ rpm}$ ($\omega = 178\text{ rad/s}$), $M_0 = 60.38\text{ Nm}$, Mômen 48 hài Fourier

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

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

  1. Xác lập định lượng tính ổn định tiệm cận của bộ truyền bánh răng: Tính toán ma trận đơn đạo của bộ truyền bánh răng: $$C = \begin{bmatrix} 0.01234 & 0.519742\times 10^{-6} \ -15.121177 & 0.012356 \end{bmatrix}$$ cho cặp nhân tử Floquet liên hợp: $$\lambda_{1,2} = 0.012348 \pm 0.0028034i \implies |\lambda_{\max}| = 0.012662$$ Do $|\lambda_{\max}| \ll 1$, hệ dao động tham số bánh răng hoàn toàn ổn định tiệm cận ở vận tốc góc $\omega_1 = 60\text{ rad/s}$.

  2. Khám phá cấu trúc phổ tần số sai số truyền động động lực $q(t)$: Phổ tần số của đạo hàm sai số truyền động $\dot{q}(t)$ bộc lộ rõ hai thành phần:

    • Các đỉnh phổ chính tại các tần số ăn khớp $n\omega_z = n \cdot z_1 \omega_1$ đặc trưng cho kích động tham số độ cứng $k(t)$.
    • Các dải phổ phụ (sidebands) bao quanh các đỉnh chính cách nhau một khoảng đúng bằng tần số quay $\omega_1 = 60\text{ rad/s}$, xuất phát từ kích động sai số bước răng $e(t)$. Kết quả tính toán số này hoàn toàn trùng khớp với dữ liệu đo đạc thực nghiệm.

  3. Tính ổn định của bộ điều chỉnh cơ cấu cam máy rèn dập: Ở tốc độ $n = 50\text{ rpm}$, tính toán ma trận đơn đạo cho cả 2 trường hợp profile cam đều cho thấy tính ổn định cao:

    • Trường hợp 1: $|\lambda_{\max}| = 0.001997 < 1$.
    • Trường hợp 2: $|\lambda_{\max}| = 0.001454 < 1$. Phổ dao động trạng thái chủ yếu bao gồm các sóng hài bậc chẵn ($2\omega, 4\omega, 6\omega$) hoặc bậc lẻ ($\omega, 3\omega, 5\omega$).

  4. Triệt tiêu sai số tích phân trong hệ trục tàu thuỷ 7 bậc tự do: Xác định chính xác véc-tơ điều kiện đầu tuần hoàn tại $\omega = 178\text{ rad/s}$: $$q(0) = [0.055281, 0.055376, 0.055538, 0.055750, 0.055844, 0.053834]^T\text{ (rad)}$$ $$\dot{q}(0) = [0.150767, 0.274259, 0.451432, 0.679615, 0.639052, 0.609877]^T\text{ (rad/s)}$$ giải quyết dứt điểm hiện tượng dao động ký sinh nhân tạo khi mô phỏng số hệ trục máy tàu.

  5. Hiện tượng rẽ nhánh nhân đôi chu kỳ trong hệ phi tuyến có khe hở: Trong hệ truyền động bánh răng có khe hở và hệ tuyến tính từng khúc, khi tăng tần số kích động hoặc giảm hệ số cản, nhân tử Floquet đi qua điểm $-1$, gây rẽ nhánh nhân đôi chu kỳ liên tiếp (Feigenbaum cascade), mở đường cho chuyển động hỗn độn (chaotic motion).

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp phương pháp luận tổng quát kết nối lý thuyết dao động tham số tuyến tính và phi tuyến mạnh, lấp đầy khoảng trống giữa lý thuyết giải tích trừu tượng và ứng dụng kỹ thuật máy.
  • Về mặt phương pháp: Giảm thiểu tới 85% chi phí tính toán so với phương pháp tích phân thời gian dò tìm xác lập thông thường nhờ thuật toán tìm điều kiện đầu $x_0$ trực tiếp.
  • Về kỹ thuật thực tiễn:
    • Đưa ra tiêu chuẩn nhận diện sớm hư hỏng bánh răng dựa trên độ lớn dải phổ phụ xung quanh tần số ăn khớp.
    • Tối ưu hóa profin cam máy dập nhằm tránh vùng cộng hưởng tham số bậc cao.
    • Tính toán chính xác dao động xoắn hệ trục động cơ diesel tàu thuỷ DTSC50 để bố trí bánh đà và bộ giảm chấn xoắn tối ưu.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án vẫn có một số giới hạn khoa học cần được nhìn nhận khách quan:

  1. Giả thiết nguồn năng lượng vô hạn (Ideal Energy Source): Tác giả giả thiết vận tốc góc của trục dẫn là hằng số ($\omega = \text{const}$), bỏ qua hiện tượng tương tác phi tuyến giữa nguồn động lực và cơ hệ dao động (hiệu ứng Sommerfeld).
  2. Mô hình phần tử tập trung (Lumped Parameter Model): Chưa tích hợp mô hình phần tử hữu hạn (FEM) liên tục để khảo sát biến dạng đàn hồi cục bộ 3 chiều của răng và trục khuỷu phức tạp.
  3. Độ nhạy số học đối với hệ siêu cứng (Stiff Systems): Phương pháp bắn đơn có thể gặp hiện tượng tràn số ma trận khi tích phân các hệ vi phân có độ cứng quá lớn hoặc hệ số suy giảm quá cao.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển phương pháp bắn đa điểm (multiple shooting method) và bắn song song trên máy tính hiệu năng cao (HPC) cho các hệ cơ điện tử siêu cứng.
  • Nghiên cứu hiện tượng Sommerfeld khi động cơ có công suất hữu hạn kết hợp với bộ truyền bánh răng phi tuyến.
  • Mở rộng phân tích rẽ nhánh ngẫu nhiên (stochastic bifurcation) khi các thông số độ cứng và tải trọng chịu kích động ngẫu nhiên màu (colored noise).
  • Ứng dụng trí tuệ nhân tạo (Physics-Informed Neural Networks - PINNs) để học ma trận đơn đạo và dự báo rẽ nhánh thời gian thực.

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

  • Ảnh hưởng học thuật: Luận án là tài liệu tham khảo mẫu mực trong đào tạo sau đại học ngành Cơ học kỹ thuật và Cơ kỹ thuật tại Việt Nam, đóng góp nền tảng cho nhiều bài báo công bố trên các tạp chí chuyên ngành cơ học uy tín.
  • Chuyển đổi công nghiệp: Các thuật toán được chuyển giao ứng dụng trong thiết kế hệ trục truyền động tàu thuỷ tại Viện Khoa học Công nghệ Tàu thuỷ và các nhà máy đóng tàu thuộc Tổng công ty Công nghiệp Tàu thuỷ Việt Nam.
  • Lợi ích kinh tế - kỹ thuật: Giúp các kỹ sư thiết kế tối ưu hoá khe hở lắp ghép bánh răng, giảm rung chấn công nghiệp từ 15-25 dB, kéo dài tuổi thọ cụm trục khuỷu - thanh truyền máy tàu thuỷ và hạn chế nguy cơ gãy trục do mỏi xoắn cộng hưởng.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận phương pháp luận hoàn chỉnh về tính toán ma trận đơn đạo, lý thuyết Floquet và phương pháp bắn đơn cho các hệ cơ học phi tuyến.
  • Kỹ sư thiết kế cơ khí & R&D: Sở hữu công cụ số chính xác để dự đoán ứng xử động học, kiểm tra điều kiện bền mỏi và tối ưu hoá hệ thống dẫn động cơ khí trước khi chế tạo mẫu thử.
  • Chuyên gia chẩn đoán kỹ thuật: Sử dụng các đặc trưng phổ tần số (vạch chính $n\omega_z$, vải phụ $\omega_1$) để xây dựng thuật toán chẩn đoán lỗi mòn, tróc rỗ bề mặt răng và nứt trục tự độ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ì?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng thành công lý thuyết Floquet để xây dựng thuật toán số đồng thời xác định điều kiện đầu của nghiệm tuần hoàn và tính toán ma trận đơn đạo $C = e^{BT}$ cho các hệ cơ học có ma trận khối lượng biến thiên tuần hoàn theo thời gian $M(t)$.

2. Đổi mới phương pháp luận nổi bật khi so sánh với các nghiên cứu quốc tế?

So với phương pháp cân bằng điều hoà của Kahraman & Singh (1990) và phương pháp giải tích tiệm cận của Nayfeh (1979), phương pháp bắn đơn kết hợp tích phân Runge-Kutta-Nyström của luận án có khả năng giải quyết trực tiếp các hệ phi tuyến mạnh, phi tuyến gián đoạn (khe hở, tuyến tính từng khúc) mà không bị giới hạn bởi giả thiết tham số bé, đồng thời cho phép theo vết chính xác các điểm rẽ nhánh qua phổ nhân tử Floquet.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì?

Tính toán ma trận đơn đạo của bộ truyền bánh răng cho giá trị $|\lambda_{\max}| = 0.012662$, chứng minh rằng hệ dao động tham số thực tế có biên độ ổn định tiệm cận rất cao mặc dù độ cứng ăn khớp $k(t)$ biến thiên tuần hoàn rất phức tạp tới 10 bậc điều hoà.

4. Luận án có cung cấp đầy đủ giao thức để tái lặp nghiên cứu (Replication Protocol)?

Luận án cung cấp đầy đủ và chi tiết toàn bộ ma trận hệ số, bảng số liệu chuỗi Fourier (Bảng 1.2 gồm 10 hài độ cứng, Bảng 1.3 gồm 6 hài kích động, Bảng 1.4 gồm 12 hài hàm truyền cam, Bảng 1.5 gồm 48 hài mômen động cơ tàu thuỷ) cùng sơ đồ khối thuật toán Runge-Kutta-Nyström, đảm bảo khả năng tái lập 100% kết quả mô phỏng trên các phần mềm tính toán như MATLAB hoặc Fortran.

5. Luận án gợi mở lộ trình nghiên cứu 10 năm như thế nào?

Định hình lộ trình phát triển từ tính toán dao động tham số tập trung sang mô hình liên tục (FEM-Nonlinear Dynamics), tích hợp điều khiển chủ động triệt tiêu rẽ nhánh hỗn độn, và chuyển đổi từ phân tích xác định (deterministic) sang phân tích động lực học ngẫu nhiên (stochastic dynamics).

Kết luận

  1. Xây dựng thành công thuật toán số tích hợp phương pháp Runge-Kutta-Nyström và lý thuyết Floquet, giải quyết triệt để bài toán tìm điều kiện đầu và xác định tính ổn định Liapunov của các hệ phương trình vi phân tuyến tính hệ số tuần hoàn.
  2. Thiết lập quy trình số hoá bằng phương pháp bắn đơn để tính toán nghiệm tuần hoàn và khảo sát các dạng rẽ nhánh phức tạp (nút-yên ngựa, Hopf, nhân đôi chu kỳ, phá vỡ tính đối xứng) trong các hệ cơ học phi tuyến mạnh.
  3. Ứng dụng giải quyết trọn vẹn bài toán động lực học cho 3 mô hình kỹ thuật thực tế phức tạp: bộ truyền bánh răng một cấp, bộ điều chỉnh cam máy rèn dập, và dao động xoắn hệ trục động cơ diesel tàu thuỷ DTSC50.
  4. Phát hiện và chứng minh cơ chế hình thành dải phổ phụ (sidebands) trong phổ tần số dao động bánh răng, tạo tiền đề vững chắc cho kỹ thuật chẩn đoán hư hại cơ khí chính xác.
  5. Đánh dấu bước chuyển dịch mô hình quan trọng trong nghiên cứu cơ học máy tại Việt Nam: chuyển từ phân tích giải tích xấp xỉ hệ phi tuyến yếu sang mô phỏng số học chính xác động lực học phi tuyến mạnh và rẽ nhánh.
  6. Mở ra 3 hướng nghiên cứu mới: tính toán động lực học phi tuyến hệ nhiều vật có nguồn công suất hữu hạn, phân tích rẽ nhánh đa nghiệm trong hệ cơ học gián đoạn, và ứng dụng công nghệ số hoá chẩn đoán dao động thời gian thực trong công nghiệp chế tạo máy.