Matlab

[invite51e65278]

Bonjour! J'ai urgemment besoin de votre aide

Je voudrais associer a chaque Klocal(pour les differentes valeurs de k) une matrice A indice K.

Donc pour K=1 j'aurai Klocal indice 1

Dans le programme suivant

for k=1: p

B(k,1)=input ('Entrer le numero du noeud de depart de la barre suivante:');
B(k,2)=input ('Entrer le numero du noeud d arrive de la barre suivante:');

E=input ('Entrer le module de Young (Kpa) de la barre concerne:');
S=input ('Entrer la section(m^2) de la barre concernee:');
I=input ('Entrer linertie(m^4)de la barre concernee:');

L=sqrt((A(B(k,2),1)- A(B(k,1),1))*(A(B(k,2),1)- A(B(k,1),1))+(A(B(k,2),2)- A(B(k,1),2))*(A(B(k,2),2)- A(B(k,1),2)));

klocal (1,1)= (E*S)/L;
klocal (1,2)= 0;
klocal (1,3)= 0;
klocal (1,4)= -(E*S)/L;
klocal (1,5)= 0;
klocal (1,6)= 0;
klocal (2,1)= klocal (1,2);
klocal (2,2)= (12*E*I)/(L*L*L);
klocal (2,3)= (6*E*I)/(L*L);
klocal (2,4)= 0;
klocal (2,5)= -(12*E*I)/(L^3);
klocal (2,6)= (6*E*I)/(L^2);
klocal (3,1)= klocal (1,3);
klocal (3,2)= klocal (2,3);
klocal (3,3)= (4*E*I)/L;
klocal (3,4)= 0;
klocal (3,5)= -(6*E*I)/(L^2);
klocal (3,6)= (2*E*I)/L;
klocal (4,1)= klocal (1,4) ;
klocal (4,2)= klocal (2,4);
klocal (4,3)= klocal (3,4);
klocal (4,4)= (E*S)/L;
klocal (4,5)= 0;
klocal (4,6)= 0;
klocal (5,1)= klocal (1,5);
klocal (5,2)= klocal (2,5);
klocal (5,3)= klocal (3,5);
klocal (5,4)= klocal (4,5);
klocal (5,5)= (12*E*I)/(L^3);
klocal (5,6)= -(6*E*I)/(L^2);
klocal (6,1)= klocal (1,6);
klocal (6,2)= klocal (2,6);
klocal (6,3)= klocal (3,6);
klocal (6,4)= klocal (4,6);
klocal (6,5)= klocal (5,6);
klocal (6,6)= (4*E*I)/L;

Je veux donc stocker tous les k local des differentes barres.Comment faire?
-----

[invite0ab6bca5]

Bonjour,
J'ai cru comprendre que vous vouliez stocker un ensemble de matrices klocal_n.
1) se méfier de tous ces "k" pas très différents des "K" et des erreurs typographiques posibles ;
2) je vois 2 possibilités pour le stockage de cette collection de matrices :
a) la solution acrobatique : utiliser un tableau tridimensionnel (voir par exemple l'aide concernant "meshgrid") ;
b) la solution classique : faire une grosse matrice par concaténation des klocal successives :
GRMAT = [] ;
for ind = 1 : indmax,
....fabrication du klocal n° ind ....
GRMAT = [GRMAT klocal]
end

On récupère la nième matrice klocal en prenant, dans GRMAT, les colonnes de (n-1)*(nb de colonnes de klocal) +1 à
(n-1)*(nb de colonnes de klocal) + (nb de colonnes de klocal).

Bonne continuation