机械系统动力学第7章

1、第七章 基于有限元法的振动分析,71 弹性力学基础： 方向变形场： 非线性的格林应变张量 线性的柯西应变张量 应力张量,弹性模量： 各向同性时： 式中，克罗尼科尔符号：,式中： E为弹性模量，v为泊松比 二维时： 一维时：,不同应力状态下弹性模量 （对称阵）的表达： 变形能 1三维应力状态 材料各向同性时： 式中： 2平面应变状态 3平面应力状态 4单轴应力状态,广义虎克定理: 则应力可表为： 或:,举例： 拉压杆的变形与应变、应力： 由边界条件 设变形场为: 解得： 变形场： （线性）柯西应变： 说明杆的应变为常数。,（非线性）格林应变： 是精确的，与端位移成非线性。 问题关建是确定变形场

2、再举拉压杆为例： 已知两端位移为 ，设 方向对应位移场为 ，故设有：,由边界条件 可得 方向位移场 则可得对应柯西张量的线性应变 和格林张量的非线性应变 ：,72 有限元基础 插值理论 （真实）精确函数 ， 插值函数 希望通过插值函数 得到精确函数 721 拉格朗日插值： 在特定点P处两函数值相等， 插值函数的构成 是形函数， 是结点函数值，如结点位移。,形函数构成 范德蒙阵， 基向量 要注意基向量的完备性，不持没有根据的偏向性。 将 代入 ，注意到 为 插值函数 在p点处的值,同时，将任一点q的位置坐标代入插值函数 的基向量 ，则可得点q处函数值 ： 因为 ，对照上式有： 即 矩阵形式表为

3、由 得形函数：,拉压杆单元为例： 杆具有两结点位移，故取基向量 的项数为2： 杆始端处 杆末端处 将 点的位置坐标代入基向量 得：,求逆后得形函数所需的范德蒙阵 ： 得形函数为： 变形场为： 应变： 应力：,6结点平面三角单元 基向量 ， 点坐标,将点坐标代入基向量 得 ，求逆后得形函数所需的范德蒙阵 ： 得形函数：,722 高阶插值 （ Hermite插值） 真实函数 和插值函数 在特征点p处不再单纯是函数值 相等，而是要满足经某种运算后仍相等： 某特定算子 插值函数的构成： 形函数构成：,平面弯曲梁单元的变形场： 梁单元边界位移 和基向量 为： 对应四个算子为：,将以上算子用于基向量 得：

4、 求逆后得形函数： 位移场为：,平面梁位移场2： 对应算子为： 将算子用于基向量 得：,求逆后得形函数： 则位移场可得：,73 有限单元方程 知位移场： 则得应变： 式中： 为线性应变， 为非线性应变。,在一般线性有限元分析中，我们只取线性应变部分，即： 在以下的分析中，我们不再使用非线性的格林应变张量，而用线性的柯应变张量： 矩阵形式的位移场、应变、应力表为：,由虚位移原理建立单元运动方程： 虚功原理：内力虚功等于外力虚功 内力虚功 单元静力平衡方程： 为简化计算，先忽略分布力，由式 得：,矩阵形式的单元静力平衡方程表为： 式中：,单元动力平衡方程： 以体积分布力的单位质量惯性力为： 单元中

5、的惯性力虚功为： 式中：,同理，假定阻尼正比于点的速度，单元中的单位阻尼力则为： 单元中的阻尼力虚功：,将式 和式 代入 得单元动力平衡方程： 或表为矩阵形式： 式中： 分别称为刚度阵、阻尼阵、质量阵，它们皆与形函数有关。,运用拉氏方程建立单元动力平衡方程： 其中拉氏函数： 动能： 取 则 位能： 耗散函数：,已知： ， ， ， 可得： 由拉氏方程得单元运动方程：,例：杆单元刚度阵K 推导： 由前已得到位移场 ：,杆单元质量阵M： 同理，可推出杆单元阻尼阵C：,梁单元运动方程：,先考虑 方向位移场 ： ， ， ， 基向量 将算子用于基向量 得： 得 方向形函数 和位移场 ：,注意到 位移场 和 位移场 的关系： 上式 对 积分得： 梁形心轴在 方向之位移场：,由梁两端点条件： 得基向量 等： 从而确定梁形心轴在 方向之位移场 至此， 位移场和 位移场皆可得到，对应的形函数为：,因为： 该梁为单轴应力状态，即又： 将形函数 和 阵等代如式 得梁单元刚度阵K、阻尼阵C、质量阵M