Arora S IOD Chapter 2
Arora S IOD Chapter 2
TO OPTIMUM
DESIGN
THIRD EDITION
JASBIR S. ARORA
The University of Iowa
College of Engineering
Iowa City, Iowa
ii
Table of Contents
Chapters
1 Introduction to Design Optimization 1-1
2 Optimum Design Problem Formulation 2-1
3
Graphical Optimization and Basic Concepts 3-1
4
Optimum Design Concepts: Optimality Conditions 4-1
5
More on Optimum Design Concepts: Optimality Conditions 5-1
6
Optimum Design with Excel Solvers 6-1
7
Introduction to Optimum Design with MATLAB 7-1
8
Linear Programming Methods for Optimum Design 8-1
9
More on Linear Programming Methods for Optimum Design 9-1
10 Numerical Methods for Unconstrained Optimum Design 10-1
11 More on Numerical Methods for Unconstrained Optimum Design 11-1
12 Numerical Methods for Constrained Optimum Design 12-1
13 More on Numerical Methods for Constrained Optimum Design 13-1
14 Practical Applications of Optimization 14-1
17 Multi-objective Optimum Design Concepts and Methods 17-1
Appendix
A
iii
1-1
1-2
Chapter 1 Introduction
1-3
CHAPTER
2
Optimum Design Problem Formulation
2.1___________________________________________________________________________
A 100 100 m lot is available to construct a multistory office building. At least 20,000 m2 total floor
space is needed. According to a zoning ordinance, the maximum height of the building can be only
21 m, and the area for parking outside the building must be at least 25 percent of the floor area. It
has been decided to fix the height of each story at 3.5 m. The cost of the building in millions of
dollars is estimated at 0.6h +0.001A, where A is the cross-sectional area of the building per floor and
h is the height of the building. Formulate the minimum cost design problem.
Solution
Step 1: Problem Statement
Shown above
Step 2: Data and Information Collection
Area of the lot =100100, m2
Area available for parking = (10000 A), m2
Total floor area = (number of floors)A =
, m2
Step 3: Definition of Design Variables
A = cross-sectional area of the building for each floor, m 2
h = height of the building, m
Step 4: Optimization Criterion
Optimization criterion is to minimize $ cost, and the cost function is defined as
Cost = (0.6h + 0.001A), million dollars
Step 5: Formulation of Constraints
Floor Space Constraint: hA/3.5 20,000
Parking Constraint: (10,000 A ) 0.25hA/3.5
Explicit Design Variable Constraints:
h 3.5, m
h 21, m
A 0
2-1
Market (bbl)
20,000
10,000
Solution
Step 1: Problem Statement
Shown above
Step 2: Data and Information Collection
Shown above
Step 3: Definition of Design Variables
A = Crude A used in barrels
B = Crude B used in barrels
Step 4: Optimization Criterion
Optimization criterion is to maximize profit, and the cost function is defined as
Profit = 200(0.6A + 0.8B) + 450(0.4A + 0.2B) 120A 150B = 180A + 100B
Step 5: Formulation of Constraints
Gasoline Market Constraint:
(0.6A + 0.8B) 20,000, bbl
Lube Oil Market Constraint:
(0.4A + 0.2B) 10,000, bbl
Explicit Design Variable Constraints:
A 20,000, bbl
B 30,000, bbl
A 0; B 0
2-2
2-3
2-4
2-5
2-6
2-7
2-8
2-9
the container
the contents
2 18 80 2 b + 4 bh
2
18 150 b 2 h
2
32 b + 2 bh 32 1 2
=
=
+
2
15 b h 15 h b
2-10
2-11
2-12
2-13
2-14
i 1
j 1
ij
x ij
ij
bi
for i = 1 to m
j 1
Demand Constraint:
ij
i 1
2-15
1
2
2H
2-16
1
2
1
2
2H ( Dt ) a
or
W s
4 H
1
2
2H
2
3
3
E D t Dt 8
FS s 2
4 H
2
3
3
E D t + Dt 8
4+H
1
2
1
2
2 H D t 250 10
2 H 210 10
t Dt
H 5, m;
H 2, m;
D/t 90;
0.1 D 2, m;
0.01 t 0.1, m
2-17
16 1 H
where d is the depth and b is the width of the beam. It is also desired that the depth of the beam shall
not exceed twice its width. Formulate the design problem for minimum cross-sectional area using
the following data: M=140 kN m, V=24 kN, =165 MPa, =50 MPa.
2-18
2-19
2-20
2-21
Truck load
Capacity
(tones)
Average truck
speed
(km/h)
Crew required
per shift
A
B
C
10
20
18
66
50
50
1
2
2
No. of hours
of operations
per day
(3 shifts)
19
18
21
Cost of each
truck($)
40,000
60,000
70,000
Solution
Step 1: Problem Statement
Shown above
Step 2: Data and Information Collection
Shown above
Step 3: Definition of Design Variables
A, B and C are the number of trucks to be purchased of the type A, B and C, respectively.
Step 4: Optimization Criterion
Optimization criterion is to maximize the capacity (tonnes kilometers), and the cost function is
defined as
Capacity = A(10 55 18) + B(20 50 18) + C(18 50 21) = 9900A + 18000B + 18900C
Transcribing into the standard form, we get
minimize f = 9900A 18000B 18900C
Step 5: Formulation of Constraints
Available Capital Constraint: A(40,000) + B(60,000) + C(70,000) 2,000,000
Limit on Available Drivers Constraint: 3A + 6B + 6C 150
Limitation on Maintenance Facility Constraint: A + B + C 30
Explicit Design Variable Constraints: A, B, C 0
2-22
2-23
k(1) = 190,000
k(2) = 240,000
k(3) = 290,000
D(1) = 330,000
D(2) = 125,000
b(1,1,2)=0.84
b(2,1,2)=0.81
b(1,2,2)=0.85
b(2,2,2)=0.77
b(1,3,2)=0.72
b(2,3,2)=0.78
Solution
Several formulations for the design problem are possible. For each formulation proper design
variables are identified. Expressions for the cost and constraint functions are derived.
2-24
g1: x i 1,200,000;
i 1
12
g2: y j 1,000,000
j 1
(2) The total tonnage of ingot stocks manufactured into products at each fabricating plant cannot
exceed its available capacity, i.e., F1 k(1); F2 k(2); F3 k(3) where F1, F2 and F3
represent the total tonnage of ingot stocks processed at three fabricating plants separately. In
2-25
( e i x i f i y i ) 330, 000;
h2 :
i 1,3,5,7,9,11
( e i x i f i y i ) 125, 000
i 2,4,6,8,10,12
where e i 's and f i 's are coefficients transferring tonnage of iron ore into products. These
coefficients are given as:
e1
e3
e5
e7
e9
e11
f1
f3
f5
f7
f9
f11
e2
e4
e6
e8
e10
e12
f2
f4
f6
f8
f10
f12
(4) There are constraints requiring that both the reduction plants and fabricating plants do not
have any inventory of their own. These constraints have been satisfied automatically since the
twenty-four design variables (paths) are chosen which satisfy these conditions.
Summarizing and transcribing into the standard model, we get
h1 = 0.3081 x1 + 0.3128 x3 + 0.2847 x5 + 0.3082 x7 + 0.2886 x9 + 0.2852 x11 + 0.3476 y1 + 0.3264 y3
+ 0.3212 y5 + 0.3216 y7 + 0.3256 y9 + 0.2976 y11 330,000 = 0
h2 = 0.3276 x2 + 0.3726 x4 + 0.3315 x6 + 0.3542 x8 + 0.2808 x10 + 0.3588 x12 + 0.3696 y2 + 0.3888 y4
+ 0.3740 y6 + 0.3696 y8 + 0.3168 y10 + 0.3744 y12 125,000 = 0
g1 = x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 + x9 + x10 + x11 + x12 1,200,000 0
g2 = y1 + y2 + y3 + y4 + y5 + y6 + y7 + y8 + y9 + y10 + y11 + y12 1,000,000 0
g3 = 0.39 x1 + 0.39 x2 + 0.46 x3 + 0.46 x4 + 0.44 y1 + 0.44 y2 + 0.48 y3 + 0.48 y4 190,000 0
g4 = 0.39 x5 + 0.39 x6 + 0.46 x7 + 0.46 x8 + 0.44 y5 + 0.44 y6 + 0.48 y7 + 0.48 y8 240,000 0
g5 = 0.39 x9 + 0.39 x10 + 0.46 x11 + 0.46 x12 + 0.44 y9 + 0.44 y10 + 0.48 y11 + 0.48 y12 290,000 0
x i 0, y i 0, i = 1 to 12
2-26
2-27
2-28
2-29
H1
w1 2
H2 1
w2
2 2
H1 H 2 1
w3
w1 2
H1
w w1
H2 1
w2 2
= 0; h4
H1
w w1
H1 H 2 1
w 3 x1
=0
2-30
Transcribe the problem into the standard design optimization model (also use
40.0 cm,
40.0 cm). Use the following data: P = 14 kN ; l = 10m; mass density, =7850kg/m3, allowable
bending stress, = 165MPa, Allowable shear stress, =50MPa.
2-31
L R o R i
2
7.85 10 1000 R
3
2
o
Ri
24.66 R
2
o
R i , kg
2
; or 4 1.4 10 4 10 3 R o R o4 R i4 1.65 10 4 0; or
R o4 R i4 1.65 10 4 0
g1 1.7825 10 7 R o
g2 : 4 P R o2 + R o R i + R i2 3 R o4 R i4
4 1.4 10
; or
R + R R + R 3 R R 5000 0 ; or
g2 5941.78 R + R R + R R R 5000 0
4
2
o
2
o
2
i
4
o
2
i
4
o
4
i
4
i
g3 R o 40 0; g 4 R i 40 0 ; g 5 R o 0 ; g 6 R i 0
2-32
2-33
2 Rt
), i.e.,
The deflection
52
g RtL
384 E R t
3
0.001 L 0,
or
gL
384 E R
10
g 3 R t 20 0 ; g 4 R 20 0 ; g 5 R 200 0
; g 6 t 1 0 ; g 7 t 10 0
or
g1 : P 2 Rt
g2 :
10
g L
384 ER
0.001 L , or
10 7.8 10
9.8 3000
384 2.1 10 R
6
0.001 3000
g 2 767.8 R
3 0
g 3 R t 20 0
g 4 20 R 0
g 5 R 200 0
g 6 1 t 0;
g 7 t 10 0
2-34