---

This paper presents a design process using a course grained parallel genetic algorithm to optimize three-dimensional steel moment frames by considering the axial force and biaxial bending moments interaction in plastic hinge formation. The objective function is to minimize the total weight of the structure subjected to the reliability constraint of the structural system. System reliability analysis is performed through the proposed Modified Latin Hypercube Simulation (M-LHS) Method. For optimization, a 3DSMF-RBO program is written in CSHARP programming language. The reliability analysis results show a large decrease in the number of simulation samples and subsequently a decrease in the execution time of optimization computation. The optimization results indicate that by considering interaction of the axial force and biaxial bending moments in plastic hinge formation rather than the only bending moment, to some extent increases the total weight of the designed structure.