Scribd
Upload
28 views

Matrix Data

This MATLAB function takes in material property values and a theta value for a composite material. It calculates and outputs two matrices - the first contains the stiffness and compliance matrices, and the second contains the reduced compliance, reduced stiffness, Q-bar and S-bar matrices for the composite. It determines these values by calculating the individual S, C, Q values from the input properties then defines and populates the various matrices that are output.

Uploaded by

Austin Stott
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
28 views

Matrix Data

This MATLAB function takes in material property values and a theta value for a composite material. It calculates and outputs two matrices - the first contains the stiffness and compliance matrices, and the second contains the reduced compliance, reduced stiffness, Q-bar and S-bar matrices for the composite. It determines these values by calculating the individual S, C, Q values from the input properties then defines and populates the various matrices that are output.

Uploaded by

Austin Stott
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

3/3/14 6:32 PM

%
%
%
%

C:\Users\Austin\Desktop\Composites Project\matrix_data.m

This function file will take in given material properties and a theta
value. It will then output 2 matrices. The first matrix will contain the
stiffness and compliance matrices. The second matrix will contain the
reduced compliance, reduced stiffness, Q_bar and S_bar matrices.

function [A,B] = matrix_data(E1,E2,E3,v23,v13,v12,G23,G13,G12,theta)


format long eng;
% Engineering Format
S11= 1/E1;
% Determine the S values
S12= -v12/E1;
S13= -v13/E1;
S22= 1/E2;
S23= -v23/E2;
S33= 1/E3;
S44= 1/G23;
S55= 1/G13;
S66= 1/G12;
S = S11*S22*S33-S11*S23*S23-S22*S13*S13-S33*S12*S12+2*S12*S23*S13;
C11= (S22*S33-S23*S23)/S;
% Determine the C values
C12= (S13*S23-S12*S33)/S;
C13= (S12*S23-S13*S22)/S;
C22= (S33*S11-S13*S13)/S;
C23= (S12*S13-S23*S11)/S;
C33= (S11*S22-S12*S12)/S;
C44=1/S44;
C55=1/S55;
C66=1/S66;
Q11 = (S22 / (S11*S22 - S12^2));
% Determine the Q values
Q12 = (-S12 / (S11*S22 - S12^2));
Q22 = (S11 / (S11*S22 - S12^2));
Q66 = 1 / S66;
m=cosd(theta);
% m&n values for T matrix
n=sind(theta);
% Calculate S_bar
S_bar_11 = (S11*(m^4))+((2*S12+S66)*(n^2)*(m^2))+(S22*(n^4));
S_bar_12 =(S11+S22-S66)*(n^2)*(m^2)+S12*(n^4+m^4);
S_bar_16 = (2*S11-2*S12-S66)*n*(m^3)-(2*S22-2*S12-S66)*(n^3)*m;
S_bar_22 = S11*(n^4)+(2*S12+S66)*(n^2)*(m^2)+S22*(m^4);
S_bar_26 = (2*S11-2*S12-S66)*(n^3)*m-(2*S22-2*S12-S66)*n*(m^3);
S_bar_66 = (4*S11+4*S22-8*S12-2*S66)*(n^2)*(m^2)+S66*((n^4)+(m^4));
% Calculate Q_bar
Q_bar_11 = Q11*(m^4)+(2*Q12+4*Q66)*(n^2)*(m^2)+Q22*(n^4);
Q_bar_12 =(Q11+Q22-4*Q66)*(n^2)*(m^2)+Q12*(n^4+m^4);
Q_bar_16 = (Q11-Q12-2*Q66)*n*(m^3)+(Q12-Q22+2*Q66)*(n^3)*m;
Q_bar_22 = Q11*(n^4)+(2*Q12+4*Q66)*(n^2)*(m^2)+Q22*(m^4);
Q_bar_26 = (Q11-Q12-2*Q66)*(n^3)*m+(Q12-Q22+2*Q66)*n*(m^3);
Q_bar_66 = (Q11+Q22-2*Q12-2*Q66)*(n^2)*(m^2)+Q66*((n^4)+(m^4));
% Define matrices
S_matrix = [S11 S12 S13 0 0 0; S12 S22 S23 0 0 0; S13 S23 S33 0 0 0;...
0 0 0 S44 0 0; 0 0 0 0 S55 0; 0 0 0 0 0 S66];
C_matrix = [C11 C12 C13 0 0 0; C12 C22 C23 0 0 0; C13 C23 C33 0 0 0;...
0 0 0 C44 0 0; 0 0 0 0 C55 0; 0 0 0 0 0 C66];

1 of 2

3/3/14 6:32 PM

C:\Users\Austin\Desktop\Composites Project\matrix_data.m

Q_matrix = [Q11 Q12 0; Q12 Q22 0; 0 0 Q66];


S_bar_matrix = [S_bar_11,S_bar_12,S_bar_16; S_bar_12,S_bar_22,S_bar_26; ...
S_bar_16,S_bar_26,S_bar_66];
Q_bar_matrix = [Q_bar_11,Q_bar_12,Q_bar_16; Q_bar_12,Q_bar_22,Q_bar_26; ...
Q_bar_16,Q_bar_26,Q_bar_66];
S_matrix_reduced = [S11 S12 0; S12 S22 0; 0 0 S66];
% Output matrices
A(:,:,1) = S_matrix;
A(:,:,2) = C_matrix;
B(:,:,1) = S_bar_matrix;
B(:,:,2) = Q_bar_matrix;
B(:,:,3) = S_matrix_reduced;
B(:,:,4) = Q_matrix;

2 of 2

You might also like