飞行器结构动力学第6章工程振动中的数值方法.pptxVIP

;;第6章 工程振动中的数值方法 ;第6章 工程振动中的数值方法 ;6.1 概 述;第6章 工程振动中的数值方法 ;6.2 结构动态特性分析;;;上式称为归一化特征向量。 ;;;? (5) 特征值的分隔性质;;6.2.2 迭代法;;; 在广义特征值问题中，质量矩阵M是对称正定，则一定存在非奇异矩阵;;;;;6.2.4 Sturm序列二分法;;;6.2.5 大、中型特征值问题的求解方法 ;;;;;;;;;;;;;;;第6章 工程振动中的数值方法 ;6.3 多自由度系统的响应分析 ; 模态迭加法是建立在展开定理基础上的一种求解动力响应的近似方法。对于一个n维自由度系统，可以通过模态坐标的解耦，确定模态坐标响应，然后通过线性变换得到物理坐标响应。 ;初始条件为;利用正交关系;用Duhamel积分求解，则强迫响应解为;;模态迭加法的计算步骤可概括如下： ①求解系统特征值问题，形成模态变换矩阵； ②通过模态坐标变换，解耦系统方程； ③求解解耦后单自由度系统的模态坐标响应解； ④进行模态坐标反变换，将模态坐标下的解变换到物理坐标。这一过程也是将各模态坐标下的响应解进行迭加，得到结构系统总的响应解。; 将加速度 和速度 用中心差分法代替;将式（6-87）和式（6-88）分别代人式（6-89）中， 可得; 由式（6-90）可知，只要利用t时刻的状态变量计算出 ，即可直接求解 。该积分过程称为显式 积分过程，有时也称为显式直接积分法。;? 中心差分积分法在计算机上实现的时间积分格式;? B. 对每一时间步长： 1. 计算在时刻t的有效荷载： 2. 求解在时刻 的位移: 3. 如果需要，计算在时刻t的加速度和速度: ;6.3.3 Houbolt法;为了得到t+Δt 时刻的解，考虑t+Δt 时刻的动平衡方程;;? Houbolt积分方法在计算机上实现的时间积分格式;? B. 对于每一个时间步长： 1. 计算在时刻 的有效荷载: 2. 求解在时刻 的位移: 3. 若需要，计算在时刻 的加速度和速度: ;6.3.4 Wilson-θ法; 令 表示时间的增量，其中 ，则加速度可表示为;将 代入上两式有;利用时刻 的平衡方程式得;; ? Wilson-θ积分方法在计算机上的实现;;? B. 对于每一个时间步长??? 1. 计算在时刻 的有效荷载： 2. 求解在时刻 的位移： 3. 若需要，计算在时刻 的加速度和速度： ;6.3.5 Newmark法; 一般取δ=1/2，α=1/4，在这种情况下，它是一种定常平均加速度法（如图6-2），是无条件稳定的。 ; ? Newmark积分法在计算机上的实现;;? B. 对于每一个时间步长： 1. 计算在时刻 的有效荷载： 2. 求解在时刻 的位移： 3. 若需要，计算在时刻 的加速度和速度： ;6.3.6 方法的选择;第6章 工程振动中的数值方法 ;6.4 有限元方法简介;6.4.1 结构的离散化 ;6.4.2 单元特性分析;梁的挠曲线微分方程为;将边界条件代入式（6-120）有; 即梁的位移 是4个结点位移的线性组合，其系数由下式确定;图6-5 梁中四个结点的位移函数;单元内位移的表达式为; 建立了单元的位移模式后，其动能、势能均可用结 点的位移表示。单元的动能为;将式（6-124）代入上式，经积分后得到 ;考虑到式（6-124）得 ;当梁上发生虚位移δqe时，作用在单元上所有外力的虚功为 ;;6.4.3 结构的综合;根据式（6-126）整个系统的动能 为;按照同样的方法可求得总刚度矩阵、总阻尼矩阵、和对应于总结点位移的广义力列阵分别为;下面进一步考虑系统的边界条件。 ;将式（6-144）和式（6-145）代入方程（6-143）中，考虑到q0=0，于是有; 例6-1 如图6-6所示，试用有限元分析长为2l的均匀悬臂梁的固有频率。 ;解：为简单起见，把梁分成两个长度均为l的单元。全梁共有3个结点，取总结点位移列阵q为;;可以确定总质量阵M和刚度矩阵K分别为 ;;;;;;6.4.4 建立分析模型考虑的主要因素 ;? 分析模型的主要影响因素 ;