Giới thiệu dự án

Trong kỷ nguyên Công nghiệp 4.0 và chuyển đổi số kỹ thuật, việc phát triển các hệ thống cơ điện tử (Mechatronics) đa miền năng lượng—kết hợp cơ học, điện từ học, thủy lực và điều khiển tự động—đang đối mặt với thách thức lớn về tính phức tạp phi tuyến và tương tác đa vật lý. Theo các báo cáo nghiên cứu từ Hiệp hội Kỹ sư Cơ khí Quốc tế (ASME), hơn 68% chi phí và thời gian phát triển các hệ thống tự động hóa công nghiệp phát sinh từ giai đoạn thử nghiệm phần cứng và khắc phục các sai số ghép nối liên ngành. Do đó, phương pháp Mô hình hóa và Mô phỏng dựa trên đồ hình liên kết (Bond Graph Modeling) kết hợp với công cụ mô phỏng chuyên dụng 20-sim là giải pháp then chốt giúp tối ưu hóa chu trình phát triển sản phẩm ảo (Virtual Prototyping).

Đồ án môn học "Mô Hình Hóa và Mô Phỏng Hệ Thống Cơ Điện Tử" do nhóm sinh viên Khóa K15, Khoa Cơ khí - Trường Đại học Công nghiệp Hà Nội (HaUI) thực hiện dưới sự hướng dẫn của TS. Lê Ngọc Duy và TS. Nguyễn Anh Tú, tập trung giải quyết bài toán cốt lõi: Thiết lập cơ sở toán học thống nhất, xây dựng đồ hình liên kết Bond Graph và tổng hợp bộ điều khiển cho 3 hệ thống cơ điện tử kinh điển:

  1. Hệ thống điều khiển vị trí động cơ điện một chiều nam châm vĩnh cửu (PMDC Motor).
  2. Hệ thống treo chủ động/bán chủ động xe bus 2 bậc tự do (Quarter-Car Suspension System).
  3. Hệ thống xe và con lắc ngược (Inverted Pendulum on Cart System).
       +-----------------------------------------------------------+
       |   Đồ án: Mô Hình Hóa & Mô Phỏng Hệ Thống Cơ Điện Tử       |
       |       Khoa Cơ Khí - Đại Học Công Nghiệp Hà Nội (HaUI)      |
       +-----------------------------------------------------------+
                                     |
         +---------------------------+---------------------------+
         |                           |                           |
         v                           v                           v
+------------------+       +------------------+       +-------------------+
|  1. Động Cơ DC   |       | 2. Hệ Treo Xe Bus|       |3. Con Lắc Ngược   |
|   Nam Châm VC    |       |   (Quarter-Car)  |       |   (Inverted Pend) |
+------------------+       +------------------+       +-------------------+
| * Mô hình hóa vật lý     | * Mô hình 2-DOF  |       | * Euler-Lagrange  |
| * Bond Graph + Causality | * Khử dao động   |       | * Tuyến tính hóa  |
| * Điều khiển P/PI/PD/PID | * Bộ điều khiển  |       | * Cân bằng vòng   |
| * Bám vị trí rad         |   PID giảm chấn  |         kín 20-sim        |
+------------------+       +------------------+       +-------------------+

