有限元静力学及动力学分析(第六章)

FiniteElementMethodandANSYS程强有限元分析及ANSYS北京工业大学机电学院二、动力学分析2.1概述2.2动力学有限元分析原理2.3模态分析2.4谐响应分析2.5瞬态分析一、静力学分析1.1静力学分析类型1.2静力学有限元分析步骤第六章结构动静力学分析1.1静力学分析类型静力分析是用来计算结构在固定不变载荷作用下的响应，如位移、应力、应变等，也就是探讨结构受到外力后变形、应力、应变的大小。与固定不变的载荷对应，结构静力分析中结构的响应也是固定不变的。静力分析中固定不变的载荷和响应是一种假定，即假定载荷和结构的响应随时间的变化非常缓慢。一般来讲，静力分析所处理的载荷通常包括位移载荷（如支座位移等）；稳定的惯性力（重力和离心力等）；外部施加的作用力(集中力、面力和体力)；温度载荷（对于温度应变）；能流载荷（对于核能膨胀）。1.2静力学分析步骤基于ANSYS进行静力分析的基本步骤与ANSYS典型分析的过程相同，一般包括建模、加载求解和检查分析结果等3个基本步骤。1．建模2．加载求解3．检查分析结果2.1动力学有限元分析原理1)动力学分析的原因2)动力学有限元分析引例3)动力学分析的定义和目的4)动力学分析类型1)动力学分析的原因静力分析也许能确保一个结构可以承受稳定载荷的条件，但这些还远远不够，尤其在载荷随时间变化时更是如此。著名的美国塔科马海峡吊桥（GallopingGertie）在1940年11月7日，也就是在它刚建成4个月后，受到风速为42英里/小时的平稳载荷时发生了倒塌。美国塔可马吊桥坍塌之谜重庆綦江彩虹桥-新彩虹桥彩虹桥始建于1994年11月5日，竣工于1996年2月16日，垮塌于1999年1月4日，建设工期1年零102天，使用寿命仅两年零222天。共造成40人死亡，其中包括18名年轻武警战士，直接经济损失628万余元。2)动力学有限元分析引例二系悬挂轨道交通车辆浮沉振动.离散系统：质量（包括转动惯量）模型只具有惯性弹簧模型只具有弹性，本身质量忽略不计阻尼模型不具有弹性，也不具有惯性，是耗能元件，相对运动时产生阻力连续系统：弹性体元件组成，典型的有杆、梁、轴、板壳等（1）结构离散2个单元3个节点1m1m2①②23（2）单元分析eeeijieeeejijiejjieitKtFtZtZkkkktFtFtZtZktFtZtZktF)}({][((((())()((())()(((）））））））ekeZj(t)Fj(t)Zi(t)Fi(t)ji（3）整体分析（4）通用运动方程1.结构离散与静力分析相同，选用适当的单元类型将连续的弹性体离散成有限多个单元和节点。3.整体分析利用各节点处的变形协调条件和动力平衡条件即达朗贝尔原理，建立整体刚度方程；)}({)}(]{[)}(]{[)}(]{[tPtKtCtMf2.单元分析从离散的弹性体中任意取出一个单元。利用给定的位移插值方式表示单元内任一点的位移{δ(t)}e,进而确定节点的速度和加速度。3）动力学分析的定义和目的1.什么是动力学分析?动力学分析是用来确定惯性（质量效应）和阻尼起着重要作用时结构或构件动力学特性的技术。2.“动力学特性”分析的目的–寻求结构振动特性（固有频率和主振型）以便更好地利用或减小振动。–分析结构的动力响应特性，以计算结构振动时的动力响应和动位移的大小及其变化规律。4）动力学分析类型（1）模态分析（2）谐响应分析（3）瞬态动力学分析（4）谱分析（1）模态分析请看下面的一些例子:–在工作中，汽车尾气排气管装配体的固有频率与发动机的固有频率相同时，就可能会被震散。那么，怎样才能避免这种结果呢?–答案：进行模态分析来确定结构的振动特性}{}]{[}]{[}]{[tPKCM–模态分析确定系统的固有特性，即与外载荷，阻尼无关。则Pf（t）=0，[C]=[0]；}0{}]{[}]{[KM什么是模态分析?模态分析是用来确定结构的振动特性（固有频率和振型）的一种技术。模态分析的好处：–使结构设计避免共振或以特定频率进行振动（例如扬声器）；–使工程师可以认识到结构对于不同类型的动力载荷是如何响应的；建议：在准备进行其它动力分析之前首先要进行模态分析。M1-28（2）谐响应分析}{}]{[}]{[}]{[tPKCM谐响应分析：用于确定横幅变频简谐激励Pf（t）下的响应。}sin{}]{[}]{[}]{[0taKCMt什么是谐响应分析？确定一个结构在已知频率的正弦（简谐）载荷作用下结构响应的技术。为什么要作谐响应分析？–确保一个给定的结构能经受住不同频率的各种正弦载荷（例如：以不同速度运行的发动机）；–探测共振响应，并在必要时避免其发生（例如：借助于阻尼器来避免共振）。21~t（3）瞬态动力学分析–一个网球排框架应该设计得能承受网球的冲击，但会稍稍发生弯曲.–解决办法：进行瞬态动力学分析来计算结构对随时间变化载荷的响应.}{}]{[}]{[}]{[tPKCM列车交会瞬态冲击压力图4轨道车辆门体位移和应力瞬态响应曲线什么是瞬态动力分析?它是确定随时间有确定变化关系的载荷作用下结构响应的技术；输入数据：作为时间函数的载荷输出数据：随时间变化的位移和其它的导出量，如：应力和应变。（4）谱分析–位于地震多发区的房屋框架和桥梁应该设计应当能够承受地震载荷要求.–太空船和飞机的部件必须能够承受持续一段时间的变频率随机载荷。–喷气发动机推力、火箭发动机振动响应解决办法：进行谱分析来确定结构对地震载荷的影响。谱分析用谱值-频率分析替代时间－历程分析，主要用于确定结构对随机载荷或随时间变化载荷(如地震、风载、海洋波浪、等等)的动力响应情况。变幅变频载荷下的幅频响应分析。第二节模态分析1.术语和概念2.模态分析方法3.ANSYS模态分析步骤4.模态分析实例通用运动方程：假定为自由振动并忽略阻尼：假定为谐运动u=u0cos(t)：1.这个方程的根是i，即特征值，i的范围从1到自由度的数目，特征值的平方根是i，它是结构的自然圆周频率（弧度/秒），并可得出自然频率fi=i/2π。2.相应的向量是{u0}，即特征向量，质点振幅列矩阵，表示结构以频率fi振动时的形状，称为振型。1.基本方程和术语(002uMK0uKuM(tFuKuCuM2.模态分析方法目前常用的求解方法有广义雅可比法、逆迭代法和子空间迭代法。（1）瑞商法（2）子空间迭代法（1）瑞商法(002uMK2010}u]{[][}u{MK方程0uKuM（2）子空间迭代法1.子空间迭代法思路通过选取m个n维向量{δi}，线性叠加为猜想振型{δ0}，从而将计算n维空间的特征值问题，转化为计算其m维子空间上的特征值问题，具体求解时采用迭代法。2.子空间迭代法的求解步骤–猜想m个向量{δi}构造初始振型{δ0}–确定初始固有频率ω0–确定迭代振型{δ0}1–确定迭代固有频率ω1–终止判断|（ω1-ω0）/ω0|e}]{[}{}]{[}{}]{[][}{0000200110MKMKTT3.ANSYS模态分析步骤前处理：定类型，画模型，设属性，分网格。求解：添约束，加载荷，查错误，求结果。后处理：列结果，绘图形，显动画，下结论。固有频率要输密度均匀网格无需添加选模态设阶次算结果列固有频率读结果绘振型制作振型动画不添约束时，前6阶振型为刚体位移，固有频率均为0.第三节谐响应分析1.定义和目的2.术语和概念3.谐响应分析求解方法4.谐响应分析步骤5.实例-弹簧质量系统谐响应分析1、定义与目的1、定义与目的1、引例-货车有阻尼受迫振动(1)力学模型(2)数学模型(3)振动方程解设波形线路简化为正弦曲线，受力分析可建立振动力学模型(图)。(1)力学模型taZpksin波形线路的简谐激励:式中，ωp=2πV/L为波形线路对车轮荷重系统的激振圆频率，V为列车速度，L为波形线路的波长，a为波形线路的波幅，高速线路a＝3～5mm。M1-69(2)数学模型由牛顿第二定律或达朗贝尔原理可建立振动方程。式中P为车体重量；M为车体质量；g为重力加速度；k为弹簧刚度；Z为车体位移；f0位弹簧静挠度。因为静态时，车体处于静力平衡状态，即P=Mg=kf0。(3)振动方程解增幅系数(受迫振动与轨面正弦曲线波幅的比值)由增幅系数与频率比η的关系曲线可见：在o＜η＜1之间，η大，即运行速度高，则受迫振动振幅大；在η＞1时，η大，即运行速度高，则受迫振动振幅小。无阻尼时，当η＝1，即自振频率等于激扰频率时，受迫振动振幅趋于无穷大，这就时共振。共振时的机车车辆速度称为共振临界速度Vc。VmaxVc，机车车辆处于亚临界速度，VVc，机车车辆处于超临界速度运行。在很小的η下，液压减振器的作用不明显；在共振时（η=1），受迫振动的振幅受到明显的抑制。这说明液压减振器除用来衰减固有振动外，还用来控制共振时的振幅。2222224)1(41DDaZ无阻尼时D=0，则211aZ式中，η为频率比，η=ωp/ω，线路激励频率ωp与结构固有频率ω之比。D为相对阻尼率，等于实际阻尼系数C与临界阻尼系数qc之比。ωp=2πV/Lc2qMK2、基本方程和术语Fmax=载荷幅值=载荷函数的相位角F1=实部,FmaxcosF2=虚部,FmaxsinFuKuCuMtitiititiieuiueeuueFiFeeFF)()(21max21max)())((21212FiFuiuKCiM通用运动方程：[F]矩阵和{u}矩阵是简谐的，频率为:谐响应分析的运动方程：umax=位移幅值=位移函数的相位角u1=实部,umaxcosu2=虚部,umaxsin1i3、ANSYS谐响应分析步骤前处理：定类型，画模型，设属性，分网格。求解：添约束，加载荷，查错误，求结果。后处理：列结果，绘图形，显动画，下结论。频率影响要输密度加谐载荷选谐响应设分析法定加载法列固有频率读结果绘振型制作振型动画求解简谐运动方程的三种方法：完整法–为缺省方法，是最容易的方法；–使用完整的结构矩阵，且允许非对称矩阵（例如：声学矩阵）。缩减法*–使用缩减矩阵，比完整法更快；–需要选择主自由度，据主自由度得到近似的[M]矩阵和[C]矩阵。模态叠加法**–从前面的模态分析中得到各模态；再求乘以系数的各模态之和；–所有求解方法中最快的。（1）选择谐响应分析建模选择分析类型和选项输入求解器，选择谐响应分析；主要分析选项是求解方法-在后面讨论；（2）选择分析选项求解方法：完整法、缩减法和模态叠加法。缺省为完整法；自由度输出格式：幅值相位或实部虚部（3）选择分析阻尼从-阻尼、-阻尼和阻尼率中选取阻尼率最常用（4）施加谐波载荷实部虚部F1maxF2max所有施加的载荷以规定的频率（或频率范围）简谐地变化规定谐波载荷时要包括：–振幅和相角–频率–阶梯载荷对线性变化载荷的说明注意：ANSYS输入振幅和相角或规定实部和虚部分量；例如，施加两个简谐力F1和F2，其相角相差:F1real=F1max(F1的振幅)F1imag=0F2real=F2maxcosF2imag=F2maxsin谐波载荷的频率：通过频率范围和在频率范围内的子步数量来规定每秒的循环次数（赫兹）；例如，在0-50频率范围内有10个子步时将给出在5，10，15...45和50Hz等频率上的解；而同一频率范围只有一个子步时，则只给出50Hz频率上的解。载荷类型：阶梯载荷（Stepped)和线性变化载荷（Rapped)。谐波载荷通常是阶梯