Benchmark D1

VMD1
Straight Cantilever Beam Under Unit Load

Overview

Reference:R. H. MacNeal, R. L. Harder, "A Proposed Standard Set of Problems to Test Finite Element Accuracy", Proceedings, 25th SDM Finite Element Validation Forum, 1984.
Analysis Type(s):Static Analysis (ANTYPE = 0)
Element Type(s):
2D 4-Node Structural Solid Elements (PLANE182)
2D 8-Node Structural Solid Elements (PLANE183)
2D 6-Node Structural Solid Elements (PLANE183)
3D 8-Node Structural Solid Elements (SOLID185)
3D 20-Node Structural Solid Elements (SOLID186)
3D 10-Node Tetrahedral Structural Solid Elements (SOLID187)
Input Listing:vmd1.dat

Test Case

A straight cantilever beam, fixed at one end, is subjected to a unit load. Determine the displacement at the end of the beam for unit loads including extension, in-plane shear, out-of-plane shear, and twist (where applicable). Examine the influence of rectangular, trapezoidal, and parallelogram element shape models on tip displacement and the percent energy error norm.

Figure 604: Straight Cantilever Beam Problem Sketch

Straight Cantilever Beam Problem Sketch

Material PropertiesGeometric PropertiesLoading
E = 10 x 106 psi
υ = 0.3
L = 6 in
d = 0.1 in
t = 0.2 in
Parameter Definitions
Θ = Element Distortion Angle
At X = 0 UX = UY = UZ = 0
At X = L Unit Load
a. Extension (FX = 1)
b. In-plane (FY = 1)
c. Out-plane (FZ = 1)
d. Twist (Equivalent FX, FY forces applied)

Representative Mesh Options

Results Comparison

Tip LoadingTip Displacement Ratio / % Error in Energy Norm
Shape/DirectionPLANE183 Triangular elementPLANE182PLANE183SOLID185SOLID187SOLID186
Rectangular
Extension.9976.9960.9995.9880.99311.9943
In-Plane Shear.98322.993240.9870.97825.96030.97117
Out-Of-Plane Shear.97327.95933.96122
Twist.89210.91024.9038
Trapezoidal (θ = 15°)
Extension.9976.9974.9995.9915.99311.9943
In-Plane Shear.98223.29366.9867.27266.95931.96918
Out-Of-Plane Shear.21567.95833.96023
Twist.85416.91024.9038
Trapezoidal (θ = 30°)
Extension.99760.999414.9936.99311.9943
In-Plane Shear.976260.10964.98217 .10064.95432.95723
Out-Of-Plane Shear.07264.95434.95425
Twist.74227.91025.9038
Trapezoidal (θ = 45°)
Extension.99760.99931.0004.9944.99311.9943
In-Plane Shear.961320.052590.96736.04759.93937.89134
Out-Of-Plane Shear.03060.94138.92136
Twist.56340.91025.9038
Parallelogram (θ = 15°)
Extension.9986.9974.9995.9915.99311.9943
In-Plane Shear.983220.81238.9881.79840.95931.97118
Out-Of-Plane Shear.74949.95833.96023
Twist.88619.91024.9038
Parallelogram (θ = 30°)
Extension.9986.9994.9995.9935.99311.9943
In-Plane Shear.98024.68051.9914.66953.95433.97221
Out-Of-Plane Shear.60865.95335.95526
Twist.86634.91025.9038
Parallelogram (θ = 45°)
Extension.99860.99931.0004.9944.99211.9943
In-Plane Shear.970290.632550.9978.62457.94038.96827
Out-Of-Plane Shear.52874.93540.94232
Twist.82051.90926.9038

Assumptions, Modeling Notes, and Solution Comments

  1. The straight cantilever beam is a frequently used test problem applicable to beam, plate, and solid elements. The problem tests elements under constant and linearly varying strain conditions. Although the problem appears rather simplistic in nature, it is a severe test for linear elements, especially when distorted element geometries are present.

  2. The fixed boundary conditions at the left edge of the beam are not representative of a "patch test." Thus, under extensional loading, the finite element solution will not agree with the beam theory solution.

  3. Element solution accuracy degrades as elements are distorted. The degradation is more pronounced for linear elements (PLANE182, SOLID185) than it is for quadratic elements (PLANE183, SOLID187, SOLID186). The degradation in performance for linear elements is most pronounced for bending loads coupled with irregular element shapes.

  4. Distorted linear elements show more pronounced locking (excessive stiffness) for trapezoidal element shapes than for parallelogram shapes. Results clearly show under predicted displacements for trapezoidal element shapes for PLANE182 and SOLID185.

  5. The quadratic elements show very good performance under all loadings and geometries. The good performance is attributed to the ability of the elements to properly handle bending and shear energy, in contrast to the linear elements.

  6. The percent error in energy norm for each test case is displayed against the tip displacement ratio to illustrate the general correlation between the two for this particular problem. The linear elements show a patterned correlation between solution accuracy and the percent error in energy norm. For the quadratic elements the correlation pattern is similar. In both cases, the results illustrate that a considerable bandwidth on norm values may exist at any desired solution accuracy level for problems under a variety of load conditions with irregular element shapes.