Scribd
Upload
33 views2 pages

All Off: For For End End

The document discusses optimizing the dimensions of an I-beam to minimize cost. It defines parameters like load, material properties, and initial beam thickness and width. Contour plots are generated showing the beam's area for different depth and height values, along with failure limits from bending and shear. The objective is to vary depth and height to find the optimal solution within the failure constraints.

Uploaded by

kwzeet
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
33 views2 pages

All Off: For For End End

The document discusses optimizing the dimensions of an I-beam to minimize cost. It defines parameters like load, material properties, and initial beam thickness and width. Contour plots are generated showing the beam's area for different depth and height values, along with failure limits from bending and shear. The objective is to vary depth and height to find the optimal solution within the failure constraints.

Uploaded by

kwzeet
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

% Opt of I-Beam for min cost.

clear
clear all
hold off
clf
% Parameters
density = 0.098; % lbf/in^3
length = 85*12; % in
w_load = 585000; % lbf
yield_strength = 2.5*10^4; %lbf/in^2
shear_strength = 1.9*10^4; % lbf/in^2
E_youngs = 10^7; % lbf/in^2
% Objective Function:
% Assum t and b in order to create plot with two variables.
% Vary t and b to optimize solution further.
t = 1.4; % in
b = 7; % in
d = 0:1:250; % in, based on data from "Mechanics of Materials" Ferdinand
h = 0:1:250; % in, based on data from "Mechanics of Materials" Ferdinand
[x1,x2] = meshgrid(d,h);
f = b*x1 - b*x2 +t*x2;
% Contour map for h vs d
figure(1)
hold on
contour(x1,x2,f,60)
% Condition for which beam is an I-beam
plot(d,h,'c--','LineWidth',3) % Plot condition d-h > 0
for i = 1:numel(d)
for j = 1:numel(h)
w_beam(i,j) = w_load + (b*d(i)-b*h(j)+t*h(j))*length*density;
end
end

% Plot limit of failure in bending


for i = 1:numel(d)
yield_limit(i) = ((3*w_load*length*d(i)-yield_strength*b*d(i)^3)/(-yield_streng
end
plot(d,yield_limit,'r-','LineWidth',2)
% Plot limit for shear
for i = 1:numel(d)
shear_limit(i) = (b*d(i)/(b-t))-w_load/(shear_strength*(b-t));
end
plot(d,shear_limit,'b-','LineWidth',2);
% Figure labeling
title('h vs. d for bending of I-beam')
legend('Area Contour','I-Beam','Bending','Shear')
xlabel('d [in]')
ylabel('h [in]')
hold off
Warning: Imaginary parts of complex X and/or Y arguments ignored

Published with MATLAB 7.9

You might also like