In this paper, we investigate parallel structural optimization methods on distributed memory MIMD machines. We have restricted ourselves to the case of minimizing a multivariate non-linear function subject to bounds on the independent variables, when the objective function is expensive to evaluate as compared to the linear algebra portion of the optimization. This is the case in structural applications, when a large three-dimensional finite element mesh is used to model the structure.This paper demonstrates how parallelism can be exploited during the function and gradient computation as well as the optimization iterations. For the finite element analysis, a ‘torus wrap’ skyline solver is used. The reflective Newton method, which attempts to reduce the number of iterations at the expense of more linear algebra per iteration, is compared with the more conventional active set method. All code is developed for an Intel iPSC/860, but can be ported to other distributed memory machines.The methods developed are applied to problems in bone remodeling. In the area of biomechanics, optimization models can be used to predict changes in the distribution of material properties in bone due to the presence of an artificial implant. The model we have used minimizes a linear combination of the mass and strain energy in the entire domain subject to bounds on the densities in each finite element.Early results show that the reflective Newton method can outperform active set methods when few variables are active at the minimum.
[1]
Thomas F. Coleman,et al.
A Reflective Newton Method for Minimizing a Quadratic Function Subject to Bounds on Some of the Variables
,
1992,
SIAM J. Optim..
[2]
A. Pothen,et al.
A Distributed Multifrontal Algorithm Using Clique Trees
,
1991
.
[3]
Thomas F. Coleman,et al.
On the convergence of interior-reflective Newton methods for nonlinear minimization subject to bounds
,
1994,
Math. Program..
[4]
Thomas H. Dunigan,et al.
Performance of the Intel iPSC/860 and Ncube 6400 hypercubes
,
1991,
Parallel Comput..
[5]
Robert A. van de Geijn.
Massively Parallel Linpack Benchmark on the Intel Touchstone Delta andIPSC/860 Systems (Progress Report)
,
1991
.
[6]
T. Coleman,et al.
On the Convergence of Reflective Newton Methods for Large-scale Nonlinear Minimization Subject to Bounds
,
1992
.
[7]
S. Nash,et al.
Assessing a search direction within a truncated-newton method
,
1990
.
[8]
G. A. Geist,et al.
A user's guide to PICL a portable instrumented communication library
,
1990
.
[9]
K G Heiple,et al.
Bone strength. The effect of screw holes.
,
1972,
The Journal of bone and joint surgery. American volume.
[10]
W. Hayes,et al.
The compressive behavior of bone as a two-phase porous structure.
,
1977,
The Journal of bone and joint surgery. American volume.