← Back to Interactive Tools

8- and 9-Node Isoparametric Elements — Curved Edges

CE/ME 532 · Section 3.5 · Switch between the Q8 serendipity and Q9 Lagrangian element. Drag the red corner nodes, the orange mid-side nodes, or the purple centre node (Q9) to bend the edges and interior. Drag the green probe to sample the shape functions and |J|.

Parent element   (ξ,η)  space

Actual element   (x,y)  space

ξ = η = x = 0.000 y = 0.000 |J| = 0.000 ΣN = 1.000

The blue mesh on the right is the image of ξ=const and η=const grid lines. With Q8/Q9 the image lines are parabolic, so each physical edge can curve. Red shading marks the region where |J| < 0. Mid-side nodes push off their chord midpoint more than the ¼–rule allows ⇒ the Jacobian can fold.

1. Node layout & shape-function families

Both elements share the same 8 boundary nodes: 4 corners and 4 mid-side nodes, one per edge of the parent square. The Q9 Lagrangian adds a 9th node at the centre (ξ,η) = (0,0).

Q8 serendipity — corners and mid-sides only:
Corners  Ni = ¼(1+ξiξ) (1+ηiη) (ξiξ+ηiη−1)
Edge ξ=0  N5 = ½(1−ξ2)(1−η),  N7 = ½(1−ξ2)(1+η)
Edge η=0  N6 = ½(1+ξ)(1−η2),  N8 = ½(1−ξ)(1−η2)
Q9 Lagrangian — tensor product of 1-D quadratic Lagrange:
L(t) = ½t(t−1),  L0(t) = 1−t2,  L+(t) = ½t(t+1)
Ni(ξ,η) = La(i)(ξ) Lb(i)(η)
Key difference. The Q9 Lagrange space contains the full biquadratic monomials {1, ξ, η, ξη, ξ2, η2, ξ2η, ξη2, ξ2η2}. Q8 is missing the ξ2η2 bubble — it is not complete to degree 2 in ξη space.

2. Live shape-function values at the probe

At the current (ξ,η):

Partition of unity check: Σ Ni should equal 1 at every point. Live value: 1.000.

3. The isoparametric map

Use the same Ni for geometry and displacement:

x = Σi Ni(ξ,η) xiy = Σi Ni(ξ,η) yi

With mid-side nodes each physical edge is the quadratic Lagrange interpolation of its three nodes — a parabola (or a straight line when the mid-side node sits on the chord midpoint).

4. Jacobian at the probe

Same definition as the 4-node case, but the derivatives Ni,ξ, Ni,η are now quadratic:

J = ,  |J| = 0.0000
Unlike Q4, |J|(ξ,η) here is a cubic polynomial in (ξ,η). It can change sign even when all nodes look reasonable.

5. Curved-edge rule: the ¼ check

Along one edge (say η=−1) the three active shape functions are the 1-D quadratic Lagrange functions of ξ. The edge in physical space is therefore a parabola through the three edge nodes.

The 1-D quadratic mapping stays monotone iff the mid-side node is within the middle half of the edge (|\Delta| < L/4 from the chord midpoint). Violate this and the edge folds on itself — the element's Jacobian goes negative near that edge.

Try the Bad mid-side preset to see the red |J| < 0 region appear.

6. Q8 vs Q9 — which to pick?

Q8 has 16 DOFs, Q9 has 18 DOFs. The extra 2 DOFs buy the full biquadratic polynomial space and, with 3×3 Gauss, exact integration on affine meshes.