newtazzmania
Inscrit le: 19 Mar 2009 |
Messages: 6 |
|
|
|
Posté le: 19 Mar , 2009 17:19 |
|
|
|
|
|
Bonjour, je rencontre un problème lors de la lecture de données d'un fichier .sauv.
Ce fichier contient les résultats d'un calcul avec la procédure STATIO, ainsi, je cherche à recupérer les contraintes concernant chaque noeud de mon maillage.
Pour cela, je commence par calculer le nombre de noeud dans mon maillage MAISTA, et ensuite, je boucle de 1 à NbNoeud, cependant, je fais une erreur en faisant cela car je boucle sur les noeuds de 1 à NbNoeud mais en fait, les NbNoeuds n'appartiennent pas forcement à ce maillage. J'ai essayé de passé par la procedure POIN qui me dit que mon OBJET est TROP COMPLEX...
Je me suis alors dit qu'il fallait que je recupère tous les noeuds en fait il y en à 39928 en tout (alors que NbNoeud = 26244), mais je ne peux voir ce nombre uniquement quand je restaure mon fichier .sauv à la ligne :
LECTURE DE 39928 OBJETS DE TYPE POINT.
Je n'ai pas trouvé comment dans TAB1, il est stocké, ni comment il est calculé.
NbNoeud = NBNO MAISTA;
impr0 = VALE IMPR ;
echo0 = VALE ECHO ;
OPTI IMPR 'res_statio.txt' ECHO 0;
sigt = TAB1 . 'CONTRAINTES' . 1 ;
XX = chan 'CHPO' MODSTA (exco 'SMXX' sigt) ;
YY = chan 'CHPO' MODSTA (exco 'SMYY' sigt) ;
ZZ = chan 'CHPO' MODSTA (exco 'SMZZ' sigt) ;
XY = chan 'CHPO' MODSTA (exco 'SMXY' sigt) ;
XZ = chan 'CHPO' MODSTA (exco 'SMXZ' sigt) ;
YZ = chan 'CHPO' MODSTA (exco 'SMYZ' sigt) ;
repe BOUCLE NbNoeud;
i=&BOUCLE;
P_i = NOEUD i;
Sxx = extr XX 'SMXX' P_i ;
Syy = extr YY 'SMYY' P_i ;
Szz = extr ZZ 'SMZZ' P_i ;
Sxy = extr XY 'SMXY' P_i ;
Sxz = extr XZ 'SMXZ' P_i ;
Syz = extr YZ 'SMYZ' P_i ;
MESS Sxx Syy Szz Sxy Sxz Syz ; ***recuperation des contraintes
MESS (COOR 1 P_i) (COOR 2 P_i) (COOR 3 P_i); ***recuperation des coordonnées de mon noeud
FIN BOUCLE;
|
|
newtazzmania
Inscrit le: 19 Mar 2009 |
Messages: 6 |
|
|
|
Posté le: 24 Mar , 2009 17:26 |
|
|
|
|
|
Voila j'ai trouvé une solution à mon problème, un petit tratement supplémentaire dans matlab et c'est bon.
S'il y en a que ça interresse faites moi signe!!!
|
|