As shown in

**FIG. 9**, the processing of Step S

**6** is performed on a three-dimensional model (analytical model) of a structure built in a virtual three-dimensional space. The analytical model is constituted by arranging a plurality of three-dimensional elements (a tetrahedron, a hexahedron, or the like). Herein, the calculation point is placed at the center point of the three-dimensional element. In

**FIG. 9**, at the respective calculation points a to e on the analysis surface

**501**, inherent strain estimated values are given. A case where the inherent strain is transferred to a transfer surface

**502** spaced from the analysis surface

**501** by θ around a z axis which is the rotation center is considered. The CPU

**111** transfers the respective estimated inherent strain values at the calculation points a to e to calculation points a′ to e′ corresponding to the calculation points a to e on the transfer surface. That is, the inherent strain values at the calculation points a′ to e′ are the same as the respective inherent strain values at the calculation points a to e. The CPU

**111** performs such processing over the transfer range (in a case of the journal shaft, in a range of 360°, and in a case of the pin shaft, in a range of 180°).