




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、固有频率在ADAMS/Linear和ADAMS/Vibration中的理解在ADAMS中,固有频率是通过本征向量计算的,为了更好的理解计算结果中各个参数的意义,解决仿真中常见的问题,在这里理论联合实际对一些基本知识在ADAMS中的应用做一基本论述。在此,不涉及ADAMS/Linear的扩展命令,所有的线性化命令实际都是在图形界面操作所得的。对于单自由度系统,如经典的弹簧质量阻尼系统,质量m的运动方程有:或 (1)这里x为质量m的位移,k为弹簧刚度系数,c为阻尼系数。根据无阻尼固有圆频率和阻尼比的定义重写等式(1): (2)这里:无阻尼固有圆频率(Undamped Natural Frequen
2、cy) (3)阻尼比(Damping Ratio) (4)可以看出,无阻尼固有圆频率只是弹簧刚度k和质量m的函数,与阻尼值无关。ADAMS/Linear实际上计算无阻尼固有圆频率的方法有所不同,它使用拉普拉斯(Laplace)在仿真运行点对模型变换为线性矩阵,再通过本征值向量(Eigenvalues)计算系统的固有圆频率和阻尼比,但计算结果与上述计算是等效的。一般,本征值由实部(Real part)和虚部(Imaginary part)两部分组成:,因此,方程式(2)可以写为: (5)本征值由下式决定:当阻尼比1, (6)当阻尼比1, (7)令:;。当系统阻尼比当1时,ADAMS/Linear
3、使用下式计算无阻尼固有圆频率与阻尼比: (8)即:,或 (9)复数平面本征值示意图显然,当实部0时,系统阻尼比0;当虚部0时,1。所以当阻尼比大于或等于1时为过阻尼或临界阻尼,此时本征值不能被计算即无振动解,阻尼值的结果均报告为1。现在,我们以一个简单的单自由度系统分别使用理论公式和ADAMS/Linear计算无阻尼固有频率和阻尼比。为简便起见,设定运动质量m1kg,弹簧刚度系数k1000N/m,阻尼系数c20N-s/m。首先,将上述参数分别代入等式(3)、(4)和(7)计算理论值:无阻尼固有圆频率(Hz)阻尼比本征值然后,启动ADAMS/View并建立一个单自由度模型,注意,为了获得振动质量
4、m的单自由度,应对其施加垂向约束。图2 单自由度模型示意图质量m弹簧与阻尼地面然后运行仿真:Static - Linear;得到的本征值信息与计算结论如图3,与理论计算结果对比得知是完全一致的。图 3 单自由度本征值计算结果表在图3中,计算结果包括4列:无阻尼固有圆频率(UNDAMPED NATURAL FREQUENCY)与阻尼比(DAMPING RATIO);其中本征值:实部(REAL)和虚部(IMAGINARY)。如果修改阻尼系数为10 N-sec/m,再次在静平衡仿真点计算本征值,结果如图4:修改后的单自由度系统本征值结果表可以得到阻尼系数的改变会影响到阻尼比、实部值和虚部值,但不会改
5、变系统的固有频率,同样与理论推导的结果同样完全一致。现在,将上述模型中的移动副删除,去处质量m的强制约束使其拥有6个自由度,执行同样的线性化命令后得到的本征值如图6所示。图5 模型核查信息图6 单物体6自由度本征值信息表其中11号(1阶非刚体)模态(有实部和虚部2部分)为质量m沿垂直轴位移的振动特性,数值的意义见前述,它是我们在工程应用中关心的模态。第711为过阻尼态,第16为无阻尼态(阻尼110-10),对于本征值虚部为0的模态,我们还会收到警告消息,如下:图 7 本征值虚部无解警告对于6自由度系统,总共有3对模态;实部不为0的第7、8、9、10组成2对模态加上第11号模态共计3对,即:(模
6、态对)2自由度数目。对于其他模型,可以用同样的方式预估模态数目。还有一种情景是当质量m的运动由数学方程(一阶微分方程、单输入,单输出传递函数、线性状态方程、通用状态方程)支配时,也会出现0值的本征值。共振频率和隔振区间的计算工程中经常发生由变质的惯性力激励产生的受迫振动,例如由于地基振动引起结构物得受迫振动,或由于转子偏心引起的受迫振动。这类受迫振动的特点是激励惯性力的振幅与频率的平方成正比。以地基振动为例,设安装质量弹簧系统的基座沿y轴方向作振幅为B、频率为的简谐振动,振动规律为则在物体上产生简谐变化的惯性力将物体相对基座的相对位移记作x1,则动力学方程可写作令相对运动的振幅为A1,振幅扩大
7、因子为1=A1/B,得到参数1,1与激励频率之间的关系为其中s=/0。由于,当激励力频率远小于固有频率是,相对运动的振幅接近于零,相位接近相同;激励频率远大于固有频率时,相对运动振幅接近于基座运动的振幅,但相位正好相反。激励频率接近固有频率时,也存在振幅急剧增大的共振现象。若将质点相对惯性坐标系的绝对位移作为响应y,则y等于相对位移y1与基座牵连位移yf之和,由于令A为绝对振动振幅,得到导出绝对运动的振幅放大因子=A/B和相位差为将地基的振动与机器设备隔离以避免将振动传至设备,隔振的效果用被动隔振系数表示,定义为隔振后设备的振幅与隔振前设备的振幅之比,导出的被动隔振系数为:再求的最大值,即求的
8、最大值:令,将其转换成则令所以当时,产生最大共振振幅。ADAMS对隔振系统的振动分析针对上述的设定好参数的弹簧阻尼质量系统的ADAMS模型,若将其地面的固定的约束解除,再给地面一个竖直方向的简谐振动位移,此时该系统就是隔振系统。地面振动通过弹簧和阻尼传递到物体上,已知m=1kg,k=1000N/m,c=10N-S/m。根据理论算出的固有频率为0=5.032921,阻尼比为=0.1581139,进一步算得当地面简谐振动频率7.117625时,隔振系统具有消振效果;当=4.916963时,物体产生最大振幅的振动。在ADAMS里,解除固定约束,给地面施加一个竖直方向的位移,并定义该位移为三角函数位移
9、(右击MOTION_7,选择定义选项modify,在弹出的定义窗口空,最初按右图图所示编辑,单机OK键)。运行仿真后,地面约束改成三角函数位移时,ADAMS计算出系统的固有频率和阻尼比保持不变。改变地面位移的参数方程,观察不同的隔振效果。(给质量块单自由度约束,添加地面位移驱动,再给地面单自由度约束)特别地,当把三角函数的方程改为10*sin(44.72*time),即频率为7.117625时,运行仿真后:分别右击地面和质量块,选择Measure观察各自的振动情况,在Measure的编辑栏里,选择Y轴的位移作为输出量,单机OK键观察其对比效果。很容易看出,此条件下,当物块的振动达到稳定时,其振
10、幅与地面振幅基本一致,这也验证了前面所计算出的理论值的正确性。 在正确性得以验证的基础上,再将竖直位移函数编辑为10*sin(30.89*time),通过运行仿真,可以进一步的观察到共振情况下,质量快的大幅振动。隔振效果曲线的绘制根据理论表达式,可以粗略的绘出不同阻尼所对应的幅频曲线:1212120.750.50.25振幅12120.750.50.25振幅其中,横坐标为/0。一般的/02时,系统有放大作用;/01时,系统共振,传递比到极大;2/03时,作用有限;3/06时,隔振能力低;6/010时,强;ADAMS专门开发了振动仿真分析模块ADAMS/Vibration,在振动仿真分析模块中可以
11、进行类似真实的实验一样对系统进行多方面的测试。它主要是在其模块建立的基础上,在输入位置定义激励,在输出位置处计算频响函数,通过对输出响应的位置分析,来决定系统振动性能方面的特性。继续研究隔振系统,可以通过这一技术的运用,绘出能够反映地面振动频率改变时,隔着系统的隔振效果。同样是先解除原有弹簧阻尼质量系统的地面固定;创建一个输入通道和激励单击菜单【Vibration】【Build】【Input Channel】【New】,在弹出的激励编辑栏中选择位移方式激励,Y方向,幅度和初相位分别为10和0,选择地面模块的中心点作为输入点。创建两个对比输出通道和激励单击菜单【Vibration】【Build】【Output Channel】【New】,在输出响应的编辑栏里选择Y轴的位移输出方式,第一个输出点选择地面模块中点,第二个输出点选择质量块中心点。仿真分析。单击菜单【Simulate】【ADAMS/Vibration】【Vibration Analysis】后,弹出运行振动分析对话框。然后选择Assembly。在Inpu_y下导入所创建的输入激励,在Output Channel下导入所创建的两个输出通道;频率Begin和End分别设置为0.1和100,步数为1000。设置
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论