Tổng quan luận án

Luận án tiến sĩ kỹ thuật với đề tài "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 nghiên cứu sinh 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, bảo vệ năm 2011 tại Trường Đại học Bách khoa Hà Nội, thuộc chuyên ngành Cơ học kỹ thuật (mã số: 62.01).

Tính cấp thiết và khoảng trống nghiên cứu

Xu hướng phát triển của kỹ thuật chế tạo máy hiện đại đòi hỏi tốc độ làm việc của máy móc ngày càng cao trong khi trọng lượng kết cấu cần được tối ưu giảm thiểu. Điều này dẫn tới sự xuất hiện của các hiện tượng dao động phức tạp. Khi máy vận hành ở chế độ chuyển động bình ổn, các mô hình dao động thực tế thường được mô tả bởi các hệ phương trình vi phân phi tuyến hoặc các phương trình vi phân tuyến tính có hệ số biến đổi tuần hoàn theo thời gian. Mặc dù các lý thuyết giải tích cổ điển đã giải quyết tốt các bài toán phi tuyến yếu hoặc dao động nhỏ, nhưng khi gặp phải các hệ cơ học phi tuyến mạnh, hệ có độ cứng ăn khớp phi tuyến, có khe hở cơ học hoặc hệ có nhiều bậc tự do, việc tìm nghiệm giải tích gặp nhiều rào cản. Khoảng trống nghiên cứu then chốt được luận án xác định là: nhu cầu ứng dụng các lý thuyết toán học hiện đại về động lực học phi tuyến (lý thuyết rẽ nhánh, lý thuyết ổn định Liapunov, lý thuyết Floquet) và phát triển các thuật toán số hiệu quả để tính toán chính xác dao động tuần hoàn, xác định điều kiện đầu của nghiệm tuần hoàn và khảo sát các dạng rẽ nhánh trong các hệ dao động máy thực tế.

Mục tiêu và nhiệm vụ nghiên cứu

Mục tiêu tổng quát của luận án là xây dựng, hoàn thiện và áp dụng các thuật toán số cùng phương pháp giải tích gần đúng để xác định nghiệm dao động tuần hoàn, đánh giá tính ổn định động lực học và phân tích hiện tượng rẽ nhánh trong các mô hình dao động máy điển hình.

Các nhiệm vụ nghiên cứu cụ thể được đánh số và triển khai tuần tự:

  1. Xây dựng thuật toán số dựa trên phương pháp Runge-Kutta bậc 4 kết hợp lý thuyết Floquet nhằm xác định điều kiện đầu, tìm nghiệm tuần hoàn và khảo sát độ ổn định của các hệ phương trình vi phân tuyến tính hệ số tuần hoàn (hệ tham số tuyến tính).
  2. Xây dựng thuật toán số dựa trên phương pháp bắn đơn (Shooting method) để giải bài toán điều kiện biên, tìm điều kiện đầu và xác định nghiệm tuần hoàn cho các hệ cơ học phi tuyến mạnh (cả hệ ôtônôm và không ôtônôm).
  3. Ứng dụng phương pháp trung bình hóa đưa hệ về dạng chuẩn Lagrange-Bogoliubov để tính toán nghiệm tuần hoàn và phân tích các dạng rẽ nhánh cục bộ của điểm cố định (rẽ nhánh nút - yên ngựa, rẽ nhánh Hopf loại 2) trong các hệ phi tuyến yếu.
  4. Xây dựng thuật toán số phân tích rẽ nhánh của nghiệm tuần hoàn dựa trên phương pháp bắn đơn thông qua sự di chuyển của các nhân tử Floquet qua vòng tròn đơn vị trong mặt phẳng phức, áp dụng cho các hệ cơ học phi tuyến mạnh và hệ tuyến tính từng khúc.

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

  • Đối tượng nghiên cứu: Các mô hình động lực học và dao động trong máy, bao gồm:
    • Hệ tham số tuyến tính: Bộ truyền bánh răng một cấp có độ cứng ăn khớp biến đổi tuần hoàn; Bộ điều chỉnh sử dụng cơ cấu cam trong máy rèn dập; Hệ trục máy tàu thủy dẫn động bởi động cơ diesel 4 xi lanh (DTSC50).
    • Hệ phi tuyến và phi tuyến mạnh: Hệ dao động Duffing, hệ Van der Pol, hệ Van der Pol-Duffing, hệ Mathieu, hệ Rôto - Móng máy, hệ cơ học tuyến tính từng khúc và bộ truyền bánh răng một cấp có xét đến ảnh hưởng của khe hở ăn khớp.
  • Phạm vi nghiên cứu: Khảo sát trạng thái chuyển động bình ổn của cơ hệ; xét các chuyển động dao động tuần hoàn quanh chuyển động cơ bản; phạm vi không gian và thời gian gắn liền với các thông số kết cấu và chế độ vận hành thực tế của các thiết bị cơ khí được khảo sát trong kỹ thuật.

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