Mục tiêu dự án

  1. Mô hình hóa giải tích: Thiết lập hệ phương trình vi phân mô tả động học và động lực học của 3 hệ thống thông qua định luật II Newton, định luật Kirchhoff và phương pháp giải tích năng lượng Euler-Lagrange.
  2. Ánh xạ Bond Graph đa miền: Chuyển đổi các phần tử vật lý (điện trở $R$, điện cảm $L$, khối lượng/quán tính $I$, độ cứng lò xo $C$, hệ số cản $R$, biến đổi cơ-điện $GY/TF$) sang cấu trúc 0-junction và 1-junction đồng nhất.
  3. Xác định quan hệ nhân quả (Causality Assignment): Thiết lập chuẩn hóa quan hệ vào/ra cho các phần tử tích trữ năng lượng (tích phân) và phần tử nguồn ($Se, Sf$).
  4. Tổng hợp bộ điều khiển PID: Thiết kế và tinh chỉnh thông số bộ điều khiển $P$, $PI$, $PD$, $PID$ nhằm triệt tiêu dao động, giảm thời gian đáp ứng quá độ và sai số xác lập.
  5. Mô phỏng và kiểm chứng trên 20-sim: Đánh giá định lượng tính ổn định, độ vọt lố ($POT < 5%$) và thời gian xác lập ($t_s < 2\text{s}$) dưới tác động của tín hiệu đặt và xung nhiễu tải.

Phạm vi và giới hạn

  • Phạm vi: Khảo sát động lực học tuyến tính hóa cục bộ quanh điểm cân bằng làm việc; mô phỏng đáp ứng miền thời gian trên phần mềm 20-sim.
  • Giới hạn: Chưa xét đến hiện tượng bão hòa từ trường phi tuyến trong lõi thép động cơ, phi tuyến hình học hành trình giảm chấn khi vượt giới hạn cơ khí, và góc lệch lớn ($\theta > 30^\circ$) của con lắc ngược chưa tích hợp bộ điều khiển phi tuyến toàn cục (Nonlinear / Sliding Mode Control).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trong kỹ thuật cơ điện tử truyền thống, kỹ sư thường mô hình hóa riêng rẽ từng miền vật lý rồi ghép nối bằng sơ đồ khối toán học (Block Diagram). Cách tiếp cận này bộc lộ nhiều nhược điểm khi xử lý các tương tác năng lượng hai chiều (Power Conjugate Variables: Effort $e$ và Flow $f$).

Tiêu chí so sánh Sơ đồ khối (Simulink / Xcos) Mã hóa phương trình (C++/Python ODE) Đồ hình Bond Graph (20-sim)
Bản chất mô hình Luồng tín hiệu một chiều (Signal-flow) Dạng giải tích số (State-space equations) Luồng công suất bảo toàn ($P = e \cdot f$)
Tính đa miền (Multi-physics) Ghép nối phức tạp, dễ sai dấu công suất Tự dẫn xuất thủ công, dễ nhầm lẫn biến Tự nhiên qua các cổng chuẩn hóa ($0/1\text{-Junction}$)
Kiểm soát tính nhân quả Cố định thủ công theo khối toán tử Không tự động phát hiện vòng lặp đại số Tự động gán nhân quả (SCAP), phát hiện suy biến
Khả năng mở rộng module Thấp (thay đổi cấu trúc phải vẽ lại khối) Trung bình (sửa hệ phương trình vi phân) Rất cao (kết nối trực tiếp qua Power Port)

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must have: Hàm truyền chuẩn xác cho 3 hệ thống; biểu đồ Bond Graph hoàn chỉnh gán đúng causality; bộ điều khiển PID đưa hệ thống về trạng thái ổn định.
  • Should have: Bảng so sánh chi tiết giữa các cấu hình P, PI, PD và PID trên cùng một đáp ứng bước (step response).
  • Could have: Tích hợp các bộ lọc nhiễu bậc cao và khảo sát tính bền vững khi tham số thay đổi $\pm 20%$.
  • Won't have: Thử nghiệm HIL (Hardware-in-the-Loop) trên mô hình phần cứng vật lý thời gian thực.

Thiết kế hệ thống

Kiến trúc mô hình hóa dựa trên cấu trúc biến liên hợp công suất (Power Bond) truyền tải đồng thời biến thế $e(t)$ (Điện áp $V$, Lực $F$, Momen $T$) và biến dòng $f(t)$ (Dòng điện $i$, Vận tốc $v$, Vận tốc góc $\omega$).

                         CẤU TRÚC BOND GRAPH CHUẨN HÓA
                
   [Nguồn Se/Sf]  --->  [0-Junction (e1=e2=e3; sum(f)=0)]  --->  [Phần tử tích C/I]
                                     |
                        [1-Junction (f1=f2=f3; sum(e)=0)]  --->  [Phần tử tiêu tán R]
                                     |
                        [Biến đổi GY / TF (Cơ <-> Điện)]

