1.一种桥上纵连板式无砟轨道疲劳应力谱分析方法，其特征在于，包括：
步骤(1)、根据气象局气象数据，结合桥上纵连板式无砟轨道温度场计算模型进行服役期间桥上纵连板式无砟轨道温度场的计算，得到服役期间时间-桥上纵连板式无砟轨道温度及温度梯度荷载时程曲线；
步骤(2)、采用APDL二次开发语言，在ANSYS通用有限元软件环境下开发高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型，采用预处理共轭梯度法求解耦合系统振动方程，在普通微机上进行高速列车-纵连板式无砟轨道-桥梁耦合系统动力仿真计算，结合步骤(1)中服役期间时间-桥上纵连板式无砟轨道温度梯度荷载时程曲线的计算结果，得到服役期间温度梯度荷载下桥上纵连板式无砟轨道各部件钢筋与混凝土应力时程曲线；
步骤(3)、采用APDL二次开发语言，在ANSYS通用有限元软件环境下开发考虑无砟轨道混凝土开裂和闭合效应及钢筋与混凝土相互作用的纵连板式无砟轨道-桥梁-桥梁墩台纵向相互作用模型，用该纵连板式无砟轨道-桥梁-桥梁墩台纵向相互作用模型进行纵向力组合作用下桥上纵连板式无砟轨道应力时程曲线计算，结合步骤(1)中服役期间时间-桥上纵连板式无砟轨道温度荷载时程曲线的计算结果，得到服役期间纵向荷载作用下桥上纵连板式无砟轨道各部件钢筋与混凝土应力时程曲线；
步骤(4)、用高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型进行列车自重荷载及轨道随机不平顺荷载下桥上纵连板式无砟轨道各部件动应力时程曲线的计算，结合服役期间时间-列车曲线，得到服役期间列车荷载作用下桥上纵连板式无砟轨道各部件钢筋与混凝土动应力时程曲线；
步骤(5)、将步骤(2)、(3)和(4)中的应力时程曲线相加，并考虑预应力的影响，得到服役期间组合荷载下桥上纵连板式无砟轨道各部件钢筋与混凝土应力时程曲线；
步骤(6)、用雨流计数法对步骤(5)中的应力时程曲线进行计数，得到服役期间组合荷载下桥上纵连板式无砟轨道各部件钢筋与混凝土二维疲劳应力谱；
步骤(7)、用Goodman曲线考虑平均应力的影响，对既包含幅值又包含均值的桥上纵连板式无砟轨道各部件钢筋与混凝土二维疲劳应力谱进行转换，得到仅含幅值的桥上纵连板式无砟轨道各部件钢筋与混凝土一维疲劳应力谱。
2.根据权利要求1所述的桥上纵连板式无砟轨道疲劳应力谱分析方法，其特征在于，所述步骤(1)具体包括：
步骤(11)、从气象局网站下载某一地区气象局相关气象资料，包括年、月、日、平均风速、平均气温、日照时数、日最低气温、日最高气温、总辐射日总量；
步骤(12)、借鉴国内外工程结构物温度场研究成果，建立桥上纵连板式无砟轨道温度场计算模型；
步骤(13)、根据气象局的气象资料，在MATLAB环境下进行二次开发，得到桥上纵连板式无砟轨道温度场计算模型边界条件，并写入一文本文件；
步骤(14)、在ANSYS环境下用APDL进行二次开发，编制桥上纵连板式无砟轨道温度场计算程序；计算程序读入MATLAB生成的桥上纵连板式无砟轨道温度场计算模型边界条件并求解，得到服役期间桥上纵连板式无砟轨道温度场；
步骤(15)、根据服役期间桥上纵连板式无砟轨道温度场的计算结果，进一步得到服役期间桥上纵连板式无砟轨道温度及温度梯度时程曲线。
3.根据权利要求1所述的桥上纵连板式无砟轨道疲劳应力谱分析方法，其特征在于，所述步骤(2)采用APDL二次开发语言在ANSYS通用有限元软件环境下开发多尺度高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型；中间桥梁及其上的无砟轨道采用精细化的建模技术，以提高桥上纵连板式无砟轨道疲劳应力谱的仿真精度，而边上桥梁及其上的无砟轨道采用比较粗糙的建模技术，以减少系统自由度，提高桥上纵连板式无砟轨道疲劳应力谱的求解速度。
4.根据权利要求3所述的桥上纵连板式无砟轨道疲劳应力谱分析方法，其特征在于，所述高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型具体建模如下：
高速列车由多节机车和车辆编组而成，机车和车辆采用在国内外得到广泛采用的10个自由度的多刚体模型，由车体、转向架、轮对、一系悬挂和二系悬挂组成；
钢轨选用梁单元进行模拟；
扣件采用弹簧-阻尼单元进行模拟；中间桥跨范围内的扣件考虑扣件的尺寸效应，将每个钢轨节点与其对应扣件尺寸范围内的轨道板节点均相连，以提高仿真精度；边跨桥梁范围内及桥梁两侧的扣件不考虑扣件的尺寸效应，以减少系统自由度，提高模型的求解速度；
轨道板和底座板采用实体单元进行模拟；中间桥梁范围其上的底座板和轨道板采用计算速度慢但精度高的20节点solid95实体单元，以提高仿真精度；边跨桥梁范围内的底座板和轨道板采用计算速度快但精度低的8节点solid45实体单元，以提高模型的求解速度；
CA砂浆采用弹簧-阻尼单元进行模拟，CA砂浆的刚度和阻尼特性根据CA砂浆的材料参数和厚度参数换算得到；
底座板与梁面之间设置的滑动层采用罚函数接触进行模拟，桥梁接缝区域的固结机构和高强度挤塑板采用线性接触刚度进行模拟；
桥梁采用实体单元进行模拟，将桥梁断面简化为矩形断面；
以竖向弹簧单元连接裂缝两侧轨道板及底座；
车轮与钢轨之间采用移动的考虑轨道不平顺影响的接触单元加以模拟。
5.根据权利要求1所述的桥上纵连板式无砟轨道疲劳应力谱分析方法，其特征在于，所述纵连板式无砟轨道-桥梁-桥梁墩台纵向相互作用模型具体建模如下：
钢轨选用梁单元进行模拟；
扣件采用非线性弹簧-阻尼单元进行模拟；
无砟轨道混凝土以梁单元进行模拟；
无砟轨道钢筋以梁单元进行模拟；
桥梁以梁单元进行模拟；
模拟无砟轨道混凝土与无砟轨道钢筋粘结滑移关系的连接单元以非线性弹簧-阻尼单元模拟；
模拟桥梁与无砟轨道混凝土间纵向相互作用的纵向连接单元以非线性弹簧-阻尼单元模拟；
在桥梁固定支座处，桥梁与无砟轨道混凝土间还有锚筋连接，以一纵向弹簧单元对此情况加以模拟；
固定支座处设置纵向线性弹簧-阻尼单元，以模拟桥梁墩台纵向刚度对系统纵向力的影响；
为了模拟无砟轨道混凝土在循环温度荷载下裂缝的开裂及闭合效应，以接触单元连接裂缝处两相邻的无砟轨道混凝土单元。
6.一种桥上纵连板式无砟轨道疲劳应力谱分析系统，其特征在于，包括：
第一处理单元、用于根据气象局气象数据，结合桥上纵连板式无砟轨道温度场计算模型进行服役期间桥上纵连板式无砟轨道温度场的计算，得到服役期间时间-桥上纵连板式无砟轨道温度及温度梯度荷载时程曲线；
第二处理单元、用于采用APDL二次开发语言，在ANSYS通用有限元软件环境下开发高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型，采用预处理共轭梯度法求解耦合系统振动方程，在普通微机上进行高速列车-纵连板式无砟轨道-桥梁耦合系统动力仿真计算，结合第一处理单元所得的服役期间时间-桥上纵连板式无砟轨道温度梯度荷载时程曲线的计算结果，得到服役期间温度梯度荷载下桥上纵连板式无砟轨道各部件钢筋与混凝土应力时程曲线；
第三处理单元、用于采用APDL二次开发语言，在ANSYS通用有限元软件环境下开发考虑无砟轨道混凝土开裂和闭合效应及钢筋与混凝土相互作用的纵连板式无砟轨道-桥梁-桥梁墩台纵向相互作用模型，用该纵连板式无砟轨道-桥梁-桥梁墩台纵向相互作用模型进行纵向力组合作用下桥上纵连板式无砟轨道应力时程曲线计算，结合第一处理单元所得的服役期间时间-桥上纵连板式无砟轨道温度荷载时程曲线的计算结果，得到服役期间纵向荷载作用下桥上纵连板式无砟轨道各部件钢筋与混凝土应力时程曲线；
第四处理单元、用于用高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型进行列车自重荷载及轨道随机不平顺荷载下桥上纵连板式无砟轨道各部件动应力时程曲线的计算，结合服役期间时间-列车曲线，得到服役期间列车荷载作用下桥上纵连板式无砟轨道各部件钢筋与混凝土动应力时程曲线；
第五处理单元、用于将第二、三、四处理单元所得的应力时程曲线相加，并考虑预应力的影响，得到服役期间组合荷载下桥上纵连板式无砟轨道各部件钢筋与混凝土应力时程曲线；
第六处理单元、用于用雨流计数法对第五处理单元所得的应力时程曲线进行计数，得到服役期间组合荷载下桥上纵连板式无砟轨道各部件钢筋与混凝土二维疲劳应力谱；
第七处理单元、用于用Goodman曲线考虑平均应力的影响，对既包含幅值又包含均值的桥上纵连板式无砟轨道各部件钢筋与混凝土二维疲劳应力谱进行转换，得到仅含幅值的桥上纵连板式无砟轨道各部件钢筋与混凝土一维疲劳应力谱。
7.根据权利要求6所述的桥上纵连板式无砟轨道疲劳应力谱分析系统，其特征在于，所述第一处理单元具体用于：
从气象局网站下载某一地区气象局相关气象资料，包括年、月、日、平均风速、平均气温、日照时数、日最低气温、日最高气温、总辐射日总量；
借鉴国内外工程结构物温度场研究成果，建立桥上纵连板式无砟轨道温度场计算模型；
根据气象局的气象资料，在MATLAB环境下进行二次开发，得到桥上纵连板式无砟轨道温度场计算模型边界条件，并写入一文本文件；
在ANSYS环境下用APDL进行二次开发，编制桥上纵连板式无砟轨道温度场计算程序；计算程序读入MATLAB生成的桥上纵连板式无砟轨道温度场计算模型边界条件并求解，得到服役期间桥上纵连板式无砟轨道温度场；
根据服役期间桥上纵连板式无砟轨道温度场的计算结果，进一步得到服役期间桥上纵连板式无砟轨道温度及温度梯度时程曲线。
8.根据权利要求6所述的桥上纵连板式无砟轨道疲劳应力谱分析系统，其特征在于，所述第二处理单元具体采用APDL二次开发语言在ANSYS通用有限元软件环境下开发多尺度高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型；中间桥梁及其上的无砟轨道采用精细化的建模技术，以提高桥上纵连板式无砟轨道疲劳应力谱的仿真精度，而边上桥梁及其上的无砟轨道采用比较粗糙的建模技术，以减少系统自由度，提高桥上纵连板式无砟轨道疲劳应力谱的求解速度。
9.根据权利要求8所述的桥上纵连板式无砟轨道疲劳应力谱分析系统，其特征在于，所述高速列车-纵连板式无砟轨道-桥梁三维有限元耦合动力学模型具体建模如下：
高速列车由多节机车和车辆编组而成，机车和车辆采用在国内外得到广泛采用的10个自由度的多刚体模型，由车体、转向架、轮对、一系悬挂和二系悬挂组成；
钢轨选用梁单元进行模拟；
扣件采用弹簧-阻尼单元进行模拟；中间桥跨范围内的扣件考虑扣件的尺寸效应，将每个钢轨节点与其对应扣件尺寸范围内的轨道板节点均相连，以提高仿真精度；边跨桥梁范围内及桥梁两侧的扣件不考虑扣件的尺寸效应，以减少系统自由度，提高模型的求解速度；
轨道板和底座板采用实体单元进行模拟；中间桥梁范围其上的底座板和轨道板采用计算速度慢但精度高的20节点solid95实体单元，以提高仿真精度；边跨桥梁范围内的底座板和轨道板采用计算速度快但精度低的8节点solid45实体单元，以提高模型的求解速度；
CA砂浆采用弹簧-阻尼单元进行模拟，CA砂浆的刚度和阻尼特性根据CA砂浆的材料参数和厚度参数换算得到；
底座板与梁面之间设置的滑动层采用罚函数接触进行模拟，桥梁接缝区域的固结机构和高强度挤塑板采用线性接触刚度进行模拟；
桥梁采用实体单元进行模拟，将桥梁断面简化为矩形断面；
以竖向弹簧单元连接裂缝两侧轨道板及底座；
车轮与钢轨之间采用移动的考虑轨道不平顺影响的接触单元加以模拟。
10.根据权利要求6所述的桥上纵连板式无砟轨道疲劳应力谱分析系统，其特征在于，所述纵连板式无砟轨道-桥梁-桥梁墩台纵向相互作用模型具体建模如下：
钢轨选用梁单元进行模拟；
扣件采用非线性弹簧-阻尼单元进行模拟；
无砟轨道混凝土以梁单元进行模拟；
无砟轨道钢筋以梁单元进行模拟；
桥梁以梁单元进行模拟；
模拟无砟轨道混凝土与无砟轨道钢筋粘结滑移关系的连接单元以非线性弹簧-阻尼单元模拟；
模拟桥梁与无砟轨道混凝土间纵向相互作用的纵向连接单元以非线性弹簧-阻尼单元模拟；
在桥梁固定支座处，桥梁与无砟轨道混凝土间还有锚筋连接，以一纵向弹簧单元对此情况加以模拟；
固定支座处设置纵向线性弹簧-阻尼单元，以模拟桥梁墩台纵向刚度对系统纵向力的影响；
为了模拟无砟轨道混凝土在循环温度荷载下裂缝的开裂及闭合效应，以接触单元连接裂缝处两相邻的无砟轨道混凝土单元。