%% Function name
%      AsBsDs

%% Revised
%      28 January 2014

%% Author
%      Robert Worcester, Trey Moore, & Autar Kaw
%      Section: All
%      Semester: Fall 2013

%% Purpose
%       Given the thickness of each ply and the extensional, coupling, and
%       bending matrices, output the extensional compliance, coupling
%       compliance, bending compliance, normalized extensional compliance,
%       normalized coupling compliance, and normalized bending compliance
%       matrices

%% Usage
%      function [As,Bs,Ds,Asn,Bsn,Dsn] = AsBsDs(A,B,D,tplies)
% Input variables
%       tplies=thickness of each ply
%       A=extensional stiffness matrix
%       B=coupling stiffness matrix
%       D=bending stiffness matrix
% Output variables
%       As=extensional compliance matrix
%       Bs=coupling compliance matrix
%       Ds=bending compliance matrix
%       Asn=normalized extensional compliance matriix
%       Bsn=normalized coupling compliance matrix
%       Dsn=normalized bending compliance matrix
% Keywords
%       extensional compliance matrix 
%       coupling compliance matrix
%       bending compliance matrix
%       normalized extensional compliance matrix
%       normalized coupling compliance matrix
%       normalized bending compliance matrix

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

%% Testing Code
clc
clear all
%% Inputs
tplies=[0.000125 0.000125 0.000125 0.000125];
A=[2.833e7,2.116e7,0;2.116e7,2.883e7,0;0,0,2.33e7];
B=[0,0,0;0,0,0;0,0,0];
D=[5.902e-1,4.408e-1,3.349e-1;4.408e-1,5.902e-1,3.349e-1;3.349e-1,3.349e-1,4.853e-1];
for i=1:1:length(tplies)
    fprintf('For Ply: %g\n',i)
    fprintf('   Ply Thickness: %g\n\n',tplies(i))
end
disp('Extensional Stiffness Matrix:')
disp('______________________________________________')
disp(A)
disp('Coupling Stiffness Matrix:')
disp('______________________________________________')
disp(B)
disp('Bending Stiffness Matrix:')
disp('______________________________________________')
disp(D)
%% Test
% Call function: AsBsDs
[As,Bs,Ds,Asn,Bsn,Dsn] = AsBsDs(A,B,D,tplies);
% Extensional compliance matrix
disp('Extensionsl Compliance Matrix:')
disp('______________________________________________')
disp(As)
% Coupling compliance matrix
disp('Coupling Compliance Matrix:')
disp('______________________________________________')
disp(Bs)
% Bending compliance matrix
disp('Bending Compliance Matrix:')
disp('______________________________________________')
disp(Ds)
% Normalized extensional compliance matrix
disp('Normalized Extensional Compliance Matrix:')
disp('______________________________________________')
disp(Asn)
% Normalized coupling compliance matrix
disp('Normalized Coupling Compliance Matrix:')
disp('______________________________________________')
disp(Bsn)
% Normalized bending compliance matrix
disp('Normalized Bending Compliance Matrix:')
disp('______________________________________________')
disp(Dsn)