Utilisation de Homard avec CREA_CHAMP

VincentMagnenet

Bonjour,

J'aimerais utiliser Homard pour affiner une zone de mon maillage,

mais j'ai l'impression qu'Aster refuse l'utilisation de MACR_ADAP_MAIL

avec un champ de type NOEU_DEPL_F. Cela est-il normal ?

(voir les fichiers joints).

Merci d'avance.

Attached file:
essai.zip, 15418kb

delmas

Oui cela est normal.

Avec CREA_CHAMP, vous n'avez défini qu'un champ de fonction. Il faut maintenant l'évaluer (toujours avec CREA_CHAMP) pour ne donner à MACR_ADAP_MAIL un champ de type *_R.

lIlIl

J'ai une erreur dans la version 17.1.9 avec CREA_CHAMP , EVAL et NOEU_DEPL_F. Je cherche à créer un champ de déplacement proche de la solution attendue pour permettre la convergence avec STAT_NON_LINE même en cas de points limite sur le chemin d'équilibre.

# Get mesh nodes (we need X, Y, Z as inputs in the FORMULE)
CHAMP_GEOM = CREA_CHAMP(OPERATION = 'EXTR',
                 TYPE_CHAM = 'NOEU_GEOM_R',
                 MAILLAGE  = mesh,
                 NOM_CHAM  = 'GEOMETRIE')

# Displacement formula (example)
# Note: the actual formula are more complex and time-dependent.
Formule_u_x = FORMULE(VALE='coeff*X1', NOM_PARA=('X1','X2','X3'), coeff=1.0)
Formule_u_y = FORMULE(VALE='coeff*X2', NOM_PARA=('X1','X2','X3'), coeff=1.0)
Formule_u_z = FORMULE(VALE='coeff*X3', NOM_PARA=('X1','X2','X3'), coeff=1.0)

# Directly create the field as 'NOEU_DEPL_F' instead of 'NOEU_NEUT_F'
champ_init_AFFE_direct = CREA_CHAMP(OPERATION = 'AFFE',
    TYPE_CHAM = 'NOEU_DEPL_F',
    MAILLAGE = mesh,
    # PROL_ZERO='OUI', # ERROR
    AFFE = (
    _F(TOUT='OUI', NOM_CMP='DX', VALE_F=Formule_u_x),
    _F(TOUT='OUI', NOM_CMP='DY', VALE_F=Formule_u_y),
    _F(TOUT='OUI', NOM_CMP='DZ', VALE_F=Formule_u_z),
    )
)

# FAILS: evaluate field
champ_init_EVAL = CREA_CHAMP(OPERATION='EVAL',
                   TYPE_CHAM='NOEU_DEPL_F',
                   CHAM_F=champ_init_AFFE,
                   CHAM_PARA=CHAMP_GEOM)

Erreur :

 ╔════════════════════════════════════════════════════════════════════════════════════════════════╗
 ║ <EXCEPTION> <MODELISA4_13>                                                                     ║
 ║                                                                                                ║
 ║ seule la grandeur NEUT_F est traitée actuellement.                                             ║
 ╚════════════════════════════════════════════════════════════════════════════════════════════════╝

J'ai tenté CREA_RESU comme dans sdld02b, wtnp110a ou wtnv141b (les seuls à mentionner 'NOEU_DEPL_F') :

CREA_RESU(OPERATION='AFFE',
TYPE_RESU='EVOL_NOLI',
NOM_CHAM='DEPL', # PROBLEM (see below)
AFFE=_F(CHAM_GD=champ_init_AFFE_direct,
INST=0.))

 ╔════════════════════════════════════════════════════════════════════════════════════════════════╗
 ║ <EXCEPTION> <SUPERVIS_4>                                                                       ║
 ║                                                                                                ║
 ║ Erreur de syntaxe dans CREA_RESU                                                               ║
 ║                                                                                                ║
 ║ For keyword NOM_CHAM: Unauthorized keyword: 'NOM_CHAM'                                         ║
 ║                                                                                                ║
 ║ Exception détaillée ci-dessous.                                                                ║
 ╚════════════════════════════════════════════════════════════════════════════════════════════════╝