Các hướng nghiên cứu trước được luận án điểm lại

Luận án đã hệ thống hóa lịch sử phát triển của lý thuyết dao động phi tuyến và các phương pháp tính toán động lực học qua hai giai đoạn chính:

  1. Giai đoạn nửa đầu thế kỷ 20: Tập trung vào các bài toán dao động và ổn định của các hệ phi tuyến yếu bằng các phương pháp giải tích gần đúng cổ điển. Tác giả điểm lại các công trình nền tảng của:

    • N. M. Krylov và N. N. Bogoliubov (1937, 1947) về phương pháp tiệm cận và phương pháp trung bình hóa;
    • I. G. Malkin (1956) về các phương pháp giải tích trong lý thuyết dao động phi tuyến;
    • A. M. Liapunov (1950) về bài toán tổng quát của sự ổn định chuyển động;
    • N. V. Butenin (1976), J. J. Stoker (1950), N. Minorsky (1962), A. H. Nayfeh và D. T. Mook (1979) về lý thuyết dao động phi tuyến cổ điển.

  2. Giai đoạn từ cuối thế kỷ 20 đến thời điểm thực hiện luận án (khoảng 30 năm): Sự phát triển của khoa học máy tính đã thúc đẩy mạnh mẽ các nghiên cứu về động lực học phi tuyến mạnh, lý thuyết rẽ nhánh, đa tạp trung tâm và chuyển động hỗn độn. Các công trình tiêu biểu được tác giả trích dẫn gồm:

    • V. I. Arnold (1974, 1983), J. Guckenheimer và P. Holmes (1983), S. Wiggins (1990), R. Seydel (1994), Yu. A. Kuznetsov (1998) về lý thuyết rẽ nhánh và hệ động lực;
    • Nguyễn Văn Khang (1997, 2004, 2007) về động lực học hệ nhiều vật và dao động phi tuyến;
    • Các phương pháp số tìm nghiệm tuần hoàn cho hệ phi tuyến mạnh như: Phương pháp bắn đơn và bắn nhiều điểm (H. B. Keller 1968; R. Seydel 1994; Stoer và Bulirsch 2002); Phương pháp cân bằng điều hòa gia lượng (Incremental Harmonic Balance - S. L. Lau và Y. K. Cheung 1981); Phương pháp ánh xạ ô (Cell mapping - C. S. Hsu 1987); Phương pháp sai phân hữu hạn và phép chiếu (U. Ascher, R. M. M. Mattheij, R. D. Russell 1995).
    • Các nghiên cứu chuyên sâu về mô hình bánh răng và hệ cơ khí cụ thể: Blankenship và Kahraman (1995), Theodossiades và Natsiavas (2000), Vaishya và Singh (2003) về dao động tham số và phi tuyến trong bộ truyền bánh răng; C. S. Choi (1999) về dao động cơ cấu cam.

Khoảng trống luận án lựa chọn giải quyết

Mặc dù lý thuyết toán học về hệ động lực phi tuyến đã đạt nhiều thành tựu, nhưng việc chuyển hóa các kết quả lý thuyết này thành các quy trình tính toán số cụ thể, áp dụng trực tiếp cho các mô hình máy móc phức tạp trong kỹ thuật tại Việt Nam còn nhiều hạn chế. Luận án chọn giải quyết khoảng trống này bằng cách kết hợp chặt chẽ giữa mô hình hóa cơ học, thuật toán số hóa tìm điều kiện đầu nghiệm tuần hoàn (phương pháp Runge-Kutta bậc 4, phương pháp bắn) và phân tích rẽ nhánh động lực học (bằng phương pháp giải tích trung bình hóa cho hệ phi tuyến yếu và phương pháp số dựa trên nhân tử Floquet cho hệ phi tuyến mạnh).

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

