平面分析单元平衡方程.pptVIP

平面问题有限元分析-单元分析2 平面问题有限元分析-单元分析2 * 平面问题有限元分析-单元分析2 * 4 平面问题有限元分析单元分析2 曹国华 4.1单元平衡方程 4.2单元刚度矩阵及其特点 4.3 等效节点载荷 在平面问题中，物体微元的应变能密度U0可表示为： 则整个物体的应变能为： * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 其中 * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 其中 为二阶子矩阵 * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 将几何矩阵和弹性矩阵代入，可得平面应力问题中的每个子块的矩阵表达式: * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 单元等效节点力 单元外力势能为 体积力的势能 表面力的势能 集中力的势能 * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 总势能泛函为 根据总势能求极值的条件为 单元平衡方程为 * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 写成普通方程 * * 平面问题有限元分析-单元分析2 4.1单元平衡方程 ，厚度为t。 例: 求图3.4所示等腰直角三角形单元的刚度矩阵，设 解 根据1（0，a）、2（0，0）、3（a，0）的坐标可得 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 单元面积 ，将以上各数据代入几何矩阵为 由于 ，则平面应力与平面应变问题的弹性矩阵 相同，为 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 应力转换矩阵，如下 则单元刚度矩阵为 从以上求得的结果可以看出，单元刚度矩阵是对称的。1、3、5列相加等于零，2、4、6列相加也等于零，单元刚度矩阵是奇异的。另外在几何矩阵中所涉及到的元素均为“坐标差量”，而与坐标的具体值无关，从而也说明了单元刚度矩阵不随单元或坐标轴的平行移动而改变。 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 1）获取单元节点信息； 2）计算单元面积 及参数 、 等； 3）计算几何矩阵 中各元素； 4）根据问题类型（平面应力、平面应变）计算弹性矩阵 ； 5）计算应力转换矩阵 ； 6）计算单元刚度矩阵 。 从单元刚度矩阵的计算公式 可以看出，在求解单元刚度矩阵时，可以不必求得应力转换矩阵 ，但在后期求解单元应力时要用到该矩阵，所以在此处一并求出，也可以根据实际情况具体考虑。 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 单元刚度矩阵具有如下特点： 1. 中的每一个元素都是一个刚度系数 假设单元的节点位移为 ，由 ，得到节点力为 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 由此可得： 表示i节点在水平方向产生单位位移时，在节点i的水平方向上需要施加的节点力； 表示i节点在水平方向产生单位位移时，在节点i的垂直方向上需要施加的节点力。选择不同的单元节点位移，可以得到单元刚度矩阵中每个元素的物理含义 因此， 中的每一个元素都是一个刚度系数，表示单位节点位移分量所引起的节点力分量， 中的每一个子块表示单元某个节点位移矢量对单元某个节点力矢量的贡献率。 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 2. 是对称矩阵 是对称矩阵。根据对称性，可以减少计算量和存储量。 3. 是奇异矩阵，即 根据这一性质，当已知单元节点位移时，可以求出单元节点力。反之，由于单元刚阵奇异，不存在逆矩阵，因此，当已知单元节点力时不能求出单元上的节点位移。 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 从物体变形的实际情况来说，单元刚阵奇异是必须的，因为单元除了产生变形外，还会产生刚体位移，即仅仅依靠节点力是无法唯一地确定刚体位移的。事实上，当单元的节点力为零时，单元仍可做刚体运动。 例如，假定单元产生了x方向的单位刚体位移即 并假设此时对应的单元节点力为零，则由 得 可以得到，在单元刚度矩阵中1，3，5列中对应行的系数相加为零，由行列式的性质可知， 。 * * 平面问题有限元分析-单元分析2 4.2单元刚度矩阵及其特点 4. 单元的刚度不随单元或坐标轴的平行移动而改变 单元的刚度取决于单元的形