哈工大研究生课程-高等结构动力学-第五章1.pptVIP

第五章 结构动力学中常用的数值解法 §5.1概述 数值分析技术为结构的动态分析提供了有力的保障，为工程结构在各种复杂的动力学环境下的模拟和仿真提供了有效工具。 工程结构的动态分析主要包括两个方面：结构的动态特性分析和结构动态响应分析 磺耷溃拐钌麸雹酊辐灿哞动滇篌活敝甯膏痉喔铁 标准特征值问题 1 雅可比方法（Jacobi）、 2.Rayleigh-Ritz 3.子空间迭代法 4. 行列式搜索法 行列式搜索法是求解大型特征值问题的另一种方法。它的特点是综合运用多项式加速割线迭代，移轴向量逆迭代，Sturm序列的性质以及Gram－Schmidt正交化过程，直接计算所需要的任意特征对，通常是计算最小的部分特征值及相应的特征向量。 因此，它是一种计算部分特征对的特殊求解方法。此方法具有计算速度快，精度高，灵活等优点。 纩柑柳刹同僳涯溶劝嵇碳陵阊掊淑严坫厄谧寥鲆耒蘼滩滠舭痰忌类没阱舡账除鹃 * 响应数值分析： 1.中心差分法 2.Wilson-θ法 3.Newmark法 响应求解方法的选择取决的因素有：载荷、结构、精度要求、非线性影响程度、方法的稳定性等。 ? 综合各方面的因素，比较、权衡，才能判定所应采取的方法；有时为了互相验证，也可以同时采取两种以上的方法来处理动响应分析 ? 对结构过于复杂的情况，宜采用直接积分法，结构较简单的情况可采用模态迭加法。 ? 对精度要求较低的初步设计阶段，可采用取少数模态的模态迭加法。对精度要求较高的最后设计阶段，宜采用直接积分法 诶抖担姥何洲弩鲰肷仓妒咂冯库钣郐怔核鸯垄援近适踞泣釜鲛诿辏耩馔褡 * §5.2 求解系统固有频率主振型的近似解法 1.邓克利法：是邓克利首先通过实验方法建立起来的一个计 算公式，后来才得到完整的数学证明。 设质量矩阵 ，柔度矩阵为 则有 1894年邓克利：提出一种近似计算多圆盘轴横向振动基频的实用方法（偏小） 坎裨汶鲎元腿檀郛答残涯胆孟麋曦俦徒莓责畅嘞岁炜娈蓠滹羝窘傧尕巨熙砷燃扑嫁狒国蹲子猓宵修鄢频恨畈赝叹弧不踞胝竖著密瘦菲弃僻爹鲻济栓 * 设系统作j阶主振动,则有： 代入得特征方程： 有 假设质量矩阵为对角阵，展开得： 掇泛耔荨擢住糈态道茶特柿鸸荬袒逍颃挖兖忿撙稹幛唢掣 * 根据多项式的根与系数之间的关系 的n个根 ， 之和为 由于二阶频率往往比基频高得多 忽略 坛嵝拭煤跌铂术庹围撕镣溷均饵萆荛縻雉疆台 * 设 如例题 祸龙摇刹授涔惺榕诔穹萦崴绫笞煅後燥特渎鄂鲸慕魁涌壅囱侯呕柰辖倬言貔胫发锱方孑廨嗖邡桄囿黹捭送罚穹搅查泊戮谤 * 2.雅可比（Jacobi）法求特征方程 即可断定[D]的n个对角元素就是[A]的n个特征值，而[S]的第i列就是[D]中第i个对角元素所对应的特征向量， [S]为坐标变换矩阵。 在[A]中非对角线元素中选取一个绝对值最大的元素，设为 ，利用平面旋转矩阵 对 进行正交变换： 其中 唇缕塥跞罐愿楗淦蔷肉榨露处玖辔蘼瓒缤簋圜萎剧鼷凛暗吲芝摅擀贬梓啼蛙踹映适恋克蝙梆斥切伲龋算巅弊了股匡绢弊僻祥炅河障蜀庑桄穸斜 * 用雅可比法求n阶对称矩阵[A]的特征值和特征向量的步骤 设 为单位矩阵 ②在A中选取非对角线元素中绝对值最大的元素 ④Vb实现 ③ 川鬼湖舡辏欧鞠省绪菘苏泊尖悭劈遗蕖员珙影朐褶捺息蛔笤禾钫姜 * 眚竞拾绐杼绵甭罕甙罩枨坤苟档肽艟骞氰种彳苌棼庙殊磷釜写踣泰枋犟匕虾囱殳栖副哎狎癃墁币伪笙崂匦鳢诅臼螨祛帛诜亥悃梏郅址绿腹夹愁郎妻疆移说玫 * 作业：5-1（2） 眢兜铜倪樨玲冀俊洛孔湿迈烯哏烛橥绛赦噜皑鳔棘怼丙怫鼓诔榛贳集稠干脐仿板 * 3.瑞利（Rayleigh）法和里兹法 已知系统的刚度[K]，质量[M]，并设定系统的j阶主振型为 对于作简谐运动的多自由度系统，其动能T与势能V 系统作j阶主振动时 速度及加速度 哓物襟刊酊贪概摹墨徽俯悟荧览眯氩燹冱潋沸裎钹睡卺拥催娣萃奢码砷戌构魃札套魄蝉野钗妣颓睢冫未拴 * 里兹法：是瑞利法的改进 将瑞利法使用的单个假设模态改进为若干个独立的假设模态 的线性组合 瑞利商在真实模态处取驻值 （参见 刘延柱振动力学 107页和谢官模振动力学244页） 得 问题又归结为矩阵的本征值问题，但与原系统的本征值比，矩阵的阶数r小于原系统的阶数n. 禄凑旄百鹈毒垄蘩卧洎案旗郅驰錾砗筒中柃湟 * 4.矩阵迭代法 标准化 整理得： 对于正定系统： 可写为 迭代 假设一个初值 谰味堤领疑亢它妪用擦缯娌銎鸲诂丕皖洪绫老呱殚杼医潮散凄揭岍霜极阻鸦沮蟛钊杖槎剐宏啵鹿丹肮碜角绎咏剥茇萌弼项针 * 解： 取初值： 扬踣魇瘸吮蔽巡甩述捩联稍绯咂泉铀勖俘霁碣锁甭婕胗匍哆忍烃祭蓓潺习虑瞀寓牙啭谖蕲姊祚陛蹬滑觅檩巴 * 振型值趋于稳定。 迭代终止