0% found this document useful (0 votes)
8 views1 page

Assignment - 2

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
0% found this document useful (0 votes)
8 views1 page

Assignment - 2

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

Department of Applied Mechanics and Biomedical Engineering, Indian Institute of Technology, Madras.

Fundamentals of Finite Element Analysis (AM 5450)


Assignment 1 – Weighted Residual Methods - I
Submission Date: 26 August 2024

Instructions:
I. Use scientific non-programmable calculator wherever necessary.
II. MATLAB/PYTHON/Python scripts used to prepare plots/figures/results should be attached with your
answer.
III. Scan all the pages of your hand-written answers, MATLAB/PYTHON scripts, output results, and figures into
a PDF file. Please submit your assignment in Google class room.

1) Find the approximate solution 𝑢̃(𝑥) of the following differential equation using the least square and Galerkin’s
methods by hand calculations. Obtain solutions by assuming polynomials of order 1 to 3. Include both domain
and boundary residual as applicable.
𝑑2 𝑢
− 6𝑥 = 0
𝑑𝑥 2
with the boundary conditions given as, 𝑢(0) = 0 and 𝑢(3) = 27.
a) Plot the solution 𝑢̃(𝑥) and its derivative 𝑢̃′(𝑥) as a function of 𝑥 in MATLAB/PYTHON for each polynomial
order 1 to 3. Also compute the exact solution using analytical method and compare its plot in
MATLAB/PYTHON with the solutions obtained by least square method.
b) Plot the residual function 𝑒(𝑥) by substituting back the obtained 𝑢̃(𝑥) for each polynomial order 1 to 3.

2) Find the approximate solution 𝑢̃(𝑥) of the following differential equation using the least square and Galerkin’s
methods in MATLAB/PYTHON. Obtain solutions by assuming polynomials of order 1 to 𝑁. Include both domain
and boundary residue as applicable.
𝑑2 𝑢
+ 𝑘 2𝑢 = 0
𝑑𝑥 2
with the boundary conditions given as, 𝑢(0) = 0 and 𝑢’(1) = 𝑢. Assume 𝑘 = 5 and 𝑢 = 10.
a) Find the highest polynomial order 𝑁 at which the assumed solution 𝑢̃(𝑥) converges to exact solution 𝑢(𝑥)
and the residual function 𝑒(𝑥) = 1E-6 (approximates to 0) for all values of 𝑥 in the domain.
b) Plot the solution 𝑢̃(𝑥) as a function of 𝑥 in MATLAB/PYTHON for each polynomial order 1 to 𝑁. Also compare
the plot of exact solution given by 𝑢(𝑥) = (2 𝑠𝑖𝑛(5𝑥))/𝑐𝑜𝑠(5) in MATLAB/PYTHON with the solutions
obtained by least square method.
c) Plot the residual function 𝑒(𝑥) by substituting back the obtained 𝑢̃(𝑥) for each polynomial order 1 to 𝑁.

Page 1 of 1

You might also like