有限元分析节直梁的有限元分析精.ppt

《有限元分析节直梁的有限元分析精.ppt》由会员分享，可在线阅读，更多相关《有限元分析节直梁的有限元分析精.ppt（33页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、有限元分析节直梁的有限元分析第1页，本讲稿共33页内容内容 Chp.2 有限元法的直接刚度法有限元法的直接刚度法 2.1 直梁的有限元分析直梁的有限元分析要求要求 了解：了解：单元单元、节点节点的概念的概念 理解：梁单元理解：梁单元节点内力节点内力、节点载荷节点载荷的概念的概念 梁梁单元刚度矩阵单元刚度矩阵的概念的概念 掌握：梁单元刚度矩阵的建立办法掌握：梁单元刚度矩阵的建立办法 利用利用节点平衡节点平衡组装整体刚度组装整体刚度作业作业 超静定悬臂梁的有限元分析超静定悬臂梁的有限元分析2有限元分析有限元分析2-1第2页，本讲稿共33页上节回顾上节回顾有限元分析的基本思想 整体离散标准单元分析单

2、元组装整体求解复杂工程结构如何离散？如何离散？分析什么？分析什么？如何组装？如何组装？求解什么？求解什么？第3页，本讲稿共33页复杂结构离散化途经 上节回顾上节回顾自然离散自然离散逼近离散逼近离散单元单元 第4页，本讲稿共33页1.直梁的离散化2.直梁的单元分析3.直梁的整体组装本节提要第5页，本讲稿共33页一个例题：MDCBAZ已知已知已知已知E E、I I、Z Z、M,AB=BC=CD=L,IM,AB=BC=CD=L,IACAC=2I,ICDCD=I=I求求求求:(1)(1)A A、D D端约束反力；端约束反力；(2)C(2)C处的挠度和转角。处的挠度和转角。处的挠度和转角。处的挠度和转角

3、。第6页，本讲稿共33页 一、直梁的有限元模型 211 1节点：节点：节点：节点：1 1，2 2单元编号单元编号单元编号单元编号 1 12 223节点：节点：节点：节点：2 2，3 32 2单元编号单元编号单元编号单元编号433 3节点：节点：节点：节点：3 3，4 43 3单元编号单元编号单元编号单元编号节点节点节点节点(nodenode):1,2,3,4):1,2,3,4单元单元单元单元(elementelement):1,2,3):1,2,3ZM43211 12 23 3Z第7页，本讲稿共33页划分单元的原则（设置节点的原则）几何形状发生改变处 外载荷规律发生改变处（含约束）边界点 计算

4、关心的位置 单元尺寸要均匀M43211 12 23 3第8页，本讲稿共33页 二、单元分析M43211 12 23 3ije eqimiqjmj截面法：截面法：单元节点内力向量：单元节点内力向量：第9页，本讲稿共33页 二、单元分析单元节点位移向量：单元节点位移向量：je efifji单元变形：单元变形：第10页，本讲稿共33页 二、单元分析je eqimiqjmj单元受力：单元受力：je efifji单元变形：单元变形：力与变形的关系：力与变形的关系：力与变形的关系：力与变形的关系：由小变形、线弹性假设，得由小变形、线弹性假设，得由小变形、线弹性假设，得由小变形、线弹性假设，得线性函数：线性

5、函数：第11页，本讲稿共33页 二、单元分析线性函数：线性函数：矩阵表达：矩阵表达：梁单元刚度矩阵梁单元刚度矩阵梁单元刚度矩阵梁单元刚度矩阵K Ke e第12页，本讲稿共33页 二、单元分析单元刚度方程：单元刚度方程：符号名称：符号名称：单元：单元节点内力节点内力向量向量 ：单元：单元节点位移节点位移向量向量 ：单元：单元刚度矩阵刚度矩阵第13页，本讲稿共33页 二、单元分析:刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义考察考察考察考察a a11 11,a a2121,a,a3131,a,a4141的含义，令的含义，令的含义，令的含义，令a11：固定

6、固定固定固定j j节点，使节点，使节点，使节点，使i i节点发生挠度为节点发生挠度为节点发生挠度为节点发生挠度为1 1，转角为，转角为，转角为，转角为0 0的节点内力的节点内力的节点内力的节点内力q qi i大小大小大小大小a21，a31，a41?i ij j第14页，本讲稿共33页 二、单元分析:刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义a a11 11 a2121 a a31 a a4141的计算：的计算：a11：固定固定固定固定j j节点，使节点，使节点，使节点，使i i节点发生挠度为节点发生挠度为节点发生挠度为节点发生挠度为1 1，转角为，

7、转角为，转角为，转角为0 0的节点内力的节点内力的节点内力的节点内力q qi i大小大小大小大小a21：固定固定固定固定j j节点，使节点，使节点，使节点，使i i节点发生挠度为节点发生挠度为节点发生挠度为节点发生挠度为1 1，转角为，转角为，转角为，转角为0 0的节点内力的节点内力的节点内力的节点内力mmi i大小大小大小大小i ij jq qi immi i由材料力学知识对悬臂梁分析，可得由材料力学知识对悬臂梁分析，可得由材料力学知识对悬臂梁分析，可得由材料力学知识对悬臂梁分析，可得第15页，本讲稿共33页i ij j考察考察考察考察a a12 12,a a2222,a,a3232,a,a

