At present, the time of easy oil and gas is over. Now, the largest part of fossil fuels is concentrated in the deepest levels of tectonic structures and in the sea shelves. One of the most cumbersome operations of their extraction is the bore-hole drilling. In connection with austere tectonic and climate conditions, their drivage every so often is associated with great and diversified technological difficulties causing emergencies on frequent occasions. As a rule, they are linked with drill string accidents. A key role in prediction of these situations should play methods of theoretical modelling. For this reason, there is a growing need for development and implementation of new numerical methods for computer simulation of critical and post-critical behavior of drill strings (DSs). In this paper, the processes of non-linear deforming of a DS in cylindrical cavity of a deep bore-hole are considered. On the basis of the theory of curvilinear flexible rods, non-linear constitutive differential equations are deduced. The effects of the longitudinal non-uniform preloading, action of torque and interaction between the DS and the bore-hole surface are taken into account. Owing to the use of curvilinear coordinates in the constraining cylindrical surface and a specially chosen concomitant reference frame, it became possible to separate the desired variables and to reduce the total order of the equation system. To solve it, the method of continuation the solution by parameter and the transfer matrix technique are applied. As a result of the completed numerical analysis, the critical states of the DS loading in the cylindrical channels of inclined bore-holes are found. It is shown that the modes of the post-critical deforming of the DS are associated with its irregular spiral curving prevailing in the zone of bottom-hole-assembly. The possibility of invariant state generation during post-critical deforming is established, condition of its bifurcation is formulated. It is shown that infinite variety of loads can correspond to one geometrical configuration of the DS. They differ each from other by contact force functions.