hola chicos, abro este hilo viejo porq estoy tratando de hacer algo parecido.
tengo varios bucles uno dentro del otro. el de afuera varia de (e) 1:3 y los dos de adentro (i,j) 1:2. por cada vez q se completa el bucle de afuera (e) se me forma una matriz de 8x8. por lo tanto voy a terminar teniendo 3 matrices de 8x8.
como hago para q me guarde esas tres matices? la parte del codigo q trata esto es esta:
k0=zeros(8,8);
pg=[1 1.2 2.2 1.8 ; -0.1 -0.5 0.0 -0.3];
r=[00 0.1]; s=[-2.5 -3.1];
for e=1:3
for i=1:2
for j=1:2
dh1r=1/4*(1+s(j)); dh1s=1/4*(1+r(i));
dh2r=-1/4*(1+s(j)); dh2s=1/4*(1-r(i));
dh3r=-1/4*(1-s(j)); dh3s=-1/4*(1-r(i));
dh4r=1/4*(1-s(j)); dh4s=-1/4*(1+r(i));
dh=[dh1r dh2r dh3r dh4r ; dh1s dh2s dh3s dh4s];
jacov=dh*pg';
invj=inv(jacov);
b1=[invj(1,1) invj(1,2) 0 0;
0 0 invj(2,1) invj(2,2);
invj(2,1) invj(2,2) invj(1,1) invj(1,2)];
b2=[dh1r 0 dh2r 0 dh3r 0 dh4r 0;
dh1s 0 dh2s 0 dh3s 0 dh4s 0;
0 dh1r 0 dh2r 0 dh3r 0 dh4r;
0 dh1s 0 dh2s 0 dh3s 0 dh4s];
b=b1*b2;
end
end
k(e)=k0+b'*b*det(jacov)
end
la idea es terminar con 3 matrices k(e) q son de 8x8, pero no logro encontrar la manera de hacerlo.
alguien sabe como?
saludos!
