作者:Sreetam Bhaduri
%% Wankel Engine Housing Profile
clear;
clc;
fprintf(‘It’‘s an no-parametric program. So please provide all the inputs in same system of measurement\n’);
%% User Inputs
e=input(‘Please provide the eccentricity magnitude\n’);
fprintf(‘K Factor=Radius of Rotor/Eccentricity\nIt’‘s value ranges in [0,4].\n’);
k_factor=input(‘What’‘s the K Factor of the engine?\n’);
%% Solver
theta=[0:pi/1000:2*pi/3];
R=e*k_factor;
if k_factor<=4
for n=1:1:3
x=(e*cos(3*theta))+(R*cos(theta+(2*n*pi/3)));
y=(e*sin(3*theta))+(R*sin(theta+(2*n*pi/3)));
%% Plotting parameters
plot(x,y,'LineWidth',2);
title('Wankel Engine Housing Profile');
xlabel('X Axis');
ylabel('Y Axis');
hold on;
end
else
disp('Learn the basics for Wankel Engine, then try again.');
end
%% End of MATLAB Program
下载地址:
http://page2.dfpan.com/fs/elc5j2f21329216ec07/
更多精彩文章请关注微信号: