Discussion :

[duhill]

Bonsoir,

Je sollicite votre aide pour le problème suivant :

Je cherche à assembler une matrice k, n fois dans une matrice initiale "a" composée de 0.
Les matrices k sont assemblées de façon décalée en ligne et colonne lorsque l'on incrémente la boucle.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
//paramètres de la matrice à assembler
l=1 
k=[12*l 6*l -12 6*l ; 6*l 4*(l^2) -6 2 ; -12 -6 12 -6*l ; 6 -2 -6 4];

//boucle pour assembler la matrice k
for j=1:3
    n=3
    
    a=zeros(2*(n+1),2*(n+1))

    a(j:j+3,j:j+3)=a(j:j+3,j:j+3)+k

//la réponse désirée de la boucle est la suivante
//a(1:4,1:4)=a(1:4,1:4)+k
//a(3:6,3:6)=a(3:6,3:6)+k
//a(5:8,5:8)=a(5:8,5:8)+k 
end

disp(a)

La réponse est la suivante:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
    0.    0.    0.     0.    0.     0.    0.    0.  
    0.    0.    0.     0.    0.     0.    0.    0.  
    0.    0.    12.    6.  - 12.    6.    0.    0.  
    0.    0.    6.     4.  - 6.     2.    0.    0.  
    0.    0.  - 12.  - 6.    12.  - 6.    0.    0.  
    0.    0.    6.   - 2.  - 6.     4.    0.    0.  
    0.    0.    0.     0.    0.     0.    0.    0.  
    0.    0.    0.     0.    0.     0.    0.    0.

Celle que je désirerais est celle ci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
    12.    6.  - 12.    6.    0.     0.    0.     0.  
    6.     4.  - 6.     2.    0.     0.    0.     0.  
  - 12.  - 6.    24.    0.  - 12.    6.    0.     0.  
    6.   - 2.    0.     8.  - 6.     2.    0.     0.  
    0.     0.  - 12.  - 6.    24.    0.  - 12.    6.  
    0.     0.    6.   - 2.    0.     8.  - 6.     2.  
    0.     0.    0.     0.  - 12.  - 6.    12.  - 6.  
    0.     0.    0.     0.    6.   - 2.  - 6.     4.

Merci de votre aide.

[Jerome Briot]

Sort les lignes suivantes de la boucle :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
n = 3
a = zeros(2*(n+1),2*(n+1))

[duhill]

Je te remercie pour ta réponse.
Le fait de sortir ces variables de la boucle permet d'obtenir un résultat différent, étonnement, mais le résultat n'est toujours pas bon.
J'ai donc bidouillé la boucle et ai ajouté une sous matrice, avec les termes sommés aux bon endroits, pour résoudre mon problème.
Je poste ma solution en espérant que sa aidera quelqu'un.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
//paramètres de la matrice à assembler
l=1 
k=[12*l 6*l -12 6*l ; 6*l 4*(l^2) -6 2 ; -12 -6 12 -6*l ; 6 -2 -6 4];

ke=[24*l 0 -12 6*l ; 0 8*(l^2) -6 2 ; -12 -6 24 -0 ; 6 -2 0 8];//matrice avec les termes sommés

n=4
a=zeros(2*(n+1),2*(n+1))
a(1:4,1:4)=k
//boucle pour assembler la matrice k
for j=3:2:(3*(n-1))
    a(j:(j+3),j:(j+3))=ke    

end
disp(k)

disp(a)

[Jerome Briot]

le résultat n'est toujours pas bon.

Car les indices doivent dépendre de j et de n :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
for j = 1:3
 
    idx = (j-1)*(n-1)+1:(j-1)*(n-1)+4
    disp(idx)
 
end

Et dans le cas où la matrice k semble fortement être la matrice de raideur de l'élément poutre en flexion 2D, je te conseille de bien vérifier les coefficients.