The contact interaction between compliant mechanisms and external objects is a common occurrence in the field of grasping and manipulation. However, accurately modeling the large deformation and contact force/stress remains a challenging task. This paper presents a contact analysis method to evaluate the deformation and normal contact force/stress of compliant mechanisms with general-shape beams. The general beams are modeled by using the chained pseudo-rigid-body model (CPRBM). On basis of CPRBM framework, linear and distributed contact springs are introduced at the interface to represent the contact force. A numerical method is established to fast calculate the boundary penetrations, and a contact stiffness formulation is proposed based on the theory of elasticity. Then, the effect of contact interaction can be evaluated by the contact potential energy of contact springs, which is part of the system’s total potential energy. By minimizing the potential energy function of the compliant system, the static equilibrium configuration and contact force/stress can be obtained easily. Numerical examples are used to verify the feasibility and accuracy of the proposed method.