本发明公开了一种钢轨表面三维几何不平顺处轮轨冲击力的模拟方法，包括如下步骤：获取钢轨表面三维几何不平顺分布数据；建立包含任意三维几何不平顺的三维钢轨几何型面；利用空间几何求解方法获得轮轨接触点和三维非对称接触几何间隙分布；建立三维非对称轮轨接触模型并计算轮轨接触解；建立车辆‑轨道耦合动力学模型，求解钢轨三维表面几何不平顺处轮轨冲击力。本发明通过离散车轮和钢轨的真实几何型面获得了三维几何不平顺处的接触点及几何间隙分布，并利用准确的三维接触模型进行求解得到接触力学行为，再带入车辆‑轨道耦合动力学中求得动态轮轨力，本发明方法能考虑三维几何不平顺，得到更为准确的计算结果。