一种三滑块单轨柔性火箭橇的垂向响应计算方法
【专利摘要】本发明公开了一种三滑块单轨柔性火箭橇的垂向响应的计算方法。该方法首先对火箭橇进行受力分析，将箭体离散为有限元梁模型；进而考虑轨道不平顺引起火箭橇与轨道结构振动，采用Hertz非线性接触模型，计算得到箭体各时刻受到的碰撞接触力；最终通过中心差分的数值方法，计算得到给定时间内三滑块单轨火箭橇系统全程的垂向动力学响应。本发明在计算过程中考虑了火箭橇橇体的柔性变形特性及火箭橇系统由于喷出燃料导致的质量折减效应，并且考虑了轨道不平稳引起的非线性接触问题，因此可以得到较为可靠的动力学响应结果。
【专利说明】
-种Ξ滑块单轨柔性火箭橋的垂向响应计算方法
技术领域
[0001] 本发明设及一种Ξ滑块单轨柔性火箭橋的垂向响应计算方法，设及采用heパz接 触模型求解碰撞接触力W及采用化wmark算法进行动力学数值计算，从而实现火箭橋在不 平顺轨道上动态接触力求解和垂向动力学响应的高效和快速预测的计算方法。
【背景技术】
[0002] 火箭橋是采用火箭发动机作为动力，沿着专口建造的轨道运行的一种试验工具。 火箭橋试验受试验件体积、外形和质量的影响较小，大至数百千克的全尺寸部件、小至数十 克的缩比模型都可W在火箭橋上进行试验。在轨道足够长的情况下，采用不同推力的发动 机或不同数量发动机组合作为橋车的动力源时，橋车的速度可W在几十米每秒至上千米每 秒的范围内进行调节，而且火箭橋试验能够很方便的观察试验结果和重复收集试验数据。 因此，在先进常规武器系统研制（如导弹碰撞和侵彻试验、子母弹抛撒试验、惯性制导系统 试验等）、飞机弹射座椅的开发、空气动力学试验等方面有极为广泛的应用。
[0003] 第二次世界大战后，美国、英国、法国和前苏联等国家为了推动武器、航空航天的 发展，相继建造了多种不同类型的火箭橋试验场，开展了包括惯性测量装置在内的一些核 屯、部件的火箭橋综合验证试验工作。火箭橋试验的显著特点是可无损回收被试惯性测量装 置，供进一步测量、检查及继续进行试验。高精度的惯性测量装置造价高，通过火箭橋试验 可重复进行多类多次的测试试验，包括环境适应性试验和精度试验，增加试验样本量，确保 飞行试验一次成功，减少试验次数，降低试验成本，加快研制周期。美国在应用火箭橋对惯 性测量装置进行试验研究方面发展最快，应用最广。从上世纪四十年代末到屯十年代初，美 国相继建造了多达25条的不同类型的火箭橋试验轨道，先后开展了大量的火箭橋试验项目 研究。
[0004] 火箭橋约束在轨道内高速运动的过程中，受到火箭推力、碰撞接触力、空气阻力和 制动力等作用，火箭橋与轨道的动态性能直接影响到车载试验件的试验环境，将直接决定 动态试验结果的准确性和可靠性。由于轨道长度与轨道不平顺的数量级相差巨大，并且橋 轨在轨道不平顺点接触所产生的附加接触变形量的精确解往往是通过多次迭代计算所求 得，导致有限元建模和计算都十分困难。
[0005] 开展高速火箭橋与轨道系统动态性能的研究，建立橋-轨动力学数学模型，有助于 提高对火箭橋、轨道动力学特性的认识，掲示系统物理本质，为研制满足要求、性能良好的 试验是被提供设计依据，最终使火箭橋及试验轨道实际更为优化可靠。对于高速火箭橋-轨 道禪合系统动力学仿真分析研究，具有重要的科学使用意义和工程应用价值。
[0006] 目前尚未有相关技术报导。
【发明内容】
[0007] 本发明解决的技术问题是:本发明克服现有技术建模困难、计算量巨大的不足，提 供一种Ξ滑块单轨柔性火箭橋的垂向响应计算方法，求解速度快，工程适用性强，在计算过 程中考虑了火箭橋橋体的柔性变形特性及火箭橋系统由于喷出燃料导致的质量折减效应， 并且考虑了轨道不平稳引起的非线性接触问题，因此可W得到较为可靠的动力学响应结 果。
[0008] 本发明采用的技术方案，一种Ξ滑块单轨柔性火箭橋的垂向响应计算方法，基于 Hertz接触力模型，实现步骤如下：
[0009] 第一步:根据Ξ滑块单轨柔性火箭橋的运动形式，将运动等效为航向的刚体平动 和垂向的振动;Ξ滑块单轨柔性火箭橋受到航向的推力和垂向接触力，根据模型受力形式， 自行构造一种5自由度的梁单元，即在欧拉梁单元的基础上附加一个航向刚体自由度，将Ξ 滑块单轨柔性火箭橋模型离散为3个5自由度的梁单元模型，得出Ξ滑块单轨柔性火箭橋模 型的动力学矩阵M、C和K，建立运动学方程为：
[0010] 編 + Cii .4· Λ'μ 二/、丫/;
[0011] M、C和Κ分别为整体质量阵、刚度阵和阻尼阵，F(t)为模型所受外力；
[0012] 第二步:根据第一步求得Ξ滑块单轨柔性火箭橋模型，然后确定模型的边界条件， 所述边界条件包括载荷与约束条件;定义迭代步长A T，已知初始时刻即Τ = 0时刻的航向位 移sx、航向速度vx及航向加速度ax和外载荷曲线，所述外载荷曲线为Ξ滑块单轨柔性火箭橋 推力、阻力和质量折减曲线，为初始已知条件，将时刻T代入外载荷曲线，时刻T由步长ΔΤ和 迭代步数η的乘积决定，得到对应时刻的外载荷与质量，根据Newmark算法，得到下一时刻T+ A T的航向位移、航向速度及航向加速度，然后根据迭代步数η求得任意时刻的航向位移、航 向速度及航向加速度；
[0013] 第Ξ步:根据第二步得到Τ时刻航向位移Sx确定Ξ滑块单轨柔性火箭橋在轨道的 距离起点即τ=0时刻所在位置，从而确定轨道不平顺信息即实测得到的轨道各测点高度与 基准面高度的差，结合对应时刻Τ的Ξ滑块单轨柔性火箭橋的垂向位移、速度和加速度，计 算接触变形与接触相对方向速度；
[0014] 第四步:滑靴相当于车轮，为火箭橋运动时与火箭橋轨道相接触的部分。确定对应 Τ的滑靴与轨道的接触状态，确定碰撞接触力大小，当滑靴在轨道上运动时，轨道与滑靴之 间存在着Ξ种状态:①滑靴与轨道上表面接触、②滑靴与轨道下表面接触、③滑靴与轨道不 接触;接触状态不同，碰撞接触力的计算公式有所不同，首先根据轨道与滑靴的初始间隙W 及当前时间步滑靴位置，确定出当前时刻的接触状态，再根据heパz接触力模型计算碰撞接 触力；
[0015] 第五步:根据第四步计算的碰撞接触力W及对应时刻ΤΞ滑块单轨柔性火箭橋的 垂向位移、速度和加速度作为输入，施加到Ξ滑块单轨柔性火箭橋模型上，采用化wmark方 法求解出下一时刻的动力学响应，即为T+ Δ T时刻的垂向位移、速度和加速度；
[0016] 第六步:将第五步得到Τ+ΔΤ动力学响应作为输入，回到第四步，确定对应时刻T+ AT的滑靴与轨道的接触状态，确定碰撞接触力大小。应用第五步中所述，将第五步得到T+ A T动力学响应施加到S滑块单轨柔性火箭橋模型上，采用化wmark方法，从而求得T+2 A T 的动力学响应，循环至时间步长达到所要求时间步为止，最终得到Ξ滑块单轨柔性火箭橋 在给定时间范围内的动力学响应，即垂向位移、垂向速度W及垂向加速度。
[0017] 本发明与现有技术相比的优点在于:本发明W轨道不平顺数值、轨道与滑靴初始 间隙、火箭橋推力曲线为输入条件，将火箭橋离散为5自由度梁单元，采用Newmark数值算 法，应用接触力模型，根据当前时间步的垂向位移和速率更新对应时间步的碰撞接触 力，最终求得各有限元节点的垂向动力学响应，具有求解速度快，工程适用性强，在计算过 程中考虑了火箭橋橋体的柔性变形特性及火箭橋系统由于喷出燃料导致的质量折减效应， 并且考虑了轨道不平稳引起的非线性接触问题，因此可W得到较为可靠的动力学响应结 果。
【附图说明】
[001引图1为本发明实现总流程图；
[0019 ]图2为本发明中的5自由度梁单元模型图；
[0020]图3为本发明中滑靴与轨道接触状态种类；
[0021 ]图4为本发明中Ξ滑块单轨柔性火箭模型；
[0022] 图5为本发明中Ξ滑块单轨柔性火箭橋简化模型；
[0023] 图6为本发明中推力曲线；
[0024] 图7为本发明中实测轨道不平顺值；
[0025] 图8为本发明中节点1垂向加速度响应(正弦不平顺值下）；
[0026] 图9为本发明中节点3垂向加速度响应(正弦不平顺值下）；
[0027] 图10为本发明中节点5垂向加速度响应(正弦不平顺值下）；
[0028] 图11为本发明中节点1垂向加速度响应；
[0029] 图12为本发明中节点3垂向加速度响应；
[0030] 图13为本发明中节点5垂向加速度响应。
【具体实施方式】
[0031] 如图1所示，本发明提出了一种基于hedz接触模型的Ξ滑块单轨柔性火箭橋垂向 响应计算方法。包括W下步骤：
[0032] (1)根据Ξ滑块单轨柔性火箭橋的运动形式，将运动等效为航向的刚体平动和垂 向的振动;Ξ滑块单轨柔性火箭橋受到航向的推力和垂向接触力，根据模型受力形式，自行 构造一种5自由度的梁单元，即在欧拉梁单元的基础上附加一个航向刚体自由度，将Ξ滑块 单轨柔性火箭橋模型离散为3个5自由度的梁单元模型，如图2所示，得出Ξ滑块单轨柔性火 箭橋模型的动力学矩阵M、C和K，建立运动学方程为：
[0033]
[0034] M、C和K分别为整体质量阵、刚度阵和阻尼阵，F(t)为模型所受外力；
[0035] 对于航向刚体运动，Ξ滑块单轨柔性火箭橋航向运动简化为杆的刚体平动，故其 运动方程可由其质屯、运动方程表达，由牛顿第二定律得到：
[0036]
[0037] 其中，m为箭体质量，Fi和F2分别为作用在5自由度梁单元两个节点处X方向上的合 力。
[0038] 对于垂向振动，梁垂向振动有限元方程可写为W下形式：
[0039]
[0040] 其中化、Ki、Ci分别为质量阵、刚度阵和阻尼阵，u=[yi目1 y2目2]τ为节点纵向位移。 F(t)外力向量。
[0041] 综合考虑梁的垂向振动和航向刚体运动，Ξ滑块单轨柔性火箭橋的运动方程可写 为：
[0042]
[0043] 其中M、K、C分别为整体质量阵、刚度阵和阻尼阵。
[0044] (2)根据第一步求得Ξ滑块单轨柔性火箭橋模型，然后确定模型的边界条件，所述 边界条件包括载荷与约束条件；定义迭代步长A T，已知初始时刻即Τ = 0时刻的航向位移 sx、航向速度vx及航向加速度ax和外载荷曲线，所述外载荷曲线为Ξ滑块单轨柔性火箭橋推 力、阻力和质量折减曲线，为初始已知条件，将时刻T代入外载荷曲线，时刻T由步长ΔΤ和迭 代步数η的乘积决定，得到对应时刻的外载荷与质量，根据化wmark算法，得到下一时刻T+ Δ T的航向位移、航向速度及航向加速度，然后根据迭代步数η求得任意时刻的航向位移、航 向速度及航向加速度；
[0045] (3)根据第二步得到Τ时刻航向位移确定Ξ滑块单轨柔性火箭橋在轨道的距离起 点即τ=ο时刻所在位置，从而确定轨道不平顺信息即实测得到的轨道各测点高度与基准面 高度的差，其中不平顺值是根据轨道实际情况每隔一定距离测定得到的一系列数据点列， 要想获得任意位置的不平顺值需要对原数据进行插值(如线性插值），轨道任意点的不平顺 值与不平顺斜率计算方法如下所示：
[0048] 其中，S为滑靴当前位置的轨道不平顺值，s/为滑靴当前位置轨道不平顺斜率，Sn+l 为滑靴位置后一个轨道监测点的监测值，Sn为滑靴位置前一个轨道监测点的监测值，h为滑 块位置距离前一个观测点的长度，L为监测点间隔。
[0049] 结合对应时刻Ξ滑块单轨柔性火箭橋的垂向位移、速度和加速度，计算接触变形 与接触相对方向速度。
[0050] (4)确定对应时刻T的滑靴与轨道的接触状态，确定碰撞接触力大小。当滑靴在轨 道上运动时，轨道与滑靴之间存在着Ξ种状态（如图3所示）：①滑靴与轨道上表面接触、② 滑靴与轨道下表面接触、③滑靴与轨道不接触。接触状态不同，碰撞接触力的计算公式有所 不同。首先根据轨道与滑靴的初始间隙W及当前时间步滑靴位置，确定出当前时刻的接触 状态，再根据hedz接触力模型计算碰撞接触力。分为如下Ξ种情况：
[0051] A.滑靴和轨道不接触。此时得到当前碰撞接触力为0;
[0052] B.滑靴与轨道下缘接触。此时，接触变形δ由当前时间步滑靴垂向位移、不平顺值 和初始间隙作差得到。接触变形速率省取为当前时间步的垂向速率。则接触力有下列公式求 得：
[0化3]
(3)
[0054] C.滑靴与轨道上缘接触。此时，接触变形δ由当前时间步滑块垂向位移、不平顺值 和初始间隙作差得到。接触变形速率取为当前时间步的垂向速率。则接触力有下列公式 求得：
[ο化5]
(4).
[0化6]其中，Η(δ)代表非线性接触力，Κδη代表的是弹性接触部分，辦r凌是冲击阻尼力部 分，梦为相对速度，μ是迟滞阻尼系数，其值与粘性阻尼系数、剪切和体积变形等有关。K和μ 的值可由试验测定或根据数值模拟方法对其进行拟合得到。
[0057] (5)根据第四步计算的碰撞接触力W及对应时刻ΤΞ滑块单轨柔性火箭橋的垂向 位移、速度和加速度作为输入，施加到Ξ滑块单轨柔性火箭橋模型上，采用化wmark方法求 解出下一时刻的动力学响应，即为T+ Δ T时刻的垂向位移、速度和加速度。
[005引（6)将第五步得到Τ+ΔΤ动力学响应作为输入，回到第四步，确定对应时刻Τ+ΔΤ的 滑靴与轨道的接触状态，确定碰撞接触力大小。应用第五步中所述，将第五步得到Τ+ΔΤ动 力学响应施加到Ξ滑块单轨柔性火箭橋模型上，采用化wmark方法，从而求得T+2 Δ T的动力 学响应，循环至时间步长达到所要求时间步为止，最终得到Ξ滑块单轨柔性火箭橋在给定 时间范围内的动力学响应，即垂向位移、垂向速度W及垂向加速度。
[0059] 实施例：
[0060] 本发明实施例采用一个Ξ滑靴的火箭橋模型，如图5所示。航向过载曲线如图6所 示，前3秒为加速段。划分五个5自由度梁单元，试验件划分1个单元，发动机划分4个单元，示 意图如图所示。各段截面参数W及碰撞参数如下表1所示，滑靴与轨道间的间隙取为1.5mm。 [0061 ] 表1基本参数
[0062]
[0063] 利用化dran语言按照图1的流程编写求解程序，轨道不平顺值采用两种情况：正 弦不平顺值和实测获得。求解得到节点垂向加速度响应如下所述：
[0064] (1)对于正弦不平顺值，即轨道不平顺按正弦规律沿轨道坐标方向变化，幅值取 1mm，波长取8m。得到各节点的加速度响应如8、9和10所示。可见Ξ滑块单轨柔性火箭橋各节 点垂向加速度响应随着时间的增加，有逐渐增大的趋势。
[0065] (2)对于实测得到的不平顺值(7)，得到各节点的加速度响应如1、图12和3所示。可 见火箭橋各节点垂向加速度响应随着时间的增加，有逐渐增大的趋势，运一点与正弦不平 顺值下的垂向响应相同，但是变化趋势不如正弦情况下均匀。运与实际测量值的趋势相符。
[0066] 综上所述，本发明提出了一种柔性火箭橋垂向响应计算方法。首先建立火箭橋模 型，将其离散为5自由度梁单元，进而组装成有限元动力学矩阵。再将不平顺值信息、推力曲 线作为输入条件，对每个求解时间步进行推力与碰撞接触力的更新，并求解下一时间步的 动力学响应。最终便得到全程的垂向动力学响应。
[0067] 本发明未详细阐述部分属于本领域技术人员的公知技术。
【主权项】
1. 一种三滑块单轨柔性火箭橇的垂向响应计算方法，其特征在于实现步骤如下： 第一步:根据三滑块单轨柔性火箭橇的运动形式，将运动等效为航向的刚体平动和垂 向的振动;三滑块单轨柔性火箭橇受到航向的推力和垂向接触力，根据模型受力形式，构造 一种5自由度的梁单元，即在欧拉梁单元的基础上附加一个航向刚体自由度，将三滑块单轨 柔性火箭橇模型离散为3个5自由度的梁单元模型，得出三滑块单轨柔性火箭橇模型的动力 学矩阵M、C和K，建立运动学方程为： Mii + Cu+Mu^-F(t) M、C和K分别为整体质量阵、刚度阵和阻尼阵，F(t)为模型所受外力； 第二步:根据第一步求得三滑块单轨柔性火箭橇模型，确定该模型的边界条件，所述边 界条件包括载荷和约束条件;定义迭代步长A T，已知初始时刻即T = 0时刻的航向位移、航 向速度及航向加速度和外载荷曲线，所述外载荷曲线为三滑块单轨柔性火箭橇推力、阻力 和质量折减曲线，为初始已知条件，将时刻T代入外载荷曲线，时刻T由步长△ T和迭代步数η 的乘积决定，得到对应时刻的外载荷与质量，根据Ne wmar k算法，得到下一时刻Τ+ △ Τ的航向 位移、航向速度及航向加速度，然后根据迭代步数η求得任意时刻的航向位移、航向速度及 航向加速度； 第三步:根据第二步得到Τ时刻航向位移，确定三滑块单轨柔性火箭橇在轨道的距离起 点即初始时刻Τ = 0时所在位置，从而确定轨道不平顺信息即实测得到的轨道各测点高度与 基准面高度的差，结合对应时刻Τ的三滑块单轨柔性火箭橇的垂向位移，计算接触变形与接 触相对方向速度； 第四步：确定对应时刻Τ的滑靴与轨道的接触状态，确定碰撞接触力大小;所述滑靴相 当于车轮，为三滑块单轨柔性火箭橇运动时与三滑块单轨柔性火箭橇轨道相接触的部分； 当滑靴在轨道上运动时，轨道与滑靴之间存在着三种状态:①滑靴与轨道上表面接触;②滑 靴与轨道下表面接触;③滑靴与轨道不接触;接触状态不同，碰撞接触力的计算公式有所不 同，首先根据轨道与滑靴的初始间隙以及当前时间步滑靴位置，确定出当前时刻的接触状 态，再根据hertz接触力模型计算碰撞接触力； 第五步:根据第四步计算的碰撞接触力以及对应时刻T三滑块单轨柔性火箭橇的垂向 位移、速度和加速度作为输入，施加到三滑块单轨柔性火箭橇模型上，采用Newmark方法求 解出下一时刻的动力学响应，即为T+ △ T时刻的垂向位移、速度和加速度； 第六步:将第五步得到Τ+Λ T动力学响应作为输入，回到第四步，确定对应时刻Τ+Δ T的 滑靴与轨道的接触状态，确定碰撞接触力大小，应用第五步中所述，将第五步得到Τ+ΔΤ动 力学响应施加到三滑块单轨柔性火箭橇模型上，采用Newmark方法，从而求得Τ+2 Δ T的动力 学响应，循环至时间步长达到所要求时间步为止，最终得到三滑块单轨柔性火箭橇在给定 时间范围内的动力学响应，即垂向位移、垂向速度以及垂向加速度。2. 根据权利要求书1中一种三滑块单轨柔性火箭橇的垂向响应计算方法，其特征在于： 所述第一步中得出三滑块单轨柔性火箭橇模型的动力学矩阵M、C和K过程如下： 将三滑块单轨柔性火箭橇模型离散为5自由度单元，5自由度梁单元的自由度为两端节 点的挠度和转角71，91，72，02以及航向的刚体平动自由度1 11，对于航向刚体平动，三滑块单轨 柔性火箭橇航向运动方程可由其质心运动方程表达，由牛顿第二定律得到： niii, = l·] + l·) 其中，m为箭体质量，Fi和F2分别为作用在5自由度梁单元两个节点处x方向上的合力； 对于垂向振动，5自由度梁单元垂向振动有限元方程写为以下形式： M.y + C.y + K.y^Fd) 其中^&分别为质量阵^度阵和阻尼阵…二卜^^^斤为节点纵向位移^⑴ 外力向量， 综合考虑梁的垂向振动和航向刚体运动，三滑块单轨柔性火箭橇的运动方程为： Mil + Cii + Ku = F(i) 其中M、K、C分别为整体质量阵、刚度阵和阻尼阵。3.根据权利要求书1中一种三滑块单轨柔性火箭橇的垂向响应计算方法，其特征在于： 第四步中，所述hertz接触力模型为： Ι!{δ) = Κδη +μ?ηδ 其中，Η(δ)代表非线性接触力，Κδη代表的是弹性接触部分，是冲击阻尼力部分j 为相对速度，μ是迟滞阻尼系数，μ值与粘性阻尼系数、剪切和体积变形有关，κ和μ的值由试 验测定或根据数值模拟方法进行拟合得到。
【文档编号】G06F17/50GK106066922SQ201610420732
【公开日】2016年11月2日
【申请日】2016年6月13日 公开号201610420732.5, CN 106066922 A, CN 106066922A, CN 201610420732, CN-A-106066922, CN106066922 A, CN106066922A, CN201610420732, CN201610420732.5
【发明人】王晓军, 田靖军, 王磊, 管闯闯, 马雨嘉
【申请人】北京航空航天大学