The objective of this study was to optimize the crawl stroke by the computer simulation. For this objective, the Particle Swarm Optimization was incorporated into our swimming human simulation model SWUM. With respect to the constrained condition, the maximum joint torque characteristics for the shoulder and elbow joints were imposed. In order to acquire these characteristics, the musculoskeletal simulation for the upper limb was employed. From the result of optimization, it was found that the optimized stoke was close to the so-called S-shape stroke.