转子动力学大作业
转子动力学大作业
学院：
姓名：
班级：
学号：
目录
一、作业题目介绍
二、转子动力学理论简介
三、参数的选择和计算
四、Ansys分析临固有频率和临界转速
五、失稳转速影响因素及计算
一、大作业题目
1、 计算临界转速；
2、 圆轴承，长颈比为0.8，油膜间隙2‰
3、 计算失稳转速
注：转子两端各一个轴承，支点在左右两端。
二、转子动力学理论知识
由于制造中的误差，转子各微段的质心一般对回转轴线有微小偏离。因此，当转子转动时，会出现横向干扰，在某些转速下还会引起系统强烈振动，出现这种情况时的转速就是临界转速。为保证系统正常工作或避免系统因振动而损坏，转动系统的转子工作转速应尽可能避开临界转速，若无法避开，则应采取特殊防振措施。这也是研究临界转速的意义。临界转速和转子不旋转时横向振动的固有频率相同，也就是说，临界转速与转子的弹性和质量分布当圆盘不装在两支撑的中点而偏于一边时，转轴变形后，圆盘的转轴线与两支点A 和B 的连线有一夹角ψ。设圆盘的自转角速度Ω，极转动惯量为p J ，则圆盘对质心o '的动量矩为p H J =Ω。它与轴线AB 的夹角也应该是ψ，见图1。当转轴有自然振动时，设其频率为n ω，则圆盘中心o '与轴线AB 所构成的平面绕AB 轴有进动角速度n ω。由于进动，圆盘的动量矩H 将不断改变方向。因此有惯性力矩
()g n n p n M H H J ωωω=-?=?=Ω?
方向与平面o AB '垂直，大小为
sin g p n M J ωψ=Ω
转子结构尺寸示意图
这一惯性力矩称为陀螺力矩或回转力矩。因夹角ψ较小，sin ψψ≈，上式可写作
g p n M J ωψ=Ω。
这一力矩与ψ成正比，相当于弹性力矩。在正进动(0/2ψπ<<)的情况下，它使转轴的变形减小，因而提高了转轴的弹性刚度，即提高了转子的临界角速度。在反进动（/2πψπ<<）的情况下，这力矩使转轴的变形增大，从而降低了转轴的弹性刚度，即降低了转子的临界角速度。通过分析，可知道陀螺力矩对转子临界转速的影响：正进动时，它提高了临界转速；反进动时，它降低了临界转速。
图 1
在大多数情况下，轴承对于转子的动力特性有很明显的影响，轴承往往是阻尼的主要来源，因而控制着转子的响应。同时，轴承的刚度和阻尼又影响着转子的临界转速和稳定性。在深入研究转子动力学问题时，因而必须考虑到轴承的作用。对于一个确定的轴承，当润滑油粘度及进油压已给定时，轴颈中心1o 的静平衡位置e 、?决定于轴颈转速Ω和静载荷W 。当载荷W 的大小或者轴颈转速Ω变化时，1o 位置也相应地变化。当铅垂载荷W 大小变化时，轴颈中心的移动在大多数情况下，并非沿铅垂方向，也即位移并不沿着载荷作用的方向。这正是油膜不同于一般机械元件的一个特点。
记x F 、y F 为油膜力在x 、y 方向的分量。我们定义油膜刚度系数为单位位移所引起的油膜力增量，即
x
xx F k x
?=?，0
y xy F k y
?=
?，0
y yx F k x
?=
?，0
y yy F k y
?=
?
定义油膜阻尼系数为单位速度所引起的油膜力增量，即
x xx F c x
?=
? ，0
x
xy F c y
?=
? ，0
y
yx F c x
?=
? ，0
y
yy F c y
?=
?
式中各系数的第一个下标代表力的方向，第二个下标代表位移或速度的方向。油膜刚度系数和阻尼系数统称为油膜动力特性系数。其中xy k ，yx k 和xy c ，yx c 分别称为交叉刚度系数和交叉阻尼系数，它们表示油膜力在两个相互垂直方向的耦合作用，交叉动力系数的大小和正
负在很大程度上影响着轴承工作的稳定性。 轴颈中心o '偏离轴承中心o 时，轴颈和轴承的间隙沿周向是不均匀的。润滑油被轴颈带动，顺着转动方向从较宽的间隙流进较窄的间隙而形成油楔，对轴颈有挤压力作用。当润滑油从较窄的间隙流到较宽的间隙时，因出现空穴而对轴有负压力。轴承的全部油膜对轴颈的总压力F 位于挤压的一侧并朝向轴颈中心o '。将F 力分解为o '点的径向力e F 和周向力F ψ。
分力e F 起支撑轴颈的作用，相当于转轴的弹性力。分力F ψ垂直于o '的向径并顺着转动方向，是
o '的速度增大，因而向径oo '增大。F ψ就是使轴颈运动失稳的力。
三、参数的选择与计算
初取0
52T C =，360/rad s Ω=，11
78000.01561060922
W Vg N ρ=
=???= /257.296(1/)N s π=Ω=,20.1D R m ==,0.80.08L D m ==,
40002210C D m -==?，选择粘度牌号为22号的润滑油，则214/mm s υ=。
动力粘度 6321410900/ 1.2610kg m Pa s ηυρ--==??=?
036301.88/()900/ 1.69210/()v c c kJ kg C kg m J m C ρ==??=??
22
2
41.261057.2960.10.080.05(()0.428609210
NDL R
S W C η--????=
=?=? 线性插值的油膜参数0.3120,0.8343, 6.4690Q P T ===
此时20=+0.8
()43018.016748.0167v R
Q
T T C c C P
ηπΩ=+= 供油工作 020max =48.0167()48.0167 4.335652.3523v R
T T T C T C c C
ηΩ+?=+
=+=工作
由于最高问题max T 跟假设的问题T 相近，故假设成立。 得轴承的动力特性系数：
1.9234,
2.6687, 6.0205, 1.6369xx xy yx yy k k k k ===-= 5.6411, 2.0471,11.8940xx xy yx yy c c c c ===-=
因为：,W W k k
c c C C ==Ω
66765.87510,8.12610, 1.810, 4.984410xx xy yx yy k k k k =?=?=-?=? 4454.771410, 1.731510, 1.00610xx xy yx yy c c c c =?==-?=? 四、Ansys 分析临固有频率和临界转速
1、 建立beam188和combination214单元，取
66765.87510,8.12610, 1.810, 4.984410xx xy yx yy k k k k =?=?=-?=? 4454.771410, 1.731510, 1.00610xx xy yx yy c c c c =?==-?=?
Campbell 图
若取刚度为87/xx
yy k k e N m ==，0/xy yx k k N m ==，13/xx yy c c e Ns m ==，
0/xy yx c c Ns m ==
可得Campbell 图
可见刚度系数对转子的固有频率及稳定性有着很大的影响。
用计算出的刚度系数，通过ANSYS计算得临界转速（r/min）
1 1831.943
2 1861.444
3 3451.431
4 3572.778
5 14317.884
6 15025.431
7 31665.300
8 33714.650
固有频率及模态振型：
一阶：
二阶：
三阶：
五、失稳转速影响因素及计算
转子运动失稳的原因除转轴的材料内阻或者圆盘与转轴配合面的摩擦以外，还有滑动轴承的油膜力。因间隙不均匀而引起的横向气动力等等。而对于大型汽轮机发电机组的转子或高速转子，轴承的油膜力引起转子运动失稳现象比较常见。
根据经验公式可以知道，失稳角速度: 11831.922=3663.886 43i c n n rpm ≈=?
机械系统动力学
机械系统动力学报告 题目：电梯机械系统的动态特性分析 姓名： 专业： 学号：
电梯机械系统的动态特性分析 一、课题背景介绍 随着社会的快速发展，城市人口密度越来越大，高层建筑不断涌现，因此，现在对电梯的提出了更高的要求，随着科技的进步，在满足客观需求的基础上，电梯向着舒适性，高速，高效的方向发展。在电梯的发展过程中，安全性和功能性一直是电梯公司首要考虑的因素，其中舒适性也要包含在电梯的设计中，避免出现速度或者加速度出现突变，或者电梯运行过程中的振动引起人们的不适。因此，在电梯的设计过程中，对电梯进行动态特性分析是十分必要的。 二、在MATLAB中编程、绘图。 通过同组小伙伴的努力，已经得到了该系统的简化模型与运动方程。因此进行编程： 该系统的微分方程：[][][]{}[]Q x k x c x M= + ? ? ? ? ? ? + ? ? ? ? ? ?? ? ? ，其中矩阵[M]、 [C]、[K]、[Q]都已知。 该系统的微分方程是一个二阶一元微分方程，在MATLAB中，提供有求解常微分方程数值解的函数，其中在MATLAB中常用的求微分方程数值解的有7个：ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb 。 ode是MATLAB专门用于解微分方程的功能函数。该求解器有变步长（variable-step）和定步长（fixed-step）两种类型。不同类型有着不同的求解器，其中ode45求解器属于变步长的一种，采用Runge-Kutta
算法；和他采用相同算法的变步长求解器还有ode23。 ode45表示采用四阶，五阶Runge-Kutta单步算法,截断误差为(Δx)^3。解决的是Nonstiff(非刚性)常微分方程。 ode45是解决数值解问题的首选方法，若长时间没结果，应该就是刚性的，可换用ode23试试。 Ode45函数调用形式如下：[T,Y]=ode45(odefun,tspan,y0) 相关参数介绍如下： 通过以上的了解，并对该微分方程进行变换与降阶，得出程序。MATLAB程序： （1）建立M函数文件来定义方程组如下： function dy=func(t,y) dy=zeros(10,1); dy(1)=y(2); dy(2)=1/1660*(-0.006*y(2)+0.003*y(4)-0.0006*y(10)-1.27*10^7*y(1)+1.27*10^7*y (3)+2.54*10^6*y(9)); dy(3)=y(4); dy(4)=1/1600*(+0.03*y(2)-0.007*y(4)+0.003*y(6)+1.27*10^7*y(1)-7.274*10^8*y(3 )+1.27*10^7*y(5)); dy(5)=y(6);
“机械动力学”课程教学大纲
“机械动力学”课程教学大纲 英文名称：Mechanical Dynamics 课程编号：MACH3441 学时：32 （理论学时：32 实验学时：课外学时：2实验） 学分：2 适用对象：机械设计、机械制造及自动化、机械电子工程、流体机械、电机、电器、材料工程等本科生高年级。 先修课程：高等数学、普通物理学、理论力学、材料力学、线性代数使用教材及参考书： [1] 石端伟主编. 机械动力学. 北京：中国电力出版社，2007. [2] 张策主编. 机械动力学.北京：高等教育出版社, 2008. [3] 倪振华主编. 振动力学. 西安交通大学出版社，1988. 一、课程性质和目的 性质：专业课 目的： 1.了解机械动力学的研究内容、发展历史以及最新研究进展。 2.培养机械系统动力学分析的基本能力。 3.了解机械系统动力学分析相关的CAE软件。 4.了解机械系统动态测试有关技术。 5.培养查阅和运用相关科技文献进行动力学分析的初步能力。 6.培养创新思维以及解决工程实际问题的能力。 7.培养科学、严谨的工作作风。
二、课程内容简介 随着现代机械装备朝着高精度、高效、大功率的方向发展，其动态性能指标的优劣越来越受到广泛关注和高度重视。机械动力学已日益成为现代机械设计与制造工程领域不可或缺的基础知识。本课程主要介绍机械系统动力分析的基本理论、分析方法、测试与控制技术以及典型机械系统动力学分析方法。通过课程的学习，培养学生能够在机械系统动力分析方面具有明确的基本概念、必要的专业基础知识、一定的机械系统动力分析能力与计算能力。 三、教学基本要求 1．了解相关机械系统动力学分析的新理论、新方法及发展趋向。 2. 掌握有关机械系统动力学分析的基本概念、基本理论与方法。 3. 了解典型机械系统动力学分析流程，具有进行工程实际问题分析的初步能力。 4. 建立正确的机械系统动力分析的思维方式，理论联系实际，具备一定的科研创新精神； 5. 课后需要查阅文献，并开展讨论，完成作业。 四、教学内容及安排 第一章：绪论 1．熟悉研究机械动力学的意义。 2．熟悉机械动力学的主要研究内容。 教学安排及教学方式
Ansys转子动力学
基于ANSYS的转子动力学分析 1、题目描述 如图1-1所示,利用有限原原理计算转子临界转速以及不平衡响应。 图 1-1 转子示意图及尺寸 2、题目分析 采用商业软件ANSYS进行分析，转子建模时用beam188三维梁单元,该单元基于Timoshenko梁理论，考虑转动惯量与剪切变形的影响。每个节点有6个(三个平动,三个转动)或7各自由度（第七个自由度为翘曲，可选）。 轴承用combine214单元模拟。该单元可以模拟交叉刚度和阻尼。只能模拟拉压刚度，不能模拟弯曲或扭转刚度。该单元如图2-1所示，其有两个节点组成，一个节点在转子上，另一个节点在基础上。
图 2-1 combine214单元 对于质量圆盘，可以用mass21单元模拟，该单元有6个自由度，可以模拟X,Y,Z 三个方向的平动质量以及转动惯性。 3、计算与结果分析 3.1 转子有限元模型 建模时，采用钢的参数，密度取37800/kg m ，弹性模量取112.1110pa ，泊松比取0.3。轴承刚度与阻尼如表1所示，不考虑交叉刚度与阻尼，且为各项同性。 Kxx Kyy Cxx Cyy 4e7N/m 4e7N/m 4e5N.s/m 4e5N.s/m 将转子划分为93个节点共92个单元。有限元模型如图3-1所示。
图 3-1 转子有限元模型 施加约束时，由于不考虑纵向振动与扭转振动，故约束每一节点的纵向与扭转自由度，同时约束轴承的基础节点。施加约束后的模型如3-2所示。 图 3-2 施加约束后的有限元模型 3.1 转子临界转速计算 在ANSYS中可以很方便的考虑陀螺力矩的影响。考虑陀螺力矩时，由于陀螺矩阵是反对称矩阵，所以求取特征值时要用特殊的方法。本文考虑陀螺力矩的影响，分析了在陀螺力矩的影响下，转子涡动频率随工作转速的变化趋势，其Campell图如图3-3所示。同时给出了转子的前四阶正进动涡动频率与反进动涡动频率以及固有频率。如表3-2所示。
matlab机电系统仿真大作业
一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环
3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置
M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置
机械动力学大作业
单自由度杆机构的Adams动力学仿真 摘要：文章分析了单自由度的铰链机构的动力学问题，已知原动件曲柄的转矩，绘制输出件摆杆的运动曲线。首先在Adams软件中构造连杆，添加三个连杆，使其成一定角度，相互连接。再在两杆之间添加转动副，并且头尾连杆与地相连。并在曲柄处加转矩，最后进行仿真，并绘出相应图表。 关键词：铰链机构；Adams仿真 1、机构模型的建立 根据题目要求，选择一个铰链四杆机构——曲柄摇杆机构为模型，其结构简图如图1所示。其中，曲柄1为原动件。 图1曲柄摇杆机构简图 在Adams软件中，建立该曲柄摇杆机构的模型如图2所示。 图2 Adams中的曲柄摇杆机构模型
曲柄摇杆机构各连杆的惯性参数参考表1。杆件的材料均选择钢材（密度ρ=7.801×10-6 kg?mm-3，杨氏模量E=2.07×105 N?mm-2，泊松比μ=0.29）。 表1 传动导杆机构各部件惯性参数 2、利用Adams软件添加约束和力矩 杆1和地之间有转动副，杆1和杆2、杆2和杆3之间有转动副，杆3和地之间有转动副。杆1为原动件，在杆1上添加转矩。转矩大小为30。 图3约束与转矩 3、进行仿真 点击仿真按钮，开始仿真，选择仿真时间为2s，可以观察到该机构各个时间的运动状态如图4和图5所示。
（a）T=0时刻（b）T=1时刻 图4仿真过程中机构模型的运动状态 （a）T=1.2时刻（b）T=2时刻 图5仿真过程中机构模型的运动状态 结论 当原动件曲柄的转矩取为30时，点击“后处理”，可以绘制出输出件摆杆的位移曲线、角速度曲线、加速度曲线分别如图10、图11和图12所示。 图10输出件摆杆的位移曲线
液压伺服 大作业
硕士学位课程考试试卷 考试科目：电液伺服控制 考生姓名：刘双龙 考生学号：20140713189 学院：机械工程学院专业：机械工程 考生成绩： 任课老师(签名) 考试日期：2014年1月20日午时至时
考试主题：电液伺服（比例）系统 考试题目： 1、为什么把液压控制阀称为液压放大元件？ 2、什么叫阀的工作点？零位工作点的条件是什么? 3、电液伺服阀由哪几部分组成?各部分的作用是什么？ 4、什么是液压固有频率?在阀控缸系统中液压固有频率与活塞位 置有关吗?为什么? 5、为什么电液伺服系统一般都要加校正装置? 6、结合自己研究领域，写一篇液压伺服系统建模、分析的论文， 字数不少于2000字。 注：要求独立完成，不允许抄袭。 交作业时间： 最迟2015年第一个学期的第一周交到7教136，交纸质档。
三自由度平台液压伺服系统建模 摘要： 我的专业是机械工程，主要方向是机械设计，所以本文选择了与我专业方向有关的一个机构进行建模。本文开始对机构进行了说明（采用已有的机构，并非自己设计），然后对其进行运动学分析，从而的到上平台和下平台的速度及加速度，和雅可比矩阵及液压缸速度。然后对驱动机构进行电液伺服系统建模。其中 一：自由度运动平台系统简介 本文所研究的三自由度运动平台类似与六自由度平台是由一个上平台(动平台)、地基(下平台)、三个支杆、三个线性作动器以及若干关节连接而成的。上平台装有负载，完成既定的位置、速度、加速度运动要求，进而实现刑于道路状况的复现。其结构示意图如图1．1所示。 图 1三自由度运动平台的结构图 该平台的结构如下：上平台与地面之间以三个支杆(strut)来约束并起支撑作用，并以三个液压缸作为驱动部件进行驱动。每个液压缸两端为关节轴承，中间为一个移动副和一个转动副连接；每根支杆两端也是采用关节轴承分别与地面和上平台相连中间一个转动副。通过计算可知每个支杆所在的支路都具有5个自由度，每个支路对上平台提供一个约束；每个液压作动器所在的支路都具有6个自由度，对于上平台没有约束。通过每个分支对上平台的约束很容易计算得出其自由度为3。因此，通过三套液压作动器的驱动，上平台能够实现对于给定运动的跟踪复现。 简单直观的对运动进行分析可得到：由于三根支杆的限制作用，上平台平动受到限制：而转动自由度相对更为自由，运动范围更大。当两竖直作动器差动动
研究生《机械系统动力学》试卷及答案
太原理工大学研究生试题 姓名： 学号： 专业班级： 机械工程2014级 课程名称: 《机械系统动力学》 考试时间: 120分钟 考试日期: 题号 一 二 三 四 五 六 七 八 总分 分数 1 圆柱型仪表悬浮在液体中，如图1所示。仪表质量为m ，液体的比重为ρ，液体的粘性阻尼系数为r ，试导出仪表在液体中竖直方向自由振动方程式，并求固有频率。（10分） 2 系统如图2所示，试计算系统微幅摆动的固有频率，假定OA 是均质刚性杆，质量为m 。（10分） 3 图3所示的悬臂梁，单位长度质量为ρ，试用雷利法计算横向振动的周期。假定梁的 变形曲线为?? ? ?? -=x L y y M 2cos 1π（y M 为自由端的挠度）。（10分） 4 如图4所示的系统，试推导质量m 微幅振动的方程式并求解θ(t)。（10分） 5 一简支梁如图5所示，在跨中央有重量W 为4900N 电机，在W 的作用下，梁的静挠度δst=，粘性阻尼使自由振动10周后振幅减小为初始值的一半，电机n=600rpm 时，转子不平衡质量产生的离心惯性力Q=1960N ，梁的分布质量略去不计，试求系统稳态受迫振动的振幅。（15分） 6 如图6所示的扭转摆，弹簧杆的刚度系数为K ，圆盘的转动惯量为J ，试求系统的固有频率。(15分) 7如图7一提升机，通过刚度系数m N K /1057823?=的钢丝绳和天轮（定滑轮）提升货载。货载重量N W 147000=，以s m v /025.0=的速度等速下降。求提升机突然制动时的钢丝绳最大张力。(15分) 8某振动系统如图8所示，试用拉个朗日法写出动能、势能和能量散失函数。(15分) 太原理工大学研究生试题纸
车辆系统动力学仿真大作业(带程序)
Assignment Vehicle system dynamics simulation 学院：机电学院 专业：机械工程及自动化 姓名： 指导教师：
The model we are going to analys: The FBD of the suspension system is shown as follow:
According to the New's second Law, we can get the equation: 2 )()(221211mg z z c z z k z m --+-=???? 221212)()(z k mg z z c z z k z m w +-----=? ??? 0)()()()(222111222111=-++--+-++--+? ? ? ? ? ? ? ?w w w w z L z k z L z k z L z c z L z c z m χχχχ 0)()()()(2222111122221111=-++----++---? ? ? ? ? ? ? ?w w w w z L z L k z L z L k z L z L c z L z L c J χχχχχ d w w w w Q z L z k z L z c z m ,111111111)()(-=------? ? ? ? ?χχ d w w w w Q z L z k z L z c z m ,222222222)()(-=-+--+-? ????χχ When there is no excitation we can get the equation: 2)()(221211mg z z c z z k z m --+-=???? 2 21212)()(z k mg z z c z z k z m w +-----=? ??? Then we substitude the data into the equation, we write a procedure to simulate the system: Date: ???? ?? ??? ??==?==?===MN/m 0.10k m 25.1s/m kN 0.20MN/m 0.1m kg 3020kg 2100kg 3250w 2l c k I m m by w b
机械系统动力学试题
机械系统动力学试题 一、 简答题： 1．机械振动系统的固有频率与哪些因素有关？关系如何？ 2．简述机械振动系统的实际阻尼、临界阻尼、阻尼比的联系与区别。 3．简述无阻尼单自由度系统共振的能量集聚过程。 4. 简述线性多自由度系统动力响应分析方法。 5. 如何设计参数，使减振器效果最佳？ 二、 计算题： 1、 单自由度系统质量Kg m 10=, m s N c /20?=, m N k /4000=, m x 01.00=, 00=? x ,根据下列条件求系统的总响应。 （a ） 作用在系统的外激励为t F t F ωcos )(0=,其中N F 1000=， s rad /10=ω。 （b ） 0)(=t F 时的自由振动。 2、 质量为m 的发电转子，它的转动惯量J 0的确定采用试验方法：在转子径向R 1的地方附加一小质量m 1。试验装置如图2所示，记录其振动周期。 a ）求发电机转子J 0。 b ）并证明R 的微小变化在R 1=（m/m 1+1）·R 时有最小影响。 3、 如图3所示扭转振动系统，忽略阻尼的影响 J J J J ===321，K K K ==21 （1）写出其刚度矩阵； （2）写出系统自由振动运动微分方程； （2）求出系统的固有频率； （3）在图示运动平面上，绘出与固有频率对应的振型图。 1 θ（图2）
（图3） 4、求汽车俯仰振动（角运动）和跳振（上下垂直振动）的频率以及振 动中心（节点）的位置（如图4）。参数如下：质量m=1000kg，回转半径r=0.9m，前轴距重心的距离l1=0.1m,后轴距重心的距离l2=1.5m,前弹簧刚度k1=18kN/m,后弹簧刚度k2=22kN/m （图4） 5、如5图所示锻锤作用在工件上的冲击力可以近似为矩形脉冲。已知 工件，铁锤与框架的质量为m1=200 Mg，基础质量为m2=250Mg，弹簧垫的刚度为k1=150MN/m，土壤的刚度为k2=75MN/m.假定各质量的初始位移与速度均为零，求系统的振动规律。
重庆大学机器人大作业
IRB 7600大功率机器人运动仿真
摘要 (2) 1引言 (3) 2机器人发展概述 (3) 2.1机器人的三大定律产生 (3) 2.2工业机器人的发展和特点 (3) 2.3工业机器人现状与前景 (5) 3 ABB机器人和大功率机器人的发展概述 (5) 3.1 ABB公司的发展 (6) 3.2 ABB工业机器人的现状 (6) 3.3简述IRB 7600机器人特点 (6) 3.４IRB 7600机器人的主要参数和应用 (7) 4. 基于ADAMS的IRB 7600大功率机器人运动学仿真 (8) 4.1 IRB 7600大功率机器人的运动学分析 (8) 4.1.1分析IRB 7600大功率机器人得到简图，建立方程 (9) 4.1.2 IRB 7600大功率机器人正向运动学解 (11) 4.2ADAMS中的的运动仿真 (12) 4.2.1在ADAMS中建立IRB 7600机器人的模型 (12) 4.2.2运动的施加 (14) 4.2.4运动结果分析 (16) 总结 (19) 参考文献 (20)
摘要 现代机器人技术飞速发展，其中工业机器人的应用也越来越广泛，成为高科技中极为重要的组成部分。本文主要针对ABB机器中的IRB 7600大功率机器人，对其运动进行仿真探究，学习机器人的一般运动方法。 ABB大功率机器人系列开辟了全新的应用领域，该机器人有多种版本，最大承重能力高达650kg。IRB 7600适合用于各行业重载场合，大转矩、大惯性、刚性结构以及卓越的加速性能等优良特性使这款市场主导产品声誉日隆。用于装配、清洁/喷涂、切割/去毛刺、研磨/抛光、机械管理、物料搬运、货盘堆跺、扳弯机管理、点焊，应用前景广。通过对IRB 7600的模型建立，基于ADAMS的点焊机器人运动学仿真，得到了机器人的仿真运动曲线和模型图。对模型和曲线分析，初步的了解到大功率机器人的运动和工作方式。 关键字：IRB 7600、ABB、ADAMS、仿真
非线性转子 动力学
航空发动机非线性转子碰磨研究 XXX （XXXX 机械工程上海200072） 摘要：综述了国内外非线性转子动力学的研究现状,讨论了非线性转子动力学研究中的7个主要问题,并引述了大量相应的国内外文献,包括:非线性转子动力学研究的一般方法;求解非线性转子动力学问题的数值积分方法;大型转子-轴承系统高维非线性动力学问题的降维求解;基于微分流形的动力系统理论方法;转子非线性动力学行为的机理研究和实验研究;高速转子-轴承系统的非线性动力学设计,最后讨论了非线性转子动力学研究中存在的问题及展望。 关键词：非线性；高速转子；数值积分法 The research for Aeroengine nonlinear rotor WANG Qing-long (Shanghai university mechainal engineering 20072 shanghai) Abstract: Reviewed the research status of nonlinear rotor dynamics both at home and abroad, discusses the seven main in the study of nonlinear rotor dynamics. To questions, and cited a large number of relevant literature both at home and abroad, include: common methods of nonlinear rotor dynamics; To solve the non-linear. Rotor dynamics problems of numerical integral method; Rotor - bearing system of large dimension reduction solution for high dimensional nonlinear dynamics; In the theory of differential dynamic system of the manifold method; Rotor nonlinear dynamics behavior of mechanism research and experiment research; High speed rotor shaft. Bearing system of the nonlinear dynamics design, and finally discusses the problems of nonlinear rotor dynamics research and prospects. Key words: nonlinear; High speed rotor; The numerical integral method. 由于旋转机械系统中各种异常振动的存在,常常引发灾难性的事故。过去研究转子-轴承-基础系统大多采用基于线性转子动力学理论。例如传统转子动力学对转子-轴承系统稳定性问题的研究,一般采用8个线性化的刚度与阻尼特性系数的油膜力模型。对于大型旋转机械中存在的油膜力、密封力、不均匀蒸汽间隙力等严重的非线性激励源,由于数学模型不够完善,以致系统中存在的许多由非线性因素引起的多种复杂动力学行为尚没有彻底搞清,不能满足现代工程设计的需要,迫切需要建立转子-轴承系统的非线性动力学理论,揭示系统存在的各种非线性动力学行为,提出转子-轴承系统的非线性动力学设计方法,研究旋转机械中存在的各种实际问题,这对提高旋转机械运行的稳定性、安全性、可靠性具有重要的现实意义和实际工程背景。 随着非线性动力学理论的发展,非线性转子动力学理论和方法也受到了关注,大量的研究成果使转子动力学面貌一新。但现有的非线性动力学理论和方法在解决高维动力系统方面还存在困难,而工程实际中的转子-轴承-基础系统是一个复杂的高维系统,从而吸引了更多的研究者从事这方面的研究,特别是现代非线性动力学理论在转子动力学中的应用,已成为当今国
机械动力学大作业
机电工程学院有限元分析及应用Ansys软件大作业 学号：S314070061 专业：机械工程 学生姓名：郭海山 任课教师：钟宇光 2014年12月18日
一.题目要求： 采用ADAMS软件或Matlab/Simulink环境，建立简单机械系统的动力学模型，借助软件进行求解计算和结果分析。 建立单自由度杆机构（有无滑块均可）动力学模型，由静止启动，选择一固定驱动力矩，，具体机构及参数自拟。 二.模型及结构分析： 利用ADAMS建立如下图1所示单自由度机构模型： 图1单自由度机构模型 结构简图如下图2： 图2 机构简图 曲柄1长度为24cm，质量为1.69kg 滑块2质量为15.6kg 导杆3长度为80cm，质量为5.19kg
部件的材料都是钢， Material Density: 7.801E-006 kg/mm**3 三.建模： 1.启动adams/view，新建模型model_1。单位设置成MMKS-mm,kg,N,s,deg。存储位置设在桌面。设置工作环境后，利用主工具箱里的基本建模工具，先后建立曲柄1、滑块2和导杆3。 2.曲柄和地面之间，曲柄和连杆之间，连杆和滑块之间，都是转动副。滑块和地面之间是移动副。在Ａ，Ｂ，Ｃ分别放，再在B点添加进行约束。 3.现在给曲柄一个匀速转动。其值如下图3所示： 图3 最后得到模型如下图4所示： 图4 四.仿真： 标签页 simulation.选择下面图标。修改仿真时间参数如下图5：
图5 完成仿真观察机构运动状况。图6为第0.97S时的仿真图像 图6 图7为第2.91S时的仿真图像 图7 图8为第8.24S时的仿真图像
转子动力学大作业
转子动力学大作业 学院： 姓名： 班级： 学号：
目录 一、作业题目介绍 二、转子动力学理论简介 三、参数的选择和计算 四、Ansys分析临固有频率和临界转速 五、失稳转速影响因素及计算
一、大作业题目 1、 计算临界转速； 2、 圆轴承，长颈比为0.8，油膜间隙2‰ 3、 计算失稳转速 注：转子两端各一个轴承，支点在左右两端。 二、转子动力学理论知识 由于制造中的误差，转子各微段的质心一般对回转轴线有微小偏离。因此，当转子转动时，会出现横向干扰，在某些转速下还会引起系统强烈振动，出现这种情况时的转速就是临界转速。为保证系统正常工作或避免系统因振动而损坏，转动系统的转子工作转速应尽可能避开临界转速，若无法避开，则应采取特殊防振措施。这也是研究临界转速的意义。临界转速和转子不旋转时横向振动的固有频率相同，也就是说，临界转速与转子的弹性和质量分布当圆盘不装在两支撑的中点而偏于一边时，转轴变形后，圆盘的转轴线与两支点A 和B 的连线有一夹角ψ。设圆盘的自转角速度Ω，极转动惯量为p J ，则圆盘对质心o '的动量矩为p H J =Ω。它与轴线AB 的夹角也应该是ψ，见图1。当转轴有自然振动时，设其频率为n ω，则圆盘中心o '与轴线AB 所构成的平面绕AB 轴有进动角速度n ω。由于进动，圆盘的动量矩H 将不断改变方向。因此有惯性力矩 ()g n n p n M H H J ωωω=-?=?=Ω? 方向与平面o AB '垂直，大小为 sin g p n M J ωψ=Ω 转子结构尺寸示意图
这一惯性力矩称为陀螺力矩或回转力矩。因夹角ψ较小，sin ψψ≈，上式可写作 g p n M J ωψ=Ω。 这一力矩与ψ成正比，相当于弹性力矩。在正进动(0/2ψπ<<)的情况下，它使转轴的变形减小，因而提高了转轴的弹性刚度，即提高了转子的临界角速度。在反进动（/2πψπ<<）的情况下，这力矩使转轴的变形增大，从而降低了转轴的弹性刚度，即降低了转子的临界角速度。通过分析，可知道陀螺力矩对转子临界转速的影响：正进动时，它提高了临界转速；反进动时，它降低了临界转速。 图 1 在大多数情况下，轴承对于转子的动力特性有很明显的影响，轴承往往是阻尼的主要来源，因而控制着转子的响应。同时，轴承的刚度和阻尼又影响着转子的临界转速和稳定性。在深入研究转子动力学问题时，因而必须考虑到轴承的作用。对于一个确定的轴承，当润滑油粘度及进油压已给定时，轴颈中心1o 的静平衡位置e 、?决定于轴颈转速Ω和静载荷W 。当载荷W 的大小或者轴颈转速Ω变化时，1o 位置也相应地变化。当铅垂载荷W 大小变化时，轴颈中心的移动在大多数情况下，并非沿铅垂方向，也即位移并不沿着载荷作用的方向。这正是油膜不同于一般机械元件的一个特点。 记x F 、y F 为油膜力在x 、y 方向的分量。我们定义油膜刚度系数为单位位移所引起的油膜力增量，即 x xx F k x ?=?，0 y xy F k y ?= ?，0 y yx F k x ?= ?，0 y yy F k y ?= ? 定义油膜阻尼系数为单位速度所引起的油膜力增量，即 x xx F c x ?= ? ，0 x xy F c y ?= ? ，0 y yx F c x ?= ? ，0 y yy F c y ?= ? 式中各系数的第一个下标代表力的方向，第二个下标代表位移或速度的方向。油膜刚度系数和阻尼系数统称为油膜动力特性系数。其中xy k ，yx k 和xy c ，yx c 分别称为交叉刚度系数和交叉阻尼系数，它们表示油膜力在两个相互垂直方向的耦合作用，交叉动力系数的大小和正
机械系统动力学
《机械系统动力学》 机械系统动力学中分析中的 仿真前沿 学院：机械工程学院 专业：机制一班 姓名：董正凯 学号：S12080201006
摘要 计算机及其相应技术的发展为建立机械系统仿真提供了一个有效的手段，机械系统动力学中的许多难题均可以采用仿真技术来解决，本文主要讲述了目前在机械系统动力学的分析中仿真技术主要的研究重点及其研究中主要存在的问题。 关键词：机械系统动力学仿真系统建模
机械系统动力学中分析中的仿真前沿 机械专业既是一个传统的专业，又是一个不断融合新技术、不断创新的专业。随着科技的发展，计算机仿真技术越来越广泛地应用在各个领域。基于多体系统动力学的机械系统动力学分析与仿真技术，从二十世纪七十年代开始吸引了众多研究者，已解决了自动化建模和求解问题的基础理论问题，并于八十年代形成了一系列商业化软件，到了九十年代，机械系统动力学分析与仿真技术更已能成熟应用于工业界。 目前的研究重点表现在以下几个方面： （1）柔性多体系统动力学的建模理论 多刚体系统的建模理论已经成熟，目前柔性多体系统的建模成了一个研究热点，柔性多体系统动力学由于本身既存在大范围的刚体运动又存在弹性变形运动，因而其与有限元分析方法及多刚体力学分析方法有密切关系。事实上，绝对的刚体运动不存在，绝对的弹性动力学问题在工程实际中也少见，实际工程问题严格说都是柔性多体动力学问题，只不过为了问题的简化容易求解，不得不化简为多刚体动力学问题、结构动力学问题来处理。然而这给使用者带来了不便，同一个问题必须利用两种分析方法处理。大多商用软件系统采用的浮动标架法对处理小变形部件的柔性系统较为有效，对包含大变形部件的柔体多体系统会产生较大仿真分析误差甚至完全错误的仿真结论。最近提出的绝对节点坐标方法，是对有限元技术的拓展和较大创新，在常规有限元中梁单元、板壳单元采用节点微小转动作为节点坐标，因而不能精确描述刚体运动。绝对节点坐标法则采用节点位移和节点斜率作为节点坐标，其形函数可以描述任意刚体位移。利用这种方法梁和板壳可以看作是等参单元，系统的质量阵为一常数阵，然而其刚度阵为强非线性阵，这与浮动标架法有截然不同的区别。这种方法已成功应用于手术线的大变形仿真中。寻求有限元分析与多刚体力学的统一近年来成为多体动力学分析的一个研究热点，绝对节点坐标法在这方面有极大的潜力，可以说绝对节点坐标法是柔性多体力学发展的一个重要进展。另外，各种柔性多体的分析方法之间是否存在某种互推关系也引起了人们的注意，如两个主要分析方法：浮动标架法、绝对节点坐标法之间是否可以互推？这些都具有重大理论意义。 另外柔性多体系统动力学中由于大范围的刚体运动与弹性变形运动相互耦合，采用浮动标架法时，即便是小变形问题，由于处于高速旋转仍会产生动力刚化现象。如果仅仅采用小变形理论，将产生错误的结论，必须计及动力刚化效应。动力刚化现象已成为柔性多体动力学的一个重要研究方面。如何利用简单的补偿方法来考虑动力刚化是问题的关键。 柔性多体系统动力学中关于柔性体的离散化表达存在三种形式：基于有限元分析的模态表达，基于试验模态分析的模态表达和基于有限元节点坐标的有限元列式。有限元列式由于大大地增加了系统的求解规模使其应用受到限制，因而一般采用模态分析方法，对模态进行模态截断、模态综合，从而缩减系统的求解规模。为了保证求解精度，同时又能提高求解速度如何进行模态截断、模态综合就成了一个关键问题。再者如何充分利用试验模态分析的结果也是一个关键性研究课题，这一方面的研究还不够深入。 柔性多体系统动力学可以计算出每一时刻的弹性位移，通过计算应变可计算计算出应力。由于一般的多柔体分析程序不具备有限元分析功能，因而柔性体的应力分析都是由有限元程序处理。由于可以计算出每个柔性体的应力的变化历
现代控制大作业
现代控制理论大作业 桥式吊车工作过程自动调节在状态空间分析中的 设计与计算 专业： 姓名： 学号： 日期：
桥式吊车工作过程自动调节在状态空间分析中的设计与计算 1. 小车-吊钩(机械)系统动力学方程 在不计小车与桥架(轨道)之间摩擦力的情况下，小车在水平(s 轴)方向上有如下作用力平衡方程： 对于吊钩，则在水平与垂直(z 轴)方向上可分别得到如下作用力平衡方程： 与上述3个力平衡方程相对应，在假定绳索长度l 不变条件下，还可得如下两个运动学方程： 为消去式(1)~(3)中的中间变量:绳索拉力p,可将式(1)、(2)两边相加得： 将(2)、(3)两边分别乘以cos θ和(-sin θ)后再相加得： 式(6)、(7)中不再含参数p ，进一步由(4)、(5)又可分别得： 最后，把式(8)、(9)代入式(6)、(7)后可分别得： 至此，小车-吊钩（机械）系统可用式(10)、(11)两个二阶非线性微分方程进行描述，显然这是一个四阶动力学系统。 解析求解式(10)、(11)是困难的，也没有必要，可以从工程角度（或通过非线性方程线性化）进行化简。 从调节（控制）技术角度讲，常可采用某种调节（控制）手段，如全状态反馈闭环调节（控制），使θ 角的变化（相对于稳态值的偏差量）控制在一个很小的范围内，例如，在此前提下，就可以进行如下近似处理，即令: 由此，式(10)、(11)可分别写为： )1( sin θp F s m A A A += )3( cos )2( sin θθp g m z m p s m B B B B B -=-= )5( cos )4( sin θθl z l s s B A B =+=)6( A B B A A F s m s m =+ )7( sin sin cos θθθg m z m s m B B B B B -=- )9( )sin cos ()8( )cos sin (22θθθθθθθθ --=+-+=l z l s s B A B )10( sin cos )(2A B B A B A F l m l m s m m =-++θθθθ (11) 0sin cos =++θθθg l s A 0sin ,1cos ,sin 2≈≈≈θθ θθθ
ANSYS 中的转子动力学计算
ANSYS 中的转子动力学计算 雷先华 安世亚太 转子动力学是固体力学的一个重要分支，它主要研究旋转机械的“转子－支承”系统在旋转状态下的振动、平衡和稳定性问题，其主要研究内容有几个方面：临界转速、动力响应、稳定性、动平衡技术和支承设计。在旋转机械研究设计中，转子动力学的性能分析是极其重要的一个方面。 传统的转子动力学分析采用传递矩阵方法进行，由于将大量的结构信息简化为极为简单的集中质量—梁模型，不能确保模型的完整性和分析的准确度；而有限元在处理转子动力学问题时，可以很好地兼顾模型的完整性和计算的效率，但多年来转子的“陀螺效应”一直是制约转子动力学有限元分析的“瓶颈”问题。ANSYS 很好地解决了动力特性分析中“陀螺效应”影响的问题，而且陀螺效应的考虑不受计算模型上的限制，使得转子动力学有限元分析变得简单高效。 本文对ANSYS 的转子动力学计算功能进行简要介绍。 1 ANSYS 转子动力学的理论基础 ANSYS 转子动力学分析中，两种参考坐标系可供选择：静止坐标系和旋转坐标系。空间点P 在静止坐标系（其原点在O ′）下的位置矢量为r ′，在旋转坐标系（其原点在O ）下的位置矢量为r 。 在静止坐标系下转子的动力方程为： [][][]{}{}([]){}{}M u C C u K u F gyr +++=&&& []C gyr 为陀螺效应矩阵。 在旋转坐标系下转子的动力方程为： [][][]{}{}([]){}([]){}M u C C u K K u F cor spin r r r +++?=&&& []C cor 为哥氏效应矩阵，[]K spin 为旋转软化效应刚度矩阵。 2 ANSYS 转子动力学的计算功能和新技术 ANSYS 转子动力学计算包含如下功能： – 无阻尼临界转速分析 – 不平衡响应分析 – 阻尼特征值分析
机械动力学大作业
机械动力学大作业 ——基于ADAMS的曲柄滑块机构动力学分析及仿真 学号： 专业：机械工程 学生姓名： 任课教师：杨恩霞教授 2013年10月21日
一、题目要求 采用ADAMS软件或Matlab/Simulink 环境，建立简单机械系统的动力学模型，借助软件进行求解计算和结果分析。建立单自由度杆机构（有无滑块均可）动力学模型，由静止启动，选择一固定驱动力矩，绘制原动件在一周内的运动关系线图，具体机构及参数自拟。 二、所选题目 我选择的机构为曲柄滑块机构，如图2-1所示，采用ADAMS 软件，建立简单机械系统的动力学模型，借助软件进行求解计算和结果分析。 具体参数如下： AB杆（曲柄）：L=0.5m W=0.2 m D=0.1 m m=20kg BC杆（连杆）：L= 2m W= 0.2 m D= 0.1 m m =40kg 滑块： L=0.6 m W=0.6 m D=0.6 m m =30kg 驱动力矩： M=200 N.m 阻力： F=100N 图2-1 三、建立模型 1.启动ADAMS 在欢迎对话框，选择New Model 项；在模型名称栏输入qubing；重力设置选择Earth Normal参数；单位设置选择MKS系统（m，kg，N，s，deg，）如图3-1
图3-1 2检查和设置建模基本环境 1）检查默认单位系统在Settings菜单中选择Units 命令，显示单位设置对话框，当前的设置应该为MKS系统。 2）设置工作栅格 ①在Settings菜单，选择Working Grid命令，显示设置工作栅格对话框； ②设置Size X=5.0，Size Y=5.0，Spacing X=0.1，Spacing Y=0.1 ③选择OK按钮。如图：
转子轴承系统动力学分析系统的设计与实现
转子轴承系统动力学分析系统的设计与实现 朱爱斌1，张锁怀2，丘大谋1，谢友柏1 (1.西安交通大学 润滑理论及轴承研究所,陕西西安 710049; 2. 上海应用技术学院,上海 200235) 摘 要: 分析了如何基于Matlab和VB开发齿轮啮合的转子轴承系统动力学分析系统的问题，介绍了系统的总体设计和具体实现途径，提出将Matlab和VB的三种集成方法混合应用，并通过实例说明系统的使用方法和计算分析内容。该系统能够有效缩短齿轮啮合的转子轴承系统的设计开发周期，优化系统的性能。关键词:转子轴承；齿轮；动力学分析；Matlab 中图分类号: TH12；TP312 文献标识码:A Design and Realization of Rotor-Bearing System's Dynamic Characteristics Analyzing System ZHU aibin1 ZHANG suohuai2 QIU damou1 XIE youbai1 (Theory of Lubrication and Bearing Institute, Xi'an Jiaotong University, Xi’an 710049, China) Abstract : Issue of how to develop dynamic characteristics analyzing system of geared rotor-bearing system with Matlab and Visual Basic was analyzed, framework design and realized approach was introduced, and method of mixed application of three integration ways between Maltab and VB was proposed. A case was given to show the computing and analyzing process. The analyzing system can efficiently short the design and development time of geared rotor-bearing system, and optimize the performance of geared rotor-bearing system. Key words : rotor bearing; gear; dynamic characteristics analyzing; matlab 齿轮耦合的转子轴承系统即多个转子－轴承 系统通过齿轮耦合联系在一起[1][2]。这种系统既保留了单个转子－轴承系统的某些动力学特性，又具有齿轮传动所引起的一些新特性。某一转子－轴承系统的动力学性能的改变，通过齿轮的耦合作用，必将影响另一转子－轴承系统的动力学性能；横向振动通过齿轮传递后，将引起转子产生扭转振动，也就是说，弯曲振动和扭转振动将同时发生，即发生弯扭耦合振动[3]；齿轮参数的改变，必将导致整个系统的动力学性能发生变化，这是该系统所独有的特性。 具有齿轮啮合的转子轴承系统在风机、压缩机、增速器等机器中广泛存在，由于齿轮的啮合作用，使原本相互独立的多个转子轴承系统联接在一起，从而使各转子轴承系统的动力特性相互影响，整个系统的动力特性与单个子系统的动力特性大不一样[4]。在齿轮耦合的转子轴承系统的研究基础上，基于Matlab和VB开发了齿轮啮合的转子轴承系统动力学分析系统，可以用于压缩机、风机等流体机械及增速器、减速器等具有齿轮传动的平行轴系的转子系统动力学分析。分析内容包括稳定性、临界转速、强迫振动响应、系统特征值及振型的计算和分析。同时本系统也能够完成转子轴承系统中任意单根转子的动力学分析。1 总体框架设计 1.1 系统设计原则 系统是面向具有转子轴承系统动力学一般知识的企业或者科研院所用户而开发的，基本的设计原则包括： (1) 建立考虑齿轮啮合因素的平行轴系的转子轴承系统的数学模型，使其计算结果能够与实际情况的误差较小； (2) 提供简单、合理和方便的使用界面，适应不同使用水平的用户； (3) 提供包括数据，图形，XML文档等多种形式的丰富的参数表示形式，给用户直观，丰富的信息； (4) 结合Matlab的数据处理，矩阵计算和图形1显示的强大功能和VB在图形用户界面开发方面的优势； 1.2 系统总体框架 收稿日期:2004 - 09 - 14 基金项目:博士学科点专项科研基金(20030698005，20050698016) https://www.360docs.net/doc/5f15356886.html,