有限元离散方法二维问题.ppt

若 讨论对象：椭圆型偏微分方程边值问题 转化为变分问题 5）约束边界处理 其中 找一个函数 满足 有限元化 试探函数空间 边界节点编号 若不在开头，总刚矩阵划去边界节点相应的行与列 有限元法内容总结 1 明确微分方程 2 转化为变分形式； 3 有限元离散； 区域剖分；确定单元基函数； 4 转化为代数方程 单元分析；总体合成；边界条件的处理； 5 解代数方程。 √ √ √ 总结：实际计算步骤 1）单元剖分 节点编号和坐标 单元编号与对应节点(i,j,m) 边界条件节点编号 2）有限元离散 基函数—单元号和单元顶点i处参数(ai,bi,ci) 计算单元刚度矩阵、单元荷载向量 合成总刚度矩阵与总荷载向量，形成代数方程组 3）根据边界调整K和F，即“掐头或去尾” 4）求解上代数方程组，得数值解 例题： 分析 解： 平面区域上无热源的定常温度场，两边给出温度值，另两边给出绝热条件 采用三角形剖分，共16个单元，15个节点 其中 1）节点坐标与编号 2）基函数与相应参数 单元 信息 １ ２ ３ … 结点号 (i,j,k) 参数值 (4,1,5) (-0.5,1) (0,-1) (0.5,0) (2,5,1) (0.5,-1) (0,1) (-0.5,0) (5,2,6) … … … 通过计算可知，奇数单元上参数值相同, 偶数单元上参数值相同