Scribd
Upload
44 views

AM5360: Composite Structures Assignment 2: Submitted By, Anilkumar P M Reg. No: CE18D755

This document contains MATLAB code that analyzes the behavior of stiffness and compliance terms for a composite laminate as it undergoes in-plane rotation. It defines material properties and the compliance matrix S for an orthotropic lamina. It then uses a transformation matrix T to determine the transformed compliance matrix S1 for different rotation angles. Several plots are generated to show how terms in S1 such as Ex/E1, νxy, and Gxy vary with the rotation angle.

Uploaded by

Anilkmar P M
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
44 views

AM5360: Composite Structures Assignment 2: Submitted By, Anilkumar P M Reg. No: CE18D755

This document contains MATLAB code that analyzes the behavior of stiffness and compliance terms for a composite laminate as it undergoes in-plane rotation. It defines material properties and the compliance matrix S for an orthotropic lamina. It then uses a transformation matrix T to determine the transformed compliance matrix S1 for different rotation angles. Several plots are generated to show how terms in S1 such as Ex/E1, νxy, and Gxy vary with the rotation angle.

Uploaded by

Anilkmar P M
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 9

AM5360: Composite Structures

Assignment 2

Submitted by,
Anilkumar P M
Reg. No: CE18D755
MATLAB CODE

%%
%%MATERIAL PROPERTIES AND S MATRIX
E1 = 145*10^9;
E2 = 20*10^9;
E3 = 12*10^9;
nu12 = 0.3;
nu21 = 0.0428;
nu13 = 0.25;
nu31 = 0.0214;
nu23 = 0.20;
nu32 = 0.12;
G12 = 6*10^9;
G13 = 5*10^9;
G23 = 4*10^9;
S= [(1/E1) (-nu12/E2) (-nu31/E3) 0 0 0; (-nu12/E1) (1/E2) (-nu32/E3) 0 0 0;
(-nu13/E1) (-nu23/E2) (1/E3) 0 0 0; 0 0 0 (1/G12) 0 0; 0 0 0 0 (1/G23) 0; 0
0 0 0 0 (1/G13)]
%%
%%TRANSFORMATION MATRIX

syms thettax1
%thettay1 = pi/2-thettax1;
%thettaz1= pi/2; thettaz2=pi/2; thettaz3=0;
%thettax2= thettax1+ pi/2; thettay2= thettay1+ pi/2;
m1= cos (thettax1);m2= -sin (thettax1); m3= 0;
n1= sin (thettax1);n2= cos(thettax1); n3= 0;
p1= 0;p2= 0; p3= 1;

T= [m1^2 n1^2 p1^2 2*n1*p1 2*m1*p1 2*m1*n1;


m2^2 n2^2 p2^2 2*n2*p2 2*m2*p2 2*m2*n2;
m3^2 n3^2 p3^2 2*n3*p3 2*m3*p3 2*m3*n3;
m2*m3 n2*n3 p2*p3 n2*p3+n3*p2 p2*m3+p3*m2 m2*n3+m3*n2;
m1*m3 n1*n3 p1*p3 n1*p3+n3*p1 p1*m3+p3*m1 m1*n3+m3*n1;
m2*m1 n2*n1 p2*p1 n2*p1+n1*p2 p2*m1+p1*m2 m2*n1+m1*n2]

S1= inv(T)*S*T

xthetta= 0:pi/48:pi/2;

figure(1)
e1= double(subs((inv(S1(1,1))/ E1), thettax1, xthetta))
%e2= double(subs((inv(S1(2,2))/E2), thettax1, xthetta))
%e3= double(subs((inv(S1(3,3))/E3), thettax1, xthetta))
xthettap= 180*xthetta/pi;
plot(xthettap, e1)
ylabel('E_x/E_1')
xlabel('Rotation angle')

figure(2)
g1= subs((S1(6,6)^-1)/G12, thettax1, xthetta)
plot(xthettap, g1)
ylabel('G_x_y/G_1_2')
xlabel('Rotation angle')
figure(3)
g2= subs((S1(4,4)^-1)/G23, thettax1, xthetta)
plot(xthettap, g2)
ylabel('G_y_z/G_2_3')
xlabel('Rotation angle')
figure(4)
g3= subs((S1(5,5)^-1)/G13, thettax1, xthetta)
plot(xthettap, g3)
ylabel('G_x_z/G_1_3')
xlabel('Rotation angle')

figure(5)
nu1= subs((S1(2,1)*S1(1,1)^-1), thettax1, xthetta)
plot(xthettap, -nu1)
ylabel('?_x_y')
xlabel('Rotation angle')
figure(6)
nu2= subs((S1(3,1)*S1(1,1)^-1), thettax1, xthetta)
plot(xthettap, -nu2)
ylabel('?_x_z')
xlabel('Rotation angle')
figure(7)
nu3= subs((S1(2,3)*S1(3,3)^-1), thettax1, xthetta)
plot(xthettap, -nu3)
ylabel('?_y_z')
xlabel('Rotation angle')

figure(8)
g1= subs((S1(1,6))*G12, thettax1, xthetta)
plot(xthettap, g1)
ylabel('?_x_y_,_x')
xlabel('Rotation angle')
figure(9)
g1= subs((S1(2,6))*G23, thettax1, xthetta)
plot(xthettap, g1)
ylabel('?_x_y_,_y')
xlabel('Rotation angle')
figure(10)
g1= subs((S1(3,6))*G13, thettax1, xthetta)
plot(xthettap, g1)
ylabel('?_x_y_,_z')
xlabel('Rotation angle')

figure(11)
g1= subs(((S1(6,1))*S1(6,6)^-1), thettax1, xthetta)
plot(xthettap, g1)
ylabel('?_z_x_,_x_y')
xlabel('Rotation angle')

You might also like