Topic 3: Finite element formulation and analysis for 2D Thermal Problem 3m 14 3 56 15 2 7 14 13 T=200C T=25C 2 3 8 13 2m 6 12 12 4 10 1 11 Consider the following finite element mesh with the boundary conditions indicated (2 sides adiabatic, 2 sides with uniform temperatures specified). Write a finite element code in ANY language that you wish (MATLAB, fortran, C, etc.) to find the temperature distribution at all the nodes. Note the following: (a) The specified input to the code should be the nodal temperatures at nodes 1-3 and 11-14, along with the heat flux q = 0 at nodes 4,7,8,10 (nodes 1,3,11,14 are corner nodes that still have to maintain the specified temperatures so the heat flux cannot be 0). (b) You can directly use the expressions for the local stiffness matrix and local load vector of a generic element, but you need to write a code to populate each local stiffness matrix and local load vector, and in turn to populate the global stiffness matrix and global load vector. (c) Assume Q = 0, and thermal conductivities kxx = kyy = 45 W/mK. (d) For ways to start the coding, I recommend looking in the notes on the chapter on 1D FEM where example codes are provided. (e) The output should be the nodal temperatures.