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
fluage maxwell [Résolu]
LaMoucheGueulante


Inscrit le: 21 Mar 2007
Messages: 10
Répondre en citant
je voudrais utiliser un modèle "mecanique elastique isotrope fluage maxwell" (à 8 cellules), je le teste pour l'instant sur un simple carrée en 2D en qua4 (compression simple, effort imposé), mais j'ai au lancement de la procédure pasapas le message d'erreur suivant: "erreur 552 dans l'opérateur COMP , le type REAL*8 de la composante EM0 du champ par élément utile ne correspond pas à celui attendu..."
J'ai consulté modflu.eso et le type attendu est bien REAL*8 (j'ai d'ailleurs tester d'autres types et le message est toujours le même).
J'ai repris un fichier avec un modèle de fluage de type norton qui tournait bien en remplaçant simplement le MATEriaux et le MODele et l'erreur apparaît également...
Comment résoudre ce problème?
Voir le profil de l'utilisateurEnvoyer un message privé
LaMoucheGueulante


Inscrit le: 21 Mar 2007
Messages: 10
Répondre en citant
J'ai trouvé! Les modules élastiques de chaque branche doivent être des évolutions des modules avec la température, et les temps propre des réels...C'était si simple Confused
Voir le profil de l'utilisateurEnvoyer un message privé
OM


Inscrit le: 05 Fév 2008
Messages: 1
Répondre en citant
Bonjour,
J'utilise un modele fluage maxwell. Pour identifier les paramètres j'utilise la procédure identi qui extrait :
1- une évolution Y1 des modules d'young en fonction du temps,
2- une liste Y2 contenant les temps de relaxation des 4 branches
3- une liste Y3 contenant les modules des 4 branches.
Mon problème est : comment introduire les paramètres et sous quelle forme dans le matériau
ci joint extrait du prog.
J'ai créé 4 listes (E1 à E4) que j'ai affecté aux EMi
Je souhaiterai savoir comment introduire les paramètres dans ce cas.

Y1 Y2 Y3 = identi maxwell jour 90. 4 EUROCODE 0.075 50 0.2 53.3;
List Y1;
List Y2;
List Y3;
prog1=prog 0.;
prog2=prog 48205.;
E0=evol manu prog1 prog2;
prog1=prog 0.1e+1;
prog2=prog 48223;
E1=evol manu prog1 prog2;
prog1=prog 0.1;
prog2=prog 48241;
E2=evol manu prog1 prog2;
prog1=prog 1.0E-2;
prog2=prog 48259;
E3=evol manu prog1 prog2;
prog1=prog 0.1E-2;
prog2=prog 48277;
E4=evol manu prog1 prog2;

matecyl = mate modecyl youn 30550E+6 nu 0.16 TR1 0.1E+1 TR2 0.1
TR3 1.0E-2 TR4 0.1E-2 EM0 E0 EM1 E1 EM2 E2 EM3 E3 EM4 E4;
Merci
OM Question
Voir le profil de l'utilisateurEnvoyer un message privé
PHAN


Inscrit le: 05 Oct 2011
Messages: 4
Répondre en citant
Bonjour,

Je voudrais faire un calcul sur le modele fluage Maxwell avec 4 braches.
Je ne sais pas comment on peut le faire?

Je viend de modifier le fichier norton_cis1.dgibi (exemple) en ajoutant les valeurs ci-dessous.

Mais je n'ai pas arrivé. Il y a ERREUR 996... "Nombre maximum de sous-pas atteint, arret de pasapas...."
Pourriez-vous m'aider le corriger.

Merci beaucoup.

PHAN

Voici la code modifié:

* fichier : norton_cis1.dgibi
* Test Norton_cis1.dgibi: Jeux de données *
* --------------------------------------- *
* *
*
**************************************************
* *
* TEST DE VALIDATION D'UNE LOI DE COMPORTEMENT *
* -------------------------------------------- *
* DE MATERIAU *
* ----------- *
* *
* LOI DE COMPORTEMENT DE FLUAGE DE : *
* NORTON *
* *
* TEST POUR DES ELEMENTS MASSIFS *
* ------------------------------ *
* *
* MAILLAGE: *
* UNE PLAQUE CARREE DE COTE L=1 M *
* *
* CHARGEMENT: *
* CISAILLEMENT PUR *
* FORCES IMPOSEES *
* DEFORMATIONS PLANES *
* *
**************************************************
opti echo 0 dime 2 elem qua4 ;

