Reference no: EM132313323
Finite Element Methods in Engineering Assignment - FE modelling of temperature distribution in a steel plate
Project Description: This project aims to solve a 2D steady state heat transfer in a plate which subjects to different boundary conditions by manually calculation and ANSYS software.
Problem: Consider steady two-dimensional heat transfer in a steel thin plate whose cross section is given in the following figure. There is no heat generation; the temperature difference through the plate thickness is negligible. The thermal conductivity of the steel plate is 50 W/(m·oC). The entire top surface 'AD' is subjected to constant heat flux q'' = 2000W/ m2. The bottom surface 'BC' is maintained at 150oC. The left side surface 'AB' of the body is insulated; the right side surface 'CD' of the body is under convection with ambient air at Tf = 20oC with a convection coefficient of h = 75 W/(m2·oC). The mesh size is Δx = 0.15m, Δy = 0.2m for manual calculation.

(1) Use the finite element method to determine the temperature for each node. You are requested to manually construct the global stiffness matrix and global load matrix.
In the report, you need to show the procedure to determine the stiffness and load matrix of each element, and how to construct of global stiffness and load matrix.
(Tips: As solving node temperature by global stiffness and load matrix is too difficult by hand, Matlab is recommended. The other program or excel also can be used for this purpose. The matlab code or excel file needs to be attached.)
(2) You are also requested to manually calculate the temperature gradient at center of each element.
In the report, you need to show the procedure how to manually calculate the temperature gradient at center of each element.
(3) You are requested construct a 2D model in ANSYS, and use ANSYS to solve the question in Item 1 and compare it with your manual calculation.
In the report, you need to show the geometry of the body, mesh of the model, boundary conditions, temperature and/or heat flux results. The comparison of manual calculation and ANSYS solution is requested.
(Tips: The mesh size in ANSYS can be different with manual calculation. But you need to determine the temperature for the locations in Item 1 for comparison).
(4) Use ANSYS to determine average reaction heat transfer rate acting on the bottom surface 'BC' and right surface 'CD'. You also need to manually determine the heat transfer rate at the top surface 'AD' and left surface 'AB'. Then you are requested to discuss the energy balance of the system.
(Tips: In ANSYS, you can use reaction probe shows the Heat transfer rate).
(5) Plot the temperatures along the line A and E for three convection coefficients acting on the right side 'CD', 10 W/(m2·oC), 75 W/(m2·oC) and 150 W/(m2·oC), and discuss the results.
(Tips: You can solve the results for h = 10 W/(m2·oC) and h = 150 W/(m2·oC) either manually or using ANSYS).