《机械结构实验模态分析》实验报告
开课实验室：汽车结构实验室 2019年月日 学院
姓名
成绩
课程 名称 机械结构实验模态分析 实验项目 名 称
机械结构实验模态分析 指导教师
教师评语
教师签名：
年 月 日
机械结构实验模态分析实验报告 一、实验目的和意义
模态分析技术是近年来在国内外得到迅速发展的一门新兴科学技术，广泛应用于航空、航天、机械制造、建筑、汽车等许多领域，在识别系统的动力学参数、动态优化设计、设备故障诊断等许多方面发挥了日益重要的作用。
本实验采用CCDS-1模态分析微机系统，对图1所示的框架结构进行分析。通过该实验达到如下目的：
212019
1817
16
1514
13121110
987
6
5
4
3
222120
20
202090
9090
90
90909090113
113
113
113
113
113
115
115
115
115 图1 框架结构图
详细了解CCDAS-1模态分析微机系统，并熟练掌握使用本系统的全过程，包括
了解测量点和激振点的选择。
了解模态分析实验采用的仪器，实验的连接、安装和调整。
1、 激励振时各测点力信号和响应信号的测量及利用这些测量信号求取传递函数，并分析影响传递
函数精度的因素。
2、 SSDAS-1系统由各测点识别出系统的模态参数的步骤。
3、 动画显示。
4、 灵敏度分析及含义。
通过CCDAS-1模态分析的全部过程及有关学习，能祥述实验模态的一般步骤。 通过实验和分析，大大提高综合分析能力和动手能力。
CCDAS-1系统模态分析的优缺点讨论并提出改进实验的意见。 二、测试及数据处理框图
加速度传感器
力传感器
脉冲锤
四个点由橡胶绳悬挂
1724 打印机
IBM PC 微型计算机 含AD板 CCMAS-1模态分析软件
双通道低 通滤波器
电荷放大器
电荷放大器
图2 测量及数据处理系统框图
三、实验模态分析的基本原理
对于一个机构系统，其动态特性可用系统的固有频率、阻尼和振型来描述，与模态质量和模态刚度一起通称为机械系统的模态参数。模态参数既可以用有限元的方法对结构进行简化得到，也可以通过激振实验对采集的振动数据进行处理识别得到。通过实验数据求取模态参数的方法就是实验模态分析。只要保证测试仪器的精度、实验条件和数据分析处理的精度就能获得高质量的模态参数。
一个线性系统，若在某一点j 施加激振力j F ，系统各点的振动响应为i X 1,2,...,i n =，系统任意两点的传递函数ij h 之间的关系可用矩阵表示如下：
11112122122212()...
0()...()...()...0n n j n n n nn x h h h x h h h F x h h h ωωωω??????
????????????
=???
???????????????
??????M M M O M (1-1)
可记为：{}{}[]X H F =
[]H 称为传递函数矩阵。其中的任意元素ij h 可以通过激振实验得到
()
()
i ij j X h F ωω=
()i X ω，()j F ω分别表示响应i X 与激振力j F 的傅立叶变换。
测量方法是给系统施加一有限带宽频率的激振力(冲击也是一有限带宽激振力)，同时测量系统的响应，将力和响应信号进行滤波，A/D 转换并离散采样，进行双通道FFT 变换，计算出激振力j F 与响应i X 之间的传递函数ij h 。
对测量的传递函数进行曲线拟和得到模态参数，一个多自由度系统曲线拟和传递函数的解析式为：*
*
1
()[]n
ijk
ijk
ij k k
k r r h S S P
S P ==
-
--∑ (1-3)
式中：S j ω=，k k k p j σω=-+
k σ—第K 阶的模态阻尼 k ω—第K 阶的固有频率 ijk r —第K 阶留数 *—表示复数共轭
由上式可知，各点传递函数的分母相同，都包含固有频率k ω和模态阻尼k σ，由任意点传递函数曲线拟和可以识别出系统的固有频率和阻尼，考虑到传递函数的测量和计算误差，将多测点传递函数相加合成一系列的传递函数，并由系统传递函数识别系统的固有频率k ω和模态阻尼k σ。
由式(1-3)可知，每一测点的传递函数曲线拟和得到的留数ijk r 是不同的，留数矩阵[]K R 与对应的固有频率振型的关系为：
1112
1111212122
2212221212...
.........[]......n n n n K n n nn n n n n r r r r r r R r r r φφφφφφφφφφφφφφφφφφ????
?????
???==????????????????M M O M M M O M (1-4)
[]K R 是对称矩阵，其中每一行，每一列都是线性相关的。其中任意元素为
ij i j r φφ= (1-5)
{}φ称为系统的模态向量，即固有振型。只要利用留数矩阵的某一行或某一列，先将对应角元
素包含的振型元素求出：
jj j j R φφ= j jj r φ= (1-6)
然后，即可求出任一点的振型元素：ijk i j r φφ=
1,...,1,1,... ij
i j
r i j j j j n φφ=
≠=-+ (1-7)
CCDAS-1模态分析微机系统要求采用某一点固定激振，逐点拾振，求出留数矩阵之一列的测
量方法，或在某一固定点拾振，求出留数矩阵之一行的方法。本实验采用后一种方法。
识别出系统的固有频率、阻尼和振型后，系统的模态阻尼和模态刚度由下式决定：
21
1 k k k
m σω=
- (1-8)
21 k k k K ωσ=- (1-9)
k m —模态质量
k K —模态刚度
k σ—模态阻尼比
至此整个系统的模态参数全部确定，其中最重要的参数是频率，阻尼和振型。振型{}k φ代表了再系统的第K 阶固有频率下，各测点位移振幅之间的比例关系。
设2
k k λω=，则系统对结构物理参数j R 的微分灵敏度为：
[][]
{}[]{}k k k k j j j
K M R R R λφλφ???=-??? (1-10) []K 、[]M 分别为系统的刚度矩阵和质量矩阵，j R 为物理参数，可为刚度和质量。振型的微
分灵敏度可表示为：
1
{}{}n
k kji i i j a R φφ=?=?∑ (1-11) 1[][]
{{}[]{}}[]
T kji i k i k i j j K M a R R φλφλλ??=
--?? (1-12)
灵敏度分析的意义在于，若你想改变系统的动态特性，对每阶振型来说，改变哪一点的质量
或刚度，其频率和振型的变化率最大；同样，若维持系统动态特性不变，改变哪一点的质量或刚度，其频率和振型的变化率最小。
四、用CCDAS-1系统对框架进行实验模态分析的内容和步骤 模态分析步骤：
1、连接仪器，并调整好各仪器的开关档位，传递函数频率分析范围为500Hz 以内，低通可取截止频率为300Hz 档。
2、选择测点，建立被测结构的几何模型。
3、测量信号的数据采集及双通道谱分析。
4、调用十个传递函数同时显示功能确定系统的前五阶固有频率。
5、求取系统多测点加权传递函数幅频特性并产生频率阻尼识别文件。
6、识别系统的频率阻尼产生频率阻尼数据文件。
7、识别系统振型，产生振型数据文件。
8、可以显示打印系统的频率阻尼表及振型表。
9、模态的动画显示观察各阶振型的特点，复模态和实模态的区别。 10、灵敏度分析
在进行框架模态分析中应注意以下问题：
1、结构测量点的选择及布置：模态分析布点的基本原则是结构的重要部分密布，次要部分稀疏，为了使活动振型的显示更直观，整个测点的连线构成的图形要保持测试系统的基本几何特点。
2、不变激振点或不变拾振点(以下简称为基本点)的选取：实验模态分析常采用单点激振，逐点拾振的方法或固定拾振点的方法。本实验采用后者。固定的激振点和固定的拾振点是非常重要的，我们应当选择能反映系统最多固有频率成分的点作为这种激振点或拾振点。用这点激振或拾振，共振频率的峰值较大，噪声对共振频率的影响较小，峰值附近相干函数值较高，另外基本点的选择不能在感兴趣的模态的振型节点处，否则这阶模态无法识别出来。
3、系统悬置问题：一般来说，测试系统的悬置方法有两种。一种是软连接，就是用橡皮绳、弹簧之类的弹性元件将测试系统悬置在某一刚性较好的支架上。测试系统自由状态下的另一种方式一般是用小功率的激振较重的系统，但在实际工作情况下，激振有时需要更大的功率，但测量的模
态参数实用意义更大一些。本实验采用软连接方法，用四根橡胶绳把测试框架结构悬置在支架上(如图2所示)。
五、实验结果及其分析
图3 第一阶扭转振型图4 第一阶弯曲振型
图5 第二阶扭转振型图6 第二阶弯曲振型
图7 第三阶扭转振型图8 第三阶弯曲振型
图9 第四阶扭转振型图10 第四阶弯曲振型
图11 传递函数图
《机械结构实验模态分析》实验报告 开课实验室：汽车结构实验室 2019年月日 学院 姓名 成绩 课程 名称 机械结构实验模态分析 实验项目 名 称 机械结构实验模态分析 指导教师 教师评语 教师签名： 年 月 日 机械结构实验模态分析实验报告 一、实验目的和意义 模态分析技术是近年来在国内外得到迅速发展的一门新兴科学技术，广泛应用于航空、航天、机械制造、建筑、汽车等许多领域，在识别系统的动力学参数、动态优化设计、设备故障诊断等许多方面发挥了日益重要的作用。 本实验采用CCDS-1模态分析微机系统，对图1所示的框架结构进行分析。通过该实验达到如下目的： 212019 1817 16 1514 13121110 987 6 5 4 3 222120 20 202090 9090 90 90909090113 113 113 113 113 113 115 115 115 115 图1 框架结构图 详细了解CCDAS-1模态分析微机系统，并熟练掌握使用本系统的全过程，包括 了解测量点和激振点的选择。 了解模态分析实验采用的仪器，实验的连接、安装和调整。 1、 激励振时各测点力信号和响应信号的测量及利用这些测量信号求取传递函数，并分析影响传递 函数精度的因素。 2、 SSDAS-1系统由各测点识别出系统的模态参数的步骤。 3、 动画显示。 4、 灵敏度分析及含义。 通过CCDAS-1模态分析的全部过程及有关学习，能祥述实验模态的一般步骤。 通过实验和分析，大大提高综合分析能力和动手能力。
CCDAS-1系统模态分析的优缺点讨论并提出改进实验的意见。 二、测试及数据处理框图 加速度传感器 力传感器 脉冲锤 四个点由橡胶绳悬挂 1724 打印机 IBM PC 微型计算机 含AD板 CCMAS-1模态分析软件 双通道低 通滤波器 电荷放大器 电荷放大器 图2 测量及数据处理系统框图 三、实验模态分析的基本原理 对于一个机构系统，其动态特性可用系统的固有频率、阻尼和振型来描述，与模态质量和模态刚度一起通称为机械系统的模态参数。模态参数既可以用有限元的方法对结构进行简化得到，也可以通过激振实验对采集的振动数据进行处理识别得到。通过实验数据求取模态参数的方法就是实验模态分析。只要保证测试仪器的精度、实验条件和数据分析处理的精度就能获得高质量的模态参数。 一个线性系统，若在某一点j 施加激振力j F ，系统各点的振动响应为i X 1,2,...,i n =，系统任意两点的传递函数ij h 之间的关系可用矩阵表示如下： 11112122122212()... 0()...()...()...0n n j n n n nn x h h h x h h h F x h h h ωωωω?????? ???????????? =??? ??????????????? ??????M M M O M (1-1) 可记为：{}{}[]X H F = []H 称为传递函数矩阵。其中的任意元素ij h 可以通过激振实验得到 () () i ij j X h F ωω= ()i X ω，()j F ω分别表示响应i X 与激振力j F 的傅立叶变换。 测量方法是给系统施加一有限带宽频率的激振力(冲击也是一有限带宽激振力)，同时测量系统的响应，将力和响应信号进行滤波，A/D 转换并离散采样，进行双通道FFT 变换，计算出激振力j F 与响应i X 之间的传递函数ij h 。 对测量的传递函数进行曲线拟和得到模态参数，一个多自由度系统曲线拟和传递函数的解析式为：* * 1 ()[]n ijk ijk ij k k k r r h S S P S P == - --∑ (1-3)
篇一：模态分析实验报告 模态分析实验报告 姓名：学号：任课教师：实验时间：指导老师：实验地点： 实验1传递函数的测量 一、实验内容 用锤击激振法测量传递函数。 二、实验目的 1) 掌握锤击激振法测量传递函数的方法； 2) 测量激励力和加速度响应的时间记录曲线、力的自功率谱和传递函数； 3) 分析传递函数的各种显示形式（实部、虚部、幅值、对数、相位）及相干函 数； 4) 比较原点传递函数和跨点传递函数的特征； 5) 考察激励点和响应点互换对传递函数的影响； 6) 比较不同材料的力锤锤帽对激励信号的影响； 三、实验仪器和测试系统 1、实验仪器 主要用到的实验仪器有：冲击力锤、加速度传感器，lms lms-scadas ⅲ测试系统，具体型号和参数见表1-1。 仪器名称 型号 序列号 3164 灵敏度 2.25 mv/n 100 mv/g 备注比利时 丹麦 b&k 数据采集和分析系统 lms-scadas ⅲ 2302-10 力锤 加速度传感器 表1-1 实验仪器 2 、测试系统 利用试验测量的激励信号(力锤激励信号)和响应的时间历程信号，运用数字 信号处理技术获得频率响应函数(frequency response function, frf)，得到系统的非参数模型。然后利用参数识别方法得到系统的模态参数。测试系统主要完成力锤激励信号及各点响应信号时间历程的同步采集，完成数字信号的处理和参数的识别。 测量分析系统的框图如图1-1所示。测量系统由振动加速度传感器、力锤和比利时lms公司scadas采集前端及modal impact测量分析软件组成。力锤及加速度传感器通过信号线与scadas采集前端相连，振动传感器及力锤为icp型传感器，需要scadas采集前端对其供电。scadas采集相应的信号和进行信号处理(如抗混滤波，a/d转换等)，所测信号通过电缆与电脑完成数据通讯。图1-1 测试分析系统框图 四、实验数据采集 1、振动测试实验台架 实验测量的是一段轴，在轴上安装了3个加速度传感器，如图1-2所示，轴由四根弹簧悬挂起来，使得整个测试统的频率很低，基本上不会影响到最终的测试结果。整个测试系统如下图所示：a1 a 测点2测点3测点4 图1-2 测试系统图
机械振动实验报告
《机械振动基础》实验报告（2015年春季学期）
专业机械设计制造及其自动化报告提交日期2015.05.07 哈尔滨工业大学
报告要求 1. 实验报告统一用该模板撰写，必须包含以下内容： (1) 实验名称 (2) 实验器材 (3) 实验原理 (4) 实验过程 (5) 实验结果及分析 (6) 认识体会、意见与建议等 2. 正文格式：四号字体，行距为1.25倍行距； 3. 用A4纸单面打印；左侧装订； 4. 报告需同时提交打印稿和电子文档进行存档，电子文档由班长收 齐，统 一发送至：Iiuyingxiang868@hit .edu .cn 5. 此页不得删除。 评语: 实验一报告正文 实验名称：机械振动的压电传感器测量及分析 教师签名: 年
二、实验器材 1、机械振动综台实验装置（压电悬臂梁）一套 2、激振器一套 3、加速度传感器一只 4、电荷放大器一台 5、信号发生器一台 6、示波器一台 7、电脑一台 & NI9215数据采集测试软件一套 9、NI9215数据采集卡一套 三、实验原理 信号发生器发出简谐振动信号，经过功率放大器放大，将简谐激励信号施加到电磁激振器上，电磁激振器振动杆以简谐振动激励安装在激振器上的压电悬臂梁。压电悬臂梁弯曲产生电流显示在示波器上，可以观测悬臂梁的振动情况；另一方面，加速度传感器安装在电磁激振器振动杆上，将加速度传感器与电荷放大器连接，将电荷放大器与数据采集系统连接，并将数据采集系统连接到计算机（PC机）上，操作NI9215数据采集测试软件，得到机械系统的振动响应变化曲线，可以观测电磁激振器的振动信号，并与信号发生器的激励信号作对比。实验中的YD64-310型压电式加速度计测得的加速度信号由DHF-2型电荷放大器后转变为一个电压信号。电荷放大器的内部等效电路如图1所示。
机械模态分析 作业:如图1所示是一个单自由系统附件一个减振器形成的的两自由振动系统,已知m 1=105kg ，m 2=7kg ，k 1=10000N/m ，k 2=410N/m ，c 2=1.15N ·m-1·s ，F 1(t)=F 1e j ωt 。求：（简化为粘性比例阻尼进行实模态分析） 1． 物理坐标下的振动微分方程； 2． 频响函数矩阵； 3． 频响函数的模态展式矩阵； 4． 脉冲相应函数； 5． 画出H 11（ω）的幅频特性曲线，相频特性曲线，实频特性曲线， 虚频特性曲线，Nyquist 图，Bode 图； 6． 固有频率，阻尼固有频率； 7． 画出振型图； 8． 模态坐标系下的振动微分方程； 9． 模态参数：复模态质量，复模态刚度，复模态阻尼。 10.按实模态系统，给出灵敏度分析。 11.集全班同学的数据（必要的话再补做不同m 2，k 2，c 2参数下的数据，画出x1的最大振幅与m 2，k 2，c 2，的变化曲线，从而分析出减振器的最佳参数。 解： 1．振动微分方程 对质量m 1、m 2绘分离体图（如图1-1），用牛二定律列分离体在铅垂方向的力平衡方程得 1221221111122122122 ()()()()F c x x k x x k x m x c x x k x x m x ???? ? ? ?? +-+--=----= （1.1） 将（1.1）整理可得： 112 2112 21122 22 2222000m x c c x k k k x F m c c k k x x x ???????? ?? -+-?????????? ????++=??????????????--?? ? ??????????? （1.2） 且m 1=105、m 2=7、k 1=10000、k 2=410、c 2=1.15，代入（1.2）得： ?? ????=????????????+???? ? ???????????+??????????????????????0 410 410-410- 104101.15 1.15- 1.15- 15.17 0 0 1051212121F x x x x x x （1.3） 可以得出此二自由度系统振动微分方程为：()M x C x Kx f t ?? ? ++= 其中M=???? ??7 0 0 105；C=?????? 1.15 1.15- 1.15- 15.1；K=??? ???410 410- 410 - 10410；f(t)= ?? ? ???0 1F 图1-1、系统的分离体图 2．频响函数矩阵 由书P25（1.4-58）公式可知，此二自由度系统频响函数矩阵为一2×2 方阵，其表达式为： 图1 两自由度振动系统
《机械振动基础》实验报告 （2015年春季学期） 姓名 学号 班级 专业机械设计制造及其自动化报告提交日期2015.05.07 哈尔滨工业大学
报告要求 1.实验报告统一用该模板撰写，必须包含以下内容： (1)实验名称 (2)实验器材 (3)实验原理 (4)实验过程 (5)实验结果及分析 (6)认识体会、意见与建议等 2.正文格式：四号字体，行距为1.25倍行距； 3.用A4纸单面打印；左侧装订； 4.报告需同时提交打印稿和电子文档进行存档，电子文档由班长收 齐，统一发送至：liuyingxiang868@https://www.doczj.com/doc/0214193162.html,。 5.此页不得删除。 评语： 教师签名： 年月日
实验一报告正文 一、实验名称：机械振动的压电传感器测量及分析 二、实验器材 1、机械振动综台实验装置(压电悬臂梁) 一套 2、激振器一套 3、加速度传感器一只 4、电荷放大器一台 5、信号发生器一台 6、示波器一台 7、电脑一台 8、NI9215数据采集测试软件一套 9、NI9215数据采集卡一套 三、实验原理 信号发生器发出简谐振动信号，经过功率放大器放大，将简谐激励信号施加到电磁激振器上，电磁激振器振动杆以简谐振动激励安装在激振器上的压电悬臂梁。压电悬臂梁弯曲产生电流显示在示波器上，可以观测悬臂梁的振动情况；另一方面，加速度传感器安装在电磁激振器振动杆上，将加速度传感器与电荷放大器连接，将电荷放大器与数据采集系统连接，并将数据采集系统连接到计算机（PC机）上，操作NI9215数据采集测试软件，得到机械系统的振动响应变化曲线，可以观测电磁激振器的振动信号，并与信号发生器的激励信号作对比。实验中的YD64-310型压电式加速度计测得的加速度信号由DHF-2型电荷放大器后转变为一个电压信号。电荷放大器的内部等效电路如图1所示。 q
实验三：简谐振动幅值测量 一、 实验目的 1、了解振动位移、速度、加速度之间的关系。 2、学会用压电传感器测量简谐振动位移、速度、加速度幅值 二、实验仪器安装示意图 三、 实验原理 由简谐振动方程：)sin()(?ω-=t A t f 简谐振动信号基本参数包括：频率、幅值、和初始相位，幅值的测试主要有三个物理量，位移、速度和加速度，可采取相应的传感器来测量，也可通过积分和微分来测量，它们之间的关系如下： 根据简谐振动方程，设振动位移、速度、加速度分别为x 、v 、a ，其幅值分别为X 、V 、A ： )sin(?ω-=t X x )cos()cos(?ω?ωω-=-==t V t X x v )sin()sin(2?ω?ωω-=--==t A t X x a 式中：ω——振动角频率 ?——初相位 所以可以看出位移、速度和加速度幅值大小的关系是：X V A X V 2ωωω===,。 振动信号的幅值可根据位移、速度、加速度的关系，用位移传感器或速度传感器、加速度传感器进行测量，还可采用具有微积分功能的放大器进行测量。 在进行振动测量时，传感器通过换能器把加速度、速度、位移信号转换成电信号，经过放大器放大，然后通过AD 卡进行模数转换成数字信号，采集到的数字信号为电压变化量，通过软件在计算机上显示出来，这时读取的数值为电压值，通过标定值进行换算，就可计算出振动量的大
小。 DASP 通过示波调整好仪器的状态（如传感器档位、放大器增益、是否积分以及程控放大倍数等）后，要在DASP 参数设置表中输入各通道的工程单位和标定值。工程单位随传感器类型而定，或加速度单位，或速度单位，或位移单位等等。 传感器灵敏度为K CH （PC/U ）（PC/U 表示每个工程单位输出多少PC 的电荷，如是力，而且参数表中工程单位设为牛顿N ，则此处为PC/N ；如是加速度，而且参数表中工程单位设为m/s 2 ，则此处为PC/m/s 2 ）； INV1601B 型振动教学试验仪输出增益为K E ；积分增益为K J （INV1601 型振动教学试验仪的一次积分和二次积分K J =1）； INV1601B 型振动教学试验仪的输出增益： 加速度：K E = 10(mV/PC) 速度：K E = 1 位移：K E = 0.5 则DASP 参数设置表中的标定值K 为： )/(U mV K K K K J E CH ??= 四、 实验步骤 1、安装仪器 把激振器安装在支架上，将激振器和支架固定在实验台基座上，并保证激振器顶杆对简支梁有一定的预压力（不要露出激振杆上的红线标识），用专用连接线连接激振器和INV1601B 型振动教学试验放大仪的功放输出接口。把带磁座的加速度传感器放在简支梁的中部，输出信号接到 INV1601B a 加速度。 2、打开INV1601B 型振动教学试验仪的电源开关，开机进入DASP2006 标准版软件的主界面，选择单通道按钮。进入单通道示波状态进行波形示波。 3、在采样参数设置菜单下输入标定值K 和工程单位m/s 2 ,设置采样频率为4000Hz ，程控倍数1倍。 4、调节INV1601B 型振动教学试验仪频率旋钮到40Hz 左右，使梁产生共振。 5、在示波窗口中按数据列表进入数值统计和峰值列表窗口，读取当前振动的最大值。 6、改变档位v (mm /s )、d (mm )进行测试记录。 7、更换速度和电涡流传感器分别测量a (m /s 2 )、v (mm /s )、d (mm )。
实验 等截面悬臂梁模态测试实验 一、 实验目的 1. 熟悉模态分析原理； 2. 掌握悬臂梁的测试过程。 二、 实验原理 1. 模态分析基本原理 理论上，连续弹性体梁有无限多个自由度，因此需要无限多个连续模型才能描述，但是在实际操作中可以将连续弹性体梁分为n 个集中质量来研究。简化之后的模型中有n 个集中质量，一般就有n 个自由度，系统的运动方程是n 个二阶互相耦合（联立）的常微分方程。这就是说梁可以用一种“模态模型”来描述其动态响应。 模态分析的实质，是一种坐标转换。其目的在于把原在物理坐标系统中描述的响应向量，放到所谓“模态坐标系统”中来描述。这一坐标系统的每一个基向量恰是振动系统的一个特征向量。也就是说在这个坐标下，振动方程是一组互无耦合的方程，分别描述振动系统的各阶振动形式，每个坐标均可单独求解，得到系统的某阶结构参数。 多次锤击各点，通过仪器记录传感器与力锤的信号，计算得到第ｉ个激励点与定响应点（例如点2）之间的传递函数 ω ，从而得到频率响应函数矩阵中的一行 频响函数的任一行包含所有模态参数，而该行的r 阶模态的频响函数 的比值，即为r 阶模态的振型。 2. 激励方法 为进行模态分析，首先要测得激振力及相应的响应信号，进行传递函数分析。传递函数分析实质上就是机械导纳，i 和j 两点之间的传递函数表示 [] ∑==N r iN r i r i r H H H 1 21 ... [] Nr r r N r r r r ir k c j m ???ωω? (2112) ∑ =++-=[]{}[] T r ir N r r iN i i Y H H H ??∑==1 21 ...
机床的模态分析方法综述 甄真 (北京信息科技大学机电工程学院，北京100192) 摘要：模态分析是研究机械结构动力特性的一种近代方法，是结构动态设计及设备的故障诊断的重要方法。机床在工作时,由于要承受各种变载荷而产生振动,其精度和寿命会受到影响。因此有必要对机床进行模态分析,了解其动态特性,以便进一步分析和改进。本文概述了模态分析的概念、研究意义及发展历史，介绍了机床模态分析的研究现状, 从理论方法与试验方法两方面指出了其关键技术以及研究发展方向。 关键词：模态分析；动态特性；机床；理论方法；实验方法 Summary of the model analysis method of machine tool ZHEN Zhen (Beijing Information Science & Technology University, Mechanical and Electrical Engineering College, Beijing, 100192) Abstract：Modal analysis is a modern method to study the dynamic characteristics of mechanical structure. It’s an important method in structure dynamic design and fault diagnosis of equipment.Its accuracy and lifetime will be affected due to withstand all kinds of variable load and vibration when the machine tool works.So it is necessary to make modal analysis and to understand the dynamic characteristics for machine tool in order to further analyze and improve. This paper summarizes the concept, significance and history of modal analysis and introduces the research status of model analysis of machine tool. It also points out the key technology and research direction in this field from two aspects of theoretical method and experimental method. Key words:model analysis; dynamic characteristics; machine tool; theoretical method; experimental method 0 引言 模态是指机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。模态分析是一种研究机械结构动力的方法，是系统辨别方法在工程振动领域中的应用。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析法搞清楚了结构物在某一个易受影响的频率范围内各阶主要模态的特性，就可预言结构在此频段内在外部或内部各种振源作用下实际响应。因此，模态分析是结构动态设计及设备的故障诊断的重要方法[1]。 模态分析将构件的复杂振动分解为许多简单而独立的振动，并用一系列模态参数来表征的过程。根据线性叠加原理，一个构件的复杂振动是由无数阶模态叠加的结果。在这些模态中。模态分析最终目标是识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。模态分析主要分为3类方法：一是，基于计算机仿真的有限元分析法；二是，基于输入(激励)输出(响应)模态试验的试验模态分析法；三是，基于仅有输出(响应)模态试验的运行模态分析法。有限元分析属结构动力学正问题，但受无法准确描述复杂边界条件、结构物理参数和部件连接状态等不确定性因素的限制难以达到很高的精度。第二、三类方法属结构动力学反问题，基于真实结构的模态试验。因而能得到更准确
研究生学院 机械工程专业硕士结课作业 课程题目：机械结构模态分析实验 指导老师： 姓名： 学号： 2015年08月23日
一、概述 模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的，则称为计算模态分析；如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数，称为试验模态分析。通常，模态分析都是指试验模态分析。 振动模态是弹性结构固有的、整体的特性。通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内的各阶主要模态的特性，就可以预言结构在此频段内在外部或内部各种振源作用下产生的实际振动响应。因此，模态分析是结构动态设计及设备故障诊断的重要方法。 机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动模态各不相同。模态分析提供了研究各类振动特性的一条有效途径。首先，将结构物在静止状态下进行人为激振，通过测量激振力与响应并进行双通道快速傅里叶变换（FFT）分析，得到任意两点之间的机械导纳函数（传递函数）。用模态分析理论通过对试验导纳函数的曲线拟合，识别出结构物的模态参数，从而建立起结构物的模态模型。根据模态叠加原理，在已知各种载荷时间历程的情况下，就可以预言结构物的实际振动的响应历程或响应谱。 模态分析的经典定义：将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵，其每列为模态振型。模态分析的最终目标是识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 模态分析技术的应用可归结为以下几个方面： 1) 评价现有结构系统的动态特性； 2) 在新产品设计中进行结构动态特性的预估和优化设计； 3) 诊断及预报结构系统的故障； 4) 控制结构的辐射噪声； 5) 识别结构系统的载荷 二、实验的基本过程 1、动态数据的采集及频响函数或脉冲响应函数分析 （1）激励方法。试验模态分析是人为地对结构物施加一定动态激励，采集各点的振动响应信号及激振力信号，根据力及响应信号，用各种参数识别方法获取模态参数。激励方法不同，相应识别方法也不同。目前主要由单输入单输出（SISO）、单输入多输出（SIMO）多输入多输出（MIMO）三种方法。以输入力的信号特征还可分为正弦慢扫描、正弦快扫描、稳态随机（包括白噪声、宽带噪声或伪随机）、瞬态激励（包括随机脉冲激励）等。 （2）数据采集。SISO方法要求同时高速采集输入与输出两个点的信号，用不断移动激励点位置或响应点位置的办法取得振形数据。SIMO及MIMO的方法则要求大量通道数据的高速并行采集，因此要求大量的振动测量传感器或激振器，试验成本较高。 （3）时域或频域信号处理。例如谱分析、传递函数估计、脉冲响应测量以及滤波、相关分析等。
模态分析实验报告 姓名： 学号： 任课教师： 实验时间： 指导老师： 实验地点：能源与动力工程学院 柴油机拆装实习一楼振动测试实验室
实验1 传递函数的测量 一、实验内容 用锤击激振法测量传递函数。 二、实验目的 1)掌握锤击激振法测量传递函数的方法； 2)测量激励力和加速度响应的时间记录曲线、力的自功率谱和传递函数； 3)分析传递函数的各种显示形式（实部、虚部、幅值、对数、相位）及相干函 数； 4)比较原点传递函数和跨点传递函数的特征； 5)考察激励点和响应点互换对传递函数的影响； 6)比较不同材料的力锤锤帽对激励信号的影响； 三、实验仪器和测试系统 1、实验仪器 主要用到的实验仪器有：冲击力锤、加速度传感器，LMS LMS-SCADAS Ⅲ测试系统，具体型号和参数见表1-1。 仪器名称型号序列号灵敏度备注 数据采集和分析系统LMS-SCADAS Ⅲ比利时力锤LC 3164 4 mV/N 加速度传感器100 mV/g 丹麦B&K 表1-1 实验仪器 2 、测试系统 利用试验测量的激励信号(力锤激励信号)和响应的时间历程信号，运用数字信号处理技术获得频率响应函数(Frequency Response Function, FRF)，得到系统的非参数模型。然后利用参数识别方法得到系统的模态参数。测试系统主要完成力锤激励信号及各点响应信号时间历程的同步采集，完成数字信号的处理和参数的识别。 测量分析系统的框图如图1-1所示。测量系统由振动加速度传感器、力锤和比利时LMS公司SCADAS采集前端及Modal Impact测量分析软件组成。力锤及加速度传感器通过信号线与SCADAS采集前端相连，振动传感器及力锤为ICP
《机械振动学》实验报告 实验名称梁的振动实验 专业航空宇航推进理论与工程 姓名刘超 学号 SJ1602006 南京航空航天大学 Nanjing University of Aeronautics and Astronautics 2017年01月06日
1实验目的 改变梁的边界条件，对比分析不同边界条件，梁的振动特性（频率、振型等）。对比理论计算结果与实际测量结果。正确理解边界条件对振动特性的影响。 2实验内容 对悬臂梁、简支梁进行振动特性对比，利用锤击法测量系统模态及阻尼比等。 3实验原理 3.1 固有频率的测定 悬臂梁作为连续体的固有振动，其固有频率为： ()1,2,.......r r l r ωλ==其中， 其一、二、三、四阶时， 1.87514.69417.854810.9955.....r l λ=、 、、 简支梁的固有频率为： ()1,2,.......r r l r ωλ==其中 其一、二、三、四阶时， 4.73007.853210.995614.1372.....r l λ=、 、、 其中E 为材料的弹性模量，I 为梁截面的最小惯性矩，ρ为材料密度，A 为梁截面积，l 为梁的长度。 试件梁的结构尺寸：长L=610mm, 宽b=49mm, 厚度h=8.84mm. 材料参数： 45＃钢，弹性模量E ＝210 (GPa), 密度ρ＝7800 (Kg/m 3) 横截面积：A ＝4.33*10-4 (m 2), 截面惯性矩：J ＝3 12 bh =2.82*10-9(m 4) 则梁的各阶固有频率即可计算出。
3.2、实验简图 图1 悬臂梁实验简图 图2简支梁实验简图
机械模态分析与实验 学院：机电工程学院 专业：机械制造及其自动化 姓名：马阳 班级：研1302 学号：2013020049
机械结构的模态分析方法研究综述 马阳 摘 要：模态分析是研究机械结构动力特性的一种近代方法，是指通过计算或实验获得机构的固有频率、阻尼比和模态振型等模态参数的过程，是结构动态设计及设备的故障诊断的重要方法。本文对模态分析的基本概念、研究目的、分类、分析方法、发展历程、发展现状和展望一一作了阐述。 关键词：模态分析；模态参数；模态识别；非线性模态 0 引言 模态是指机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。模态分析是一种研究机械结构动力的方法，是系统辨别方法在工程振动领域中的应用[1]。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析法搞清楚了结构物在某一个易受影响的频率范围内各阶主要模态的特性，就可预言结构在此频段内在外部或内部各种振源作用下实际响应。因此，模态析是结构动态设计及设备的故障诊断的重要方法[2]。 1 模态分析概述 1.1 模态分析定义 模态分析的经典定义是：将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵，其每列为模态振型[3]。 模态分析将构件的复杂振动分解为许多简单而独立的振动，并用一系列模态参数来表征的过程。根据线性叠加原理，一个构件的复杂振动是由无数阶模态叠加的结果。在这些模态中。模态分析最终目标是识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊 断和预报以及结构动力特性的优化设计提供依据。 1.2 模态分析基础 1.2.1 无阻尼的情况 实际的机械结构在计算的过程中常常会被简化成多自由度系统。该系统为线性时不变系统。多自由度振动系统可以写成如下的耦合方程形式，可以用矩阵来表示如下： 上式中X 为系统每个自由度的位移向量，对应着系统的各个自由度，M 为系统的质量矩阵，C 为系统的阻尼阵，K 为系统的刚度矩阵，F(t)为系统所受的外力。 为了求解方便，首先考虑没有阻尼的特殊情况（C=0）。 []{()}{()}{()}M X t K X t f t ??+= （1） 设,st X e s jw ψ== 为了求解上式，考虑外力为0时的自由振动齐次解。得到特征方程： 2[][]0Det w M K -+= 由特征方程可以求得特征值（固有频率），与固有频率一一对应可以求得满足式（1）等于0的{}{}{}12,,n ψψψ值，即为求得的特征向量。振动模态之间有正交性，可以证明： ()M X C X KX F t ???++=
模态分析的应用及它的试验模态分析 模态分析是研究结构动力特性一种近代方法，是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得，这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的，则称为计算模记分析；如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数，称为试验模态分析。通常，模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性，就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此，模态分析是结构动态设计及设备的故障诊断的重要方法。 模态分析最终目标是在识别出系统的模态参数，为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 模态分析技术的应用可归结为以下几个方面： 1) 评价现有结构系统的动态特性； 2) 在新产品设计中进行结构动态特性的预估和优化设计； 3) 诊断及预报结构系统的故障； 4) 控制结构的辐射噪声； 5) 识别结构系统的载荷。 机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动千姿百态、瞬息变化。模态分析提供了研究各种实际结构振动的一条有效途径。首先，将结构物在静止状态下进行人为激振，通过测量激振力与胯动响应并进行双通道快速傅里叶变换（FFT）分析，得到任意两点之间的机械导纳函数（传递函数）。用模态分析理论通过对试验导纳函数的曲线拟合，识别出结构物的模态参数，从而建立起结构物的模态模型。根据模态叠加原理，在已知各种载荷时间历程的情况下，就可以预言结构物的实际振动的响应历程或响应谱。 近十多年来，由于计算机技术、FFT分析仪、高速数据采集系统以及振动传感器、激励器等技术的发展，试验模态分析得到了很快的发展，受到了机械、电力、建筑、水利、航空、航天等许多产业部门的高度重视。已有多种档次、各种原理的模态分析硬件与软件问世。在各种各样的模态分析方法中，大致均可分为四个基本过程： （1）动态数据的采集及频响函数或脉冲响应函数分析 1）激励方法。试验模态分析是人为地对结构物施加一定动态激励，采集各点的振动响应信号及激振力信号，根据力及响应信号，用各种参数识别方法获取模态参数。激励方法不同，相应识别方法也不同。目前主要由单输入单输出（SISO）、单输入多输出（SIMO）多输入多输出（MIMO）三种方法。以输入力的信号特征还可分为正弦慢扫描、正弦快扫描、稳态随机（包括白噪声、宽带噪声或伪随机）、瞬态激励（包括随机脉冲激励）等。 2）数据采集。SISO方法要求同时高速采集输入与输出两个点的信号，用不断移动激励点位置或响应点位置的办法取得振型数据。SIMO及MIMO的方法则要求大量通道数据的高速并行采集，因此要求大量的振动测量传感器或激振器，试验成本较高。 3）时域或频域信号处理。例如谱分析、传递函数估计、脉冲响应测量以及滤波、相关分析等。（2）建立结构数学模型根据已知条件，建立一种描述结构状态及特性的模型，作为计算及识别参数依据。目前一般假定系统为线性的。由于采用的识别方法不同，也分为频域建模和时域建模。根据阻尼特性及频率耦合程度分为实模态或复模态模型等。
机械振动实验报告 1.测量简支梁的固有频率和振型 1.1实验目的 用激振法测量简支梁的固有频率和固有振型。掌握多自由度系统固有频和振型的简单测量方法。 1.2实验原理 共振法测量振动系统的固有频率是比较常用的方法之一。共振是指当激振频率达到某一特定值时，振动量的振动幅值达到极大值的现象。本次试验主要利用调整激振频率使简支梁达到位移振动幅值的方法来测量简支梁的一阶，二阶以及三阶固有频率以及从计算机上读取其当时的振型！ 1.3实验内容与结果分析 (1)将激振器通过顶杆连接到简支梁上（注意确保顶杆与激振器的中心线在一直线上），激振点位于简支梁中心偏左50mm 处（已有安装螺孔），将信号发生器输出端分别与功率放大器和数据采集仪的输入端连接，并将功率放大器与激振器相连接。 (2) 用双面胶纸（或传感器磁座）将加速度传感器A 粘贴在简支梁上5#测点（实验时固定不动，用于与其他测点比较相位），将加速度传感器连接，将电荷放大器输出端与数据采集仪的输入端连接。 (3)将信号发生器和功率放大器的幅值旋钮调至最小，打开所有仪器电源。打开控制计算机，打开做此次试验所需的测试软件，进入页面设置好各项参数。通过调节激振频率，观察简支梁位置幅值振动情况。可以通过放在简支梁上的装有一定量塑质小球的小型透明容器直观的观察里面小球的振动情况，小球振动越厉害，也就说明简支梁振动的位移幅值越大；还可以通过分辨简支梁在不同激振频率下的发出的振动声音，声音越大，说明振动幅值越大！ （4）通过（3）中的方法，可以测量出在简支梁在某一激振频率范围内的振动幅值，则此激振频率就是我们需要测量的一阶，二阶以及三阶固有频率，在测出固有频率的同时将计算机上画出的各阶振型的图像保存，以便结果的分析。（5）在完成所有的试验内容之后，通过记录下的实验数据分析实验的结果。所得的实验结果如下：测得的简支梁的一阶、二阶以及三阶的固有频率为 1=35.42HZ ω, 2=131.54HZ ω, 3=258.01HZ ω。振型如下图所示：）
精品资料 悬臂梁各阶固有频率及主振形的测定试验 一、实验目的 1、用共振法确定悬臂梁横向振动时的前五阶固有频率； 2、熟悉和了解悬臂梁振动的规律和特点； 3、观察和测试悬臂梁振动的各阶主振型，分析各阶固有频率及其主振型的实测值与理论计算值的误差。 二、仪器和设备 悬臂梁固定支座；脉冲锤1个；圆形截面悬臂钢梁标准件一个；加速度传感器一个；LMS振动噪声测试系统。 三、实验基本原理 瞬态信号可以用三种方式产生,分述如下： 一是快速正弦扫频法.将正弦信号发生器产生的正弦信号,在幅值保持不变的条件下,由低频很快地连续变化到高频.从频谱上看,该情况下,信号的频谱已不具备单一正弦信号的特性,而是在一定的频率范围内接近随机信号. 二是脉冲激励.用脉冲锤敲击试件,产生近似于半正弦的脉冲信号.信号的有效频率取决于脉冲持续时间τ,τ越小则频率范围越大. 三是阶跃激励.在拟定的激振点处,用一根刚度大、重量轻的弦经过力传感器对待测结构施加张力,使其产生初始变形,然后突然切断张力弦,相当于给该结构施加一个负的阶跃激振力. 用脉冲锤进行脉冲激振是一种用得较多的瞬态激振方法,它所需要的设备较
少,信号发生器、功率放大器、激振器等都可以不要,并且可以在更接近于实际工作的条件下来测定试件的机械阻抗. 四、实验结果记录 前五阶固有频率表 阶数固有频率（Hz） 1 8.491 2 54.216 3 154.607 4 304.354 5 494.691 实验测得的前五阶振型图如下： 1阶振型图
2阶振型图 3阶振型图 4阶振型图
5阶振型图 五、理论计算悬臂梁固有频率 圆截面悬臂钢梁有关参数可取：Pa E 11101.2?=,7850=ρkg/3 m 。用直尺测 量悬臂梁的梁长L=1000mm 、梁直径D=12mm 。计算简支梁一、二、三、四阶固有频率和相应的振型，并将理论计算结果填入表。 悬臂梁的振动属于连续弹性体的振动，它具有无限多自由度及其相应的固有频率和主振型，其振动可表示为无穷多个主振型的叠加。对于梁体振动时，仅考虑弯曲引起的变形，而不计剪切引起的变形及其转动惯量的影响，这种力学分析 模型称为欧拉-伯努利梁。 运用分离变量法，结合悬臂梁一端固定一端自由的边界条件，通过分析可求得均质、等截面悬臂梁的频率方程 1 L Lch cos -=ββ （5-1） 式中：L ——悬臂梁的长度。 梁各阶固有频率为 4 2 2(Al EI l f i i ρπ β）= （5-2） 悬臂梁固有圆频率及主振型函数
第十章周期对称结构的模态分析 ANSYS的周期对称分析支持静力（Static）分析和模态（Modal）分析。静力分析支持线性和大变形非线性；模态分析支持带有预应力的模态分析和不带有预应力的两种，关于带有预应力的模态分析本书第九章有专门讲述。本章只讲述不带有预应力的模态分析。在静力分析和模态分析这两种分析类型中，关于模型建立部分的要求是一致的，不同的是在进行模态分析时需要指定求解的节径数以及指定对于每个节径数的求解的模态阶数。对于每个节径，ANSYS均将其作为一个载荷步。ANSYS将周期对称边界条件施加于每一载荷步，并且每求解一个载荷步（即节径）后，都将构成周期对称边界条件的约束方程删除（保留任何用户自定义的约束方程）。在静力分析中ANSYS只求解零节径，而在模态分析中默认将求解全部节径。 本章中介绍的实例依然是第7章的轮盘，包括模型和边界条件。 10．1 问题描述 某型压气机盘，见7.1节的对其描述。要求查看其低阶频率结构和振动模态。 10．2 建立模型 在周期对称分析中，在建立模型后，划分网格之前，需要指定周期对称选项。 10．2．1 设定分析作业名和标题 在进行一个新的有限元分析时，通常需要修改数据库文件名（原因见第二章），并在图形输出窗口中定义一个标题用来说明当前进行的工作内容。另外，对于不同的分析范畴（结构分析、热分析、流体分析、电磁场分析等）ANSYS6.1所用的主菜单的内容不尽相同，为此我们需要在分析开始时选定分析内容的范畴，以便ANSYS6.1显示出跟其相对应的菜单选项。 （1）选取菜单路径Utility Menu >File >Change Jobname，将弹出修改文件名(Change Jobname)对话框，如图10.1所示。
题目：完成一个综合作业(What I hear, I forgot. What I see, I remember. What I do, I understand.) 作业:如图所示的两自由振动系统,已知m 1=100kg ，m 2=5kg ，k 1=10000N/m ， k 2=500N/m ，c 2=1N ·m-1·s ，F 1(t)=F 1e j ωt 。求： 1． 物理坐标下的振动微分方程； 2． 频响函数矩阵； 3． 频响函数的模态展式矩阵； 4． 脉冲相应函数； 5． 画出H 11（ω）的幅频特性曲线，相频特性曲线，实频特性曲线，虚频特性曲线，Nyquist 图，Bode 图； 6． 固有频率，阻尼固有频率； 7． 画出振型图； 8． 模态坐标系下的振动微分方程； 9． 模态参数：复模态质量，复模态刚度，复模态阻尼。 解： 1．振动微分方程 对质量m 1、m 2绘分离体图（如图1-1），用牛二定律列分离体在铅垂方向的力平衡方程得 12212211111 22122122 ()()()()F c x x k x x k x m x c x x k x x m x ? ? ?? ? ? ?? +-+--=----= （1.1） 将（1.1）整理可得： 1122112 21122 22 222200 0m x c c x k k k x F m c c k k x x x ??????????-+-?????????? ????++=????? ?????????--?? ?????????? ?? （1.2） 且m 1=100、m 2=5、k 1=10000、k 2=500、c 2=1，代入（1.2）得： 111122************ 5000 511500 5000x x x F x x x ??????????--?????????? ????++=??????? ???????--???????? ???? ?? （1.3） 可以得出此二自由度系统振动微分方程为：()M x C x K x f t ?? ? ++= 其中100 00 5M ??=? ???；111 1C -??=??-??；10500 500500 500K -??=??-??；1()0F f x ?? =???? 图1-1、系统的分离体图 2．频响函数矩阵 由书P25（1.4-58）公式可知，此二自由度系统频响函数矩阵为一2×2方阵，其表达式为： 2 1 ()() H K M j C ωωω-=-+，其中100005M ??=? ???；1 111C -?? =? ?-?? ；10500 500500 500K -?? =? ?-?? （2.1） 写成矩阵形式：