*
* Maillage
*
p1 = 0. 0. ;
p2 = 1. 0. ;
p3 = 1. 1. ;
p4 = 0. 1. ;
li1 = d p1 1 p2 ;
li2 = d p2 1 p3 ;
li3 = d p3 1 p4 ;
li4 = d p4 1 p1 ;
mail1 = 'SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN' ;
*
* Modele
*
AF1 = 1.E-45 ;
AF2 = 5. ;
AF3 = 1.2 ;
SMAX = 2.1E8;
youn0 = 1.5E+4;
nu0 = 0.3 ;

prog1=prog 0.;
prog2=prog 1.5E+4;
E0=evol manu prog1 prog2;

prog1=prog 0.1E+1;
prog2=prog 1.11E+6;
E1=evol manu prog1 prog2;

prog1=prog 0.1;
prog2=prog 1.12E+6;
E2=evol manu prog1 prog2;

prog1=prog 1.0E-2;
prog2=prog 3.43E+5;
E3=evol manu prog1 prog2;

prog1=prog 0.1E-2;
prog2=prog 4.14E+5;
E4=evol manu prog1 prog2;

TRR1 = 0.46;
TRR2 = 0.45;
TRR3 = 3.62;
TRR4 = 45.1;

mod0 = 'MODE' mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE'
'FLUAGE' 'MAXWELL';
mat0 = 'MATE' mod0 'YOUN' youn0 'NU' nu0 'EM0' E0
'EM1' E1 'EM2' E2 'EM3' E3 'EM4' E4
'TR1' TRR1 'TR2' TRR2 'TR3' TRR3 'TR4' TRR4;
*
* Conditions aux limites
*
cl1 = 'BLOQ' mail1 'UY' ;
cl2 = 'BLOQ' li1 'UX' ;
cl0 = cl1 'ET' cl2 ;
*
* Chargement
*
depmax = 7.E-1 ;
chp1 = 'COOR' 2 mail1 ;
chp1 = chp1 * depmax / 1. ;
chp1 = 'NOMC' chp1 'UX' ;
titr ' Champ de forces impose' ;
*trace chp1 mail1 ;
chp2 = 'MANU' 'CHPO' mail1 1 'UY' 0. ;
chp0 = chp1 + chp2 ;
dep1 = 'FORCE' li1 'FX' ( -1. * depmax ) 'FY' 0. ;
dep2 = 'FORCE' li3 'FX' depmax 'FY' 0. ;
dep0 = dep1 'ET' dep2 ;
ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-4 1000. )
y ( 'PROG' 0. 1. 1. ) ;
cha0 = 'CHAR' 'MECA' dep0 ev0 ;
*
dt0 = 0.020 ;
dt1 = 2. * dt0 ;
dt2 = 5. * dt1 ;
dt3 = 5. * dt2 ;
tfin0 = 20. ;
*
* Appel à PASAPAS
*
ta1 = 'TABLE' ;
ta1 .'MODELE' = mod0 ;
ta1 .'CARACTERISTIQUES' = mat0 ;
ta1 .'CHARGEMENT' = cha0 ;
ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' 1.E-4 2.E-4 ) 'ET'
( 'PROG' 3.E-4 dt0 ) 'ET'
( 'PROG' dt1 'PAS' dt0 dt2 ) 'ET'
( 'PROG' ( dt2 + dt2 ) 'PAS' dt2 dt3 ) 'ET'
( 'PROG' ( dt3 + dt3 ) 'PAS' dt3 tfin0 ) ;
*
list ta1.'TEMPS_CALCULES';
TMASAU=table;
ta1. 'MES_SAUVEGARDES'=TMASAU;
TMASAU.'DEFTO'=VRAI;
TMASAU.'DEFIN'=VRAI;
*
PASAPAS ta1;
*
* Post traitement
*
dim0 = 'DIME' ta1 .'TEMPS' ;
eps10 = 'PROG' 0. ;
sigm10 = 'PROG' 0. ;
def10 = 'PROG' 0. ;
def_th2 = 'PROG' 0. ;
sig_th2 = 'PROG' 0. ;
err_e2 = 'PROG' 0. ;
err_s2 = 'PROG' 0. ;
temp2 = 'PROG' 0. ;
rac0 = ( 3. ) ** ( .5 ) ;
H0 = youn0 / ( 2. * ( 1. + nu0 ) ) ;
*
opti echo 1;
prob0 = vrai ;
'REPETER' bloc0 ( dim0 - 1 ) ;
i1 = &bloc0 ;
*
* Résultats calculés
*
def0 = ta1 .'DEFORMATIONS_INELASTIQUES'.i1 ;
depl0 = ta1 .'DEPLACEMENTS'.i1 ;
sigm0 = ta1 .'CONTRAINTES'.i1 ;
temp0 = ta1 .'TEMPS'.i1 ;
var0 = ta1 .'VARIABLES_INTERNES'.i1 ;
eps0 = 'EPSI' mod0 depl0 ;
eps1 = 'EXTR' ( 'EXCO' eps0 'GAXY' 'SCAL' ) 'SCAL' 1 1 1 ;
sigm1 = 'EXTR' ( 'EXCO' sigm0 'SMXY' 'SCAL' ) 'SCAL' 1 1 1 ;
def1 = 'EXTR' ( 'EXCO' def0 'GIXY' 'SCAL' ) 'SCAL' 1 1 1 ;
eps10 = eps10 'ET' ( 'PROG' eps1 ) ;
sigm10 = sigm10 'ET' ( 'PROG' sigm1 ) ;
sigeq0 = 'VMIS' sigm0 mod0 ;
sigeq1 = 'EXTR' sigeq0 'SCAL' 1 1 1 ;
def10 = def10 'ET' ( 'PROG' def1 ) ;
temp2 = temp2 'ET' ( 'PROG' temp0 ) ;
var1 = 'EXTR' ( 'EXCO' var0 'EPSE' 'SCAL' ) 'SCAL' 1 1 1 ;
sig0 = 'CHAN' 'CHPO' eps0 mod0 ;
sigx = 'EXCO' sig0 'EPXX' ;
sigy = 'EXCO' sig0 'EPYY' ;
sigxy = 'EXCO' sig0 'GAXY' ;
*titr ' Contraintes suivant X' ;
trace sigx mail1 ;
*titr ' Contraintes suivant Y' ;
trace sigy mail1 ;
*titr ' Contraintes de cisaillement' ;
trace sigxy mail1 ;
*
* Résultats analytiques
*
y2 = ( ( rac0 * depmax ) ** ( AF2 ) ) ;
y2 = y2 * ( temp0 ** AF3 ) ;
x2 = y2 * AF1 * rac0 ;
*
def_th1 = x2 ;
sig_th1 = H0 * ( eps1 - x2 ) ;
*
def_th2 = def_th2 'ET' ( 'PROG' def_th1 ) ;
sig_th2 = sig_th2 'ET' ( 'PROG' sig_th1 ) ;
*mess ' sig calculée =' sigm1 ;
*mess ' sig theorique =' sig_th1 ;
*mess ' def calculée =' def1 ;
*mess ' def théorique =' def_th1 ;
*
* Erreurs
*
'SI' ( def1 '>' ( 1.E-1 * eps1 ) ) ;
err_e1 = ( 'ABS' ( def1 - def_th1 ) ) / def1 ;
err_e1 = err_e1 * 100. ;
'SINON' ;
err_e1 = 0. ;
'FINSI' ;
'SI' ( sigm1 '>' 0. ) ;
err_s1 = ( 'ABS' ( sigm1 - sig_th1 ) ) / sigm1 ;
err_s1 = err_s1 * 100. ;
'SINON' ;
err_s1 = 0. ;
'FINSI' ;
*
'SI' ( ( err_e1 '>' .25 ) 'OU' ( err_s1 '>' .25 ) ) ;
prob0 = faux ;
mess 'err_e1 err_s1' err_e1 err_s1 ;
'FINSI' ;
err_e2 = err_e2 'ET' ( 'PROG' err_e1 ) ;
err_s2 = err_s2 'ET' ( 'PROG' err_s1 ) ;
*
'MESS' ' Fin du calcul de restitution au temps ' temp0 ;
*mess ' ' ;
*
'FIN' bloc0 ;
*
defo0 = 'DEFO' mail1 depl0 0. ;
defo1 = 'DEFO' mail1 depl0 1. 'ROUG' ;
*titr ' deformee finale' ;
*trace ( defo0 'ET' defo1 ) ;
*
'SI' ( prob0 'EGA' faux ) ;
'ERRE' 5 ;
'SINON' ;
'ERRE' 0 ;
'FINSI' ;
*
* Tracés
*
evo1 = 'EVOL' 'MANU' 'epsp' def10 'sigm' sigm10 ;
evo2 = 'EVOL' 'MANU' 'epsp' def_th2 'sigm' sig_th2 ;
evo3 = 'EVOL' 'MANU' 'epsp' def10 'erreur_%' err_e2 ;
evo4 = 'EVOL' 'MANU' 'epsp' def10 'erreur_%' err_s2 ;
evo5 = 'EVOL' 'MANU' 'temps' temp2 'epsp' def10 ;
evo6 = 'EVOL' 'MANU' 'temps' temp2 'epsp' def_th2 ;
TAB = 'TABLE' ;
TAB. 1 = 'MARQ CROI ' ;
TAB. 2 = 'MARQ PLUS ' ;
TAB.'TITRE' = 'TABLE' ;
TAB.'TITRE'. 1 = MOT 'CALCUL' ;
TAB.'TITRE'. 2 = MOT 'THEORIE' ;
'DESSIN' ( evo1 'ET' evo2 )
'TITR' ' Comparaison calcul theorie'
'LEGE' TAB ;
TAB = 'TABLE' ;
TAB. 1 = 'MARQ CROI ' ;
*TAB. 2 = 'MARQ PLUS ' ;
TAB.'TITRE' = 'TABLE' ;
TAB.'TITRE'. 1 = MOT 'DEFORMATION' ;
*TAB.'TITRE'. 2 = MOT 'CONTRAINTE' ;
'DESSIN' ( evo3 )
'TITR' ' Evolution des erreurs (en pourcentage) '
'LEGE' TAB ;
TAB = 'TABLE' ;
TAB. 1 = 'MARQ CROI ' ;
TAB. 2 = 'MARQ PLUS ' ;
TAB.'TITRE' = 'TABLE' ;
TAB.'TITRE'. 1 = MOT 'CALCUL' ;
TAB.'TITRE'. 2 = MOT 'THEORIE' ;
'DESSIN' ( evo5 'ET' evo6 )
'TITR' ' Comparaison calcul theorie, deformation-temps'
'LEGE' TAB ;
*
'FIN' ;



Ex Fluage.dgibi
 Description:

Télécharger
 Nom du fichier:  Ex Fluage.dgibi
 Taille du fichier:  7.09 Ko
 Téléchargé:  402 fois

Voir le profil de l'utilisateurEnvoyer un message privé
Fluage
imen


Inscrit le: 01 Mai 2013
Messages: 2
Localisation: Tunisie
Répondre en citant
Bonjour,

Je voudrai faire un calcul de fluage dans Castem avec la loi de double puissance J(t)=(1/E0)x(1+phi x t'^(-m)x (t-t')^(n)
Je ne vois pas comment faire. Je n'ai rien trouvé sur le site à ce propos.

Est-ce que quelqu'un aurait un exemple d'un tel calcul? Prière de m'aidez.....

Et Merci.
Voir le profil de l'utilisateurEnvoyer un message privé
fluage maxwell [Résolu]
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