Discussion :

[anita.crus]

super, ça marche merci

Donc la pour calculer les modes propres, j'ai trouvé ce code :
Mise en œuvre sous MATLAB :

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
close all
 
% matrice de masse
matm = diag([m1 m2 m3]) ;
% matrice de raideur
matk = [k1+k2 -k2    0 ; ...
        -k2  k2+k3 -k3 ; ...
         0   -k3    k3 ] ;
% RECHERCHE DES MODES PROPRES
[V,D]=eig(matk,matm) ;
% fréquences propres
fqce = sqrt(diag(D))/(2*pi)
% vecteurs propres normés en composante maximale unitaire
V1 = V(:,1)/max(abs(V(:,1)))*sign(V(find(abs(V(:,1))==max(abs(V(:,1)))),1))
V2 = V(:,2)/max(abs(V(:,2)))*sign(V(find(abs(V(:,2))==max(abs(V(:,2)))),2))
V3 = V(:,3)/max(abs(V(:,3)))*sign(V(find(abs(V(:,3))==max(abs(V(:,3)))),3))

La question est !
Est-ce que MATLAB pourra calculer les matrices de rigidités élémentaires car la j'ai plus de 100 barres et ça va être compliqué et long de les calculer à la main. croyez vous qu'il existe un moyen??

[Jerome Briot]

Les matrices élémentaires ne dépendent que de la position des nœuds, de l'orientation et de la section des barres, ainsi que du matériau.

Tu as donc tout ce qu'il te faut pour faire une boucle qui va te créer les matrices locales et l'assemble de la matrice globale.

[anita.crus]

Est ce que tu peux m'expliquer le code que t'avais fait. Certes que ça marche super bien mais il y a quelques lignes que je ne comprends pas.

[Invité]

il y a quelques lignes que je ne comprends pas.

Et... pourrait-on savoir lesquelles ?

[anita.crus]

Oui bien sur:

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
20
21
22
23
% Recuperation des coordonnees des noeuds
idx = ~cellfun('isempty', strfind(X{1},'geompy.MakeVertex'));
subX = X{1}(idx);
 
points = struct('coord', zeros(numel(subX),3), 'name', {});
 
lignes = struct('point1',{}, 'point2',{}, 'name', {});
 
for n = 1:numel(subX)
    lignes(1).point1{n} = sscanf(subX{n}, '%*s = geompy.MakeLineTwoPnt(%s%*s');
    lignes(1).point2{n} = sscanf(subX{n}, '%*s%*s%*s%s');
    lignes(1).name{n} = sscanf(subX{n}, '%s%*s%*s%*s');
end
 
lignes.point1 = cellfun(@(x) x(1:end-1), lignes.point1, 'UniformOutput',false);
 
 
 
for n = 1:numel(lignes.point1)
 
    idx = strcmp(points.name, lignes.point1{n}) | strcmp(points.name, lignes.point2{n});
 
    if strncmp(lignes.name{n}, 'BARR', 4)

[Invité]

Je te suggère la lecture de la FAQ Tableaux de cellules, et plus particulièrement Comment appliquer une fonction à chaque cellule ?

[anita.crus]

Je te suggère la lecture de la FAQ Tableaux de cellules, et plus particulièrement Comment appliquer une fonction à chaque cellule ?

X=textscan(fid, '%s', 'delimiter', '\n') : lire le fichier fid, sépare les chaines de caractère et le mettre dans une seule cellule.
Est ce que c'est correct??

[Jerome Briot]

Non pas vraiment, dans ce cas, textscan stocke chaque ligne du fichier dans une cellule du tableau de cellules X

C'est très utile pour lire un fichier en une seule fois quand son contenu change à chaque ligne.

Sinon, il faudrait faire :

Code :

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

1
2
3
4
5
6
7
8
9
fid = fopen(..., 'rt')
 
while 1
 
   str = fgetl(fid);
 
end
 
fclose(fid);

[anita.crus]

Merci beaucpu !!
donc

subX=X{1}(idx) :: stocke les valeurs vrai de X{1} ??

[Jerome Briot]

Oui, selon les valeurs prises par la variable idx pour la condition qui précède.
Voir la FAQ : Qu'est-ce que l'indexage logique ?

[anita.crus]

Je te suggère la lecture de la FAQ Tableaux de cellules, et plus particulièrement Comment appliquer une fonction à chaque cellule ?

J'ai 3 groupes d'éléments avec 3 paramètres différents. Comment je pourrais attribuer ces paramètres à ces groupes de barres.

[Invité]

Pourrais-tu expliquer ton problème plus en détails ? Je ne sais absolument pas à quoi correspond chaque élément mentionné.

[anita.crus]