Ngăn xếp công nghệ và thông số chuẩn định

  • Môi trường mô phỏng: Controllab 20-sim Professional Simulation Software (Hỗ trợ ngôn ngữ SIDOPS+).
  • Thuật toán tích phân số: Runge-Kutta 4th Order (RK4) / Vode tích phân biến bước giải hệ phương trình vi phân đại số (DAE).
  • Thông số thực nghiệm trích xuất từ đề tài:
System_1_DC_Motor:
  Inductance_L: 0.002 # Henry (H)
  Resistance_R: 0.8   # Ohm
  Damping_b: 0.0076   # N.m.s/rad
  Inertia_J: 0.1      # N.m.s^2/rad
  Torque_Constant_K: 0.3 # N.m/A (V.s/rad)
  Input_Voltage: 24.0 # V

System_2_Bus_Suspension:
  Sprung_Mass_m1: 2500.0  # kg (Thân xe)
  Unsprung_Mass_m2: 320.0 # kg (Cụm bánh xe)
  Suspension_Stiffness_k1: 60000.0 # N/m
  Tire_Stiffness_k2: 400000.0      # N/m
  Suspension_Damping_b1: 350.0     # N.s/m
  Tire_Damping_b2: 15020.0         # N.s/m
  Disturbance_Force: 100.0         # N (t=1s -> 2s)

System_3_Inverted_Pendulum:
  Cart_Mass_M: 0.6    # kg
  Pendulum_Mass_m: 0.1 # kg
  Length_l: 0.3       # m
  Inertia_I: 0.006    # kg.m^2
  Friction_b: 0.1     # N.s/m

Phương pháp luận nghiên cứu (Methodology)

Quy trình nghiên cứu áp dụng mô hình phát triển hướng mô hình (Model-Based Design - MBD) 5 giai đoạn:

  1. Giai đoạn 1: Phân tích cơ sở lý thuyết vật lý: Thiết lập các phương trình cân bằng lực động học D'Alembert và định luật vòng điện áp Kirchhoff.
  2. Giai đoạn 2: Xây dựng cấu trúc Bond Graph: Xác định các nút đẳng thế (0-junction) và các nút đẳng dòng/vận tốc (1-junction), ghép nối các phần tử $R, C, I, GY, TF$.
  3. Giai đoạn 3: Thuật toán gán nhân quả (Sequential Causality Assignment Procedure - SCAP): Xác định quan hệ nhân quả tích phân cho toàn bộ các phần tử $I$ và $C$.
  4. Giai đoạn 4: Thiết kế mạch vòng điều khiển hồi tiếp: Xây dựng thuật toán điều khiển PID trong miền liên tục $u(t) = K_p e(t) + K_i \int e(t)dt + K_d \frac{de(t)}{dt}$.
  5. Giai đoạn 5: Tinh chỉnh tham số và đánh giá đáp ứng: Khảo sát độ nhạy và đánh giá chất lượng điều khiển theo các chỉ tiêu quá độ ($t_r, t_s, POT, e_{ss}$).

Implementation và kết quả

Quá trình phát triển và chi tiết giải thuật

1. Động cơ điện một chiều nam châm vĩnh cửu (PMDC Motor)

Áp dụng định luật II Newton cho phần cơ và định luật Kirchhoff cho mạch phần ứng:

$$\begin{cases} J \cdot \ddot{\theta}(t) + b \cdot \dot{\theta}(t) = K \cdot i(t) \ L \cdot \frac{di(t)}{dt} + R \cdot i(t) + K \cdot \dot{\theta}(t) = V(t) \end{cases}$$

