Heat Transfer Problem
Heat Transfer Problem
Heat Transfer Problem
PROJECT REPORT ON
Submitted to:
Dr. Dinesh Kumar
By:
Bajak Rahul (2014PDE5064)
ACKNOWLEDMENT
I take immense pleasure to express my deep and sincere
gratitude to my guide Prof. Dr. Dinesh Kumar Sir (Department Of
Mechanical Engineering, MNIT Jaipur) for his valuable guidance and
suggestions. He helped me a lot in developing project of this magnitude.
Also while working on this project it helped me a lot in understanding
the subject.
I would also like to thank all my friends for helping me during
the project work and for inevitable suggestions whenever required.
ABSTRECT
CONTENT
1. Introduction :
1.1.
Finite element method
1.2.
Problem statement
2. Development of the FEM model
2.1.
Weak form formulation
2.2.
Finite element model
2.3.
Interpolation function and element matrices
2.4.
Assembly of the element matrices
2.5.
Imposition of the boundary condition
2.6.
Solution
3. Analysis
3.1.
software name and version
3.2.
Analysis using Ansys mechanical APDL
3.3.
Analysis using Ansys workbench 15.0
4. Results and Discussions
4.1.
Comparisons of analytical and FEM modeling results :
5. Conclusion
: Comparisons of results
Chapter 1
Introduction :
This term project is intended to enhancing the knowledge of the
Advanced finite element method which is base of various analysis
software.
Chapter 2
Development of the FEM model
2.1 Weak form formulation:
Above problem is of 2-D steady state heat conduction through plate.
Steady state heat conduction problem is govern by the following
differential equation:
T
(k x ) (k y ) = f(x, y)
x
x
y
y
= temperature (0C)
kx
T
T
nx + k y
n + (T T ) = q n
x
y y
K eij = k (
e
x x
w T
= k (
e
x x
+
+
w T
y y
w T
y y
wf) dxdy w (k x
T
x
nx + k y
T
y
ny ) ds
=B (w, T) l(T)
Where, B(w ,T) and l(T) are the bilinear and linear form respectively,
w T
B(w ,T) = k (
e
x x
w T
y y
wf) dxdy + wT ds
l(t) = wf dx dy + wT ds + w q n ds
e
T = Tje ej (x , y)
j=1
Where ,
e
ei j
e
K ij = k (
x x
e
ei j
y y
) dxdy
Hije = ei ej ds
Pie = ei T ds
he23
he34
e
Hije = 12
ei ej ds + e23 ei ej ds + e34 ei ej ds
0
0
he41
+ e41 ei ej ds
0
he12
he23
he34
e
12
23
34
Pie = 12
T
ei ds + e23 T
ei T ds + e34 T
ei T ds
0
he41
41
+ e41 T
ei T ds
0
x+xi
a
So that ,
x
1e =(1 ) (1 )
x
e2 =( ) (1 )
x
e3 =( ) ( )
x
e4 =(1 ) ( )
) (1
y+yi
b
e
e
12
h12
2
1
[
0
0
1
2
0
0
e34 he34
6
0 0
0 0 0 0
e
e
23 h23 0 2 1 0
0 0
]+
[
]
0 0
0 1 2 0
6
0 0
0 0 0 0
0 0 0 0
2 0 0
e
e
41 h41 0 0 0
0 0 0 0
[
]+
[
0 0 2 1
0 0 0
6
0 0 1 2
1 0 0
1
0
]
0
2
1
0
0
e
e
e
e
23
34
23 h23 T 1
34 h34 T 0
1
[ ]+
[ ]+
[ ]
0
1
1
2
2
2
0
0
1
1
e
e
41
41 h41 T 0
+
[ ]
0
2
1
e
e
12
12
h12
T
A plate shown in the fig. describe the discretization of the plate in the
small element .Every element of the discretization model is having
length of a and height of b.
In our problem the width of the plate is 4a = 0.4 m and height of the
plate 2b = 0.2 m. And because of similarity of the all the elements of the
mesh the stiffness metrix of all elements are same.
And stiffness metrix obtain from above weak form and shape function
(interpolation function) is given by,
2
k b 2
K e =
[
6 a 1
1
2
2
1
1
1
1
2
2
K e =
1
2 1
1] + k a [ 1
2
2
6 b 1 2
2
2 1
kb e
ka e
[X ij ] +
[Y ]
6a
6 b ij
kb
ka
[ X] +
[ Y]
6a
6b
1
2
2
1
2
1]
1
2
2
0
4 2
2
4
0 2
0
0
1
0
2 1
1
2
0 1
0
0
0
0
0
0
0
0
0
0
0
0
2
1
0
1
4
1
0
1
4
0
0
1
0
0
0
2 1
0
1 4 1
[Y] = 0 1 4
0
0 1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
[ 0
0
0
0
0
1 1
0
0
0
0
0
0
0
0
0 1
2 1
0
0
0
0
0
0
2
0
0 1
2 1
0
0
0
0
0
4 2
0
0 1
2 1
0
0
0
0
2
2
0
0
0 1
1
0
0
0
0
0
0
4 4
0
0
0
1 1
0
0
0
0 4
8 4
0
0 1
2 1
0
1
0
0 4
8 4
0
0 1
2 1
2 1
0
0 4
8 4
0
0 1
2
1
1
0
0
0 4
4
0
0
0 1
0
0
1 1
0
0
0 2
2
0
0
0
0 1
2 1
0
0
2
4 2
0
0
0
0 1
2 1
0
0 2
4 2
0
0
0
0 1
2 1
0
0 2
4
0
0
0
0
0 1
1
0
0
0 2
0
0 2 1
0
0
0 1 4 1
1
0
0 1 4
4
1
0
0 1
1
2
0
0
0
0
0
4
2
0
0
0
2
8
2
1
0
0
2
8
4 1
0
0
2
1 2
0
0
0
0
0 2 1
0
0
0 1 4 1
0
0
0 1 4
0
0
0
0 1
0
0
0
0
0
0
0
1
4
1
0
0
2
8
2
0
0
1
4
1
Where,
[X] and [Y] are in global.
Now in our case, for simplicity I have taken,
Kx ,ky = k = 30 W/ m C
423 = 823 =18 W/ m2 C
0
0
0
0
0
0
0
0
1
1
0
0
0
2
2]
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
2
0
0
0
0
0
0 2 1
0
0
0
0 1 4 1
0
0
0
0 1 4 1
0
2
0
0 1 4 1
4
0
0
0 1 2
0
2
1
0
0
0
0
1
4
1
0
0
0
0
1
4
1
0
1
0
0
1
4
1
2
0
0
0
1
2]
[H] =
e23 he23
6
0
0
0
0
0
0
0
0
0
0
0
0
0
0
[0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
2
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
4
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
2
0]
{[K] + [H]} =
20 5
5
40
0
5
0
0
0
0
5 10
10 10
0 10
0
0
0
0
0
0
0
0
0
0
0
0
[ 0
0
0
0
5
0
40
5
5
40
0
5
0
0
10
0
10 10
10 10
0
10
0
0
0
0
0
0
0
0
0
0
0
5 10
0
10 10
0
0
10
5
0
0
20.6
0
0
0
40 10
0
10 80
0
0
10
10
0
0
4.7
0
0
0
5 10
0
10 10
0
0
10
0
0
0
0
0
0
0
10
10
10
0
0
10
80
10
0
0
10
10
10
0
0
0
10
10
0
0
0
10
80
10
0
0
10
10
10
0
0
0
0
0
0
0
0
0
0
0
0
10
0
0
0
4.7
0
0
0
0
5 10
0
0
10 10 10
0
0
10 10
10
0
0
10
41.2
0
0
0
0
0
15
0
0
5
40
5
0
0
5
40
9.7
0
0
5
5.0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
10
0
10 10
10 4.7
0
0
0
0
5
0
40 4.4
5
20 ]
2.6 Solution:
So that we can obtain finally,
20 5
0
0
5
40
5
0
0
5
40
5
0
0
5
40
0
0
0
5
5 10
0
0
10 10 10
0
0 10 10 10
0
0 10 10
0
0
0
10
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
[ 0
0
0
0
Q1
1
0
Q2
2
0
Q3
0
2
Q4
0
2
Q5
1
9
0
2
0
0
0
4
= 4 + 0 + 0
0
4
0
0
18
2
Q
1
0
11
Q12
2
0
Q13
0
2
2
0
Q14
[1] [ 9 ] [Q15 ]
40
10
0
0
[ 0
10
0
80 10
10 80
0
10
0
0
0
0
0
5
20.6
0
0
0
10
4.7
0
0
0
0
0
5
10
0
0
0
40
10
0
0
0
5
10
0
0
0
10
10
10
0
0
10
80
10
0
0
10
10
10
0
0
0
10
10
10
0
0
10
80
10
0
0
10
10
10
0
0
0
10
10
0
0
0
10
80
10
0
0
10
10
10
T6
0
0
1502
T7
0
0
3004
T
10
0
8 = 3004
80 10 T9
3004
]
[
10 41.2 [T10 ]
1454]
0
0
0
10
4.7
0
0
0
10
41.2
0
0
0
9.7
5.0
0
0
0
0
0
5
10
0
0
0
0
5
0
0
0
0
0
0
0
0
10
10
10
0
0
15
40
5
0
0
0
0
0
0
0
0
10
10
10
0
0
5
40
5
0
0
0
0
0
0
0
0
10
10
10
0
0
5
40
5
0
0
0
0
0
0
0
0
0
0
T6
0
T7
0
T8
0
10 T9
4.7 T10
0
100
0
100
0
100
4.4 100
20 ] [100]
By using the simple matlab program we can find the value of the
temperature
and
heat
flux
at
the
various
nodes,
T7 = 50.0610 C
T8 = 50.0230 C
T9 =49.7229 C
T10 = 47.3599 C
And heat flux at node 1 to 5 and node 11 to 15 are
Q1 = -751.9 W/m2
Q2 = -1503.5 W/m2
Q3 = -1500.1 W/m2
Q4 = -1473.1 W/m2
Q5 = -230.6 W/m2
Q11 = 748.1 W/m2
Q12 = 1496.5 W/m2
Q13 = 1499.9 W/m2
Q14 = 1601.1 W/m2
Q15 = 774 W/m2
Chapter 3
Analysis
3.1 software name and version:
Ansys 15.0 workbench
Ansys 15.0 mechanical APDL
Modeling software: Creo parametric 2.0
Modeling :
2. Temperature distribution:
3. Heat flux:
Chapter 4
Results and Discussions
4.1 Comparisons of analytical and FEM modeling
results:
Analytical result of the problem compare with the result of software.
Node no.
6
7
8
9
10
Analytical result of
Temperature(C)
50.0653
50.0610
50.0230
49.7229
47.3599
Table 2 Comparisons of results
Software results of
temperature(C)
47
46
46.5
45
44.44
Chapter 5
Conclusion
Here we note that the result of analysis of problem using software is
much closer to the value found analytically. There is amount of error
that taken into consideration. This error in the Temperature distribution
is due to our law accuracy of calculation with limited amount of mesh
taken into consideration while computer software is work with the high
degree of accuracy in calculation. If we take more number of element for
given problem our answer will go closer to the actual one.