舰炮抽壳系统柔体动力学分析方法研究.docxVIP

舰炮抽壳系统柔体动力学分析方法研究

张建，唐文献

（江苏科技大学江苏镇江212003）

摘要：抽壳系统的动态特性直接影响到舰炮的可靠性，由于在实际工作工程中存在多重非线性，难以用传

统方法研究抽壳机构的力学特性能。为此，基于现代设计理论与方法，提出一种抽壳系统柔体动力学分析

方法。应用显式积分法求解非线性方程组，三维非协调实体单元离散几何模型，通用接触算法模拟部件之

间的动态接触，Rayleigh阻尼和Mises屈服面定义材料特性，失效连接单元模型弹壳抽壳特性。为了检验此

分析方法的正确性，以某舰炮的抽壳系统为实际案例，使用PROE软件建立实体模型，HyperMesh软件进

行网格划分，电子拉伸试验机测试材料的力学性能，然后采用ABAQUS\Explicit进行数值计算，最后通过

动应力测试验证模型的正确性。结果表明，本文所提出的方法有效可靠，可以为舰炮抽壳系统的研究提供

理论基础。

关键词：舰炮抽壳系统柔体动力学分析方法有限元分析

中图分类号：TJ303文献标识码：A

舰炮抽壳系统的主要作用是抽出药筒（弹进行系统性多柔体分析的研究很少，因为存在如下

壳），其质量状况的好坏直接影响抽壳过程和后续三个难点：(1)由于系统的大位移和大转动引起的

炮弹的进膛。在对某新开发的舰炮物理样机进行几何非线性；(2)在工作过程中，材料会产生塑性

试验时，抽壳系统各部件存在很多问题，如抽筒变形，存在明显的物理非线性；(3)系统各部件之

子时常发生塑性变形甚至断裂，抽壳模板部位磨间存在非线性接触。整个系统的多柔体有限元分析

损比较严重导致无法抽出药筒，从而失去战斗力。模型只有充分考虑这些因素才能有效可靠。

因此，有必要研究抽壳系统的力学特性。本文针对这些关键问题所涉及动力学方程理

试验方法用于研究舰炮的力学特性，是一种论、材料本构、阻尼、接触等问题展开分析，提出

最直接且最准确的方式，但测试周期比较长且耗抽壳系统柔体动力学分析方法。为了验证该方法的

费的经费大。随着计算机技术、先进计算方法以可行性，以某的抽壳系统为工程实例，根据结构尺

及力学理论的发展，计算机辅助工程（Computer寸用PROE软件建立三维实体模型，用HyperMesh

AidedEngineering）应用于预测舰炮系统的力学

软件划分网格，测试了的炮钢材料的力学特性，最

特性成为可能，这种基于虚拟技术的现代设计理后定义载荷和边界条件进行数值计算，并通过动态

论与分析方法在军工行业的广泛应用[1][2][3]，可发

应力试验验证了模型的正确性。结果表明，本文提

现及早结构设计缺陷，减少物理样机的试制次数，出的多柔体动力学分析方法是切实可行的，为舰炮

有效地缩短产品开发周期及降低开发费用。然而，抽壳机构的研制开发提供了一种较为经济实用的分

大多数研究仅局限于静态有限元分析、多刚体动析方法。

力学分析以及单个零件动应力校核上，而对舰炮1柔性体动力学分析理论

1.1柔性体动力学方程抽壳系统是由抽筒子、抽筒子轴、摇臂、摇臂

柔性体的动力学运动方程是建立在广义坐标轴、抽壳模板等组成复杂系统。对抽壳系统进行有

基础上，它不仅要能够反映柔性体的大范围或非限元建模时，传统做法采用结构单元来离散这些部

线性刚体位移，而且还需体现柔性体微小弹性变件，这必然要对结构作较大简化，并为了使不同类

形。运用广义坐标推导的基于拉格朗日方程的控型单元边界处相互协调，需要进行较多的人工干预，

制性动力微分方程的形式如下：劳动量大且很难反映系统的真实情况；对抽壳系统

TT

1?M(1)

`

??????

M?????????????K??f?D?????F

M

????

2????

g

????

采用统一的三维实体单元来划分网格，单元的一致

式中：?、??、???是自由体广义坐标及其对

性较好，网格质量较高，分析工程师的劳动量也大

为减少。

时间的导数,M是柔性体的质量矩阵,M?是柔性

实体单元可分为普通单元和非协调单元，目前，

体质量矩阵一阶导数,?M??是质量矩阵关于广

CAE软件中大都采用等参单元，而改善等参单元性

义坐标偏导数,K是广义刚度矩阵,fg是广义重

能的一个重要步骤是引入非协调模式，此时单元位

移模式就产生了附加自由度，可以使用较少单元得

力,D是模态阻尼矩阵,?是代数约束方程,?是

到较高的分析精度。在三维分析中的某些情况下，

约束的拉格朗日乘子,F是广义的激励力[4][5]。

[4][5]。

利用非协调模式可以大大减少分析的成本[8]。综上

1.2方程解法探讨

所述，使用六面体非协调单元对抽壳系统进行