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
écrouissage mode pas a pas
sam


Inscrit le: 29 Jan 2010
Messages: 3
Répondre en citant
bonjour j ai un petit probleme :
je simule une pression variable sur une plaque trouée
j applique une pression de 100 puis de 0 puis de 100 etc ....
qd je trace la déformée au moment ou j applique la pression nulle je devrais retrouver ma plaque trouée d origine ms la je la retrouve avec qques endroits du maillage qui sont incohérents, comme vous pouvez le voir en executant le prgme. Qqn saurait il m aider s il vous plait .
merci d avance

mon premier fichier que j appelle ecrou.txt
opti dime 2 echo 1 elem QUA8;

lplaque = 10.; r = 1.;
p1 = 0. 0.;
p2 = lplaque 0.;
p3 = lplaque lplaque;
p4 = 0. lplaque;
a2 = r 0.;
a4 = 0. r;
a3 = (r/1.4142) (r/1.4142);
nelem = 15; nel_bis = -25;
d12 = droite nel_bis a2 p2
dini 0.25 dfin 1.5 ;
d23 = droite nelem p2 p3;
d34 = droite nelem p3 p4;
d41 = droite nel_bis p4 a4
dini 1.5 dfin 0.25 ;
bis = droite nel_bis a3 p3
dini 0.3 dfin 2.2 ;
arc23 = cerc nelem a2 p1 a3;
arc34 = cerc nelem a3 p1 a4;
dom1 = dall
d12 d23 (inve bis) (inve arc23);
dom2 = dall (inve d41) (inve d34)
(inve bis) arc34;
dom = dom1 et dom2;
trac qual dom;

E1 = 210000. ;
nu1 = 0.2 ;
*
* Limite élastique MPa
Sy1 = 420. ;




*
* Comportement en traction (PLASTIQUE ISOTROPE)
*
ListE1 = PROG 0.0000 0.002 0.202 ;
*
ListC1 = PROG 0. 420. 741.643 ;
*
EvolMat1 = EVOL 'MANU' 'def' ListE1 'Cont' ListC1 ;

DESS EvolMat1 ;



MOD = MODEL dom MECANIQUE ELASTIQUE PLASTIQUE ISOTROPE QUA8;

MAT = MATE MOD 'YOUN' E1 'NU' nu1 'RHO' 7.85E-9 'TRAC' EvolMat1;


KK = rigi mod mat;




*CONDITIONS LIMITES ET FORCE
Ad12 = bloquer uy d12;
Ad41 = bloquer ux d41;

CL3= Ad12 ET Ad41;






PRES0 = PRESSION MASS MOD 1. (d34);
CHAR3 = PRES0 ;




*--- RESOLUTION:
RIG3 = RIGI MOD MAT ;
RIGCL3 = RIG3 ET CL3 ;
DEP3 = RESOU RIGCL3 CHAR3 ;
*
*--- POST TRAITEMENT
*DEFORMEE
DEF0 = DEFO DOM DEP3 0. BLAN ;
DEF1 = DEFO DOM DEP3 ROUG ;
TITR 'DEFORMEE APRES CHARGEMENT PRESSION' ;
TRAC (DEF0 ET DEF1) ;
*
*CALCUL DES CONTRAINTES TOTALES ET MECANIQUES
SIGTOT3 = SIGMA MOD MAT DEP3 ;
SIGM3 = SIGTOT3 ;

*
S_XX = (EXCO SMXX SIGM3) ;
TITR ' CONTRAINTES S_XX
(MAXI=' (MAXI S_XX) 'MPa)' ;
TRAC MOD S_XX ;
S_YY = (EXCO SMYY SIGM3) ;
TITR ' CONTRAINTES S_YY (MAXI=' (MAXI S_YY) 'MPa)' ;
TRAC MOD S_YY ;
*


CPNPLA=CHAN CHPO mod SIGM3;
* évolution le lond du ligament
EVSXX=EVOL CHPO CPNPLA SMXX arc23;
* tracé de cette évolution
dessin EVSXX;


*-----------------------------CHARGEMENT----------------------------*
*


LI1 = PROG 0. 1. 2. 3. 4. 5. 6. ;
LI2 = PROG 0. -100. 0. -100. 0. -100. 0. ;
EV = EVOL MANU 'Temps' LI1 'PRESSION' LI2 ;
CHA1 = CHAR MECA CHAR3 EV ;
DESS EV TITRE 'évolution pression en fonction du temps' ;
*
*--------------------------CALCUL PAS A PAS-------------------------*
*

