文档库 最新最全的文档下载
当前位置:文档库 › 什么是模态分析

什么是模态分析

什么是模态分析
什么是模态分析

什么是模态分析

你不是第一个要求我用通俗易懂的语言解释模态分析的人,这样一来,任何人都能明白模态分析到底是怎样一个过程。简单地说,模态分析是根据用结构的固有特征,包括频率、阻尼和模态振型,这些动力学属性去描述结构的过程。那只是一句总结性的语言,现在让我来解释模态分析到底是怎样的一个过程。不涉及太多的技术方面的知识,我经常用一块平板的振动模式来简单地解释模态分析。这个解释过程对于那些振动和模态分析的新手们通常是有用的。

考虑自由支撑的平板,在平板的一角施加一个常力,由静力学可知,一个静态力会引起平板的某种静态变形。但是在这儿我要施加的是一个以正弦方式变化,且频率固定的振荡常力。改变此力的振动频率,但是力的峰值保持不变,仅仅是改变力的振动频率。同时在平板另一个角点安装一个加速度传感器,测量由此激励力引起的平板响应。

现在如果我们测量平板的响应,会注意到平板的响应幅值随着激励力的振动频率的变化而变化。随着时间的推进,响应幅值在不同的频率处有增也有减。这似乎很怪异,因为我们对此系统仅施加了一个常力,而响应幅值的变化却依赖于激励力的振动频率。具体体现在,

当我们施加的激励力的振动频率越来越接近系统的固有频率(或者共振频率)时,响应幅值会越来越大,在激励力的振动频率等于系统的共振频率时达到最大值。想想看,真令人大为惊奇,因为施加的外力峰值始终相同,而仅仅是改变其振动频率。

时域数据提供了非常有用的信息,但是如果用快速傅立叶变换(FFT)将时域数据转换到频域,可以计算出所谓的频响函数(FRF)。这个函数有一些非常有趣的信息值得关注:注意到频响函数的峰值出现在系统的共振频率处,注意到频响函数的这些峰出现在观测到的时域响应信号的幅值达到最大时刻的频率处。

如果我们将频响函数叠加在时域波形之上,会发现时域波形幅值达到最大值时的激励力振动频率等于频响函数峰值处的频率。因此可以看出,既可以使用时域信号确定系统的固有频率,也可以使用频响函数确定这些固有频率。显然,频响函数更易于估计系统的固有频率。

许多人惊奇结构怎么会有这些固有特征,而更让人惊奇的是在不同的固有频率处,结构呈现的变形模式也不同,且这些变形模式依赖于激励力的频率。

现在让我们了解结构在每一个固有频率处的变形模式。在平板上均匀分布45个加速度计,用于测量平板在不同激励频率下的响应幅值。如果激励力在结构的每一个固有频率处驻留,会发现结构本身存在特定的变形模式。这个特征表明激励频率与系统的某一阶固有频率相等时,会导致结构产生相应的变形模式。我们注意到当激励频率在第一阶固有频率处驻留时,平板发生了第1阶弯曲变形,在图中用蓝色表示。在第2阶固有频率处驻留时,平板发生了第1阶扭转变形,在图中用红色表示。分别在结构的第3和第4阶固有频率处驻留时,平板发生了第2阶弯曲变形,在图中用绿色表示,和第2阶扭转变形,在图中用红紫红色表示。这些变形模式称为结构的模态振型。(从纯数学角度讲,这种叫法实际上不完全正确,但在这儿作为简单的讨论,从实际应用角度讲,这些变形模式非常接近模态振型。)

我们设计的所有结构都具有各自的固有频率和模态振型。本质上,这些特性取决于确定结构固有频率和模态振型的结构质量和刚度分布。作为一名设计工程师,需要识别这些频率,并且当有外力激励结构时,应知道它们怎样影响结构的响应。理解模态振型和结构怎样振动有助于设计工程师设计更优的结构。模态分析有太多的需要讲解的地方,但这个例子仅仅是一个非常简单的解释。

现在我们能更好地理解模态分析主要是研究结构的固有特性。理解固有频率和模态振型(依赖结构的质量和刚度分布)有助于设计噪声和振动应用方面的结构系统。我们使用模态

分析有助于设计所有类型的结构,包括机车、航天器,宇宙飞船、计算机、网球拍、高尔夫球杆……这些清单举不胜举。

我希望这次简明的介绍有助于解释什么是模态分析。我用上面的例子向我母亲解释模态分析,她第一次真正明白了我到底在做什么。从此以后,她一直用一系列非常像模态分析的词语向她的朋友讲解模态分析,而她称这种分析为傻瓜式的分析……当然,这又是另一个故事了。

模态分析意义

模态分析意义模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的,则称为计算模态分析;如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动千姿百态、瞬息变化。模态分析提供了研究各种实际结构振动的一条有效途径。首先,将结构物在静止状态下进行人为激振,通过测量激振力与胯动响应并进行双通道快速傅里叶变换(FFT)分析,得到任意两点之间的机械导纳函数(传递函数)。用模态分析理论通过对试验导纳函数的曲线拟合,识别出结构物的模态参数,从而建立起结构物的模态模型。根据模态叠加原理,在已知各种载荷时间历程的情况下,就可以预言结构物的实际振动的响应历程或响应谱。近十多年来,由于计算机技术、

