StructuralMind

Standard tests, NAFEMS

Updated:

On this page

Results.

ExampleReferenceRecordsLargest deviation
NAFEMS LE1, elliptic membranePublished source82.42%
NAFEMS FV12, free thin square platePublished source60.906%
NAFEMS FV32, cantilevered tapered membranePublished source60.724%

Problem. The NAFEMS problems are standard tests used to compare finite element programs. LE1 measures a stress concentration in an elliptic membrane, FV12 the natural frequencies of a free thin plate and FV32 those of a tapered membrane cantilever.

Reference document. The reference values are the target values published by NAFEMS. The Standard NAFEMS Benchmarks (1990) is used for LE1, and Selected Benchmarks for Natural Frequency Analysis (1987) for FV12 and FV32.

Numerical model. LE1 is a quarter model solved in plane stress, and the stress at point D is read at the corner of the element. In the FV12 plate the in-plane translations and the in-plane rotation are held at every node. The FV32 cantilever is held along its root edge and only in-plane vibration is solved.

NAFEMS LE1, elliptic membrane

Published source. NAFEMS, The Standard NAFEMS Benchmarks (1990), LE1

QuantityConditionReferenceComputedDeviation
σyy\sigma_{yy} at point D [MPa]mesh = 4×892.791.28521.53%
Nyy=σyy tN_{yy} = \sigma_{yy}\,t at point Dmesh = 4×89.279.128521.53%
σyy\sigma_{yy} at point D [MPa]mesh = 8×1692.794.88172.35%
Nyy=σyy tN_{yy} = \sigma_{yy}\,t at point Dmesh = 8×169.279.488172.35%
σyy\sigma_{yy} at point D [MPa]mesh = 16×3292.794.94632.42%
Nyy=σyy tN_{yy} = \sigma_{yy}\,t at point Dmesh = 16×329.279.494632.42%
σyy\sigma_{yy} at point D [MPa]mesh = 32×6492.794.18341.60%
Nyy=σyy tN_{yy} = \sigma_{yy}\,t at point Dmesh = 32×649.279.418341.60%

Note. The stress is read at the corner of the element touching point D. Because the corner value is extrapolated from inside the element, on the finer meshes it stays 1.5 to 2.5% above the target. An independent four-node element solution gives the same values within 1% on the same meshes. Since the model shows section forces, the NyyN_{yy} row is the same corner value multiplied by the thickness and is read in the shell detail window of the element at D.

NAFEMS FV12, free thin square plate

Parameters. mesh = 32×32

Published source. NAFEMS, Selected Benchmarks for Natural Frequency Analysis (1987)

QuantityReferenceComputedDeviation
Mode 4 frequency [Hz]1.6221.621790.0130%
Mode 5 frequency [Hz]2.362.362090.0886%
Mode 6 frequency [Hz]2.9222.927370.184%
Mode 7 frequency [Hz]4.2334.194650.906%
Mode 8 frequency [Hz]4.2334.194650.906%
Mode 9 frequency [Hz]7.4167.39580.272%

Note. As the plate is free, the first three modes are rigid body motions with zero frequency. The comparison starts at the fourth mode, and modes 7 and 8 are a pair with the same frequency.

NAFEMS FV32, cantilevered tapered membrane

Parameters. mesh = 32×16

Published source. NAFEMS, Selected Benchmarks for Natural Frequency Analysis (1987)

QuantityReferenceComputedDeviation
Mode 1 frequency [Hz]44.62344.58570.0835%
Mode 2 frequency [Hz]130.03130.0910.0472%
Mode 3 frequency [Hz]162.7162.7060.0039%
Mode 4 frequency [Hz]246.05246.8520.326%
Mode 5 frequency [Hz]379.9382.650.724%
Mode 6 frequency [Hz]391.44391.9080.120%