En fait mes éléments sont des barres ; et j'ai 3 groupes de barres qui ont des modules de young et des sections différentes.

[Invité]

"barre" = graphique ?
Tu aurais un schéma à présenter, que ce soit plus clair ?

[anita.crus]

"barre" = graphique ?
Tu aurais un schéma à présenter, que ce soit plus clair ?

Oui !! en fait ma structure est celle la .module.py

[anita.crus]

et j'ai utilisé ce code pour la générer sur MATLAB :

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
% fid = fopen('module.py', 'r');
fid = fopen('treilli.py', 'r');
X = textscan(fid, '%s', 'delimiter', '\n');
fclose(fid);
 
% Recuperation des coordonnees des noeuds
idx = ~cellfun('isempty', strfind(X{1},'geompy.MakeVertex'));
subX = X{1}(idx);
 
points = struct('coord', zeros(numel(subX),3), 'name', {});
 
for n = 1:numel(subX)
    points(1).coord(n,:) = sscanf(subX{n}, '%*s = geompy.MakeVertex(%f,%f,%f)');
    points(1).name{n} = sscanf(subX{n}, '%s%*s%*s%*s%*s');
end
 
% Recuperation de la connectivite des barres
idx = ~cellfun('isempty', strfind(X{1},'geompy.MakeLineTwoPnt'));
subX = X{1}(idx);
 
lignes = struct('point1',{}, 'point2',{}, 'name', {});
 
for n = 1:numel(subX)
    lignes(1).point1{n} = sscanf(subX{n}, '%*s = geompy.MakeLineTwoPnt(%s%*s');
    lignes(1).point2{n} = sscanf(subX{n}, '%*s%*s%*s%s');
    lignes(1).name{n} = sscanf(subX{n}, '%s%*s%*s%*s');
end
 
lignes.point1 = cellfun(@(x) x(1:end-1), lignes.point1, 'UniformOutput',false);
lignes.point2 = cellfun(@(x) x(1:end-1), lignes.point2, 'UniformOutput',false);
 
figure
plot3(points.coord(:,1), points.coord(:,2), points.coord(:,3), 'ro', 'markerfacecolor', 'r');
hold on
 
for n = 1:numel(lignes.point1)
 
    idx = strcmp(points.name, lignes.point1{n}) | strcmp(points.name, lignes.point2{n});
 
    if strncmp(lignes.name{n}, 'BARR', 4)
        h_barr(n) = line(points.coord(idx,1), points.coord(idx,2), points.coord(idx,3) ...
            , 'color', 'g', 'linewidth', 2);
    else
        h_geomobj(n) = line(points.coord(idx,1), points.coord(idx,2), points.coord(idx,3) ...
            , 'color', 'b', 'linewidth', 2);
    end
 
end
axis equal vis3d

Le code marche super bien et c'est DUT qui me l'a filé

[anita.crus]

Donc je voudrais affecter des paramètres à mes barres pour pouvoir calculer la matrice de rigidité.

[Jerome Briot]

Vu que tu débutes, je ne peux que te conseiller la lecture de cette ouvrage : MATLAB Guide to Finite Elements: An Interactive Approach par Peter Kattan

[anita.crus]

Merci pour le lien, je vais le faire. Mais ce que j'ai à résoudre en ce moment c'est très urgent.c'est pour ça que je cherche les réponses un peu partout.

[anita.crus]

J'ai effectué ce code pour un élément barre dans un treillis :

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 
%%%% exemple pour une barre %%
 
% Calcul de la matrice raideur Ke 
% pour un élément (iel) d'une structure treillis
%
% 
% en entrée iel : numéro de l'élément
% en sortie Ke  : matrice raideur élémentaire 
 
 
function [Ke] = barre_ke(iel)
 
 
Nprop=[1;1];              % pour chaque élément N° de la propriété
Prop=[ 2.e10 0.5];        % tableau des différentes valeurs de E et S 
 
E=Prop(Nprop(iel),1);
S=Prop(Nprop(iel),2);
X  = Coord(Connec(iel,:),:);    % ----- Coordonnées des 2 noeuds de l'élément
dX = X(2,:) - X(1,:);           % ----- x2-x1  et y2-y1
if ndim == 1
    L = abs(dX);
    Ke = (ES/L)*[  1  -1
                  -1   1 ];
elseif ndim == 2
    L = sqrt(dX(1)^2 + dX(2)^2);
    c = dX(1)/L; s = dX(2)/L;   % ----- Cosinus directeurs de l'élément
    cc = c*c; ss = s*s; cs = c*s;
    Ke = (E*S/L)*[  cc  cs -cc -cs
                   cs  ss -cs -ss
                  -cc -cs  cc  cs
                  -cs -ss  cs  ss];
 end
return