提出了一种有限元和虚拟关节法相结合的协作机器人非线性刚度建模方法,该方法不仅包含连杆和关节模块线性及耦合刚度,还考虑了机器人传动系统非线性刚度、负载作用下雅可比矩阵变化率以及负载和自重引起的位姿变化对机器人刚度非线性的影响。通过有限元子结构法提取机器人各模块的结构刚度矩阵,与传动系统非线性刚度模型相结合,建立机器人各模块的非线性综合刚度模型;采用数值计算方法计算在外力和自重工况下的静力平衡位姿及各模块受力情况;基于虚拟关节法构建机器人的非线性刚度模型,通过非线性有限元仿真验证了该建模方法的正确性。