Skip to content

Commit 2eba581

Browse files
committed
updating ouput metrics in md file
1 parent 2326360 commit 2eba581

2 files changed

Lines changed: 27 additions & 28 deletions

File tree

.vscode/settings.json

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,4 @@
1+
{
2+
"python-envs.defaultEnvManager": "ms-python.python:conda",
3+
"python-envs.defaultPackageManager": "ms-python.python:conda"
4+
}

docs/benchmarks/linear elasticity/plate-with-hole.md

Lines changed: 23 additions & 28 deletions
Original file line numberDiff line numberDiff line change
@@ -94,45 +94,40 @@ B(\boldsymbol u,\boldsymbol v) &= \int_{\Omega} \boldsymbol\varepsilon(\boldsymb
9494
$$
9595

9696

97-
In order to solve the weak formulation, the finite-element method (FEM) can be used. This method discretizes the domain $\Omega$ into so called finite elements that can for example be triangles or quadrilaterals in 2D. On these elements, ansatz functions are defined such that they are continous on the boundaries between elements. These functions form a basis for the solution space for an approximate solution $\boldsymbol{u}_h$ of the problem.
97+
In order to solve the weak formulation, the finite-element method (FEM) can be used. This method discretizes the domain $\Omega$ into so called finite elements that can for example be triangles or quadrilaterals in 2D. On these elements, ansatz functions are defined such that they are continuous on the boundaries between elements. These functions form a basis for the solution space for an approximate solution $\boldsymbol{u}_h$ of the problem.
9898

99-
## Solution metrics
99+
## Output metrics
100100

101-
The approximate solution is compared with the analytical solution using the $L_2$ norm which is defined as
101+
- The approximate solution is compared with the analytical solution using the $L_2$ norm which is defined as
102102

103-
$$
104-
\Vert \boldsymbol{f}\Vert_{L_2} = \sqrt{\int_\Omega \left|\boldsymbol{f}(\boldsymbol{x})\right|^2 \mathrm d \boldsymbol x}
105-
$$
106-
107-
This norm is computed for the error between the finite element solution and the analytical solution of displacements i.e., $\Vert \boldsymbol{u}-\boldsymbol{u}_h\Vert_{L_2}$.
108-
109-
The maximum displacement error is computed at the nodes of the finite element mesh with respect to the analytical solution.
110-
111-
$$
112-
e_{\max}
113-
=
114-
\max_{i \in \mathcal{N}}
115-
\left\|
116-
\mathbf{u}(\mathbf{x}_i)
117-
-
118-
\mathbf{u}_h(\mathbf{x}_i)
119-
\right\|
120-
$$
121-
122-
where $\mathcal{N}$ denotes the set of nodes of the finite element mesh.
103+
$$
104+
\Vert \boldsymbol{f}\Vert_{L_2} = \sqrt{\int_\Omega \left|\boldsymbol{f}(\boldsymbol{x})\right|^2 \mathrm d \boldsymbol x}
105+
$$
123106

124-
Max Von-Mises stress
107+
This norm is computed for the error between the finite element solution and the analytical solution of displacements i.e., $\Vert \boldsymbol{u}-\boldsymbol{u}_h\Vert_{L_2}$.
125108

126-
Displacement at the top-right corner of the plate.
109+
- The maximum displacement error is computed at the nodes of the finite element mesh with respect to the analytical solution.
127110

128-
The reaction force at the left boundary.
111+
$$
112+
e_{\max}
113+
=
114+
\max_{i \in \mathcal{N}}
115+
\left\|
116+
\mathbf{u}(\mathbf{x}_i)
117+
-
118+
\mathbf{u}_h(\mathbf{x}_i)
119+
\right\|
120+
$$
129121

130-
Number of Degrees of Freedom which supports in comparing the outputs between different meshes.
122+
where $\mathcal{N}$ denotes the set of nodes of the finite element mesh.
131123

124+
- Max Von-Mises stress
132125

126+
- Displacement at the top-right corner of the plate.
133127

128+
- The reaction force at the left boundary.
134129

135-
With these metrices, we can perform a convergence analysis for different approximations $\boldsymbol{u}_h$ which differ in the element size $h$. Plotting the error over the used element-size in a log-log plot lets us determine the convergence order of the approximation.
130+
- The number of Degrees of Freedom in the finite element mesh.
136131

137132
## Table of parameters
138133

0 commit comments

Comments
 (0)