ISOPARAMETRIC ELEMENT Bruce Irons 1924 1983 Quadrilateral Shape

ISOPARAMETRIC ELEMENT Bruce Irons, 1924 -1983 • Quadrilateral Shape – – – Most commonly used element (irregular shape) Generalization of rectangular element Use mapping to transform into a square (Reference element). The relationship between (x, y) and (s, t) must be obtained. All formulations are done in the reference element. y v 3 v 4 4 3 t u 3 u 4 4 (-1, 1) v 1 1 3 (1, 1) s u 1 v 2 2 Actual Element (x, y) u 2 1 (-1, -1) 2 (1, -1) x Reference Element (s, t) 1

ISOPARAMETRIC MAPPING • Definition – the same interpolation method is used for displacement and geometry. • Procedure – Construct the shape functions N 1, N 2, N 3, and N 4 at the reference element t 4 (-1, 1) 3 (1, 1) s 1 (-1, -1) 2 (1, -1) 2

HOW DOES IT LOOK? y 4 t 3 4 (-1, 1) 3 (1, 1) s 1 1 (-1, -1) 2 2 (1, -1) x Two questions: For a given (x, y), find corresponding (s, t). For a given (s, t), find corresponding (x, y). Linearity of transformation on boundary is important for continuity 3

THE ISO PART – Use the shape functions for interpolating displacement and geometry. – For a given value of (s, t) in the parent element, the corresponding point (x, y) in the actual element and displacement at that point can be obtained using the mapping relationship. Displacement interpolation Geometry interpolation 4

EXAMPLE 6. 6 • Find mapping point of A in the physical element – At point A, (s, t) = (0. 5, 0. 5) – Physical coord y 3 (2, 4) 4 (− 1, 1) 2 (4, 4) t 3 (1, 1) A (. 5, . 5) s B (1, 2) 4 (0, 0) 1 (6, 0) x 1 (− 1, − 1) 2 (1, − 1) 5

EXAMPLE cont. • Find mapping point of B in the reference element – At point B, (x, y) = (1, 2) y 3 (2, 4) 4 (− 1, 1) 2 (4, 4) t – Thus, (s, t) = (0, 1) 3 (1, 1) A (. 5, . 5) s B (1, 2) 4 (0, 0) 1 (6, 0) x 1 (− 1, − 1) 2 (1, − 1) 6

Quiz-like questions • Consider the element in Example 6. 4. We are given the u displacement at the four nodes as u 1=1, u 2=3, u 3=0, u 4=2. Find the equation of u along the 3 -4 edge of the quadrilateral as function of x or y. Why is it important that it is linear? • What point in the real element corresponds to the center of the reference element (s=t=0)? • What is the equation of the line y=x in the s, t plane? • Answers in the notes page. y 3 (2, 4) 4 (− 1, 1) 2 (4, 4) t 3 (1, 1) A (. 5, . 5) s B (1, 2) 4 (0, 0) 1 (6, 0) x 1 (− 1, − 1) 2 (1, − 1) 7

JACOBIAN OF MAPPING Shape functions are given in (s, t). But, we want to differentiate w. r. t. (x, y) in order to calculate strain and stress. Use chain rule of differentiation. In Matrix Form Matrix [J] is called the Jacobian matrix of mapping. How to calculate matrix [J]? 8

JACOBIAN OF MAPPING cont. • Derivatives of shape functions w. r. t. (x, y) coordinates: • Determinant |J|: Jacobian • What happen if |J| = 0 or |J| < 0? – shape function derivative cannot be obtained if the |J| = 0 anywhere in the element – Mapping relation between (x, y) and (s, t) is not valid if |J| = 0 or |J| < 0 anywhere in the element (– 1 ≤ s, t ≤ 1). 9

MEANING OF NON-POSITIVE JACOBIAN • Every point in the reference element should be mapped into the interior of the physical element • When an interior point in (s, t) coord. is mapped into an exterior point in the (x, y) coord. , Jacobian becomes negative • If multiple points in (s, t) coordinates are mapped into a single point in (x, y) coordinates, the Jacobian becomes zero at that point • It is important to maintain the element shape so that the Jacobian is positive everywhere in the element • Jacobian IS negative when an angle exceeds 180 o y x 10

EXAMPLE 6. 7 Jacobian must not be zero anywhere in the domain (-1 ≤ s, t ≤ 1) y • Nodal Coordinates 3(2, 2) 4(0, 1) • Iso-Parametric Mapping 1(0, 0) 2(1, 0) x • Jacobian Matrix 11

JACOBIAN • Jacobian – It is clear that |J| > 0 for – 1 ≤ s ≤ 1 and – 1 ≤ t ≤ 1. Constant s t 4 (-1, 1) 3 (1, 1) s Constant t 1 (-1, -1) 2 (1, -1) 12

Quiz-like problem • Find the equation of the determinant of the Jacobian if Node 3 of the quadrilateral element of Slide 11 is moved the midpoint of Node 2 and Node 4, that is to (0. 5, 0. 5), so that the quadrilateral is actually a triangle. 13