This paper presents a numerical analysis technique for application to shape optimization problems of linear elastic structures subject to multiple loading conditions. The problems dealt with here are a mean compliance minimization problem in relation to individual load cases and a fully stressed design problem. The proposed technique is based on the traction method which analyzes the domain variation. A shape optimization system was developed and applied to fundamental problems in two and three dimensions. The computed results confirmed the validity and usefulness of the proposed technique.