Problem 5 002
Problem 5 002
Problem 5 002
EXAMPLE 5-002
SOLID – STRAIGHT BEAM WITH STATIC LOADS
PROBLEM DESCRIPTION
In this example, a straight cantilever beam, modeled with solid objects, is
subjected to unit forces at the tip in the three orthogonal directions and unit
moments at the tip about the three orthogonal directions, each in a different load
case. The tip displacements in the direction of the load are compared with hand-
calculated results.
The basic geometry, properties and loading are as described in MacNeal and
Harder 1985. The cantilever beam is 6 inches long, 0.2 inch wide parallel to the
Z direction and 0.1 inch wide parallel to the Y direction. Three different models
are created, each with a different element shape. Models A, B and C use
rectangular-, trapezoidal- and parallelogram-shaped elements, respectively.
It is important to note that this example is an extreme case presented for testing
and verification of the solid object. Solid objects are not in general intended for
use in modeling a beam with an nx1x1 mesh.
Six load cases are created for each model. Load cases 1 through 3 apply unit
forces at the cantilever tip in the X, Y and Z directions, respectively. Load cases
4 through 6 apply unit moments at the cantilever tip about the X, Y and Z
directions, respectively. The moments are applied as couples.
The independent solution is derived using elementary beam theory that assumes
no local Poisson’s effect occurs at the support. The beam is modeled in SAP2000
with all joints fully restrained at the fixed end and with all joints at the free end
assigned to a body constraint. Thus, the SAP2000 model is slightly different
from the hand calculations because in the SAP2000 model there is a slight local
Poisson’s effect at the beam ends.
EXAMPLE 5-002 - 1
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
GEOMETRY
1 - Joint number
Model A – Rectangular Shaped Elements
1 - Solid object number
0.2"
8, 22 9, 23 10, 24 11, 25 12, 26 13, 27 14, 28
Z 1 2 3 4 5 6
Y 1, 15 2, 16 3, 17 4, 18 5, 19 6, 20 7, 21
6 @ 1" = 6"
X
1,15 - Joint numbers
Model B – Trapezoidal Shaped Elements 1 - Area object number
0.9" 1.2" 0.8" 1.2" 0.8" 1.1"
0.2"
8, 22 9, 23 10, 24 11, 25 12, 26 13, 27 14, 28
1 2 3 4 5 6
1, 15 2, 16 3, 17 4, 18 5, 19 6, 20 7, 21
1.1" 0.8" 1.2" 0.8" 1.2" 0.9"
6"
0.2"
8, 22 9, 23 10, 24 11, 25 12, 26 13, 27 14, 28
1 2 3 4 5 6
1, 15 2, 16 3, 17 4, 18 5, 19 6, 20 7, 21
0.9" 1" 1" 1" 1" 1.1"
6"
PROPERTIES
E = 10,000,000 lb/in2
ν = 0.3
G = 3,846,154 lb/in2
EXAMPLE 5-002 - 2
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
LOADING
The following table defines the loading applied to each model.
Px
1 Axial Fx = +0.25 lb at joints 7, 14, 21 and 28
Extension
Vz and My
2 Shear and Fz = +0.25 lb at joints 7, 14, 21 and 28
bending
Vy and Mz
3 Shear and Fy = +0.25 lb at joints 7, 14, 21 and 28
bending
EXAMPLE 5-002 - 3
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
RESULTS COMPARISON
The SAP2000 results are presented separately for models with and without the
incompatible bending modes option. The independent results are hand calculated
using the unit load method described on page 244 in Cook and Young 1985. In
addition, the torsional stiffness of the section, J, is calculated using item 4 in
Table 20 on page 290 in Roark and Young 1975. Independent results are also
published in MacNeal and Harder 1985.
EXAMPLE 5-002 - 4
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 5
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 6
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 7
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 8
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 9
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
For all shapes of elements, the bending behavior of the models is improved by
meshing along the length of the beam.
For all shapes of elements, the results for twist are improved by using a 4 x 8
mesh in the cross-sectional plane of the beam.
In general the models that do not use the incompatible bending modes option
seem to be more sensitive to element aspect ratio than models that use the
incompatible bending modes option. In addition, the models that do not use the
EXAMPLE 5-002 - 10
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
incompatible bending modes option appear to require a much more refined mesh
to achieve acceptable results.
CONCLUSIONS
The model with a 6x1x1 mesh, rectangle-shaped elements, and the incompatible
bending modes option has acceptable results for all loading types except twist.
All other models with a 6x1x1 mesh have unacceptable results for all loadings
except axial extension. The results for all models are improved by refining the
mesh.
Bending results are improved by meshing along the length of the beam.
Acceptable twisting results are achieved with at least four elements in each
direction of the cross-sectional plane of the beam and with the aspect ratio of the
elements in the cross-sectional plane of the beam approaching 1:1.
EXAMPLE 5-002 - 11
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
HAND CALCULATION
EXAMPLE 5-002 - 12
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 13
Software Verification
PROGRAM NAME: SAP2000
REVISION NO.: 2
EXAMPLE 5-002 - 14