%% Function name
%       QandS

%% Revised:
%       16 January 2014

%% Author
%       Sotirios Chatzisavvas, Trey Moore, & Autar Kaw
%       Section: All
%       Semester: Fall 2012

%% Purpose
%       Given the elastic modulii, Poisson's ratio, and shear modulus
%       of a unidirectional lamina, output the Reduced Stiffness Matrix,and
%       the Compliance Matrix.

%% Usage:
%       function [Q_S]= QandS(moduli)
%   Input Variables
%       moduli=vector with four elastic modulii of unidirectional lamina
%       [moduli]=[E1 E2 nu12 G12]
%       E1=longitudinal elastic modulus
%       E2=transverse elastic modulus
%       nu12=major Poisson's ratio
%       G12=in-plane shear modulus
%   Output Variables
%       [Q]=vector of reduced stiffness matrix
%       [S]=compliance matrix
%   Keyword
%       reduced stiffnes matrix
%       compliance matrix

%% License Agreement
%       http://www.eng.usf.edu/~kaw/OCW/composites/license/limiteduse.pdf

%% Testing code
clc
clear all
%% Inputs
% Material: Graphite/Epoxy
moduli=[181 10.3 0.28 7.17];
fprintf('\nLongitudinal Elastic Modulus: %g ',moduli(1))
fprintf('\nTransverse Elastic Modulus: %g ',moduli(2))
fprintf('\nPoisson''s Ratio: %g ',moduli(3))
fprintf('\nShear Modulus: %g \n\n\n',moduli(4))
%% Test
% Testing for reduced stiffnes matrix and compliance matrix
% Call function: Q,S
[Q,S]= QandS(moduli);
disp('               Reduced Stiffness Matrix')
disp('__________________________________________________')
disp(Q)
disp(' ')
disp('               Compliance Matrix')
disp('__________________________________________________')
disp(S)