节直梁的有限元分析

李建宇天津科技大学有限元分析FiniteElementAnalysis内容有限元法的直接刚度法1直梁的有限元分析要求了解：单元、节点的概念理解：梁单元节点内力、节点载荷的概念

梁单元刚度矩阵的概念掌握：梁单元刚度矩阵的建立办法利用节点平衡组装整体刚度作业超静定悬臂梁的有限元分析2有限元分析2-1上节回顾有限元分析的基本思想

整体离散标准单元分析单元组装整体求解复杂工程结构如何离散？

分析什么？

如何组装？

求解什么？

复杂结构离散化途经

上节回顾自然离散逼近离散单元

直梁的离散化直梁的单元分析直梁的整体组装本节提要一个例题：MDCBAZ已知E、I、Z、M,AB=BC=CD=L,IAC=2I,ICD=I求:(1)A、D端约束反力；

(2)C处的挠度和转角。Ansys演示…E=3.0e11I=6.5e-7,A=6.8e-4L=1,M=720,Z=3000

一、直梁的有限元模型

211节点：1，2单元编号1223节点：2，32单元编号433节点：3，43单元编号节点(node):1,2,3,4单元(element):1,2,3ZM4321123Z划分单元的原则（设置节点的原则）几何形状发生改变处外载荷规律发生改变处（含约束）边界点计算关心的位置单元尺寸要均匀M4321123

二、单元分析M4321123ijeqimiqjmj截面法：单元节点内力向量：

二、单元分析单元节点位移向量：jefifji单元变形：

二、单元分析jeqimiqjmj单元受力：jefifji单元变形：力与变形的关系：由小变形、线弹性假设，得线性函数：

二、单元分析线性函数：矩阵表达：梁单元刚度矩阵Ke

二、单元分析单元刚度方程：符号名称：

：单元节点内力向量：单元节点位移向量：单元刚度矩阵

二、单元分析:刚度矩阵元素的力学含义考察a11,

a21,a31,a41的含义，令a11：固定j节点，使i节点发生挠度为1，转角为0的节点内力qi大小a21，a31，a41?ij

二、单元分析:刚度矩阵元素的力学含义a11a21a31a41的计算：a11：固定j节点，使i节点发生挠度为1，转角为0的节点内力qi大小a21：固定j节点，使i节点发生挠度为1，转角为0的节点内力mi大小ijqimi由材料力学知识对悬臂梁分析，可得ij考察a12,

a22,a32,a42的含义，令：a12：固定j节点，使i节点发生转角为1，挠度为0的节点内力qi大小a22，a32，a42?

二、单元分析:刚度矩阵元素的力学含义

二、单元分析:刚度矩阵元素的力学含义a12a22a32a42的计算：a12：固定j节点，使i节点发生转角为1，挠度为0的节点内力qi大小a22：固定j节点，使i节点发生转角为1，挠度为0的节点内力mi大小由材料力学知识对悬臂梁分析，可得ijqimi

二、单元分析单元刚度矩阵的具体表达对称矩阵

二、单元分析单元分析小结：单元变形分析单元内力分析单元刚度方程单元分析的任务：单元内力用单元变形表达

三、单元组装利用单元分析技术，对各单元有：单元1：21q11m11q21m211单元2：32q22m22q32m322单元3：43q33m33q43m433

三、单元组装组装原理：位移协调条件节点平衡条件位移协调条件：各单元共享节点位移相等

三、单元组装节点平衡条件：各节点满足平衡条件M4321123节点1Z1M1m11q11节点2节点4Z4M4m43q43Z2M2m22q22q21m21Z3M3m33q33q32m32节点3节点1Z1M1m11q11节点1平衡方程：节点2平衡方程：Z2M2m22q22q21m21同理，可列出节点3、4平衡方程。节点4Z4M4m43q43Z3M3m33q33q32m32节点3整体刚度方程

将节点1、2、3、4的平衡方程整合为一个矩阵方程结构整体刚度方程

其中，为整个结构的节点载荷向量（外载、约束力）

为整个结构的节点位移向量

为结构的整体刚度矩阵，也称总刚度矩阵对称，稀疏，奇异，主对角元恒正四、总体刚度方程的求解力的边界条件位移边界条件四、总体刚度方程的求解对称，稀疏，非奇异，主对角元恒正

五、回代求解未知载荷求解得解将节点位移代入各节点平衡方程，求得未知外载荷。节点位移向量本节课小结有限元法求直梁基本过程第一步，对直梁进行离散化，划分为有限个单元。第二步，对各结点和单元进行编码。第三步，进行单元分析，形成单元刚度矩阵。第四步，进行整体分析，形成整体刚度矩阵。第五步，引入边界条件。边界条件的引入可以使问题具有解的唯一性。第六步，求解方程组，计算结构的整体结点位移。第七步，求单元内力，未知外载荷。

作业MC