Abstract: Contact analysis and stability of squirm analysis were discussed in terms of bellows with internal pressure load by building axis-symmetrical plane and three-dimensional solid model.Numerical experimental comparisons were conducted with regard to the specific values of friction coefficient,contact stiffness and ply distance related to contact effect among plies.At last,initial defection in nonlinear squirm analysis and pressure value needed to load to model in elastic-plastic analysis was described and the FEM analysis on multi-plies bellows was discussed aimed at providing some guidance for the development and design of bellows.