Flow over bicycle

PART1:

%Program to claculate the drag force

clear all
close all
clc

% Inputs
% Drag coefficient
c_d=0.8

% Area m^2
A=1

% Density kg/m^3
rho=1.2

% velocity m/s
v= [1:50];

drag_force=0.5*v.^2*rho*c_d*A;

plot(v,drag_force)
xlabel(\'Velocity\')
ylabel(\'Drag Force\')

Part 2:

% drag force for various shapes and their drag coefficient

clear all
close all
clc

% Inputs
% Area m^2
A=1

% Density kg/m^3
rho=1.2

% velocity m/s
v= [1:50]

%streamline body
c1=0.04

%half sphere
c2=0.43

%sphere
c3=0.6

%cone
c4=0.5

%long cylinder
c5=0.82

%cube
c6=1.05

drag_force=0.5*v.^2*rho*c3*A;

plot(c3,drag_force)
xlabel(\'cd\')
ylabel(\'Drag Force\')
grid on