ASSIGNMENT

MEHNAZ RASOOL ENTRYNO. 2012AMZ8406

COURSE:MECHANICS OF COMPOSITES
function compositefailure
sigma1ft= 1500e6;
sigma1fc= 1200e6;
sigma2ft= 50e6;
sigma2fc= 100e6;
sigma12f= 90e6;
E1=170e9;
E2=10e9;
G12=13e9;
m12=0.3;
alpha11=-0.9e-6;
alpha22=27e-6;
beta11=150e-6;
beta22=4800e-6;
Tc=175;
Tr=20;
dM=0.5;
n=3;
h=0.000125;
theta=[30 0 45];
Nx=1;
Ny=0;
Nxy=0;
Mx=0;
My=0;
Mxy=0;
[sigmax1,sigmay1,sigmaxy1]= ...
mehnazfailure1(E1,E2,G12,m12,alpha11,alpha22,beta11,beta22,...
Tc,Tr,dM,n,h,theta,0,0,0,0,0,0);
[sigmax2,sigmay2,sigmaxy2]= ...
mehnazfailure1(E1,E2,G12,m12,0,0,0,0,0,0,0,n,h,theta,Nx,Ny,Nxy,Mx,My,Mxy);
for i=1:n
fprintf('\nFor ply-%d at %d degrees:',i,theta(i));
m= cosd(theta(i));
n= sind(theta(i));
inT1= [m^2 n^2 -2*m*n; n^2 m^2 2*m*n; m*n -m*n m^2-n^2];
Sft= inT1*[sigma1ft; sigma2ft; sigma12f];
Sxft= Sft(1);
Syft= Sft(2);
Sxyf= Sft(3);
inT2= [m^2 n^2 -2*m*n; n^2 m^2 2*m*n];
Sfc= inT2*[sigma1fc; sigma2fc; sigma12f];
Sxfc= Sfc(1);
Syfc= Sfc(2);
%At bottom,
Res1t_b= (Sxft-sigmax1(2*i-1))/sigmax2(2*i-1);
Res2t_b= (Syft-sigmay1(2*i-1))/sigmay2(2*i-1);
Res12_b= (Sxyf-sigmaxy1(2*i-1))/sigmaxy2(2*i-1);
Res1c_b= (-Sxfc-sigmax1(2*i-1))/(-sigmax2(2*i-1));
Res2c_b= (-Syfc-sigmay1(2*i-1))/(-sigmay2(2*i-1));
%At top,

Res1t_t= (Sxft-sigmax1(2*i))/sigmax2(2*i);
Res2t_t= (Syft-sigmay1(2*i))/sigmay2(2*i);
Res12_t= (Sxyf-sigmaxy1(2*i))/sigmaxy2(2*i);
Res1c_t= (-Sxfc-sigmax1(2*i))/(-sigmax2(2*i));
Res2c_t= (-Syfc-sigmay1(2*i))/(-sigmay2(2*i));
Resultant=[abs(Res1t_b);abs(Res2t_b);abs(Res12_b);abs(Res1c_b);...
abs(Res2c_b);abs(Res1t_t);abs(Res2t_t);abs(Res12_t);...
abs(Res1c_t);abs(Res2c_t)];
Res(i)= min(Resultant);
end
Res_fail=Res(1);
index=1;
for j=2:n
if Res(j)<Res_fail
Res_fail=Res(j);
index=j;
end
end
fprintf('\n\n\nMinimum value of Nx for the laminate= %f',Res_fail);
end

function [sigmax,sigmay,sigmaxy]= ...

mehnazfailure1(E1,E2,G12,m12,alpha11,alpha22,beta11,beta22,...
Tc,Tr,dM,n,h,theta,Nx,Ny,Nxy,Mx,My,Mxy)
A= zeros(3,3);
B= zeros(3,3);
D= zeros(3,3);
Nt=zeros(3,1);
Mt=zeros(3,1);
Nm=zeros(3,1);
Mm=zeros(3,1);
dT= Tr-Tc;
z(1)= -(n*(h/2));
for i=1:n
z(i+1)= z(i) + h;
m21= m12*(E2/E1);
Q= zeros(3,3);
Q(1,1)= E1/(1-m12*m21);
Q(2,2)= E2/(1-m12*m21);
Q(3,3)= G12;
Q(1,2)= m21*Q(1,1);
Q(2,1)= m21*Q(1,1);
m= cosd(theta(i)); n= sind(theta(i));
T= [m^2 n^2 2*m*n; n^2 m^2 -2*m*n; -m*n m*n m^2-n^2];
R= [1 0 0; 0 1 0; 0 0 2];
Qxy(:,:,i)= inv(T)*Q*R*T*inv(R);
A= A + (z(i+1)-z(i)).*Qxy(:,:,i);
B= B + (0.5)*((z(i+1))^2-(z(i))^2).*Qxy(:,:,i);
D= D + (1/3)*((z(i+1))^3-(z(i))^3).*Qxy(:,:,i);

