clear all;
clc;

fid11=fopen('Q11.txt','w');
fid12=fopen('Q12.txt','w');
fid13=fopen('Q13.txt','w');
fid14=fopen('Q14.txt','w');
fid15=fopen('Q15.txt','w');
fid16=fopen('Q16.txt','w');
fid21=fopen('Q21.txt','w');
fid22=fopen('Q22.txt','w');
fid23=fopen('Q23.txt','w');
fid24=fopen('Q24.txt','w');
fid25=fopen('Q25.txt','w');
fid26=fopen('Q26.txt','w');
fid31=fopen('Q31.txt','w');
fid32=fopen('Q32.txt','w');
fid33=fopen('Q33.txt','w');
fid34=fopen('Q34.txt','w');
fid35=fopen('Q35.txt','w');
fid36=fopen('Q36.txt','w');
fid41=fopen('Q41.txt','w');
fid42=fopen('Q42.txt','w');
fid43=fopen('Q43.txt','w');
fid44=fopen('Q44.txt','w');
fid45=fopen('Q45.txt','w');
fid46=fopen('Q46.txt','w');
fid51=fopen('Q51.txt','w');
fid52=fopen('Q52.txt','w');
fid53=fopen('Q53.txt','w');
fid54=fopen('Q54.txt','w');
fid55=fopen('Q55.txt','w');
fid56=fopen('Q56.txt','w');
fid61=fopen('Q61.txt','w');
fid62=fopen('Q62.txt','w');
fid63=fopen('Q63.txt','w');
fid64=fopen('Q64.txt','w');
fid65=fopen('Q65.txt','w');
fid66=fopen('Q66.txt','w');



 E1=sym('E1');
 E2=sym('E2');
 E3=sym('E3');
 G12=sym('G12');
 G13=sym('G13');
 G23=sym('G23');
 v12=sym('v12');
 v21=sym('v21');
 v13=sym('v13');
 v31=sym('v31');
 v23=sym('v23');
 v32=sym('v32');
 mx=sym('mx');
 nx=sym('nx');
 my=sym('my');
 ny=sym('ny');
 mz=sym('mz');
 nz=sym('nz');
 
 class(E1);
 class(E2);
 class(E3);
 class(G12);
 class(G13);
 class(G23);
 class(v12);
 class(v21);
 class(v13);
 class(v31);
 class(v23);
 class(v32);
 class(mx);
 class(nx);
 class(my);
 class(ny);
 class(mz);
 class(nz);
 
Q(1,1)=-(E1^2*(E3 - E2*v32^2))/(E1^2*v31^2 + 2*E2*E1*v12*v31*v32 + E2*E1*v32^2 - E3*E1 + E2*E3*v12^2);
Q(1,2)=-(E1*E2*(E3*v12 + E1*v31*v32))/(E1^2*v31^2 + 2*E2*E1*v12*v31*v32 + E2*E1*v32^2 - E3*E1 + E2*E3*v12^2);
Q(1,3)=-(E1*E3*(E1*v31 + E2*v12*v32))/(E1^2*v31^2 + 2*E2*E1*v12*v31*v32 + E2*E1*v32^2 - E3*E1 + E2*E3*v12^2);
Q(1,4)=0;
Q(1,5)=0;
Q(1,6)=0;
Q(2,1)=Q(1,2);
Q(2,2)=-(E1*E2*(E3 - E1*v31^2))/(E1^2*v31^2 + 2*E2*E1*v12*v31*v32 + E2*E1*v32^2 - E3*E1 + E2*E3*v12^2);
Q(2,3)=-(E1*E2*E3*(v32 + v12*v31))/(E1^2*v31^2 + 2*E2*E1*v12*v31*v32 + E2*E1*v32^2 - E3*E1 + E2*E3*v12^2);
Q(2,4)=0;
Q(2,5)=0;
Q(2,6)=0;
Q(3,1)=Q(1,3);
Q(3,2)=Q(2,3);
Q(3,3)=-(E3^2*(E1 - E2*v12^2))/(E1^2*v31^2 + 2*E2*E1*v12*v31*v32 + E2*E1*v32^2 - E3*E1 + E2*E3*v12^2);
Q(3,4)=0;
Q(3,5)=0;
Q(3,6)=0;
Q(4,1)=Q(1,4);
Q(4,2)=Q(2,4);
Q(4,3)=Q(3,4);
Q(4,4)=G23;
Q(4,5)=0;
Q(4,6)=0;
Q(5,1)=Q(1,5);
Q(5,2)=Q(2,5);
Q(5,3)=Q(3,5);
Q(5,4)=Q(4,5);
Q(5,5)=G13;
Q(5,6)=0;
Q(6,1)=Q(1,6);
Q(6,2)=Q(2,6);
Q(6,3)=Q(3,6);
Q(6,4)=Q(4,6);
Q(6,5)=Q(5,6);
Q(6,6)=G12;
 
