Skip to main content

Iso/super/sub-parametric; Hermite shape function

Short answer

Isoparametric: same shape functions for geometry and displacement. Superparametric: geometry interpolation higher-order than displacement. Subparametric: geometry lower-order. Hermite beam shape functions (length L): N1=1-3x^2/L^2+2x^3/L^3, N2=x-2x^2/L+x^3/L^2, N3=3x^2/L^2-2x^3/L^3, N4=-x^2/L+x^3/L^2.

Define iso-parametric, super parametric and sub parametric elements. Derive the expression for the Hermite shape function used for the interpolation in beam elements.

Types of parametric elements

  • Isoparametric: the same shape functions (same number of nodes) are used for the geometry and for the displacement, x=∑Nixix=\sum N_ix_i, u=∑Niuiu=\sum N_iu_i.
  • Superparametric: the geometry is described by a higher-order interpolation (more nodes) than the displacement, e.g. quadratic geometry and linear displacement.
  • Subparametric: the geometry is described by a lower-order interpolation than the displacement, e.g. straight sides with a quadratic displacement.

Hermite shape functions for a beam element

A beam element of length LL has two nodes with two DOFs each: deflection vv and slope θ=dv/dx\theta=dv/dx. Both vv and θ\theta must be continuous between elements (C1C^1 continuity), so a cubic (4 constants) is assumed:

v(x)=α1+α2x+α3x2+α4x3v(x)=\alpha_1+\alpha_2x+\alpha_3x^2+\alpha_4x^3 θ(x)=dvdx=α2+2α3x+3α4x2\theta(x)=\frac{dv}{dx}=\alpha_2+2\alpha_3x+3\alpha_4x^2

Apply the four nodal conditions: v(0)=v1v(0)=v_1, θ(0)=θ1\theta(0)=\theta_1, v(L)=v2v(L)=v_2, θ(L)=θ2\theta(L)=\theta_2:

{v1θ1v2θ2}=[100001001LL2L3012L3L2]{α1α2α3α4}\begin{Bmatrix}v_1\\\theta_1\\v_2\\\theta_2\end{Bmatrix}=\begin{bmatrix}1&0&0&0\\0&1&0&0\\1&L&L^2&L^3\\0&1&2L&3L^2\end{bmatrix}\begin{Bmatrix}\alpha_1\\\alpha_2\\\alpha_3\\\alpha_4\end{Bmatrix}

Solving: α1=v1\alpha_1=v_1, α2=θ1\alpha_2=\theta_1,

α3=3(v2−v1)L2−2θ1+θ2L,α4=2(v1−v2)L3+θ1+θ2L2\alpha_3=\frac{3(v_2-v_1)}{L^2}-\frac{2\theta_1+\theta_2}{L},\qquad \alpha_4=\frac{2(v_1-v_2)}{L^3}+\frac{\theta_1+\theta_2}{L^2}

Substituting and collecting the coefficients of v1,θ1,v2,θ2v_1,\theta_1,v_2,\theta_2 gives v=[N1 N2 N3 N4]{d}v=[N_1\ N_2\ N_3\ N_4]\{d\} with

N1=1−3x2L2+2x3L3N2=x−2x2L+x3L2N3=3x2L2−2x3L3N4=−x2L+x3L2\begin{aligned} N_1&=1-\frac{3x^2}{L^2}+\frac{2x^3}{L^3}\\ N_2&=x-\frac{2x^2}{L}+\frac{x^3}{L^2}\\ N_3&=\frac{3x^2}{L^2}-\frac{2x^3}{L^3}\\ N_4&=-\frac{x^2}{L}+\frac{x^3}{L^2} \end{aligned}

In terms of s=x/Ls=x/L: N1=1−3s2+2s3N_1=1-3s^2+2s^3, N2=L(s−2s2+s3)N_2=L(s-2s^2+s^3), N3=3s2−2s3N_3=3s^2-2s^3, N4=L(−s2+s3)N_4=L(-s^2+s^3).

Check: N1(0)=1N_1(0)=1, N1(L)=0N_1(L)=0, N1′(0)=N1′(L)=0N_1'(0)=N_1'(L)=0; N2′(0)=1N_2'(0)=1, N2(0)=N2(L)=N2′(L)=0N_2(0)=N_2(L)=N_2'(L)=0, and so on for N3,N4N_3,N_4. Also N1+N3=1N_1+N_3=1 (rigid-body translation is represented).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Computational Techniques in Civil Engineering chapter 4 questions →