Scribd
Upload
8 views

MCM Assignment 1

The document presents code to calculate various elastic properties and ultimate strength properties of composite materials with varying volume fractions of fiber (Vf), using empirical, semi-empirical and analytical approaches. The properties calculated include longitudinal and transverse Young's modulus (E1, E2), Poisson's ratios (u12, u21), shear modulus (G12) and ultimate tensile, compressive and shear strengths. Plots of each property against Vf are also generated.

Uploaded by

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

MCM Assignment 1

The document presents code to calculate various elastic properties and ultimate strength properties of composite materials with varying volume fractions of fiber (Vf), using empirical, semi-empirical and analytical approaches. The properties calculated include longitudinal and transverse Young's modulus (E1, E2), Poisson's ratios (u12, u21), shear modulus (G12) and ultimate tensile, compressive and shear strengths. Plots of each property against Vf are also generated.

Uploaded by

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

MCM Assignment 1A: Elastic Properties Using SOM Approach

Name: Sakshi Vinayak Tonde


C.No. : C22021552707

Code
%Sakshi Tonde
%C22021552707
%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 = ');

%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);
a=(Em*Vf(i))+(Ef*Vm);
E2(i)=(Ef*Em)/a;
u12(i)=(Vf(i)*vf)+(Vm*vm);
u21(i)=(u12(i)*E2(i))/E1(i);
b=(Gm*Vf(i))+(Gf*Vm);
G12(i)=(Gf*Gm)/b;
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

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 3.7611 4.2079 4.7753 5.5195 6.5385


8.0189 10.3659 14.6552 25.0000 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.0944 0.0597 0.0462 0.0398 0.0370


0.0368 0.0394 0.0469 0.0683 0.2000

G12 =
1.3077 1.4471 1.6197 1.8390 2.1271 2.5223
3.0977 4.0132 5.6971 9.8152 35.4167

>>

Graph (Empirical Approach)

Fig1: Plot of elastic moduli properties w.r.t volume fraction of fiber using
SOM approach
MCM Assignment 1B: Elastic Properties Using Halpin-Tsai Approach

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)

Fig2: Plot of Elastic moduli properties w.r.t volume fraction of fiber


using Halpin-Tsai approach
MCM Assignment 1C: Ultimate Strength Properties

Name: Sakshi Vinayak Tonde


C.No. : C22021552707

Code
%Sakshi Tonde
%C22021552707
%Ultimate Tensile Strength

%Inputs
close all;clear all;clc;
Ef=input('Enter Ef = ');
Em=input('Enter Em = ');
vf=input('Enter mu f = ');
vm=input('Enter mu m = ');
sigma_FT_ult = input('Enter sigma FT ult = ');
sigma_MT_ult = input('Enter sigma MT ult = ');
sigma_MC_ult = input('Enter sigma MC ult = ');
tau_F_ult = input('Enter tau F ult = ');
tau_M_ult = input('Enter tau M ult = ');
d_by_s = input('Enter d/s = ');

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

for i=1:n

Vm=1-Vf(i);
E1(i)=(Ef*Vf(i))+(Em*Vm);
u12(i)=(Vf(i)*vf)+(Vm*vm);
a=(Em*Vf(i))+(Ef*Vm);
E2(i)=(Ef*Em)/a;

%sigma_1C_ult
ep1(i) = epslon_MT_ult*(1-sqrt(Vf(i)));
ep2 = epslon_MT_ult*(d_by_s*((Em/Ef)-1)+1);
if ep1(i)<ep2
epslon_2T_ult(i)=ep1(i);
else
epslon_2T_ult(i)=ep2;
end
sigma_1C_ult(i) = (E1(i)*epslon_2T_ult(i))/u12(i);

%sigma_1T_ult
sigma_1T_ult(i) = (sigma_FT_ult*Vf(i))+(epslon_FT_ult*Em*(1-
(Vf(i))));

%sigma_2C_ult
epslon_2C_ult = epslon_MC_ult*((d_by_s*(Em/Ef))+(1-d_by_s));
sigma_2C_ult(i) = E2(i)*epslon_2C_ult;

