Scribd
Upload
56 views1 page

Nonlinear Finite Element Analysis: Mema/Cses 648 Spring 2002

This document provides an assignment to develop a computer program to solve nonlinear transient heat conduction problems using finite element analysis. Students are asked to revise their subroutine to account for time approximation and solve a specific problem of heat conduction in a rectangular isotropic medium with given boundary and initial conditions. Key aspects include: - Developing a program to solve a nonlinear partial differential equation governing heat conduction - Accounting for linear and quadratic elements in the program - Validating the program by solving a problem of heat conduction in a medium with temperature-dependent thermal conductivity - Applying given boundary, initial conditions, material properties to the problem for simulation

Uploaded by

stefan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
56 views1 page

Nonlinear Finite Element Analysis: Mema/Cses 648 Spring 2002

This document provides an assignment to develop a computer program to solve nonlinear transient heat conduction problems using finite element analysis. Students are asked to revise their subroutine to account for time approximation and solve a specific problem of heat conduction in a rectangular isotropic medium with given boundary and initial conditions. Key aspects include: - Developing a program to solve a nonlinear partial differential equation governing heat conduction - Accounting for linear and quadratic elements in the program - Validating the program by solving a problem of heat conduction in a medium with temperature-dependent thermal conductivity - Applying given boundary, initial conditions, material properties to the problem for simulation

Uploaded by

stefan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

MEMA/CSES 648 Spring 2002

Nonlinear Finite Element Analysis

Date: March 28, 2002 Due: April 5, 2002

ASSIGNMENT No. 7 (Program No. 6)


Nonlinear Transient Analysis of
Heat Conduction Type Problems

Consider the equation

∂u ∂ ∂u ∂ ∂u
² ³ ² ³
ρC1 − axx − ayy + a00 u = f (x, y, t) in Ω (1)
∂t ∂x ∂x ∂y ∂y

where axx and ayy are, in general, functions of x, y, and u. Write a computer program to study
linear as well as nonlinear solutions of Eq. (1). Assume that C1 , a00 , axx , and ayy are of the form

C1 = C10 + C1X ∗ x + C1Y ∗ Y


a00 = A00 + A0X ∗ x + A0Y ∗ Y
axx = AX0 + AXX ∗ x + AXY ∗ Y + AXU ∗ U
ayy = AY 0 + AY X ∗ x + AY Y ∗ Y + AY U ∗ U (2)

Your program should account for linear and quadratic rectangular Lagrange family of elements.
In particular, revise subroutine ELKMFR to account for the time approximation, and validate
the program by solving the following problem (corresponds to heat conduction in a rectangular,
isotropic medium):

Problem. The domain is of dimensions a = 0.18m and b = 0.1m along the x and y coordinates,
respectively; internal heat generation is zero f = 0; the conductivity k = a11 = a22 is of the form

k = k0 (1 + βT ) (3)

where k0 is the constant thermal conductivity and β is the temperature coefficient of thermal
conductivity. Take k0 = 0.2 W/(m.◦ C) and β = 2 × 10−3 (◦ C−1 ). The boundary conditions are

∂T
T (0, y, t) = 500◦ C , T (a, y, t) = 300◦ C , = 0 at y = 0, b (4)
∂y

The initial condition is


T (x, y, 0) = 0◦ C (5)
Take ρC1 = 1.0, and use a time step (you may have to estimate ∆t) that is reasonable for the
problem.

You might also like