4.2多自由度系统的固有频率与主振型

合集下载

第六章 多自由度系统固有频率和主振型的两种近似解法.doc

第六章  多自由度系统固有频率和主振型的两种近似解法.doc

第六章 多自由度系统固有频率和主振型的两种近似解法从多自由度问题的精确解的求解过程可知,求振系的固有频率及主振型是一项必不可少的过程,当自由度较少时,可直接求固有频率及主振型,但当自由度较多时,关于固有频率的求解就很复杂,如一个16自由度的振动问题,仅为展开频率方程的行列式,就需要进行720次计算,当然这些计算可借助计算机解决,但关于固有频率的近似计算及其计算思想,在实际应用及理论研究中仍具有一定的意义。

本章主要介绍求固有频率的两种方法:矩阵迭代法及传递矩阵法。

6-1矩阵迭代法矩阵迭代法适合于自由度较多的复杂系统,该法可以同时计算出系统的固有频率和相应的主振型,当自由度很多,但只要计算出低阶的几个频率时,矩阵迭代法很为适用,其大量的计算可由计算机完成。

在第五章已经介绍过,多自由度无阻尼系统的振动微分方程有两种形式,一种是用刚度矩阵建立的,其固有频率和主振型可由下式求,[]{}[]{}{}02=-A M p A K或写成[]{}[]{}A M p A K 2= (6-1)另一种是用柔度矩阵建立的,其固有频率和主振型可由下式求出{}[][]{}{}012=-A M R A p 或写成{}[][]{}A M R A p=21(6-2) 用[]1-M 前乘(6-1)式,得[]1-M []{}{}A p A K 2= (6-3)方程(6-2)(6-3)可写成如下统一的形式[]{}{}A A D λ= (6-4)(6-4)式称为特征值问题的标准形式,即矩阵迭代法的基本迭代公式。

式中[]D 称为动力矩阵,λ则是矩阵[]D 的特征值,当[]D 是按刚度矩阵形成时,即[][][]K M D 1-=,则λ表示固有频率的平方,λ=p 2,而当[]D 是按柔度矩阵形成时,即[][][]K R D =,则λ表示固有频率的平方的倒数,λ=1/p 2。

显然,任一阶固有频率和主振型都是(6-4)式的精确解。

下面介绍从(6-4)式出发进行迭代的基本过程:1) 某个经过基准化了的初始迭代向量{}1A (所谓基准化就是选取迭代向量的某个分量为基准值1),现选取{}1A ,使其第一个元素A 1,1为基准值1,并作[]{}1AD =运算,运算得到一个新的列阵{}1B ,再将{}1B 基准化,即将新的列阵{}1B 中的各元素均除以B 1,1,可得[]{}{}{}21,111A B B A D ==2) 与{}2A ,如果{}1A ≠{}2A ,则重复上述步骤,以[]D 乘{}2A ,得[]{}{}{}32,122A B B A D ==3) 比{}2A 与{}3A ,如果{}3A ≠{}2A ,则继续重复上述步骤,以[]D 乘{}3A ,…,直到第k 次迭代[]{}{}{}1,1+==k k k k A B B A D ,当式中{}k A ={}1+k A 时停止,这时,特征值1λ=B 1,k ,而相应的特征向量就等于{}k A 。

振动力学[PDF]

振动力学[PDF]

