close all;clear all; clc;

%Rozmìry experimentu
l=300;
hh=0.7;
hj=39.5;
hd=0.7;
b=91.5;
F=1e3;
dT=[28.7246;28.7246;28.7246];

%Parametry materiálù
Eh=2e5;
Gh=77e3;
alfah=17.3e-6;
vh=(Eh-2*Gh)/2/Gh;

Ej=12.6864;
alfaj=4.5e-5;
vj=0.5682;
Gj=Ej/(2*vj+2);

Ed=73.1e3;
Gd=28e3;
alfad=23.2e-6;
vd=(Ed-2*Gd)/2/Gd;

%Matice tuhostí
%Horní potah
C11=Eh*(1-vh)/(1+vh)/(1-2*vh);
C12=Eh*vh/(1+vh)/(1-2*vh);

Ch=[C11    C12   C12    0               0      0;
    C12    C11   C12    0               0      0;
    C12    C12   C11    0               0      0;
     0      0     0   1/2*(C11-C12)     0      0;
     0      0     0     0        1/2*(C11-C12) 0;
     0      0     0     0               0  1/2*(C11-C12);];
 
%Jádro
C11=Ej*(1-vj)/(1+vj)/(1-2*vj);
C12=Ej*vj/(1+vj)/(1-2*vj);

Cj=[C11    C12   C12    0               0      0;
    C12    C11   C12    0               0      0;
    C12    C12   C11    0               0      0;
     0      0     0   1/2*(C11-C12)     0      0;
     0      0     0     0        1/2*(C11-C12) 0;
     0      0     0     0               0  1/2*(C11-C12);];
 
%Dolní potah
C11=Ed*(1-vd)/(1+vd)/(1-2*vd);
C12=Ed*vd/(1+vd)/(1-2*vd);

Cd=[C11    C12   C12    0               0      0;
    C12    C11   C12    0               0      0;
    C12    C12   C11    0               0      0;
     0      0     0   1/2*(C11-C12)     0      0;
     0      0     0     0        1/2*(C11-C12) 0;
     0      0     0     0               0  1/2*(C11-C12);];
 
%Vzdálenosti
h1=-hh-hj/2;
h2=-hj/2;
h3=hj/2;
h4=hd+hj/2;
 
%Matice A
A11=Ch(1,1)*(h2-h1)+Cj(1,1)*(h3-h2)+Cd(1,1)*(h4-h3);
A12=Ch(1,2)*(h2-h1)+Cj(1,2)*(h3-h2)+Cd(1,2)*(h4-h3);
A16=Ch(1,6)*(h2-h1)+Cj(1,6)*(h3-h2)+Cd(1,6)*(h4-h3);
A21=Ch(2,1)*(h2-h1)+Cj(2,1)*(h3-h2)+Cd(2,1)*(h4-h3);
A22=Ch(2,2)*(h2-h1)+Cj(2,2)*(h3-h2)+Cd(2,2)*(h4-h3);
A26=Ch(2,6)*(h2-h1)+Cj(2,6)*(h3-h2)+Cd(2,6)*(h4-h3);
A61=Ch(6,1)*(h2-h1)+Cj(6,1)*(h3-h2)+Cd(6,1)*(h4-h3);
A62=Ch(6,2)*(h2-h1)+Cj(6,2)*(h3-h2)+Cd(6,2)*(h4-h3);
A66=Ch(6,6)*(h2-h1)+Cj(6,6)*(h3-h2)+Cd(6,6)*(h4-h3);

A=[A11     A12     A16;
   A21     A22     A26;
   A61     A62     A66;];

 
%Matice AT
AT11=alfah*Ch(1,1)*(h2-h1)+alfaj*Cj(1,1)*(h3-h2)+alfad*Cd(1,1)*(h4-h3);
AT12=alfah*Ch(1,2)*(h2-h1)+alfaj*Cj(1,2)*(h3-h2)+alfad*Cd(1,2)*(h4-h3);
AT16=alfah*Ch(1,6)*(h2-h1)+alfaj*Cj(1,6)*(h3-h2)+alfad*Cd(1,6)*(h4-h3);
AT21=alfah*Ch(2,1)*(h2-h1)+alfaj*Cj(2,1)*(h3-h2)+alfad*Cd(2,1)*(h4-h3);
AT22=alfah*Ch(2,2)*(h2-h1)+alfaj*Cj(2,2)*(h3-h2)+alfad*Cd(2,2)*(h4-h3);
AT26=alfah*Ch(2,6)*(h2-h1)+alfaj*Cj(2,6)*(h3-h2)+alfad*Cd(2,6)*(h4-h3);
AT61=alfah*Ch(6,1)*(h2-h1)+alfaj*Cj(6,1)*(h3-h2)+alfad*Cd(6,1)*(h4-h3);
AT62=alfah*Ch(6,2)*(h2-h1)+alfaj*Cj(6,2)*(h3-h2)+alfad*Cd(6,2)*(h4-h3);
AT66=alfah*Ch(6,6)*(h2-h1)+alfaj*Cj(6,6)*(h3-h2)+alfad*Cd(6,6)*(h4-h3);

AT=[AT11     AT12     0;
    AT21     AT22     0;
    AT61     AT62     0;];

