Chapter 4 Optimum Design Concepts: Optimality Conditions
There are 16 cases because there are 4 inequality constraints. The case which yields a solution is given as
. The solution is
12
0 4267, 0 00658,uu= =..
. It is
seen from the graph for the problem that the solution is indeed a minimum point and the constraints on
maximum b and h are active.
MATLAB Code for Exercise 4.92
clear all
[b,h]=meshgrid(–3:0.1:15 , –3:0.1:23);
f=(32/15)*(1./h+2./b);
g1=b–10;
cv1=[0:0.03:0.5];
const1=contour(b,h,g1,cv1,‘g’);
cv1=[0 0.001];
const1=contour(b,h,g1,cv1,‘k’);
cv4=[0:0.05:0.9];
const4=contour(b,h,g4,cv4,‘g’);
cv4=[0 0.01];
const4=contour(b,h,g4,cv4,‘k’);
fv=[0.5 0.545185 1 2];
fs=contour(b,h,f,fv,‘b’);