Matlab Codes For FEM Problems
Uploaded by
Abid YusufMatlab Codes For FEM Problems
Uploaded by
Abid Yusuf`
Reg.No: 20MMF0004_______________________
Name : Abid
Yusuf____________________________
Slot : L3-4
SCHOOL OF MECHANICAL ENGINEERING
FINAL ASSESEMENT TEST 2020
COURSE NAME : FEM in Manufacturing SEMESTER : FALL
COURSE CODE : MEE 5001 PROGRAMME : M.Tech (MFG)
MAXIMUM MARKS : 40 DURATION : 90 minutes.
CLASS NUMBER : VL2020210104223 DATE : 09.11.2020 TIME : 10.00 to 11.30 AM
COMPONENT MARKS
1 Algorithm & Code 15
2 Solution and Output 10
3 Viva 15
40
Total
A tapered bar is of circular cross section and is of length “L”. It has a larger diameter of 80
mm and a smaller diameter of 40 mm. The 80 mm diameter side is fixed and a tensile load of “P” is
applied axially at the 40 mm side. The Youngs modulus of the material is “E”. Find
1. Assuming 2 elements show the element stiffness Matrix.
2. Solve the problem using 2, 20, 200, 2,000 and 20,000 elements.
3. Compare the free end deflection solutions using different number of elements with the closed of the
4𝑃𝐿
given problem given by the expression .
𝜋𝐸𝐷𝑑
4. Provide your inferance
X is the last digit of your registration number
Take L = 100 * (X+1) in mm =500mm
P = 10 * (X+1) in kN =50kN
E = 200 GPa
My CODE
% Axial deformation in a tapered bar
L=500;
E=200000;
D=80;
d=40;
Load=50000;
N=2000;
l=L/N;
K= zeros(N+1,N+1);
for i = 1:N
di= D-((D-d)*(2*i-1)/(2*N));
aebyl= (pi*di*di*E)/(4*l);
ki= [aebyl,-aebyl;-aebyl,aebyl];
K(i,i)= K(i,i)+ki(1,1);
K(i,i+1)=K(i,i+1)+ki(1,2);
K(i+1,i)=K(i+1,i)+ki(2,1);
K(i+1,i+1)=K(i+1,i+1)+ki(2,2);
end
for a = 1:N
for k= 1:N
KK(a,k)=K(a+1,k+1);
end
end
P=zeros(N,1);
P(N,1)=P(N,1)+Load
U=KK\P
Solutions:
1: Element Stiffness Matrix for two elements
K=
1.0e+06 *
3.0788 -3.0788 0
-3.0788 4.6496 -1.5708
0 -1.5708 1.5708
2: Output:
No. of Elements Solution Refinement at free end (Deflection)
2 0.0481
20 0.0497
200 0.0497
2000 0.0497
20000 0.0497
3: The closed form of the solution gives the value: 0.049715909
Our Finite Element analysis of this bar has returned an exactly the same solution upto four decimal
(quite accurate) places at 20 elements itself, our study up until 20,000 elements just strengthens our finding.