Skip to main content

Linear shape functions for one-dimensional element

Short answer

For a two-node bar element N1=1-x/L and N2=x/L. They satisfy Kronecker delta, sum to unity, and give C0 displacement with constant strain. Strain-displacement matrix is [B]=(1/L)[-1 1].

Derive the shape functions for a two-node one-dimensional (linear) element of length LL. State their properties and use them to write the strain-displacement matrix.

Displacement function

Two nodal values u1u_1 (at x=0x = 0) and u2u_2 (at x=Lx = L) allow a linear polynomial:

u(x)=a1+a2xu(x) = a_1 + a_2x
   u1                u2
   o-----------------o
  x=0               x=L

Apply the nodal conditions:

u(0)=a1=u1,u(L)=a1+a2L=u2⇒a2=u2−u1Lu(0) = a_1 = u_1,\qquad u(L) = a_1 + a_2L = u_2 \Rightarrow a_2 = \frac{u_2 - u_1}{L}

Substitute back:

u=u1+u2−u1Lx=(1−xL)u1+xL u2u = u_1 + \frac{u_2 - u_1}{L}x = \left(1 - \frac{x}{L}\right)u_1 + \frac{x}{L}\,u_2

So u=N1u1+N2u2u = N_1u_1 + N_2u_2 with

N1=1−xL,N2=xLN_1 = 1 - \frac{x}{L},\qquad N_2 = \frac{x}{L}

Properties

  • Ni=1N_i = 1 at node ii and 00 at the other node (Kronecker delta property).
  • N1+N2=1N_1 + N_2 = 1 everywhere (partition of unity); this allows rigid-body motion.
  • Each NiN_i is linear, so the displacement is continuous between elements (C0C^0), but the strain is constant per element.
  • In natural coordinate ξ=2x/L−1\xi = 2x/L - 1: N1=1−ξ2N_1 = \dfrac{1 - \xi}{2} and N2=1+ξ2N_2 = \dfrac{1 + \xi}{2}.

Strain-displacement matrix

ε=dudx=dN1dxu1+dN2dxu2=1L[−1  1]{u1u2}\varepsilon = \frac{du}{dx} = \frac{dN_1}{dx}u_1 + \frac{dN_2}{dx}u_2 = \frac{1}{L}\left[-1\ \ 1\right]\begin{Bmatrix} u_1 \\ u_2 \end{Bmatrix} [B]=1L[−1  1][B] = \frac{1}{L}\left[-1\ \ 1\right]

✓ Answer checked by Claude Opus 5.5 (high effort)

All Finite Element Method chapter 5 questions →