Chuyển đổi sang miền Laplace và rút ra hàm truyền vòng hở từ điện áp đặt $V(s)$ tới vận tốc góc $\dot{\theta}(s)$:

$$P(s) = \frac{\dot{\theta}(s)}{V(s)} = \frac{K}{(J s + b)(L s + R) + K^2}$$

Thay các thông số thực tế ($L=0.002\text{ H}, R=0.8\ \Omega, b=0.0076\text{ Nms/rad}, J=0.1\text{ Nms}^2/\text{rad}, K=0.3$):

$$P(s) = \frac{0.3}{(0.1s + 0.0076)(0.002s + 0.8) + 0.09} = \frac{0.3}{0.0002 s^2 + 0.0800152 s + 0.09608}$$

// Mã nguồn cấu trúc SIDOPS mô phỏng PMDC Motor trên 20-sim
parameters
    real L = 2.0e-3;
    real R = 0.8;
    real J = 0.1;
    real b = 7.6e-3;
    real K = 0.3;
variables
    real u, i, omega, theta, e_back, torque;
equations
    // Mạch điện phần ứng
    L * ddt(i) = u - R * i - e_back;
    e_back = K * omega;
    
    // Ghép nối biến đổi cơ-điện qua Gyrator
    torque = K * i;
    
    // Chuyển động quay rotor
    J * ddt(omega) = torque - b * omega;
    ddt(theta) = omega;

2. Hệ thống treo xe bus 2 bậc tự do (Quarter-Car Model)

Phương trình vi phân mô tả chuyển động của thân xe ($x_1$) và bánh xe ($x_2$) khi chịu mấp mô mặt đường ($w$) và lực điều khiển chủ động ($u$):

$$\begin{cases} m_1 \ddot{x}_1 + b_1(\dot{x}_1 - \dot{x}_2) + k_1(x_1 - x_2) = u \ m_2 \ddot{x}_2 - b_1(\dot{x}_1 - \dot{x}_2) - k_1(x_1 - x_2) + b_2(\dot{x}_2 - \dot{w}) + k_2(x_2 - w) = -u \end{cases}$$

Ma trận hàm truyền hệ thống biểu diễn độ êm dịu thân xe $(x_1 - x_2)$ theo lực điều khiển $u(s)$ khi $w(s)=0$:

$$G_1(s) = \frac{X_1(s) - X_2(s)}{U(s)} = \frac{(m_1 + m_2)s^2 + b_2 s + k_2}{\Delta(s)}$$

Với định thức ma trận đặc trưng:

$$\Delta(s) = (m_1 s^2 + b_1 s + k_1)(m_2 s^2 + (b_1 + b_2)s + (k_1 + k_2)) - (b_1 s + k_1)^2$$

3. Hệ thống xe và con lắc ngược (Inverted Pendulum)

Áp dụng phương pháp hàm Lagrange $L = T - V$ với vector tọa độ suy rộng $q = [x, \theta]^T$:

$$T = \frac{1}{2}(M + m)\dot{x}^2 - m l \dot{x} \dot{\theta} \cos\theta + \frac{1}{2}(I + m l^2)\dot{\theta}^2$$ $$V = m g l \cos\theta$$

Phương trình Euler-Lagrange phi tuyến:

$$\begin{cases} (M + m)\ddot{x} + m l \dot{\theta}^2 \sin\theta - m l \ddot{\theta}\cos\theta + b\dot{x} = F \ (I + m l^2)\ddot{\theta} - m l \ddot{x}\cos\theta - m g l \sin\theta = 0 \end{cases}$$

Tuyến tính hóa quanh vị trí thẳng đứng cân bằng không bền ($\theta \approx 0 \Rightarrow \cos\theta \approx 1, \sin\theta \approx \theta, \dot{\theta}^2 \approx 0$):

