We apply a Hamiltonian particle method, which is one of the particle methods, to simulate seismic wave propagation in a cracked medium. In a particle method, traction-free boundaries can readily be implemented and the spatial resolution can be chosen in an arbitrary manner. Utilisation of the method enables us to simulate seismic wave propagation in a cracked medium and to estimate effective elastic properties derived from the wave phenomena. These features of the particle method would bring some advantages of numerical efficiencies (e.g. calculation time, computational memory) and the reduction of time for pre-processing.We first describe our strategy for the introduction of free surfaces inside a rock mass, i.e., cracks, and to refine the spatial resolution in an efficient way. We then model a 2D cracked medium which contains randomly distributed, randomly oriented, rectilinear, dry and nonintersecting cracks, and simulate the seismic wave propagation of P-and SV-plane waves through the region. We change the crack density in the cracked region and determine the effective velocity in the region. Our results show good agreement with the modified selfconsistent theory, which is one of the effective-medium theories. Finally, we investigate the influence of the ratio of crack length to particle spacing on the calculated effective velocities.The effective velocity obtained becomes almost constant when the ratio of crack length to particle spacing is more than approximately 20. Based on this result, we propose to use more than 20 particles per crack length.