MATLAB 三层框架的屈曲特征值和特征向量求解解析.pptVIP

三层框架的屈曲特征值和特征向量求解 题目： 轴向荷载作用下的3层框架如图所示，除AB、CD横梁长为10米外，其他梁、柱长度均为5米。 底层框架立柱的截面抗弯刚度EI1为8千牛/米，第二层和第三层框架立柱的截面抗弯刚度EI2和EI3为4千牛/米，第一层、第二层和第三层横梁截面抗弯刚度值大小为2千牛/米。 把框架的每根梁和柱作为梁单元进行分析，每个结点含有竖向位移、横向位移和扭转3个自由度。 一、几个基本概念 1. 单元刚度矩阵 将杆端位移和杆端力联系起来的矩阵，称为单元刚度矩阵。 单元干度矩阵中任一元素Kij的物理意义：当单元的第j个结点位移为单位位移而其他结点位移为零时，需在单元第i个结点位移方向上施加的结点力的大小。 2. 刚度矩阵集成和整体刚度矩阵 3. 几何刚度矩阵 梁单元在受到轴向力和弯曲荷载的作用时，刚度系数由于轴向力的存在得到修正。对应于刚度系数Kij的修正系数就是几何刚度Kgij。 它定义为由结构轴向力以及结点j处的单位位移引起的结点i处的力。这些系数可以通过虚功原理进行计算。（推导简要见文章） 4. 坐标转换问题与坐标转换矩阵 将局部坐标系下的表达量转换成整体坐标系下表达的量，或相反，称为坐标转换。 二、综合编程 1. 整体弹性刚度矩阵与整体几何刚度矩阵的建立 单元弹性刚度矩阵和单元几何刚度矩阵 框架结构的离散化 数据化（结构各参数） 整体坐标 结点编码 单元编号 位移编码 弹性刚度矩阵和几何刚度矩阵集成 施加约束条件 求解 2. 定位向量 3. 如何进行边界条件处理 4. 如何求解 普通单元的刚度矩阵 单元几何刚度矩阵 3）单元刚度矩阵由局部坐标系转化为整体坐标系 7. 屈曲特征值的求解 将单刚集结成总刚后，得到弹性屈曲的有限元求解方程： END * 平面框架结构 单元划分 1 2 3 4 5 6 7 8 9 10 11 平面框架结构 三、结构矩阵分析的基本思路 化整为零 （单元分析） 集零为整 （结点力平衡、位移协调） 整体分析：研究整体的平衡条件、平衡方程的组成规律和求解方法。 编制程序：根据矩阵位移法的分析原理，绘制程序运行框图并选择一种计算机语言给予实现，又称为程序设计。 离 散 化：确定座标、单元编码、结点编码（总体码和局部码）、位移编码（总体码和局部码） 单元分析：研究单元的力学特性，建立单元杆端力和杆端位移的关系。 1. 拟解决的问题 ? 2.平面弯曲自由式单元刚度方程 e ? ? 3.单元刚度矩阵的性质 系数的物理意义：单元刚度矩阵是结点力与杆端位移之关系；矩阵的阶数与位移分量数相等； 表示 引起的杆端力 的大小。 单刚一般具有奇异性：受力角度存在刚体位移；数学角度向量相关。 单刚具有对称性：由功的互等定理可知 ? 4. 座标转换问题 在讨论了单元在局部座标系中的刚度方程的基础上，由于体系中的各单元方位不尽相同，要考虑结点的位移协调和受力平衡，应该有一个统一的标准，故引入结构的整体座标系。 两种座标系间物理量的转换过程称座标转换。（Transformation of coordinates） ? 1） 平面自由式单元位移的座标转换 x y e ? 2）座标转换矩阵 5. 整体分析 1）利用定位向量集成结构刚度矩阵 事实上，结构刚度矩阵的元素是由单元刚度矩阵的元素组成的，只要确定了单元刚度矩阵各元素在结构刚度矩阵中的位置，就可以由单元刚度矩阵元素直接集成结构刚度矩阵。 位 移 码：按先起点后终点将杆件位移顺序排列的单元位移序号称单元局部位移码；而按整体结点码的顺序将结点位移顺序排列的位移序号称为整体位移码。 定位向量：由单元局部位移码所对应的结构整体位移码所组成的向量称为单元的定位向量。 2）定位向量 集成的规则 1、先求出单元在整体座标系中的单元刚度矩阵； 2、单元刚度矩阵中元素kij所在的行码 i、列码 j 就分别与单元定位向量对应的结构刚度矩阵中一个位置分量行码 s，列码 r 相匹配； 3、单元刚度矩阵中元素 kij 分别送到结构刚度矩阵中的的对应位置，从而实现按位叠加。 ? 3）结构刚度方程集成举例 采用“对号入座”的方法集成 1(1,2,3) 2(4,5,6) 3(7,8,9) 1 2 FP2 FP1 M ? 1 2 3 1 2 3 ? 结构原始刚度矩阵 ? 6.边界约束的处理 在总刚方程的形成过程中要反映边界约束的影响，这种影响可在形成总刚前引入,亦可在总刚生成后引入，因而有先处理法和后处理法之分 我们发现结构刚度矩阵中包含边界处的已知位移，综合结点荷载列阵中包含未知的支座反力。因此刚度方程还需要做边界处理才能求解； 是后处理法还是先处理法，我们可以在定义单元定位向量的时候予以体现。 ?