Khung lý thuyết và khái niệm cốt lõi

  1. Lý thuyết ổn định Liapunov: Xét hệ phương trình trạng thái $\dot{y} = f(t, y)$ với $y \in \mathbb{R}^n$. Nghiệm $y^* = \varphi(t)$ được gọi là ổn định theo nghĩa Liapunov nếu với mọi $\varepsilon > 0$, tồn tại $\delta(\varepsilon, t_0) > 0$ sao cho mọi nghiệm $y(t)$ thỏa mãn $|y(t_0) - \varphi(t_0)| < \delta$ sẽ dẫn tới $|y(t) - \varphi(t)| < \varepsilon$ với mọi $t \ge t_0$. Nghiệm ổn định tiệm cận nếu thỏa mãn thêm điều kiện $\lim_{t \to \infty} |y(t) - \varphi(t)| = 0$.
  2. Lý thuyết Floquet cho hệ tuyến tính hệ số tuần hoàn: Xét hệ thuần nhất $\dot{x} = P(t)x$ với $P(t+T) = P(t)$. Ma trận cơ bản của hệ có dạng $\Phi(t) = A(t)e^{Bt}$, trong đó $A(t+T) = A(t)$, $B$ là ma trận hằng số. Ma trận $C = e^{BT} = \Phi(T)\Phi^{-1}(0)$ được gọi là ma trận đơn đạo (Monodromy matrix). Phương trình đặc trưng là $\det(C - \lambda E) = 0$.
    • Các nghiệm $\lambda_k$ ($k=1, \dots, n$) là các nhân tử đặc trưng (nhân tử Floquet).
    • Các số $\rho_k = \frac{1}{T}\ln \lambda_k$ là các số mũ đặc trưng.
    • Tiêu chuẩn ổn định Floquet: Hệ ổn định tiệm cận khi và chỉ khi tất cả các nhân tử Floquet nằm trong vòng tròn đơn vị trên mặt phẳng phức, tức $|\lambda_k| < 1$ (hoặc $\text{Re}(\rho_k) < 0$); hệ không ổn định nếu tồn tại ít nhất một $|\lambda_k| > 1$.
  3. Lý thuyết rẽ nhánh: Khảo sát sự thay đổi định tính về cấu trúc tôpô của nghiệm (điểm cân bằng hoặc chu trình giới hạn) khi một hoặc nhiều tham số điều khiển biến thiên. Bao gồm các dạng rẽ nhánh điểm cố định (nút - yên ngựa, chuyển qua tới hạn, hình dĩa, Hopf) và các dạng rẽ nhánh của nghiệm tuần hoàn (phá vỡ tính đối xứng, nếp gấp - chu trình, nhân đôi chu kỳ, rẽ nhánh Hopf loại 2 / Neimark-Sacker).
          [Khung phân tích động lực học của luận án]
[Hệ tham số tuyến tính]                     [Hệ phi tuyến mạnh]
  - Khai triển Taylor quanh q_R(t)            - Phương pháp bắn đơn (Shooting)
  - Tìm x(0) qua Runge-Kutta 4                - Ma trận Jacobi sai phân
  - Tích phân Runge-Kutta-Nyström             - Tính nhân tử Floquet từ Monodromy
  - Tính ma trận đơn đạo & Floquet            - Khảo sát rẽ nhánh (vượt vòng đơn vị)

Phương pháp nghiên cứu và công cụ tính toán

  • Phương pháp số tìm điều kiện đầu nghiệm tuần hoàn của hệ tham số tuyến tính: Chia chu kỳ $T$ thành $m$ khoảng bằng nhau với bước $h = T/m$. Ứng dụng công thức Runge-Kutta bậc 4 để thiết lập hệ thức đại số tuyến tính liên kết giữa các trạng thái $x_i = A_{i-1}x_{i-1} + b_{i-1}$. Kết hợp điều kiện tuần hoàn $x(0) = x(T)$, hệ được giải trực tiếp để tìm vector điều kiện đầu $x_0 = x(0)$ mà không cần lặp qua nhiều chu kỳ.
  • Phương pháp Runge-Kutta-Nyström: Áp dụng tích phân số trực tiếp cho hệ phương trình vi phân cấp hai dạng $\ddot{x} = f(t, x, \dot{x})$ sau khi đã biết điều kiện đầu $x(0), \dot{x}(0)$, giúp tăng tốc độ tính toán và giảm thiểu sai số tích lũy so với việc chuyển về hệ cấp một.
  • Phương pháp lặp Von-Mises: Dùng để xác định trị riêng có môđun lớn nhất $\lambda_{\max}$ của ma trận đơn đạo $C$, phục vụ đánh giá nhanh điều kiện ổn định.
  • Phương pháp bắn đơn (Shooting method): Chuyển bài toán biên hai điểm tìm nghiệm tuần hoàn chu kỳ $T$ ($x(T, x_0) - x_0 = 0$) thành bài toán tìm nghiệm của hệ phương trình phi tuyến đối với giá trị đầu $x_0$ bằng thuật toán Newton-Raphson kết hợp tích phân số xác định ma trận Jacobi.
  • Phương pháp trung bình hóa: Đưa phương trình vi phân phi tuyến về dạng chuẩn Lagrange-Bogoliubov, tính giá trị trung bình theo thời gian để khử biến nhanh, từ đó khảo sát tính ổn định và rẽ nhánh của các nghiệm dừng.

