模态分析理论
- 格式:docx
- 大小:358.19 KB
- 文档页数:9
模态分析的基础理论模态分析是一种研究系统中不同模式的分布、生成和演化规律的方法。
在这个理论中,模态是指系统中不同状态或形式的存在形式,例如质量分数、温度、湿度等。
模态分析的基础理论包括概率论、统计学和模态分析技术等。
概率论是模态分析的基础之一、它研究随机事件的发生概率和规律。
在模态分析中,我们可以利用概率论来描述不同模态出现的概率分布,并通过分析系统中的模式,得出不同模态的生成规律。
通过概率论的方法,我们可以预测不同模态的变化趋势,从而指导系统的优化设计和运行管理。
统计学也是模态分析的基础理论之一、统计学研究如何收集、处理、分析和解释数据,通过对大量数据的统计分析,揭示数据背后的规律和趋势。
模态分析中,统计学的方法可以用于分析模态数据的分布情况,寻找模态之间的相关性和影响因素,并建立相应的模型来预测和优化系统的运行情况。
在模态分析技术方面,主要包括聚类分析、主成分分析和模态分析方法等。
聚类分析是一种将相似的对象分组的方法,通过对模态数据进行聚类分析,我们可以将相似的模态归为一类,从而描述系统中的不同模态分布情况。
主成分分析是一种降维技术,它可以将高维的模态数据降低到低维,并保留大部分信息。
这可以帮助我们更好地理解系统模态之间的关系和重要性。
模态分析方法包括有限元模态分析、频响函数法和模态参数识别等。
通过这些方法,我们可以对系统的模态进行分析,包括振型、频率和阻尼等,并找出模态的摄动源和分布规律。
模态分析的基础理论对于理解和优化系统具有重要意义。
通过对模态的分析和研究,我们可以了解系统的特性和不同模态之间的关系,从而指导系统的设计和运行。
同时,模态分析也可以帮助我们发现和解决系统中存在的问题,提高系统的稳定性和可靠性。
因此,深入理解和应用模态分析的基础理论对于各个领域的研究和实践具有重要价值。
模态叠加法一.思想要点是在积分运动方程以前,利用系统自由振动的固有振型将方程组转换为n 个相互不耦合的方程,对这种方程可以解析或数值地进行积分。
对于每个方程可以采用各自不同的时间步长,即对于低阶振型可采用较大的时间步长。
当实际分析的时间历程较长,同时又只需要少数较低阶振型的结果时,采用振型叠加法将是十分有利的。
求解步骤:1.求解系统的固有频率和振型2.求解系统的动力响应二.求解固有频率与振型(求解不考虑阻尼影响的振动方程) ..()(){0}M a t Ka t += 解可假设为:0sin ()a t t φω=-φ是n 阶向量,ω是向量φ的振动频率,t 是时间变量,0t 是由初始条件确定的时间常数。
代入振动方程,得到一个广义特征值问题:20K M φωφ-=求解可得n 个特征解221122(,),(,),ωφωφ···2,(,)n n ωφ120ωω≤<<···n ω< 特征向量12,,φφ···,n φ代表系统的n 个固有振型,幅度可按以下要求规定T i i M φφ=1(i=1,2,···,n ),这样规定的固有振型又称正则振型。
将22(,)(,)i i j j ωφωφ代回特征方程,得:2i i i K M φωφ= 2j j j K M φωφ=前式两边前乘以j φT,后式两边前乘以i φT ,得:2j i i j i K M φφωφφTT = 2i j i i jK M φφωφφT T = 由()TTj i j i i j K K K φφφφφφT T==得:22i j i j i j M K ωφφωφφT T =,推出22()0i j j i M ωωφφT-=当i j ωω≠时,有0j i M φφT =这表明固有振型对于矩阵M 是正交的,可表示为:1 ()0 ()i j i j M i j φφT=⎧=⎨≠⎩得:2 ()0 ()i i j i j K i j ωφφT ⎧==⎨≠⎩如果定义123n [ ]φφφφΦ=K21222 0 0 n ωωω⎡⎤⎢⎥⎢⎥Ω=⎢⎥⎢⎥⎢⎥⎣⎦O则特征解的性质可表示成:M K T T ΦΦ=I ΦΦ=Ω原特征值问题可表示为:K M Φ=ΦΩ三.求解动力响应1.位移基向量的变换引入变换()()1ni i i a t x t x φ==Φ=∑其中()[]12 n x t x x x =L代入运动方程,并两边前乘以T Φ,可得:()()()()()...x t C x t x t Q t R t T T +ΦΦ+Ω=Φ= 初始条件相应地转换成:..0000 x x Ma M a T T =Φ=Φ 阻尼为振型阻尼,则:()()2 i=j 0 i j i i ij C ωξφφT ⎧⎪=⎨≠⎪⎩ 或11222 0 2 0 2n n C ωξωξωξT ⎡⎤⎢⎥⎢⎥ΦΦ=⎢⎥⎢⎥⎣⎦O 其中i ξ(i=1,2,···,n )是第i 阶振型阻尼比,可得n 个相互不耦合的二阶常微分方程()()()()...22i i i i i i i x t x t x t r t ωξω++= (i=1,2,···,n )若C 是Rayleigh 阻尼,即C M K αβ=+根据试验或相近似结构的资料已知两个振型的阻尼比i ξ和j ξ,可得22222()()2()()i j j i i j j i j j i i j i ξωξωαωωωωξωξωβωω-=--=-2.求解单自由度系统振动方程在振动分析中常常采用杜哈美(Duhamel )积分,又称叠加积分,其基本思想是将任意激振力()i r t 分解为一系列微冲量的连续作用,分别求出系统对每个微冲量的响应,然后根据线性系统的叠加原理,将它们叠加起来,得到系统对任意激振的响应。
模态分析理论1模态分析简介1.1 模态简介模态是结构固有的振动特性,每一个模态具有一个特定的固有频率、阻尼比和模态振型。
这些模态参数可以由分析软件分析取得,也可以经过试验计算获得,这样一个软件或者试验分析过程称为模态分析。
这个分析结果如果是由有限元计算的方法取得的,则称为计算模态分析;如果结果是通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。
模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。
1.2 固有频率简介固有频率是物体的一种物理特性,由它的结构、大小、形状等因素决定的。
这种物理特征不以物体是否处于振动状态而转移。
当物体在多个频率上振动时会渐渐固定在某个频率上振动,当他受到某一频率策动时,振幅会达到最大值,这个频率就是物体的固有频率。
1.3 振型简介振型是指体系的一种固有的特性。
它与固有频率相对应,即为对应固有频率体系自身振动的形态。
每一个物体实际上都会有无穷多个固有频率,每一阶固有频率相对应物体相对应的形状改变我们称之为振型。
理论上来说振型也有无穷多个,但是由于振型阶数越高,阻尼作用造成的衰减越快,所以高振型只有在振动初期才较明显,以后则衰减。
因此一般情况下仅考虑较低的几个振型.1.4模态分析的目的模态分析技术从上世纪60年代开始发展至今,已趋于成熟。
它和有限元分析技术一起,已成为结构动力学中的两大支柱。
到目前,这一技术已经发展成为解决工程振动问题的重要手段,在机械、航空航天、土木建筑、制造化工等工程领域被广泛的应用。
我国在这一方面的研究,在理论上和应用上都取得了很大的成果,处于世界前列。
模态分析的最终目标就是识别出系统的模态参数,为结构系统的振动特性的分析、振动故障的诊断和检测以及结构的优化提供依据。
模态分析技术的应用可归结为以下几个方面:1) 评价所求结构系统的动态特性;2) 在新产品设计中进行结构特性的预估,优化对结构的设计;3) 诊断及预报结构系统中的故障;4) 识别结构系统的载荷。
模态分析理论1模态分析简介1.1 模态简介模态是结构固有的振动特性,每一个模态具有一个特定的固有频率、阻尼比和模态振型。
这些模态参数可以由分析软件分析取得,也可以经过试验计算获得,这样一个软件或者试验分析过程称为模态分析。
这个分析结果如果是由有限元计算的方法取得的,则称为计算模态分析;如果结果是通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。
模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。
1.2 固有频率简介固有频率是物体的一种物理特性,由它的结构、大小、形状等因素决定的。
这种物理特征不以物体是否处于振动状态而转移。
当物体在多个频率上振动时会渐渐固定在某个频率上振动,当他受到某一频率策动时,振幅会达到最大值,这个频率就是物体的固有频率。
1.3 振型简介振型是指体系的一种固有的特性。
它与固有频率相对应,即为对应固有频率体系自身振动的形态。
每一个物体实际上都会有无穷多个固有频率,每一阶固有频率相对应物体相对应的形状改变我们称之为振型。
理论上来说振型也有无穷多个,但是由于振型阶数越高,阻尼作用造成的衰减越快,所以高振型只有在振动初期才较明显,以后则衰减。
因此一般情况下仅考虑较低的几个振型.1.4模态分析的目的模态分析技术从上世纪60年代开始发展至今,已趋于成熟。
它和有限元分析技术一起,已成为结构动力学中的两大支柱。
到目前,这一技术已经发展成为解决工程振动问题的重要手段,在机械、航空航天、土木建筑、制造化工等工程领域被广泛的应用。
我国在这一方面的研究,在理论上和应用上都取得了很大的成果,处于世界前列。
模态分析的最终目标就是识别出系统的模态参数,为结构系统的振动特性的分析、振动故障的诊断和检测以及结构的优化提供依据。
模态分析技术的应用可归结为以下几个方面:1) 评价所求结构系统的动态特性;2) 在新产品设计中进行结构特性的预估,优化对结构的设计;3) 诊断及预报结构系统中的故障;4) 识别结构系统的载荷。
精心整理模态分析指的是以振动理论为基础、以模态参数为目标的分析方法。
首先建立结构的物理参数模型,即以质量、阻尼、刚度为参数的关于位移的振动微分方程;其次是研究其特征值问题,求得特征对(特征值和特征矢量),进而得到模态参数模型,即系统的模态频率、模态22¨330m 0z k 2k k z 000m 0k k z 0z +--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦(9) 定义主振型由于是无阻尼系统,因此系统守恒,系统存在振动主振型。
主振型意味着各物理坐标振动的相位角不是同相(相差0o )就是反相位(相差180o ),即同时达到平衡位置和最大位置。
主振型定义如下:()i i j ωt+i i sin ωt+=Im(e )φφi mi mi z =z z (10)其中为第i 阶频率下,各自有度的位移矢量,为第i 个特征矢量,表示第i 阶固有频率下的振型,i ω为第i 阶频率下的第i 个特征值,i φ为(去除项化简得以矩阵的形式展开得:2i 2i mi 2i k-ωm -k 0-k 2k-ωm -k z =00-k k-ωm ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(15)有非零解,则2i 2i 2i k-ωm -k 0-k 2k-ωm -k =00-k k-ωm ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(16)即()234222ω-m ω+4km ω-3k m =0(17)阶固有频率,每一个特征根对应一个特征矢量,表示对应模态下该由式3i i 21=z k 如果设定了1z 值,则就可以求出三个特征根值下,2z 和3z 相对于1z 的位移。
假设m=k=1, 一阶模态,1ω=0:21z =1z ,31z =1z ,即;二阶模态,223kω=m :21z=0z,31z=-1z,即;三阶模态,23kω=m :21z=-2z,31z=1z,即。
运动方程的解耦图错误!未指定顺序。
运动方程解耦过程在进行坐标变换之前需对刚度矩阵和质量矩阵进行归一化。
模态分析指的是以振动理论为基础、以模态参数为目标的分析方法。
首先建立结构的物理参数模型,即以质量、阻尼、刚度为参数的关于位移的振动微分方程;其次是研究其特征值问题,求得特征对(特征值和特征矢量),进而得到模态参数模型,即系统的模态频率、模态矢量、模态阻尼比、模态质量、模态刚度等参数。
特征根问题以图3所示的三自由度无阻尼系统为例,设123m =m =m =m ,123k =k =k =k ,图 1 三自由度系统其齐次运动方程为:(8)其中分别为系统的质量矩阵和刚度矩阵,123m 00m 00m=0m 0=0m 000m 00m ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦,11212221k -k 0k -k 0k=-k k +k -k =-k 2k -k 0-k k 0-k k ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦,则运动方程展开式为: ¨11¨22¨33z m 00k k 0z 00m 0z k 2k k z 000m 0k k z 0z ⎡⎤⎢⎥-⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥+--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦(9)定义主振型由于是无阻尼系统,因此系统守恒,系统存在振动主振型。
主振型意味着各物理坐标振动的相位角不是同相(相差0o )就是反相位(相差180o ),即同时达到平衡位置和最大位置。
主振型定义如下:()i ij ωt+i i sin ωt+=Im(e)φφi mi mi z =z z (10)其中为第i 阶频率下,各自有度的位移矢量,为第i 个特征矢量,表示第i 阶固有频率下的振型,i ω为第i 阶频率下的第i 个特征值,i φ为初始相位。
对于三自由度系统,在第i 阶频率下,等式可以写成1m1i 2m2i i i 3m3i z z z =z sin(ωt+)z z φ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(11)mki z 表示第k 个自由度在第i 阶模态下的模态矩阵。
各种模态分析方法总结与比较一、模态分析模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。
模态分析的理论经典定义:将线性定常系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,以便求出系统的模态参数。
坐标变换的变换矩阵为模态矩阵,其每列为模态振型。
模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。
模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。
这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。
这个分析过程如果是由有限元计算的方法取得的,则称为计算模记分析;如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。
通常,模态分析都是指试验模态分析。
振动模态是弹性结构的固有的、整体的特性。
如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。
因此,模态分析是结构动态设计及设备的故障诊断的重要方法。
模态分析最终目标是在识别出系统的模态参数,为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。
二、各模态分析方法的总结(一)单自由度法一般来说,一个系统的动态响应是它的若干阶模态振型的叠加。
但是如果假定在给定的频带内只有一个模态是重要的,那么该模态的参数可以单独确定。
以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。
在给定的频带范围内,结构的动态特性的时域表达表示近似为:()[]}{}{T R R t r Q e t h rψψλ= 2-1而频域表示则近似为:()[]}}{{()[]2ωλωψψωLR UR j Q j h r tr r r -+-= 2-2 单自由度系统是一种很快速的方法,几乎不需要什么计算时间和计算机内存。
模态分析理论范文模态分析理论的核心理念是,人们在特定的社会和文化情境下会表现出不同的态度和行为。
它认为,我们的态度、信念和行为不仅受到个体心理因素的影响,还受到社会和文化环境的影响。
因此,要全面了解一个人的态度或行为,就需要考虑到这个人所处的情境。
首先,模态指的是人们在特定情境下所采取的态度、信念和行为。
它可以通过探究个体的思考方式、观点和意见来理解。
例如,一些人可能会对一些产品持有积极的态度,这可能是因为他对产品的特点和功能有较高的认同。
其次,資源指的是人们在模态形成过程中所依赖的信息和知识。
在分析模态时,人们使用各种不同的资源来评估和形成自己的态度。
这些资源可以是个体的经验、心理特征、社会身份或文化价值观。
通过了解人们所依赖的资源,我们可以更好地理解他们的态度和行为。
最后,情境是指人们所处的社会和文化环境。
情境对个体的态度和行为具有重要的影响。
在不同的情境下,人们可能表现出不同的态度和行为。
例如,同一个人在工作时可能持有不同的观点和做法,而在家庭生活中可能又是另一种态度和行为。
模态分析理论的应用非常广泛。
在广告和市场营销领域,模态分析理论被用于理解消费者的态度和行为,从而更好地设计和推广产品。
在政治和公共政策领域,模态分析理论可以帮助政治家和政策制定者了解公众的意见和需求,有针对性地制定政策和决策。
在社会学和心理学领域,模态分析理论可以用来研究群体行为和态度的变化,揭示社会和文化因素对个体的影响。
然而,模态分析理论也存在一些限制。
首先,因为人们的态度和行为是受多个因素的影响,所以模态分析理论不能解释所有的情况。
其次,模态分析理论强调了情境对个体行为的影响,但情境本身也是由个体创造和改变的,所以情境也会受到个体行为的影响。
最后,模态分析理论对于一些复杂的社会和文化现象可能无法提供充分的解释,因为这些现象涉及多个层面和多个因素的交互作用。
总之,模态分析理论是一种有用的社会科学研究方法,可以帮助我们理解人们在特定情境下的态度和行为。
Tyler & Sofrin 模态分析理论非定常来流与叶片干涉产生的声波在风扇或涡轮中并非是任意形态存在的。
Goldstein 在假定平均流场有势的前提下,建立起了平均流场中任意一点的扰动量与远前方来流扰动量之间的相互关系,给出了下面的方程,()()00002000111I D D Dt c Dt ϕρϕρρρ⎛⎫-∇∇=∇ ⎪⎝⎭u (0-1)由以上方程可知,在非均匀平均流的情况下,来流扰动不仅通过边界条件与声扰动相互作用,而且在传播过程中也会与声扰动耦合,并形成如(2-2)右边所示声源[68]。
在航空发动机叶轮机内部,最重要的边界条件就是管道效应,由于管道边界的限制,声波在其中只能以特定的形态出现,也就是我们常说的模态。
在均匀平均流中,考虑一个环形管道,硬壁条件,对小扰动有下面的对流波动方程[12],2222222110i M p p x x r r r r ωυ⎛⎫∂∂∂∂∂⎛⎫+-+++= ⎪ ⎪∂∂∂∂∂⎝⎭⎝⎭ (0-2)波动方程描述的特征值问题是可解的,环形管道中我们可以将它的一般解展开为傅里叶-贝塞尔形式的模态()()()1,,m m ik x ik x im m m m m p x r A e B e U r e μμθμμμμθ+-∞∞---=-∞==+∑∑ (0-3)这里径向模态和径向、轴向波数分别满足()22222210m m m m m m m m m U U U r r Mk k k μμμμμμμμααω±⎛⎫'''++-= ⎪⎝⎭=--= (0-4) 其中,径向特征模态()m U r μ以贝塞尔函数的形式出现,m 和μ分别表示周向和径向模态数。
满足上述波动方程的声波解在环形或圆形管道中会以图2-4所示的螺旋波形式出现和传播。
Tyler 和Sofrin 是最早研究叶轮机内部叶片非定常气动力旋转模态特征的学者,他们的研究结果已经成为当代航空燃气涡轮发动机气动声学设计的主要理论基础之一。
模态分析理论Document number:WTWYT-WYWY-BTGTT-YTTYU-2018GT模态分析指的是以振动理论为基础、以模态参数为目标的分析方法。
首先建立结构的物理参数模型,即以质量、阻尼、刚度为参数的关于位移的振动微分方程;其次是研究其特征值问题,求得特征对(特征值和特征矢量),进而得到模态参数模型,即系统的模态频率、模态矢量、模态阻尼比、模态质量、模态刚度等参数。
特征根问题以图3所示的三自由度无阻尼系统为例,设123m =m =m =m ,123k =k =k =k ,图三自由度系统其齐次运动方程为:mz̈+kz =0(8)其中m ,k 分别为系统的质量矩阵和刚度矩阵,123m 00m 00m=0m 0=0m 000m 00m ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦,11212221k -k 0k -k 0k=-k k +k -k =-k 2k -k 0-k k 0-k k ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦,则运动方程展开式为:¨11¨22¨33z m 00k k 0z 00m 0z k 2k k z 000m 0k k z 0z ⎡⎤⎢⎥-⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥+--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦(9) 定义主振型由于是无阻尼系统,因此系统守恒,系统存在振动主振型。
主振型意味着各物理坐标振动的相位角不是同相(相差0o )就是反相位(相差180o ),即同时达到平衡位置和最大位置。
主振型定义如下:()i i j ωt+i i sin ωt+=Im(e )φφi mi mi z =z z (10)其中z i 为第i 阶频率下,各自有度的位移矢量,z mi 为第i 个特征矢量,表示第i 阶固有频率下的振型,i ω为第i 阶频率下的第i 个特征值,i φ为初始相位。
对于三自由度系统,在第i 阶频率下,等式可以写成1m1i 2m2i i i 3m3i z z z =z sin(ωt+)z z φ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(11)mki z 表示第k 个自由度在第i 阶模态下的模态矩阵。
特征值对式(10)二次求导,得2i i i =-ωsin(ω+)φ¨i mi z z (12)代入齐次运动方程得m [−ωi 2z mi sin (ωi +i )]+k [z mi sin (ωi +i )]=0(13)去除sin (ωi +i )项化简得 (k −ωi 2m )z mi =0(14)以矩阵的形式展开得:2i 2i mi 2i k-ωm -k 0-k 2k-ωm -k z =00-k k-ωm ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(15) z mi 有非零解,则2i 2i 2i k-ωm -k 0-k 2k-ωm -k =00-k k-ωm ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦(16)即()234222ω-m ω+4km ω-3k m =0(17)方程解如下:1ω=0,2ω=±,3ω=±。
三个解对应该系统的前三阶固有频率,每一个特征根对应一个特征矢量z i ,表示对应模态下该系统的振型。
特征矢量由式 (k −ωi 2m )z mi =0得矩阵展开形式:2i m1i 2i m2i 2i m3i k-ωm -k 0z -k 2k-ωm -k z =00-k k-ωm z ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(18) 展开第一行和第二行,忽略下脚标m 和i ,得()()2i1221i3k-ωm z -kz =0-kz 2k-ωm kz+-=(19)得22i 124223ii21z k-ωm =z k z m ω-3km ω+k =z k (20)如果设定了1z 值,则就可以求出三个特征根值下,2z 和3z 相对于1z 的位移。
假设m=k=1,一阶模态,1ω=0:21z =1z ,31z =1z ,即z 1=[111];二阶模态,223k ω=m :21z =0z ,31z =-1z ,即z 2=[10−1];三阶模态,23kω=m :21z =-2z ,31z =1z ,即z 3=[1−21]。
模态矩阵所谓模态矩阵就是指各列由各阶模态特征矢量构成的矩阵,如图4所示。
图模态矩阵对于前面提到的三自由度系统,模态矩阵如下:z m =[11110−21−11]运动方程的解耦对于一个复杂的系统,在物理坐标系统中建立的运动方程之间存在耦合关系,因此求解起来比较麻烦,因此需要进行坐标系转化,将耦合的运动方程变为非耦合的运动方程,再将求得的结果转化为物理坐标系下的结果,运动方程解耦过程如下图5:图运动方程解耦过程在进行坐标变换之前需对刚度矩阵和质量矩阵进行归一化。
任意上面的三自由度系统为例,由式m [−ωi 2z mi sin (ωi +i )]+k [z mi sin (ωi +i )]=0得ωi 2mz mi =kz mi (21)ωj 2mz mj =kz mj (22)对式(21)左乘z mj T 得ωi 2z mj T mz mi =z mj Tkz mi (23)又因为ωj 2z mj T m T =z mj T k T因为系统对称所以,m T =m ,k T =k ,则:ωj 2z mj T m =z mj Tk (24)对式(24)右乘z miωj 2z mj T mz mi =z mj Tkz mi (25)则式(23)—式(25)得(ωi 2−ωj 2)z mj T mz mi =0(26)当(ωi 2−ωj 2)≠0时,则z mj Tmz mi =m ji =0(27)当(ωi 2−ωj 2)=0,即i =j ,则z mj T mz mi 可以为任何值,令z mj Tmz mi =m ii (28)则对质量矩阵和刚度矩阵的归一化结果如下:m n =z m T mz m (29)k n =z m Tkz m (30)特征矢量的归一化由于特征矢量只是位移之比,而不是绝对振幅,因此可以对其进行归一化处理。
令z ni T mz ni =1.0,其中z ni =z mi[z mi T mz mi ]12=z mi q i(31)q i =[∑z mji (∑m jk z mki n k=1)n j=1]12(32)对于对角质量矩阵q i =[∑m k z mki 2n j=1]12(33)则三自由度系统:z m =[11110−21−11](34)=00n z (35) 则归一化的质量矩阵为100010001⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦Tn n n m =z mz (36) 同理归一化后的刚度矩阵为000k =010m003⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦n k (37) 可以看出归一化后的刚度矩阵对角线上的各项就是各阶模态固有频率的平方。
运动方程解耦将物理坐标系下的运动方程¨11¨22¨33z m 00k -k 0z 0 0m 0z +-k 2k -k z =000m 0-k k z 0z ⎡⎤⎢⎥⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦⎢⎥⎣⎦按照前面介绍的归一化方法转化为主坐标系下的运动方程,其结果如下:¨p1p1¨p2p2¨p3p30z 00z 0k 00z +-k z =0m 00z 03k z 0-km 001101⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦38) 可以看出在主坐标系中的运动方程之间没有耦合关系,分别单独描述各阶模态的运动特性。
初始条件和激励的坐标转换物理坐标系中的非齐次运动方程为..mz+kz =F (39)做如下变形..T-1T -1Tn n nn n n n z mz z z+z kz z z =z F (40)其中T nn z mz ,Tn n z kz 就是前面介绍的质量和刚度矩阵的对角化。
令Tp n n m =z mz ,主坐标质量矩阵; T p n n k =z kz ,主坐标刚度矩阵; ....-1p nz z =z ,主坐标系加速度矢量;-1n p z z =z ,主坐标系位移矢量; T n p z F =F ,主坐标系激励矢量。
同样的关系也适用于初始位移和速度:-1op n o ..-1op noz =z z z =z z (42)两种坐标系的对比动模态下的位移和速度。
由主坐标系转变为物理坐标系前面介绍了物理坐标系与主坐标系之间的关系为-1n p z z =z (43)对式(41)左乘n z ,变为=-1n n n p z z z =z z z (44)同理p =..n z z z (45)非参数模型传递函数传递函数由系统的本质特性所决定,与系统的输入输出无关。
知道了系统的传递函数就可以根据输入求输出或根据输出求输入。
以图2的单自由度粘性阻尼系统为例,图单自由度系统则该系统的运动方程为:...m z +c z +kz=F (1)其中m 为质量,c 为阻尼系数,k 为刚度系数,z ,ż,z̈分别为位移、速度和加速度。
对二阶微分方程进行拉普拉斯变换,其中二阶导数项的拉普拉斯变换为:ℒ{z̈(t)}=s 2z (s )−sz (0)−ż(0)(2)假设初始位移和速度都为零,则ℒ{z̈(t)}=s 2z (s )(3)则经过拉普拉斯变换后的运动方程为:ms 2z (s )+csz (s )+kz (s )=F(s)(4)求解拉氏方程得传递函数:22z(s)11/m==c kF(s)ms +cs+k s +s+m m(5)其中定义2n kω=m为非阻尼系统的固有频率,rad/sec;cr c =界阻尼值,ζ为阻尼比,一般为阻尼与临界阻尼的比值,crc=c ζ,则n c 2ω=mζ。
则传递函数又可以写成:22n n z(s)1/m=F(s)s +2ωs+ωζ(6) 频响函数FRF用“j ω”代替s ,得系统的频响函数,其中j 是虚数项:()()22n n 22n n z(j ω)1/m=F(j ω)j ω+2ζωj ω+ω1/m=-ω+2ζωωj+ω(7)其中n kω=m,ζ,则频响函数可以写成 2z(j ω)1=F(j ω)-m ω+j ωc+k(8) 质量、阻尼、刚度对FRF 的影响刚度增大导致共振频率的增大,并且降低FRF 在低频段的幅值。
增加阻尼会使共振频率略微减小,但它的主要作用是减小频响函数在共振点的幅值,同时使相位的改变较为平缓。
如果阻尼为零,在共振点振动振幅将趋于无穷大,相位会突变180o 。
增大质量会降低共振频率,同时也降低FRF 在高频段的幅值。