在 FEM 上应用边界条件

计算科学 有限元 边界条件 计算物理学 微分方程
2021-12-18 00:39:30

我有一个偏微分方程,如下所示。

ddx((1+x)du(x)dx)=0
具有以下边界条件。
u(0)=0,u(3)=10
为了使用 FEM 解决它,我将 PDE 与残差函数相乘,然后将其转换为弱形式。

然后使用全局边界内的三个元素构造全局矩阵。

请记住,要求解矩阵系统,我们需要求解以下方程。

[A][u]=[b]+[D]
[D]对应于包含边界条件信息的矩阵。

由于我们在第一个方程的右手边没有常数,所以我们没有[b]矩阵。

我构建了[D]矩阵,如下图。

D=((1+x)B00(1+x)B)

在这一点之后我该怎么办?应该怎样1+x术语会影响该系统的强加 Dirichlet BCs 吗?

0个回答
没有发现任何回复~