%Matice B
B11=Ch(1,1)*(h2^2-h1^2)+Cj(1,1)*(h3^2-h2^2)+Cd(1,1)*(h4^2-h3^2);
B12=Ch(1,2)*(h2^2-h1^2)+Cj(1,2)*(h3^2-h2^2)+Cd(1,2)*(h4^2-h3^2);
B16=Ch(1,6)*(h2^2-h1^2)+Cj(1,6)*(h3^2-h2^2)+Cd(1,6)*(h4^2-h3^2);
B21=Ch(2,1)*(h2^2-h1^2)+Cj(2,1)*(h3^2-h2^2)+Cd(2,1)*(h4^2-h3^2);
B22=Ch(2,2)*(h2^2-h1^2)+Cj(2,2)*(h3^2-h2^2)+Cd(2,2)*(h4^2-h3^2);
B26=Ch(2,6)*(h2^2-h1^2)+Cj(2,6)*(h3^2-h2^2)+Cd(2,6)*(h4^2-h3^2);
B61=Ch(6,1)*(h2^2-h1^2)+Cj(6,1)*(h3^2-h2^2)+Cd(6,1)*(h4^2-h3^2);
B62=Ch(6,2)*(h2^2-h1^2)+Cj(6,2)*(h3^2-h2^2)+Cd(6,2)*(h4^2-h3^2);
B66=Ch(6,6)*(h2^2-h1^2)+Cj(6,6)*(h3^2-h2^2)+Cd(6,6)*(h4^2-h3^2);

B=[B11     B12     B16;
   B21     B22     B26;
   B61     B62     B66;];

B=1/2*B;

%Matice BT
BT11=alfah*Ch(1,1)*(h2^2-h1^2)+alfaj*Cj(1,1)*(h3^2-h2^2)+alfad*Cd(1,1)*(h4^2-h3^2);
BT12=alfah*Ch(1,2)*(h2^2-h1^2)+alfaj*Cj(1,2)*(h3^2-h2^2)+alfad*Cd(1,2)*(h4^2-h3^2);
BT16=alfah*Ch(1,6)*(h2^2-h1^2)+alfaj*Cj(1,6)*(h3^2-h2^2)+alfad*Cd(1,6)*(h4^2-h3^2);
BT21=alfah*Ch(2,1)*(h2^2-h1^2)+alfaj*Cj(2,1)*(h3^2-h2^2)+alfad*Cd(2,1)*(h4^2-h3^2);
BT22=alfah*Ch(2,2)*(h2^2-h1^2)+alfaj*Cj(2,2)*(h3^2-h2^2)+alfad*Cd(2,2)*(h4^2-h3^2);
BT26=alfah*Ch(2,6)*(h2^2-h1^2)+alfaj*Cj(2,6)*(h3^2-h2^2)+alfad*Cd(2,6)*(h4^2-h3^2);
BT61=alfah*Ch(6,1)*(h2^2-h1^2)+alfaj*Cj(6,1)*(h3^2-h2^2)+alfad*Cd(6,1)*(h4^2-h3^2);
BT62=alfah*Ch(6,2)*(h2^2-h1^2)+alfaj*Cj(6,2)*(h3^2-h2^2)+alfad*Cd(6,2)*(h4^2-h3^2);
BT66=alfah*Ch(6,6)*(h2^2-h1^2)+alfaj*Cj(6,6)*(h3^2-h2^2)+alfad*Cd(6,6)*(h4^2-h3^2);

BT=[BT11     BT12     0;
    BT21     BT22     0;
    BT61     BT62     0;];

BT=1/2*BT;

%Matice D
D11=Ch(1,1)*(h2^3-h1^3)+Cj(1,1)*(h3^3-h2^3)+Cd(1,1)*(h4^3-h3^3);
D12=Ch(1,2)*(h2^3-h1^3)+Cj(1,2)*(h3^3-h2^3)+Cd(1,2)*(h4^3-h3^3);
D16=Ch(1,6)*(h2^3-h1^3)+Cj(1,6)*(h3^3-h2^3)+Cd(1,6)*(h4^3-h3^3);
D21=Ch(2,1)*(h2^3-h1^3)+Cj(2,1)*(h3^3-h2^3)+Cd(2,1)*(h4^3-h3^3);
D22=Ch(2,2)*(h2^3-h1^3)+Cj(2,2)*(h3^3-h2^3)+Cd(2,2)*(h4^3-h3^3);
D26=Ch(2,6)*(h2^3-h1^3)+Cj(2,6)*(h3^3-h2^3)+Cd(2,6)*(h4^3-h3^3);
D61=Ch(6,1)*(h2^3-h1^3)+Cj(6,1)*(h3^3-h2^3)+Cd(6,1)*(h4^3-h3^3);
D62=Ch(6,2)*(h2^3-h1^3)+Cj(6,2)*(h3^3-h2^3)+Cd(6,2)*(h4^3-h3^3);
D66=Ch(6,6)*(h2^3-h1^3)+Cj(6,6)*(h3^3-h2^3)+Cd(6,6)*(h4^3-h3^3);

D=[D11     D12     D16;
   D21     D22     D26;
   D61     D62     D66;];

D=1/3*D;

%Matice N
N=[0;
   0;
   0 ];
%Matice M
M=[F*l/4/b;
   0;
   0 ];

disp('Laminatova teorie s vlivem teploty');
%Prùhyb bez smyku
disp('bez smyku');

X=inv(B*inv(A)*-B+D)*(BT*dT-B*A^-1*AT*dT);
w_ft=-X(1,1)/2*l^2/4+X(1,1)*l^2/4;

X=inv(B*inv(A)*-B+D);
w_f=F*l^3/(48*b)*X(1,1)+w_ft

%Prùhyb se smykem
disp('se smykem')
F45=0;
F54=0;
F55=Ch(5,5)*hh+Cj(5,5)*hj+Cd(5,5)*hd;
F44=Ch(4,4)*hh+Cj(4,4)*hj+Cd(4,4)*hd;
F55hv=F44/F44/F55;

w_fs=F*l^3/48/b*X(1,1)*(1+12*F55hv/X(1,1)/l^2)+w_ft


%Matice teplotní roztažnosti
%alfa=[alfa;