%Transformationsmatrix der Spannungen bei Drehung um die X-Achse
TSig_alpha=[1  0     0      0         0      0;
    0  mx^2  nx^2   2*mx*nx   0      0;
    0  nx^2  mx^2  -2*mx*nx   0      0;
    0 -mx*nx mx*nx  mx^2-nx^2 0      0;
    0  0     0      0         mx    -nx;
    0 -mx*nz mx*nx  0         nx     mx];


%Transformationsmatrix der Spannungen bei Drehung um die Y-Achse
TSig_beta=[ my^2   0 ny^2  0   2*my*ny    0;
    0      1 0     0   0          0;
    ny^2   0 my^2  0  -2*my*ny    0;
    0      0 0     my  0         -ny;
    -my*ny  0 my*ny 0   my^2-ny^2  0;
    0      0 0     ny  0          my];


%Transformationsmatrix der Spannungen bei Drehung um die Z-Achse
TSig_gamma=[mz^2   nz^2  0 0   0   2*mz*nz;
    nz^2   mz^2  0 0   0  -2*mz*nz;
    0      0     1 0   0   0;
    0      0     0 mz -nz  0;
    0      0     0 nz  mz  0;
    -mz*nz mz*nz 0 0   0   mz^2-nz^2]; 
 
TSig_Gesamt=TSig_gamma*TSig_beta*TSig_alpha;%Transformationsmatrix eines Bar i in dem VE ve
TSig_Gesamt_inv=TSig_Gesamt^-1;
TSig_Gesamt_invtra=TSig_Gesamt_inv';
Qtrans=TSig_Gesamt_inv*Q*TSig_Gesamt_invtra;

a=char(Qtrans(1,1));
b=char(Qtrans(1,2));
c=char(Qtrans(1,3));
d=char(Qtrans(1,4));
e=char(Qtrans(1,5));
f=char(Qtrans(1,6));
g=char(Qtrans(2,1));
h=char(Qtrans(2,2));
i=char(Qtrans(2,3));
j=char(Qtrans(2,4));
k=char(Qtrans(2,5));
l=char(Qtrans(2,6));
m=char(Qtrans(3,1));
n=char(Qtrans(3,2));
o=char(Qtrans(3,3));
p=char(Qtrans(3,4));
q=char(Qtrans(3,5));
r=char(Qtrans(3,6));
s=char(Qtrans(4,1));
t=char(Qtrans(4,2));
u=char(Qtrans(4,3));
v=char(Qtrans(4,4));
w=char(Qtrans(4,5));
x=char(Qtrans(4,6));
y=char(Qtrans(5,1));
z=char(Qtrans(5,2));
aa=char(Qtrans(5,3));
ab=char(Qtrans(5,4));
ac=char(Qtrans(5,5));
ad=char(Qtrans(5,6));
ae=char(Qtrans(6,1));
af=char(Qtrans(6,2));
ag=char(Qtrans(6,3));
ah=char(Qtrans(6,4));
ai=char(Qtrans(6,5));
aj=char(Qtrans(6,6));


fprintf(fid11,'%s',a);
fprintf(fid12,'%s',b);
fprintf(fid13,'%s',c);
fprintf(fid14,'%s',d);
fprintf(fid15,'%s',e);
fprintf(fid16,'%s',f);
fprintf(fid21,'%s',g);
fprintf(fid22,'%s',h);
fprintf(fid23,'%s',i);
fprintf(fid24,'%s',j);
fprintf(fid25,'%s',k);
fprintf(fid26,'%s',l);
fprintf(fid31,'%s',m);
fprintf(fid32,'%s',n);
fprintf(fid33,'%s',o);
fprintf(fid34,'%s',p);
fprintf(fid35,'%s',q);
fprintf(fid36,'%s',r);
fprintf(fid41,'%s',s);
fprintf(fid42,'%s',t);
fprintf(fid43,'%s',u);
fprintf(fid44,'%s',v);
fprintf(fid45,'%s',w);
fprintf(fid46,'%s',x);
fprintf(fid51,'%s',y);
fprintf(fid52,'%s',z);
fprintf(fid53,'%s',aa);
fprintf(fid54,'%s',ab);
fprintf(fid55,'%s',ac);
fprintf(fid56,'%s',ad);
fprintf(fid61,'%s',ae);
fprintf(fid62,'%s',af);
fprintf(fid63,'%s',ag);
fprintf(fid64,'%s',ah);
fprintf(fid65,'%s',ai);
fprintf(fid66,'%s',aj);
fclose(fid11);
fclose(fid12);
fclose(fid13);
fclose(fid14);
fclose(fid15);
fclose(fid16);
fclose(fid21);
fclose(fid22);
fclose(fid23);
fclose(fid24);
fclose(fid25);
fclose(fid26);
fclose(fid31);
fclose(fid32);
fclose(fid33);
fclose(fid34);
fclose(fid35);
fclose(fid36);
fclose(fid41);
fclose(fid42);
fclose(fid43);
fclose(fid44);
fclose(fid45);
fclose(fid46);
fclose(fid51);
fclose(fid52);
fclose(fid53);
fclose(fid54);
fclose(fid55);
fclose(fid56);
fclose(fid61);
fclose(fid62);
fclose(fid63);
fclose(fid64);
fclose(fid65);
fclose(fid66);