- Đạo hàm do vận tốc góc quanh trục Zb ,r
4.3 Chương trình tính toán
Từ lý thuyết đã xét và xây dựng sơ đồ khối giải thuật như trên chúng ta có thể xây dựng một chương trình được lập trình bằng ngôn ngữ Matlab để có thể giải quyết vấn đề nêu trên cho một máy bay mà các thông số ban đầu về hình học cũng như là khí động có thể thay đổi được, giúp chúng ta có thể so sánh đưa ra những nhận xét kết luận về sự ổn định của một máy bay hạng nhẹ.
Dưới đây là giao diện của chương trình “Tính toán ổn định của máy bay”, bằng một ngôn ngữ khá mạnh là Matlab chương trình có thể dễ dàng đưa ra các hình vẽ, đồ thị mà chúng ta có thể xem xét một cách trực quan:
- Giao diện mà ta bắt gặp lúc đầu là phần nhập các thông số hình học, khí động ... của máy bay mà ta quan tâm, ở đây các thông số ban đầu là mặc định cho máy bay thông dụng loại nhỏ là máy bay NAVION.
- Ở phía dưới chúng ta có thể xem xét tính ổn định của máy bay ở chế độ dọc hoặc ngang bằng cách nhấn chuột vào nút “Kết quả ổn định dọc hoặc ngang”.
- Trong chế độ chuyển động dọc thì ta lại có 2 loại là ổn định dọc với chu kỳ dài, tần số thấp và ổn định dọc chu kỳ ngắn, tần số cao.
- Trong chế độ chuyển động ngang thì ta có 3 loại là Roll mode, Dutch Roll mode và Spiral mode. Chúng ta có thể xem các đáp ứng về vận tốc, vận tốc góc và các góc Euler bằng kết quả đạt được trên đồ thị.
Hình * Đáp ứng theo thời gian của các đại lượng chuyển động theo trục dọc với chu kỳ dài
Hình * Đáp ứng theo thời gian của các đại lượng chuyển động theo trục dọc với chu kỳ ngắn
Từ đồ thị này chúng ta có thể rút ra nhận xét là máy bay này ổn định ở chế độ đang xét vì các đại lượng vật lý trên theo diễn biến thời gian có biên độ giảm dần tới giá trị cân bằng ban đầu trong cả chuyển động với chu kỳ ngắn lẫn chu kỳ dài.
Dưới đây là kết quả tính toán đối với một máy bay lớn hơn là máy bay CONVAIR 880
Đồ thị biểu diễn sự biến thiên các đại lượng ổn định dọc trong chế độ chu kỳ dài
Đồ thị biểu diễn sự biến thiên các đại lượng ổn định ngang trong chế độ Roll mode
Đồ thị biểu diễn sự biến thiên các đại lượng ổn định ngang trong chế độ Dutch Roll mode
Đồ thị biểu diễn sự biến thiên các đại lượng ổn định ngang trong chế độ Spiral mode