有限元讲义6-动力分析有限元.doc

第6章 结构动力分析有限元法 此前述及的问题属于静力分析问题，即作用在结构上的荷载是与时间无关的静力。由此求得的位移、应力等均与时间无关。实际工程中的大部分都可简化成静力问题。但当动载与静载相比不容忽略时，一般应进行动力分析。如地震作用下的房屋建筑，风荷载作用下的高层建筑等，都应计算动荷载作用下的动力反应。研究课题中以动力问题为主。 解决动力问题有两大工作要做：一是动荷载的模拟和计算，二是结构反应分析。本章将讨论如何用有限元来解决动力计算问题。 6.1 结构动力方程 一．单元的位移、速度和加速度函数 设单元的位移函数为； 6—1—1 式中： 单元位移函数列阵、结点位移函数列阵均是时间t的函数。 由6-1-1可求得单元的速度、加速度函数： 6—1—2 6—1—3 二．单元的受力分析 设图示三角形单元，当它处于运动状态时，其上的荷载一般应包括： 单元上的荷载；单元对结点的作用力, 单元内部单位体积的： 惯性力： 6—1—4 阻尼力（设正比于运动速度）： 6—1—5 干扰力（已知的条件）： 根据达朗贝尔原理，上述四力将构成一瞬时平衡力系，使单元处于动平衡状态。 为此寻求四者之间的关系； 三．结点力与结点位移、速度和加速度之间的关系 用虚功原理推导： 令单元结点发生任意可能的虚位移，它满足单元所定义的位移场，即虚位移场成立。 作用在单元上的外力所作的外力虚功： 单元内部应力在由于虚位移所引起的虚应变上所做的内力虚功： 根据虚功原理（T=W），若将惯性力，阻尼力用上面的6—1—4，6—1—5代替，得： 由于虚位移的任意性，可从等式两边各项中消去 ，得： 简写为： 6—1—6 式中： 单刚（第一项为弹性恢复力） 单元阻尼矩阵（第二项为阻尼力） 质量矩阵（第三项为惯性力） 包括由作用在单元上的干扰力转化成的等效结点荷载 6—1—6即为单元结点力之间的关系式。 四、结构的动力方程 有了上述单元力关系式，象在静力问题中对每个结点建立平衡方程一样，根据达朗贝尔原理，对每个结点建立动平衡方程后，即可得到结构的动力方程组： 6—1—7 式中：分别为总质量、总阻尼、总刚度矩阵。 为外力（结点力），实为干扰力（当不考虑静载时） 当不计阻尼影响时，上式成为： 6—1—8 若干扰力为零，得： 6—1—9 即结构的无阻尼自由振动微分方程组。由此可求得结构的自振特性（频率，振型）。 由上可见，动力问题首先要解决及的形成和方程组的求解问题，下面逐一加以讨论。 6.2 单元质量矩阵和单元阻尼矩阵 由于单元质量矩阵表达式中包含了推导单元刚度矩阵时相同的形函数，因此常将按此式形成的称为协调质量矩阵（或一致质量矩阵），下面对此讨论。 一、几种常见的协调质量矩阵 梁单元（如图） 设梁单元位移函数： 式中形函数 设单元的质量沿梁的长度方向均匀分布，则有： 6—2—1 二、集中质量矩阵 从上可知，单元的协调质量矩阵和单刚具有相同的阶数，因此，从总质量矩阵的阶数也与总刚相同。或者说采用协调质量矩阵后，结构的振动自由度和结构的 静力自由度是相同的，动力问题的这种做法，其求解是很费时的： 形成质量矩阵的工作量等同于总刚 特征值的求解 但是工程实际和试验证明，在某种干扰力作用下，结构的动力反应是有明显主次之分的。因为工程上通常把单元的分布质量集中到各结点而成为集中推聚质量，这样可使问题得到很大简化，且计算经验表明，二者给出的计算精度相差无几（如果小而重的物体放置在一个轻型结构的结点上集中质量公式是近乎精确的） 质量集中按静力等效原则，且常忽略转动惯量的影响，上述各单元的质量矩阵简化为： 梁单元： 常应变三角单元： 矩形薄板单元： 由此可见，采用集中质时，集中质量即采用“就近堆积”的原则。 三、阻尼矩阵 在结构动力分析中，较多采用的是粘滞阻尼理论，即假定阻尼力与速度成正比，由此得到单元的阻尼矩阵： ① 这样做，可给方程组的求解带来方便。但这个假定并不能很好的符合结构的实际情况。因此在实际应用中也常采用M和K的线性组合： ② 或 ③ 也可写成更一般的形式： ④ 当取S=0、1两项时，即为式② 式中系数由下式假定： 由于阻尼矩阵依赖于质量矩阵和刚度矩阵，故可通过、而获得。 6.3 结构的自由振动和特征值问题 一、特征值方程 这是一个大家都很熟悉的问题，故着重讨论程序设计上的一些处理方法。 结构作无阻尼自由振动的微分方程： 6—1—9 设结构作简谐振动： 将其代回6—1—9得： 求解特征方程 6—3—1，即可获得几个运动自由度所对应的频率和振型。 还有另