第四章多自由度系统的振动4.1多自由度系统运动方程的建立4.2 耦合与坐标变换4.3 固有频率和主振型4.4振型矩阵、主坐标和正则坐标4.5 固有频率相等的情况4.6 固有频率为零的情况4.7 无阻尼系统对初始条件的响应4.8 无阻尼系统对任意激励的响应4.9 多自由度系统的阻尼4.10 有阻尼系统的响应4.11 一般粘性阻尼系统的响应一般粘性阻尼系统的响应i nj nj j j i j j i nj j j i ==•=••111i n j n j j j i j j i n j j j i Q q k q c q m =++∑∑∑==•=••111•••[][]{}[]{}{}Q q k q c q m =++⎭⎬⎫⎩⎨⎧•••nn n n n n 212222111211212222111211nn n n n n 212222111211212222111211nn n n n n 212222111211212222111211••••••n 11•••n 11n 11n 11inj nj j j i j j i nj j j i ==•=••111i n j n j j j i j j i n j j j i P x k x c x m =++∑∑∑==•=••111•••[][]{}[]{}{}P x k x c x m =++⎭⎬⎫⎩⎨⎧•••n 2121•i ••i1m 12m 23m 3ii •i i Q q Dq T =∂∂+∂∂•1m 12m 23m 32222)2221k +2222•2221⎟⎠⎞+•x c jjj j q W δδ11x δ11P 111x P δ1x δ22P 33P 2⎟⎠⎞21221212212111••••122121221211123323212332321222•••••2332321233232122233323332333••••33323332333•••••••••321333322221321321321333322221321321000000003213213333222213213213333222210•1•1θv ••2θv •1•=1θ•2•=2θ22θ−+mg l 22θl +k Oθ222yk+=•••[][]{}[]{}{}P x k x c x m =++⎭⎬⎫⎩⎨⎧•••••••••n n i j j i i 1111•••nn i j j i i 1111in n i j j i i 1111i j刚度影响系数k i j 若系统各自由度的广义速度和广义加速度为零,除j i i j i 。

多自由度系统的振动模态分析

多自由度系统的振动模态分析

多自由度系统的振动模态分析振动是物体在受到外界作用力或受到初始扰动后产生的周期性运动。

在工程领域中,多自由度系统的振动模态分析是一项重要的研究内容。

本文将介绍多自由度系统的振动模态分析的基本原理和方法。

一、多自由度系统的定义多自由度系统是指由多个相互连接的质点组成的系统。

每个质点都可以在三个坐标方向上自由运动,因此系统的自由度就是质点的个数乘以每个质点的自由度。

多自由度系统的振动模态分析可以帮助我们了解系统的固有振动特性,为工程设计和结构优化提供依据。

二、振动模态的概念振动模态是指多自由度系统在固有频率下的振动形态。

每个固有频率对应一个振动模态,振动模态的数量等于系统的自由度。

振动模态分析可以帮助我们确定系统在不同频率下的振动特性,从而预测系统的响应和寻找可能的共振点。

三、振动模态分析的方法1. 模态分析方法模态分析是一种通过数学方法求解系统的固有频率和振动模态的方法。

常用的模态分析方法包括有限元法、模态超级位置法等。

有限元法是一种基于离散化的方法,将系统分割成有限个小单元,通过求解每个单元的振动特性,最终得到整个系统的振动模态。

模态超级位置法是一种基于物理原理的方法,通过测量系统在不同频率下的振动响应,推导出系统的振动模态。

2. 模态参数的计算模态参数是指描述振动模态特性的参数,包括固有频率、振型、振幅等。

模态参数的计算可以通过实验测量和数值模拟两种方法。

实验测量是通过激励系统,测量系统在不同频率下的振动响应,并通过信号处理和频谱分析等方法计算出模态参数。

数值模拟是通过建立系统的数学模型,利用计算机仿真软件求解系统的振动模态。

四、振动模态分析的应用振动模态分析在工程领域有广泛的应用。

首先,振动模态分析可以帮助工程师了解系统的固有振动特性,从而优化设计和改善结构。

其次,振动模态分析可以用于故障诊断和预测,通过对系统的振动模态进行监测和分析,可以判断系统是否存在异常或潜在故障。

此外,振动模态分析还可以应用于声学工程、航天工程、汽车工程等领域。

机械振动运动学第四章 多自由度系统振动(改)

