偏微分方程的有限元法.ppt

HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法例5.2.1一个边长为1的二维正方形静电场域，电位函数为φ(x,y)，边界条件如图所示，试用有限元法确定二维静电场域的电位分布。解：该二维静电场域的电位函数φ(x,y),可以用下列第一类边界条件的偏微分方程描述：HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法按照右图进行三角形单元剖分，单元编号按照从左到右，从下到上的顺序编号。节点编号：1(0,0)2(0,1)3(0.5,0.5)4(1,0)5(1,1)三角形单元编号：e(j,k,l)单元内顶点按逆时针编号1(1,3,2)2(1,4,3)3(2,3,5)4(3,4,5)第一类边界条件HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法对于三角元1(1,3,2)HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104扩展到全部节点5.2基于变分原理的有限元法HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法同样，对于三角元2(1,4,3)HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104扩展到全部节点5.2基于变分原理的有限元法HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法同样，对于三角元3(2,3,5)HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104扩展到全部节点5.2基于变分原理的有限元法HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法同样，对于三角元4(3,4,5)HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104扩展到全部节点5.2基于变分原理的有限元法HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104全部节点K5.2基于变分原理的有限元法HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法添加标题有限元方程01添加标题可以采用迭代法和直接解法，求解此线性代数方程组。02HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法迭代法：迭代公式为代入初值HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104直接解法：需要处理第一类边界条件5.2基于变分原理的有限元法HarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法A差分法得到差分递推公式BHarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/1045.2基于变分原理的有限元法差分网格A利用边界条件BHarbinInstituteofTechnologyYangkunkyang@hit.edu.cn*/104国外软件ANSYS、ADINA、A