A unilateral impact double pendulum model with hinge links is constructed to detect subharmonic bifurcation for the high dimensional non-smooth system. The non-smooth and nonlinear coupled factors lead a barrier for high dimensional conventional nonlinear techniques. By introducing reversible transformation and energy time scale transformation, the system is expressed as a smooth decoupling form of energy coordinates. Thus, the concept of subharmonic Melnikov function is extended to high-dimensional nonsmooth systems, and the influence of impact recovery coefficient on the existence of subharmonic periodic orbits of double pendulum is revealed. The efficiency of the theoretical results is verified by phase portraits, time process portraits and Poincaré section.