probleme matlab Résolu/Fermé

Question

Bonjour, je dois créer un programme matlab afin de déterminer la matrice de rigidité d'après la méthode des éléments fini.

Je dois créer une matrice C, composé de la façon suivante:

N est le nombres de barres 

   u = eye(2*N);             
for i = 1:1:N
      for j = 1:1:2
         u(e(i,j)*2-1,:)
         u(e(i,j)*2,:)

      end
end
 
Le problème c'est que je reçois la réponse sous forme de ligne, alors que je souhaiterai pouvoir mettre les réponses sous la forme d'une matrice de la taille:

m(1:4;2*N)

Merci de votre aide.

Sacabouffe · Answer

Salut
Ton ptit bout de code veut rien dire du tout... Du moins tel quel sans les détails, ça veut rien dire, donc si tu pouvais détailler un peu on pourrait peut-être t'aider.
A plus

jeremy124 · Answer

bonjour, je copie donc mon code/

clear all
clc
[n, e,A,Emod,BC] = bar_input; %le fichier [n, e,A,Emod,BC] comprend des données sous formes de matrices
[N,M] = size(n);

%lenght
for i = 1:1:N
l(i) = sqrt((n(e(i,2),1)-n(e(i,1),1))^2+(n(e(i,2),2)-n(e(i,1),2))^2)
end

%angle
for i = 1:1:N
if ((n(e(i,2),1)-n(e(i,1),1))==0)&&((n(e(i,2),2)-n(e(i,1),2))>0)
alp(i)=pi()/2;
elseif ((n(e(i,2),1)-n(e(i,1),1))==0)&&((n(e(i,2),2)-n(e(i,1),2))<0)
alp(i)=3*pi()/2;
elseif ((n(e(i,2),1)-n(e(i,1),1))>0)&&((n(e(i,2),2)-n(e(i,1),2))==0)
alp(i)=0;
elseif ((n(e(i,2),1)-n(e(i,1),1))<0)&&((n(e(i,2),2)-n(e(i,1),2))==0)
alp(i)=pi();
elseif ((n(e(i,2),1)-n(e(i,1),1))>0)&&((n(e(i,2),2)-n(e(i,1),2))>0)
alp(i)=atan((n(e(i,2),2)-n(e(i,1),2))/((n(e(i,2),1)-n(e(i,1),1))));
elseif ((n(e(i,2),1)-n(e(i,1),1))>0)&&((n(e(i,2),2)-n(e(i,1),2))<0)
alp(i)=pi()+atan((n(e(i,2),2)-n(e(i,1),2))/((n(e(i,2),1)-n(e(i,1),1))));
end
end

%matrice k
for i = 1:1:N

K(1:4,1:4,i)=(Emod(i)*A(i)/l(i))*[cos(alp(i))*cos(alp(i)) cos(alp(i))*sin(alp(i)) -cos(alp(i))*cos(alp(i)) -cos(alp(i))*sin(alp(i)); cos(alp(i))*sin(alp(i)) sin(alp(i))*sin(alp(i)) -cos(alp(i))*sin(alp(i)) -sin(alp(i))*sin(alp(i)); -cos(alp(i))*cos(alp(i)) -cos(alp(i))*sin(alp(i)) cos(alp(i))*cos(alp(i)) cos(alp(i))*sin(alp(i)); -cos(alp(i))*sin(alp(i)) -sin(alp(i))*sin(alp(i)) cos(alp(i))*sin(alp(i)) sin(alp(i))*sin(alp(i))];

end

%matrice C
u = eye(2*N); % we create the identity matrix
for i = 1:1:N
for j = 1:1:2
u(e(i,j)*2-1,:)
u(e(i,j)*2,:)

end
end

la reponse de la matrice C est de la forme /

ans =

1 0 0 0 0 0

ans =

0 1 0 0 0 0

ans =

0 0 1 0 0 0

ans =

0 0 0 1 0 0
/merci de votre aide/
bonne journée

Sacabouffe · Answer

Salut
C'était juste ça?
Donc si j'ai bien compris, effectivement tout le code servait pas à grand chose pour te répondre... désolé... :-D
Tu veux juste mettre les lignes des ans dans la matrice C?
Par contre il y a un truc que je comprends pas, t'as pas écrit toutes les lignes de la sortie écran dans ton message? Parce qu'apparemment ton N=3, donc il devrait y avoir 12 lignes à la sortie, non? Et dans ce cas pourquoi le m(1:4;2*N) de ton premier message?

