DarylL - loganAFirstCourse Bar Convergence-1
DarylL - loganAFirstCourse Bar Convergence-1
DarylL - loganAFirstCourse Bar Convergence-1
and qpp AL E ^1x þ 2d^2x Þ ÿ f^ ¼ 0
¼ ðÿ2 d 2x
qd^2x 2 L2
( ) ( )
d^1x f^1x
qpp AE 1 ÿ1 0
¼ ÿ ¼ ð3:10:27Þ
^
qfdg L ÿ1 1 d^2x f^
2x
0
As expected, Eq. (3.10.28) is identical to the stiffness matrix obtained in Section 3.1.
Finally, instead of the cumbersome process of explicitly evaluating pp , we can
use the matrix differentiation as given by Eq. (2.6.12) and apply it directly to Eq.
(3.10.19) to obtain
qpp
¼ AL½BT ½D½Bfdg
^ ÿ f f^g ¼ 0 ð3:10:29Þ
^
qfdg
where ½DT ¼ ½D has been used in writing Eq. (3.10.29). The result of the evaluation
of AL½BT ½D½B is then equal to ½^k given by Eq. (3.10.28). Throughout this text, we
will use this matrix differentiation concept (also see Appendix A), which greatly sim-
plifies the task of evaluating ½^k.
To illustrate the use of Eq. (3.10.20a) to evaluate the equivalent nodal loads for a
bar subjected to axial loading traction T^x , we now solve Example 3.12.
Example 3.12
A bar of length L is subjected to a linearly distributed axial loading that varies from
zero at node 1 to a maximum at node 2 (Figure 3–28). Determine the energy equiva-
lent nodal loads.
Using Eq. (3.10.20a) and shape functions from Eq. (3.10.9), we solve for the
energy equivalent nodal forces of the distributed loading as follows:
8 9
x^ >
>1 ÿ >
( ) ðð >
>
f^1x
ð L<
L=
>
T
f f^0 g ¼ ¼ ½N fT^x g dS ¼ fC x^g d x^ ð3:10:30Þ
f^2x S1
0 >
>
> x^ > >
>
: ;
L
8 2 9L
C x^ C x^3 >
ÿ
>
> >
> >
< 2 3L =
¼
>
>
> C x^3 >
>
>
: ;
3L 0
8 9
2
>
> CL > >
< 6 >
> =
¼ ð3:10:31Þ
>
>
> CL 2 > >
>
: ;
3
where the integration was carried out over the length of the bar, because T^x is in units
of force/length.
Note that the total load is the area under the load distribution given by
1 CL 2
F ¼ ðLÞðCLÞ ¼ ð3:10:32Þ
2 2
Therefore, comparing Eq. (3.10.31) with (3.10.32), we find that the equivalent nodal
loads for a linearly varying load are
1
f^1x ¼ F ¼ one-third of the total load
3
ð3:10:33Þ
^ 2
f2x ¼ F ¼ two-thirds of the total load
3
In summary, for the simple two-noded bar element subjected to a linearly varying
load (triangular loading), place one-third of the total load at the node where the dis-
tributed loading begins (zero end of the load) and two-thirds of the total load at the
node where the peak value of the distributed load ends. 9
Example 3.13
For the rod loaded axially as shown in Figure 3–29, determine the axial displacement
and axial stress. Let E ¼ 30 10 6 psi, A ¼ 2 in. 2 , and L ¼ 60 in. Use (a) one and (b)
two elements in the finite element solutions. (In Section 3.11 one-, two-, four-, and
eight-element solutions will be presented from the computer program Algor [9].
3.10 Potential Energy Approach to Derive Bar Element Equations d 117
where Tx is a line load in units of pounds per inch and f^0 ¼ F 0 as x ¼ x^. Therefore,
using Eq. (3.1.4) for ½N in Eq. (3.10.34), we obtain
8 9
>
> x
ð L< 1 ÿ > >
L
=
fF0 g ¼ fÿ10xg dx ð3:10:35Þ
0 >
> x > >
: ;
L
ÿ10ð60Þ 2 >
8 9 8 9 8 9
> ÿ10L 2 10L 2 > ÿ10L 2 >
þ
> >
> > >
> >
> >
> >
> >
>
6
F1x < 2 3 = < 6 = < =
or ¼ ¼ ¼
F2x >
> ÿ10L 2 >
> > ÿ10L 2 >
> > > ÿ10ð60Þ 2 >
> >
>
: >
; > : >
; > : >
;
3 3 3
or F1x ¼ ÿ6000 lb F2x ¼ ÿ12;000 lb ð3:10:36Þ
Using Eq. (3.10.33), we could have determined the same forces at nodes 1 and 2—that
is, one-third of the total load is at node 1 and two-thirds of the total load is at node 2.
118 d 3 Development of Truss Equations
We first obtain the element forces. For element 2, we divide the load into a uni-
form part and a triangular part. For the uniform part, half the total uniform load is
placed at each node associated with the element. Therefore, the total uniform part is
ð30 in:Þðÿ300 lb=in:Þ ¼ ÿ9000 lb
and using Eq. (3.10.33) for the triangular part of the load, we have, for element 2,
( ð2Þ ) ( )
ÿ½12 ð9000Þ þ 13 ð4500Þ
f2x ÿ6000 lb
¼ ¼ ð3:10:40Þ
f
ð2Þ
3x
ÿ½12 ð9000Þ þ 23 ð4500Þ ÿ7500 lb
3.10 Potential Energy Approach to Derive Bar Element Equations d 119
For element 1, the total force is from the triangle-shaped distributed load only and is
given by
1
2 ð30 in:Þðÿ300 lb=in:Þ ¼ ÿ4500 lb
On the basis of Eq. (3.10.33), this load is separated into nodal forces as shown:
( ð1Þ ) ( )
1
3 ðÿ4500Þ
f1x ÿ1500 lb
ð1Þ
¼ 2
¼ ð3:10:41Þ
f 3 ðÿ4500Þ
ÿ3000 lb
2x
where the boundary condition d3x ¼ 0 has been substituted into Eq. (3.10.45). Now,
solving equations 1 and 2 of Eq. (3.10.45), we obtain
d1x ¼ ÿ0:006 in:
ð3:10:46Þ
d2x ¼ ÿ0:00525 in:
The element stresses are as follows:
Element 1
d1x ¼ ÿ0:006
1 1
sx ¼ E ÿ
30 30 d2x ¼ ÿ0:00525
¼ 750 psi ðTÞ ð3:10:47Þ
120 d 3 Development of Truss Equations
Element 2
( )
d2x ¼ ÿ0:00525
1 1
sx ¼ E ÿ
30 30 d3x ¼ 0
¼ 5250 psi ðTÞ ð3:10:48Þ
9
or 5L 3
C1 ¼ ÿ ð3:11:4Þ
3AE
Substituting Eq. (3.11.4) into Eq. (3.11.3) makes the final expression for displacement
5
u¼ ðx 3 ÿ L 3 Þ ð3:11:5Þ
3AE
3.11 Comparison of Finite Element Solution d 121
Figure 3–32 Comparison of exact and finite element solutions for axial displacement
(along length of bar)
Substituting A ¼ 2 in.2 , E ¼ 30 106 psi, and L ¼ 60 in. into Eq. (3.11.5), we obtain
u ¼ 2:778 10ÿ8 x3 ÿ 0:006 ð3:11:6Þ
The exact solution for axial stress is obtained by solving the equation
PðxÞ 5x 2
sðxÞ ¼ ¼ ¼ 2:5x 2 psi ð3:11:7Þ
A 2 in 2
Figure 3–32 shows a plot of Eq. (3.11.6) along with the finite element solutions
(part of which were obtained in Example 3.13). Some conclusions from these results
follow.
1. The finite element solutions match the exact solution at the node
points. The reason why these nodal values are correct is that the
element nodal forces were calculated on the basis of being energy-
equivalent to the distributed load based on the assumed linear
displacement field within each element. (For uniform cross-sectional
bars and beams, the nodal degrees of freedom are exact. In general,
computed nodal degrees of freedom are not exact.)
2. Although the node values for displacement match the exact solution,
the values at locations between the nodes are poor using few elements
(see one- and two-element solutions) because we used a linear
displacement function within each element, whereas the exact solution,
Eq. (3.11.6), is a cubic function. However, because we use increasing
122 d 3 Development of Truss Equations
Figure 3–33 Comparison of exact and finite element solutions for axial stress (along
length of bar)
X
Fx ¼ 0: ÿ 12 ð300 lb=in:Þð30 in:Þ þ 4500 lb ¼ 0
For other kinds of elements (other than beams), this adjustment is ignored in practice.
The adjustment is less important for plane and solid elements than for beams. Also,
these adjustments are more difficult to formulate for an element of general shape.