FFT 分析仪、高速数据采集系统以及振动传感器、激励器等技术的发展,试验模态分析得到了很快的发展,受到了机械、电力、建筑、水利、航空、航天等许多产业部门的高度重视。已有多种档次、各种原理的模态分析硬件与软件问世。在各种各样的模态分析方法中,大致均可分为四个基本过程:(1)动态数据的采集及频响函数或脉冲响应函数分析1)激励方法。试验模态分析是人为地对结构物施加一定动态激励,采集各点的振动响应信号及激振力信号,根据力及响应信号,用各种参数识别方法获取模态参数。激励方法不同,相应识别方法也不同。目前主要由单输入单输出(SISO)、单输入多输出(SIMO)多输入多输出(MIMO)三种方法。以输入力的信号特征还可分为正弦慢扫描、正弦快扫描、稳态随机(包括白噪声、宽带噪声或伪随机)、瞬态激励(包括随机脉冲激励)等。2)数据采集。SISO 方法要求同时高速采集输入与输出两个点的信号,用不断移动激励点位置或响应点位置的办法取得振形数据。SIMO 及MIMO 的方法则要求大量通道数据的高速并行采集,因此要求大量的振动测量传感器或激振器,试验成本较高。3)时域或频域信号处理。例如谱分析、传递函数估计、脉冲响应测量以及滤波、相关分析等。(2)建立结构数学模型根据已知条件,建立一种描述结构状态及特性的模型,作为计算及识别参数依据。目前一般假定系统为线性的。由于采用的识别方法不同,也分为频域建模和时

模态分析和频率响应分析的目的

有限元分析类型 一、nastran中的分析种类 (1)静力分析 静力分析是工程结构设计人员使用最为频繁的分析手段,主要用来求解结构在与时间无关或时间作用效果可忽略的静力载荷(如集中载荷、分布载荷、温度载荷、强制位移、惯性载荷等)作用下的响应、得出所需的节点位移、节点力、约束反力、单元内力、单元应力、应变能等。该分析同时还提供结构的重量和重心数据。 (2)屈曲分析 屈曲分析主要用于研究结构在特定载荷下的稳定性以及确定结构失稳的临界载荷,NX Nastran中的屈曲分析包括两类:线性屈曲分析和非线性屈曲分析。 (3)动力学分析 NX Nastran在结构动力学分析中有非常多的技术特点,具有其他有限元分析软件所无法比拟的强大分析功能。结构动力分析不同于静力分析,常用来确定时变载荷对整个结构或部件的影响,同时还要考虑阻尼及惯性效应的作用。 NX Nastran的主要动力学分析功能:如特征模态分析、直接复特征值分析、直接瞬态响应分析、模态瞬态响应分析、响应谱分析、模态复特征值分析、直接频率响应分析、模态频率响应分析、非线性瞬态分析、模态综合、动力灵敏度分析等可简述如下: ?正则模态分析 正则模态分析用于求解结构的固有频率和相应的振动模态,计算广义质量,正则化模态节点位移,约束力和正则化的单元力及应力,并可同时考虑刚体模态。 ?复特征值分析 复特征值分析主要用于求解具有阻尼效应的结构特征值和振型,分析过程与实特征值分析类似。此外

Nastran的复特征值计算还可考虑阻尼、质量及刚度矩阵的非对称性。 ?瞬态响应分析(时间-历程分析) 瞬态响应分析在时域内计算结构在随时间变化的载荷作用下的动力响应,分为直接瞬态响应分析和模态瞬态响应分析。两种方法均可考虑刚体位移作用。 直接瞬态响应分析 该分析给出一个结构随时间变化的载荷的响应。结构可以同时具有粘性阻尼和结构阻尼。该分析在节点自由度上直接形成耦合的微分方程并对这些方程进行数值积分,直接瞬态响应分析求出随时间变化的位移、速度、加速度和约束力以及单元应力。 模态瞬态响应分析 在此分析中,直接瞬态响应问题用上面所述的模态分析进行相同的变换,对问题的规模进行压缩,再对压缩了的方程进行数值积分,从而得出与用直接瞬态响应分析类型相同的输出结果。 ?随机振动分析 该分析考虑结构在某种统计规律分布的载荷作用下的随机响应。例如地震波,海洋波,飞机超过建筑物的气压波动,以及火箭和喷气发动机的噪音激励,通常人们只能得到按概率分布的函数,如功率谱密度(PSD)函数,激励的大小在任何时刻都不能明确给出,在这种载荷作用下结构的响应就需要用随机振动分析来计算结构的响应。NX Nastran中的PSD可输入自身或交叉谱密度,分别表示单个或多个时间历程的交叉作用的频谱特性。计算出响应功率谱密度、自相关函数及响应的RMS值等。计算过程中,NX Nastran不仅可以像其他有限元分析那样利用已知谱,而且还可自行生成用户所需的谱。 ?响应谱分析 响应谱分析(有时称为冲击谱分析)提供了一个有别于瞬态响应的分析功能,在分析中结构的激励用各个小的分量来表示,结构对于这些分量的响应则是这个结构每个模态的最大响应的组合。 ?频率响应分析 频率响应分析主要用于计算结构在周期振荡载荷作用下对每一个计算频率的动响应。计算结果分实部和虚部两部分。实部代表响应的幅度,虚部代表响应的相角。 直接频率响应分析 直接频率响应通过求解整个模型的阻尼耦合方程,得出各频率对于外载荷的响应。该类分析在频域中主要求解两类问题。第一类是求结构在一个稳定的周期性正弦外力谱的作用下的响应。结构可以具有粘性阻尼和结构阻尼,分析得到复位移、速度、加速度、约束力、单元力和单元应力。这些量可以进行正则化以获得传递函数。 第二类是求解结构在一个稳态随机载荷作用下的响应。此载荷由它的互功率谱密度定义。而结构载荷由上面所提到的传递函数来表征。分析得出位移、加速度、约束力或单元应力的自相关系数。该分析也对自功率谱进行积分而获得响应的均方根值。 模态频率响应 模态频率响应分析和随机响应分析在频域中解决的两类问题与直接频率响应分析解决相同的问题。

模态试验分析系统