$$G_{\theta}(s) = \frac{\Theta(s)}{U(s)} = \frac{\frac{m l}{q} s}{s^3 + \frac{b(I + m l^2)}{q}s^2 - \frac{(M + m)m g l}{q}s - \frac{b m g l}{q}}$$

Trong đó hệ số quán tính tương đương: $q = (M + m)(I + m l^2) - (m l)^2$.

Thử nghiệm và đánh giá chất lượng điều khiển

Các bài toán mô phỏng được thực hiện trên phần mềm 20-sim để so sánh định lượng chất lượng điều khiển giữa các cấu hình bộ điều khiển.

       ĐÁP ỨNG QUÁ ĐỘ ĐIỀU KHIỂN VỊ TRÍ ĐỘNG CƠ DC (Step = 1 rad)
  1.4 |               ... PI (Vọt lố cao)
      |              .   .
  1.2 |             .     .
  1.0 |------------+-------+--------------------- Target (1.0 rad)
      |           / \     / \
  0.8 |          /   .....   \____ PID Tối ưu (Kp=1, Ki=100, Kd=0.85)
      |         /                  (t_s < 1s, POT < 5%)
  0.4 |        /
      |       /   --- P Thuần túy (Đáp ứng chậm, t_s ~ 7s)
  0.0 +------+------+------+------+------+------+---> Time (s)
      0      1      2      3      4      5      6

Bảng kết quả thực nghiệm mô phỏng trên 20-sim

Hệ thống khảo sát Cấu hình bộ điều khiển Tham số hiệu chỉnh ($K_p, K_i, K_d$) Độ vọt lố ($POT$) Thời gian xác lập ($t_s$) Sai số xác lập ($e_{ss}$)
Động cơ DC Vòng hở (Open-loop) Không điều khiển ($V=24\text{V}$) Không xác lập Tăng tuyến tính $\infty$
Bộ điều khiển P $K_p = 1$ $0%$ $7.0\text{ s}$ Nhỏ
Bộ điều khiển PI $K_p = 1, K_i = 100$ $18.5%$ $3.2\text{ s}$ $0%$
Bộ điều khiển PD $K_p = 1, K_d = 0.85$ $0%$ $1.0\text{ s}$ $1.2%$
Bộ điều khiển PID $\mathbf{K_p = 1, K_i = 100, K_d = 0.85}$ $\mathbf{< 4.2%}$ $\mathbf{0.95\text{ s}}$ $\mathbf{0.0%}$
Hệ treo xe bus Không điều khiển $u = 0$ (Xung lực $100\text{N}$) Biên độ đỉnh $55.0$ $8.0\text{ s}$ Dao động tắt dần
PID thử nghiệm 1 $K_p = 10000, K_i = 5, K_d = 100$ Dao động mạnh $> 6.0\text{ s}$ Mất ổn định đầu
PID thử nghiệm 2 $K_p = 100000, K_i = 100, K_d = 5$ Biên độ đỉnh $1.2$ $4.0\text{ s}$ Ổn định dần
PID Tối ưu hóa $\mathbf{K_p = 100000, K_i = 500, K_d = 0.2}$ Biên độ đỉnh 2.8 $\mathbf{< 1.8\text{ s}}$ Triệt tiêu nhanh
Con lắc ngược Vòng hở $u = 0$ Mất ổn định Rơi tức thì $\infty$
PID chưa chỉnh $K_p = 1, T_i = 1, T_d = 1$ Dao động lớn Không ổn định Không hội tụ
PID Tinh chỉnh $\mathbf{K_p = 1000, T_i = 1, T_d = 0.05}$ $\mathbf{< 3.8%}$ $\mathbf{1.2\text{ s}}$ $\mathbf{0.0^\circ}$

