杆和梁结构的有限元法教学幻灯片讲义.pptVIP

§2.3.1 简单梁单元 再设： 同理，由梁的变形公式和平衡条件可求得刚度矩阵的第二列元素： 梁单元变形 由刚度方程可得： §2.3.1 简单梁单元 同样的方法可以求出其余2列元素，从而求出单元刚度矩阵： 显然，与弹簧和杆单元一样，该梁单元的刚度矩阵具有如下性质： 1）对称性； 2）奇异性； 3）主对角元素恒正。 刚度矩阵求得后，单元特性就完全确定： 单元虚应变能： 对杆单元应用虚功原理，得： 考虑到 的任意性，立刻得到： 这就是刚度矩阵的一般形式，可推广到其他类型的单元。 ——杆单元刚度矩阵 §2.2.1 等截面杆单元 对于上面的杆单元： 与前面直接法得到的公式相同！ §2.2.1 等截面杆单元 （三）关于杆单元的讨论 1）在单元坐标系下，每个节点一个未知位移分量，单元共有 2个自由度。 2）单元刚度矩阵元素的物理意义： 刚度方程中令： 单元刚度方程 §2.2.1 等截面杆单元 所以，单元刚度矩阵的第i（i=1,2)列元素表示当维持单元的第i个自由度位移为１，其它自由度位移为０时，施加在单元上的所有节点力分量。 ３）单元刚度矩阵对称、奇异、主对角元素恒正。 §2.2.1 等截面杆单元 （四）举例 例1：求图示２段杆中的应力。 解：系统分为２个杆单元，单元之间在节点２连接。 单元刚度矩阵分别为： §2.2.1 等截面杆单元 参考弹簧系统的方法，装配系统的有限元方程（平衡方程）： 引入边界位移约束和载荷： 系统平衡方程化为： §2.2.1 等截面杆单元 上述方程组中删除第１，３个方程，得到： 位移解： 单元1应力： 解得： §2.2.1 等截面杆单元 单元2应力： 提示： 1）本例中单元应力的计算采用了材料力学中的方法,与采用有限元单元应力公式 的结果相同。 2）对锥形杆，单元截面积可用平均值。 3）求应力之前需要求出节点位移——有限元位移法。 §2.2.1 等截面杆单元 例题2： 已知： 求：杆两端的支反力 解 §2.2.1 等截面杆单元 §2.2.2 二维杆单元 （一）2-D空间中杆单元（平面桁架） 1-D空间杆单元 2- D空间杆单元 坐标变换 原来1-D空间中的杆坐标系作为局部坐标系 局部 总体 每节点一个dof 每节点2个dof §2.2.2 二维杆单元 1. 节点位移向量的坐标变换： §2.2.2 二维杆单元 向量的坐标变换矩阵为： 显然是正交阵，即： 2. 单元节点位移向量的变换式 或 3. 单元节点力向量的变换式： §2.2 二维杆单元 4.刚度矩阵的坐标变换 局部坐标系下杆单元的刚度方程为： 扩充到4自由度形式： 写成矩阵符号形式： §2.2.2 二维杆单元 利用前面的向量坐标变换式，得： 考虑到变换矩阵的正交性，得： 总体坐标系中的杆单元刚度矩阵为： 用单元刚度矩阵装配系统刚度矩阵的方法与1-D情况相同，按节点号对子块重新排列。 §2.2.2 二维杆单元 5. 单元应力计算： 即： §2.2.2 二维杆单元 （二）例题 平面桁架由2根相同的杆组成（E，A，L）。求： 1）节点2位移 2）每根杆应力 解： 1. 求出每个单元在总体坐标下的刚度矩阵： §2.2.2 二维杆单元 单元1：1-2 §2.2.2 二维杆单元 单元2：2-3 §2.2.2 二维杆单元 2. 将单元1，2的刚度矩阵扩大到系统规模（6阶）叠加得到总刚度矩阵，再列出系统平衡方程： §2.2.2 二维杆单元 3. 引入边界约束和载荷： 则上面6阶有限元方程凝聚为： 4. 解出未知位移： §2.2.2 二维杆单元 5. 按公式计算杆应力： 得到： §2.2.2 二维杆单元 平面内一般梁单元 简单梁单元 （弯曲变形） 三维空间梁单元简介 2.3.1 2.3.2 2.3.3 结构总刚度矩阵及其性质 梁单元的单元特性 梁单元的单元刚度矩阵 离散结构的整体分析 平面刚架的整体分析 单元与节点 局部坐标系下的平面梁单元 单元刚度矩阵的坐标变换 三维空间梁单元刚度矩阵 2.3 梁单元和平面刚架 §2.3.1 简单梁单元 一、离散化，节点位移与节点载荷 对图(a)直梁，根据结构和载荷情况，分为3段，每段为一个单元。单元之间和端点是节点。梁单元节点的物理模型是“焊接”。 梁上任一节点i处有2个位移分量： 挠度 及转角 。 §2.3.1 简单梁单元 一个节点位移用列阵表示为： 称为节点i的节点位移。 对应节点位移分量，梁上任一节点i的载荷也有2项： 横向力 和弯矩 ，称为广义力。 §2.3.1 简单梁单元 梁上若有分布载荷，可近似地等效到节点上。 称为节点i的节点载荷