Nguồn dữ liệu và tư liệu

Dữ liệu sử dụng để tính toán kiểm chứng số học trong luận án được thu thập từ các công trình nghiên cứu cơ học công bố quốc tế và số liệu thiết kế thực tế của các động cơ, thiết bị cơ khí (bộ truyền bánh răng thí nghiệm theo tài liệu của Vaishya và Singh; cơ cấu cam máy rèn dập theo tài liệu của C. S. Choi và Nguyễn Văn Khang; hệ trục tàu thủy trang bị động cơ DTSC50 theo tài liệu của Bộ môn Cơ học Trường Đại học Hàng hải Việt Nam).

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

Chương 1: Tính toán dao động tuần hoàn và ổn định động lực của một số hệ cơ học tuyến tính hệ số tuần hoàn

Chương 1 thiết lập cơ sở toán học cho phương pháp số tìm điều kiện đầu nghiệm tuần hoàn và tính toán ổn định Floquet, sau đó áp dụng chi tiết vào 3 mô hình cơ học:

1. Bộ truyền bánh răng một cấp

Mô hình gồm 2 bánh răng quay với vận tốc góc không đổi $\omega_1, \omega_2$, kể đến biến thiên độ cứng ăn khớp $k(t)$ và sai số bước răng $e(t)$. Phương trình vi phân dao động quy về biến sai số truyền động $q(t)$: $$m\ddot{q} + c\dot{q} + k(t)q = k_0 q_0 + [k(t) - k_0]e(t) + c\dot{e}(t)$$ Trong đó các thông số tính toán gồm: mômen quán tính $J_1 = 9{,}3 \cdot 10^{-2}\text{ kg}\cdot\text{m}^2$, $J_2 = 0{,}272\text{ kg}\cdot\text{m}^2$; bán kính cơ sở $r_1 = 30{,}46 \cdot 10^{-3}\text{ m}$; số răng $z_1 = 14$, $z_2 = 39$; vận tốc góc $\omega_1 = 60\text{ rad/s}$; sai số tĩnh $q_0 = 1{,}2 \cdot 10^{-5}\text{ m}$; tần số riêng không cản trung bình $\omega_0 = 5462\text{ s}^{-1}$ ($f_0 = 869{,}3\text{ Hz}$); độ cản Lehr $D = 0{,}024$.

Bảng 1: Các hệ số Fourier của độ cứng ăn khớp $k(t)$ (10 hài) và hàm kích động $e(t)$ (6 hài)

$n$ $k_n$ (N/m) $\psi_n$ (rad) $i$ $e_i$ (mm) $\alpha_i$ (rad)
0 $8{,}1846 \cdot 10^6$ — 1 0,0015 -0,0490
1 $1{,}5286 \cdot 10^6$ -0,8696 2 0,0035 -1,7661
2 $0{,}8369 \cdot 10^6$ 1,2280 3 0,0027 -0,7286
3 $0{,}5142 \cdot 10^6$ -2,8931 4 0,0011 -0,5763
4 $0{,}3645 \cdot 10^6$ -0,7539 5 0,0005 -0,7810
5 $0{,}2756 \cdot 10^6$ 1,4025 6 0,0013 1,8172

Kết quả tính toán ma trận đơn đạo cho: $$C = \begin{bmatrix} 0{,}01234 & 0{,}519742 \cdot 10^{-6} \ -15{,}121177 & 0{,}012356 \end{bmatrix}$$ Các nhân tử Floquet: $\lambda_{1,2} = 0{,}012348 \pm 0{,}0028034i \implies |\lambda_{\max}| = 0{,}012662 < 1$. Nghiệm tuần hoàn ổn định tiệm cận. Phổ tần số của sai số động lực chỉ ra các đỉnh tại tần số ăn khớp $n\omega_z$ và các dải điều hòa phụ cách nhau đúng bằng tần số quay $\omega_1$, hoàn toàn phù hợp với kết quả đo thực nghiệm.

2. Bộ điều chỉnh bằng cơ cấu cam trong máy rèn dập

Hệ 3 bậc tự do $(\varphi_0, \varphi_1, q_2)$ được tuyến tính hóa quanh chuyển động cơ bản $\varphi_0 = \omega t$, đưa về hệ 2 bậc tự do $q = [q_1, q_2]^T$ dạng ma trận: $M(t)\ddot{q} + C(t)\dot{q} + K(t)q = d(t)$. Thông số: $n = 50\text{ vòng/phút}$ ($\omega = 5{,}236\text{ rad/s}$), $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/rad}$. Hàm truyền bậc nhất $U(\omega t)$ xấp xỉ bằng chuỗi Fourier 12 số hạng với 2 trường hợp hệ số $a_k$.

