2020-5-29有限元从0开始第九天(有限元的等价性)

继续给出参考资料镇楼
Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007.
2020-5-27有限元从0开始第八天(插值算子)

上次内容简要回忆

在前面博客中给出了有限元局部插值算子, 以及全局插值算子的定义, 例子以及一些相关的性质, 在本节主要接招有限元的等价性.

3.4 有限元的等价性

在全局插值的应用中, 找到插值${\mathcal{I}}_{\mathcal{T}}$ 在需要范数意义下的一致有界性(不依赖于网格) 是非常重要的, 为此需要在所有单元上比较局部插值算子,建立有限元之间的等价性在做这一点的时候是非常重要的.

定义:(仿射等价有限元) 假设 $(K,\mathcal{P},\mathcal{N})$ 为有限元, 记 $F(x)=Ax+b$ 为非退化仿射变换, 称有限元 $\left(\hat{K},\hat{\mathcal{P}},\hat{\mathcal{N}}\right)$ 和有限元 $(K,\mathcal{P},\mathcal{N})$ 之间是仿射等价的, 如果下面几条成立:

• $F\left(K\right)=F\left(\hat{K}\right)$
• $F^{\ast }\hat{\mathcal{P}}=\mathcal{P}$
• $F_{\ast }\mathcal{N}=\mathcal{N}$

其中 $F^{\ast }(\hat{f})=\hat{f}\circ F$ 称为 pull-back, $\left(F_{\ast }\mathcal{N}\right)(\hat{f})=\mathcal{N}\left(F^{\ast }(\hat{f})\right)$ 称为 push-forward.

性质: 仿射等价性为等价关系.

对称性显然

利用 push-forward 以及 pull-pack 的性质, 可以知道传递性显然

下面证明自反性质, 即假设 $(K,\mathcal{P},\mathcal{N})$ 和 $\left(\hat{K},\hat{\mathcal{P}},\hat{\mathcal{N}}\right)$ 是仿射等价的, 则存在非退化仿射变换 $F=Ax+b$, 使得

• $F\left(K\right)=F\left(\hat{K}\right)$
• $F^{\ast }\hat{\mathcal{P}}=\mathcal{P}$
• $F_{\ast }\mathcal{N}=\mathcal{N}$

下面定义 $F$ 的逆:
$$F^{-1}(\hat{x})=A^{-1}(\hat{x}-b)$$
显然, 容易验证有 $F^{-1}(\hat{K})=F^{-1}(K)$
再有:
$$\begin{array}{r}(F^{-1}{)}^{\ast }f(\hat{x})=f\circ F^{-1}(\hat{x})=f(x)\end{array}$$
$$F_{\ast }^{-1}\left(\hat{\mathcal{N}}\right)(f)=\hat{\mathcal{N}}\left((F^{-1}{)}^{\ast }f\right)=\hat{\mathcal{N}}(\hat{f})$$

仿射等价有限元的例子:

• 记 $K$ 为三角形单元, $\mathcal{P}={\mathcal{P}}_1$, N = { K 在 顶 点 处 的 取 值 } \mathcal{N} = \left\lbrace K 在顶点处的取值\right\rbrace N={K在顶点处的取值}.
• 记 $K$ 为三角形单元, $\mathcal{P}={\mathcal{P}}_2$, N = { K 在 顶 点 处 以 及 每 条 边 的 中 点 处 的 取 值 } \mathcal{N} = \left\lbrace K 在顶点处以及每条边的中点处的取值\right\rbrace N={K在顶点处以及每条边的中点处的取值}.

存在某种节点的选取, 使得任意给定次数的 Lagrange 有限元是仿射等价的.

定义: (插值等价的有限元)

称有限元 $(K,\mathcal{P},\mathcal{N})$ 和 $\left(K,\mathcal{P},\tilde{\mathcal{N}}\right)$ 是插值等价的, 如果对一切充分光滑的函数 $f$, 成立:
$${\mathcal{I}}_{\mathcal{N}}f={\mathcal{I}}_{\tilde{\mathcal{N}}}f,$$
其中, ${\mathcal{I}}_{\mathcal{N}}$ 即为前面所定义的插值算子 ${\mathcal{I}}_K$, 即:
$${\mathcal{I}}_{\mathcal{N}}f=\sum _{i=1}^n{N}_i(f){\varphi }_i$$
性质: $(K,\mathcal{P},\mathcal{N})$ 和 $\left(K,\mathcal{P},\tilde{\mathcal{N}}\right)$ 为有限元, 则这两个有限元是插值等价充要条件为 $\mathcal{N}$ 中的,每一个节点变量, 可以由 $\tilde{\mathcal{N}}$ 中的节点变量进行线性表出.

