RDforum.free.fr/Cast3M   Index du Forum
Site Web Cast3M
RDforum.free.fr/Cast3M
Forum d'entraide des utilisateurs de Cast3M
Répondre au sujet
Problème élastique simple convergence avec PASAPAS
Vincent G.


Inscrit le: 21 Oct 2009
Messages: 1
Répondre en citant
Bonjour,

J'ai un problème de non-convergence avec PASAPAS sur un modèle élastique extrèmement simple (code ci dessous).
Il s'agit ici tout simplement d'un bloc élastique en compression, le cas que je veux traiter au final est évidemment plus complexe (géométrie différente, viscosité et chargement variable au cours du temps ...), mais vu que j'étais confronté à ces problèmes de convergence je l'ai simplifié progressivement pour comprendre ce qui se passait et je suis finalement arrivé à celui-ci...
Je me doute bien qu'il doit s'agir d'une erreur évidente mais je n'arrive pas à mettre le doigt dessus.

Merci d'avance,
Vincent


Code:

OPTI DIME 2 ELEM tri3;

*maillage
DENSITE 0.02;
P1  =    0.   0.;
P2  =    1.   0.;
P3  =    1.   1.;
P4  =    0.   1.;

P1P2 = P1 DROI P2;
P2P3 = P2 DROI P3;
P3P4 = P3 DROI P4;
P4P1 = P4 DROI P1;

C1=P1P2 ET P2P3 ET P3P4 ET P4P1;
S1= SURF C1;

*modele, materiaux
MO1=MODE S1 MECANIQUE 'ELASTIQUE' 'ISOTROPE';
MAT1='MATERIAU' MO1 YOUN 20.e9 NU 0.25 RHO 2900.;

*temps de calcul
N_cal=10;
D_cal=200;
dt_cal=D_cal/N_cal;         
LX= PROG 0. PAS dt_cal NPAS  N_cal;
LY=PROG 'LINE' 'A' (1./D_cal) 'B' 0. LX;

*chargement
F1=MANU 'CHPO' P3P4  2 FX 0. FY -10. NATURE DIFFUS;
ev='EVOL' 'MANUEL' 'temps' LX 'force' LY;
CHA1=CHAR 'MECA' F1 ev;

*rigidite
R1=RIGI MO1 MAT1;
Cl1=BLOQ P1P2  'UY';
Cl2=BLOQ P2P3  'UX';
Cl3=BLOQ P4P1  'UX';
RIG1=R1 ET Cl1 ET Cl2 ET Cl3;


*lancement pasapas
TAB1=TABLE;
TAB1.'BLOCAGES_MECANIQUES'=RIG1;
TAB1.'MODELE'= MO1;
TAB1.'CHARGEMENT'=CHA1;
TAB1.'CARACTERISTIQUES'= MAT1;
TAB1.'TEMPS_CALCULES'=LX;
PASAPAS TAB1;

Voir le profil de l'utilisateurEnvoyer un message privé
Problème élastique simple convergence avec PASAPAS
Vous ne pouvez pas poster de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas voter dans les sondages de ce forum
Vous pouvez joindre des fichiers
Vous pouvez télécharger des fichiers
Toutes les heures sont au format GMT + 2 Heures  
Page 1 sur 1  

  
  
 Répondre au sujet