Bảng 2: Điều kiện đầu và nhân tử Floquet của bộ điều chỉnh cơ cấu cam | Trường hợp | Điều kiện đầu $q(0)$ | Điều kiện đầu $\dot{q}(0)$ | $|\lambda_{\max}|$ | Trạng thái ổn định | | :--- | :--- | :--- | :---: | :---: | | Trường hợp 1 | $[-0{,}002897; 0{,}808339 \cdot 10^{-5}]^T$ | $[-0{,}555743; 0{,}015551]^T$ | 0,001997 | Ổn định tiệm cận | | Trường hợp 2 | $[-0{,}043726; 0{,}001019]^T$ | $[-0{,}280632; 0{,}008297]^T$ | 0,001454 | Ổn định tiệm cận |

3. Dao động xoắn hệ trục máy tàu thủy

Khảo sát động cơ 4 xi lanh DTSC50 kết nối bánh đà, hộp số và chân vịt (7 bậc tự do, chuyển thành hệ 6 tọa độ tương đối $q_1, \dots, q_6$). Tốc độ quay chân vịt $n = 1700\text{ vòng/phút}$ ($\omega = 178\text{ rad/s}$), độ cứng xoắn $k_1 = k_2 = k_3 = 1{,}2 \cdot 10^6\text{ Nm/rad}$, $k_5 = 0{,}1166 \cdot 10^6\text{ Nm/rad}$, hệ số cản $c_1 = \dots = c_4 = 1{,}484\text{ kg}\cdot\text{m}^2/\text{s}$, mômen xoắn trung bình $M_0 = 60{,}38\text{ Nm}$ với chuỗi Fourier 48 điều hòa. Tìm được chính xác vector điều kiện đầu nghiệm tuần hoàn:

  • $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)}$

Chương 2: Tính toán dao động tuần hoàn của một số hệ cơ học phi tuyến mạnh bằng phương pháp bắn

Chương 2 trình bày thuật toán bắn đơn (Shooting method) để giải bài toán biên tìm nghiệm tuần hoàn chu kỳ $T$ của hệ phi tuyến:

  1. Hệ không ôtônôm: $\dot{x} = f(t, x)$ với $f(t+T, x) = f(t, x)$. Bài toán biên đặt ra là tìm giá trị đầu $x_0$ sao cho $F(x_0) = x(T, x_0) - x_0 = 0$. Thuật toán lặp Newton-Raphson được sử dụng: $$x_0^{(k+1)} = x_0^{(k)} - \left[ \frac{\partial x(T, x_0^{(k)})}{\partial x_0} - E \right]^{-1} \left( x(T, x_0^{(k)}) - x_0^{(k)} \right)$$ Trong đó ma trận đạo hàm riêng (ma trận Jacobi) được tính toán bằng cách tích phân đồng thời hệ phương trình biến phân hoặc dùng công thức sai phân số.
  2. Hệ ôtônôm: $\dot{x} = f(x)$, chu kỳ dao động $T$ chưa biết trước và trở thành một ẩn số cần tìm cùng với vector điều kiện đầu $x_0$. Luận án bổ sung điều kiện pha (chẳng hạn cố định một thành phần vận tốc ban đầu bằng 0) để khép kín hệ phương trình.
  3. Ứng dụng tính toán cụ thể:
    • Xác định nghiệm tuần hoàn và chu trình giới hạn của hệ Van der Pol một bậc tự do và hệ dao động Duffing.
    • Tính toán dao động tuần hoàn của hệ Rôto - Móng máy chịu lực kích động mất cân bằng.
    • Tính toán hệ cơ học có đặc tính đàn hồi tuyến tính từng khúc (Piecewise linear system).
    • Tính toán dao động tuần hoàn trong bộ truyền bánh răng một cấp có xét đến độ cứng ăn khớp thay đổi theo thời gian và khe hở ăn khớp giữa các răng (mô hình phi tuyến gián đoạn).

Chương 3: Tính toán và khảo sát rẽ nhánh dao động tuần hoàn của một số hệ cơ học phi tuyến bằng phương pháp trung bình hoá

