结构坐标单元刚度矩阵.ppt

3.7、罚函数法 当结构的某个自由度给定位移值： 当结构的某个自由度给定位移值： 修改后的第 j 个方程为 2、组集结构刚度矩阵基本方法示例 l （1） （2） X Y P ○ ○ ○ ③ ① ② 1 2 3 4 5 6 图3-5 图（3-5）示出了一简单桁架。 有三个节点，编号为：①、②、③ 二个单元，编号为：（1）、（2） 受集中力P作用。 结构共有6个自由度： 6个节点位移， 6个节点力。 结构坐标系如图中XY。 （1）单元刚度矩阵 1 2 ① ○ ○ 3 4 ② （1） X Y 3 4 ③ ② 5 6 ○ ○ （2） 设，对于结构坐标系而言，单元刚度矩阵已知。（这一过程将在后续章节中涉及，此处暂从略） 5 6 3 4 5 6 3 4 （2）结构刚度矩阵 ① 膨胀单元刚度矩阵 1 2 3 4 1 2 3 4 l （1） （2） X Y P ○ ○ ○ ③ ① ② 1 2 3 4 5 6 图3-5 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 l （1） （2） X Y P ○ ○ ○ ③ ① ② 1 2 3 4 5 6 图3-5 * 第二章： 单元分析的主要任务是： 一、组集单元刚度矩阵； 二、组集单元等价节点力矩阵。 以平面应力三角形单元为例： （2-39） （2-40） （2-41） ，组合见式（2-40） 等价节点力 （2-45） （2-46） （2-35） 第三章 系统分析 3.1 关于结构坐标概念 下面，转入对系统（结构整体）进行分析。 结构整体（系统）分析时，必须在统一的坐标系内进行，如图（3-1）所示。 X Y X Y P ○ ○ ○ ○ ○ P 图（3-1） (a) 平面桁架(杆件单元) 悬臂深梁 (平面三角形单元) 理由： 系统分析的基本内容是把被离散各单元集合成结构整体，研究结构整体的平衡和变形协调。显然，这一过程只能在统一的坐标系中完成。此坐标系称为“结构坐标”或“整体坐标”。 X Y X Y P ○ ○ ○ ○ ○ P 然而，单元分析时所取的坐标系只有少数情况才和结构坐标系一致，如常应变三角形单元。许多时候，单元分析时的坐标系和结构坐标不一致。此时，所有与坐标系有关的量：单元位移、单元力、单元刚度矩阵等，都应明确加上“单元坐标”字样给予限定，如单元坐标单元位移、单元坐标单元力、单元坐标单元刚度矩阵。 系统分析前，必须先把单元在单元坐标系中的量转化为结构坐标系中的量。加“结构坐标”字样界定，如结构坐标单元位移 、结构坐标单元力 、结构坐标单元刚度矩阵 等。 本节下面分析中假设单元的所有量都已经是结构坐标系中的量。至于如何从单元坐标转化为结构坐标将在以后各章中针对不同类型单元详述。 3.2 系统能量分析 1、系统应变能 系统应变能（Uf）是各单元应变能之和： （3-1） ??? 、 [k] 分别是结构坐标单元位移列向量、结构坐标单元刚度矩阵。 [?]1×n1——系统节点位移向量组成的行阵。 3、系统总势能? 2、系统外力势能V {P}n1×1——作用在系统节点上的外荷载向量 （3-2） 1、变换（3-3）中的应变能部分 作如下2项变换： 3.3 结构刚度矩阵及其组集 （3-3） 令各单元位移 中的元素和系统节点位移 中的相应元素位置相同、大小相等，如图3-2中的?1与①单元?3和 ②单元?1相等，以保证将各离散单元集合成结构系统后在系统结点处的变形满足变形协调条件。 ?3 l/2 3 1 ① ?1 ?2 ?4 l/2 1 2 ② ?1 ?2 ?3 ?4 ① ② ③ 图3-2 ?1 ?2 ?3 ?4 1 2 3 4 x y ?5 ?6 ?7 ?8 （1）变换单元位移列阵??? 实施步骤： 建立一个与结构节点位移列阵体积一样的 n1 ×1 阶列阵 （膨胀后的单元位移列阵）； 将单元位移列阵 元素放入 中的结构节点位移 相应的位置； 再将元素换成结构的节点位移。 1 2 3 4 5 6 7 8 ?3 (1) ?2 (1) ?1 (1) ?4 (1) ?1 ?2 ?