模态试验分析系统 系统简介 模态试验与分析系统是指通过数据采集系统获得激励(和响应)数据,经动态信号分析与模态参数识别,确定机械结构的固有频率、阻尼比、振型和模态参与因子等揭示结构动态特性的参数。模态实验广泛应用于振动排故、状态检测、故障诊断和结构健康监测,以及动态响应预报、结构动态修改、有限元模型修正、动态分析与设计、振动控制等。 系统特点 ★快速几何建模 1、集成交互式几何建模模块,实现节点、连线、多边形、3D对象的交互式选择、移动、旋转、放大、删除、修改等功能 2、可定义总体坐标和局部坐标,具有笛卡尔、柱、以及球等三种坐标系统,各种坐标系统间转换方便 3、可实现线段、直线、矩形、梯形、扇面、椭圆、圆台、球体等规则3D对象的快速建模,还可自

定义三维单元库 4、除了交互式几何建模,模型几何信息也可通过配置信息界面直接进行修改、添加、删除等操作 ★快速、易用的信号分析功能 1、向导式的信号处理参数设置,实现趋势去除、时域抽取、快速傅立叶变换(FFT)、加窗函数等功能 2、 FFT长度:基2整数,根据实测数据自由可选;重叠:0%~83%,可从下拉列表中选择;平均次数:用户自定义;窗函数:矩形窗、汉窗、海明窗、平顶窗、指数窗、力窗、指数窗等;分析频率范围:采样频率的1/2或1/2.56 3、功率谱估计:自谱、互谱、功率谱密度矩阵、半功率谱密度矩阵 4、单输入多输出(SIMO)的频率响应函数(FRF)估计:H1、H2估计 5、多输入多输出(MIMO)的频率响应函数估计及相干函数估计 6、多线程支持的信号处理过程,并可采用不同设置参数重复进行 ★灵活的二维\三维图形显示、控制和输出 1、提供专用的二维曲线与三维图形控制面板,以及鼠标、快捷键、菜单等多种控制方式 2、多种曲线表达方式,诸如频率响应函数的幅值(线性、对数、dB坐标)、相位、展开相位、实部、虚部、奈奎斯特图等 3、方便灵活的二维曲线显示与控制,网格、图例等元素可显示或隐藏,并能提供相应曲线的完善测量信息(测量节点、方向,是否原点测量等) 4、缩放(具有不同缩放状态的记忆能力)、选段、寻峰寻谷等实用功能 5、方便灵活的三维图形显示与控制,节点号、输入/输出标记、坐标轴等元素可显示或隐藏,并能轻易实现平移、缩放、旋转等功能 6、提供三维图形的俯仰、左右、前后等各向视图,能实现结构的框架线显示或着色面渲染 7、二维曲线和三维图形的各元素颜色均可自定义 8、基于OpenGL的三维图形动画控制,实现播放、暂停、帧播放、幅度控制、速度控制等功能 9、各种二维曲线和三维图形均可复制到操作系统剪贴板中,亦可一键存储为BMP或JPG文件 10、振型动画和ODS可直接输出成AVI文件 ★先进、准确、可靠的模态分析技术 1、EMA : 基于输入(激振力)、输出(响应)测量的试验模态分析技术 (1)单输入/多输出(SIMO)的全局模态识别技术,可识别得到全局模态参数 (2)多点激振的多输入/多输出(MIMO)模态识别技术,具有识别高密度或重频模态的能力,是大型、复杂结构试验模态分析的理想方法 (3)单参考点和多参考点锤击法(MRIT)模态识别技术。 2、OMA: 环境激励下仅有输出(响应)可测量的运行模态分析技术,可以对桥梁、建筑、汽车、飞机、旋转机械等机械结构在运行状态进行试验与分析,无须人工激振,只需测量响应 (1)不仅简单可行,同时还可获得结构在真实运行状态下的动态特性,且天然具备多参考点特性,具有解耦密集模态的能力 (2)基于全功率谱密度矩阵的窄带模态参数识别方法(频域空间域分解法,FSDD),方便易用,结

模态试验及分析的基本步骤

模态试验及分析的基本步骤 1.动态数据的采集及响应函数分析 首先应选取适当的激励方式。激励方式可以是正弦、随机或瞬态中的任何一种。激励方式不同,相应的模态参数识别方法也不同。目前主要有单输入单输出、单输入多输出和多输入多输出三种方法。然后进行数据采集。对于单输入单输出方法要求同时高速采集输入与输出两个点的信号,用不断移动激励点位置或响应点位置的办法取得振型数据;单输入多输出及多输入多输出的方法要求大量通道数据的高速采集,因此要求大量的振动测量传感器或激振器,试验成本极高。在采集信号数据以后,还要在时域或频域对信号进行处理,例如谱分析、传递函数估计、脉冲响应测量以及滤波、相关分析等。 2.建立结构数学模型 根据己知条件,建立一种描述结构状态及特性的模型,作为计算及参数识别的依据,目前一般假定系统为线性的。由于采用的识别方法不同,数学建模可分为频域建模和时域建模。根据阻尼特性及频率藕合程度又可分为实模态和复模态等。 3.参数识别 按识别域的不同可分为频域法、时域法和混合域法。激励方式不同,相应的识别参数方法也不尽相同。并非越复杂的方法识别的结果越可靠。对于目前能够进行的大多数不是十分复杂的结构,只要取得了可靠的频响数据,用简单的识别方法也可能获得良好的模态参数;反之,即使用最复杂的数学模型、最高级的拟合方法,如果频响测量数据不可靠,识别的结果也不会理想。 4.振型动画 参数识别的结果得到了结构的模态参数模型,即一组固有频率、模态阻尼以及相应各阶模态的振型。但是由于结构复杂,由许多自由度组成的振型的数组难以引起对振动直观的想象,所以必须采用振型动画的办法,将放大的振型叠加到原始的几何形状上。

