Scribd
Upload
32 views

MCM Assignment 2

The document presents the results of a semi-empirical approach to calculate the Young's modulus, Poisson's ratio, and in-plane shear modulus based on the volume fraction of fibers in a composite material. It provides the code to input material properties and calculate the results, which are then plotted against volume fraction for analysis.

Uploaded by

2707SAKSHI TONDE
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
32 views

MCM Assignment 2

The document presents the results of a semi-empirical approach to calculate the Young's modulus, Poisson's ratio, and in-plane shear modulus based on the volume fraction of fibers in a composite material. It provides the code to input material properties and calculate the results, which are then plotted against volume fraction for analysis.

Uploaded by

2707SAKSHI TONDE
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

MCM Assignment 1: Semi-Empirical Approach to Find Youngs

Modulus, Poisson’s Ratio & In-Plane Shear Modulus

Name: Sakshi Vinayak Tonde


C.No. : C22021552707

Code
%Sakshi Tonde
%C22021552707
%Semi-Empirical Approach to find E1, E2, u12, u21, G12 wrt Vf

%Inputs
close all;clear all;clc;
Ef=input('Enter Ef = ');
Em=input('Enter Em = ');
vf=input('Enter mu f = ');
vm=input('Enter mu m = ');
zeta1=input('Enter zeta1 = ');

%Calculations
n=11;
Vf=[0:0.1:1]
Gf=Ef/(2*(1+vf))
Gm=Em/(2*(1+vm))

for i=1:n
Vm=1-Vf(i);
E1(i)=(Ef*Vf(i))+(Em*Vm);
eta1=((Ef/Em)-1)/((Ef/Em)+zeta1);
E2(i)= Em*((1+(zeta1*eta1*Vf(i))))/(1-(eta1*Vf(i)));
u12(i)=(Vf(i)*vf)+(Vm*vm);
u21(i)= (u12(i)*E2(i))/E1(i);
if (Vf(i)>0.3)
zeta2=1+(40*(Vf(i)^10));
else
zeta2=1;
end
eta2=((Gf/Gm)-1)/((Gf/Gm)+zeta2);
G12(i)= Gm*(1+(zeta2*eta2*Vf(i)))/(1-(eta2*Vf(i)));
end

%Results
E1
E2
u12
u21
G12

%Plots
%E1 VS Vf
subplot(3,2,1)
plot(Vf,E1,'b')
grid on;
xlabel('Volume fraction of fibre Vf')
ylabel('Longitudinal Youngs Modulus GPa ')
legend('E1 vs Vf')
title('Variation wrt Youngs Modulus')
hold on;

%E2 VS Vf
subplot(3,2,2)
plot(Vf,E2,'r')
grid on;
xlabel('Volume fraction of fibre Vf')
ylabel('Transverse Youngs Modulus GPa')
legend('E2 vs Vf')
title('Variation wrt Youngs Modulus')
hold on;

%u12 vs Vf
subplot(3,2,3)
plot(Vf,u12,'g')
grid on;
xlabel('Volume fraction of fibre Vf');
ylabel('Minor Poissons Ratio');
legend('u12 vs Vf')
title('Variation wrt Poissons Ratio')
hold on;

%u21 vs Vf
subplot(3,2,4)
plot(Vf,u21,'b')
grid on;
xlabel('Volume fraction of fibre Vf');
ylabel('Major Poissons Ratio');
legend('u21 vs Vf')
title('Variation wrt Poissons Ratio')
hold on;

%G12 VS Vf
subplot(3,2,5.5)
plot(Vf,G12,'g')
grid on;
xlabel('Volume fraction of fiber Vf')
ylabel('In plane Shear modulus Gpa');
legend('G12 vs Vf')
title('Variation wrt Shear Modulus')

Code Output
Enter Ef = 85

Enter Em = 3.4

Enter mu f = 0.2

Enter mu m = 0.3

Enter zeta1 = 2
Vf =

0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000


0.8000 0.9000 1.0000

Gf =

35.4167

Gm =

1.3077

E1 =

3.4000 11.5600 19.7200 27.8800 36.0400 44.2000 52.3600 60.5200


68.6800 76.8400 85.0000

E2 =

3.4000 4.3951 5.6054 7.1091 9.0276 11.5600 15.0571 20.2000


28.5077 44.2000 85.0000

u12 =

0.3000 0.2900 0.2800 0.2700 0.2600 0.2500 0.2400 0.2300


0.2200 0.2100 0.2000
u21 =

0.3000 0.1103 0.0796 0.0688 0.0651 0.0654 0.0690 0.0768


0.0913 0.1208 0.2000

G12 =

1.3077 1.5755 1.9044 2.3180 2.8566 3.6136 4.9274 8.1293


16.2293 27.6906 35.4167

>>

Graph (Semi-Empirical Approach)

Fig: Plot of E1/E2 Vs Vf, u12/u21 Vs Vf & G12 Vs Vf (Semi-Empirical


approach)

You might also like