TAB1 = TABLE ;
TAB1.'MODELE' = mod ;
TAB1.'CARACTERISTIQUES' = MAT ;
TAB1.'CHARGEMENT' = CHA1 ;
TAB1.'BLOCAGES_MECANIQUES' = CL3 ;
TAB1.'TEMPS_CALCULES' = LI1 ;
TAB1.'TEMPS_SAUVES' = LI1 ;
PASAPAS TAB1 ;



OPTI SAUV '\ecrou.txt';
SAUV TAB1;
FIN;
mon 2ieme fichier pr le post traitement :


OPTI REST '\ecrou.txt';
REST;
*
*--- point de chargement
*

PCHAR = 0;
*

SIGM1 = TAB1.CONTRAINTES.PCHAR ;
SIGEQ1 = VMIS MOD SIGM1 ;
TEMPO = tab1.temps.PCHAR ;
TITRE 'CONRAINTE DE VON MISES AU TEMPS' PCHAR;
TRAC MOD SIGEQ1 ;


PCHAR = 1;
*

SIGM1 = TAB1.CONTRAINTES.PCHAR ;
SIGEQ1 = VMIS MOD SIGM1 ;
TEMPO = tab1.temps.PCHAR ;
TITRE 'CONRAINTE DE VON MISES AU TEMPS' PCHAR;
TRAC MOD SIGEQ1 ;


PCHAR = 2;
SIGM1 = TAB1.CONTRAINTES.PCHAR ;
SIGEQ1 = VMIS MOD SIGM1 ;
TEMPO = tab1.temps.PCHAR ;
TITRE 'CONRAINTE DE VON MISES AU TEMPS' PCHAR;
TRAC MOD SIGEQ1 ;


PCHAR=1;

DEPPLA=TAB1.'DEPLACEMENTS'. PCHAR;
* déformée
DEFPLA0=defo DEPPLA dom 0. vert;
DEFPLA1=defo DEPPLA dom rouge;
* tracé de la déformée
TITRE 'deformée au temps' PCHAR;
trace (DEFPLA0 ET DEFPLA1);

PCHAR=2;

DEPPLA=TAB1.'DEPLACEMENTS'. PCHAR;
* déformée
DEFPLA0=defo DEPPLA dom 0. vert;
DEFPLA1=defo DEPPLA dom rouge;
* tracé de la déformée
TITRE 'deformée au temps' PCHAR;
trace (DEFPLA0 ET DEFPLA1);


PCHAR=3;

DEPPLA=TAB1.'DEPLACEMENTS'. PCHAR;
* déformée
DEFPLA0=defo DEPPLA dom 0. vert;
DEFPLA1=defo DEPPLA dom rouge;
* tracé de la déformée
TITRE 'deformée au temps' PCHAR;
trace (DEFPLA0 ET DEFPLA1);

PCHAR=4;

DEPPLA=TAB1.'DEPLACEMENTS'. PCHAR;
* déformée
DEFPLA0=defo DEPPLA dom 0. vert;
DEFPLA1=defo DEPPLA dom rouge;
* tracé de la déformée
TITRE 'deformée au temps' PCHAR;
trace (DEFPLA0 ET DEFPLA1);
fin;

**
Voir le profil de l'utilisateurEnvoyer un message privé
lguenhael
Administrateur

Inscrit le: 25 Nov 2006
Messages: 146
Localisation: Laboratoire Roberval - UTC
Répondre en citant
Bonjour,

Lorsque l'on trace une déformée, si aucun coefficient d'amplification de cette déformée n'est défini Cast3m en choisi un lui-même tel que le déplacement des nœuds soit d'un ordre en dessous de celui de la dimensions de la fenêtre.

Ainsi, lorsque la déformée que vous souhaitez afficher est proche de 0, Cast3M choisi un coefficient d'amplification qui va théoriquement tendre vers l'infini! (il est précisé en haut à droite)

Pour éviter cela, il vous suffit donc d'imposer vous même ce coefficient de déformée en choisissant par exemple 3.0E2.

Cordialement.
Voir le profil de l'utilisateurEnvoyer un message privéEnvoyer un e-mail
sam


Inscrit le: 29 Jan 2010
Messages: 3
Répondre en citant
vs avez raison, merci beaucoup pr votre réponse , et de m'avoir répondu aussi vite

sam.
Voir le profil de l'utilisateurEnvoyer un message privé
écrouissage mode pas a pas
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