1.一种变轴力双向支座摩阻力计算方法，其特征在于，包括以下步骤：获取结构外加荷载 F 、刚度矩阵 K 和质量矩阵 M 、以及支座参数，所述支座参数包括支座轴向弹性刚度 K n 、支座摩擦系数 μ ；根据所述结构外加荷载 F 、刚度矩阵 K 和质量矩阵 M 建立 t 时刻结构变形增量迭代公式；根据所述支座参数建立变轴力双向支座内力 和切线刚度 公式；将所述变轴力双向支座内力 和切线刚度 公式代入所述 t 时刻结构变形增量迭代公式，计算所述变轴力双向支座摩阻力；所述 t 时刻结构变形增量迭代公式为： ；式中， ；式中， F 为结构外加荷载， 、 和 分别为位移、速度和加速度增量， 为支座内力函数对速度 的导数， 为支座内力函数对位移 的导数， 为荷载残量， 为t时刻外加荷载， 为当前迭代步支座内力， 为当前迭代步位移， 为当前迭代步速度， 为当前迭代步加速度；所述变轴力双向支座内力 公式为： ；式中， ； ；式中， 为支座轴向内力， 为支座水平 x 方向摩阻力， 为支座水平 y 方向摩阻力， 为当前迭代步支座轴向位移， k 为支座滑动速度系数， 、 为前一时刻支座两个方向摩阻系数， 、 为当前迭代步支座两个方向摩阻系数， 为动力分析时间间隔， 为当前时刻支座水平 x 方向速度， 为前一时刻支座水平 x 方向速度， 为当前迭代步支座水平 y 方向速度， 为前一时刻支座水平 y 方向速度；所述变轴力双向支座切线刚度 公式为： ；式中，当 时， ；当 时， ；式中， 为支座轴向切线刚度， 为支座水平 x 方向切线刚度， 为支座 x、y 方向耦合的切线刚度， 为支座水平 y 方向切线刚度。
2.如权利要求1所述的变轴力双向支座摩阻力计算方法，其特征在于，所述结构外加荷载 F 、刚度矩阵 K 和质量矩阵 M 、以及支座参数通过有限元模型获取。
3.如权利要求1所述的变轴力双向支座摩阻力计算方法，其特征在于，所述内力 包括支座轴向压力、支座水平方向的摩阻力。
4.如权利要求1所述的变轴力双向支座摩阻力计算方法，其特征在于，所述支座滑动速度系数 k≥ 10000 。
5.如权利要求1所述的变轴力双向支座摩阻力计算方法，其特征在于，所述变轴力双向支座内力 和切线刚度 公式以数组和矩阵形式代入所述 t 时刻结构变形增量迭代公式。
6.如权利要求5所述的变轴力双向支座摩阻力计算方法，其特征在于，以数组和矩阵形式代入时， 。
7.一种采用如权利要求1所述的变轴力双向支座摩阻力计算方法的变轴力双向支座摩阻力计算装置，其特征在于，其包括：参数获取模块，其用于获取结构外加荷载 F 、刚度矩阵 K 和质量矩阵 M 、以及支座参数，所述支座参数包括支座轴向弹性刚度 K n 、支座摩擦系数 μ ；第一公式建立模块，其用于根据所述结构外加荷载 F 、刚度矩阵 K 和质量矩阵 M 建立 t 时刻结构变形增量迭代公式；第二公式建立模块，其用于根据所述支座参数建立变轴力双向支座内力 和切线刚度 公式；计算模块，其用于将所述变轴力双向支座内力 和切线刚度 公式代入所述 t 时刻结构变形增量迭代公式，计算所述变轴力双向支座摩阻力。
