第八章 弹性连续体振动的有限元解.docVIP

第八章 弹性连续体振动的有限元解 上一章讨论的弹性连续体振动问题都只是在简单的特殊边界情况下才能得到精确解，而对于复杂弹性连续体的振动，通常无法得到精确解。因此，只能采用近似解，近似解方法很多，其要旨在于将无限自由度系统（连续体）变换成为有限多自由度系统（离散系统）来处理。下面仅介绍现代工程中常用的有限元分析法（简称有限元）。 有限元的基本思想是将一个复杂结构（连续系统）看成是有限个基本元素（单元）在有限个结点彼此相联结的组合结构。每个单元都是一个弹性体。有限元法通常是采用位移法，即以结点处的位移作为基本未知量，单元的位移是用结点位移的插值函数表示，单元以至整个结构的一切参数包括位移、应变、应力等都通过结点位移表示出来。从振动问题来看，最后是将一个连续体的振动问题变成了一个以有限个结点位移为广义坐标的多自由度系统的振动问题。由于平面刚架所采用的单元综合有受轴向力杆及横向弯曲的两种作用，因此本节着重通过平面刚架（杆系结构）的振动为例来阐明有限元的基本概念与分析步骤。 有限单元法分析过程基本上可分为结构离散化、单元分析、整体分析三个步骤。 一. 结构的离散化 如图8-1所示平面刚架，它由三根杆件组成，杆件之间通过刚结点2、3相互连接，若在某些结点上受有作用在刚架平面内的激扰力，这就是一个平面刚架的振动问题。结构的离散化就是将结构分成若干个单元。在此例中，取三个杆各为一个单元。共分成三个单元，在图中分别以①、②、③标明各单元的标号。各单元之间以及它们与基础的联结点称为结点。结点编号分别为1、2、3、4。结点1、4是杆件和基础固定联结的结点，它们不能运动，称为固定结点。2、3点是杆件间的联结点，是可以运动的，称为可动结点。可动结点的数目乘以每个结点的基本未知量的数目就是整个结构的自由度数。对于工程结构来说，结构划分的单元愈多，则结构的自由度数愈多，计算工作量就愈大，但计算结果愈有可能趋于精确解，因此在结构离散化过程中，划分单元的多少，应综合考虑精确度与计算工作量的因素。 二. 单元分析 单元分析的目的在于建立结点力与结点位移的函数关系，此步是有限元分析中最关键的一步。在单元分析中，首先将作用在单元上跨间载荷按照静力等效原理移到该单元的结点上去而成为结点载荷，单元的运动可看作是在结点力作用下的运动。在平面刚架中所采用单元的动力性质将是纵向振动杆和横向振动梁这两种性质的综合。在线性弹性范围内这两种性质是可以简单迭加的，即它们之间是不耦合的。现取一代表性单元如图8-2所示，在单元分析中，取坐标原点在单元左端，x轴沿单元轴线，这种和单元连结在一起的坐标系称为局部坐标系，相应地表示整体位置的坐标系称为总体坐标系，如此建立局部坐标系的目的是使单元特性分析较简单。设单元的长度为L，截面面积为A，惯性矩为J，单位体积质量为，材料的弹性模量为E及泊松比为，单元的左端、右端编号分别用i和j表示。 当杆同时作横向振动和纵向振动时，每个结点的广义坐标是轴向位移U（t）、横向位移V（t）和转角，相应的广义力是轴向力X（t）、横向力Y（t）和弯矩M（t），每个结点的广义位移和广义力分别如图8-2所示，规定位移和力按图中所示的方向取为正。以q表示位移，F表示力，则结点的位移和力的矩阵表示形式如下， ， （8-1） 其中 ，， ， 以下先分别研究只发生纵向振动的杆单元和只发生横向振动的杆单元的特性，然后再把它们迭加，得到这两种振动同时发生时的单元特性。 1. 纵向振动杆单元的特性分析 建立单元位移和结点位移的关系 （1）纵向振动杆单元 杆纵向振动时的基本未知量是杆各个截面的轴向位移u（x，t），而有限元分析中杆单元纵向振动的基本未知量是杆的结点位移ui、uj，杆任一x截面的轴向位移u（x，t）则是利用结点位移ui、uj经插值得到，即 （8-2） 式中 Niu（x）、Nju（x）为位移插值函数（也叫作形函数）。 插值函数实质上是一种假设模态，由于每个单元的尺寸很小，故插值函数的形式一般取得很简单，常取为多项式形式，且总是按常数项、一次项……的次序来增加项次（取常数项可描述单元刚体运动，取一次项可描述单元的常应变），对现研究的均匀杆在常值结点力作用下的静位移情况，由于杆单元的轴向位移状态是由两个位移参数所决定，因此杆单元的位移模式可设为包含两个待定常数的多项式 (a) 待定常数可根据结点轴向位移条件来确定。结点轴向位移条件是：x＝0时，u＝ui；x=L时，u＝uj。据此由（8-2）式及（a）式可求得，及 Niu（x）＝1－x/L, Nj