Ozymandias
Well-Known Member
- Joined
- Oct 21, 2014
- Messages
- 544
- Reaction score
- 0
I got to do a pretty cool project for my Computing for Engineers class. I wrote two simple scripts that will calculate the shape for HAACK series nose cones. Below is the script for the popular LD-HAACK (aka. Von Karmann Ogive). I have it set up for a 15" long, 3" diameter nose. If you would like a nose of a different size, simply change L (nose length) and R (nose radius @ base) accordingly. This script will work in MATLAB as well as GNU Octave. You could also plot it on a graphing calculator with a little reverse engineering. I hope this helps someone out there. Script follows:
%LD-HAACK Cone Calculator
%Written by David A. Roy
C = 0;
L = 15;
R = 1.5;
x = linspace(0, 15);
theta = acos(1-(2*x)/(L));
y = R*sqrt((theta-sin(2*theta)/2+C*(sin(theta)).^3))/(sqrt(pi));
plot(x,y)
hold on
plot(x,-y)
axis('equal');
xlabel('length')
ylabel('radius')
title('Nose Profile')
%LD-HAACK Cone Calculator
%Written by David A. Roy
C = 0;
L = 15;
R = 1.5;
x = linspace(0, 15);
theta = acos(1-(2*x)/(L));
y = R*sqrt((theta-sin(2*theta)/2+C*(sin(theta)).^3))/(sqrt(pi));
plot(x,y)
hold on
plot(x,-y)
axis('equal');
xlabel('length')
ylabel('radius')
title('Nose Profile')