机械振动运动学第四章  多自由度系统振动(改)
(4.19)
或简写成
上式还可以简写成:
(4.21)
(4.20)
上式表明,在动力作用下系统产生的位移等于系统的柔 度矩阵与作用力的乘积。它也可写成:
(4.22) 柔度矩阵与刚度矩阵之间转换关系为:
(4.23)
上式说明,对于同一个机械振动系统,若选取相同的广 义坐标,则机械振动系统的刚度矩阵和柔度矩阵互为逆矩矩 阵。
可用矩阵形式表达为:
(4.48)
(4.49)
(4.50) (4.51) 将式(4.50)和式(4.51)代入式(4.48)和式(4.49) 中,得到机械系统的动能T和势能V的表达式分别为:
(4.52)
故得
(4.53) (4.54)
(4.55)
单自由度无阻尼系统在作自由振动时,其动能T和势能V (4.57) (4.58)
现在选取以下三组不同的广义坐标来分别写出振动系统 的运动作用力方程。
①取C点的垂直位移 yc和刚杆绕C点的转角c为广义坐标。 如图4.6(b)所示。
图4.6(b) 刚体振动系统广义坐标示意图 应用达朗伯原理,得出振动系统的运动方程式:
(4.62)
将上式写成矩阵形式:
(4.63)
上式中,刚度矩阵是非对角线矩阵,反映在方程组中,即 为两个方程通过弹性力项互相耦合,故称为弹性耦合。
为使系统的第 j坐标产生单位位移,而其它坐标的位移 为零时,在第i 坐标上所需加的作用力大小。
现以图4.1所示的三自由度系统为例,说明确定影响系数和 系数矩阵的方法。
1、确定 及[k] 设 x₁ 1, x₂ 0,x₃ 0 则得到系统的刚度矩阵
2、确定 及[C] 设 设 设
得 C₁₁ C₁ C₂, C₂₁ C₂, C₃₁ ; 得 C₂₂ C₂ C₃;C₁₂ C₂;C₃₂ C₃ 得C₃₃ = C₃; C₂₃ = C₃; C₁₃ = 0

第4章多自由度系统的振动题解

第4章多自由度系统的振动题解

习 题4-1 在题3-10中,设m 1=m 2=m ,l 1=l 2=l ,k 1=k 2=0,求系统的固有频率和主振型。

