有限元方法及应用_01基本概念.pptVIP

sdustzhu 一、有限元示例 几何模型：h1=100mm，h2=50mm，b1=200mm，b2=100mm，厚度t=6mm。 材料特性：材料的弹性模量E=2×1011Pa，泊松比ν=0.25。 边界条件：左端固定，右端受到力F的作用，F=5000N。 问题：对此受力平板进行力学分析。 有限元分析基本步骤 单元划分 单元分析 整体分析 边界条件处理 求解及分析 问题分析 问题类型：弹性力学平面应力问题 边界条件：混合边界条件 初步方案：3节点三角形线性单元 1. 单元划分 2. 单元分析 节点位移 单元位移 单元应变 单元应力 节点受力 单元刚度矩阵 2.1 节点位移 2.2 单元位移 利用节点位移求解待定系数β1～ β6 将系数β1～ β6带入单元位移方程得到节点位移表示的单元内的位移分布 2.3 单元应变 2.4 单元应力 虚位移原理：如果在虚位移产生之前物体处于平衡状态，则当产生虚位移时，外力所做的虚功等于物体的虚应变能。 2.5 节点受力 节点虚位移 节点外力 外力虚功 单元虚位移 单元虚应变 单元应力 单元虚应变能 对于单元厚度为t，面积为A的平面应力问题 2.6 单元刚度矩阵 节点位移与节点受力之间的关系 单元刚度矩阵的特性 对称性 奇异性 主元横正 3 整体分析 例如展开 整体刚度矩阵的直接集成 整体刚度矩阵的特性 对称性 奇异性：须施加边界条件才能求解 主元横正 稀疏性 4 边界条件 应力边界条件的处理 确定各节点受力，消除整体刚度矩阵的奇异性，得到以节点位移为未知数的矩阵方程。 位移边界 条件的处理 5 求解及分析 单元位移 单元应变 单元应力 其它应力：von Mises应力、主应力、最大剪应力等 其它参数：应力梯度、应变及梯度等 可视化输出 其它参数分析及输出 编程计算 输入节点坐标和单元节点编号 输入材料特性 计算单元刚度矩阵 形成整体刚度矩阵 引入边界条件修正整体刚度矩阵和节点力矩阵 计算节点位移 其它参数计算 可视化输出 节点编号 x/m y/m 1 0 0.1 2 0 0 3 0.1 0.1 4 0.12 0 5 0.2 0.05 6 0.2 0 单元编号 节点i 节点j 节点m 1 1 2 3 2 2 4 3 3 3 4 5 4 4 6 5 材料特性：弹性模量E=2×1011Pa，泊松比ν=0.25 输入参数 单元刚度矩阵 k1=1.0e+11 * 8.8000 -4.0000 -2.4000 1.6000 -6.4000 2.4000 -4.0000 8.8000 2.4000 -6.4000 1.6000 -2.4000 -2.4000 2.4000 2.4000 0 0 -2.4000 1.6000 -6.4000 0 6.4000 -1.6000 0 -6.4000 1.6000 0 -1.6000 6.4000 0 2.4000 -2.4000 -2.4000 0 0 2.4000 k2=1.0e+11 * 5.4133 0.6667 -4.9333 0.9333 -0.4800 -1.6000 0.6667 2.2133 1.7333 -0.9333 -2.4000 -1.2800 -4.9333 1.7333 7.3333 -3.3333 -2.4000 1.6000 0.9333 -0.9333 -3.3333 7.3333 2.4000 -6.4000 -0.4800 -2.4000 -2.4000 2.4000 2.8800 0 -1.6000 -1.2800 1.6000 -6.4000 0 7.6800 k3=1.0e+11 * 3.4844 -1.7778 -0.3556 -0.1778 -3.1289 1.9556 -1.7778 5.2178 0.6222 -5.0222 1.1556 -0.1956 -0.3556 0.6222 4.4444 2.2222 -4.088