机械零件的有限元分析.pptVIP

* 第二章 梁单元 如图所示，为一简单平面纯弯曲梁单元（只考虑弯曲，不考虑轴向变形） 位移边界条件（挠度、转角）： 定义位移函数： 将边界条件代入可得： 写成矩阵形式为： 求解上述矩阵方程为： 则位移函数可表示为： 因此，定义形函数为： 简化后得： 则位移函数可简化为： 梁单元在局部坐标系统中有： 形函数可变为： 引入一量纲变量： 称为自然坐标 则形函数变为： 那么，位移函数可进一步表示为： 梁的曲率为： 令： 称为应变矩阵 则： 梁的截面弯矩为： 正应力： 梁单元的应变能： 其中： 外力做功为： 所以，梁单元刚度矩阵为： 将 B 代入上式可得： 因此，梁单元刚度方程为： 例题1：如图所示，求中间结点的挠度和转角以及两端点的反力与反力矩。 显然，将整个梁分成两个梁单元。 单元1与2的刚度矩阵分别为： 整体刚度矩阵为： 整体有限元方程为： 载荷和位移边界条件： 简化方程可得： 求解可得： 则反力及反力矩为： 例题2：如图所示，求挠度、转角以及反力。 可以分成两个梁单元，一个弹簧单元进行求解。 已知： 弹簧单元矩阵： 两个梁单元矩阵为： 整体有限元方程为： 边界条件： 边界条件代入有限元方程，简化后求解得： 再将求解结果代入有限元方程，可得反力。 如图所示，为一考虑轴向变形的平面弯曲梁单元。 这种梁单元相当于等截面直杆单元与不考虑轴向变形的纯弯曲梁单元的合成。 其单元刚度矩阵 为两