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

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

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

宽屏显示 收藏 分享

资源预览需要最新版本的Flash Player支持。
您尚未安装或版本过低,建议您

专业文档,值得下载专业文档,值得珍藏文章编号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值的增大而变化不大。
编号:201311161020074308    大小:377.00KB    格式:DOC    上传时间:2013-11-16
  【编辑】
5
关 键 词:
精品文档
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 人人文库网仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
  人人文库网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

当前资源信息

4.0
 
(2人评价)
浏览:27次
dingyx0101上传于2013-11-16

官方联系方式

客服手机:13961746681   
2:不支持迅雷下载,请使用浏览器下载   
3:不支持QQ浏览器下载,请用其他浏览器   
4:下载后的文档和图纸-无水印   
5:文档经过压缩,下载后原文更清晰   

相关资源

相关资源

相关搜索

精品文档  
关于我们 - 网站声明 - 网站地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2015-2017 人人文库网网站版权所有
苏ICP备12009002号-5