摘要

基于虚功原理,推导了三维接触面单元的刚度矩阵,并引进预留单元的概念,避免了接触过程中由于接触面变化节点和单元需要重新编号的麻烦.文中算例验证了计算方法的有效性