Notes For Finite Element Method Thermal Analysis
Notes For Finite Element Method Thermal Analysis
Vable
du dx
Work Potential:
L m
W A = p x u dx +
0
Fq u ( xq )
q=1
WT =
qx T dx + Qq T ( xq )
0 q=1
Functional
Potential Energy A = UA WA
n
T = UT WT
n
Ci fi ( x )
i=1
T(x) =
Ci fi ( x )
i=1 L
df j df k = EA dx d x d x
0 L m
K jk =
L
k d x d x dx
0 m
df j df k
Rj =
px fj
0
dx +
Fq fj ( xq )
q=1
R j = q x f j dx +
0
Qq fj ( xq )
q=1
Other Applications:
Flow through pipes; Flow through porous media; Electrostatics By non-dimensionalizing the problem a software can be used to solve all the above applications.
M. Vable
Linear Elements:
T = T1 L1 ( x ) + T2 L2 ( x ) [K
(e) (e) (e)
Q1
(e)
(e)
Q2
(e)
(e)
k ] = --------) 1 1 (e L 1 1
(e) (e)
{R
[K
(e)
k A ] = ------------------ 1 1 (e) 1 1 L
{R
(e)
Quadratic Elements:
T = T1 L1 ( x ) + T2 L2 ( x ) + T3 L3 ( x ) 7 8 1 k (e) [ K ] = ------------ 8 16 8 (e) 3L 1 8 7
(e) (e) (e) (e)
Q1
(e) (e)
Q3
(e) (e) 1 qx L
Class Problem 1
Heat qo is being added at point B at a constant rate as shown in Fig. 1. Using two linear elements determine the temperature at point B and the heat flowing out at A and C in terms of k, L, and qo for the following two cases: (a) The ends of the bar are maintained at a constant zero temperature. (b) The ends of the bars are maintained at a constant temperature To.
qo A L B C 2L A L 3qo B L qo C L D
Fig .1
Fig .2
Home Problem 1
Heat 3qo is being added at point B and qo is being taken out at a constant rate at point C as shown in Fig. 2. Using three linear elements determine the temperature at points B and C and the heat flowing out at A and D in terms of k, L, and qo for the following two cases: (a) The ends of the bar are maintained at a constant zero temperature. (b) The ends of the bars A and D are maintained at a temperature To and 2To, respectively. ANS: TB = ( 5qo L ) 3k T C = ( q o L ) 3k q A = ( 5q o 3 )
M. Vable
Thermal Stresses
Stress-strain curve with temperature effects
Normal Stress
O T
O1
-E
Normal Strain
= -- + T = -- + o E E
where, = linear coefficient of thermal expansion. o = Initial strain= Thermal strain unconstrained body due to uniform temperature changes.
U AT
1 = -- xx ( xx o ) dV = 2
V L L
-- E ( xx o ) 2
0 L
A dx or
L L
U AT
1 1 du 2 du 1 2 2 = -- EA xx dx EA xx o dx + -- EA o dx = -- EA dx EA o dx + U o d x d x 2 2 2
0 L 0 m 0 0 0
WA =
px ( x )u ( x ) dx + Fq u ( xq )
0 L q=1
A = U AT W A
1 du 2 du = -- EA dx EA o dx + U o d x d x 2
0 0
px ( x )u ( x ) dx + Fq u ( xq )
0 q=1
M. Vable
1 du 2 = -- EA dx + U o d x 2
0 L m
px ( x )u ( x ) dx + Fq u ( xq ) + EAo d x dx
du
0 L q=1 0
= U A W AT
W AT =
L
px ( x )u ( x ) dx + Fq u ( xq ) + EAo d x dx
du
0 q=1 L L 0 0
du I = EA o dx = EA o u ( x ) d x
0
d o EA u dx = EATu ( x ) dx
0 L
L L 0
d EA T u dx dx
0
W AT =
where,
F T1 EAT 1 = F T2 EAT 2 qx = 0
EA p xT = ----------- q x k
Thermal loads are added only at the element ends. Linear Element:
[K
(e)
EA ] = ------- 1 1 L 1 1
{R
(e)
Quadratic Element:
(e)
[K
7 8 1 EA ] = ------- 8 16 8 3L 1 8 7
F( e ) T1 0 + F( e ) T3
M. Vable
Class Problem 2
Heat qo is being added at point B at a constant rate as shown in Fig. 3. Using two quadratic elements determine the displacement of node B, reaction force at A and the axial stress just before B in terms of E, A, , k, L, and qo. Assume that the entire bar was at zero temperature and the ends of the bars are maintained at zero temperature.
qo A L B C 2L A L 3qo B L qo C L D
Fig .3
Fig .4
Home Problem 2
Heat 3qo is being added at point B and qo is being taken out at a constant rate at point C as shown in Fig. 4. Using three linear elements determine the displacements of points B and C, the reaction force at A, and the axial stress just before B in terms of E, A, , k, L, and qo. Assume that the entire bar was at zero temperature and the ends of the bars are maintained at zero temperature. ANS: u B
5q o L = -----------------18k
2
uC
7q o L = -----------------18k
5EAq o L R A = -----------------------9k
M. Vable
or
k T = q v
t dx dy
Isotropic material:
[k] = k
---scaler quantity.
Element Approximation:
(e) T1 n T(e ) 2 (e) T ( x ) = T i f i ( x, y ) = f 1 f 2 f n i=1 (e) Tn
UT =
-- { d 2
(e) T
} [B
(e) T
] [k
(e)
][B
(e)
]{d
(e)
M. Vable
[B
A
(e) T
] [k
(e)
][B
(e)
]t dx dy
where,
n= direction of the unit normal to the boundary. nx, ny = direction cosines of the unit normal to the boundary. qn = specified heat flow in the n-direction on the boundary.
where,
= the boundary of the element. s= tangential coordinate along the element boundary.
T = h ( Tf T ) n
(e)
where,
h = convection heat transfer coefficients Tf = Temperature of the surrounding fluid h depends upon many factors: velocity of fluid, viscosity of fluid, density of fluid, and other properties of fluid. It also depends upon the surface roughness and surface geometry. Addition to Element Matrix: Addition to Element Right Hand Side Vector:
K ij = Ri
(e) (e)
(e)
f i f j t ds f i T f t ds
(e)
(e)
(e)
If
is an element boundary in the interior, then there is no convection there and hence no addition to the matrix or the element RHS vector.
(e)
M. Vable
where,
BC is the proportionality constant. Tr is the temperature of the other radiating body. Temperatures Tr and T are in absolute degrees i.e., oK.
For two infinite parallel black bodies (planes) it is called the Boltzmann
constant. For regular bodies BC depends upon the emissivity of the bodies, the geometry and other factors. Radiation boundary conditions lead to non-linear thermal problem. A general approach is:
k T 2 2 = BC ( Tr + T ) ( Tr + T ) ( Tr T ) = hr ( Tr T ) n
In the iteration process, at each step treat the radiation boundary condition like a convection term with coefficient dependent upon the temperature at a particular step.