ANSYS中的模态分析与谐响应分析

ANSYS中的模态分析与谐响应分析 作者:未知时间:2010-4-15 8:59:49 模态分析是分析结构的动力特性,与结构受什么样的荷载没有关系,只要给定了质量、弹性模量、泊松比等材料参数,并施加了边界约束就可以得到此状态下的各阶自振频率和振型(也称为模态)。 谐响应分析是分析结构在不同频率的简谐荷载作用下的动力响应,是与结构所受荷载相关的,只是结构所受荷载的都是简谐荷载,而且荷载频率的变化范围在谐响应分析时要给出来。 比如,在ANSYS谐响应分析中要给出这样的语句 FK,3,FX,7071,7071 !指定点荷载的实部和虚部(或者幅值和相位角) HARFRQ,0,2.5, !指定荷载频率的变化范围,也就是说只分析结构所受频率从0到2.5HZ之间的荷载 NSUBST,100, !指定频率从0到2.5之间分100步进行计算 这样,结构所受的这个点荷载的表达式实际上是 F=(7071+i*7071)*exp(i*omiga*t) !式中omiga从0到2.5*2*3.1415926变化 分析得到结果是各点物理量随频率变化的,但物理量的值一般为复数,包括实部的虚部,这可以从后处理LIST结点值看出来。 个人认为进行谐响应分析并不一定要先进行模态分析(也叫振型分析、振型分解等),而直接进行谐响应分析后查看结构的物理量随频率变化曲线时也会看到在结构的自振频率处响应会放大(共振)。如果已经进行过模态分析的话,会发现谐响应分析时的共振频率和模态分析提到的自振频率是一致的。但有些时候模态分析中得到的有些频率在谐响应分析的频响曲线里可能很不明显。因此,只能说在谐响应分析前进行一下模态分析可以对结构的自振特性有个了解,以便验证谐响应分析结果是否合理。 另外,谐响应分析应该是频域分析方法的一个部分。对于相地震那样的时间过程线,直接进行时域分析(ANSYS里用暂态分析)可得到结构随时间的响应。而如果进行频域分析,就应该通过傅立叶变换把时域地震曲线变为由多个简谐荷载的叠加,然后再以此简谐荷载做为谐响应分析时的荷载进行谐响应分析,最后再对谐响应分析得到的结果进行傅立叶逆变换得到时域的结果。不知道这种理解是否正确,我也没有用ANSYS这样做过。如果正确的话,时域分析和频域分析的结果应该是一致的。 模态分析的应用及它的试验模态分析 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模

结构模态分析方法

模态分析技术的发展现状综述 摘要:本文首先系统的介绍了模态分析的定义,并以模态分析技术的理论为基础,查阅了大量的文献和资料后,介绍了三种模态分析技术在各领域的应用,以及国内外对于结构模态分析技术研究的发展现状,分析并总结三种模态分析技术的特点与发展前景。 关键词:模态分析技术发展现状 Modality Analysis Technology Development Present Situation Summary Abstract:This article first systematic introduction the definition of modality analysis,and based on modal analysis theory,after has consulted the massive literature and the material.Introduced application about three kind of modality analysis technology in various domains. At home and abroad, the structural modal analysis technology research and development status quo.Analyzes and summarizes three kind of modality analysis technology characteristic and the prospects for development. Key words:Modality analysis Technology Development status 0 引言 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。模态分析的过程如果是由有限元计算的方法完成的,则称为计算模态分析;如果是通过试验将采集的系统输入与输出信号经过参数识别来获得模态参数的,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备故障诊断的重要方法。 1 数值模态分析的发展现状 数值模态分析主要采用有限元法,它是将弹性结构离散化为有限数量的具体质量、弹性特性单元后,在计算机上作数学运算的理论计算方法。它的优点是可以在结构设计之初,根据有限元分析结果,便预知产品的动态性能,可以在产品试制出来之前预估振动、噪声的强度和其他动态问题,并可改变结构形状以消除或抑制这些问题。只要能够正确显示出包含边界条件在内的机械振动模型,就可以通过计算机改变机械尺寸的形状细节。有限元法的不足是计算繁杂,耗资费时。这种方法,除要求计算者有熟练的技巧与经验外,有些参数(如阻尼、结合面特征等)目前尚无法定值,并且利用有限元法计算得到的结果,只能是一个近似值。 正因如此,大多数数学模拟的结构,在试制阶段常应做全尺寸样机的动态试验,以验证计算的可靠程度并补充理论计算的不足,特别对一些重要的或涉及人身安全的结构,就更是如此。 70 年代以来,由于数字计算机的广泛应用、数字信号处理技术以及系统辨识方法的发展 , 使结构模态试验技术和模态参数辨识方法有了较大进展,所获得的数据将促进产品性能的改进、更新[1] 。在硬件上,国外许多厂家研制成功各种类型的以FFT和

模态分析软件操作

模态分析软件操作说 明及实例 东方振动和噪声技术研究所 1999.3.16 目录 一模态分析的步骤 (2) 1.确定分析的方法 (2) 2.测点的选取、传感器的布置 (2) 3.仪器连接 (3) 4.示波 (3) 5.输入标定值 (3) 6.采样 (4) 7.传递函数分析 (4) 8.进行模态分析 (4) 二模态分析实例 (5)

