sam
Inscrit le: 29 Jan 2010 |
Messages: 3 |
|
|
|
Posté le: 09 Fév , 2010 19:54 |
|
|
|
|
|
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;
**
|
|
lguenhael
Administrateur
Inscrit le: 25 Nov 2006 |
Messages: 146 |
Localisation: Laboratoire Roberval - UTC |
|
|
Posté le: 09 Fév , 2010 21:30 |
|
|
|
|
|
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.
|
|
sam
Inscrit le: 29 Jan 2010 |
Messages: 3 |
|
|
|
Posté le: 09 Fév , 2010 21:41 |
|
|
|
|
|
vs avez raison, merci beaucoup pr votre réponse , et de m'avoir répondu aussi vite
sam.
|
|