Pneumatic artificial muscles (PAMs) have been widely applied to robotic systems, especially assistive devices, which could benefit from PAMs’ intrinsic viscoelasticity. However, nonlinearity and hysteresis make it challenging to achieve high-accuracy control. Using pre-tensioned springs to replace one side PAMs in an antagonistic structure is one efficient way to improve control accuracy. However, in the previous study, the selection method for spring, especially the range of spring constant has not been sufficiently investigated. This study designed a method to find the optimal spring constant based on gradient descent algorithm; and verified it with prototype and simulation experiments. Results showed that the group with optimal spring constant leads to a faster response time, smaller overshoot, and lower steady-state error, which indicated that the proposed method could find optimal spring constant. Moreover, the sensitivity of PAMs dynamic model parameters to the optimal spring stiffness was investigated. This study provides insights to practical use of soft actuators.