一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法

1.一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，包括，其特征在于：包括以下步骤：步骤A，铁路简支梁桥墩任意截面输入制定截面输入规则，将桥墩沿着高度方向分段输入，每个节段均输入顶面和底面的截面信息；步骤B，铁路简支梁桥墩节点离散及标序对铁路简支梁桥墩，从墩顶到墩底截面，按一定间距划分节点，并标注节点号，两节点之间构成一个梁单元；步骤C，计算单位力作用下桩基水平位移α和转角位移β在桩基承台底部中心作用单位水平力，得到桩基承台的水平位移α和转角位移β；步骤D，计算单位力作用下墩身弹性变形在桥墩节点上作用一个单位力，计算出作用节点处墩身的弹性水平位移和转角位移；步骤E，形成桩基耦合柔度矩阵F综合考虑桩基变形和墩身弹性变形的影响，分别计算出墩身受力节点的水平位移及其它节点的水平位移，并将两类水平位移组建成桩基耦合柔度矩阵F；步骤F，形成质量矩阵M按照桥墩的节点划分，将节点上下两段质量之和的一半计入该节点的质量，构成质量矩阵的对角阵，质量矩阵中其它矩阵元素为零，对于需要考虑列车、梁体、支座质量的工况，分别将列车、梁体、支座的质量计入质量矩阵的对角阵中；步骤G，求解各阶振型及特征值利用上述桩基耦合柔度矩阵和质量矩阵列出特征方程，利用迭代算法计算出各阶振型和特征值；步骤H，求解地震力利用步骤G得到的各阶振型和频率所对应的桥墩各节点的地震力，将各阶振型的地震力进行叠加得到最终地震力。 2.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤A中截面输入规则适用于任意一种形状的桥墩截面，截面输入规则具体如下：将截面的边按照逆时针方向输入，每条边数据中间用分号隔开，每条边包含五项信息，每项信息之间用逗号隔开，五项信息如下：a.轮廓编号，其中，外轮廓大于零，内轮廓小于零；b.起始端点X坐标；c.起始端点Z坐标；d.半径R，其中，数据等于零代表直线段，数据大于零表示圆弧的半径值；e.圆弧标示，其中数据等于零代表整圆，数据大于零表示劣弧，数据小于零代表优弧。 3.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤B中对桥墩沿着竖直方向进行节点划分，从墩顶到墩底进行标号，墩身可划分为NP段，基础可划分为NC段，考虑列车、简支梁梁体及支座的影响，分顺桥向、横桥向有车、横桥向无车三种工况进行节点标号，各工况的节点数为：横桥向有车工况：NP+NC+1+3横桥向无车工况：NP+NC+1+2顺桥向工况：NP+NC+1+1以上，NP、NC分别为墩身和基础的节段划分数量。 4.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤C中计算单位力作用下桩基水平位移α和转角位移β，按照以下公式进行计算：以上，α、β分别为承台板底面的水平位移和绕坐标原点的转角，γ αα 、γ ββ 、γ αβ 、γ βα 参数为桩基计算参数，H、M分别为承台板底面原点所承受的水平力和弯矩，此处H取单位力，M取零。 5.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤D中计算单位力作用下墩身弹性变形，按照以下公式进行计算：其中，以上，y 0 (n)、 分别为单位力作用下墩身弹性变形受荷节点处位移和转角，E i 为第i个单元的弹性模量，l i 、l j 为第i、j个单元的单元长度。 6.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤E中，形成桩基耦合柔度矩阵F的计算过程如下：首先，单位力作用下考虑桩基耦合效应的受荷节点处位移计算公式为：CF(i,i)＝α+β·L i +y 0 (n)以上，CF(i,i)为单位力作用下墩身和桩基础共同影响下的受荷节点处位移，L i 为第i个节点到承台板底面的距离；然后，单位力作用下考虑桩基耦合效应的各节点的转角位移计算公式为：以上，fai_1为单位力作用下各节点的转角，再后，单位力作用下考虑桩基耦合效应的受荷节点以外各节点的位移，表达式为：CF(j,i)＝CF(i,i)+fai_1·L j以上，CF(j,i)为单位力作用下墩身和桩基础共同影响下的受荷节点以外各节点的位移，L j 为第j个节点到承台板底面的距离，最后，将CF(i,i)和CF(j,i)组合而形成柔度矩阵F。 7.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤F中形成质量矩阵M的具体过程如下：按照步骤B对墩台基础进行节点划分后，两节点之间的任一个单元均可计算出质量，可记作m 0i ，对于节点处的集中质量m i ，取节点上下两段质量之半进行计算，公式为：m i ＝(m 0i +m 0i+1 )/2对于桥墩墩顶的节点，暂且记作0号节点，其公式为：m 0 ＝m 00 /2对于基础底层的节点，暂且记作NP+NC号节点，其公式为：m NP+NC ＝m 0Np+NC /2对于墩顶以上的节点，其分别代表支座+垫石质量、梁体质量、列车活载质量，分别将其竖向质量计入相应的质量矩阵位置中，以上，m i 为墩身中间节点的质量，m 0 为墩顶节点的质量，m NP+NC 为基础底面节点的质量,m 0i 、m 0i+1 分别为第i、i+1节段的质量，m 00 、m 0Np+NC 分别为墩顶位置和基础底位置节段的质量。 8.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤G求解各阶振型及特征值，具体计算过程如下：首先，利用柔度矩阵和质量矩阵，列出通用特征方程：FMX＝λX然后，利用迭代法求解各阶振型频率，过程如下：首先假设向量 代入公式通用特征方程中，得到当按上式迭代求出第一振型 和第一特征值λ 1 后，第二振型向量 可下式计算按上式采用迭代方式收敛到第二振型和第二特征值的结果，同理，求出第n阶振型和特征值的求解迭代公式以上，F为柔度矩阵，M为质量矩阵，λ 1 、λ 2 ……λ n 为各阶特征值， 为各阶振型，T s1 、T s2 ……T sn-1 为各阶振型的消去矩阵。 9.根据权利要求1所述的一种基于墩身任意截面的铁路简支梁桥墩地震力计算方法，其特征在于：步骤H中求解地震力的具体计算过程如下：首先，根据桥梁所在场地条件因素，查询桥梁的重要性系数C i 、水平地震基本加速度α，然后，对于任一振型j，计算该振型的动力放大系数β j ，当结构自振周期小于2s，且阻尼比ζ＝0.05时，动力放大系数β可如下取值：当T＜0.1 当0.1≤T≤T g β j ＝2.25当T g ＜T≤5T g 当T＞5T g β j ＝0.45当T≥2s时或者阻尼比ζ≠0.05时，另行研究，以上，T为第j阶的自振周期,T g 为地震动反应谱特征周期；再后，对于任一振型j，按以下公式，计算该振型的参与系数，再后，对于任一振型j，按以下公式，计算任一节点i处由自身节点质量引起的水平地震力:F ijE ＝C i ·α·β j ·γ j ·x ij ·m i再后，计算任一截面位置地震力之和及地震弯矩之和，任一截面位置以节点号表示，以节点i位置为例：节点i处的水平地震力，需要对该节点即i点及以上各点的水平地震力求和，记做Pi：节点i处的水平地震力产生的弯矩：最后，各振型地震作用效应的加权平均地震作用效应包括地震力、地震弯矩和位移参数，对于多阶振型，应按照如下公式进行加权平均：以上，F ijE 为j振型i点的水平地震力，C i 为桥梁的重要性系数，α为水平地震基本加速度，β j 为j振型动力放大系数，γ j 为j振型参与系数，x ij 为j振型在第i段桥墩质心处的振型坐标，S i 为地震作用效应，可代表位移、水平地震力、地震弯矩参数。