Preliminary Vehicle Structure Design Application: Michael J. Scott and Erik K. Antonsson
Preliminary Vehicle Structure Design Application: Michael J. Scott and Erik K. Antonsson
Proceedings of the
10th International Conference on Design Theory and Methodology
ASME, Paper Number DETC98/DTM-5646, (September, 1998)
Abstract The Method of Imprecision, or MoI, is a formal method for incorporating im-
precise information into a design process. This methodology has been exercised
on a problem in preliminary vehicle structure design in collaboration with VW
Wolfsburg. Results show that the method is useful in trading off multiple con-
flicting attributes, including styling preferences and engineering requirements.
Keywords:
Industrial Applications of DTM; Vehicle Structure Design; Design Meth-
ods and Models; Design Representations; Computational Methods of Design;
Fuzzy Sets
Introduction
Preliminary design is inherently imprecise [3, 4, 12, 40], and many prelim-
inary design decisions are made informally. Preliminary design has enormous
economic importance, as much of the cost of a design is determined by these
(often informal) preliminary decisions [37]. A further complication is the dif-
ficulty of communicating imprecise information between different members or
groups involved in the design process. Many “interface” decisions are made
after design analysis is complete; these post hoc decisions can result in costly
redesigns.
The Method of Imprecision, or MoI [39, 1] has been developed to formally
incorporate imprecise information into the engineering design process. In the
summer of 1997, an application of the MoI to preliminary vehicle structure
183
184 IMPRECISION IN ENGINEERING DESIGN
µp
1
crisp
fuzzy
0
500 km range
~ and µp (~
(where µd (d) p) are themselves aggregations of their constituent pref-
erences), and candidate designs can be compared on the basis of this overall
preference.
Further research on the MoI developed techniques for including noise [26]
and adjustments, or tuning parameters [25], in the imprecision calculations,
and placed an axiomatic framework on the calculations [23]. Implementa-
186 IMPRECISION IN ENGINEERING DESIGN
tion of the MoI continued with the development of a computational tool [19].
The applicability of the method was seen to be limited by large computational
requirements, so the inclusion of Design of Experiments (DOE) approxima-
tions [21] and other computational innovations [18] followed. The need for
more than two aggregation functions to model different trade-off levels was
recognized, and a family of aggregations introduced [32, 33]. The MoI is re-
viewed in more detail, and compared to other methods, in [1].
Other researchers have applied fuzzy methods to design optimization prob-
lems [7, 8, 28, 29, 6]. Related research includes: chemical process synthesis
[9]; fuzzy constraint propagation applied to manufacturing [10]; fuzzy schedul-
ing [11]; application of fuzzy methods to windturbine design [13]; multiob-
jective scheduling [14]; management of uncertain knowledge in engineering
design [15]; engineering design optimization [16]; imprecise calculations in
engineering design [5]; evaluation of design alternatives [17]; fuzzy MADM
methods in system design [22]; fuzzy evaluations [27]; multiobjective fuzzy
optimization [30, 34]; fuzzy ratings and utility analysis in preliminary design
evaluation of multiple attributes [35]; scheduling system design [36]; and fuzzy
multi-criteria decision making [44, 45, 43]. In addition, there is increasing in-
terest in related work in evolutionary algorithms [46], including the combina-
tion of evolutionary algorithms with the MoI [31].
1 The demonstration here was conducted using a subset of the design variables; the method can be applied
3. Weight (kg)
The stated design problem was to achieve 10% improvements over the refer-
ence model in the three measured performances. In addition, it was understood
that that the design must not be difficult to manufacture, and that this year’s
model should have a somewhat longer and sleeker look.
1.0 1.0
0.8 0.8
0.6 0.6
µ µ
0.4 0.4
0.2 0.2
0.0 0.0
2000 3000 4000 5000 6000 7000 8000
Bending Stiffness (N/mm) Torsional Stiffness (N-m/deg)
1.0
0.8
0.6
µ
0.4
0.2
0.0
120 140 160 180
Weight (kg)
ent “optimization” problem that includes the imprecise information that would
be left to the negotiation stage in traditional design.
The calculated performance requirements on bending stiffness, torsional
stiffness, and weight were originally expressed as targets of 10% improve-
ments over the reference model. As was discussed above, this is unrealistically,
and indeed unproductively, precise. In place of these hard targets, imprecise
performance requirements were specified with a linear interpolation between
two points. In the implementation of the MoI, it is common to name the cus-
tomer as the source of the performance preferences; in fact, it is more likely to
be a manager, perhaps informed by market research, serving as the customer’s
proxy. To specify these imprecise requirements, the manager must answer two
simple questions: “What is the lowest performance you can live with (where is
µ = 0)? What performance would satisfy you completely (where is µ = 1)?”
These bounds are clearly dependent on a number of factors, including the target
market and the performance of competitors’ products; we have found that en-
gineering managers can answer these two questions with little more effort than
is needed to settle on the initial crisp target. Figure 8.5 shows the imprecise
requirements on stiffnesses and weight.
Preliminary Vehicle Structure Design Application 191
1. The sheet steel for stamping the A-pillar is only available in certain in-
crements, so this plot is discrete rather than continuous. The manufac-
turing engineer has a higher preference for thinner sheets, since they are
easier to form; this is a design preference for manufacturability.
2. The B-pillar thickness is continuous and more complicated than the lin-
ear performance preferences. This preference does not indicate that the
physical B-pillar might be 1.113 or 1.114 mm thick; rather it means that
the designer knows that the finite element model is simplified, and that a
high number for B-pillar thickness means that more reinforcing features
will need to be added to the B-pillar. The designer would like to keep
the B-pillar as simple as possible.
3. The floor pan thickness is preferred thicker by the designer for ease of
attachments and for durability.
5. The design preference for B-pillar location comes from the stylists, and
captures the directive for a longer, sleeker look for this year’s model. It
has been specified differently from the other design preferences, using
α-cuts [1], so that the stylists have given a range of perfectly accept-
able values, a range of barely acceptable values, and a range of values
that fall in the middle. This method of specifying preferences can have
computational advantages.
1.0 1.0
0.8 0.8
0.6 0.6
µ µ
0.4 0.4
0.2 0.2
0.0 0.0
0.6 0.7 0.8 0.9 1.0 1.1 0.8 1.0 1.2 1.4 1.6
A-pillar thickness (mm) B-pillar thickness (mm)
1.0 1.0
0.8 0.8
0.6 0.6
µ µ
0.4 0.4
0.2 0.2
0.0 0.0
0.9 1.0 1.1 1.2 1.3 1.4 0.7 0.8 0.9 1.0 1.1 1.2 1.3
Floor thickness (mm) Floor rail thickness (mm)
1.0
0.8
0.6
µ
0.4
0.2
0.0
-50 0 50 100 150 200
B-pillar location (mm aft of nominal)
The computation of µo (d) ~ for a single design point d~ is limited by the finite
element stiffness calculation, which takes about a minute on a Sun Ultra1-
170MHz workstation; the calculations of weight and preference aggregation
are of negligible cost, regardless of the trade-off strategies employed. Even in
this relatively modest problem, where there are only five design dimensions,
an exhaustive calculation of preferences over the design space is prohibitively
expensive. The MoI exploits the structure of the problem to speed the search
for preferred solutions: the internal calculations linearize where possible, ef-
fectively reducing the dimension of the search space, and employ Powell’s
method to locate internal extrema [20].
Results
The design problem, including all imprecise information, was solved in
two different ways. First, in order to demonstrate the method, the finite el-
ement analysis was run 3125 (= 55 ) times to provide a coarse but complete
check of the entire design space. The point of peak overall preference of
µo = 0.50 was found at d~ = (1.0, 0.9, 0.9, 1.0, 50), where the design pref-
erences µd are (0.6, 1.0, 1.0, 0.5, 1.0); the stiffnesses and weight at this point
were p~ = (2832, 5836, 147), with preferences (0.23, 0.14, 0.62). The max-
imum achievable stiffnesses are 3365 N/mm (µp = 0.77) in bending and
6029 N-m/degree (µp = 0.25) in torsion, but the corresponding weight of
170 kg is unacceptable. Similarly, a weight of 144 kg (µp = 0.78) is achiev-
able, but stiffnesses drop to 2803 N/mm (µp = 0.20) and 5730 N-m/degree
(µp = 0.08). The combined overall preference µo also takes into account the
design preferences on style, manufacturability, and the like.
The power of the method lies not in an ability to find a single overall “best”
point, but in the information it contains of how the total combined preference
µo varies with each of the design variables. Although it is impossible to dis-
194 IMPRECISION IN ENGINEERING DESIGN
play all five dimensions varying at once, a tool was written that uses a com-
mercial package (Matlab) to display results interactively. Using the tool, the
designer can see the change in preference that would occur by varying each
design variable independently from a chosen beginning point. Results can be
seen on five simultaneous plots in two dimensions (see Figure 8.7), or on a
three-dimensional surface plot (see Figure 8.8) with the remaining design vari-
ables set to nominal values. The interpretation of these graphical results is
discussed in greater detail below.
Approximations
Naturally, the exhaustive evaluation of points in the design space would not
be performed on a real design problem. It was performed here only for com-
parison purposes. An approach that utilizes Design of Experiments (DOE) [2]
to approximate the finite element calculations for bending and torsional stiff-
nesses reached substantially similar results in only 21 runs (or approximately
20 minutes). The average difference (from the exhaustive evaluation) in bend-
ing stiffness was approximately 1%, with a maximum difference of less than
4%, while the average and maximum differences for torsional stiffness were
both less than 1%.
In some cases, the nonlinearities of the analysis function f will defeat a
linear or even polynomial approximation, but in many cases, such as the ex-
ample presented here, these simple approximations can drastically reduce the
required computation. Since precise answers are not required for preliminary
design, it is sensible to exploit approximation tools when possible. If more
computation can be justified, a more thorough calculation can be made.
Discussion
Engineering analysis usually requires some judgement on the part of the
designer. Unless a full-scale exact prototype is to be built and tested, the ac-
curacy of any calculated performance measure depends on the fidelity of the
model employed. Even when exact data are available for some attributes, fi-
nal decisions about a design incorporate other, unmodelled concerns, such as
manufacturing and styling.
The MoI constructs a model of the entire decision process, expressing the
calculated overall performance µo (d)~ as a function of the design variables. It
depends on many factors: the function f for calculating measurable perfor-
mances, the specification of design preferences µd and performance prefer-
ences µp , the weighting of these preferences, and the specification of trade-off
aggregations between attributes. A change in any of these will affect the shape
of the function µo in design space, and thus affect the decision. The analysis
f is here relatively expensive to compute, and changes in the finite element
Preliminary Vehicle Structure Design Application 195
0 0 0
0.7 0.8 0.9 1 1.1 0.9 1 1.1 1.2 0.8 0.9 1 1.1
A pillar thickness B pillar thickness Floor sill thickness
1 1
0.8 0.8
0.6 0.6
µo µo
0.4 0.4
0.2 0.2
0 0
1 1.2 1.4 −50 0 50 100 150
Floor pan thickness B pillar location
µ*o = 0.40159
1
0.8
0.6
µo
0.4
0.2
0
1.4
1.3 150
1.2 100
50
1.1
0
1 −50
Floor Thickness B pillar location
model are costly to propagate to overall preference. Changes in the other fac-
~
tors, on the other hand, are easily incorporated, as finite element results f (d)
are stored so that the same design point need never be analyzed more than once.
This allows the MoI to support an iterative decision process, when the informa-
tion from the first round of calculations inspires a change in the preference
structure.
In the example shown in this paper, the shape of the function µo (d)~ is sensi-
tive to changes in the styling preference µd (d5 ), which is not surprising, since
this preference is accorded a large weight. This and other features of the de-
sign problem can be seen in the advanced interface shown in Figure 8.7. The
vertical dashed (red) lines indicate the selected values of the design variables.
In this figure, d~ = (0.8, 1, 0.9, 1.2, 50), a point which is representative, not
optimal. The solid (blue) lines indicate how the overall preference µo would
change by varying that design variable while holding the other four fixed at
their current values. For instance, decreasing the floor pan thickness d4 will
result in a more preferred design. The dashed (black) lines show the maximum
achievable µo for each value of each design variable. Finally, the solid (red)
lines joining the circles are the specified preferences on design variables.
In this example, the desired improvements in performance were achieved
by small changes to the sheet-metal thicknesses and B-pillar location. Visually
the improved structure would appear quite similar to that shown in Figure 8.3.
No change in vehicle structure configuration was required here.
Figure 8.7 shows that the overall preference µo varies qualitatively, though
not quantitatively, with four of the five design variables; the exception is d4 ,
floor pan thickness. The variation of the overall preference µo with respect to
d4 shows a conflict between the calculated stiffness and weight requirements
µp and the provisions for attachment and durability captured by the designer
preference µd (d4 ). This indicates d4 as a likely candidate for change in a
potential redesign. The resolution of the conflict would be achieved by a choice
of d4 that provides the best overall trade-off between the competing attributes.
Designers can interact with the preference display to examine trends in the
structure of the overall preference.
Conclusions
In preliminary vehicle structure design, as in preliminary engineering design
in general, many important decisions are made informally on the basis of im-
precise information. Concerns of styling and manufacturability, for instance,
can carry great weight in the design process although they are not modelled by
any formal analysis. The MoI is a tool to formally incorporate such imprecise
information into the design process, and thus to make decisions on a sound
basis. In a demonstration of the MoI prepared for VW Wolfsburg, concerns of
198 IMPRECISION IN ENGINEERING DESIGN
Acknowledgements
Volkswagen Wolfsburg, and particularly Mr. W. Lange, and Dr. O. Tegel
(presently Chief Engineer of the Institut für Maschinenkonstruktion at the Tech-
nical University in Berlin), were indispensable in the development of this project.
The authors would like express their appreciation to Zee Khoo and Juan
Nuño for their many hours of work stripping and load-testing the Rabbit, and
building the solid model, and to Mike Gerfen, manager of Caltech’s Central
Engineering Shop, for his generous help in providing space and equipment.
This material is based upon work supported, in part, by the National Science
Foundation under NSF Grant Number DMI-9523232. Any opinions, findings,
conclusions, or recommendations expressed in this publication are those of the
authors and do not necessarily reflect the views of the sponsors.
Preliminary Vehicle Structure Design Application 199
References
[1] A NTONSSON , E. K., AND OTTO , K. N. Imprecision in Engineering De-
sign. ASME Journal of Mechanical Design 117(B) (Special Combined
Issue of the Transactions of the ASME commemorating the 50th anniver-
sary of the Design Engineering Division of the ASME.) (June 1995), 25–
32. Invited paper.
[2] BARKER , T. Quality by Experimental Design. Marcel Dekker, Inc., New
York, 1985.
[3] B ECKER , J. M. A Structural Design Process Philosophy and Methodol-
ogy. PhD thesis, University of California, Berkeley, 1973.
[4] B LOCKLEY, D. I. The Nature of Structural Design and Safety. Ellis
Norwood Limited, Chichester, 1980.
[5] C HEN , J. E., AND OTTO , K. N. A tool for imprecise calculations in
engineering design. In Proceedings of the Fourth IEEE International
Conference on Fuzzy Systems (FUZZ-IEEE/IFES’95) (Mar. 1995), vol. 1,
IEEE, pp. 389–394.
[6] D HINGRA , A. K. Multiple objective design optimization. In Fuzzy Sets
in Engineering Design and Configuration, H.-J. Sebastian and E. K. An-
tonsson, Eds. Kluwer, 1996, pp. 53–88.
[7] D IAZ , A. R. Fuzzy set based models in design optimization. In Advances
in Design Automation - 1988 (New York, NY, Sept. 1988), S. S. Rao, Ed.,
vol. DE-14, ASME, pp. 477–485.
[8] D IAZ , A. R. A strategy for optimal design of hierarchical systems us-
ing fuzzy sets. In The 1989 NSF Engineering Design Research Con-
ference (College of Engineering, University of Massachusetts, Amherst,
June 1989), J. R. Dixon, Ed., NSF, pp. 537–547.
[9] D JOUAD , S., F LOQUET, P., D OMENECH , S., AND P IBOULEAU , L.
Fuzzy set methods in separation process synthesis. In Fuzzy Informa-
tion Engineering: A Guided Tour of Applications, D. Dubois, H. Prade,
and R. R. Yager, Eds. John Wiley & Sons, 1996.
[10] D UBOIS , D. Tuning the cutting edge conditions on a machining transfer
line using fuzzy constraint propagation. In Fuzzy Information Engineer-
ing: A Guided Tour of Applications, D. Dubois, H. Prade, and R. R. Yager,
Eds. John Wiley & Sons, 1996.
[11] FARGIER , H. Fuzzy scheduling: Principles and experiments. In Fuzzy
Information Engineering: A Guided Tour of Applications, D. Dubois,
H. Prade, and R. R. Yager, Eds. John Wiley & Sons, 1996.
[12] G AVIN , J R ., J. G. Fly me to the Moon: An interview with Joseph G.
Gavin, Jr. Technology Review 97, 5 (July 1994), 61–68.
200 IMPRECISION IN ENGINEERING DESIGN
Fitting to y = mx + c:
m = 4960.74 N-m/deg
c = −10.16 N-m
Fitting to y = mx + 0:
m = 4917.04N-m/deg
1800.00
1600.00
1400.00
1200.00
moment (N-m)
1000.00
800.00
200.00
0.00
0.00000 0.05000 0.10000 0.15000 0.20000 0.25000 0.30000 0.35000
twist (deg)
Bending
Load (N) Deflection (mm)
0.00 0.00000
284.67 0.101060
551.55 0.17780
836.22 0.33020
2001.60 0.78740
2286.27 0.91440
2837.82 1.16840
Fitting to y = mx + c:
m = 2426.16N/mm
c = 51.79N
Fitting to y = mx + 0:
m = 2484.80N/mm
3500.00
3000.00
2500.00
2000.00
load (N)
1500.00
data points
1000.00 fitted line
500.00
0.00
0.00000 0.20000 0.40000 0.60000 0.80000 1.00000 1.20000 1.40000
deflection (mm)