飞行器结构动学_第7章.ppt

建立在展开定理基础上的一种求解动态响应的近似方法。 对于一个n-DOF系统，可以通过模态坐标的解耦，确定模态坐标响应，然后通过线性变换得到物理坐标响应。 只适用于线性系统 模态迭加法的计算步骤： ①求解系统特征值问题，形成模态变换矩阵； ②通过模态坐标变换，解耦系统方程； ③求解解耦后单自由度系统的模态坐标响应解； ④模态坐标反变换，将模态坐标下的解变换到物理坐标。这一过程也是将各模态坐标下的响应解进行叠加，得到总响应。 7.2.2 直接积分法 将加速度和速度用中心差分法代替 只要利用t时刻的状态变量计算出 ，即可直接求解 Ut+Dt 。 该积分过程称为显式积分过程，有时也称为显式直接积分法。 ? 中心差分积分算法 2 Houbolt法 与中心差分方法的区别是它采用向后差分格式 ? Houbolt法的实现 3 Wilson-θ法 ? Wilson-θ法的实现 4 Newmark法 线性加速度法的推广 采用近似表达式 ? Newmark法的实现 7.2.3 方法的选择 决定因素：载荷、结构、精度、非线性程度、方法稳定性等 7.2.3 方法的选择 决定因素：载荷、结构、精度、非线性程度、方法稳定性等 思想：将连续体动力学问题转化为多自由度系统动力学问题。代数方程代替偏微分方程。 优点 （1）可满足任意（不规则）的边界 （2）可在单元内采用低阶的插值函数（型函数），由于单元很小，因此可用低阶多项式来满足函数的可导性。 分析步骤 结构离散化。把复杂结构分割成若干个彼此之间只在结点处相互连接的单元，每个单元都是一个弹性体 单元分析。单元内位移用结点位移插值函数表示。对每个单元，由位移插值函数和动力学基本原理确定刚度矩阵、质量矩阵及其它特征矩阵。 结构综合。形成系统总刚度（质量、阻尼、载荷）矩阵，施加边界条件，得到系统运动方程。 数值求解。 7.3.1 结构离散化 1）把梁分成S个单元， 2）分别对单元和节点进行编号。 7.3.2 单元特性分析 梁单元坐标原点在单元左端，x轴沿单元轴线。 与单元相联系的坐标系也称为局部坐标系。 设单元长度为l，抗弯刚度为El(x)，单位长度的质量为ρ， 由于单元长度一般较小，El(x) 、ρ均当作常量。 挠曲线方程 单元内位移为 单元动能 7.3.3 结构综合 例：用有限元分析长为2l的均匀悬臂梁的固有频率。 7.3.4 建立分析模型考虑的主要因素 选择合适的单元 一/二/三维、直线边单元、曲线边单元 满足精度条件下，选用简单单元 划分网格时应考虑： 几何形状、边界条件、外力等变化剧烈的部分，网格相对较密； 二维和三维单元在各方向的长度应尽量接近，避免畸形单元； 相邻单元的尺寸不宜相差太大； 选取不同单元组合时，应考虑结点上不同单元自由度的协调； 单元疏密视需要而定，尽量减小计算规模； 应充分利用结构的对称性，以减小计算规模； 单元的力学特性与结构的力学特性相符合。如，板单元模拟长度远大于厚度的结构；梁单元模拟长度远大于截面尺寸的结构等。 7.3 有限元方法简介 7.3 有限元方法简介 7.4 子结构模态综合法简介 在结构静力分析中，对于大型复杂结构问题往往采用子结构技术，即将结构划分为若干个子结构，先进行局部分析，然后综合组集，再作整体分析。这种先局部后整体的分析方法是科学研究的普遍方法。实际上有限元法本身也就是这种分析方法的具体应用。 7.4 子结构模态综合法简介 ? 再通过各子结构的界面连接条件，作第二次坐标变换，消去不独立的模态坐标，即对整个结构的模态坐标进行独立坐标变换，得到一组用独立的各子结构模态坐标组成的描述整个系统运动的独立广义坐标。 7.4 子结构模态综合法简介 分析过程可归纳为两个基本步骤： 7.4 子结构模态综合法简介 解：将两端固定的梁划分为α和β两个子结构，并将每个子结构的物理坐标集U分为内部坐标集Ui与界面坐标集Uj 。对α与β两个子结构，用矢量形式表示的物理坐标为 7.4 子结构模态综合法简介 7.4 子结构模态综合法简介 7.4 子结构模态综合法简介 7.4 子结构模态综合法简介 7.4 子结构模态综合法简介 7.4 子结构模态综合法简介 7.4 子结构模态综合法简介 单元运动方程还不能求解，因为它没有反映各单元间的联系和结构的约束条件，因此必须建立整个系统的运动方程。 令总结点位移向量 给出q与qe的关系矩阵Ai i:单元编号，s:单元数目, Ai:长方阵(行数=单元自由度，列数=系统自由度) （7.3-18） 7.3 有限元方法简介 逞踞蚤持整耽碧踊卤慌王禁巢胺劣浊凌患响狭启靠哆讲钒稻泻媳镀儒勿磁飞行器结构动