Donc... si effectivement je suis pas à côté de la plaque, personnellement j'écrirais ça :
% matrice C
C = [];
u = eye(2*N); % we create the identity matrix
for i = 1:1:N
for j = 1:1:2
C=[C;u(e(i,j)*2-1,:);u(e(i,j)*2,:)]
end
end 
Ou sous forme matricielle 
% matrice C
u = eye(2*N); % we create the identity matrix
C=u(reshape([reshape((e.^2).'-1,1,numel(e));reshape((e.^2).',1,numel(e))],2*numel(e),1),:);
Dis moi si je suis carrément à côté parce qu'en fait je vois pas pourquoi il y a que 4 lignes à ta sortie écran...
A plus

jeremy124 · Answer

Bonjour, 
Je te remercie de ta patience, comme tu l as remarqué je n'ai pas copié toutes lignes de sortie d'écran, 
j ai bien 12 lignes, le problème c'est que je souhaite avoir en sortie des matrices de la forme  m(1:4;2*N), 
Soit 3 matrices pour ce cas là, de la taille m(4,6)
En gros une matrice par barres c'est à dire par éléments étudiés.

Merci de ton aide

jeremy124 · Answer

le problème c'est que cela doit être general, est s appliquer pour tout les cas possibles, par exepemple si j'ai 100 éléments, je doi avoir 100 matrice différentes et c'est cela qui me pose problème, je n'arrive pas à créer une fonction pour avoir N matrice.

Sacabouffe · Answer

Ben pour ça je pense pas qu'il existe de solution simple.
En fait c'est un ensemble de matrices paramétrées par un entier p &#1028; [1,N] que tu veux donc ce que tu peux faire c'est définir un tableau à 3 dimensions par exemple ou utiliser des cellules.
En partant de la matrice C (de taille 4*N x 2*N) du message 3 ça peut donner un truc du genre :

