% Script til Oblig 1 for MEK-4540 Komposittmaterialer h?st 2008
% Brian Hayman 2008
%
% Materialegenskapene samt laminatoppbyggingen m? byttes ut med de riktige dataene.
%
clear all;

% Materialdata
EL=40000;
ET=7700;
vLT=0.25;
GLT=2800;


R=[45 0 -45]*3.14159/180;   % orientering av lag i radianer
h=[-.96 -.64 .64 .96];  % z-koordinater for lagenes topp/bunnflater som definert i fig. 6-5


% Compliance (S-matrix)
S(1,1)=1/EL;
S(2,2)=1/ET;
S(1,2)=-vLT/EL;
S(3,3)=1/GLT;
S(2,1)=S(1,2);

% Stiffness (Q-matrix)
Q=inv(S);

Q(3,3)=2*Q(3,3);        % for tensorielle t?ynigner ved rotasjon


A(1:3,1:3)=0;
B(1:3,1:3)=0;
D(1:3,1:3)=0;

for i=1:length(R)                 % Beregner T og Q for hvert lag 'k'. Summerer A, B og D 

    % Transformation matrix (T-matrix)
    T(1,1)=cos(R(i))^2;
    T(2,2)=cos(R(i))^2;
    T(1,2)=sin(R(i))^2;
    T(2,1)=sin(R(i))^2;
    T(3,1)=-sin(R(i))*cos(R(i));
    T(3,2)=sin(R(i))*cos(R(i));
    T(1,3)=2*sin(R(i))*cos(R(i));
    T(2,3)=-2*sin(R(i))*cos(R(i));
    T(3,3)=cos(R(i))^2-sin(R(i))^2;
    
    Qk(:,:,i)=inv(T)*Q*T;             % Stivhetsmatrise for tensorielle t?yninger
    Qk(:,3,i)=Qk(:,3,i)/2;            % G?r tilbake til ingeni?rt?yninger etter rotasjon
    A=A+Qk(:,:,i)*(h(i+1)-h(i));                % Beregner A-matrisen
    B=B+.5*Qk(:,:,i)*(h(i+1)^2-h(i)^2);         % Beregner B-matrisen
    D=D+1/3*Qk(:,:,i)*(h(i+1)^3-h(i)^3);        % Beregner D-matrisen
end

A
B
D

% Total stivhetsmatrise for platen gir sammenheng mellom spenninger og ingeni?rt?yninger
Ktot=[A B
    B D]