Chương 3 áp dụng lý thuyết giải tích gần đúng và lý thuyết rẽ nhánh cục bộ của điểm cố định:

  • Phương pháp trung bình hóa: Chuyển hệ phương trình phi tuyến về dạng chuẩn Lagrange-Bogoliubov: $\dot{y} = \varepsilon f(t, y)$. Áp dụng toán tử trung bình hóa theo thời gian đối với trường hợp $f(t,y)$ tuần hoàn và trường hợp tổng quát để thu được hệ phương trình tự trị trung bình $\dot{\xi} = \varepsilon f_0(\xi)$.
  • Phân loại rẽ nhánh một thứ nguyên và rẽ nhánh Hopf: Khảo sát các điểm cân bằng của hệ trung bình (tương ứng với nghiệm tuần hoàn của hệ ban đầu) và tính chất ổn định thông qua các giá trị riêng của ma trận Jacobi:
    • Rẽ nhánh nút - yên ngựa (Saddle-node bifurcation): Một giá trị riêng thực đi qua gốc tọa độ ($\lambda = 0$).
    • Rẽ nhánh chuyển qua tới hạn (Transcritical bifurcation) và rẽ nhánh hình dĩa (Pitchfork bifurcation).
    • Rẽ nhánh Hopf (Hopf bifurcation): Một cặp giá trị riêng thuần ảo vượt qua trục ảo ($\text{Re}(\lambda_{1,2}) = 0$), dẫn đến sự sinh ra hoặc mất đi của một chu trình giới hạn (tương ứng với dao động hai tần số / nghiệm gần tuần hoàn trong hệ gốc).
  • Khảo sát rẽ nhánh trên các hệ dao động cụ thể:
    • Rẽ nhánh nút - yên ngựa trong hệ dao động Duffing chịu kích động điều hòa.
    • Rẽ nhánh nút - yên ngựa trong hệ phương trình tham số Mathieu phi tuyến.
    • Rẽ nhánh nút - yên ngựa và rẽ nhánh Hopf loại 2 trong hệ Van der Pol chịu lực kích động cưỡng bức ngoài.
    • Rẽ nhánh nút - yên ngựa và rẽ nhánh Hopf loại 2 trong hệ tự chấn Van der Pol-Duffing kết hợp.

Chương 4: Tính toán và khảo sát rẽ nhánh dao động tuần hoàn trong các hệ cơ học phi tuyến bằng phương pháp bắn

Chương 4 mở rộng phương pháp bắn đơn để khảo sát rẽ nhánh của nghiệm tuần hoàn trong các hệ phi tuyến mạnh thông qua cơ chế toán học của các nhân tử Floquet:

                  [Cơ chế rẽ nhánh qua nhân tử Floquet]
[Rẽ nhánh nhân đôi chu kỳ]   [Rẽ nhánh nút-yên ngựa / Gập]   [Rẽ nhánh Hopf loại 2 / Neimark]
   λ = -1 vượt vòng tròn        λ = +1 vượt vòng tròn        Cặp λ = e^(±iθ) rời vòng tròn
   (Xuất hiện nghiệm 2T)        (Nghiệm biến mất/hợp nhất)     (Xuất hiện chuyển động xuyến)
  1. Lý thuyết rẽ nhánh của nghiệm tuần hoàn:
    • Rẽ nhánh phá vỡ tính đối xứng (Symmetry-breaking bifurcation): Nghiệm tuần hoàn mất tính đối xứng ban đầu.
    • Rẽ nhánh nếp gấp - chu trình (Fold-cycle / Limit point bifurcation): Một nhân tử Floquet thực đi ra khỏi vòng tròn đơn vị tại điểm $+1$ ($\lambda = +1$).
    • Rẽ nhánh chuyển qua giới hạn (Transcritical bifurcation của chu trình).
    • Rẽ nhánh nhân đôi chu kỳ (Period-doubling / Flip bifurcation): Một nhân tử Floquet đi ra khỏi vòng tròn đơn vị tại điểm $-1$ ($\lambda = -1$), chu kỳ của dao động chuyển từ $T$ thành $2T, 4T, \dots$ (con đường dẫn tới hỗn độn Feigenbaum).
    • Rẽ nhánh Hopf loại hai (Neimark-Sacker bifurcation): Một cặp nhân tử Floquet phức liên hợp rời khỏi vòng tròn đơn vị $|\lambda_{1,2}| = 1$ ($\lambda_{1,2} = e^{\pm i\theta}$), dẫn tới sự xuất hiện của mặt xuyến 2 chiều (chuyển động gần tuần hoàn hai tần số).
  2. Thuật toán số phân tích rẽ nhánh: Sử dụng phương pháp bắn đơn kết hợp kỹ thuật kéo dài tham số (Path following/continuation) để vẽ đường cong đáp ứng biên độ - tần số và tính toán vết của ma trận đơn đạo dọc theo nhánh nghiệm.
  3. Các kết quả tính toán số học tiêu biểu:
    • Phân tích rẽ nhánh trong hệ Van der Pol-Duffing phi tuyến mạnh.
    • Phân tích rẽ nhánh nhân đôi chu kỳ trong hệ cơ học có độ cứng tuyến tính từng khúc.
    • Khảo sát biểu đồ rẽ nhánh và sự biến đổi của nhân tử Floquet trong bộ truyền bánh răng một cấp khi có kể đến khe hở ăn khớp, phát hiện các vùng mất ổn định và quá trình phân nhánh phức tạp của sai số truyền động khi vận tốc góc thay đổi.

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