%sigma_2T_ult
epslon_2T_ult = epslon_MT_ult*((d_by_s*(Em/Ef))+(1-d_by_s));
sigma_2T_ult(i) = E2(i)*epslon_2T_ult;

%tau_12_ult
b=(Gm*Vf(i))+(Gf*Vm);
G12(i)=(Gf*Gm)/b;
Gma_12_M_ult = tau_M_ult/Gm;
Gma_12_ult = Gma_12_M_ult*((d_by_s*(Gm/Gf))+(1-d_by_s));
tau_12_ult = G12*Gma_12_ult;

end

%Results
sigma_1C_ult
sigma_1T_ult
sigma_2C_ult
sigma_2T_ult
tau_12_ult

%Plots

%sigma_1C_ult VS Vf
subplot(3,2,1)
plot(Vf,sigma_1C_ult,'b')
grid on;
xlabel('Volume fraction of fibre Vf','FontSize',8)
ylabel('Longitudinal Compressive Strength','FontSize',8)
legend('sigma1Cult vs Vf')
title('Longitudinal compressive strength Vs Vf')
hold on;

%sigma_1T_ult VS Vf
subplot(3,2,2)
plot(Vf,sigma_1T_ult,'r')
grid on;
xlabel('Volume fraction of fibre Vf','FontSize',8)
ylabel('Longitudinal tensile strength','FontSize',8)
legend('sigma1Tult vs Vf')
title('Longitudinal tensile strength Vs Vf')
hold on;

%sigma_2C_ult VS Vf
subplot(3,2,3)
plot(Vf,sigma_2C_ult,'g')
grid on;
xlabel('Volume fraction of fibre Vf','FontSize',8);
ylabel('Transverse compressive strength','FontSize',8);
legend('sigma2Cult vs Vf')
title('Transverse compressive strength Vs Vf')
hold on;

%sigma_2T_ult VS Vf
subplot(3,2,4)
plot(Vf,sigma_2T_ult,'b')
grid on;
xlabel('Volume fraction of fibre Vf','FontSize',8);
ylabel('Transverse tensile strength','FontSize',8);
legend('sigma2Tult vs Vf')
title('Transverse tensile strength Vs Vf')
hold on;

%tau_12_ult VS Vf
subplot(3,2,5.5)
plot(Vf,tau_12_ult,'g')
grid on;
xlabel('Volume fraction of fiber Vf','FontSize',8)
ylabel('Inplane shear strength','FontSize',8);
legend('tau12ult vs Vf')
title('In plane shear strength Vs Vf')

Code Output
Enter Ef = 85000
Enter Em = 3400
Enter mu f = 0.2
Enter mu m = 0.3
Enter sigma FT ult = 1550
Enter sigma MT ult = 72
Enter sigma MC ult = 102
Enter tau F ult = 34
Enter tau M ult = 35
Enter d/s = 0.9441
Vf =

Columns 1 through 9

0 0.1 0.2 0.3 0.4


0.5 0.6 0.7 0.8

Columns 10 through 11

0.9 1

Gf =

35417

Gm =

1307.7

epslon_MT_ult =

0.021176

epslon_FT_ult =

0.018235

epslon_MC_ult =

0.03

sigma_1C_ult =

Columns 1 through 9

22.479 79.065 139.69 204.81 274.94


350.68 432.73 521.91 619.2

Columns 10 through 11

397.63 0
sigma_1T_ult =

Columns 1 through 9

62 210.8 359.6 508.4 657.2


806 954.8 1103.6 1252.4

Columns 10 through 11

1401.2 1550

sigma_2C_ult =

Columns 1 through 9

9.5537 10.568 11.824 13.418 15.509


18.373 22.532 29.127 41.18

Columns 10 through 11

70.248 238.84

sigma_2T_ult =

Columns 1 through 9

6.7438 7.46 8.3463 9.4716 10.948


12.969 15.905 20.56 29.068

Columns 10 through 11

49.587 168.6

tau_12_ult =

Columns 1 through 9

3.1766 3.5151 3.9344 4.4673 5.1671


6.1269 7.5247 9.7487 13.839

Columns 10 through 11

23.843 86.032

>>
Graph (Ultimate Strength Properties)

Fig3: Plot of ultimate strength properties w.r.t volume fraction of fiber

You might also like