python实现三角，矩形有限元求解二维泊松方程

已于 2023-04-28 20:15:25 修改
阅读量1.6k 收藏 1
点赞数
分类专栏： 传统数值方法 文章标签： 线性代数
于 2021-01-22 22:26:44 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/113005110
版权
传统数值方法 专栏收录该内容
35 篇文章 18 订阅 ¥199.90 ¥299.90
订阅专栏
本文介绍了使用Python通过C0类型矩形和三角形有限元方法解决二维泊松方程的问题。讨论了有限元的构造、数值结果、Dirichlet边界条件和变系数Dirichlet问题。文章提供了代码实现，并对比了矩形和三角有限元在不同范数下的收敛性。
摘要由CSDN通过智能技术生成

重点强调：可以通过本人知乎主页了解博客内容添加链接描述

问题提出

{ − Δ u = f , x ∈ Ω , u =

了解本专栏 订阅专栏 解锁全文
Galerkin码农选手
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
订阅专栏
python：二维有限元方法求解，泊松方程
走向CTO的路上...
12-13 1468
在求解泊松方程时，二维有限元方法将求解域分割成多个小区域（单元），并在每个单元内使用一个局部函数（有限元形函数）逼近原方程的解。《The Finite Element Method: Linear Static and Dynamic Finite Element Analysis》（作者：Thomas J. R. Hughes）- 这本书介绍了有限元方法在结构力学中的应用，包括泊松方程的求解。计算物理学：二维有限元方法在计算物理学中也有广泛的应用，如量子力学中的薛定谔方程、流体力学中的势流方程等。
三角形单元有限元程序
05-10
有限元三角形单元程序，很经典！希望大家能够互相学习与分享。
离散域下的泊松方程求解（python实现）
qq_57731481的博客
01-07 6849
目录 一、背景 二、原理 1、离散Laplace算子介绍 2、Laplace卷积 3、Possion方程解法介绍 三、验证 四、Python下的算法实现 a、DCT求解 1、定义函数calMSE计算误差Mean Square Error 2、导入原图并记下大小 3、拓展原图 4、卷积并求DCT变换 5、除以分母 6、逆变换、拉回低频并显示。 7、打印误差与耗时 b、DST求解 五、误差与结果 六、源码 一、背景 数字化之后的图...
矩形区域的泊松方程，深度学习模拟差分法
forrestguang的博客
10-13 1028
差分法是核心 假设我们已经对区域进行了对应的网格剖分： xi=i∗δx,i=0,…,M−1x_{i}=i*\delta x,i=0,\ldots,M-1xi​=i∗δx,i=0,…,M−1 yj=j∗δy,i=0,…,N−1y_{j}=j*\delta y,i=0,\ldots,N-1yj​=j∗δy,i=0,…,N−1 考虑差分法的运算过程： -(ui+1,j−2ui,j+ui−1,jδx2+ui,j+1−2ui,j+ui,j−1δy2\frac{u_{i+1,j} - 2u_{i,j}+u_{i-1,j
二维有限元方程matlab,有限元法求解二维Poisson方程的MATLAB实现
weixin_28850485的博客
03-17 3030
有限元法求解二维 Poisson 方程的 MATLAB 实现 陈 莲a ，郭元辉b ，邹叶童a ( 西华师范大学 a． 数学与信息学院; b． 教育信息技术中心，四川南充 6437009) 摘 要: 文章讨论了圆形区域上的三角形单元剖分、有限元空间，通过变分形式离散得到有限元方程． 用 MATLAB 编程求得数值解，并进行了误差分析． 关键词: Poisson 方程; 有限元方法; MATLAB ...
二维泊松_泊松方程_poisson_二维泊松_
10-03
一个用于解决二维泊松方程的求解，用有限元中的三角形剖分
poisson:使用有限元方法求解任意二维域上的泊松方程
06-12
在本主题中，我们将聚焦于如何使用有限元方法（Finite Element Method, FEM）在Python环境下求解任意二维域上的泊松方程，这涉及到数值计算和编程技术。 有限元方法是一种将连续区域划分为离散单元的数值方法，通过...
二维泊松方程：迭代求解二维泊松方程，使用 5 点有限差分模板-matlab开发
05-30
在MATLAB环境中，我们可以利用数值方法来求解二维泊松方程。本教程将深入探讨如何使用5点有限差分模板，结合迭代方法在MATLAB中解决这个问题。 二维泊松方程通常表示为： \[ \nabla^2 u = \frac{\partial^2 u}{\...
用FFT法与SOR法求解二维泊松方程的比较
12-14
### 用FFT法与SOR法求解二维泊松方程的比较 #### 摘要 本文探讨了快速傅立叶变换（FFT）方法与逐次超松弛迭代法（SOR）在解决二维泊松方程方面的应用，并对其精确度、速度等方面进行了全面比较。研究结果显示，优化...
C 代码 求解矩形上的二维泊松方程， 使用有限元法， 和分段二次三角形单元.rar
最新发布
05-27
在本压缩包中，我们关注的是使用C语言解决二维...总之，这个C代码实现了有限元方法求解二维矩形上的泊松方程，通过分段二次三角形单元提高了精度。理解这种算法和代码的细节对于学习数值方法和C编程都具有很高的价值。
有限元矩形网格剖分
04-07
本程序采用MATLAB编写，可以自动进行矩形网格的剖分，以及图形的生成。
泊松方程差分求解
04-13
利用matlab对珀松方程作有限差分计算
galerkin有限元法matlab实现,有限元法求解二维Poisson方程的MATLAB实现
weixin_35225144的博客
03-18 1072
有限元法求解二维Poisson方程的MATLAB实现陈莲a，郭元辉b，邹叶童a【摘要】文章讨论了圆形区域上的三角形单元剖分、有限元空间，通过变分形式离散得到有限元方程. 用MATLAB编程求得数值解，并进行了误差分析.【期刊名称】洛阳师范学院学报【年(卷),期】2018(037)005【总页数】4【关键词】 Poisson方程；有限元方法； MATLAB编程；三角形单元剖分0 引言热学、流体力学、...
python有限元传热求解_Python进行有限元编程-平面应力问题（三节点三角形单元）...
weixin_39688856的博客
12-20 2366
参考书籍是：《有限元方法基础教程》(国际单位制)(第五版)章节为：第6章 建立平面应力和平面应变刚度方程重要的事情继续强调(有限元的基本计算流程)：Step 1: 选择单元类型。Step 2: 选择位移函数。Step 3: 定义应变/位移和 应力/应变关系。Step 4: 推导单元刚度矩阵和方程。Step 5: 组合单元刚度方程得出总体方程并引进边界条件。Step 6: 求解未知自由度。Step ...
二维三角元有限元方法matlab,有限元C++编程实践.doc
weixin_26854555的博客
03-23 629
有限元C++编程实践.doc 基于有限元算法的编程实践学号：2011043010031 姓名：廖校毅电子科技大学 物理电子学院【摘要】本文就有限元算法在电磁场分析中的应用展开研究与实践，从电磁场的 Maxwell 方程出发 ,根据电磁场的边值问题及变分公式建立了有限元方程组。具体在实践中，将这些知识运用到 C++语言和 Matlab中，并将这两种语言有机结合，编程并实现二维 FEM。程序最后通过计...
PINN学习记录（2）
weixin_45805559的博客
11-27 6104
PINN学习记录（2） PINN基于解物理的方程的应用，所以我自己学习了一段时间，参考了网上很多的开源项目，末尾会贴出一些，自己总结了一下思路 解微分方程 1、ODE f′(x)=f(x)f'(x)=f(x)f′(x)=f(x) f(0)=1f(0)=1f(0)=1 网络构造 这里说明一下，之后用nn.module，来解决，这只是建立一个通用网络 import torch import torch.nn as nn import numpy as np class Net(nn.Module):
深度ritz方法求解高维泊松方程
forrestguang的博客
05-03 1075
此时样本数目过于庞大，一定会导致训练特别耗时。为此可以采取的策略是，每训练一定次数，就重新撒点，对于边界点也是类似的。举例子，如果我们在区域内部采集1000个样本点，扩散到每个维度，那么每个维度最多采集了2个点，因为。，这样导致样本数目过小，如果一开始每个维度采集20个点，那么需要采集的样本数目就是。长度因子需要满足在边界上取值为0，在内部区域恒大于0.，利用netg训练边界数据，构建均方误差函数，训练完以后netg满足。表示区域维度，这里我们使用PFNN算法来求解这个问题。即可，这个思想非常简单。
PINN学习记录（1）
weixin_45805559的博客
11-21 6912
PINN学习记录（1） 安装Deepdxe 官方链接安装 pip install deepxde 直接在evns 环境下安装即可 在git hub 上下载zip，然后set.up 链接:deepxde 注意： 默认backend 是tensorflow_compact_v1 我自己用pytorch ,所以需要更改 在C:\Users\admin\.deepxde下更改即可 学习例子 1、buegers equations (完整代码） (部分代码) def pde(x, y): dy
深度学习求解微分方程系列四：一种基于自适应激活函数的PINN求解方法—burger方程逆问题
山里娃的博客
11-10 5203
下面我将介绍内嵌物理知识神经网络（PINN）求解微分方程。首先介绍PINN基本方法，并基于Pytorch的PINN求解框架实现求解含时间项的一维burger方程逆问题。
给我有限元求解二维泊松方程的matlab参考程序
03-31
以下是一个简单的 MATLAB 程序，用于求解二维泊松方程： ```matlab % 定义网格尺寸和边界条件 N = 50; % 网格数量 L = 1; % 区域长度 h = L/N; % 网格尺寸 x = (h/2:h:L-h/2)'; % 网格节点 [X,Y] = meshgrid(x,x); % 网格 f = -8*pi^2*sin(2*pi*X).*sin(2*pi*Y); % 要求解的泊松方程右侧 % 定义边界条件 u = zeros(N,N); u(:,1) = sin(2*pi*x); % 左边界 u(:,end) = -sin(2*pi*x); % 右边界 u(1,:) = 0; % 上边界 u(end,:) = 0; % 下边界 % 构造系数矩阵和右侧向量 A = zeros(N^2,N^2); b = zeros(N^2,1); for i = 1:N for j = 1:N k = j + (i-1)*N; if i==1 || i==N || j==1 || j==N % 边界节点 A(k,k) = 1; b(k) = u(i,j); else % 内部节点 A(k,k) = -4/h^2; A(k,k-1) = 1/h^2; A(k,k+1) = 1/h^2; A(k,k-N) = 1/h^2; A(k,k+N) = 1/h^2; b(k) = f(i-1,j-1); end end end % 求解线性方程组 U = A\b; % 将解向量转换为网格形式 U = reshape(U,N,N); % 绘制解 surf(X,Y,U) xlabel('x') ylabel('y') zlabel('u') ``` 在这个程序中，我们首先定义网格尺寸和边界条件。然后，我们构造系数矩阵和右侧向量，并使用 MATLAB 中的“\”运算符求解线性方程组。最后，我们将解向量转换为网格形式，并绘制解的表面。注意，这个程序中的边界条件是非零的，因此我们需要将边界节点的值插入到系数矩阵和右侧向量中。如果边界条件为零，则无需这样做。 这个程序只是一个简单的例子，实际应用中可能需要更复杂的技术，如多重网格方法、预处理等。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Galerkin码农选手

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值