Lesson2第2章有限元法的直接刚度法-1梁单元汇编.pptVIP

2）将整体刚度矩阵写成分块形式（分块矩阵的阶数等于结构的节点数） （2-36） 3）叠加形成整体刚度矩阵 （2-37） 2.1直梁的有限元分析 2.1.4 求解有限元基本方程 分析图2.1所示梁结构的边界约束条件，可以得出结构的位移边界条件为： ， 载荷边界条件为： ， ， 。 将边界约束条件代入有限元基本方程（2-30），得 （2-39） 对式（2-39）求解，可以求得节点的未知位移值： 将求得的节点未知位移值回代入式（2-38），可求得节点的未知载荷值： 2.1直梁的有限元分析 2.1.5 举例 例：如图2.7（a）所示梁，已知 、 、 、 ，求 截面的转角及约束力。 （a） 直梁模型 图2.7 直梁 2.1直梁的有限元分析 解：（1）划分单元，对单元和节点进行编号，如图2.7 b 所示 b 直梁的有限元分析模型 图2.7 2.1直梁的有限元分析 （2）分析每一个单元，得到单元刚度矩阵，并写成分块形式。 直梁单元刚度矩阵 的公式前面已求得，见式（2-22），由此可以得出1号单元和2号单元的单元刚度矩阵 和 （2-40） （2-41） 2.1直梁的有限元分析 （3）叠加法形成整体刚度矩阵。 （2-42） 2.1直梁的有限元分析 （4）求解有限元基本方程。有限元基本方程为： （2-43） 2.1直梁的有限元分析 分析结构的边界约束条件 ，位移边界条件为： ，载荷边界条件为： 。将结构的边界约束条件代入有限元基本方程（2-43），得 （2-44） 在上式中划掉已知位移所在的行和列，得 （ 2-45） 解之，可得 （2-46） 2.1直梁的有限元分析 将求得的 回代入有限元基本方程（2-44），求得各节点的未知载荷，即约束力为： （2-47） 2.1直梁的有限元分析 2.1.6 总结 通过以上推导、分析和举例，我们可以总结出有限元求解问题的基本过程为： （1）结构的离散化：划分单元，对单元和节点进行编号。 （2）单元分析，得到每一个单元的刚度矩阵，并将单元刚度矩阵写成分块形式。 （3）求出结构的整体刚度矩阵，得到有限元基本方程。 求出结构的整体刚度矩阵，有两种方法：①通过建立节点平衡方程式得到整体刚度矩阵。即对整个结构的每个节点进行受力分析，并列节点平衡方程式，得到整个结构的有限元基本方程，从而得到整体刚度矩阵。②通过叠加法得到整体刚度矩阵。即将单元刚度矩阵直接叠加形成整体刚度矩阵。 （4）分析结构的边界约束条件，求解有限元基本方程。 将位移边界条件和载荷边界条件代入有限元基本方程，并降阶处理（划掉已知位移所在的行和列），求解降阶后的基本方程，求得节点的未知位移值，并将结果回代入没有降阶的有限元基本方程，求得节点的未知载荷值等。 知识点： 直梁和平面刚架的直接刚度法 重点： 梁单元杆和刚架单元的自由度 单元的坐标变换 难点：直接刚度法的计算过程与物理意义 Ⅰ. 关于梁和弯曲的概念 受力特点： 杆件在包含其轴线的纵向平面内，承受垂直于轴线的横向外力或外力偶作用。 变形特点： 直杆的轴线在变形后变为曲线。 梁——以弯曲为主要变形的杆件称为梁。 工程实例 F2 F1 起重机大梁 图1 钻床的摇臂 图5 悬臂梁 2 梁的基本形式 简支梁 外伸梁 纵向对称面 对称弯曲——外力作用于梁的纵向对称面内，因而变形后梁的轴线 挠曲线 是在该纵对称面内的平面曲线。 非对称弯曲——梁不具有纵对称面 例如Z形截面梁 ，因而挠曲线无与它对称的纵向平面；或梁虽有纵对称面但外力并不作用在纵对称面内，从而挠曲线不与梁的纵对称面一致。 本章讨论对称弯曲时梁的内力和应力。 1、梁整体上处于静止的平衡状态 2、研究变形比较小的情况 3、简单载荷下的直梁的扰度、转角计算公式请参考《材料力学》。 对称弯曲时和特定条件下的非对称弯曲时，梁的挠曲线与外力所在平面相重合，这种弯曲称为平面弯曲。 使微段梁有顺时针转动趋势的剪力、外力为正，反之为负；使微段梁产生向下凸变形的弯矩、外力矩为正，反之为负。 符号规定： Q＞0 Q＜0 M＞0 M＜0 剪力、外力的正负符号规定 弯矩、外力与外力偶矩正负符号规定 2.1直梁的有限元分析 a 直梁模型 图2.1 直梁 b 直梁的有限元模型 以直梁为例来说明有限元法的直接刚度法。 如图2.1 a 所示直梁，已知E、I、Z、M，AB BC CD l，IAC 2I，ICD I。 2.1.1划分单元 两个节点之间的杆件构成一个单元，杆件结构的节点可按以下原则选取： 1、杆件的交点一定要选为节点。 2、阶梯形杆截面变化处一定要取为节点。 3、支承点和自由端要取为节点。 4、集中载荷作用处要取为节点。 5、欲求位移的点要取为节点。 6、单元长度不要相差太多。 按照杆件结构划分单元的原则，对图2.1