介绍了常用的数值解法及接触算法的运算过程,详细地分析了各解法的优缺点。对目前在应力迭代算法中显式和隐式2种算法的适用范围进行了分析,指出在钻柱系统动力学研究中应该采用隐式算法。重点对阻尼模型和钻头岩石互作用模型进行了研究,借鉴海洋波浪力学中的成熟理论解决了钻柱在钻进过程中受到的钻井液阻尼问题,并建立了简化的钻头-岩石互作用模型。实例数值模拟计算表明,采用的方法可以对复杂井眼轨迹的1000m钻柱进行动态工作行为分析。研究结果对石油工程中采用数值方法分析超细长杆管柱在复杂井身结构及拟卖工况下的工作行为具有一定的指导意义。