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
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;
|