解:由题3-10的结果22121111)(l g m l g m m k k +++=,2221l gm k -=,2212l g m k -=,22222l gm k k += 代入m m m ==21,021==k k ,l l l ==21 可求出刚度矩阵K 和质量矩阵M⎥⎦⎤⎢⎣⎡=m m M 00;⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=l mg lmg l mg l mg K 3 由频率方程02=-M p K ,得0322=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=mp l mgl mg lmgmp lmg B 0242222242=+-∴l g m p l g m p m l g p )22(1-=∴ ,lgp )22(2+= 为求系统主振型,先求出adjB 的第一列⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=l mg mp lmg adjB 2分别将频率值21p p 和代入,得系统的主振型矩阵为⎥⎦⎤⎢⎣⎡-=112)1(A⎥⎦⎤⎢⎣⎡+=112)2(A题4-1图4-2 题4-2图所示的均匀刚性杆质量为m 1,求系统的频率方程。

解:设杆的转角θ和物块位移x 为广义坐标。

利用刚度影响系数法求刚度矩阵k 。

设0,1==x θ,画出受力图,并施加物体力偶与力2111,k k ,由平衡条件得到,222111a k b k k +=, a k k 221-=设1,0==x θ,画出受力图,并施加物体力偶与力2212,k k ,由平衡条件得到, 12k a k 2-=, a k k 222= 得作用力方程为⎥⎦⎤⎢⎣⎡=⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡--++⎭⎬⎫⎩⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡0000312222221221x a k a k a k a k b k x m a m θθ由频率方程02=-M K p ,得031222222212221=----+p m a k ak a k p a m a k b k4-3 题4-3图所示的系统中,两根长度为l 的均匀刚性杆的质量为m 1及m 2,求系统的刚度矩阵和柔度矩阵,并求出当m 1=m 2=m 和k 1=k 2=k 时系统的固有频率。

机械振动学中的固有频率与振型分析

机械振动学中的固有频率与振型分析

机械振动学中的固有频率与振型分析机械振动学是研究机械系统在受到外界激励作用下产生振动现象的一门学科。

在机械系统中,固有频率与振型分析是非常重要的内容,可以用来描述系统的动态特性和振动行为。

本文将介绍机械振动学中固有频率与振型分析的基本概念和应用。

一、固有频率固有频率是指机械系统在没有外界激励下自由振动的频率。

对于一个简单的振动系统,其固有频率可以通过运动方程的解析解求得。

固有频率是系统的固有特性之一,可以用来描述系统的动态响应特性和结构的刚度、质量、阻尼等参数。

在实际工程应用中,固有频率的计算对于系统结构设计和振动控制至关重要。

通过对系统的固有频率进行分析,可以避免共振现象的发生,减小系统动态响应,提高系统的稳定性和可靠性。

二、振型分析振型分析是指对机械系统的振动模式和振动幅值进行分析和描述。

振型是指系统在特定频率下的振动模式,可以通过振动实验和有限元分析等方法得到。

振型分析可以提供系统的模态形式和振动幅值信息,有助于分析系统的受力情况和结构设计。

振型分析在工程实践中具有广泛的应用,可以用于评估系统的结构健康状况、辅助设计优化和振动控制。

通过对系统的振型进行分析,可以找到系统的薄弱环节和潜在问题,及时进行改进和优化,提高系统的性能和可靠性。

三、结语固有频率与振型分析是机械振动学中重要的内容,对于机械系统的设计和性能评估具有重要意义。

通过对系统的固有频率和振型进行分析,可以优化系统的结构设计,降低系统的动态响应,提高系统的稳定性和可靠性。

希望本文的介绍能够帮助读者更好地理解机械振动学中固有频率与振型分析的相关知识。

第2章——多自由度系统的振动——多自由度方程的建立

第2章——多自由度系统的振动——多自由度方程的建立

船体振动基础1第2章多自由度系统的振动第章多自由度系统的振一、引言二、两自由度系统的振动2上节课内容的回顾1.几个重要概念主振型第阶主振型第二阶主振型多自由度系统主振型,第一阶主振型,第二阶主振型基频,第一阶固有频率,第二阶固有频率,……主振动,模态个自度系自上节课内容的回顾2.两个自由度系统的自由振动(P37)⎬⎫=++−=−++00)(2212111x k k x k xm x k x k k xm &&&&⎭)(2321222个自度系自上节课内容的回顾2.两个自由度系统的自由振动(P41-43)m &&⎭⎬⎫=++−=−++0)(0)(23212222212111x k k x k xm x k x k k x&&①假设简谐形式的解振动时,两个质量按相同频率和相位角作简谐振动。

()()⎭⎬⎫+=+=θωθωt A x t A x n n sin sin 2211上节课内容的回顾将简谐振动解代入运动方程式上节课内容的回顾解特征方程式的根,可以得到:上节课内容的回顾将特征值代入②的振幅A1和振幅A2,得到对应于和的振幅A1和振幅A2之间的两个确定的比值:21ω上节课内容的回顾⑥主振动的确定。

z 系统以某一阶固有频率按其相应的主振型作振动,z 称为系统的主振动(1)(1)⎫第一阶主振动为:()1111(1)(1)(1)22111111sin()sin()sin x A t xA t A t ωθωθβωθ=+⎪⎬=+=+⎪⎭第二阶主振动为:(2)(2)1122sin()x A t ωθ⎫=+⎪()(2)(2)(2)22222122sin()sin x A t A t ωθβωθ⎬=+=+⎪⎭z 系统作主振动时,各点同时经过静平衡位置和到达最大偏离位置,z 以确定的频率和振型作简谐振动。

上节课内容的回顾⑦一般情况下自由振动的通解。

并非在任何情况下系统都会作主振动形式的运动,一般情况下系统运动方程的通解为上述两种主振动的叠加:o在一般情况下,系统的自由振动是两种不同频率的主振动的线性组合,其结果不一定是简谐振动。

多自由度系统振动理论及应用

多自由度系统振动理论及应用
多自由度系统的作用力方程
对一些较简单的问题,用牛顿定律来建立振动微分方程是简便的.
图4-1所示为无阻尼三自由度弹簧质量系统,可参照二自由度系统的方
法,写出其微分方程:

下一页
返回
4.1
多自由度系统的振动微分方程

或更一般地写成

该式可简单地写成

式(4-2)称为用矩阵符号表示的作用力方程,它可以代表许多种运动方程
种心灵的孤独。
2. 与 个 别 人 难 以 相 处
一些学生能够与多数人保持良好的关系,但与个别人交往
不 良 。 因 此 ,常 会 影 响 情 绪 ,如 鲠 在 喉 。
上一页 下一页
返回
任 务 一了解自己与人交往的现状

3. 与 他 人 交 往 平 淡

一些学生虽然能与他人交往,但多属点头之交,没有关系
人际关系新起点
1
任 务 一 了解自己与人交往的现状
2
任 务 二 调整不良交际心态
任 务 一了解自己与人交往的现状







任 务 提 出 :了 解 自 己 与 人 交 往 的 现 状 。
任 务 目 标 :了 解 自 己 与 人 交 往 的 现 状 ,激 发 学 习 热 情 ,明 确 努
力方向。
喜欢独来独往。

(3) 嫉 妒 心 理 。 部 分 大 学 生 不 能 正 确 对 待 别 人 的 长 处 和 优
点,看到别人冒尖心里嫉妒,对比自己水平高的同学采取
讽 刺 、 挖 苦 、 打 击 、 嘲 笑 等 不 当 方 式 ,给 别 人 造 成 伤 害 ,严
重影响了同学之间的沟通。
上一页
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

4.2 多自由度系统的固有频率与主振型自由振动微分方程:
设它为:
可得如下主振型方程
(4-11)与(4-12)化成具有相同的形式,对(4-11)式两端乘以,可得
可化为
)可化为
)就有着相同的形式。

