The inherent low stiffness in soft robots makes them preferable for working in close proximity to humans. However, having this low stiffness creates challenges when operating in terms of control and sensitivity to disturbances. To alleviate this issue, soft robots often have built-in stiffness tuning mechanisms that allow for controlled increases in stiffness. Additionally, redundant pneumatic manipulators can utilize antagonistic pressure to achieve identical positions under increased stiffness. In this paper, we develop a model to predict the stiffness and configuration of a pneumatic soft manipulator under different pressure inputs and external forces. The model is developed based on the physical characteristics of a soft manipulator while enabling efficient parameter estimation and computation. The efficacy of the modeling approach is supported via experimental results.