-Avec un tableau à 3 dimensions : C=permute(reshape(C.',2*N,4,N),[2 1 3]); Pour tout p &#1028; [1,N], C(:,:,p) sera la p-ième matrice.

-Avec des cellules : C=mat2cell(C,4*ones(1,N),2*N); Pour tout p &#1028; [1,N], C{p} sera la p-ième matrice.

Ça te convient ce petit bidouillage?

jeremy124 · Answer

Je te remercie, c’est exactement ce que je recherché….
Merci de ta patience et de tes explications qui mon permis de comprendre.

jeremy124 · Answer

bonjour, je reviens toujours a la recherche d aide.
j ai un probleme pour reduire une matrice.
j ai une matrice de taille k= [2*N,2*N]
une autre,ma matrice de condition aux limites de la forme BC=[2,N],elle est reöpli de 1 et de 0

exemple

% Boundary Conditions
%     X1 Y1  X2 Y2 X3 Y3
BC = [0  0 ;  0   0;  1   1 ]

%matrice k
k =

  1.0e+007 *

    0.7425   -0.7425   -0.0000   -0.0000   -0.7425    0.7425
   -0.7425    2.8425   -0.0000   -2.1000    0.7425   -0.7425
   -0.0000   -0.0000    2.1000    0.0000   -2.1000         0
   -0.0000   -2.1000    0.0000    2.1000         0         0
   -0.7425    0.7425   -2.1000         0    2.8425   -0.7425
    0.7425   -0.7425         0         0   -0.7425    0.7425

je souhaiterai pouvoir reduire la matrice de sorte qu il ne me reste que les valeurs ou les conditions limites ns sont pas nul, donc ds cette exemple avoir seulement une matrice 2*2 de cette forme 

 2.8425   -0.7425
 -0.7425    0.7425

merci de votre aide

Sacabouffe · Answer

Salut jeremy124
Je jette un œil ce we ou ce soir à ton message.
A plus tard

Sacabouffe · Answer

Re salut jeremy124
Pas facile de te suivre des fois :-D
Faudrait que tu détailles un peu plus (et telle qu'elle est écrite, ta matrice BC a plutôt une taille Nx2)
Mais bon... Si j'ai bien compris, ce que tu demandes est un truc dans ce genre :
k=k(find(reshape(BC.',1,numel(BC))),find(reshape(BC.',1,numel(BC))));
Si c'est pas ça c'est que j'ai pas compris. Dans ce cas détaille un peu plus... Comment tu veux sélectionner les éléments de k en fonction de BC?
Tiens moi au courant...
A plus

jeremy124 · Answer

Merci beaucoup, c’est exactement ce que je souhaitai.
Je te remercie encore une fois pour ton aide précieuse.
Surement à bientôt pour d’autres questions ….
A plus

jeremy124 · Answer

bonjour

J´ai de nouveau un probleme, je copie mon code pr plus de precision:

f = [f1x f1y; f2x f2y; f3x f3y];

z=reshape((flipud(rot90(f))),1,2*N); % passage en vecteur de la matrice des forces
inconnues=z(find(reshape(BC.',1,numel(BC)))) %prend les valeurs ds C suivant les "1" ds BC
      
f=inv(kfinal)*inconnues'

j aimerai pouvoir former un nouveau vecteur avec les valeurs de f et des zero pour le reste (cad la ou BC = 0 )

jeremy124 · Answer

je donne un exemple si  BC = [0 0 ; 1 0; 1 1 ]  alors  f=   -0.0524
                                                                                -0.0524                                                                                                                   
                                                                                -0.2529

 et je souhaiterai avoir un vecteur de la forme vecteur=[ o,o,-0.0524,o,-0.0524, -0.2529]

merci d avance

Sacabouffe · Answer

Salut l'ami
Je pense qu'un truc du genre devrait faire l'affaire :

fext=zeros(1,numel(BC));
fext(logical(reshape(BC.',1,numel(BC))))=f;

Dis-moi si c'est ça.
A plus

jeremy124 · Answer

salut

merci encore pour ton aide, j ai de nouveau un probleme, ca n arrete jamais.....


for i = 1:1:N 
kt(:,:,i)=(Emod(i)*A(i)/l(i))*[C(:,:,i)'*K(1:4,1:4,i)*C(:,:,i)];
end 

ktotal=sum(kt,N); % we create summary matrix kt 

kfinal=ktotal(find(reshape(BC.',1,numel(BC))),find(reshape(BC.',1,numel(BC)))); %matrice k en fonction des BC

z=reshape((flipud(rot90(f))),1,2*N); % passage en vecteur de la matrice des forces

inconnues=z(find(reshape(BC.',1,2*N))) %prend les valeurs ds C suivant les "1" ds BC


force=inv(kfinal)*inconnues'

et pour les condition suivante 
 
A =400; % mm^2
E =210000; % MPa
% 
f1x = 0;
f1y = 0;
f2x = 90;
f2y = 100;
f3x = 10;
f3y = 110;

% Nodes
n =[2 6; 2 2; 6 2];

% Elements
e =[1 2; 2 3; 1 3];

% Areas
A =[A; A; A];

% Modulus of Elasticity
Emod =[E; E; E];

% Forces
f = [f1x f1y; f2x f2y; f3x f3y];

% Boundary Conditions
BC =[0 0 ; 1 1; 1 1];

ca me donne cela

inconnues =

    90   100    10   110

Warning: Matrix is singular to working precision.
> In solveode at 63

force =

   Inf
   Inf
   Inf
   Inf


displacement =

     0
     0
   Inf
   Inf
   Inf
   Inf


je ne comprend pas mon erreur, donc voila mon probleme.....
merci de ton aide

jeremy124 · Answer

finalement j ai resolu mom probleme 
merci de ton aide
mon programme semble fini
bon WK

jeremy124 · Answer

j ai de nouveau un probleme, j'ai un message d erreur qui s'affiche

??? Error using ==> reshape
To RESHAPE the number of elements must not change.

Error in ==> solveode at 33
C=permute(reshape(C.',2*N1,4,N1),[2 1 3]); % we create the different matrix c

mon code est ceci:

C = [];
u = eye(2*N); % we create the identity matrix

for i = 1:1:N1
    for j = 1:1:2
  C=[C;u(e(i,j)*2-1,:);u(e(i,j)*2,:)];
    end
end 
C=permute(reshape(C.',2*N1,4,N1),[2 1 N1]);

la taille de la matrice e change suivant les problemes d'elements finis a resoudre.
merci de votre aide

bibou12 · Answer

salut 
je suis débutante en matlab et je doit faire une connection avec une base de donnée
si vous pouvez m'aider pour la syntaxe
merci d'avance

jeremy124 · Answer

hej 
alors j'ai de nouveau un probleme, je souhaiterai pouvoir afficher une matrice dans une fonctions sprintf, cela est il possible?
 merci d avance
jeremy

medi2 · Answer

insupportable ce forum ,j'ai envoye un message deux fois ei il a ete effacé !!
je peux comprendre pourquoi.

Probleme matlab

20 réponses

Discussions similaires