证明: 先证明充分性, 任取 $N_i\in \mathcal{N}$,由条件可以知道, 存在 $c_j,j=1,2,\cdots ,n$, 使得:

$$N_i=\sum _{j=1}^n{c}_j{\tilde{N}}_j$$
下面证明: ${\mathcal{I}}_{\mathcal{N}}f={\mathcal{I}}_{\tilde{\mathcal{N}}}f$ 对于一切充分光滑的 $f$ 成立. 思路是去证明 $N_i\left({\mathcal{I}}_{\tilde{\mathcal{N}}}f\right)=N_i(f)=N_i({\mathcal{I}}_{\mathcal{N}}f)$.

$$\begin{array}{rl}{N}_i\left({\mathcal{I}}_{\tilde{\mathcal{N}}}f\right)& =\sum _{j=1}^n{c}_j{\tilde{N}}_j({\mathcal{I}}_{\tilde{N}}f)=\sum _{j=1}^n{c}_j{\tilde{N}}_j(f)=N_{i}f=N_i({\mathcal{I}}_{\mathcal{N}}f)\end{array}$$
,从而马上可以得到 $N_i({\mathcal{I}}_{\tilde{\mathcal{N}}}f-{\mathcal{I}}_{N}f)=0$, 利用 $i$ 的任意性,以及 $\mathcal{N}$ 的定义, 立马知道: ${\mathcal{I}}_{\mathcal{N}}f={\mathcal{I}}_{\tilde{\mathcal{N}}}f.$

再证明必要性, 假设 $(K,\mathcal{P},\mathcal{N})$ 和 $\left(K,\mathcal{P},\tilde{\mathcal{N}}\right)$ 是插值等价的有限元, 即有:
$${\mathcal{I}}_{\mathcal{N}}f=\sum _{i=1}^n{N}_i(f){\varphi }_i=\sum _{j=1}^n\widetilde{N_j}(f){\varphi }_j={\mathcal{I}}_{\tilde{\mathcal{N}}}f$$
下面去证明
$$N_i=\sum _{\alpha =1}^n{N}_i({\varphi }_{\alpha }){\tilde{N}}_{\alpha }$$
注意到只需证明上式对一切的 { ϕ j } j = 1 n \left\lbrace\phi_{j}\right\rbrace_{j=1}^{n} {ϕj​}j=1n​成立即可, 利用插值的性质:
$$\begin{array}{rl}{N}_i({\varphi }_j)& =N_i({\mathcal{I}}_{\mathcal{N}}({\varphi }_j))=N_i\left({\mathcal{I}}_{\tilde{\mathcal{N}}}({\varphi }_j)\right)\\ & =N_i\left(\sum _{\alpha =1}^n{\tilde{N}}_{\alpha }({\varphi }_j){\varphi }_{\alpha }\right)\\ & =\sum _{\alpha =1}^n{N}_i({\varphi }_{\alpha }){\tilde{N}}_{\alpha }({\varphi }_j).\end{array}$$

插值等价有限元的例子:

• 参考资料第76页给出的两个 Hermite 有限元具有插值等价性.

定义: 假设有限元 $(K,\mathcal{P},\mathcal{N})$ 和 有限元 $\left(\hat{K},\hat{\mathcal{P}},\hat{\mathcal{N}}\right)$ 之间仿射等价, 并且 $\left(\hat{K},\hat{\mathcal{P}},\hat{\mathcal{N}}\right)$ 和 $\left(\tilde{K},\tilde{\mathcal{P}},\tilde{\mathcal{N}}\right)$ 之间插值等价, 则称有限元 $(K,\mathcal{P},\mathcal{N})$ 和 有限元 $\left(\tilde{K},\tilde{\mathcal{P}},\tilde{\mathcal{N}}\right)$ 之间仿射插值等价.

插值插值等价有限元的例子:

• 显然, 根据定义, 所有仿射等价的有限元都是仿射插值等价的. 取 $\left(\left(\hat{K},\hat{\mathcal{P}},\hat{\mathcal{N}}\right)=\left(\tilde{K},\tilde{\mathcal{P}},\tilde{\mathcal{N}}\right)\right)$
• 对 Hermite 有限元选取合适的边以及内点, 可以做到仿射插值等价.
• 所有的 Argyris 元不具有仿射插值等价性.

性质: 假设有限元 $(K,\mathcal{P},\mathcal{N})$ 和 $\left(\tilde{K},\tilde{\mathcal{P}},\tilde{\mathcal{N}}\right)$ 之间具有仿射插值等价性, 则: $\mathcal{I}\circ F^{\ast }=F^{\ast }\circ \mathcal{I}$, 其中 $F$ 是 $K$ 到 $\tilde{K}$ 之间的仿射变换.

后面几天会进行有限元的编程训练, 因此文档可能会停个几天, 等到把程序清理一下再继续写这个文档.