例一自由梁的模态分析实例 (5) 例二楼房的模态分析实例 (15) 模态分析是一种参数识别的方法,因为模态分析法是在承认实际结构可以运用所谓“模态模型”来描述其动态响应的条件下,通过实验数据的处理和分析,寻求其“模态参数”。 模态分析的关键在于得到振动系统的特征向量(或称特征振型、模态振型)。试验模态分析便是通过试验采集系统的输入输出信号,经过参数识别获得模态参数。具体做法是:首先将结构物在静止状态下进行人为激振(或者环境激励),通过测量激振力与振动响应,找出激励点与各测点之间的“传递函数”,建立传递函数矩阵,用模态分析理论通过对试验导纳函数的曲线拟合,识别出结构的模态参数,从而建立起结构物的模态模型。 东方所研制的模态分析系统,自推出以来参与了许多重大的科研项目如大型航空航天设备(长征火箭、通信卫星、大型雷达、火箭发射平台等)、大桥(火箭激振钱塘江大桥、锤击法激振乌海黄河铁路大桥属国内首次)、大楼、大坝、、机车(汽车)车辆和大型港口机械等,分析精度高、操作简便,尤其是变时基模态分析及高速模态三视图动画技术更是在国内外处于领先地步。 一、模态分析的步骤 1. 确定分析的方法 DASP中提供的模态分析方法有多输入单输出法、单输入多输出法和多输入多输出方法。一般采用较多的是多输入单输出或单输入多输出方法,在这两种方法中选取时,视哪一种方法简便而定,如激励装置大、不好移动但传感器移动方便就选取单输入多输出方法(即单点激励、多点移步拾振);如传感器移动不方便但激励装置小、容易移动就选取多输入单输出方法(即单点拾振、多点移步激励)。 有时结构因为过于巨大和笨重,以至于采用单点激振时不能提供足够的能量,将我们所感兴趣的模态都激励出来;其次,结构在同一频率时可能有多个模态,这样单点激振就不能把它们分离出来,这时就要采取两个甚至多个激励来激发结构的振动,即采取多输入多输出方法。 在DASP中进行模态分析时,由于采用了高弹性聚能力锤和先进的变时基传递函数分析技术,对于象大型铁路桥、火箭发射平台这样的大型结构用力锤敲击就能分析出结构的模态;对于大型的混凝土结构(如大楼)可以以天然脉动作为激励信号进行模态分析。所以在大多数情况下,采取单输入多输出或多输入单输出方法就可完全满足工程需要。 2. 测点的选取、传感器的布置 选择好分析方法后,就要根据结构的特点和试验目的确定测点的数目和布置,以及传感器的安装方法等。

各种模态分析方法总结与比较

各种模态分析方法总结与比较 一、模态分析 模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。 模态分析的理论经典定义:将线性定常系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵,其每列为模态振型。 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的,则称为计算模记分析;如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过 AHA12GAGGAGAGGAFFFFAFAF

模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。 模态分析最终目标是在识别出系统的模态参数,为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 AHA12GAGGAGAGGAFFFFAFAF

