Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2019-11-13 09:39:40

Mahions
Member
Registered: 2019-05-13
Posts: 26

How to define a temperature field along the direction of z axis?

1. I want to impose a temperature field along z axis (not a function, just specific value), but I don't know how to define it?

2. Look at the follow message, what is the meaning of VALE_REF? initial temperature, is that right? Can it be defined by a form of field?

temp1 = CREA_CHAMP(
  AFFE=_F(
    NOM_CMP=('TEMP', ),
    TOUT='OUI',
    VALE=(300.0, )
  ),
  MAILLAGE=mesh,
  OPERATION='AFFE',
  TYPE_CHAM='NOEU_TEMP_R'
)
fieldma = AFFE_MATERIAU(
  AFFE=_F(
    MATER=(mater, ),
    TOUT='OUI'
  ),
  AFFE_VARC=_F(
    CHAM_GD=temp1,
    NOM_VARC='TEMP',
    VALE_REF=0.0
  ),
  MAILLAGE=mesh
)


3. I tried to use LIAISON_ELEM to connect DKT and POU_D_E,but failed. The wrong message :


# Commande No :  0010            Concept de type : char_meca
  # ------------------------------------------------------------------------------------------
  load = AFFE_CHAR_MECA(MODELE=model,
                        LIAISON_ELEM=_F(GROUP_NO_2='connect1',
                                        GROUP_MA_1=('load', ),
                                        OPTION='COQ_POU',
                                        AXE_POUTRE=(0.0, 0.0, 1.0),
                                        CARA_ELEM=elemprop,
                                        NUME_LAGR='NORMAL',
                                        ANGL_MAX=1.0,),
                        DDL_IMPO=_F(LIAISON='ENCASTRE',
                                    GROUP_NO=('bcpoint', ),),
                        FORCE_NODALE=_F(GROUP_NO=('bc', ),
                                        FZ=100.0,),
                        INFO=1,
                        VERI_NORM='OUI',)


Le système linéaire à résoudre a 18606 degrés de liberté:
   - 18606 sont des degrés de liberté physiques
     (ils sont portés par 3101 noeuds du maillage)
   - 0 sont les couples de paramètres de Lagrange associés
     aux 0 relations linéaires dualisées.
   
   !--------------------------------------------------------------------------------!
   ! <EXCEPTION> <CALCUL_30>                                                        !
   !                                                                                !
   !  Erreur de programmation :                                                     !
   !  on n'a pas pu extraire toutes les composantes voulues du champ global associé !
   !  au paramètre: PCACOQU  (option: CARA_SECT_POUT3  type_élément: MEBODKT )      !
   !                                                                                !
   ! --------------------------------------------                                   !
   ! Contexte du message :                                                          !
   !    Option         : CARA_SECT_POUT3                                            !
   !    Type d'élément : MEBODKT                                                    !
   !    Maillage       : mesh0                                                      !
   !    Maille         : M101                                                       !
   !    Type de maille : SEG2                                                       !
   !    Cette maille appartient aux groupes de mailles suivants :                   !
   !       load                                                                     !
   !    Position du centre de gravité de la maille :                                !
   !       x=5.000000 y=0.000000 z=0.000000                                         !
   !--------------------------------------------------------------------------------!
   
   Maille concernée : M101
   Option de calcul : CARA_SECT_POUT3
   Commentaire:
     Paramètre d'entrée : PCACOQU
      Commentaire:
PCACOQU  :  PROPRIETES COQUES
PCACOQU : CARACTERISTIQUE DE COQUE
PCACOQU : CARACTERISTIQUE DE COQUE
PCACOQU : CARACTERISTIQUE DE COQUE,
NECESSITE DE FOURNIR LE CONCEPT PRODUIT PAR AFFE_CARA_ELEM
Champ de caracteristiques pour les coques. Concept CARA_ELEM
PCACOQU : CARACTERISTIQUES DE COQUE,
NECESSITE DE FOURNIR LE CONCEPT PRODUIT PAR AFFE_CARA_ELEM
CARACTERISTIQUES DES ELEMENTS DE COQUE
PCACOQU : CARACTERISTIQUE DE COQUE
PCACOQU : CARACTERISTIQUE DE COQUE,
NECESSITE DE FOURNIR LE CONCEPT PRODUIT PAR AFFE_CARA_ELEM
SHELLS: CARACTERISTIQUE DE COQUE

     Grandeur associée au champ : CACOQU
      Commentaire:
CACOQU  Type:R  Caracteristiques des coques
EP      : epaisseur de la coque
SECT_L  : somme des sections d'armatures dans la direction L
ALPHA   : 1er angle de l'axe de reference dans le plan transverse
BETA    : 2eme angle de l'axe de reference dans le plan transverse
Ces deux angles definissent par rapport au repere de reference
le vecteur a projeter sur le plan tangent de l'element afin
d'y definir le repere (T,N,L).
La direction L est perpendiculaire au plan tangent.
KAPPA   : coefficient de cisaillement transverse
C_METR  : coefficient de modification metrique pour l'element coque
DIST_N  : excentrement de la nappe d'armatures vs la maille support
CTOR    : constante multiplicative de la rigidite en rotation associee a
la direction perpendiculaire au plan tangent de l'element
EXCENT  : excentrement du feuillet moyen
INERTIE : prise en compte de l'inertie de rotation
TENS    : (réservé aux membranes) application d'une contrainte
initiale pour faire converger à la première iteration.
Elle disparait aux increments suivants


What does  PCACOQU  and CARA_SECT_POUT3 mean? I can't find it in document. In addition, there are unkown elements (named MEBODKT)in the model.


model = AFFE_MODELE(MAILLAGE=mesh0,
                      AFFE=(_F(PHENOMENE='MECANIQUE',
                               TOUT='OUI',
                               MODELISATION=('DKT', ),),
                            _F(GROUP_MA=('beam', ),
                               PHENOMENE='MECANIQUE',
                               MODELISATION=('POU_D_E', ),),
                            ),
                      INFO=1,
                      VERI_JACOBIEN='OUI',
                      DISTRIBUTION=_F(PARTITIONNEUR='METIS',
                                      METHODE='SOUS_DOMAINE',),
                      )


Sur les 3460 mailles du maillage mesh0, on a demandé l'affectation de 3460, on a pu en affecter 3460
Modélisation     Type maille  Élément fini     Nombre
                  SEG2         MEBODKT          420
POU_D_E          SEG2         MECA_POU_D_E     40
DKT              QUAD4        MEDKQU4          3000


Your help is urgently needed.

Offline