Matrix Finite Element
Methods By Energy Minimization
· Write Potential Energy as a
function of nodal displacements: V=V(u)
· Minimize V(u):
· Solve for u
To calculate V(u):
Write the strain and axial force of each member as a function of nodal displacements.
For
Element e connected to joints (a) and (b):
Write the Strain Energy for
each member:
and
are the element
stiffness matrix and element displacement vector.
is symmetric!!!
In 3D
(Same as it ever
was…)
The Global Stiffness Matrix
The total strain
energy of the truss may be computed by adding together the strain energy of
each element:
It is more
convenient to express W in terms of the global displacement vector, u
Consider a
Two-Member (Two-Dimensional) Truss:
W(1) and W(2)
in terms of the global displacement vector u:
[K] is the Global Stiffness Matrix. It is the sum
of all the element stiffness matrices.
Because the element stiffness matrix is symmetric, the global stiffness
matrix must also be symmetric.
CAUTION:
SOME ASSEMBLY REQUIRED!!!
(Batteries not included)
Bookkeeping:
The displacement of node #b
in the ei-direction, ui(b) (i=1,2)
is the (2(b-1)+i)th element of the global displacement vector, u:
The degree of freedom number assigned to the displacement
of node b in the i-direction is
(2(b-1)+i).
In 3D: and the degree of freedom number assigned
to the displacement of node b in the i-direction is
(3(b-1)+i).)
The global stiffness positioning (2D):
Where do the element stiffness components land?
Externally Applied Forces
These are specified by prescribing the node on which
the loading acts and the force vector P.
Potential energy
due this load is then
The total
contribution to the potential energy due to prescribed forces on all of the
loaded nodes is
Global residual force vector: r
But of course, we have to
write Q in terms of the global force vector u.
Q=-r×u
ASSEMBLY REQUIRED!!!
Recall
the degree of freedom number assigned to the displacement
of node b in the i-direction is
(2(b-1)+i).
3D degree of freedom number assigned to
the displacement of node b in the i-direction is
(3(b-1)+i).
Stay Tuned!!
Still
to Come:
·
Minimize
V[u]
·
Take
boundary conditions into account
·
Solve
for u.
·
Go
back to get element forces and strains.