牵引炮全发射过程非线性有限元建模与分析.docxVIP

牵引炮全发射过程非线性有限元建模与分析 有限源法用于对火炮结构动力学的模拟研究，采用隐式时间积分法和显式时间积分法。由于含大位移滑动接触等非线性因素,目前国内外在这方面的研究多采用显式方法 1 非线性结构分析的基本理论 1.1 节点内力向量的变化 在n+1时间步,结构动态问题的有限元离散动量方程为: 式中:d 动力学问题的隐式时间积分可以通过引入一个非常有效的参数α对方程(1)进行如下变换: 该式与式(1)比较,主要变化在于驱动节点力的位移计算: 对于一个性线系统,节点内力向量的定义为f 应用Newmarkβ法,假设更新的位移和速度为 式中: 式中: 其中 隐式时间积分程序通过控制参数α来控制由于数值方法引进的阻尼:当α=0时积分器没有阻尼,当α=-1/3时积分器有较大的阻尼。应用参数α,可以在没有过多降低精度的前提下,有效地改善高频的数值耗散。 更新的加速度可以通过求解方程(4)得到 将方程(6)代入方程(2)得 方程(7)是节点位移d 1.2 lagrange乘子法 接触问题属于最困难的非线性问题之一,结构动力学有限元离散方程包括一般运动方程和在接触界面上要满足的面力条件和不可侵彻性条件。目前ABAQUS/Standrard处理接触界面约束的方法是Lagrange乘子法。 Lagrange乘子法是通过Lagrange乘子施加接触体必须满足的非穿透约束条件的带约束极值问题的描述方法。在接触问题的离散化中,在接触界面上乘子必须是近似的,乘子必须满足法向面力是压力的约束,这种方法是把约束条件加在一个系统中最完美的数学描述。接触问题的离散运动方程和不可侵彻条件为 式中:A是Lagrange乘子场,v是接触物体的速度场,G=∫ 如果I在A上如果I在B上,A和B表示两个 相互接触的物体,Γ Lagrange乘子法没有用户设定的参数,并且当节点相邻时,接触约束(不可侵彻性条件)几乎可以精确地得到满足。该方法增加了未知量的数目,并使系统矩阵主对角线元素为零,使得在数值方案的贯彻中需要处理非正定系统,数学上将发生困难。在对于高速碰撞,Lagrange乘子法常常导致非常不平顺的结果,因此Lagrange乘子法更适合于静态和低速的接触问题。在本文中,火炮身管与摇架导轨间的大位移滑动接触属于低速问题,在接触计算时,可以预先知道接触发生的确切部位,以便通过Lagrange乘子法施加界面接触约束 2 火炮结构动态有限预算的提出 2.1 有限元模型建立 为了建立有限元模型,首先应用I-DEAS软件建立某牵引火炮的实体模型。在建立火炮实体模型时,抑制了对火炮质量、质心位置、转动惯量及结构强度等问题影响不大的零部件,如瞄准机构、防危板等,该火炮简化为5部分:火炮身管、摇架、上架、下架、大架。 进行动力学分析需要耗费大量的计算资源,为合理减少有限元模型的规模,应用壳单元模拟火炮架体中大量的焊接钢板,用梁单元来模拟火炮身管。由于本文重点关注大架在火炮发射过程中的应变、应力,所以应用减缩积分二阶壳单元(S8R)对其进行网格划分,其他部件应用一阶壳单元(S4R)。本文应用实体单元、壳单元和梁单元相结合的方法建立了该炮的有限元模型,该模型由317个实体单元、13 563个壳单元、12个梁单元,共32 557个节点组成,如图1所示。图中火炮身管处于xy平面内,垂直向上的方向为y轴方向,水平正前方向为x方向,身管与x轴成70°射角,z轴按右手法则确定。 2.2 反后坐装置非线性连接模拟的可行性 火炮各部件间的连接非常复杂,往往具有高度非线性特性,需要按照经济、有效、实用的原则来处理其间的非线性连接。从火炮整体来看,能否准确地模拟反后坐装置非线性连接特性及炮身与摇架导轨之间的大滑移接触是全炮结构动力学建模和数值模拟的关键,需要重点考虑。对其他连接,则采用不同形式的简化模型进行模拟。 2.2.1 刚性接触模拟 该火炮发射过程中,身管上的滑块沿着摇架导轨进行后坐-复进运动,并产生不可忽视的接触力。能否准确模拟这一过程是全炮结构动力学分析的一个关键问题。由于滑块尺寸较小,本文将滑块接触面进行刚性处理,在刚性面与摇架导轨接触面之间定义大滑移接触。 火炮复进机力是后坐位移的非线性函数,用一个非线性弹簧单元模拟。弹簧特性从火炮设计说明书中复进机力与后坐位移的关系确定,如图2所示。 驻退机节流孔的平方阻尼系数是其面积的函数,而该面积随后坐位移变化。可以用一个具有复杂特性的阻尼器单元去模拟驻退机特性。也可以在特定炮膛合力和后坐部分质量下,在炮尾与摇架的驻退机力连接点上加一对共线且反向的力作为载荷,力的大小等于设计计算或实验测得的总后坐阻力中扣除复进机力和后坐部分重力分量的成分。本文采用后一种方法。 2.2.2 耦合y方向自由度约束 摇架与上架耳轴之间的连接用自由度耦合来定义,只释放绕耳