Kết quả đạt được

  1. Kiểm soát vị trí động cơ PMDC: Cấu hình PID ($K_p=1, K_i=100, K_d=0.85$) giúp góc quay đạt chính xác vị trí đặt $1\text{ rad}$ chỉ sau $0.95\text{s}$, độ vọt lố được khống chế dưới $4.2%$, triệt tiêu hoàn toàn sai số xác lập.
  2. Dập tắt dao động hệ treo xe bus: Với bộ điều khiển PID ($K_p=100000, K_i=500, K_d=0.2$), biên độ dịch chuyển thân xe khi gặp xung kích $100\text{N}$ giảm từ $55.0$ xuống còn $2.8$ (giảm chấn hiệu quả đạt $94.9%$), thời gian dập tắt dao động rút ngắn từ $8\text{s}$ xuống dưới $1.8\text{s}$ (giảm $77.5%$).
  3. Cân bằng con lắc ngược: Giữ vững con lắc ở phương thẳng đứng $(\theta = 0\text{ rad})$ khi có tác động lực đẩy, thời gian quá độ đạt $1.2\text{s}$.

Đổi mới và đóng góp

  • Thống nhất phương pháp luận Bond Graph cho hệ đa vật lý: Thay vì giải quyết rời rạc từng hệ phương trình vi phân phức tạp, đồ án đã chuẩn hóa việc biểu diễn động cơ điện (Điện - Từ - Cơ quay), hệ treo (Cơ tịnh tiến - Giảm chấn thủy lực) và con lắc ngược (Cơ học phi tuyến) trên cùng một ngôn ngữ đồ hình liên kết.
  • Quy trình gán nhân quả tự động hóa: Áp dụng triệt để quy tắc SCAP giúp loại bỏ hoàn toàn các biến trạng thái phụ thuộc (dẫn đến đạo hàm số gây mất ổn định giải thuật), đảm bảo $100%$ các phần tử tích năng ở dạng tích phân thuần nhất ($I$-integral, $C$-integral).
  • Chiến lược tối ưu hóa tham số PID thực nghiệm: Xác lập được quy luật điều chỉnh tham số cho hệ thống treo xe bus tải trọng lớn: $K_p \gg K_i > K_d$ ($K_p = 10^5, K_i = 500, K_d = 0.2$), mang lại hiệu quả giảm chấn vượt trội so với các phương pháp chỉnh định Ziegler-Nichols cổ điển.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong thực tế

  • Hệ thống lái và trợ lực điện tử (EPS / Steer-by-Wire): Thuật toán điều khiển vị trí động cơ PMDC được ứng dụng trực tiếp trong các cơ cấu chấp hành tay lái xe điện thông minh, đảm bảo góc quay chính xác và phản hồi mô-men tức thì.
  • Hệ thống treo khí nén chủ động xe bus đô thị (Active Air Suspension): Mô hình hệ treo 2-DOF giúp tối ưu van điều tiết thủy lực/khí nén, nâng cao độ êm dịu chuyển động (Ride Comfort) và ổn định thân xe khi vào cua hoặc phanh gấp.
  • Robot tự cân bằng hai bánh (Self-Balancing Transporter): Cơ sở toán học của hệ con lắc ngược là nền tảng trực tiếp phát triển các dòng xe Segway, robot hai bánh vận chuyển hàng hóa trong nhà kho thông minh (AGV).
                            LỘ TRÌNH TRIỂN KHAI PHẦN CỨNG (ROADMAP)
  
  [ Mô hình hóa 20-sim ] ---> [ Tự động xuất mã C ] ---> [ Triển khai Vi điều khiển ]
        (Bond Graph)                 (ANSI C Code)            (STM32F4 / TMS320F28335)
             |                            |                              |
             v                            v                              v
     Kiểm chứng mô phỏng           Hardware-in-the-Loop           Thử nghiệm cơ cấu thực tế
      (Milestones: M1)              (Milestones: M2)                 (Milestones: M3)

