动力响应理论.doc

第2章 动力响应理论2.1引言机柜结构动力响应的计算机仿真分析是以设备动力响应理论为基础的，是进行设备结构动力响应研究的一种有效手段。论文中主要研究设备动力响应两个方面的内容：设备结构固有特性分析和结构在地震波作用下的响应分析。固有特性分析可以得到结构的固有频率和固有振型，是进行响应分析的基础；地震波响应分析将得到设备响应的时间历程变化。在使用有限元工具对结构进行建模、分析之前必须掌握结构动力响应的理论和相关的有限元基本原理。因此，本章重点叙述了与设备结构动力响应相关的机械振动学理论及其有限元仿真技术。2.2结构动力响应分析相关理论2.2.1结构固有特性分析理论机柜设备结构的固有特性包括固有频率和振型，是响应分析的基础。通过进行结构的固有特性分析可以使设计有效地避开结构的共振频率。机柜设备是一个复杂振动系统，在理论分析过程中，常常可以把机柜设备简化为多自由度集中参数系统。一般，多自由度系统的自由振动方程可以写成如下形式： 式中：， 和分别为系统的质量矩阵、阻尼矩阵和刚度矩阵;、分别为系统的位移列向量、速度列向量和加速度列向量。而多自由系统的无阻尼自由振动方程可以写成如下形式: 通常系统的自由振动是简谐振动，所以可以假设式 的解为： 式中：为系统的振幅列向量;为系统的自由振动频率。将代入，就可以得到系统的振型方程，其具体形式如下： 可以看到，式是一个齐次线性方程组，根据线性代数知识，它具有非零解的充分必要条件为系数矩阵的行列式为零，亦即有下式成立。 式 称为系统的特征方程或频率方程，他是关于的次代数方程。解之，我们可以得到个解。具体如下： 若系统为正定系统，则有 式中：为系统的第阶固有圆频率。将分别代入系统的振型方程（2-3）中，可以解得与之对应的振幅列向量，也称为系统的主模态或主振型。注意：（1）系统固有频率和模态的数目与系统拥有的自由度数相同，n自由度系统必有n阶固有频率和模态（主振型）。（2）系统的固有频率仅与系统的质量矩阵和刚度矩阵有关，而与外界干扰力无关，它们是系统本身的固有性质。（3）系统的主振型表示系统以频率作自由振动时系统的各个自由度的振幅的相对比值。2.2.2结构冲击响应分析理论 通过系统的自由振动方程，可以解得系统得固有频率和主振型。在冲击激励作用下，多自由度系统振动将作受迫振动。利用激励响应理论可以解得系统受迫振动过程中系统的响应。2.2.2.1 多自由度系统的受迫振动方程多自由度系统的受迫振动方程可以写成如下形式： 式中：为系统受到的激励力向量，如果系统受到基础加速度激励则激励力向量就是该加速度系统的等效惯性力；其他符号的意义同式（2-1）。通常情况下，系统的各个自由度之间存在耦合，方程（2-7）中的、和不是对角阵，方程难以求解。为了得到系统的响应必须选择适当的坐标系将（2-7）式解耦。2.2.2.2振动方程的解耦（1） 正则振型与正则振型矩阵 在对振动方程解耦之前先介绍系统的正则振型与正则振型矩阵的概念。通过式（2-3）与式（2-4）已经求得系统的各阶固有频率和相应的主振型。令 使得 式中：为第i阶正则振型；为待定系数。将（2-8）带入（2-9），可以得到 将（2-10）带入（2-8）可得到各阶正则振型，进而可以写出正则振型矩阵，具体如下： （2）振动方程解耦 由于各阶主振型关于、是正交的，各阶正则振型也是关于、正交的，可以利用正则振型矩阵来对振动方程进行解耦。选取一组广义坐标又称为主坐标，令 将（2-12）代入（2-7）中，并且将方程两端左乘，则可以得到解耦后的振动方程： 式中，为单位阵；伪对角阵，其对角限元素为各阶固有频率的平方；为对角化的阻尼阵，采用比例阻尼假设，其对角线上的第i个元素为： 为外力对角阵，可进一步表示为： 式（2-13a）是一个非齐次线性方程组，其中第i个方程可以写成： 式中，为第i个主坐标。2.2.2.3单自由度系统对任意激励的响应 式（2-14）在形式上与单自由度系统的振动方程一致，可以将上式看作单自由系统的受迫振动方程。求解方程（2-14）即计算单自由度系统对任意激励的响应。单自由度系统对任意激励的响应，可以看作系统对一连串单位脉冲响应的叠加。单自由度系统对单位脉冲的响应可以用下式表述： 式（2-15）表示的是单自由度系统对发生在t=0时刻的单位脉冲的响应，单自由度系统对发生在任意时刻的单位脉冲的响应可以写成如下形式： 则方程（2-14）所示的单自由度系统对任意激励的响应可用杜哈美积分表示为： 式（2-15）、（2-16）、（2-17）中，为式（2-14）表示的单自由度系统的阻尼比；为考虑阻尼时系统的第i阶固有频率。式（2-15），（2-16）称为单自由度系统的脉冲响应函数。至此，第i个自由度的响应己经求得。式（2-17）表述的为振动方程的全解，包含了瞬态解和稳态解两部分。2.2.2.4多自由度系统任意激励的响应 上节求得了第i个自由度的响应解。同理可以解得其他自由度的响应解，全部自由度的解可以写成向量，具体形式如下： 式中，由式确定。上式表示的响应为系统在主坐标下的响应。要得到原物理坐标下的响应只要将式（2-18）代入式（2-13）。具体如下： 也可以写成和的形式： 式中，为第i阶正则振型，可由式（2-8）确定，为杜哈美积分，可由式（2-17）确定。2.3 结构动力响应的有限元仿真技术2.3.1有限元分析的一般步骤 有限元方法是求解微分方程的一种非常有效的数值方法，其基本思想是用分片函数（单元形状函数）来逼近原函数，也就把无限自由度问题转化为有限自由度问题求解，通过求解一个线性方程组，得到原问题的近似解。有限元方法的名称是Clough在处理平面弹性问题时，第一次提出使用的。有两种通常与有限元方法相关的方法。一种是力法或称作柔度法，它使用内力作为问题的基本未知量;另一种是位移法或称作刚度法，它使用节点位移为问题的基本未知量。这两种方法在分析中得出不同的未知量（力或位移），并得出与其公式相关的不同的矩阵（柔度矩阵或刚度矩阵）。目前绝大多数的通用有限元程序，如ANSYS， HYPERMESH，ABAQUS等，采用的都是位移法(或刚度法)。有限元方法涉及用相互连接的、称作有限单元的小单元来模拟连续结构。使用有限元工具进行分析计算通常有三个步骤：第一，把对象进行合理的简化建立其有限元模型；第二，按照实际的工况对模型施加约束和激励并进行必要的求解设置，然后开始求解；第三，待求解结束后，综合运用理论、经验或实验结果对计算结果进行分析、论证。2.3.1.1有限元分析模型的建立 在有限元分析过程中，从结构对象的实际模型到计算模型的简化和物理特征提取的过程称为结构对象的特征建模。特征建模的过程也就是对结构实际模型进行简化的过程。实际工程中的结构对象可以分为两种：一种是离散体结构，另一种是连续体结构。杆梁结构体系是最常见的离散体结构，由于其本身存在有自然的连接关系即自然节点，一般可以直接基于这些节点进行单元划分和离散。而连续体离散化过程需要人为地在连续体内部和边界上划分节点，以分片（单元）连续的形式来逼近原来复杂的几何形状。一般情况下，连续体的有限元分析特征模型和结构对象的实际模型不是直接对应的，需要更好地进行问题的特征建模。它需要分析人员具有相应的数学和力学基础、工程分析经验，以及对软件的熟练操作，以便做出准确合理的简化，得到既能反映实际结构对象的特征，又具有合理离散方案的有限元计算模型。在机柜设备结构有限元分析中，一个重要的问题就是如何建立逼真的有限元模型。近40年来实体造型理论和技术的成熟以及大批实用实体造型软件的出现，现在可以很轻松的建立分析对象的几何模型。这样，有限元模型的逼真程度主要取决于物理特征建模，即选定单元类型、设定材料特性、确定边界条件、约束和载荷。2.3.1.2有限元模型的求解 有限元动力响应分析的求解阶段，程序根据前处理阶段建立的网络模型、模型的材料数据、约束和激励的信息，通过单元位移插值函数和材料的本构关系计算模型中各个单元的质量矩阵、阻尼矩阵和刚度矩阵，然后将各个单元的单元矩阵组装成模型的整体矩阵，再求解模型的固有频率及振型、位移响应、应力响应等结果数据。在将建立的有限元模型提交给计算机计算之前，需要指定一些分析计算的控制信息，如采用何种数值算法、需要提取那些分析结果以及模型的动力参数等等。进行合理的求解设置有助于减少计算需要的时间、提高计算效率，在一定程度上还可能影响得到的分析结果。2.3.1.3有限元分析结果的处理由于有限元方法的近似性，分析结果不可避免的存在误差。有限元模型的网格密度会影响有限元分析结果的精度。在一定范围内，适当加大网格的密度，有限元分析的精度会有所提高。但是，网格精度对有限元分析结果精度的影响不是绝对的，当结果的精度提高到一定程度后，再加大网格的密度，计算结果的精度不仅不会提高反而会降低。正由于有限元法是一种工程近似求解的方法，存在误差，所以在用有限元软件进行工程分析计算时，必须对结果进行分析评估。通常以下方式来判定有限元分析结果正确与否：（1） 将分析结果与理论计算结果进行比较。（2） 将分析结果与实验结果进行比较。例如，模态实验可以作为判别有限元模型正确与否的标准。（3） 根据己有的经验判断分析结果。 了解有限元分析结果的特点和产生误差的原因对理解和解释有限元分析结果具有重要的意义。只有经过检验为正确的分析结果才是真正对工程应用有意义的。2.3.2 ANSYS用于机柜设备结构动力响应分析的可行性ANSYS软件是结构、热、流体、电磁、声学于一体的大型CAE通用有限元分析软件，可以广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、电子、造船、轻工等一般工业及科学研究。ANSYS软件具有丰富和完善的单元库、材料模型库和求解器，从而保证了其出色的分析计算能力。ANSYS的结构分析计算功能非常强大。它可以完成结构静力分析、结构动力学分析和结构非线性分。其中，ANSYS可以进行的结构动力学分析类型包括:模态分析、瞬态动力学分析、谐波响应分析及随机振动响应分析。本论文涉及到ANSYS软件的动力学分析功