1D Spring Systems Finite Element Analysis
1D Spring Systems Finite Element Analysis
1D-Spring Systems
Left Wall
Right Wall
Q0
Rod Area = A
Q0
k0= AE/L
k0= AE/L
Idealizations:
- rod properties, A, E and L, and load, Q0, are known exactly.
- rod has uniform cross-section and material properties
- steel stress-strain relationship is linear (small strains)
- load is applied at rod centroid, and rod is perfectly straight (no bending)
- rod does not buckle (can't be too slender)
CEE 379
1D-Spring Systems
D1
Q2
Q1
Q3
2
1
2
Unknown Displacements (Du): D1
3
Known Displacements (Dk): D2 = 0., D3 = 0.
Q1 = Q0
(applied loads)
3.
D3
Q2 , Q3
(reactions)
dF
qN
qF
k0
N
The forces at the ends of the spring (qN ,qF) are related to the displacements at the ends of
the springs (dN ,dF) as follows:
qN1 = k0 ( dN1 - dF1)
This relationship is the same for spring #1 and spring #2, because the spring stiffness is
the same.
CEE 379
1D-Spring Systems
4. Express Nodal Equilibrium in Terms of Kinematic Degrees of Freedom (D1, D2, D3)
qN1
qF1
qN2
qF2
Q2
2
Q1 = Q0
qN1
qF1
Q3
qN2
qF2
Nodal Equilibrium:
Node 1:
qF1
Node 2:
qN1
+ qN2
= Q1
(= Q0 )
= Q2
Node 3:
qF2
= Q3
= Q1= Q0
Node 2:
( k0 dN1 - k0 dF1)
= Q2
(-k0 dN2 + k0 dF2) = Q3
Node 3:
Enforce Compatibility:
Element Connectivity:
dN1 = D2 ;
dF1 = D1 ;
= Q1= Q0
Node 2:
( k0 D2 - k0 D1)
= Q2
(-k0 D1 + k0 D3) = Q3
D2 = 0.;
Boundary Conditions:
D3 = 0.
Node 1:
k0 D1) +
Node 2:
-k0 D1)
(k0 D1
= Q1= Q0
= Q2
(-k0 D1
Node 3:
= Q3
6.
dF2 = D3 ;
Node 1:
Node 3:
5.
dN2 = D1 ;
D1 = Q0 / 2k0
Q2 = -Q0 / 2
Node 3:
Q3 = - Q0 / 2
CEE 379
7.
1D-Spring Systems
From Step 3:
Spring #1
qN1 = - Q0 / 2
qF1 = Q0 / 2
Spring #2
qN2 = Q0 / 2
qF2 = - Q0 / 2
CEE 379
1D-Spring Systems
Same as before
2.
Same as before.
3.
dF
qN
qF
k0
N
q N1
=
q F1
k1
k
1
k1
k1
qN 2
=
qF 2
k2
k
2
k2
k 2
d N1
d F1
q1 = k1 d1
dN2
dF2
q2 = k2 d2
4. Express Nodal Equilibrium in Terms of Kinematic Degrees of Freedom (D1, D2, D3)
qN1
qF1
qN2
qF2
Q2
2
Q1 = Q0
qN1
qF1
Q3
qN2
qF2
qF1
Node 2:
qN1
+ qN2
= Q1
= Q2
(= Q0 )
CEE 379
1D-Spring Systems
Node 3:
qF2
= Q3
Node1
Node2
Node3
Node1
Node2
Node3
5.
k1 + k 2
k
1
k 2
k1
k1 + k 2
k
1
k 2
k1
k1
0
k1
0
k 2 D1 Q1
0 D2 = Q2 or KD = Q
k 2 D3 Q3
k 2 D1 Q0
0 0. = Q2
k 2 0. Q3
We only need the first equation to solve for the unknown displacement, D1. Better yet,
we only to consider the first stiffness term of the first equation (k1+k2), because the other
terms are multiplied by displacements that we already know to be equal to 0.0. In
practice we benefit greatly from noting that it is not necessary to determine all of the
stiffness coefficients of K.
From the first equation, D1 = Q0/(k1+k2)
6.
From Step 3:
q
Spring #1, q1 = k1 d1 = N 1 =
q F1
q
Spring #2, q2 = k2 d2 = N 2 =
qF 2
k1
k
1
k1Q0
0
.
k1 Q k1 + k 2
=
0
(tension)
k1 k + k k1Q0
1
2
k +k
1
2
k Q
Q0 2 0
k2
k +k
k1 + k 2 = 1k Q2 (compress.)
k2
2 0
0. k + k
1
2
k2
k
2
2.
3.
CEE 379
1D-Spring Systems
q2 = k2 d2
Consider nodal equilibrium at all degrees of freedom. Take into account connectivity
compatibility as well as the boundary conditions.
KD = Q
To distinguish among various components of the stiffness matrix, the stiffness equation,
KD = Q, can be partitioned as follows.
K 12 D u Q k
K
= Q
=
KD = 11
K 21 K 22 D k Q u
where Qk and Dk are the known external loads and displacements, and Qu and Du are the
known external loads and displacements.
In our example, the components of the partitioned matrix are:
K 11 = [k1 + k 2 ]
K 12 = [ k1
k2 ]
k
K 21 = 1
k 2
k
K 22 = 1
0
0
k 2
Qk = K11 Du + K12Dk
(Eq. 1)
CEE 379
Qu = K21 Du + K22 Dk
5.
1D-Spring Systems
(Eq. 2)
In practice we benefit greatly from noting that it is not necessary to determine all of the
stiffness coefficients of K. From the first equation, one can solve the unknown
displacements,
Du = [K11]-1 [Qk]
Note that in our example, K11 = k1+k2, so [K11]-1 = 1/(k1+k2)
6.
Use matrix equation #2, if necessary to determine the equations. If you have not
assembled the K21 and K22 matrices, you can also get the reactions by considering the
ends of the members, which is what is usually done in practice.
7.
From Step 3:
Spring #1, q1 = k1 d1
Spring #2, q2 = k2 d2
Repeat for all springs
IMPLEMENTATION NOTES
It is often cumbersome to assemble the full stiffness matrix, K, particularly in problems
with many degrees of freedom, and in which the specified displacements are 0.0. In this
case, the K12 components multiply 0.0. In solving your homework problems, only
assemble the full K matrix if necessary to solve the problem, or required by the problem
statement.
The following pages demonstrate the solution of a three-spring problem with two free
DOF using a spreadsheet program.