Yêu cầu triển khai phần cứng

  • Bộ vi điều khiển khuyến nghị: Vi điều khiển 32-bit DSP hoặc ARM Cortex-M4/M7 (ví dụ: STM32F407VGT6 hoặc TI TMS320F28335) hỗ trợ đơn vị tính toán số thực FPU phần cứng.
  • Cảm biến và cơ cấu chấp hành:
    • Động cơ DC: Cảm biến Encoder quang học độ phân giải tối thiểu 1024 xung/vòng (PPR).
    • Hệ treo: Cảm biến gia tốc MEMS 3 trục (IMU) và cảm biến khoảng cách LVDT đo hành trình nhún.
    • Con lắc ngược: IMU 6-DOF (MPU6050/ICM-20948) kết hợp bộ lọc Kalman ước lượng góc nghiêng $\theta$.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Mô hình tuyến tính hóa cục bộ: Hệ con lắc ngược chỉ cân bằng tốt ở dải góc nhỏ ($|\theta| < 15^\circ$); khi góc nghiêng vượt quá giới hạn, hệ thống mất kiểm soát do ảnh hưởng của các thành phần phi tuyến $\sin\theta, \cos\theta, \dot{\theta}^2$.
  2. Hiện tượng trễ và bão hòa tín hiệu (Actuator Saturation): Mô phỏng chưa tích hợp giới hạn điện áp nguồn cấp thực tế ($V_{max} = 24\text{V}$) và dòng điện cực đại của mạch cầu H (H-Bridge Driver).
  3. Nhiễu đo lường thực tế: Chưa khảo sát ảnh hưởng của nhiễu trắng cảm biến (Sensor White Noise) và trễ truyền thông CAN/PWM.

Hướng phát triển tiếp theo

  • Tích hợp bộ điều khiển hiện đại: Phát triển thuật toán điều khiển tối ưu tuyến tính toàn phương LQR (Linear Quadratic Regulator) và điều khiển trượt (Sliding Mode Control - SMC) cho bài toán con lắc ngược với góc khởi tạo lớn.
  • Điều khiển thích nghi (Adaptive / Fuzzy-PID): Tự động thích ứng tham số $K_p, K_i, K_d$ của hệ thống treo khi tải trọng thân xe bus thay đổi từ không tải ($2500\text{kg}$) sang đầy tải ($4500\text{kg}$).
  • Kiểm nghiệm trên mô hình thực nghiệm (Hardware-in-the-Loop - HIL): Kết nối 20-sim với hệ thống dSPACE hoặc NI CompactRIO để kiểm thử thời gian thực với phần cứng điều khiển.

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

+-------------------------------------------------------------------------------+
|                           ĐỐI TƯỢNG HƯỞNG LỢI                                 |
+-------------------------------------------------------------------------------+
| * Sinh viên ngành Cơ điện tử / Tự động hóa:                                   |
|   Nắm vững công cụ Bond Graph, làm chủ phần mềm 20-sim, nâng cao kỹ năng      |
|   mô hình hóa hệ thống vật lý đa miền.                                        |
+-------------------------------------------------------------------------------+
| * Kỹ sư R&D và Lập trình nhúng:                                              |
|   Sở hữu bộ thông số chuẩn định và thuật toán PID đã được kiểm chứng số học,  |
|   rút ngắn thời gian phát triển firmware điều khiển động cơ và hệ thống treo. |
+-------------------------------------------------------------------------------+
| * Doanh nghiệp sản xuất ô tô và thiết bị công nghiệp:                         |
|   Tối ưu hóa chi phí thử nghiệm thực địa thông qua quy trình Virtual Prototyping,|
|   nâng cao độ tin cậy và tiện nghi của sản phẩm thương mại.                   |
+-------------------------------------------------------------------------------+
| * Giảng viên và Nhà nghiên cứu:                                               |
|   Tài liệu tham khảo học thuật trực quan cho các môn học Mô hình hóa & Mô phỏng, |
|   Lý thuyết điều khiển tự động và Động lực học hệ nhiều vật.                  |
+-------------------------------------------------------------------------------+

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

1. Yêu cầu kỹ thuật phần mềm và phần cứng để chạy mô phỏng đồ án là gì?