8、4242的含义，令：的含义，令：的含义，令：的含义，令：a12：固定固定固定固定j j节点，使节点，使节点，使节点，使i i节点发生转角为节点发生转角为节点发生转角为节点发生转角为1 1，挠度为，挠度为，挠度为，挠度为0 0的节点内力的节点内力的节点内力的节点内力q qi i大小大小大小大小a22，a32，a42?二、单元分析:刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义第16页，本讲稿共33页 二、单元分析:刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义刚度矩阵元素的力学含义a a12 12 a a2222 a a3232 a

9、a4242的计算：的计算：a12：固定固定固定固定j j节点，使节点，使节点，使节点，使i i节点发生转角为节点发生转角为节点发生转角为节点发生转角为1 1，挠度为，挠度为，挠度为，挠度为0 0的节点内力的节点内力的节点内力的节点内力q qi i大小大小大小大小a22：固定固定固定固定j j节点，使节点，使节点，使节点，使i i节点发生转角为节点发生转角为节点发生转角为节点发生转角为1 1，挠度为，挠度为，挠度为，挠度为0 0的节点内力的节点内力的节点内力的节点内力mmi i大小大小大小大小由材料力学知识对悬臂梁分析，可得由材料力学知识对悬臂梁分析，可得i ij jq qi immi i第17

10、页，本讲稿共33页 二、单元分析单元刚度矩阵的具体表达对称矩阵对称矩阵第18页，本讲稿共33页 二、单元分析单元分析小结：单元分析小结：单元变形分析单元变形分析单元内力分析单元内力分析单元刚度方程单元刚度方程单元分析的任务：单元内力用单元变形表达单元分析的任务：单元内力用单元变形表达第19页，本讲稿共33页 三、单元组装利用单元分析技术，对各单元有：单元单元1：21 1q11m11q21m211单元单元2 2：32 2q22m22q32m322单元单元3：43 3q33m33q43m433第20页，本讲稿共33页 三、单元组装组装原理：组装原理：位移协调条件位移协调条件节点平衡条件节点平衡条件

11、位移协调条件：各单元共享节点位移相等第21页，本讲稿共33页 三、单元组装节点平衡条件：各节点满足平衡条件M43211 12 23 3节点节点节点节点1 1Z1M1m11q11节点节点节点节点2 2节点节点节点节点4 4Z4M4m43q43Z2M2m22q22q21m21Z3M3m33q33q32m32节点节点节点节点3 3第22页，本讲稿共33页节点节点节点节点1 1Z1M1m11q11节点节点1平衡方程：平衡方程：第23页，本讲稿共33页节点节点2平衡方程：平衡方程：Z2M2m22q22q21m21第24页，本讲稿共33页同理，可列出节点同理，可列出节点3、4平衡方程。平衡方程。节点节点节

12、点节点4 4Z4M4m43q43Z3M3m33q33q32m32节点节点节点节点3 3第25页，本讲稿共33页整体刚度方程 将节点将节点1、2、3 3、4 4的平衡方程整合为一个矩阵方程的平衡方程整合为一个矩阵方程结构整体刚度方程结构整体刚度方程结构整体刚度方程结构整体刚度方程 其中，为整个结构的节点载荷向量（外载、约束力）为整个结构的节点位移向量第26页，本讲稿共33页 为结构的整体刚度矩阵，也称总刚度矩阵为结构的整体刚度矩阵，也称总刚度矩阵对称，稀疏，奇异，主对角元恒正第27页，本讲稿共33页四、总体刚度方程的求解力的边界条件力的边界条件力的边界条件力的边界条件位移边界条件位移边界条件位移

13、边界条件位移边界条件第28页，本讲稿共33页四、总体刚度方程的求解对称，稀疏，非奇异，主对角元恒正第29页，本讲稿共33页 五、回代求解未知载荷求解得解得解得解得解将节点位移代入各节点平衡方程，求得未知外载荷。将节点位移代入各节点平衡方程，求得未知外载荷。将节点位移代入各节点平衡方程，求得未知外载荷。将节点位移代入各节点平衡方程，求得未知外载荷。节点位移向量节点位移向量节点位移向量节点位移向量第30页，本讲稿共33页本节课小结有限元法求直梁基本过程第一步，对直梁进行离散化，划分为有限个单元。第二步，对各结点和单元进行编码。第三步，进行单元分析，形成单元刚度矩阵。第四步，进行整体分析，形成整体刚度矩阵。第五步，引入边界条件。边界条件的引入可以使问题具有解的唯一性。第六步，求解方程组，计算结构的整体结点位移。第七步，求单元内力，未知外载荷。第31页，本讲稿共33页作作 业业MCBAZ2 2、已知、已知、已知、已知E E、I I、Z、M,AB=BC=L有限元法求有限元法求有限元法求有限元法求:(1)(1)A A、B B端约束反力；端约束反力；端约束反力；端约束反力；(2)C(2)C处的挠度和转角。处的挠度和转角。处的挠度和转角。处的挠度和转角。1 1 1 1、自学教材P17P17，例2-12-1。第32页，本讲稿共33页再再 见见 第33页，本讲稿共33页