This research tackles topology optimization challenges with large deformation in hyperelastic materials. Utilizing nonlinear response characteristics and significant deformation capabilities, we introduced a novel sensitivity analysis method. large deformation finite element analysis was combined with an updating method using a reaction-diffusion equation (RDE) and SIMP (Solid Isotropic Material with Penalization) method. The optimization targeted the minimization of shape compliance. By assuming linear elasticity for sensitivity analysis to reduce computational cost, the optimal topology was determined, validating the effectiveness of the proposed method. This study contributes to the structural optimization design in hyperelastic materials and offers new paths for future applications.