README.adoc

NAFEMS LE1 Benchmarck

This benchmark is extract from the Abaqus Benchmarks Manual.

Definition

We focus on the LE1 benchmarks in particular.

Running the case

The command line to run this case is

mpirun -np 4 feelpp_toolbox_solid --case "github:{path:toolboxes/solid/NAFEMS-LE1}"

Copy command line to clipboard Copy case option to clipboard

Data files

The case data files are available in Github here

• CFG file - [Edit the file]

• JSON file - [Edit the file]

Geometry

The geometry is given here by :

Boundary conditions

We set:

• \$u_y = 0\$ on DC

• \$u_x = 0\$ on AB

• \$\bar{\bar{\varepsilon}}\cdot\mathbf{n}=1e^7\$ on BC.

Inputs

We have the following parameters:

Table 1. Inputs

Name	Value
\$E\$	\$210\, GPa\$
\$\nu\$	\$0.3\$
\$\rho\$	\$7800\, kg/m^2\$

Outputs

We want to compare the value of \$\sigma_{yy}\$ at the point D. The reference value is \$92.7\, MPa\$.

Results

The value of \$\sigma_{yy}\$ at the point D is \$94.09\, MPa\$ for \$32 000\$ dofs, which is \$1.49%\$ higher than the target.

One possibility to get a more accurate output is to use a mixed formulation, where the stress tensor would also be an unknown.