Lập trình vẽ biên dạng cánh NACA 4 kí tự bằng MATLAB

Trần Quang Khôi | MSSV : 1411875 | Khí động lực học I/ Cấu trúc chương trình VẼ BIÊN DẠNG CÁNH NACA KÍ TỰ BẰNG PHẦN MỀM MATLAB Page Page II/ Ý nghĩa thông số họ biên dạng cánh NACA kí tự - Họ biên dạng cánh NACA kí tự xác định hoàn toàn phương trình đường cong phương trình phân bố độ dày Trong tên biên dạng cánh ( kí hiệu kí tự ) có ý nghĩa khác giúp xác định số cho phương trình đề cập, cụ thể sau : + Số xác định độ cong tối đa đường cong biên dạng cánh tính theo phần trăm dây cung cánh Ví dụ : NACA 2112 cho biết độ cong tối đa biên dạng 2% độ dài dây cung cánh, tức biên dạng cánh có độ dài 1m độ cong tối đa 2cm + Số thứ xác định vị trí theo trục x có độ cong tối đa tính theo phần mười dây cung cánh Ví dụ : NACA 2112 cho biết vị xmc = 1/10c - trí có độ cong tối đa 2% nằm vị trí có tọa độ , hay cụ thể 10cm tính từ đầu mép cánh cánh có độ dài dây cung 1m +2 số cuối cho biết độ dày tối đa cánh tính theo phần trăm dây cung cánh Ví dụ : NACA 2112 cho biết biên dạng cánh có độ dày tối đa 12% dây cung cánh, cụ thể 12cm cánh có độ dài dây cung 1m Như nhu cầu đặt chương trình phải trích xuất giá trị từ input người dùng tên biên dạng cánh với kí tự Chương trình giải vấn đề cách cho phép người dùng nhập vào chuỗi kí số (vd: 2112) sau trích xuất số vào biến khởi tạo Minit, số thứ vào biến khởi tạo Pinit số cuối vào biến khởi tạo Tinit thông qua hàm str2double matlab, có tác dụng chuyển dạng liệu từ chuỗi sang giá trị số thực Sau tính toán giá trị số theo tỷ lệ ý nghĩa kí số đề cập bên III/ Các phương trình sử dụng III.1/ Phương trình đường cong biên dạng cánh (camber line) Page - yc Phương trình đường cong biên dạng cánh giúp xác định tọa độ tập hợp điểm tạo thành đường cong (camber line) theo x Từ vẽ đường cong thông qua ma trận tọa độ [ yc ] Trong : m : độ cong tối đa biên dạng cánh p : vị trí có độ cong tối đa ( xmc ) m = M init /100 p = Pinit /10 Code chương trình cụ thể sau : yc = ones(gridPts,1); if (x(i) >= && x(i) < P) yc(i) = (M/P^2)*((2*P*x(i))-x(i)^2); elseif (x(i) >=P && x(i)