内容发布更新时间 : 2025/2/24 7:58:51星期一 下面是文章的全部内容请认真阅读。

基于Matlab 地发动机悬置系统地固有频率和主振型计算.txt看一个人地地心术,要看他地眼神；看一个人地身价,要看他地对手；看一个人地底牌,要看他地朋友.明天是世上增值最快地一块土地,因它充满了希望.基于Matlab 地发动机悬置系统地固有频率和主振型计算<二）

3 运用MATLAB 对动力总成悬置系统固有特性地计算

3.1 理论计算

动力总成系统固有特性地计算, 即计算系统地固有频率和振型.动力总成悬置系统无阻尼地自由振动微分方程:

式中: M——对称正定惯性矩阵。

K——对称正定刚度矩阵.

求多自由度振动系统地固有频率, 从数学上讲就是求特征值地问题:

设式<13） 地解为: X=Xsin<ωt+a）

代入式<13） 化简后得: KX=ω2MX

左乘M- 1 得: M- 1KX=ω2X <14）

令M- 1K=A, 则: AX=ω2X <15）

ω2 即为A 阵地特征值, X 为其特征向量.

由于M 对称正定, K 也是对称阵, 因而式<13） 是广义特征值问题.可用广义特征值地方法求得特征值及特征向量, 所求特征值即为系统地固有频率.b5E2RGbCAP

3.2 MATLAB 计算过程

Matlab 是Matrix Laboratory <矩阵实验室） 地缩写, 它是由美国Mathwork 公司于1967 年推出地软件包, 已发展为一种功能强大地计算机语言, 特别适合于科学与工程计算.p1EanqFDPw

<1） 将动力总成系统质量参数代入式<6） 可得惯性矩阵M.

<2） 将各悬置点地位置参数及悬置块地主刚度代入, 可得EiBiDi.再根据式<12）求得

总体地刚度矩阵K.DXDiTa9E3d

<3） 编制Matlab 程序, 由上述<1）、<2） 得到矩阵M, K, 由式<14）、<15）即可求得A.RTCrpUDGiT

<4） 由式<15）, 通过Matlab 命令eig

4 振动占优方向地判定

在系统定坐标系中, 根据系统地质量矩阵[M] 及振型矩阵, 可以求出系统在做各阶主振动时地能量分布, 将它写成矩阵形式, 定义为能量分布矩阵[EG]j.当系统以第j 阶固有频率振动时,此矩阵地

式中:

[M]kl——质量矩阵地

{u< j）}k——第j 阶振型列阵地第k 个元素。

{u< j）}1——第j 阶振型列阵地第l 个元素。

ωj——为第j 阶固有频率.

根据此矩阵各行元素值总和地大小便可以判别出在系统以第j 阶固有频率振动时地占优方向.

在确定了系统地固有频率和固有振型以及相应地振动占优方向后, 便可以从避开共振频率这一理论来初步评价系统地隔振性能.xHAQX74J0X

目录 文尾口口口 m工程地质计算机应用 2004年 1 期 总33期

--------------------------------------------------------------------------------LDAYtRyKfE

应用MATLAB计算结构自振频率和振型地一种方法

关文阁 杨黎萌 魏翠玲<河北工程学院土木工程系 河北邯郸 056038）

【摘要】 本文从质量归一化原理出发,采用MATLAB工具箱中eig函数,推导出一个新地求解多质点弹性体系自振频率和振型地方法.算例表明该方法比传统地雅可比方法更简捷更易于应用.Zzz6ZB2Ltk

【关键词】自振频率 振型 雅可比法 特征值法 Matlab

当采用地震反应谱方法计算水平地震作用标准值时,需要求得模型结构地多个主振型及其相应地自振频率<周期）.因此,计算多质点弹性体系地自由振动<包括自振周期、振型等）是进行结构抗震设计地必要步骤.dvzfvkwMI1

1 基本理论

多质点弹性体系地无阻尼自由振动方程

(1>

式中： 为体系刚度矩阵； 为体系质量矩阵.

方程<1）左乘 整理后得

(2>

令 则有

(3a>

或 (3b>

这是一个求特征值和特征向量问题.该方程非零解地充要条件,是它地系数行列式等于零,即 (4> rqyn14ZNXI

式<4）称为方阵 地特征方程. 称为方阵 地特征值或特征根.将所求得地 个 逐个回代到式<3b）,解出 , 称为方阵 与 相对应地特征向量,也就是所要求解地第 振型.因此,求解体系地自振频率与振型问题也就是求解方阵 地特征值和特征向量问题.再利用公式(4>求出体系地频率.EmxvxOtOco

(5>

1.1雅可比法

雅可比(Jacobi>法可以求解实对称阵地特征值和特征向量.对式(4>首先需要把 作对称化处理.

令 (6>