二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现

有限元法求解二维 Poisson 方程的 MATLAB 实现陈莲a ，郭元辉b ，邹叶童a ( 西华师范大学 a．数学与信息学院; b．教育信息技术中心，四川南充 6437009) 摘要: 文章讨论了圆形区域上的三角形单元剖分、有限元空间，通过变分形式离散得到有限元方程．用 MATLAB 编程求得数值解，并进行了误差分析．关键词: Poisson 方程; 有限元方法; MATLAB 编程; 三角形单元剖分中图分类号: O241． 8 文献标识码: A 文章编号: 1009 －4970( 2018) 05 －0015 －04 0 引言热学、流体力学、电磁学、声学等学科中的相关过程，都可以用椭圆型方程来描述．最为典型的椭圆型方程就是泊松方程［1］．泊松方程在许多工程领域中有着广泛的应用，如它可以用来模拟半导体器件的静电过程．在实际应用中，求解区域往往是不规则的，因此要想求出问题的解析解是非常困难的，很有必要研究其数值计算方法．求解偏微分方程数值解常用的方法为有限差分法和有限元方法．有限元方法是由 20 世纪 50 年代初的工程师们提出，其最初的数学思想是康托提出在三角形网格上用逐片线性函数去逼近 Dirichlet 边值问题［2］．有限元方法［3 －4］是基于传统的Ｒitz －Galerkin 方法和有限差分方法优点的一种数值方法，相对有限差分方法它更擅长处理各种复杂区域，是求解偏微分方程数值解的一种有效方法．基于此，结合 MATLAB 强大的数值计算能力和绘图功能，本文用有限元方法求解圆形区域上的泊松方程，通过 MATLAB 编程实现求解过程．讨论二维 Poisson 方程－ Δu = f( x，y) ， ( x，y) ∈Ω， u =0， ( x，y) ∈Ω{ ， ( 1) 其中 Δu = 2u x2 + 2u y2 ，f( x，y) ∈L2( Ω) ，Ω 为Ｒ2 中的有界凸区域，区域 Ω = x，( )y x2 + y2{ }＜1 ． 1 二维 Poisson 方程的有限元法 1． 1 有限元方法的基本原理和步骤有限元法是基于变分原理和剖分技术的一种数值计算方法，把微分方法的定解问题转化为求解一个等价的“变分问题”．具体步骤如下． 1) 根据微分方程的边界条件给出对应定解问题的泛函及其等价的变分问题; 2) 选定单元的性状，对求解区域进行单元剖分，构造单元基函数; 3) 对多元函数的泛函求极值，形成有限元方程; 4) 处理边界插值函数，求解有限元方程; 5) 收敛性及误差分析． 1． 2 Poisson 方程的变分形式对 u ∈ H2 0( Ω) ， v ∈ H1 0( Ω) ，在方程－ Δu = f( x，y) 两边同乘以v，在区域Ω上积分，由 Green 积分公式［5］得，  Ω u x v x + u y v  ( ) y dxdy =  Ω fvdxdy ( 2) 收稿日期: 2017 － 09 － 18 基金项目: 西华师范大学创新团队项目( CXTD2015 － 2) ; 西华师范大学英才基金( 17YC180) 作者简介: 陈莲( 1993—) ，女，四川德阳人，硕士．研究方向: 数值计算．通讯作者: 郭元辉( 1962—) ，男，重庆巴县人，博士，教授，硕士生导师．研究方向: 工程模拟与数值计算方法． ·51· 2018 年 5 月第 37 卷第 5 期洛阳师范学院学报 Journal of Luoyang Normal Unive

二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现相关推荐

galerkin有限元法matlab实现,有限元法求解二维Poisson方程的MATLAB实现
有限元法求解二维Poisson方程的MATLAB实现陈莲a,郭元辉b,邹叶童a [摘要]文章讨论了圆形区域上的三角形单元剖分.有限元空间,通过变分形式离散得到有限元方程. 用MATLAB编程求得数值 ...
泊松方程 matlab,MATLAB编程求解二维泊松方程
<MATLAB编程求解二维泊松方程>由会员分享,可在线阅读,更多相关<MATLAB编程求解二维泊松方程(3页珍藏版)>请在人人文库网上搜索. 1. 真解 u=sin(pi*x) ...
【生产优化】基于matlab遗传算法求解帐篷工序问题【含Matlab源码 2145期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源: [生产优化]基于matlab遗传算法求解帐篷工序问题[含Matlab源码 2145期] 点击上面蓝色字体,直接付费下载,即可. 获取代码 ...
MATLAB指令求解二维PDE三维PDE
偏微分方程工具箱用有限元法求解偏微分方程. 偏微分方程工具箱可以解椭圆型.双曲型.抛物型.以及本征值方程. PDE工具箱解二维PDE,不仅可以用APP,还可以用指令,而且指令还能解三维PDE. 解题步 ...
一节双曲型方程基于MATLAB的求解,二维双曲型方程的分组并行格式及其数值实验...
第 28卷第 2期 2010年 6月湖北民族学院学报 (自然科学版 ) Journal of Hubei University for Nationalities(Natural Science E ...
matlab迭代求解泊松方程,MATLAB编程求解二维泊松方程
%%%% 真解 u=sin(pi*x)*sin(pi*y) %%% %%%% 方程 -Laplace(u)=f %%%%%% %%%% f=2*pi^2*sin(pi*x)*sin(pi*y) %%% ...
matlab有限元法求解转子固有频率,多支承转子系统辛空间传递矩阵法及应用
[1] 孟光. 转子动力学研究的回顾与展望[J]. 振动工程学报, 2002, 15(01): 5-13. Meng Guang. Retrospect and prospect to the re ...
【指派问题】基于matlab遗传算法求解指派优化问题【含Matlab源码 2292期】
⛄一.遗传算法求解指派优化问题简介 1 遗传算法 1.1 遗传算法简介遗传算法是模拟生物在自然环境中的遗传和进化过程而形成的一种自适应全局优化概率搜索算法,它是一种多学科融合交叉的产物.遗传算法通过 ...
matlab求洛伦兹方程的解,[转载]用Matlab求解洛伦兹方程
1. 洛伦兹方程求解本文说明用Matlab工具箱求解洛伦兹方程的过程,并给出吸引子的三维动态图象.洛伦兹方程如下: (1)这是一个自洽的方程组,求解过程如下: (1) 建立自定义函数 functio ...

二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现

二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现相关推荐

最新文章

热门文章