Theo công bố của tác giả trong văn bản luận án, các kết quả và đóng góp mới bao gồm:

Đóng góp về mặt lý luận và phương pháp

  1. Xây dựng và hoàn thiện thuật toán số dựa trên phương pháp Runge-Kutta bậc 4 để tìm chính xác điều kiện đầu của nghiệm tuần hoàn trong hệ phương trình vi phân tuyến tính hệ số tuần hoàn, khắc phục nhược điểm tích phân lặp tốn thời gian và sai số lũy tích.
  2. Thiết lập quy trình số kết hợp phương pháp bắn đơn với lý thuyết Floquet để tính toán điều kiện đầu, xác định chu trình dao động tuần hoàn và phân tích rẽ nhánh cho các hệ cơ học phi tuyến mạnh, hệ ôtônôm và hệ có độ cứng gián đoạn từng khúc.
  3. Hệ thống hóa và làm sáng tỏ mối liên hệ giữa lý thuyết rẽ nhánh điểm cố định (thông qua phương pháp trung bình hóa) và lý thuyết rẽ nhánh của nghiệm tuần hoàn (thông qua chuyển dịch nhân tử Floquet trên vòng tròn đơn vị) trong các bài toán cơ kỹ thuật.

Đóng góp về mặt thực tiễn kỹ thuật

  1. Ứng dụng thành công các thuật toán số đã xây dựng vào việc tính toán dao động và kiểm tra ổn định động lực cho các hệ thống kỹ thuật thực tế: bộ truyền bánh răng một cấp, bộ điều chỉnh cơ cấu cam trong máy rèn dập và hệ trục xoắn máy tàu thủy 4 xi lanh (DTSC50).
  2. Làm rõ cơ chế sinh dao động trong bộ truyền bánh răng: phổ dao động gồm các hài chính của tần số ăn khớp $n\omega_z$ và các dải điều hòa phụ cách nhau đúng bằng tần số quay $\omega_1$, cung cấp cơ sở định lượng để chẩn đoán, đánh giá và định vị hư hại của bánh răng dựa trên phân tích phổ rung động.
  3. Cung cấp bộ số liệu tính toán chi tiết về vector điều kiện đầu, ma trận đơn đạo và nhân tử Floquet cho các kỹ sư thiết kế máy nhằm kiểm tra độ ổn định động lực và tránh các vùng rẽ nhánh nguy hiểm.

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

Dựa trên các giả thiết và phạm vi mô hình hóa được tác giả trình bày trong luận án:

  • Hạn chế:
    • Các bài toán dao động tham số mới chỉ khảo sát ở chế độ chuyển động bình ổn với giả thiết vận tốc góc của khâu dẫn hoặc trục máy là hằng số tuyệt đối ($\omega = \text{const}$), chưa xét đến tương tác phi tuyến ngược chiều giữa dao động và nguồn động lực (hiện tượng Sommerfeld).
    • Khai triển Taylor tuyến tính hóa trong Chương 1 bỏ qua các số hạng bậc cao, chỉ áp dụng chính xác cho lân cận nhỏ của chuyển động cơ bản.
    • Hàm truyền của cam, mômen xoắn và độ cứng ăn khớp được xấp xỉ bằng chuỗi Fourier với số phần tử hữu hạn ($K=12$, $N=10$, $m=48$), có thể bỏ sót các thành phần tần số rất cao.
  • Hướng nghiên cứu tiếp:
    • Mở rộng thuật toán bắn đơn và phân tích rẽ nhánh cho các hệ cơ học phi tuyến nhiều bậc tự do không bình ổn (vận tốc góc khâu dẫn biến thiên).
    • Khảo sát sâu hơn các trạng thái chuyển động hỗn độn (Chaos) sau chuỗi rẽ nhánh nhân đôi chu kỳ trong các hệ cơ cấu có khe hở và va đập.

Giá trị tham khảo

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ học kỹ thuật, Kỹ thuật cơ khí: Tham khảo phương pháp thiết lập phương trình Lagrange loại 2 cho hệ nhiều vật đàn hồi, phương pháp chuyển đổi tọa độ tương đối và quy trình lập trình các thuật toán số (Runge-Kutta-Nyström, Shooting method, thuật toán Floquet).
  • Giảng viên đại học: Sử dụng làm tài liệu giảng dạy các chuyên đề cao học về Dao động phi tuyến, Động lực học máy và Hệ động lực.
  • Kỹ sư tính toán thiết kế máy, kỹ sư tàu thủy: Ứng dụng các công thức và bảng số liệu phân tích dao động xoắn hệ trục động cơ nhiều xi lanh, dao động cơ cấu cam và dao động ăn khớp bánh răng để tối ưu hóa kết cấu, hạn chế tiếng ồn và phục vụ công tác chẩn đoán kỹ thuật qua phổ dao động.

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

