Scribd
Upload
46 views9 pages

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
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
46 views9 pages

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
We take content rights seriously. If you suspect this is your content, claim it here.
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