动力学问题的有限元法.pdf

资源描述

动力学问题的有限元法.pdf

《动力学问题的有限元法.pdf》由会员分享，可在线阅读，更多相关《动力学问题的有限元法.pdf（22页珍藏版）》请在冰豆网上搜索。

动力学问题的有限元法.pdf

第七章动力学问题的有限元法第七章动力学问题的有限元法结构动力学是研究动载荷作用下结构动力反应规律的学科，讨论结构在动力荷载作用下反应的分析方法，寻找结构固有动力特性、动力荷载和结构反应三者间的相互关系。

研究结构在动力荷载作用下的反应规律，能够为结构的动力可靠性（安全、舒适）设计提供依据。

前面介绍的静力学问题的研究对象是受不随时间变化的载荷作用。

而动力学问题的对象受随时间而变的载荷的作用，从而使在结构中产生的位移、速度、应力和应变都随时间而变。

当结构受随时间变化的载荷作用，且这种载荷的作用对结构的变形和应力的产生起主要作用，以致影响设备的安全性，或舒适性。

这时就要进行动力学分析，充分认识其规律性，从设计阶段就抑制这种不利状况的发生。

例如，有时虽然动载荷不大，但结构在交变力的作用下，其某些固有频率与激励力的作用频率相接近时，就会引起很大的振动、变形或应力，这时，就必须对结构作动力学分析。

又如，要利用结构在周期性作用力驱动下的定向振动，例如利用这种运动输送产品，这时，就必须巧妙地设计结构，使其具有某些与激励频率一致的固有频率，并且使结构对激励具有适当的响应能力。

总之，不管是利用振动，还是抑制振动，都需要进行结构动力学分析。

当前结构动力学的研究内容有三类。

第一类问题：

反应分析（结构动力计算），第二类问题：

参数（或称系统）识别，第三类问题：

荷载识别。

第一类问题是已知系统动态特性和动载荷作用部位及大小，求出系统的响应随时间变化的位移，速度，加速度和应力等。

第二类问题是已知系统的输入输出特性，分析系统固有的动态特性，结构模态分析就属于这一类问题。

第三类问题是在已知系统动态特性的条件下,通过测量系统的响应,或由响应准则预先给出响应要求,以此识别对响应的外载荷。

三类结构动力学研究内容的载荷、结构和响应之间的关系如图7-1所示。

动载荷种类大致分类如图7-2所示。

图7-1结构动力学研究的内容图7-2动载荷种类1本章主要介绍结构动力学分析的基础知识，并主要介绍系统固有特性的有限元分析方法有限元模态分析。

主要知识点：

1.导出有限元的结构动力学方程和模态分析有限元方程。

2.介绍模态分析方程中质量矩阵与阻尼矩阵的计算。

3.把计算结构的固有频率与振型的问题归纳为一个求特征值的问题。

7.1结构动力学方程及有限元方程7.1结构动力学方程及有限元方程系统的受迫振动的微分方程：

m（）xcxkxRt+=?

其中：

质量，阻尼系数，弹性系数。

mck,xxx?

分别为加速度，速度和位移。

1.单元的动力方程：

（用“虚功原理”求解）结构动力学有限元分析与静力学的有限元分析一样，首先要建立单元的有限元方程。

为求解起见，要引入达朗贝尔原理。

质点的达朗贝尔原理可表述为：

当非自由质点运动时，主动力，约束反力和惯性力构成一个动平衡力系。

惯性力是一种虚构的力，由质点的加速度引起。

mFNS即：

0fns+=现考察单元的动力学方程的建立。

1）.准备工作：

将位移函数表达成近似函数。

单元内任一点位移（,）（）eNNxyz=et，是单元内位移的近似函数,这儿N只是位置的函数，与时间无关。

而e与时间相关。

（,）（）eeNNxyzt=eeNN=?

2）.分析单元上的作用力。

图7-3单元分析主动力：

集中力，面力，体积力：

psRRRRV=+阻尼力：

与速度相关，并与速度方向相反，两者的比例关系为阻尼系数：

ceRCCN=?

2惯性力：

与加速度，质量相关：

单位体积惯性力meRN=?

阻尼力、惯性力与主动力方向相反。

密度，单位体积质量。

3）.建立动力学平衡方程：

利用虚功原理和达朗贝尔原理，系统处于动平衡，首先将载荷移置到节点上去。

利用虚功原理。

外力做的虚功：

11*11*TeTTppTeTTSsSsTeTTVVVVTeTTVCVCTeTTVmVmRNRRdSNRdSWRdVNRdVNRdVRdVNRdV=集中力：

面力：

体积力：

阻尼力：

惯性力：

RdV内力做的虚功：

*TeTTeVVUdVBDB=dV0虚功原理：