1. Tại sao cần xác định điều kiện đầu của nghiệm tuần hoàn thay vì tích phân số thông thường từ một điều kiện đầu bất kỳ?

Trong các hệ thực tế, việc tích phân số từ một điều kiện ban đầu tùy ý đòi hỏi hệ phải trải qua giai đoạn dao động quá độ rất dài trước khi đạt tới trạng thái dừng tuần hoàn. Quá trình này làm tích lũy sai số tính toán của thuật toán số và tốn nhiều thời gian máy tính. Việc tìm trực tiếp vector điều kiện đầu $x(0)$ bằng thuật toán Runge-Kutta bậc 4 hoặc phương pháp bắn giúp xác định ngay lập tức nghiệm tuần hoàn chính xác trong đúng một chu kỳ $T$.

2. Tiêu chuẩn ổn định Floquet được áp dụng như thế nào trong luận án?

Tính ổn định của nghiệm tuần hoàn được xác định thông qua ma trận đơn đạo $C = \Phi(T)$. Bằng phương pháp số hoặc phép lặp Von-Mises, luận án tìm các giá trị riêng $\lambda_k$ (nhân tử Floquet). Nếu nhân tử có môđun lớn nhất thỏa mãn $|\lambda_{\max}| < 1$, nghiệm tuần hoàn chắc chắn ổn định tiệm cận; nếu $|\lambda_{\max}| > 1$, dao động sẽ mất ổn định.

3. Phổ dao động của sai số truyền động bánh răng trong Chương 1 phản ánh điều gì?

Phổ tần số chỉ ra hai thành phần đặc trưng:

  • Các đỉnh biên độ chính tại tần số ăn khớp $\omega_z$ và các bội số $n\omega_z$, phản ánh kích động tham số do độ cứng ăn khớp biến thiên $k(t)$.
  • Các dải điều hòa phụ bao quanh các đỉnh chính, cách nhau đúng một khoảng bằng tần số quay $\omega_1$ của trục bánh răng, phản ánh kích động sai số bước răng $e(t)$. Độ lớn của các dải phụ này là căn cứ để chẩn đoán mức độ hư hỏng bề mặt răng.

4. Rẽ nhánh nhân đôi chu kỳ trong Chương 4 xảy ra theo cơ chế toán học nào?

Rẽ nhánh nhân đôi chu kỳ xảy ra khi tham số điều khiển biến thiên làm cho một nhân tử Floquet thực của nghiệm tuần hoàn rời khỏi vòng tròn đơn vị qua điểm $-1$ ($\lambda = -1$). Tại thời điểm này, nghiệm chu kỳ $T$ ban đầu mất ổn định và sinh ra một nghiệm tuần hoàn mới có chu kỳ $2T$.

5. Sự khác biệt giữa phương pháp trung bình hóa (Chương 3) và phương pháp bắn đơn (Chương 4) khi khảo sát rẽ nhánh là gì?

Phương pháp trung bình hóa là phương pháp giải tích gần đúng, áp dụng cho các hệ dao động phi tuyến yếu có chứa tham số bé $\varepsilon$, đưa hệ về các điểm cố định dừng để khảo sát rẽ nhánh cục bộ. Ngược lại, phương pháp bắn đơn là phương pháp số thuần túy, áp dụng được cho các hệ phi tuyến mạnh và phi tuyến gián đoạn (tuyến tính từng khúc, có khe hở), khảo sát rẽ nhánh trực tiếp trên chu trình tuần hoàn thông qua các nhân tử Floquet.

Kết luận

Luận án của tác giả Hoàng Mạnh Cường đã giải quyết bài toán tính toán dao động tuần hoàn, đánh giá độ ổn định và khảo sát rẽ nhánh cho cả hai lớp mô hình: hệ tham số tuyến tính và hệ phi tuyến mạnh trong máy. Bằng việc kết hợp phương pháp Runge-Kutta bậc 4, phương pháp bắn đơn, lý thuyết Floquet và phương pháp trung bình hóa, công trình đã xây dựng hệ thống thuật toán số hoàn chỉnh và kiểm chứng thành công trên các mô hình kỹ thuật thực tế như bộ truyền bánh răng, cơ cấu cam máy rèn dập và hệ trục xoắn máy tàu thủy. Các kết quả đạt được cung cấp cơ sở khoa học và công cụ tính toán phục vụ thiết kế động lực học và chẩn đoán trạng thái máy trong ngành cơ khí.