X=0 x = 1/2 x= L u U2 Uz (a) Trial solution for a 1-D quadratic elastic bar element can be written as follows: ū(x) = [N]{u} where, [N] = [N1 N2 N3] and {u} u2 13 1 and Ni L2 L2 [N] and {u} are known as interpolation function matrix and nodal displacement, respectively. (272 – 3L + L´), N= = (22- La), Ns = 12 (2=– LE) Derive the expression for element stiffness matrix, (Kelem) and element force vector {F} using following expression and exact integration procedure. dN NT dr Here, (B) is known as strain-displacement matrix. (b) Use finite element formulation derived in (a) to generate numerical values of element stiffness matrix and force vector for the following bar with uniformly distributed load. The bar is fixed at r = 0 and a point load, P = 10 kN/m is applied. First, solve for displacement for the problem using 1-element model. 9(x) = 5 kN/m E = 120 Gpa PA = 1×10-3 m² L = 1 m (c) Now repeat the problem using uniformly meshed 2-element model as shown below. Generate element stiffness matrix and force vector for both elements, assemble them and find the nodal displacement. Compare displacement solution from (b) and (c) by plotting against the length (you can use MATLAB or MS Excel to generate the plot). ez u uz us