基于有限元的某多管火箭炮模态分析

基于有限元的某多管火箭炮模态分析
弹箭与制导第3O卷
基于有限元的某多管火箭炮模态分析
徐立黄,张贵林,张波涛,谢家松
(1南京理工大学机械工程学院,南京210094;2湖北江山重工有限责任公司,湖北襄樊441003)
摘要:为掌握某多管火箭炮的发射动力学特性,文中利用CAD软件Pro/E建立某多管火箭炮的简化实体
模型,然后将其导人ANSYS有限元软件建立有限元模型并计算,得到了详细的模态固有频率和振型.为某多
管火箭炮的振动特性分析和结构动力特性的优化设计提供了依据.
关键词:多管火箭炮;有限元;模态参数
中图分类号:TJ393文献标志码:A TheModalAnalysisofaMLRSBasedonFEM
XUI.ihuang,ZHANGGuilin,ZHANGBotao.,XIEJiasong
(1SchoolofMechanicalEngineering,NUST,Nanjing210094,China;
2JiangshanHeavyIndustriesCo.Itd,HubeiXiangfan441003.China)
Abstract:TodeterminethedynamicpropertiesofaM1RS,thethree—dimensionalCADmodelofthesystemwascreatedinPro/E andthentheformerwasintroducedintoANSYS,inwhichitcreatedtheparticularfiniteeleme ntmodeiofthesystemandthemo—dalnaturalfrequenciesandvibrationmodewereobtained.Theresultsareveryusefulforopti maldesign.
Keywords:MIRS;finiteelementmethod;modalparameter
IJ5l昌
火箭炮是现代战争中重要的压制性武器,但受射
击精度影响,限制了其应用范围.影响射击精度的主
要因素之一是发射系统的振动,设法减少火箭炮发射瞬问的振动,将有效的减小火箭弹起始段的扰动,提高火箭炮的射击密集度.要减小发射系统振动对射击精度的影响,甚至利用发射系统的振动补偿其他扰动因素造成的不利影响,首先必须了解和掌握其固有振动特性,以利于火箭炮结构的优化设计和改进.通过ANSYS对火箭炮进行模态分析,得到结构的固有频率和振型,能为某多管火箭炮的振动特性分析和结构动力学特性的优化设计提供依据.
l模态分析的基本理论
模态分析理论吸取了振动理论,信号分析,数据
处理,数理统计及自动控制理论中的有关"营养",结合自身内容的发展,形成了一套独特的理论.自动控制中的传递函数概念的引入,对模态分析理论的发展起着很大的推动作用.
模态分析的核心内容是确定描述系统动态特性
的参数.对于一个有N个自由度的线性系统,其运
动微分方程为:
』++KX=F(f)(1)
对于自由振动(F(f)=0),上式可写为:
M+CX+KX=0(2)
若忽略阻尼影响,式(2)可进一步简化为:
^+c:k=0(3)
式中:M为质量矩阵;c为阻尼矩阵;K为刚度矩阵;X 为位移向量;F(f)为作用力向量;t为时间.
自由振动时,结构上各点作简谐振动,各节点位
移为:
X一一e(4)
南式(3),式(4)得:
(K一.M)=0
由该式可求出"个特征值,CO…,和相对
应的"个特征向量,,…,…其中特征值的平
方根,(i—l,2,…,)就是结构的第阶固有频率,
特征向量,就是结构的第阶模态振型,它反映了结
构按频率co,振动时各自由度方向振幅问相对比例关
系.
*收稿日期:2009—1214
作者简介:徐立黄(1984一),男,江苏泗阳人,硕十:研究生,研究方向:兵器发射理论与技术.
第6期徐立黄等:基于有限元的某多管火箭炮模态分析
2建立有限元模型
由于火箭炮结构复杂,在ANSYS中直接建模难
度较大,文中选用Pro/E建立实体模型,然后直接导
人ANSYS.合理简化后的实体模型见图1,有限元
模型见图2.
图1物理模型图2有限元模型
建立物理模型时,需对火箭炮进行合理简化,其
原则是简化不能影响结构的实际动力学特性.遵循
这一原则,忽略结构中的一些圆角,倒角,螺纹孔等.
考虑射击时回转与俯仰方向均被锁定,因此将箱体与
回转体之间用于调整俯仰角的扇形齿轮简化为刚性
连接,耳轴与回转体之间的轴承连接也简化为刚性
连接.
整个模型采用由20个节点定义的SOIID95单
元,它能够吸收不规则形状的单元而精度没有损失.
SOIID95单元有可并立的位移形状并且对于曲线边
界的模型能很好的适应,因此非常适合于复杂模型和
由CAD导入的模型.模态分析对网格质量的要求与
静力学分析并无大的区别,但可以忽略应力集中处对网格的细化处理,更加注重网格的均匀性,这是因为固有频率和振型主要取决于结构质量分布和刚度分布,不存在类似应力集中的现象,采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差太大,可减小数值计算误差.
定向管为薄壁长圆管,若采用自由网格划分则产
生的单元数太多,计算耗时,甚至导致一般计算机无法进行计算.为减少单元和节点,定向管采用扫掠法划分网格,为保证与支撑板部分网格正确连接,将定向管沿支撑板面切割后再划分网格.其余部分结构不规则,均采用自由网格划分.对2O节点单元退化为四面体单元的,将其转化为1O节点的SOlID92单元.通过划分网格,共产生45460个单元,166075个节点.
整炮均为钢材,其密度为7830kg/m.,弹性模量
为209GPa,泊松比为0.3.模拟发射角为45..这里
需要强调的是ANSYS中单位制(Units)问题,AN—SYS中有一个只能从命令行输入的命令:"/UNIT, SI",它的作用仅是标记作用,让用户有个地方做标记,它没有任何单位转换的功能.这要求单位由用户去统一,这类似于振动实验中的标定问题.如果在Pro/E中建模的长度单位为mm,这与密度单位
kg/m.不符合.若在材料属性中直接输入7830,必将
导致模型质量太大,模态固有频率太小.有人提出对密度单位进行换算.由牛顿第二定律F—Ma可知,
1N一1kg*m/s,那么质量单位应该用吨,相应的密
度单位就应该是吨每立方毫米了.但问题并未因此
得到完美的解决,因为这时的弹性模量也需要调整, 否则结果依旧会失常.文中采用在Pro/E中建立模
型的时候把单位改成国际单位制,同时,在装配时调入模型之前先把装配的模板的单位改成国际单位制, 这样,调入ANSYS中的单位就没有问题了.
3模态计算及分析
有限元模型建好后,进入ANSYS求解器,选择
分析类型为Modal,在求解选项(GUI—AnalysisOp—tions)中指定模态提取方法为BlockLanzcos.这种
模态提取方法采用Lanzcos算法,用一组向量来实现Lanzcos递归计算,同Subspace法相比,计算精度相当,但速度要快,存储要求较低,只是所需内存比Subspace法多50%.无论用EQSIV命令指定过哪
一
种求解器进行求解,BlockLanzcos法都将自动采
用稀疏矩阵方程求解器.当计算某系统特征值谱所
包含的一定范围的固有频率时,采用这种方法特别有效.计算时,求解从频率谱中间位置到高频范围内的固有频率时的求解收敛速度和求解低阶频率时基本相同,因此当采用频移频率FREQB来提取从FREQB起始的n阶模态时,该法提取大于FREQB 的n阶模态和提取n阶低阶模态的速度基本相同. 同Subspace法一样,BlockIanzcos法特别适用于大
型对称特征值求解问题.点击确定按钮,在随后弹出的对话框中设置提取模态数No.ofmodestoextract
为12.为了在后处理器中查看振型,必须进行模态
扩展(也就是将振型写入结果文件).将Expansion PassOn/0ff设置为ON.在NumberofModestO Expand中指定要扩展的模态数.这里将其设置为
12(与提取模态数相同),设置完成后,通过Main Menu-Solution—CurrentLS对所建立的有限元模型
进行计算.
(下转第115页)
第6期王晓方等:某型电子时间引信电参数综合测试系统设计?115? 一图5测试系统丁作过程中
5结论
文中在分析了引信检测电参数的基础之上,设计
了一种基于单片机的引信电参数检测系统,系统结构
简单,具有很好的人机交互性能,使用灵活方便,实现
(上接第1l1页)
通过计算,得到了模型的前12阶模态,由于从第
6阶模态起出现局部模态,且为集中在定向器上的高
阶模态,因此这里只对前6阶主要模态进行分析,详
见表1和图3～图8,图中位移仅表示振型,非实际位
移尺寸.
表1固有频率和振型
■一同3一阶模态振型图4二阶模态振型
一■同5阶模态振型罔6四阶模态振型
了生产过程中引信电子部件电参数的现场检测,具有
较高的经济价值.
参考文献:
EG徐爱钧.智能化测量控制仪表原理与设计[M].北京:北
京航空航天大学出版社,1995.
E23李彦学,智墩旺.无线电与电子时间引信EM].北京:兵器
工业出版社,1996.
[3]李学海.PIC单片机实用教程[M].北京:北京航空航天
大学出版社,2006.
[4]向红军,雷彬.便携式表面粗糙度测试仪设计[J].计算
机测量与控制,2006,14(5):692—694.
图7五阶模态振型图8六阶模态振型
4结论
1)该模态分析结果与同类武器(参考文献E4-I)相
符合,证明计算结果的正确性.因此可作为谐谱响应
分析,频谱分析等其他动力学分析的基础.可以为该
火箭炮结构动力特性的优化设计提供理论依据,从而缩短研发周期,降低研发费用.
2)火箭炮发射时间间隔一般约为o.5s,即频率为
2Hz左右,该发射系统主振频率较好的避开了2Hz
及其邻近倍频.但考虑偏航和俯仰振动对射击精度
的敏感性,建议加强摇架的扭转刚度和箱体俯仰连接刚度.如可将箱体底部加厚钢板改为钢架结构,从而
减轻质量,提高固有频率.
参考文献:
[1]傅志方.振动模态分析与参数辨识[M].北京:机械工业出版社,1990.
E2]刘相新盂宪颐.ANSYS基础与应用教程[M].北京:科学出版社,2005.
[3]张琪,刘莉.导弹固有特性的有限元分析[J].弹箭与制导,2008,28(2):6—7.
E4]潘宏侠,王福明.某火箭炮模态试验分析及动态修改[J].兵工武器分册.1993(3):58—64.。