某炮发射过程的动力学分析.docxVIP

某炮发射过程的动力学分析 采用元方法的试验过程 火炮是一个非常复杂的非线性动力学系统。 基于动力学理论, 应用非线性有限元方法, 将全炮各构件均处理为柔体, 模拟了全炮结构 (含大位移滑动接触) 后坐复进过程中炮口位移和速度的动态响应并获得了发射过程中关键部位的动态刚强度。计算结果可为该舰炮的设计和优化提供一定的理论依据。 1 线性结构分析问题的提出 结构的非线性问题指结构的刚度和阻尼随其变形而改变的分析问题。实际上所有的物理结构均为非线性的, 而线性分析只是一种方便的近似。在结构力学模拟中有三种非线性的来源:材料非线性、边界非线性、几何非线性。 1.1 动力学平衡方程 非线性动力学 其中:[M]、[C]和[K]分别代表质量矩阵、阻尼矩阵和刚度矩阵;[I]=[C]{u}+K{u}为粘性效应项, 考虑阻尼、粘塑、粘弹等效应;[P]为外部激励作用;{u}、{u}和{u}为节点位移、速度、加速度向量; 对应t, t+?t时刻, 动力学平衡方程可表示为如下形式: 为了便于求解, 将式 (1) 改写为增量形式: Newmark法得到t+?t时刻的位移和速度向量可以表示为: 其中:权重因子β∈ (0, 1/2) , γ∈ (0, 1]。 对于Newmark时间积分算法而言, 当取γ=1/2, β=1/4时, Newmark方法即为梯形积分算法 (或称为平均加速度法) , 针对梯形算法无法考虑阻尼影响, Hilber 其中权重因子β、γ如下式: 此时任意给定一个权重因子α∈[-1/3, 0], 就可得到一组对应的权重因子β、γ, 然后通过式 (5) 和 (6) 可求出t+?t的位移和速度向量。 1.2 lagrange乘子法在非零情况下的极值问题描述 接触碰撞问题属于最困难的非线性问题之一, 动力学有限元离散方程包括一般运动方程和在接触界面上要满足的面力条件和不可侵彻性条件。目前ABAQUS/Standard处理接触界面约束的方法是Lagrange乘子法 Lagrange乘子法是通过Lagrange乘子施加接触体必须满足的非穿透约束条件的带约束极值问题的描述方法。在接触问题的离散化中, 接触界面上乘子必须是近似的, 乘子必须满足法向面力是压力的约束, 这种方法是把约束条件加在一个系统中最完美的数学描述。接触问题的离散运动方程和不可侵彻条件为 式中:λ是Lagrange乘子场, ν是接触物体的速度场, A和B表示两个相互接触的物体, Γ Lagrange乘子法没有用户设定的参数, 并且当节点相邻时, 接触约束几乎可以精确的得到满足。该方法增加了未知量数目, 并使系统矩阵主对角线元素为零, 使得在数值方案的贯彻中需要处理非正定系统, 数学上将发生困难。在对于高速碰撞, Lagrange乘子法常常导致非常不平顺的结果, 因此Lagrange乘子法更适于静态和低速的接触问题。在本文中, 火炮身管与摇架导轨间的大位移滑动接触属于低速问题, 在接触计算时, 可以预先知道接触发生的确切部位, 以便通过Lagrange乘子法施加界面接触约束。 2 海洋结构的动态建模 2.1 积分单元划分 为了建立结构动力学模型, 首先用Unigraphic NX4软件建立某舰炮三维实体模型。将火炮质量、质心、转动惯量及结构强度等问题影响不大的零部件去除, 将火炮简化为四大部分:身管、摇架、托架、回转支撑。利用ABAQUS软件分别采用非协调实体积分单元 (C3D8I) 和减缩积分壳单元 (S3R) 对身管、摇架、托架、回转支撑进行了网格划分, 共计约五万个单元, 如图1所示。图中沿身管轴向为z轴, 依右手定则沿身管的横向和垂向分别为x和y轴。 2.2 经济实用的原则 火炮各部件间的连接是非常复杂的, 具有高度的非线性, 需要按照经济实用的原则来处理其间的非线性连接。从火炮整体出发, 会发现能否准确的模拟后坐部分与摇架之间的非线性连接特性及导轨间的大滑动接触是全炮结构动力学仿真的关键, 其它连接采用不同形式的简化。 2.2.1 摇架的后坐过程 驻退机可以在特定炮膛合力和后坐部分质量下, 在炮尾与摇架的驻退机力连接点上加一对共线且反向的力作为载荷, 力的大小等于设计计算或实验测得的总后坐阻力。 该炮摇架为组合型摇架, 即筒形和槽形的混合结构, 接触定义比较复杂。发射过程中, 身管与摇架相对滑动, 在身管与摇架和炮尾与摇架之间都定义了软接触。由于身管与摇架之间的摩擦阻力很难精确给定, 且摩擦阻力相对于其它载荷来说非常小, 在分析中可以忽略不计, 因此将身管与摇架之间的切向行为定义为无摩擦阻力, 以简化问题的求解。法向接触属性结合模型的实际情况定义为指数接触模式。身管与筒型摇架的接触定义如图2和图3所示。 2.2.2 耦合自由度原则 摇架与托架耳轴之间的连接用自由度耦合来定义,