有限元第七讲.pptVIP

第七章 动力问题的有限单元法 1、一致质量矩阵 三、质量矩阵 3）平面常应变单元 第七章 动力问题的有限单元法 2、集中质量矩阵 三、质量矩阵 将分布质量按某种原则换算成结点集中质量，按单元动力自由度的顺序放入相应位置形成的单元质量矩阵，称作单元集中质量矩阵。 当质量均匀分布时，常按照结点所分担的线段、面积和体积确定该结点集中质量的大小。 因为假设集中质量集中成质点，故没有转动惯量，与转动自由度相对应的质量为零。 1）平面杆单元 第七章 动力问题的有限单元法 2、集中质量矩阵 三、质量矩阵 2）平面常应变单元 第七章 动力问题的有限单元法 三、质量矩阵 1）单元质量矩阵和总体质量矩阵均为对称矩阵。如果是一致质量矩阵，则是正定的，如果是集中质量矩阵，则至少是半正定的。 3、质量矩阵的性质及特点 2）单元一致质量矩阵为满阵，数值计算费时，总体一致质量矩阵具有总体刚度矩阵同样的带宽。集中质量矩阵为对角阵，占用内存较少，计算简单，节约机时。 3）用一致质量矩阵算得的频率是结构真实频率的上限，用集中质量矩阵算得的频率是结构真实频率的下限，但两者相差不大。 工程上采用集中质量矩阵计算的情况居多 第七章 动力问题的有限单元法 三、质量矩阵 3、质量矩阵的性质及特点 例：分别采用一致质量矩阵和集中质量矩阵有限元法分析等截面悬臂梁的自振频率。 将长度为L的悬臂梁用N个长度相等的单元离散，为分析单元数量的影响，分别取N=1,2,3,4,5 对于一致质量矩阵，每个结点有两个（动力）自由度，而对于集中质量矩阵一个结点只有一个动力自由度。 第七章 动力问题的有限单元法 三、质量矩阵 3、质量矩阵的性质及特点 从例题可以看出： 1）一致质量矩阵法的上限收敛性和集中质量矩阵法的下限收敛性。 2）增加单元的数量可以很快提高分析精度。 3）有限元法分析低阶自振频率的精度高于高阶频率。 4）对此问题，虽然同等条件下一致质量矩阵法的精度高于集中质量矩阵，但其计算量是后者的两倍。 第七章 动力问题的有限单元法 四、阻尼矩阵 1、单元阻尼矩阵 在运动方程的推导过程中已得到单元的阻尼矩阵为： 式中右侧第一项是假定阻尼力正比于质点运动速度的结果，为粘性阻尼，如果单元阻尼系数μ1,为常数，则此项比例于单元质量矩阵： 第七章 动力问题的有限单元法 四、阻尼矩阵 1、单元阻尼矩阵 在运动方程的推导过程中已得到单元的阻尼矩阵为： 式中右侧第二项是假定阻尼力正比于应变速度的结果，为材料阻尼，如果单元阻尼系数μ2为常数，则此项比例于单元刚度矩阵： 所以： 第七章 动力问题的有限单元法 四、阻尼矩阵 2、总体阻尼矩阵 单元阻尼矩阵： 由于μ1、μ2难以确定，所以常采用下面的公式确定总体阻尼矩阵： —瑞利阻尼 式中： 第七章 动力问题的有限单元法 1、特征值问题 无阻尼自由振动方程为： 五、结构的自振特性分析 是一个常系数齐次线性常微分方程组，其解的形式为： 带入自由振动方程得： 上式是齐次线性方程组，有非零解的条件是： 第七章 动力问题的有限单元法 1、特征值问题 五、结构的自振特性分析 如果[K]和[M]的阶数是n，则 是 的n次方程， 称其为自由振动特征方程，通过它可解出n个特征值，将这些特征值 再带入 可解出n个特征向量 第 i 个 、 合称第i个特征对， 为结构的第i个固有频率, 为结构的第i个振型 第七章 动力问题的有限单元法 1、特征值问题 五、结构的自振特性分析 将 按从小到大的顺序排列： 其中 称作基本频率，相应的振型 称作基本振型。 对于上式可以采用广义雅克比法，拟迭代法、子空间迭代法等数值方法直接求出特征值和相应的特征向量。 式 在数学上称为广义特征值问题，常记作： 第七章 动力问题的有限单元法 2、特征值和特征向量的性质 五、结构的自振特性分析 1）当[K]和[M]是实系数对称矩阵时，其特征值一定是实数，且特征向量也是实向量。 2）如果[K]为正定矩阵，则特征值一定是正实数，如果[K]为半正定，则特征值为非负实数，且特征值为零的个数等于结构刚体位移自由度的个数。如果集中质量矩阵为半正定，其对角线上有r个零元素，则n个特征值的最后r个为无穷大。 3）振型的规格化（三种方法）： a）以第一个元素为1规格化； b）以振型中的最大元素为1规格化； 第七章 动力问题的有限单元法 2、特征值和特征向量的性质 五、结构的自振特性分析 3）振型的规格化（三种方法）： a）以第一个元素为1规格化； b）以振型中的最大元素为1规格化； c）以矩阵[M]、[K]进行规格化，使振型满足： 这时振型向量的各个元素应除以 第七章 动力问题的有限单元法 2、特征