AHA12GAGGAGAGGAFFFFAFAF 二、各模态分析方法的总结 (一)单自由度法 一般来说,一个系统的动态响应是它的若干阶模态振型的叠加。但是如果假定在给定的频带内只有一个模态是重要的,那么该模态的参数可以单独确定。以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。在给定的频带范围内,结构的动态特性的时域表达表示近似为: ()[]}{}{T R R t r Q e t h r ψψλ= 2-1 而频域表示则近似为: ()[]}}{ {()[]2ωλωψψωLR UR j Q j h r t r r r -+-= 2-2 单自由度系统是一种很快速的方法,几乎不需要什么计算时间和计算机内存。 这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正确的。然而实际情况通常并不是这样的,所以就需要用包含若干模态的模型对测得的数据进行近似,同时识别这些参数的模态,就是所谓的多自由度(MDOF)法。 单自由度算法运算速度很快,几乎不需要什么计算和计

模态分析与振动测试技术

模态分析与振动测试技术 固体力学 S0902015 李鹏飞

模态分析与振动测试技术 模态分析的理论基础是在机械阻抗与导纳的概念上发展起来的。近二十多年来,模态分析理论吸取了振动理论、信号分析、数据处理数理统计以及自动控制理论中的有关“营养”,结合自身内容的发展,形成了一套独特的理论,为模态分析及参数识别技术的发展奠定了理论基础。 一、单自由度模态分析 单自由度系统是最基本的振动系统。虽然实际结构均为多自由度系统,但单自由度系统的分析能揭示振动系统很多基本的特性。由于他简单,因此常常作为振动分析的基础。从单自由度系统的分析出发分析系统的频响函数,将使我们便于分析和深刻理解他的基本特性。对于线性的多自由度系统常常可以看成为许多单自由度系统特性的线性叠加。 二、多自由度系统模态分析 对于多自由度系统频响函数数学表达式有很多种,一般可以根据一个实际系统来讨论,给出一种形式;也可根据问题的要求来讨论,给出其他不同的形式。为了课程的紧凑,直接联系本课程的模态分析问题,我们就直接讨论多自由度系统通过频响函数表达形式的模态参数和模态分析。即多自由度系统模态参数与模态分析。 多自由度系统模态分析将主要用矩阵分析方法来进行。 我们以N个自由度的比例阻尼系统作为讨论的对象。然后将所分析的结果推广到其他阻尼形式的系统。 设所研究的系统为N个自由度的定常系统。其运动微分方程为: (2—1) ++= M X CX KX F ?)阶式中M,C,K分别为系统的质量、阻尼及刚度矩阵。均为(N N 矩阵。并且M及K矩阵为实系数对称矩阵,而其中质量矩阵M是正定矩阵,刚度矩阵K对于无刚体运动的约束系统是正定的;对于有刚体运动的自由系统则是半正定的。当阻尼为比例阻尼时,阻尼矩阵C为对称矩阵(上述是解耦条件)。 N?阶矩阵。即 X及F分别为系统的位移响应向量及激励力向量,均为1

模态分析与谐响应分析区别联系

模态分析是分析结构的动力特性,与结构受什么样的荷载没有关系,只要给定了质量、弹性模量、泊松比等材料参数,并施加了边界约束就可以得到此状态下的各阶自振频率和振型(也称为模态)。 谐响应分析是分析结构在不同频率的简谐荷载作用下的动力响应,是与结构所受荷载相关的,只是结构所受荷载的都是简谐荷载,而且荷载频率的变化范围在谐响应分析时要给出来。 比如,在ANSYS谐响应分析中要给出这样的语句 FK,3,FX,7071,7071 !指定点荷载的实部和虚部(或者幅值和相位角) HARFRQ,0,2.5, !指定荷载频率的变化范围,也就是说只分析结构所受频率从0到2.5HZ之间的荷载NSUBST,100, !指定频率从0到2.5之间分100步进行计算 这样,结构所受的这个点荷载的表达式实际上是 F=(7071+i*7071)*exp(i*omiga*t) !式中omiga从0到2.5*2*3.1415926变化 分析得到结果是各点物理量随频率变化的,但物理量的值一般为复数,包括实部的虚部,这可以从后处理LIST结点值看出来。 个人认为进行谐响应分析并不一定要先进行模态分析(也叫振型分析、振型分解等),而直接进行谐响应分析后查看结构的物理量随频率变化曲线时也会看到在结构的自振频率处响应会放大(共振)。如果已经进行过模态分析的话,会发现谐响应分析时的共振频率和模态分析提到的自振频率是一致的。但有些时候模态分析中得到的有些频率在谐响应分析的频响曲线里可能很不明显。因此,只能说在谐响应分析前进行一下模态分析可以对结构的自振特性有个了解,以便验证谐响应分析结果是否合理。 另外,谐响应分析应该是频域分析方法的一个部分。对于相地震那样的时间过程线,直接进行时域分析(ANSYS里用暂态分析)可得到结构随时间的响应。而如果进行频域分析,就应该通过傅立叶变换把时域地震曲线变为由多个简谐荷载的叠加,然后再以此简谐荷载做为谐响应分析时的荷载进行谐响应分析,最后再对谐响应分析得到的结果进行傅立叶逆变换得到时域的结果。不知道这种理解是否正确,我也没有用ANSYS这样做过。如果正确的话,时域分析和频域分析的结果应该是一致的。 模态分析的应用及它的试验模态分析 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的,则称为计算模记分析;如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。 模态分析最终目标在是识别出系统的模态参数,为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 模态分析技术的应用可归结为一下几个方面: 1) 评价现有结构系统的动态特性; 2) 在新产品设计中进行结构动态特性的预估和优化设计; 3) 诊断及预报结构系统的故障; 4) 控制结构的辐射噪声; 5) 识别结构系统的载荷。 机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动千姿百态、瞬息变化。模态分析提供了研究各种实际结构振动的一条有效途径。首先,将结构物在静止状态下进行人为激振,通过测量激振力与胯动响应并进行双通道快速傅里叶变换(FFT)分析,得到任意两点之间的机械导纳函数(传递函数)。用模态分析理论通过对试验导纳函数的曲线拟合,识别出结构物的模态参数,从而建立起结构物的模态模型。根据模态叠加原理,在已知各种载荷时间历程的情况下,就可以预言结构物的实际振动的响应历程或响应

最新模态试验及分析的基本步骤

模态试验及分析的基本步骤 1 1.动态数据的采集及响应函数分析 2 首先应选取适当的激励方式。激励方式可以是正弦、随机或瞬态中的任何一种。激3 励方式不同,相应的模态参数识别方法也不同。目前主要有单输入单输出、单输入多4 输出和多输入多输出三种方法。然后进行数据采集。对于单输入单输出方法要求同时5 高速采集输入与输出两个点的信号,用不断移动激励点位置或响应点位置的办法取得6 振型数据;单输入多输出及多输入多输出的方法要求大量通道数据的高速采集,因此要7 求大量的振动测量传感器或激振器,试验成本极高。在采集信号数据以后,还要在时8 域或频域对信号进行处理,例如谱分析、传递函数估计、脉冲响应测量以及滤波、相9 关分析等。 10 2.建立结构数学模型 11 根据己知条件,建立一种描述结构状态及特性的模型,作为计算及参数识别的依 12 据,目前一般假定系统为线性的。由于采用的识别方法不同,数学建模可分为频域建13 模和时域建模。根据阻尼特性及频率藕合程度又可分为实模态和复模态等。 14 3.参数识别 15 按识别域的不同可分为频域法、时域法和混合域法。激励方式不同,相应的识别参16 数方法也不尽相同。并非越复杂的方法识别的结果越可靠。对于目前能够进行的大多17 数不是十分复杂的结构,只要取得了可靠的频响数据,用简单的识别方法也可能获得18 良好的模态参数;反之,即使用最复杂的数学模型、最高级的拟合方法,如果频响测量19 数据不可靠,识别的结果也不会理想。 20 4.振型动画 21 参数识别的结果得到了结构的模态参数模型,即一组固有频率、模态阻尼以及相应22 各阶模态的振型。但是由于结构复杂,由许多自由度组成的振型的数组难以引起对振23

运用hypermesh进行模态分析