（）WUWU=，最小势能原理根据虚功原理，应有WU=1*1*eTeTeTTTpSsVVeTeTeTTTTeVCVmVNRNRdSNRdVNRdVNRdVBDBdV+=由*eT的任意性，可约去，且由阻尼力，惯性力的定义可知：

cmeeRCCNRN=?

，代入上式，有11TTTpSsVVeeTTTeVVVNRNRdSNRdVBDBdVCNNdVNNdV+=+?

=eRi令：

11TTTepSsVVTeVTeVTeVNRNRdSNRdVRKBDBdVcCNNdVmNNdV+=上式可表达为：

，这就是单元的动力学有限元方程。

eeemcK?

+=+=2.结构整体动力学有限元方程与静力学的集成原理相同：

根据局部单元节点编号与整体的节点编号的对应关系，利用，即一个节点上的节点力是该节点所在的所有单元的相应节点（同这个节点）上的节点力的总和来集成。

eeiiRR=3位移唯一性。

得MCK+=?

ReMm=，eCc=eKk=其中，当0C=，称为无阻尼受迫振动，此时MKR+=?

当0,0CR=，称为无阻尼自由振动，此时0MK+=?

7.2单元质量矩阵和阻尼矩阵的表达式7.2单元质量矩阵和阻尼矩阵的表达式在动力学有限元方程中除刚度矩阵之外，还有两个重要的系数矩阵：

即质量矩阵和阻尼矩阵。

这一节中将介绍单元质量矩阵和单元阻尼矩阵的表达式。

一.单元质量矩阵单元质量矩阵有二种表达形式：

一致质量矩阵和集中（堆积）质量矩阵。

1.一致质量矩阵：

用公式计算的单元质量矩阵，称为一致质量矩阵，是因其中的TeVmNN=dVN就是位移模式形函数N而得名。

或者说，推导质量矩阵采用的形函数与采用推导单元刚度矩阵时采用的形函数是一致的。

以平面应力问题为例：

三角单元，ijmNNININI=m（位移模式中的000000ijmijNNNNNNN=N）图7.3三节点三角形单元则三角形单元的一致质量矩阵em：

（以单位厚度计）ii22jj2j0N0N0N00N00N00jmijmijmeTijNNNNNNNNmNNdxdyNN=?

j2mm2mmNN0N000N0N0mijmijmdxdyNNNNNNN?

由于m，利用数值积分：

2iNii00N0N,iijjmNLNLNL=4!

（2）abcijmabcLLLdxdyabc=!

，有22!

24!

62（112）!

12iijNdxdyNNdxdy=+从而得到20210212010210102010102em=对称2.集中质量矩阵：

将单元质量假想地集中到节点上去，也就是说只在节点上有质量。

这样，一个节点的加速度就不会影响其它节点的惯性力（相互之间不受影响）具体处理方法是将质量平均分到所有节点上去。

例如三角形三个节点单元，每个节点分到1/3的质量，即1m3it=,t为厚度，为三角形面积。

111111111331111ewmt=资料表明，集中质量矩阵计算出来的频率与一致质量矩阵的计算结果相差无几。

在单元相同的情况下，集中质量矩阵计算的结果低于一致质量矩阵的计算频率，而振型则是一致质量矩阵的计算结果精确。

由于集中质量矩阵式对角线矩阵，所以计算中常采用。

3.附加质量当有附加的质量在节点上时，要在计算整体质量矩阵时，加上这些附加质量cm即ecMmm=+二.单元阻尼矩阵根据阻尼的成因可以分为粘滞阻尼和结构阻尼。

1粘滞阻尼（比例阻尼）前面得到的阻尼矩阵是将阻尼看作是正比于节点的运动速度，称为粘滞阻尼。

它消耗振动的动能（能量）。

从可知，这个阻尼矩阵与质量矩阵成正比。

TeCCNNd=VdVeTmNN=eeCCm=，/Cc=如果,c均为常数的话。

所以称为比例阻尼。

（相对于质量矩阵而言）NTeCCNd=V2结构阻尼比例于结构的变形速度?

，由于材料内部的摩擦引起的，称为结构阻尼。

这时的阻尼力可以简5化为cRD=?

，这样可得到单元的结构阻尼矩阵使其具有的形式，所以结构阻尼正比于刚度矩阵：

BeTkCCBDd=VeekCk=3Rayleigh阻尼在实际分析中，精确决定阻尼矩阵是相当困难的，通常允许将整体阻尼矩阵简化为质量矩阵M和刚度矩阵K的线性组合，这样可用来解耦有阻尼的振动方程。

M+KC=，称为Rayleigh阻尼。

其中j22i222（-）2（-）ijiijjijjiji=其中22iii+=称为结构的第i阶振型的振型阻尼比振型阻尼比。

2i+称为结构的第i个校准振型的振型阻尼常数。

工程中一般取i=0.020.24，它与结构的类型，材料的能量消耗特性及振型有关。

可通过实测，从系统的固有频率和辐频特性曲线获得。

质疑质疑：