机械振动数值分析.pptVIP

雅可比变换法定义阀值收敛条件第94页，共148页，星期日，2025年，2月5日雅可比变换法收敛性检验考虑特征值问题，其中由式（47）可计算得第95页，共148页，星期日，2025年，2月5日雅可比变换法经第一次循环，对非对角元（1，2）消零，得由式（52）可计算得第96页，共148页，星期日，2025年，2月5日雅可比变换法类似的，分别对非对角元（1，3）和（2，3）消零，得第97页，共148页，星期日，2025年，2月5日雅可比变换法雅可比变换法的一般步骤给定当前循环的阀值，一般取第m次迭代的阀值为10-2m；对上三角块的非对角元按式（48）计算，若计算结果大于给定的阀值，则按式（47）进行旋转变换；检查式（49）式是否满足，若满足则由迭代最终步骤计算出的矩阵给出特征值，计算特征向量，否则转向第1步，开始新一轮循环。第98页，共148页，星期日，2025年，2月5日广义雅可比变换法对雅可比变换矩阵进行修改，得转换矩阵第99页，共148页，星期日，2025年，2月5日广义雅可比变换法系数的选取应使刚度矩阵和质量矩阵在变换后其第i行j列的元素同时为零，即若，可取，否则第100页，共148页，星期日，2025年，2月5日广义雅可比变换法其中第101页，共148页，星期日，2025年，2月5日广义雅可比变换法例1，用雅可比法求解图示系统的固有频率第102页，共148页，星期日，2025年，2月5日广义雅可比变换法第103页，共148页，星期日，2025年，2月5日广义雅可比变换法精确解第104页，共148页，星期日，2025年，2月5日单元内部积分二维和三维问题第62页，共148页，星期日，2025年，2月5日约束条件的处理直接代入法特点方程阶数降低第63页，共148页，星期日，2025年，2月5日约束条件的处理对角元素置1法特点适用于零位移情况第64页，共148页，星期日，2025年，2月5日约束条件的处理对角元素置大数法特点相当于罚函数法第65页，共148页，星期日，2025年，2月5日等参变换等参变换的基本思想就是局部系向总体系的映射局部坐标系向物理坐标系的映射类似于构造形函数，可构造物理坐标和局部坐标之间的关系。仍然以三角形平面单元为例第66页，共148页，星期日，2025年，2月5日等参变换需要注意的地方等参变换第67页，共148页，星期日，2025年，2月5日等参变换对面积分来说三维问题？？？等参变换第68页，共148页，星期日，2025年，2月5日等参变换后的单元刚度矩阵、单元质量矩阵格式等参变换单元等效结点载荷的一般形式第69页，共148页，星期日，2025年，2月5日第三章大型系统特征值问题系统运动方程对特征值问题，可建立系统运动方程如下自由振动解将式（2）代入式（1），得第70页，共148页，星期日，2025年，2月5日系统运动方程由式（3）得特征方程解得n个特征根及相应的特征向量n次代数方程另一种形式第71页，共148页，星期日，2025年，2月5日特征向量正交性取结合式（6），得第72页，共148页，星期日，2025年，2月5日特征向量正交性进一步得正交性关系一个重要的关系第73页，共148页，星期日，2025年，2月5日刚度矩阵为半正定矩阵质量矩阵为半正定矩阵半正定矩阵的特征值问题特征向量特征值2自由度系统的例子第74页，共148页，星期日，2025年，2月5日Rayleigh商对给定的振型向量，可定义Rayleigh商Rayleigh商的上下界考虑式（12），可将式（15）表示为第75页，共148页，星期日，2025年，2月5日Rayleigh商将式（16）展开，得物理意义——假设模态相当于增加了约束，使得系统刚度提高第76页，共148页，星期日，2025年，2月5日误差估计设和是标准特征值问题的近似解，定义残差考虑下式得第77页，共148页，星期日，2025年，2月5日误差估