%Script til Oblig 1 for MEK-4540 Komposittmaterialer v?r 2007
%Terje Glomsaker 2007
%
clear all;

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


R=[45 0 -45]*3.1415/180;   %orientering av lag i radianer
h=[-.96 -.64 .64 .96];    %tykkelse av lagene


%Compliance (S-matirx)
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

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