一种基于非线性接触关系的新型二维轮轨耦合单元及OpenSees实现
基于轮轨竖向非线性接触关系，提出一种新型通用的二维轮轨耦合单元模型，并在有限元OpenSees软件平台上实现。所提单元由轮节点和与之接触的梁单元节点组成，通过建立和求解轮轨作用力的一元三次方程，得到轮轨之间的接触力，继而计算由于轮轨相互作用产生的耦合单元各节点力，定义为单元内力，其中不包括梁单元内力。通过与文献中计算结果对比，验证了该单元模型的准确性和可靠性，并基于此模型分析了高速列车通过桥梁时在轨道不平顺激励和地震作用下的动力响应。此耦合单元模型易于集成到有限元计算平台中，和已有的列车模型、轨道和桥梁等模型联合使用，能够考虑轨道不平顺和轮轨脱离等情况，可用来分析复杂的竖向车桥耦合系统的动力问题。
由图5所示，基于本文轮轨耦合滑动单元计算所得的质块竖向位移和简支梁跨中挠度时程曲线与文献3阶模态解析解（标注为解析解）和曾庆元（标注为数值解）数值解高度吻合；质块的加速度时程曲线与文献数值解在趋势和数值上均吻合的很好，但是在末段，本文耦合单元模型数值解结果与3阶数值解有一定偏差，考虑到解析解仅考虑了3阶模态，其误差几乎是可以忽略，完全满足工程需求。
基于轮轨非线性接触关系提出了一种新型二维轮轨耦合单元。所提单元由轮节点和与之接触的梁单元节点组成，通过建立和求解轮轨作用力的一元三次方程，得到轮轨之间的接触力，继而计算由于轮轨相互作用产生的耦合单元各节点力，定义为单元内力，其中不包括梁单元内力。由于此耦合单元采用有别于传统有限元概念的单元内力和位移的定义方式，大大降低了加入现有的通用有限元框架的难度，使得该单元非常实用且能够方便的加入现有的通用有限元框架。并体现与其他方法不具备的优势：（1）能够同时模拟弹性接触和密贴接触两种接触关系，即可考虑车辆跳轨情况；（2）能够分析轨道不平顺激励的影响。