et si je retire 'NOM_CHAM':

 ╔════════════════════════════════════════════════════════════════════════════════════════════════╗
 ║ <EXCEPTION> <SUPERVIS_4>                                                                       ║
 ║                                                                                                ║
 ║ Erreur de syntaxe dans CREA_RESU                                                               ║
 ║                                                                                                ║
 ║ 'Keyword NOM_CHAM is mandatory'                                                                ║
 ║                                                                                                ║
 ║ Exception détaillée ci-dessous.                                                                ║
 ╚════════════════════════════════════════════════════════════════════════════════════════════════╝

Je ne vois pas comment avancer à partir de là.

basile-simvia

Bonjour, la syntexe de CREA_RESU a changé un tout petit peu dans les versions 17, l'argument NOM_CHAM est à mettre dans le AFFE maintenant par exemple dans votre cas :

CREA_RESU(OPERATION='AFFE',
                    TYPE_RESU='EVOL_NOLI',
                    AFFE=_F(NOM_CHAM='DEPL', CHAM_GD=champ_init_AFFE_direct,
                    INST=0.))

et pour l'histoire du EVAL il me semble que votre champ doit être de type NOEU_NEUT_F pour être évaluable.

lIlIl

Merci Basile.

Voici le code que j'ai réussi à faire tourner :

# Get the geometry field (X,Y,Z)
CHAMP_GEOM = CREA_CHAMP(OPERATION = 'EXTR', # extraire
                     TYPE_CHAM = 'NOEU_GEOM_R',
                     MAILLAGE = mesh,
                     NOM_CHAM = 'GEOMETRIE')

# Define your transformation functions

def approximation_of_u_x(X,Y,Z, INST, max_force):
    return ... # define your own

def approximation_of_u_y(X,Y,Z, INST, max_force):
    return ... # define your own

def approximation_of_u_z(X,Y,Z, INST, max_force):
    return ... # define your own

# Time-dependent formulas

Formule_zero = FORMULE(VALE='0', NOM_PARA=('X','Y','Z', 'INST'))

params_formule = {
    'VALE': 'f(X,Y,Z, INST, max_force)',
    'NOM_PARA': ('X','Y','Z', 'INST'),
    'max_force': 1500 # Newton
}
Formule_u_x = FORMULE(f=approximation_of_u_x, **params_formule)
Formule_u_y = FORMULE(f=approximation_of_u_y, **params_formule)
Formule_u_z = FORMULE(f=approximation_of_u_z, **params_formule)

# Create a displacement field
DEPL_init_AFFE_direct = CREA_CHAMP(OPERATION = 'AFFE',
    TYPE_CHAM = 'NOEU_DEPL_F',
    MAILLAGE = mesh,
    # PROL_ZERO='OUI', # ne fonctionne pas avec NODE_DEPL_F
    AFFE = (
        _F(
            GROUP_NO = 'most_of_the_nodes',
            NOM_CMP = ('DX','DY','DZ'),
            VALE_F = (Formule_u_x, Formule_u_y, Formule_u_z)),
        _F(
            GROUP_NO = 'some_nodes_that_need_to_stay_at_zero',
            NOM_CMP = ('DX','DY','DZ'),
            VALE_F = (Formule_zero, Formule_zero, Formule_zero)),
    )
)

# Create a result (RESULTAT/EVOL_NOLI)
DEPL_init_RESU = CREA_RESU(OPERATION = 'AFFE',
    TYPE_RESU = 'EVOL_NOLI',
    AFFE = _F(CHAM_GD = DEPL_init_AFFE_direct,
    NOM_CHAM = 'DEPL',
    LIST_INST = lst_reels)
)

# # then:
# result = STAT_NON_LINE( ETAT_INIT=_F(EVOL_NOLI=DEPL_init_RESU), ... )

Malheureusement ETAT_INIT n'a pas eu le comportement que j'espérais : la résolution avec STAT_NON_LINE est initialisée avec le dernier instant de DEPL_init_RESU, puis reprend ce dernier état convergé pour calculer les paliers de charge suivants (il ignore tous les instants de mon résultat "guide"). J'aurais voulu que STAT_NON_LINE se base sur mon champ approximé pour initialiser chaque pas de temps et accompagner la résolution non linéaire.

basile-simvia

lIlIl peux tu me donner le cas complet .comm, .export et .med ?

AsterO'dactyle

lIlIl

PREDICTION='DEPL_CALCULE' dans le mot-clef facteur NEWTON devrait faire ça