Problems with uniqueness and high parametric sensitivity of the solution of equations of motion, encountered in the static friction regime, are addressed. Friction in joints of a multibody system with closed-loop kinematic chains is discussed. Three different models of friction are studied: the discontinuous Coulomb model with stiction regime represented in terms of additional constraints; the approximate Coulomb model, smoothed in the vicinity of zero relative velocity; and the LuGre model with presliding displacements represented in terms of auxiliary state variables. Firstly, a rigid body model is investigated. It is shown that in the case of constraint addition approach, problems with uniqueness of solution emerge in the static friction regime. In the case of continuous models of friction, the solution in the stiction regime and its vicinity is highly sensitive to some hardly measurable or arbitrarily chosen parameters of the model of friction. Origins of nonuniqueness and high sensitivity are investigated, and the questionable credibility of the stiction regime simulation results is discussed. Secondly, a simplified model of body and joint elasticity is introduced to investigate the impact of flexibility on the mechanism frictional behavior. It is shown that taking the flexibility into consideration may eliminate the uniqueness and sensitivity problems. Moreover, the quantities that represent flexibility may be regarded as the key factors influencing the results of stiction regime simulation. Five examples are provided to illustrate the presented considerations.