1 引言 系统的模态参数(模态频率、模态阻尼、振型)对系统的动态分析和优化设计具有实用价值。通常由试验模态分析和计算模态分析两种方法。但由于受实验条件和时间的限制,组织实施往往比较困难,而且在测量次数,测量数据的处理准确性方面也难以得到充分的保证,在设计阶段难以实现。基于虚拟样机技术的虚拟实验方法在履带车辆箱体类零部件模态参数测量方面在设计阶段就能为方案优化提供指导,缩短产品开发周期,节省费用。因此,开展在虚拟环境下测试箱体类零部件的模态参数研究与探讨并扩展其应用具有重要意义。本文以某型履带车辆传动箱设计为例,应用HyperMesh为前处理软件,对其进行了有限元网格的划分,进而对箱体的模态进行了分析。 2 箱体有限元模型的建立及模态分析 首先依据传动箱体的尺寸,建立箱体的三维实体模型。利用HyperMesh对传动箱体的实体模型进行有限元网格划分,箱体的材料为铝合金,其密度为 2.66e33kg/m3,泊松系数为0.31,杨氏模量为7.7e72N/m2,强度极限为176.4MPa。整个箱体共划分76151个4面体单元,22262个节点。在此过程中,还必须考虑到箱体有限元模型建立后与各传动轴之间的连接,即柔性体与刚体间的连接。传动箱各轴都是通过轴承与箱体连接的,笔者在有限元模型中应用多点约束(MPC,Multi-point Constraint)来模拟轴承的作用。所谓多点约束是将某节点的依赖自由度定义为其他若干节点独立自由度的函数。多点约束可以用于不相容单元间的载荷传递,表征一些特定的物理现象,比如刚性连接、铰接、滑动等。笔者在箱体有限元模型中各轴孔的中心点处建立一个虚拟杆单元,如图1所示。轴孔内表面各节点的自由度则依赖于对应的虚拟杆单元。各传动轴与箱体间的约束也是在对应的虚拟单元处建立,各传动轴上的作用力则通过相应的虚拟杆单元和多点约束作用于箱体之上。文中建立的包括轴承模型的传动箱箱体有限元模型如图2所示。

试验模态分析的两种方法

试验模态分析的两种方法 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。模态分析最终目标在是识别出系统的模态参数,为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。 试验模态分析主要有以下两种方法,OROS模态分析软件MODEL 2 完全具备了这两种常用的模态方 法。 锤击法模态测试 用于满足锤击法结构模态试验,以简明、直观的方法测量和处理输入力和响应数据,并显示结果。提供两种锤击方法:固定敲击点移动响应点和固定响应点移动敲击点。用力锤来激励结构,同时进行加速度和力信号的采集和处理,实时得到结构的传递函数矩阵。能够方便地设置测量参数,如触发量级、测量带宽和加窗类型,同时对最优的设置提供建议指导。 激振器法模态测试 主要是通过分析仪输出信号源来控制激振器,激励被测试件,输出信号有先进扫频正弦,随机噪声,正弦,调频脉冲等信号。支持单点激励(SIMO)与多点同时激励法(MIMO)。 1)几何建模 结构线架模型生成,节点数和部件数没有限制,测量点DOF自动加到通道标示;建立几何模型,以3维方式显示测量和分析结果。结构模型可以作为单个部件的装配,及采用不同的坐标系(直角、圆柱、球体坐标系),要求除点的定义外,还可定义线和面,真实的显示试验结构。结构线架模型生成,节点数和部件数没有限制,测量点自由度自动加到通道标示。

模态分析在工程中的应用概述

模态分析在工程中的应用概述 学号:XXXXXX 姓名:XXX 模态分析是研究结构动力特性的一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的,则称为计算模态分析(FEA);如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为实验模态分析(EMA)。通常,模态分析都是指实验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一个易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。 模态分析所寻求的最终目标在于改变机械结构系统由经验与类比和静态设计为动态、优化设计方法;在于借助试验与理论分析相结合的方法,对已有结构系统进行识别、分析和评价,从中找出结构系统在动态性能上所存在的问题,确保工程结构能安全可靠及有效地工作;在于根据现场测试的数据来这段及预报振动故障和进行噪声控制。通过这些方法为老产品的改进和新产品的设计提供可靠的依据。[1] 模态分析是一项综合性技术,可以应用于各个工程部门及各种工程结构。机器、建筑物、航天航空飞行器、船舶、汽车等的实际振动千姿百态、瞬息万变。模态分析提供了研究各种实际结构振动的一条有效途径。首先,将结构物在静止状态下进行人为激振,通过测量激振力与响应并进行双通道快速Fourier 变换(FFT)分析,得到任意两点之间的机械导纳函数(传递函数)。用模态分析理论通过对实验导纳函数的曲线拟合,识别出结构物体的模态参数,从而建立起结构物体的模态模型。根据模态叠加原理,在已知各种载荷时间历程的情况下,就可以预言结构物体的实际振动的响应历程或响应谱。[2] 模态分析技术的应用可以归纳为以下几个方面:评价现有结构系统的动态特性,在新产品设计中进行结构动态特性的预估及优化设计,诊断及预报机构系统的故障,控制结构的辐射噪声,识别结构系统的载荷。[1] 下面对近几年国内模态分析在工程中各个方面的应用分别进行概述。 1.评价现有结构系统的动态特性 在处理结构的振动问题时,必须对其动态特性有全面的了解,而其动态特性

模态分析理论应用实际的讨论