柔度矩阵之间存在着互逆关系,即有

两种系统矩阵之间有着互逆关系:
、柔度矩阵以及质量矩阵一般都是对称矩阵,但是其系统矩阵和一般已不再是对称矩阵。

振型问题。

鉴于方程(4-15)与(4-17)属于同一形式,故只需讨论其中之一。

征方程。

对它展开的结果,可得一个关于的次代数方程:
(4-20)的特征根,亦称矩阵的特征值。

特征值与系统固有频率之间有如下关系:
可以是单根,也可以是重根;可以是实数,也可以是复数。

但是,在我们所考虑的情形中,由于系统质量矩阵是正定的实对称阵,刚事实上,由正定与半正定的条件,对于任何非零的,有
上述结论。

半正定)的系统,称为正定系统(或半正定系统)。

所以,上述结论可改述为:正定系统的特征值都是正的,而半正定系统的特征
18),可求得各个相应的,他们称为系统的主振型(或固有振型),亦称为矩阵的特征矢量。

这样,对于任何一个自由度系征矢量)。

统的特征矢量也可以从的伴随矩阵得出。

事实上,按逆阵的表示,有
中各列与充其量只相差一个常数乘子。

统中,有,。

求系统的主振型。

(4-23)中第三列正是取为基准的主振型:
23),可得

形式:
,是特征值。

在大多数算法中还假设是对称阵。

17)都具有(4-24)式的形式。

不过无论是还是一般都不是对称阵。

为了将它们化为对称阵,可进行如下坐标变换。

的实对称阵,它总可以分解为
转置阵。

)与(4-17),得
是对称阵,并且与系统矩阵有着相同的特征值。

事实上,是由对称阵经变换(4-27)得来的,按矩阵乘积的的转置规则,
(4-30)
阵经变换(4-30)得来的。

因为
有着相同的特征值。

不过,这时的特征矢量已不同于特征矢量,但如果求出后,就不难通过逆变换:
,对式(4-25)求逆,有
,式(4-32)可改写为
系统矩阵已是对称阵,且有
着相同得特征值。

换将大为简化。

这时有
为中对应元素的平方根;的元素分别为中对应元素的倒数。

从而有
统中有,,。

试将系统矩阵化为对称阵。



数的一个专题。

已经发展了许多有效的算法来求解各种形式的矩阵的特征值问题。

关于这一问题的详细论述,请读者参阅有关专著
交性。

这种正交性表现为关于质量矩阵与刚度矩阵的加权正交性,即当时,有

个与第个主振型,因而有

二式前乘以,然后两者相减,可得
)中的第一式得证;同理可证明其第二式。

取何有限值,式(4-38)恒成立,因而可取


包含一个任意常数乘子,因而可以选取归一化振型,使式(4-39)中各个都等于1,这时各个就等于。

/jp2009/04/wlkc2/dd04/4-2.htm。

相关文档
最新文档