平面问题有限元法.pptVIP

第六章 弹性平面问题有限元 6.1有限元的基本思想 6.1有限元的基本思想 6.1有限元的基本思想 6.1有限元的基本思想 6.1有限元的基本思想 四、利用变分原理(弱解形式)(加权余量法)形成有限元求解方程组 单元方程 总装 总体方程(首先建立单元方程) 线弹性问题 线性代数方程组(关于结点未知场函数值的方程组) 五、引入位移强制边界条件(消除系数矩阵的奇异性) 六、解线性代数方程组：结点未知场函数值 七、计算应力,应变。由结点位移计算单元的应力,应变 八、其它要求 (进行其他工程上的要求计算) 6.2 平面问题的常应变单元 一、单元形式 整体结点编号：1，2，3，…， i,j,…m,n,…, 单元结点编号：1，2，3， X 编号顺序：逆时针方向， 对应于右手坐标系，次序不能任意。 y 单元内某一点位移 节点位移 二、单元位移模式及插值函数 1.线性位移模式 位移模式 （1） 关于坐标x,y的多项式 待定系数。称之为广义坐标 2.求解广义坐标? 在(1)式中，代入i,j,m结点坐标得到： (a) 及 (b) 首先求解 方程(a) 系数矩阵 这里A是三角形面积 其中， ， ， 2)广义坐标可用结点位移表示，将其代入位移模式，可以得到单元位移场用单元结点位移表示 其中： for i, j, m 又： 代入上式，得： 表示结点i有位移 时，对单元内应变的贡献 即： ** 三结点三角形单元的应变矩阵是常数阵，单元位移确定后，单元内各点应变相等，称常应变单元 2、应力矩阵 ， 称为应力矩阵，它建立了与单元内应力的关系。 又： ， ， 可见 在三结点三角形单元中是常数阵。 四、变分法建立有限元格式 由最小位能原理： 系统总位能 应变能 外力势能 单位体力， 已知边界面力， 为边界线，板厚 t=1 虚功原理：如果内力与外力在物体发生任意位移或虚位移过程，所作的内功和外功相等，则物体上所有的点处于平衡状态。 离散体系的总位能 == 各单元总位能之和 离散系统的总位能还可以表述为： 为单元结点位移，对每一个不同的单元，表示不同的位移，不能简单求和。 ** 这里引入结构总的位移列阵 和单元结点位移自由度转换矩阵 (总结点为n ,实际上就相当于前面杆系结构的单元贡献矩阵) 平面问题：三结点三角形单元有3个结点，6个自由度。 因此，结构总的结点位移列阵有n个结点2n个自由度 显然对于每个单元是不同的，仅与单元结点编号有关，将其代入泛函中， 令单元刚度矩阵： 则总刚为： 单元体积力载荷列阵： ，单元面力载荷列阵： 总体等效结点载荷为： 因为： 因此总位能可以写作： 这个位能泛函的未知量是结点位移 变分原理： 这样就得到了有限元方法的求解方程。 五、单元刚度矩阵 1.单元刚度矩阵的计算 * 对于三结点三角形 应变矩阵， 弹性矩阵， 均为常量阵。t 为单元厚度等于常数。 *这样都可以提到积分号外，剩下的积分为三角形面积。 这里 其中： 所以我们说 对称的，并且 对角线上的子阵对称。 这一点是很容易被证明的。 2.单元刚度矩阵的物理意义 同整体刚度方程一样，对单元应用变分原理得： 或 如令 其他位移为零则有： 单元各结点施加的结点力应满足平衡， 由于 的对称性，所以还应有： 推广到一般情况： 的物理意义( i, j=1,2,3,…6) 当第 j 个自由度上位移为1，其他自由度位移为零，需要在 i 自由度方向施加的结点力的大小。 同样有： 每列 每行 3.单元刚度矩阵的特性 1) 对称性：计算单刚时仅需上三