模态分析理论应用实际的讨论 模态分析在结构设计中的应用认识小结 在结构设计中,我们通常要运用模态分析的方法来辅助设计,提高结构设计的合理性和科学性。模态参数获取有两种方法: 一种是有限元法,一般的FEA软件都可以计算,WB当然也没有问题拉; 一种是测试的方法,比如用LMS https://www.wendangku.net/doc/df9905456.html,b来测试。 这两种方法对于测试简单的结构是没有问题的,分析结果和试验结果很吻合。但是对于复杂的装配体结构,FEA软件就显得无能为力了,因为装配体有令人讨厌的结合面,对于结合面的分析,据我所知目前还没有比较好的办法(就算是最高的CAE高手恐怕也算不准)。所以复杂装配体的模态一般用测试的方法解决。当然CAE工程师可以用实验数据得到的结合面刚度阻尼值来修正自己的有限元模型。 一般模态分析的结果中,最受关注的是固有频率值及其振型。固有频率主要用以对照结构外的激振频率,看是否出现共振,共振出现的后果很严重,它会使设备的加工精度降低很多,另外固有频率值是衡量结构动静刚度的标杆,如果我想提高结构的动静刚度,不断改变自己设计的结构一般就能实现,当然设计水平也很重要;而通过观察振型我可以判断这个振型是否影响我设备的加工精度,如果影响的话,我会考虑将改变这个振型的频率,避免实际生产中出现加工精度降低的情况。模态分析在CAE中应该很简单,算出固有频率和振型也很轻松。但是如何在设计中运用好这个工具其实有很多学问。对于振型而言,可能不同的领域关注的焦点可能会不一样。以机床为例,如果计算机床的床身模态振型,可能振型有弯曲,扭转等众多振型,如果存在机床进刀、加工方向的振型,那么有可能这些振型会影响机床的加工精度。那么在设计阶段就必须对结构进行调整,比如修改结构内部的肋板分布,提高影响加工精度振型的固有频率,减少发生共振进而影响机床加工精度的可能性。我的看法是,振型模态分析要和结构强度刚度分析结合在一起,强度分析结果的高应力区如果和某一阶模态振型位移较大区域重合,就可认为结构是偏危险的,这些高应力区域有可能就是疲劳裂纹的萌生位置,而实际中的连续结构体振型应该是无穷多的,经典理论认为实际工程中能够对结构安全产生影响的往往只是低阶的频率振型,所以只要结构避开低阶共振区就能安全运行,然而随着结构形式运行条件等因素的不断变化,现代机械的振动形式也越来越复杂,除了静态强度刚度,动态强度刚度也越来越重要,在水中的湿模态分析,目前似乎还没有完美简洁的解决办法,计算分析所采用的模型和计算条件与实际运行中结构之间的差异会直接影响计算结果的精度,所以如何减小这个差异,或者说如何使分析过程更加接近实际是一直以来我们的目标。 模态分析中经常遇到的问题就是当分析对象为装配体的时候。装配体模态计算的正确性绝不仅仅在熟悉产品这么简单,尤其是类似于螺栓结合面、导轨结合面的地方,关于结合面的研究老早就到了一个瓶颈了,由于结合部特性参数的影响因素众多,如结合面材料、加工方法和表面质量,结合面介质及其性质,结合面几何形状及法面压力大小等,特别是在结合部作用机理尚未被真正揭示之前,要在理论上精确获得结合部的特性参数及其分析计算表达式非常困难,故用有限元法识别精度还有待验证。 结合部动力学参数识别问题的确是个技术性难题。目前解决好这一问题的手段是:测试+仿真,建立混合模型。另外对于产品的认知度问题是个值得讨论的问题,比如加强劲板形状的设计就是个问题。你是否已经能够罗列出各种简单振动模式下最好的结构形式?首先列一张表,然后你会心里有数些。但产品并非那么简单,所以需要设计复杂结构。那么,仅仅凭借模态测试是不够的,需要做结构形式的优化,那我们现有的优化技术中,拓扑优化是解决这一问题的好帮手。 曾经拿一家公司的产品,测试和计算发现他们的产品第一阶模态就到了300Hz以上,而同形式的产品,国内仅能到70几Hz.这个差距是何等的大?想办法把我们的产品也做到这样,那你就牛了。 这里谈到结构优化,我就插一句,ANSYS Workbench在分析或者说验证方面很不错,但是要涉及到拓扑优化和形貌优化则比较差,几乎不能应用到实际工程中,最多使用的尺寸优化。如果大家要做结构优化的话,建议使用一下HyperWorks/Optistruct,这个在结构优化上可以说是绝对领先的。. 还有就是共振的实际分析

模态测试与分析报告基本概念

模态测试与分析基本概念 1.模态假设:线性假设、时不变假设、互易性假设、可观测性假设 线性假设:结构的动态特性是线性的,就是说任何输入组合引起的输出等于各自输出的组合,其动力学特性可以用一组线性二阶微分方程来描述。 时不变性假设:结构的动态特性不随时间变化,因而微分方程的系数是与时间无关的常数。 可观测性假设:这意味着用以确定我们所关心的系统动态特性所需要的全部数据都是可以测量的。 互易性假设:结构应该遵从Maxwell互易性原理,即在q点输入所引起的p点响应,等于在p点的相同输入所引起的q点响应。 2.EMA、OMA、ODS 试验模态分析(Experimental Modal Analysis, EMA) 力锤激励EMA技术 激振器激励EMA技术 工作模态分析(Operational Modal Analysis, OMA) 工作变形模态(Operational Deflection Shape, ODS) 3.SISO、SIMO、MIMO SISO:设置1个响应测点,力锤激励遍历所有测点,也称为SRIT SIMO:设置若干响应测点,力锤激励遍历所有测点,也称为MRIT;用一个激振器固定在某测点处激励结构,测量所有测量自由度的响应,经FFT快速测量计算FRF MIMO:用多个激振器激励结构,测量所有测量自由度的响应,经FFT快速测量计算MIMO-FRFs,输入能量均匀,数据一致性好,能分离密集和重根模态,在大型复杂或轴对称结构模态试验尤为重要 4.模态分析基本步骤 建立模型:确定测量自由度、生成几何、确定各类参数:BW,参考点、触发等 测量:FRF,(时域数据可选) 参数估计:曲线拟合、参数提取 验证:MAC、MOV、MP等

相关文档
相关文档 最新文档