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
Paramètre matériel en fonction du temps ( PERSO1 )
François131


Inscrit le: 10 Juin 2010
Messages: 1
Répondre en citant
Bonjour,

Je souhaite réaliser une évolution d'un paramètre matériel (TAB1.CARACTERISTIQUES) en fonction du temps ( et non en fonction d'un chargement, thermique par exemple ).
L'évolution en fonction du temps est justifié par le calcul de la réduction de mon matériel.

Utiliser un objet évolution ne suffit pas, puisqu'en abscisse, un chargement est demandé ( hors, je veux un temps à la place ).

J'ai donc pensé à changer mes propriétés matérielles en utilisant la procédure PERSO1 pour pouvoir intervenir au temps de calcul qui m'interesse.
J'arrive à faire changer TAB1.CARACTERISTIQUES ( en le vérifiant après PASAPAS avec list TAB1.CARACTERISTIQUES; ) mais le calcul n'est pas modifié...

Je joins ici la procédure PERSO1 que j'ai utilisé :
Je précise que je travaille sur un exemple simple afin de vérifier le calcul ( cad le chargement thermique est simplifié ).
Le changement de parametre se fait à température constante.


-----------------

*Le chargement thermique suit cette liste
LIST_DT = PROG 0. (T_ambia-T_fritt) (T_ambia-T_fritt+1) 0.;
*Je vous passe l'écriture du chargement en lui-même.
( en gros je reviens au point de départ et je compare avec l'arrivée, or je suis en calcul élastique => je devrais voir une modification en changeant un paramètre )

debp perso1 TAB_P*'TABLE' ;

*RANG DU TEMPS_CALCULES
n_rang = TAB_P.cpt_tp ; *cpt_tp est un compteur, qui fournit le numéro du rang (n_rang)
TAB_P.cpt_tp = (TAB_P.cpt_tp) + 1 ;

si (n_rang EGA 2); * je change le CTE ( alpha si vous préférer )
mess 'Changement de CTE' ;
t1=(T_ambia-T_fritt) ; * afin qu'il y ait un changement , je crée une évolution de CTE pour un chargement de température fictif (delta=1)
t2=(T_ambia-T_fritt+1);
L_T_f = PROG t1 t2;
L_CTE_f= PROG CTE_CO CTE_CR; *CTE_Cermet oxydé et CTE_cermet réduit
EV_CTE = EVOL 'MANU' 'T' L_T_f 'ALPH' L_CTE_f;

*PROPRIETES MATERIELLES
MAT_OA_P = MATER MOD_OA 'YOUN' E_CO 'NU' MU_CO 'ALPH' EV_CTE;
MAT_AB_P = MATER MOD_AB 'YOUN' E_EM 'NU' MU_EM 'ALPH' CTE_EM;

MAT_TOTP=MAT_OA_P ET MAT_AB_P;

TAB_P.'CARACTERISTIQUES'=MAT_TOTP;
finsi;

finp TAB_P ;

( J'ai simplifie ici en enlevant une étape supplémentaire : au rang suivant, je remplace le CTE par une valeur fixe ( CTE_CR ) au lieu de l'objet évolution )


Voilà, si quelqu'un a déjà rencontrer ce problème ( Evolution de TAB.CARACTERISTIQUES en fonction du temps)...

Merci.

PS : j'ai essayé d'utiliser WTAB, ou PRECED et MMMM ( cf procedure PASAPAS sur le site de CASTEM ) mais je n'y arrive pas..
Voir le profil de l'utilisateurEnvoyer un message privé
Paramètre matériel en fonction du temps ( PERSO1 )
cueffgui


Inscrit le: 05 Jan 2012
Messages: 8
Répondre en citant
Bonjour,

Le sujet date de 2010 mais on ne sait jamais. Avez vous trouvé une solution pour le problème de mise à jour des paramètres du matériau à chaque itération ?

Pour ma part je souhaite mettre à jour le module d'élasticité d'un matériau en fonction de l'avancement d'un réaction.

Je sais qu'il est possible de mettre à jour les paramètres thermiques à chaque pas de temps via la procédure PARA_THER mais elle ne semble pas permettre la mise à jour des caractéristiques mécaniques.

PERSO1 fonctionne-t-elle pour cela?

Merci d'avance.

G.CUEFF
Voir le profil de l'utilisateurEnvoyer un message privé
Paramètre matériel en fonction du temps ( PERSO1 )
Vincecro


Inscrit le: 01 Fév 2015
Messages: 1
Répondre en citant
Bonjour,

Je rencontre le même problème, J'ai un calcul dans lequel je souhaite modifier des rigidités additionnelles et l'amortissement d'une structure en fonction de l'indice de mon calcul PASAPAS.

Mon code est le suivant :

********************************
DATA=TABLE;
DATA.'MODAL'=TABMOD;
DATA.'MAS'=MASSETOT;
DATA.'CONNEX'=CONNEX;
DATA.'AMO1'=MATAMO1;
DATA.'AMO2'=MATAMO2;
DATA.'RIG1'=CLADD1;
DATA.'RIG2'=CLADD2;
**********************************
TAB1 = TABLE ;
TAB1. 'BLOCAGES_MECANIQUES' = KCL;
TAB1. 'MODELE' = MODETOT;
TAB1. 'CHARGEMENT' = CHARTOT;
TAB1. 'CARACTERISTIQUES' = MATETOT;
TAB1. 'PROCEDURE_PERSO1' = VRAI ;
TAB1. 'TEMPS_CALCULES' = LISCALC;
TAB1. 'MESDONNEES' = DATA;
TAB1.'CHARGEMENT' = CHARTOT;
TAB1.'DYNAMIQUE' = VRAI;
TAB1.'TEMPS_CALCULES' = SINETEMP;
****************************************************
DEBP PERSO1 TTU*'TABLE' ;
*NUMERO DU PAS ET INDIQUATION DE L'ETAT DU SYSTEME
NUM = DIME( TTU . 'TEMPS' );
A = TTU . 'MESDONNEES' . 'CONNEX' ;
B = EXTR A NUM;
LISTE NUM ;
WTAB = TTU . 'WTABLE' ;
LISTE B;
*CHOIX DES PARAMETRES D'AMORTISSEMENT ET DE RIGIDITE
SI (B EGA 1)
MESS ;
MESS ' *****************';
MESS ' RAIDEUR MIN' ;
MESS ' *****************' ;
MESS ;
WTAB.'AMORTISSEMENT'= TTU.'MESDONNEES'.'AMO2';
WTAB.'RIGIDITE_CONSTANTE'=TTU.'MESDONNEES'.'RIG2';
SINO;
MESS ;
MESS ' *****************';
MESS ' RAIDEUR MAX' ;
MESS ' *****************' ;
MESS ;
WTAB.'AMORTISSEMENT'= TTU.'MESDONNEES'.'AMO1';
WTAB.'RIGIDITE_CONSTANTE'=TTU.'MESDONNEES'.'RIG1';
FINS;
FINP ;
PASAPAS TAB1 ;



J'ai essayé pas mal de façon différente pour essayé d'appliquer le changement de rigidité et d'amortissement mais je tombe toujours sur le même message d'erreur (ci-dessous).
Ce message indique un problème de conditions au limites pourtant celles-ci ne sont pas modifié lors des changement que j'applique. J'ai essayé de faire oublier les paramètres de rigidité et d'amortissement avant d'en appliquer de nouveaux mais rien de change. Ce calcul fonctionne sans problème sans la procédure PERSO

Quelqu'un aurait il une idée ?

Merci d'avance,

Vincent



erreur.PNG
 Description:

Télécharger
 Nom du fichier:  erreur.PNG
 Taille du fichier:  12.28 Ko
 Téléchargé:  39 fois

Voir le profil de l'utilisateurEnvoyer un message privéEnvoyer un e-mail
Paramètre matériel en fonction du temps ( PERSO1 )
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