Understanding the optimal designs in nature is critical in bionics. This paper presents a method for designing the configuration of fusiform muscle with a maximum contractile displacement based on topology optimization methods. A nearly incompressible continuum constitutive model of skeletal muscle is utilized. The contractile displacement from the relaxed state to the contracted state is regarded as the objective function. To handle the numerical difficulties that result from the existence of element density, an energy interpolation equation is employed, and a modification of the constitutive model of skeletal muscle is proposed. Several numerical examples are given to demonstrate the reasonability of the proposed method.