SUMMARYLagrangian particle method has been widely used in computer physics and graphics; however, numerically solving the partial differential physical equation on a great number of particles is a computationally complex task. In this paper, a unified particle method on graphics processing units is proposed to simulate fluid-solid interaction with large density ratio interactively. Motivated by microscopic molecular dynamics, we consider the solid object as a particular fluid limited to solid motions; therefore, fluid-solid interaction as well as solid-solid interaction could be solved directly using multiphase weakly compressible smoothed particle hydrodynamics solvers. And then, we present a momentum-conserving particle collision handling scheme to prevent fluid penetrating into solid objects. In the simulation, a measure of particle densities is used to handle density discontinuities at fluid-solid interfaces, and consequently, new formulations for density-weighted inter-particle pressure and viscous forces are derived. Moreover, to realistically simulate various small-scale interaction phenomena such as water droplets flowing on solids' surfaces, a surface tension model that uses density-weighted color gradient and can obtain a stable and accurate surface curvature is employed to capture the interfacial fluid-solid tensions. Because all of the computation is carried out on graphics processing unit and no CPU processing is needed, the proposed algorithm can exploit the massive computational power of graphics processing unit for interactive simulation with a higher particle resolution. The experiment results show that our method can simulate realistic fluid-solid couplings at interactive frame-rates even for up to 126 k particles.