matlab单自由度系统,单自由度系统,响应的matlab程序.doc

单自由度系统,响应的matlab程序.doc

1FUNCTIONFHANNING_IMPT,TC,A2FZEROSSIZET3FTTCA/21COS2PITTTC/TC4END下面是计算单自由度系统响应的MATLAB程序，计算传递函数，画实频、虚频、幅频、相位、导纳图PLAINVIEWPLAINCOPY1M1002K10003C10045NUM16DENMCK7SYSTFNUM,DEN8DT0000019FS1/DT采样频率HZ100HZ实际并不需要这么高的采样频率，但是如果采样时间太小，HANNING脉冲不完整10为了得到准确的响应DT一定要小，否则做出的相位可能不对1112T0DT20013TC000114A1015UHANNING_IMPT,TC,A16YLSIMSYS,U,T17YIMPULSESYS,T18YY 1920NLENGTHU21FYFFTY22FUFFTU23FTFY/FU24F0N1FS/N25FT_RREALFT26FT_IIMAGFT2728PARTF302930FIGURE NAME , 实频 31PLOTFPART,FT_RPART32FIGURE NAME , 虚频 33PLOTFPART,FT_IPART34FIGURE NAME , 幅值 35PLOTFPART,ABSFTPART36FIGURE NAME , 相位 37PLOTFPART,PHASEFTPART180/PI38FIGURE NAME , 导纳圆 39PLOTFT_RPART,FT_IPART, 40AXISEQUAL41HOLDON42XK000011654300002857448297E5450000486946000068247XK1000080174800003082490000285750000119451000145152YK100028175300030955400031855500028695600026335758YK000314959000318560000317961000316562000312163K564A65SUMXK2SUMXKYKSUMXK66SUMXKYKSUMYK2SUMYK67SUMXKSUMYKK68B69SUMXK3XKYK270SUMXK2YKYK371SUMXK2YK272RLTAB73X0RLT10574Y0RLT20575RSQRTRLT12/4RLT22/4RLT376FAI00012PI77XX0RCOSFAI78YY0RSINFAI79PLOTX,Y, R 实验模态分析非数学公式的简单概述之二分类模态空间译文|标签锤击法激振器窗函数曲线拟合201107022300阅读3040评论0为何只需获得频响函数矩阵的一行或一列理解从可能得到频响函数矩阵的不同元素中得到模态振型对我们来说是非常重要的。在这不涉及数学层面的知识，让我们来讨论这个问题。首先考虑频响函数矩阵的第三行，并且只关注第1阶模态，留意频响函数虚部的峰值振幅，很容易就能得出结构的第1阶模态振型，如图8A所示。因此，从测量数据中提取模态振型似乎相当直观。一种快速但又粗略的方法就是在不同的测点处仅仅测量频响函数虚部的峰值振幅。图8A从频响函数矩阵第三行得到的1阶模态接着考虑频响函数矩阵的第二行，并且只考察第1阶模态，如图8B所示。留意频响函数虚部的峰值振幅，从这一行也易于得到第1阶模态振型。图8B从频响函数矩阵第二行得到的1阶模态我们同样可以从频响函数矩阵的第一行得到这一阶模态振型。这是理论所表达的一种简单示意性描述。我们可以使用频响函数任一行得到系统的模态振型。故很显然，这些测量包含有与系统模态振型相关的信息。现在再考虑频响函数矩的阵第三行，并且只考察第2阶模态，如图8C所示。还是留意频响函数的虚部的峰值振幅，很容易得到第2阶模态振型。图8C从频响函数矩阵第三行得到2阶模态而观察频响函数矩阵的第二行，并且只考察第2阶模态。此时会有点奇怪，因为这一行没有第2阶模态可用的幅值，如图8D所示。这是我不希望发生的，但是如果我们考察第2阶模态振型，那么很快就会发现位置2是第2阶模态的节点。此时参考点位于模态节点上。图8D从频响函数矩阵第二行得到2阶模态这就指明了模态分析和实验测量中一个非常重要的方面参考点不能位于某阶模态的节点上，否则该阶模态在频响函数中将不可见，并且得不到该阶模态。在这我们仅用了3个测点去描述该悬臂梁的模态。如果我们增加更多的输入输出测点，就能得到更光顺的模态振型，如图9所示。图9显示了15个频响函数，其中前面讨论的3个测点的频响函数高亮突出显示。显示的15个频响函数用瀑布图式样绘出。利用这种方式绘图，通过频响函数的虚部峰值连线能更容易确定模态振型。目前为止，我们所讨论的测量是从锤击法测试中得到的，如果我们使用激振器测试，那么测量的频响函数会是什么样的呢图9瀑布图显示悬臂梁频响函数锤击法测试和激振器测试有什么不同之处从理论角度看，频响函数是由激振器测试得到还是由锤击法测试得到，并没有什么区别。图10A和10B给出了由锤击法测试和激振器测试得到的频响函数。锤击法测试通常测量频响函数矩阵中的一行，而激振器测试通常测量频响函数矩阵中的一列。因为描述系统的频响函数矩阵是对称的方阵，故互易性是成立的。例如，对于上面已讨论的情况，频响函数矩阵的第三行和第三列是完全相同的。理论上讲，激振器测试和锤击法测试两者没有差异，但那仅仅是理论观点。假如我可以对结构施加一个纯外力，外力与结构二者之间没有任何相互作用，并且用一个无质量的传感器测量响应，要求该传感器对结构没有任何影响，那么上面所讲的是正确的。但是事实并非如此，结果又将怎样呢现在我们从现实角度出发，考虑实际测试中存在的不同之处。模态测试过程中，关键在于激振器和响应传感器通常对结构确实有影响。需要注意最主要的一点是被测结构已不是你想得到模态参数的那个结构。因为在结构上已附加了与数据采集过程有关的东西结构支承条件、安装的传感器的质量、激振器推力杆的刚度影响等等。因此虽然理论告诉我们，锤击法测试和激振器测试不存在任何差异，但现实中却因数据采集方面导致二者存在差异。激振器测试过程中，最明显的差异是由移动加速度计引起的。加速度计的质量相对于结构的总质量可能非常小，但是它的质量相对于结构不同部分的有效质量可能又非常大。特别是多通道测试系统，这个问题更加突出，为了获得所有频响函数，有多个加速度计在结构上移动。特别是轻质结构，这个问题尤为突出。纠正此问题的方法之一是在结构上安装所有的加速度计，即使一次只用到少数几个加速度计。另一个方法是在非测量位置上安装与加速度计质量相等的质量哑元，这将能消除移动加速度计带来的影响。图10A移动力锤测试过程图10B移动响应传感器的测试过程另一个差异在于激振器推力杆带来的影响。本质上，结构的模态受激振器附属装置的质量和刚度的影响。虽然我们试图将这部分影响减少到最低程度，但是它们仍然是存在的。激振器推力杆的作用是分离激振器对结构的影响，然而，多数结构，激振器附属装置的影响仍然很大。因为锤击法测试不存在这些问题，所以得到的结果不同于激振器测试得到的结果。所以虽然理论上讲激振器测试和锤击法测试二者不存在差异，但一些非常基本的现实情况却会引起一些差异。为了计算频响