本发明属于高速铁路无砟轨道领域，涉及一种有限元仿真分析方法，尤其涉及一种高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法。
背景技术：
高速铁路在当下是一个热门的话题，在国家大力倡导之下，我国的高速铁路工程也得到了大力发展，列车行驶在传统的有渣轨道上容易引起道砟飞溅、碎石被碾碎等病害，这些病害可能会导致路基沉降从而影响钢轨的不平顺性，轻则限制列车行驶速度，严重情况下会影响列车在高速行驶过程中的安全性和平稳性。综合上述缺点，无砟轨道技术逐步从国外引进，到现在我国有了具有自主知识产权的crtsiii型板式无砟轨道。目前，crtsiii型板式无砟轨道技术发展只历经短短几年时间，该领域的统一标准也未出台，因此亟待对crtsiii型板式无砟轨道技术进行优化和创新。
无砟轨道相对有砟轨道来说具有刚度大、沉降小等特点，列车行驶在刚度较大的无砟轨道结构上必然会引起剧烈振动，在何种振动频率以及加速度下会使得列车能够安全地行驶，在何种横垂向位移变化量下会使钢轨变得不平顺等等问题。解决这些问题通常可以归结为动力响应问题。要解决这些问题通常可以通过试验的方法，确定相关实测数据然后进行改进设计，但进行一次试验耗费的时间资金也是巨大的，最广为有效的方法是建立一套与实际情况相似的车桥三维实体耦合有限元模型。
目前，国内绝大多数研究单位对于车辆轨道耦合动力问题的研究都建立在空间动态耦合的状态下。考虑车辆系统为两系悬挂系统，充分考虑车辆的自由度；考虑轨道系统的各结构层，但现有的研究仍存在以下不足：对于车辆轨道耦合动力学的研究，绝大多数集中在对于车辆振动和桥梁结构动力响应的研究，忽略了无砟轨道结构自身的振动；对于无砟轨道结构模型进行大量的简化，典型的是忽略无砟轨道结构层厚度，采用梁单元或板壳单元对无砟轨道各结构层进行模拟，仅对无砟轨道表层动力响应进行研究，忽略其内部的动力响应状况；模型中还忽略了实际的无砟轨道内部配筋情况，现有的无砟轨道模型显然对其实际结构进行了大量的简化，使得无砟轨道计算结果与实际动力响应存在更大的偏差。综合以上不足，建立一套模拟列车在桥上无砟轨道高速行驶的列车-crtsiii型板式无砟轨道-桥梁结构空间耦合动力有限元模型显得尤为必要，其结果可以有效地指导国内crtsiii型板式无砟轨道设计。
技术实现要素：
为了解决背景技术中存在的上述技术问题，本发明提供了一种计算结果与实测数据对比发现吻合度较高、可有效地指导国内crtsiii型板式无砟轨道的设计、施工以及养护的高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法。
为了实现上述目的，本发明采用如下技术方案：
一种高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法，所述高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法包括以下步骤：
1)运用abaqus有限元软件建立高铁无砟轨道桥梁结构空间耦合动力有限元模型；
2)对步骤1)建立得到的高铁无砟轨道桥梁结构空间耦合动力有限元模型进行动力仿真分析。
作为优选，本发明所采用的步骤1)的具体实现方式是：
1.1)建立高速铁路车辆耦合模型；
1.2)建立crtsiii型板式无砟轨道模型；
1.3)建立桥梁模型；
1.4)定义接触方式；
1.5)定义荷载；
1.6)定义边界条件。
作为优选，本发明所采用的步骤1.1)的具体实现方式是：
所述高速铁路车辆耦合模型包括轮对、转向架以及车体；所述轮对、转向架以及车体均定义为刚体，轮对与转向架连接方式设置一系悬挂，转向架与车体之间设置二系悬挂，一系悬挂以及二系悬挂均考虑纵横垂向刚度和阻尼，车轮的踏面采用铁道部标准tb1967-87规定的lm型87磨耗型踏面；通过定义轮对、转向架和车体刚体约束耦合在一起；所述高速铁路车辆耦合模型的尺寸符合铁道部铁道部标准tb/t1010-2005规定；所述车辆耦合是通过在轮对、转向架和车体形心处建立约束参考点，然后与分别与自身进行刚体约束，通过mpc梁多点约束将车体、转向架以及轮对联系在一起。
作为优选，本发明所采用的步骤1.2)的具体实现方式是：
所述crtsiii型板式无砟轨道模型的结构自上而下分别为：钢轨、弹簧扣件、轨道板、自密实混凝土填充层、底座以及桥梁；所述桥梁上是单元式板式无砟轨道，轨道板与轨道板之间存在一定间距；所述crtsiii型板式无砟轨道模型采用abaqus中的实体单元对钢轨、轨道板、自密实混凝土填充层以及底座进行模拟，采用abaqus中的桁架单元对普通钢筋、预应力钢筋和u型钢筋进行模拟；所述crtsiii型板式无砟轨道模型中的各部件尺寸均与实际尺寸相同，根据实际的材料属性对crtsiii型板式无砟轨道模型各部件进行定义。
作为优选，本发明所采用的步骤1.3)的具体实现方式是：
通过abaqus中的实体单元来模拟铁路上常见的32米简支梁桥，桥梁模型尺寸与实际尺寸相同。
作为优选，本发明所采用的步骤1.4)的具体实现方式是：
轮轨接触方式采用表面与表面接触，主面为轮对表面，从面为钢轨顶面和内侧面；填充层与底座表面接触方式同样定义为表面与表面接触；轮轨接触包含切向接触和法向接触；切向采用罚函数设定摩擦系数来定义，所述罚函数的函数表达式是：
f＝[uk+(us-uk)e-βv]·p(t)
其中：
us为静摩擦系数；
uk为动摩擦系数；
β为衰减系数；
v为主从面的相对滑移速度；
t为时间；
p(t)表征为某一时刻对应的轮轨法向接触力；
所述切向采用通过公式轮轨间法向接触力采用“herz”非线性接触理论，即不同压力下所产生的过盈量，所需参数为接触压力和过盈量；普通钢筋、预应力钢筋和u型钢筋均采用embedded约束，其他各结构层之间接触方式均采用tie约束。
作为优选，本发明所采用的步骤1.5)的具体实现方式是：
荷载参数根据和谐号动车组列车质量参数来定义，并考虑轨道板中预应力荷载；预应力荷载的添加是对轨道板中预应力钢筋的模拟采用等效降温法来实现；
通过公式
作为优选，本发明所采用的步骤1.6)的具体实现方式是：
考虑对车体、轮对、转向架的转动约束，考虑钢轨两侧的竖向约束和底部的转动约束，根据实际简支梁桥支座的布置方式来对支座边界条件进行模拟。
作为优选，本发明所采用的步骤2)的具体实现方式是：
通过运行高铁无砟轨道桥梁结构空间耦合动力有限元模型计算得到车辆动力学安全性指标；所述车辆动力学安全性指标包括轮轨接触力、车体横垂向振动加速度、脱轨系数以及轮轴减载率；
通过运行高铁无砟轨道桥梁结构空间耦合动力有限元模型计算得到板式无砟轨道结构动力学性能指标；所述板式无砟轨道结构动力学性能指标包括钢轨横垂向振动加速度、钢轨横垂向位移、轨道板横垂向振动加速度、轨道板横垂向位移、自密实混凝土横垂向振动加速度、自密实混凝土横垂向位移、底座横垂向振动加速度以及底座横垂向位移；
通过运行高铁无砟轨道桥梁结构空间耦合动力有限元模型计算得到动力学分析结果；所述动力学分析结果包括桥梁的横垂向振动加速度、桥梁横垂向位移以及梁端转角。
本发明的优点是：
本发明提供了一种高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法，将传统的列车-轨道-桥梁耦合分析方法大多做了很多简化，如简化为梁、板、壳等结构，旨在改变现有技术所得结果并不能真实地反映列车高速行驶时无砟轨道和桥梁的状态。本发明建立了三维空间实体精细化耦合分析模型，将列车、crtsiii型板式无砟轨道和桥梁三者耦合在一起综合考虑，计算结果与实测数据对比发现吻合度较高。通过采用本发明建立的高铁无砟轨道桥梁结构空间耦合动力有限元模型，计算结果趋于准确化，模型也更加精细化，对于研究列车在高速行驶时crtsiii型板式无砟轨道和桥梁结构各细部区域的动力响应结果有一定的价值意义。采用本发明建立的高铁无砟轨道桥梁结构空间耦合动力有限元模型，其结构的参数和材料的参数可以依据铁路线路实际的结构和材料(如轨道板尺寸变化、预应力钢筋配筋率变化等)实际情况而变化，可以有效地指导国内crtsiii型板式无砟轨道的设计、施工以及养护。
具体实施方式
下面对本发明作更进一步的描述：
本发明提出了一种高铁无砟轨道桥梁结构空间耦合动力有限元仿真分析方法，包括以下步骤：
1)运用abaqus有限元软件建立高铁无砟轨道桥梁结构空间耦合动力有限元模型：
1.1)建立高速铁路车辆耦合模型，列车模型包括轮对、转向架和车体，均定义为刚体，通过定义轮对、转向架和车体刚体约束耦合在一起。轮对与转向架连接方式设置一系悬挂，转向架与车体之间设置二系悬挂，一系二系悬挂均考虑纵横垂向刚度和阻尼。
1.2)采用实体单元对钢轨进行模拟，钢轨的截面面积与实际钢轨横截面面积相等，对钢轨赋予密度、弹性模量、泊松比等材料参数。
1.3)对弹簧扣件的模拟采用abaqus中的“cartesian”进行模拟，弹簧扣件假定为线性的，考虑纵横垂三个方向的刚度和阻尼。
1.4)采用实体单元对iii型板的轨道板进行模拟，其几何尺寸来源于实际无砟轨道铁路线路上的轨道板尺寸，考虑密度、弹性模量、泊松比等材料参数。
1.5)采用实体单元对无砟轨道填充层自密实混凝土进行模拟，底部含有两个凸块对无砟轨道结构层进行限位，其几何尺寸来源于实际无砟轨道铁路线路上的自密实混凝土尺寸，考虑密度、弹性模量、泊松比等材料参数。
1.6)采用实体单元对底座进行模拟，底部含有两个凹槽，凹槽尺寸略比凸块大，其几何尺寸来源于实际无砟轨道铁路线路上的底座的尺寸，考虑密度、弹性模量、泊松比等材料参数。
1.7)采用桁架单元对无砟轨道中的钢筋进行模拟，钢筋包括轨道板预应力钢筋、轨道板u型钢筋和结构层普通钢筋，其中预应力的模拟要在边界条件中采用等效降温法赋予给预应力钢筋，考虑密度、弹性模量、泊松比等材料参数。
1.8)采用实体单元对轨下基础桥梁结构的模拟，桥梁采用简支梁桥模拟，对支座的模拟通过分割梁体底面实际的支座面积大小，然后添加参考点，将参考点与分割的面耦合在一起，通过对参考点的约束以此来达到对支座的模拟，考虑密度、弹性模量、泊松比等材料参数。
1.9)轮对与钢轨之间的接触方式定义为切向接触和法向接触，其中，切向接触采用采用滑动摩擦力予以表征，具体公式：f＝μ·p(t)；u＝uk+(us-uk)e-βv，公式中us为静摩擦系数，uk为动摩擦系数，β为衰减系数，v为主从面的相对滑移速度，在abaqus中输入静摩擦系数、动摩擦系数和衰减系数即可完成对切向接触的定义。法向接触采用赫兹非线性接触模型，定义法向接触时还要考虑压力过盈，及钢轨在不同压力下所允许发生的变形量。
1.10)通过对以上各部件装配然后定义约束和边界条件，即完成高铁无砟轨道桥梁结构空间耦合动力有限元模型的建立。
2)对列车一crtsiii型板式无砟轨道-桥梁结构空间耦合动力有限元模型进行动力仿真分析，具体是：
通过运行该模型计算得到轮轨接触力、车体横垂向振动加速度、脱轨系数、轮轴减载率等车辆动力学安全性指标；通过运行该模型计算得到钢轨横垂向振动加速度、钢轨横垂向位移、轨道板横垂向振动加速度、轨道板横垂向位移、自密实混凝土横垂向振动加速度、自密实混凝土横垂向位移、底座横垂向振动加速度、底座横垂向位移等iii型板式无砟轨道结构动力学性能指标；通过运行该模型计算得到桥梁的横垂向振动加速度、桥梁横垂向位移、梁端转角等动力学分析结果。
实施例1
本实施例通过模拟列车在时速200km/h、250km/h、300km/h、350km/h的速度下，运行模型，得到轮轨横垂向力的仿真结果对比如下表1。
表1模型在不同列车时速下轮轨力最大值的仿真结果对比
通过模拟列车在时速200km/h、250km/h、300km/h、350km/h的速度下运行模型，得到在以上工况下，钢轨、轨道板、自密实混凝土层、底座和桥梁在同一节点处时程曲线垂向加速度最大值的仿真结果对比如下表2。
表2各结构层在不同列车时速下同一节点时程加速度最大值的仿真结果对比
通过模拟列车在时速200km/h、250km/h、300km/h、350km/h的速度下运行模型，得到在以上工况下，钢轨、轨道板、自密实混凝土层、底座和桥梁在同一节点处时程曲线垂向位移最大值的仿真结果对比如下表3。
表3各结构层在不同列车时速下同一节点时程位移最大值的仿真结果对比
通过模拟列车在时速200km/h、250km/h、300km/h、350km/h的速度下运行模型，得到在以上工况下，钢轨、轨道板、自密实混凝土层、底座和桥梁在同一节点处时程曲线拉应力最大值的仿真结果对比如下表4。
表4各结构层在不同列车时速下同一节点时程拉应力最大值的仿真结果对比
结合表1、表2和表3中数据可以看出，轮轨横垂向力、各结构层垂向位移、垂向加速度随着列车速度的增大而不断增大；从表4中可以看出，轨道板的垂向拉应力随着列车速度的增大而逐渐减小，自密实混凝土垂向拉应力、底座纵向拉应力和桥梁纵向拉应力随着列车速度增加而逐渐增加。
实施例2
同crtsiii型板式无砟轨道结构中ca砂浆结构层一样，填充层能够起到良好的变形协调作用，不同客运专线所选用的材料均存在差异，本实施例通过调整模型中无砟轨道结构层中自密实混凝土的材料参数，以此来得到自密实混凝土弹性模量的变化对无砟轨道结构的动力响应值的影响规律。因此，通过将自密实混凝土弹性模量设定为25.5gpa、30gpa、32.5gpa、34.5gpa、36gpa共五种工况，列车速度取350km/h，运行模型，得到无砟轨道各结构层车辆安全性指标动力响应最大值的仿真结果对比如下表5。
表5不同自密实混凝土层厚度下车辆安全性指标动力响应最大值
通过将自密实混凝土弹性模量设定为25.5gpa、30gpa、32.5gpa、34.5gpa、36gpa共五种工况，列车速度取350km/h，运行模型，得到无砟轨道各结构层垂向位移动力响应最大值的仿真结果对比如下表6。
表6不同自密实混凝土层厚度下无砟轨道各结构层垂向位移动力响应最大值
通过将自密实混凝土弹性模量设定为25.5gpa、30gpa、32.5gpa、34.5gpa、36gpa共五种工况，列车速度取350km/h，运行模型，得到无砟轨道各结构层应力指标动力响应最大值的仿真结果对比如下表7。
表7不同自密实混凝土层厚度下无砟轨道各结构层应力动力响应最大值
从表5中可以看出，对于车辆安全性指标来说，自密实混凝土弹性模量的增加会引起车体垂向加速度的增加，而对轮轨垂向力、脱轨系数以及轮轴减载率影响不大。从表6中可以看出，对于垂向位移指标，钢轨、轨道板、自密实混凝土、底座和桥梁的垂向位移均随着自密实混凝土弹性模量的增加而减小，这是由于填充层的刚度变化，故无砟轨道结构层的垂向位移也相应减小。从表7可以看出，自密实混凝土弹性模量的增加会引起轨道板、自密实混凝土和底座的横向拉应力增加，拉应力过大会使得混凝土结构开裂破坏，故在实际环境中，应使用合理强度的混凝土标号。
应当知道，以上所述实施例仅为对本发明的一种说明起着示意性的作用，而并非限制性的作用。对于本领域内的工程技术人员在不脱离本发明的原理下，可以做若干的改进和替换，而这些修改和替换并不使相应技术方案的本质脱离本发明的保护范围。