会员注册 | 登录 | 微信快捷登录 QQ登录 微博登录 | 帮助中心 人人文库renrendoc.com美如初恋!
站内搜索 百度文库

热门搜索: 直缝焊接机 矿井提升机 循环球式转向器图纸 机器人手爪发展史 管道机器人dwg 动平衡试验台设计

   首页 人人文库网 > 资源分类 > DOC文档下载

基于Matlab 挠性转子动力学分析的软件设计.doc

  • 资源星级:
  • 资源大小:377.00KB   全文页数:5页
  • 资源格式: DOC        下载权限:注册会员
您还没有登陆,请先登录。登陆后即可下载此文档。
  合作网站登录: 微信快捷登录 支付宝快捷登录   QQ登录   微博登录
友情提示
2:本站资源不支持迅雷下载,请使用浏览器直接下载(不支持QQ浏览器)
3:本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

基于Matlab 挠性转子动力学分析的软件设计.doc

专业文档,值得下载专业文档,值得珍藏文章编号16718178(2003)03006205基于Matlab挠性转子动力学分析的软件设计1程胜文2杨光(1湖北职业技术学院,湖北孝感4320002武汉理工大学,湖北武汉430074)摘要文章主要利用传递矩阵算法对复杂挠性转子系统进行动力学分析的相关理论,在MATLAB平台上,开发了一个挠性转子动力特性分析软件,而且利用MATLAB的GUI编辑多个人机界面,可方便实现结果的可视化。关键词挠性转子动力学分析软件设计中图分类号TH113.2文献标识码A1引言迄今为止,转子动力学的研究已有100多年的历史。早期的转子动力学研究中,转子模型较为简单,而且转子的转速长期限制在一阶临界转速以下工作1。50年代以来,旋转机械的转子越来越柔、功率越来越大、转速越来越高,甚至达到了三、四阶临界以上。这就使得挠性转子系统动力学的研究显得尤为重要。常用的对转子进行动力学分析的方法有传递矩阵法和有限元法。传递矩阵法具有占用储存空间小,计算速度快,矩阵维数恒定,能计算至任意高阶固有频率,无需预知振型等优点,笔者仅限于利用传递矩阵法对挠性转子动力学分析的软件开发。MATLAB是一个强大的数学计算和仿真工具,利用它我们可以避免复杂的数学计算编程,并且借助其绘图函数,可方便实现计算的可视化和仿真。与Basic、Fortran以及C语言相比较,MATLAB的语法规则更简单,编程特点更贴近人的思维方式。用MATLAB写程序有如在便笺上列公式和求解,因而MATLAB被称为科学便笺式的科学工程计算语言。22传递矩阵法求解各阶临界转速和振型22.1系统的动力模型本文中笔者研究的是多圆盘、多支承的挠性转子动力学系统。对于这类复杂的挠性转子系统,要进行机械系统精确的弹性动力分析是十分困难的。因而,必须对实际的工程问题进行简化和抽象,做一个繁简适度的力学描述和数学描述,即建立其动力模型。简化后建立图1所示模型。将连续的挠性转子系统简化为薄圆盘和无质量的弹性轴段。按照当量轴计算出来的结果与实际情况会有所差异。分段多一些可以接近实际情况。但是过多的分段计算误差会增加,而且当分段达到一定数目后,由于其他参数而产生的影响,可能会超过由于分段造成的影响。所以分段并非越多越好。程胜文杨光基于Matlab挠性转子动力学分析的软件设计63在笔者设计的软件中将考虑由于圆盘偏摆和动刚度的变化对转子系统临界转速的影响。上图2所示模型中有陀螺力矩npnngJHHM,方向与平面O′AB垂直,大小为sinnpgJM由于陀螺力矩的存在,当正进动40ψπ/2的情况下,它使转轴的变形减小,因而提高了转轴的弹性刚度,即提高了转子的临界角速度。在反进动π/2ψπ的情况下,这力矩使转轴的变形增大,从而降低了转轴的弹性刚度,即降低了转子的临界速度。轴承支承刚度是由轴承座的刚度和油膜的刚度的结合。总刚度k为22bbpbbpmkkmkkk(1)其中pk为油膜刚度,bk,bm分别为轴承座的等效刚度和质量。2.2传递矩阵法对于第i个部件,如其左右两端截面的编号分别为i与1i,则由截面i的挠度iy,斜率i,弯矩iM及剪力iQ所组成的列阵,称为该截面的状态向量iz。即TiiQMyz,,,(2)任一部件两端截面的状态向量总存在一定的关系,即iiizTz1(3)iT称为该部件的传递矩阵。当状态向量有r个元素时,iT是rr阶方阵。其中iT为idpdpdpiKmlJJKmlEIlEIlJJEIlKmEIlEIlEIlJJEIllKmEIlT1001212622612222222322223(4)式(4)中dJlm,,,pJ分别为圆盘质量、轴段长度、极转动惯量、径转动惯量。对于图1的模型,若边界条件设置为1M=01Q=05若设有矩阵121TTTAiii(i2,3,4,N1)6因此对于任一截面ii2,3,,N1有4131211111}{}{aaaazAzii42322212aaaa43332313aaaaiaaaa44342414100y41312111aaaa142322212iaaaa1y7式中ikai,k1,2,3,,N1表示第i段质心到第k个支承之间的距离,a即为转子的总跨度。在利用(6)式计算iA时,在各矩阵的连乘过程中,图1转子的当量模型图2圆盘的偏摆湖北职业技术学院学报2003年第3期第6卷第3期64只需计算第一、二列的元素。对于末端截面N1,由(7)式知有1NQM4131aaNaa42321y8在末端截面N1的边界条件是0,011NNQM9我们称△(ω2)4131aaNaa4232010为剩余量。由(8)式知,当边界条件(9)式满足时,这一齐次式有解的条件是剩余量为零上式就是转子系统的频率主程。当某阶临界转速求得后,由(7)式可得Niyaaaaaayaaaayiii,,3,2122211211112221121111由此即可求得各截面的y与的比例解,亦即对应于该阶临界转速的模态振型。3软件的设计3图7是该挠性转子动力分析软件的程序设计框图。该软件是在Matlab平台上设计开发的。根据前面的计算思路,利用MATLAB编制了一个求解多圆盘、多支承转子临界速度的通用计算分析软件。利用该软件计算转子的临界速度时,在界面上最多显示6阶的数值。图3是软件的主界面,主界面右上角显示余量2的关系图。图4、图5分别为总体参数和轴段参数输入界面。限于篇幅的限制,本文中不给出具体的程序。图3程序主界面图4基本数据输入界面图5轴段参数输入界面图6单圆盘转子系统简化图程胜文杨光基于Matlab挠性转子动力学分析的软件设计65图7传递矩阵法计算框图同时软件有一个简单的分析功能分析轴的刚度和支承的刚度对转子临界转速的影响。现以单圆盘转子为例分析系统刚度对系统临界转速的影响。模型如图6,现分析,当支承刚度不变时,轴的抗弯刚度EI对转子系统临界转速的影响。设K150000000N/m。编程分析如下图8所示。由图可以看出,随着轴的刚度的增大,转子系统的临界转速也随着增加。而当轴的刚度设为不变时,图8反映了支承刚度的变化对转子临界转速的影响。图8是设轴的刚度EI8000N.m2时的变化曲线图。由图可以看出,当支承刚度K由一个较小值变到2000000N/m时,转子固有频率随K值的变化发生很大的变化,当K值到达一定值时,转子的固有频率随支承刚度K值的增大而变化不大。

注意事项

本文(基于Matlab 挠性转子动力学分析的软件设计.doc)为本站会员(dingyx0101)主动上传,人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知人人文库网([email protected]),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。

copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5