Robotics with flexible legs have attracted significant attention. Engineers often design and analyze the motion of legged robots from kinematic and biomimetic perspective. However, the influence of flexibility of the feet on robot locomotion is often not given sufficient considerations, which is also very crucial to the motion posture of the flexible legged robots, especially as the soft robots design becomes increasingly popular. The mainly difficulties lies in the traditional numerical methods in handling the dynamic motion analysis with both large rigid motion and large deformation. In this paper, the finite particle method (FPM) is used to simulate the motion and deformation coupled problems of the flexible six-legged robot. A shell-based particle model of a six-leg robot and the contact model between legs and ground are built. Without iterative and modification of the FPM analytical framework, structural nonlinearity is efficiently handled after eliminating rigid body motions by a fictitious reverse motion. The motion and deformation of a single leg with varying leg thickness, locomotion speed, and leg-to-ground friction coefficients were simulated. By analyzing the stress distribution in the leg and the number of contact points with the ground, the mechanical leg was optimized in design. Furthermore, the motion and deformation of the entire six-legged robot were simulated using the FPM. The numerical results' feasibility was validated through comparison with experimental data obtained from robot walking tests. The proposed method effectively simulates the motion and deformation of flexible robots, providing significant insights for the design of soft robots.