alpha_xx(i)= alpha11*(cosd(theta(i)))^2 + alpha22*(sind(theta(i)))^2;

alpha_yy(i)= alpha11*(sind(theta(i)))^2 + alpha22*(cosd(theta(i)))^2;
alpha_xy(i)= 2*sind(theta(i))*cosd(theta(i))*(alpha11-alpha22);
alphaxy(:,:,i)= [alpha_xx(i); alpha_yy(i); alpha_xy(i)];
Nt= Nt+ (z(i+1)-z(i))*dT.*(Qxy(:,:,i)*alphaxy(:,:,i));
Mt= Mt+ (0.5)*((z(i+1))^2-(z(i))^2)*dT.*(Qxy(:,:,i)*alphaxy(:,:,i));
beta_xx(i)= beta11*(cosd(theta(i)))^2 + beta22*(sind(theta(i)))^2;
beta_yy(i)= beta11*(sind(theta(i)))^2 + beta22*(cosd(theta(i)))^2;
beta_xy(i)= 2*sind(theta(i))*cosd(theta(i))*(beta11-beta22);
betaxy(:,:,i)= [beta_xx(i); beta_yy(i); beta_xy(i)];
Nm= Nm+ (z(i+1)-z(i))*dM.*(Qxy(:,:,i)*betaxy(:,:,i));
Mm= Mm+ (0.5)*((z(i+1))^2-(z(i))^2)*dM.*(Qxy(:,:,i)*betaxy(:,:,i));
end
d= inv(D-B*inv(A)*B);
b= -inv(A)*B*d;
a= inv(A)-b*B*inv(A);
ep= [a, b; b', d]*([Nx; Ny; Nxy; Mx; My; Mxy] + [Nt; Mt] + [Nm; Mm]);
exo= ep(1); eyo= ep(2);
exyo= ep(3);
kx= ep(4);
ky= ep(5);
kxy= ep(6);
for i=1:(n+1)
ex(i)= exo+z(i)*kx;
ey(i)= eyo+z(i)*ky;
exy(i)=exyo+z(i)*kxy;
end
for i=1:n
q1(1:3)= Qxy(1,1:3,i);
q2(1:3)= Qxy(2,1:3,i);
q3(1:3)= Qxy(3,1:3,i);
sigmax_bottom= q1(1:3)*( [ex(i);ey(i);exy(i)] - dT.*alphaxy(:,:,i) dM.*betaxy(:,:,i) );
sigmax_top= q1(1:3)*( [ex(i+1);ey(i+1);exy(i+1)] - dT.*alphaxy(:,:,i) dM.*betaxy(:,:,i) );
sigmay_bottom= q2(1:3)*( [ex(i);ey(i);exy(i)] - dT.*alphaxy(:,:,i) dM.*betaxy(:,:,i) );
sigmay_top= q2(1:3)*( [ex(i+1);ey(i+1);exy(i+1)] - dT.*alphaxy(:,:,i) dM.*betaxy(:,:,i) );
sigmaxy_bottom=q3(1:3)*( [ex(i);ey(i);exy(i)] - dT.*alphaxy(:,:,i) dM.*betaxy(:,:,i) );
sigmaxy_top=q3(1:3)*( [ex(i+1);ey(i+1);exy(i+1)] - dT.*alphaxy(:,:,i) dM.*betaxy(:,:,i) );
sigmax(2*i-1)=sigmax_bottom;
sigmax(2*i)=sigmax_top;
sigmay(2*i-1)=sigmay_bottom;
sigmay(2*i)=sigmay_top;
sigmaxy(2*i-1)=sigmaxy_bottom;
sigmaxy(2*i)=sigmaxy_top;
end

end

Results:For ply-1 at 30 degrees:

Minimum value of stress resultant for ply-1= 443261.477321
For ply-2 at 0 degrees:
Minimum value of stress resultant for ply-2= 503660.279814
For ply-3 at 45 degrees:
Minimum value of stress resultant for ply-3=

Minimum value of stress resultant for the laminate= 185405.976097