Structural Engineering & Statics
where the load goes
Concept

Shape function — where it appears

The assumed variation of displacement within a finite element, from which its stiffness and its load vector are both derived. It defines the space the solution is found in, so an element is exact for any behaviour its shape functions contain and approximate for everything else.

Named by 3 essays across 2 fields — each of them below, with the objects they name alongside it.

One member's stiffness, scattered into the freedoms it touches. A member's own six-by-six stiffness matrix relates the forces at its two ends to the displacements there, and it is written in the member's own axes. Assembly is two operations and no physics: rotate it into the structure's axes, then add each of its thirty-six entries into the row and column of the global freedom that entry belongs to. Every member does the same, and the sum is the structure. The shaded rows and columns are the six freedoms this one member reaches; every other entry it contributes is exactly zero, and that is the whole reason a global stiffness matrix is sparse. Nothing here is an approximation — the result is the same equilibrium and the same compatibility a hand method writes, in an order a machine can follow.

The answer that depends on how it was divided

Every computed answer in this collection came out of a structure chopped into pieces — elements, strips, stations, trial positions. The chopping is invisible in the result and it is not neutral: some divisions give the exact answer, some give one that is always too stiff, and one of them changes nothing but the cost of getting there.

deflection · Discretisation Two load sets with the same resultant and different work. The two ways of putting a uniform load of 10 kN/m onto a beam element 6.00 m long. Both put 30.0 kN at each node, so both have the same resultant and the same moment about any point — they are equivalent for a rigid body. The consistent set adds a couple of 30.0 kN·m at each end, in opposite senses, which is what makes it do the same virtual work over the element's shape functions as the real load does. The couples cancel in the resultant, which is exactly why the resultant cannot see them, and they are the whole difference between an exact answer and one that is a third out.

Equivalent in work, not in resultant

Two force systems with the same resultant and the same moment about every point are interchangeable — for a rigid body. A finite element is not a rigid body, and substituting one for the other on a beam element leaves the tip of a cantilever a third too low with no warning of any kind.

equilibrium · Force couple A downward pressure that pulls the corners up. The nodal forces that do the same work as a uniform pressure on one eight-node serendipity quadrilateral, as fractions of the whole load, each node's share being the integral of its shape function over the element. The shares are −1/12, −1/12, −1/12, −1/12, 1/3, 1/3, 1/3, 1/3 at its eight nodes, adding up to one. The four corner forces are negative: the pressure pushes down and the work-equivalent forces at the corners pull up.

The corners the pressure pulls up

Replace a uniform pressure on a beam element by its work-equivalent nodal loads and the ends acquire couples the resultant cannot see. Do the same on an eight-node plate element and the corners acquire forces pointing the wrong way: a pressure pushing down is represented by the corners being pulled up, a twelfth of the load each. It is correct, it is the vector that makes the solution the best one available, and it is the reason a printout of nodal forces is not a picture of where a load goes.

equilibrium · Force couple

Named alongside it

The objects these essays reach for when they reach for this one.

Virtual workConvergenceDiscretisationEquivalent loadFinite-elementResultantBandwidthCompatibilityDeflectionDegrees of freedomDistributed loadEigenvalue

All concepts