Hệ thống yêu cầu cài đặt phần mềm Controllab 20-sim (phiên bản 4.8 hoặc 5.0 trở lên) trên hệ điều hành Windows 10/11 hoặc Linux (qua Wine). Cấu hình phần cứng tối thiểu: CPU Intel Core i3 thế hệ 8 trở lên, RAM 4GB, đồ họa tích hợp hỗ trợ OpenGL để hiển thị mô hình 3D Animation.

2. Ưu điểm nổi bật của Bond Graph so với sơ đồ khối Simulink trong mô hình hóa động cơ DC là gì?

Bond Graph mô tả chính xác tính hai chiều của năng lượng qua phần tử Gyrator ($GY$): Dòng điện sinh ra mô-men quay ($T = K \cdot i$) và vận tốc góc đồng thời sinh ra sức phản điện động ($e_{back} = K \cdot \omega$). Simulink đòi hỏi lập hai khối tách rời và dễ gây lỗi vòng lặp đại số (Algebraic Loop) nếu không xử lý đúng biến trạng thái.

3. Tại sao trong điều khiển hệ thống treo xe bus, hệ số $K_p$ phải chọn rất lớn ($100000$)?

Khối lượng thân xe bus rất lớn ($m_1 = 2500\text{kg}$) kết hợp độ cứng lò xo cao ($k_1 = 60000\text{N/m}$). Để cơ cấu chấp hành chủ động sinh đủ lực $u(t)$ đối kháng lại xung lực mặt đường và giữ ổn định thân xe trong dải hành trình ngắn, hệ số khuếch đại tỉ lệ $K_p$ phải đủ lớn để tạo phản lực tức thời.

4. Có thể chuyển trực tiếp mô hình từ 20-sim sang vi điều khiển như thế nào?

20-sim cung cấp công cụ C-code Generation. Sau khi mô hình hóa và kiểm chứng thuật toán điều khiển PID thành công, phần mềm có thể tự động biên dịch sơ đồ khối sang mã nguồn ANSI-C tiêu chuẩn. Lập trình viên có thể tích hợp trực tiếp file .c và .h này vào dự án Keil C, STM32CubeIDE hoặc Code Composer Studio.

5. Chi phí triển khai thử nghiệm hệ thống treo chủ động thực tế khoảng bao nhiêu?

Chi phí cho một cụm thử nghiệm Quarter-Car thực tế quy mô phòng thí nghiệm dao động từ 30 - 50 triệu VNĐ (bao gồm động cơ servo chấp hành hoặc van thủy lực tỉ lệ, loadcell, cảm biến gia tốc, vi điều khiển công nghiệp và bộ nguồn công suất). Việc mô phỏng ảo trên 20-sim giúp loại trừ rủi ro hỏng hóc phần cứng, tiết kiệm hơn $80%$ chi phí R&D ban đầu.

Kết luận

Đồ án môn học "Mô Hình Hóa và Mô Phỏng Hệ Thống Cơ Điện Tử" của nhóm sinh viên Trường Đại học Công nghiệp Hà Nội đã thực hiện thành công việc kết nối chặt chẽ giữa lý thuyết giải tích động lực học cổ điển với phương pháp luận đồ hình liên kết Bond Graph hiện đại trên phần mềm 20-sim. Các kết quả nghiên cứu đã giải quyết trọn vẹn 3 bài toán kỹ thuật trọng điểm: điều khiển vị trí động cơ PMDC đạt độ chính xác cao ($POT < 5%$), giảm chấn hệ thống treo xe bus đạt hiệu quả dập tắt dao động $94.9%$, và giữ cân bằng bền vững cho hệ con lắc ngược phi tuyến. Đây là nền tảng học thuật và thực tiễn vững chắc phục vụ công tác thiết kế, phát triển các hệ thống điều khiển thông minh trong công nghiệp ô tô, robot tự hành và cơ điện tử chính xác.