A. Tafreshi et Rt. Fenner, GENERAL-PURPOSE COMPUTER-PROGRAM FOR SHAPE OPTIMIZATION OF ENGINEERING STRUCTURES USING THE BOUNDARY-ELEMENT METHOD, Computers & structures, 56(5), 1995, pp. 713-720
Optimization of the shape of structural components is often concerned
with minimizing stress concentration effects to ensure a longer fatigu
e life. A general purpose computer program, STRESOPT, for shape optima
l design of two-dimensional structures in order to smooth stress peaks
is presented. The program has three main parts: the stress analyzer,
design sensitivity analyzer and optimizer. The stress analyzer and des
ign sensitivity analyzer use the boundary element method with isoparam
etric quadratic boundary elements. The boundary element method is very
suitable for shape optimization and in comparison with the finite ele
ment method needs much less data, related only to the boundary of the
structure being considered. Because a differentiated form of the bound
ary integral equation (on which the boundary element method is based)
can be used directly to determine the derivatives of the objective and
constraint functions, the accuracy of computation is very high. The n
umerical optimization method used in the program is the extended penal
ty function approach, using the BFGS variable metric for unconstrained
minimization, together with the Golden Section method for the one-dim
ensional search. Initial mesh preparation and regeneration of the boun
dary elements during the iterative process of optimization is both str
aightforward and fast. Hermitian cubic splines are well suited for the
boundary shape representation, and complex geometries can be describe
d in a very compact way by a small number of design variables. Applica
tions of the program to the optimum shape design of fillets and holes
in plates and bars are presented.