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

You are not logged in.

#1 Re: Code_Aster usage » What's wrong with "Saturated file"? » 2020-04-28 04:27:19

Well, I added "- max_base 250000 ", then astk run. Besides, I have to change temoraray directory for Aster executions(rep_trav).
Actually, it took 6 days to complete the calculation.How to speed-up it with hundreds of contact pairs?And there is a warning in mess file.

!----------------------------------------------------------------------------------------------------------------------------------!
   ! <A> <CONTACT_22>                                                                                                                 !
   !                                                                                                                                  !
   ! Contact méthode continue.                                                                                                        !
   !   La méthode ALGO_CONT='PENALISATION  autorise une pénétration des zones de contact. On estime cette pénétration à 1.565956 pour !
   ! cent par rapport à la plus petite maille du maillage.                                                                            !
   !   COEF_PENA_CONT doit être suffisamment grand pour empêcher une trop grande pénétration. Il faut donc augmenter la valeur de     !
   ! COEF_PENA_CONT (et de COEF_PENA_FROT )                                                                                           !
   !   de sorte à inférieur à cinq pour cent de pénétration.                                                                          !
   !   Conseils :                                                                                                                     !
   !   -------                                                                                                                        !
   !   Il n'y a pas de méthode de référence pour le choix des coefficients de pénalisation.                                           !
   !   - une estimation empirique serait de multiplier COEF_PENA_CONT par (1+1.565956)/100*module d'Young du corps le plus dur        !
   !   - Il n'y a pas d'estimation empirique pour COEF_PENA_FROT.                                                                     !
   !                                                                                                                                  !

What's the meaning of this warning? Should I choose PENALISATION method?

#2 Re: Code_Aster usage » Compound mesh VS. Liaison_Mail » 2020-04-28 04:01:37

I have same problem. I think it's OK for shell and 3D elements which share their nodes(compound mesh). However, something is wrong when I simulate a simple expansion case like this.
!------------------------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <FACTOR_10>                                                                                                      !
   !                                                                                                                              !
   ! Problème : la matrice est singulière ou presque singulière :                                                                 !
   !   Lors de la factorisation de la matrice, on a rencontré un problème                                                         !
   !   (pivot nul ou presque nul) à la ligne 2850 qui correspond au degré de liberté donné ci-dessus.                             !
   !   Le nombre de décimales "perdues" pour cette ligne est : 13.                                                                !
   !   Ce nombre de décimales perdues est lié au degré de singularité de la matrice.                                              !
   !   Plus il est grand, plus le système est singulier.                                                                          !
   !   Quand on a perdu plus de 8 décimales, on estime que l'on a perdu 50% de la précision                                       !
   !   des nombres réels (qui ont 15 décimales environ).

#3 Re: Code_Aster usage » What's wrong with "Saturated file"? » 2020-04-21 02:51:39

jeanpierreaubry wrote:

how are you running the calculation astk or as_run?
how much memory did you allocate?

Hello,
I run the case by astk version 2017. And the version of code-aster is stable, lower than 13.4.The total memory is 20480 MB.
As I attached, the model has only three laps,actually 16 laps. I don't know how to speed up it and whether my computer can solve complete model.

#4 Re: Code_Aster usage » What's wrong with "Saturated file"? » 2020-04-20 11:12:19

Hi,
Thank you very much. I will try it.
Actually, I would like to speed up the calculation.Do you have any suggestion?
Here is med file.

#6 Code_Aster usage » What's wrong with "Saturated file"? » 2020-04-20 10:42:19

Mahions
Replies: 6

Hi all,
I want to simulate a contact problem caused by the change of temperature field. The problem includes hundreds of contact pairs. To save times, I select part of the model for analysis. However, the cost is still high, more than three days. And it reaches the maximum number of record of the base GLOBALE.

!----------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <JEVEUX_42>                                                                            !
   !                                                                                                    !
   !  Saturated file, the maximum number of record 62914 of the base GLOBALE is reached                 !
   !  it is necessary to start again computation while passing a basic maximum size on the command line !
   !  argument" - max_base" followed by the Mo value.                                                   !
   !----------------------------------------------------------------------------------------------------!

1. What's the meaning of that? Just input " - max_base" by Arguments on astk version 2017.0?
2. How to perform complete problem quickly? This is very important.

Thanks for your suggestions.

#7 Salome-Meca usage » How to create steel wire ropes by SALOME-MECA? » 2020-03-26 02:52:37

Mahions
Replies: 1

Hi all,
   Do you know how to create wire ropes in GEOMETRY MODULE?

Thank you for your help.

#8 Re: Code_Aster usage » [SOLVED] Augmentez le temps CPU » 2019-12-12 02:15:18

The default maximum CPU time is 23 hours. How to increase it?

#9 Re: Code_Aster usage » How to import node temperature? » 2019-12-06 02:35:13

Thank you very much.
I think DEFI_NAPPE can solve my problem, but I can't find parameter 'radian'. How to define local cylindrical coordinate system? So ,temperature depends on two parameter 'Z’ and 'radian'(which is not in the parameter).

#10 Re: Code_Aster usage » How to import node temperature? » 2019-12-04 10:21:42

Hello,
Here is my model. The temperature field is as follows:
height(mm)(Z axis)     mark1(K)   mark2   mark3    ...      mark12
    -380                               **             **            **        ...         **     
     588                               **             **            **        ...         **
     ...                                                                 

How to import such a temperature field for thermal expansion?
Thank you very much.

#11 Code_Aster usage » How to import node temperature? » 2019-12-02 09:22:07

Mahions
Replies: 3

Hi,
I want to import many node temperature values for thermal expansion.Do you have any cases about that?
Thanks in advance for your help.

#12 Code_Aster usage » Does Code-Aster have a limit on quantity of grid? » 2019-12-02 07:03:20

Mahions
Replies: 0

Hi, all
Considering the computational efficiency, how many grids can Code-Aster calculate at most? Is there any way to speed up the solution?

#13 Re: Code_Aster usage » Il n'y a pas de mailles a projeter. » 2019-11-20 11:10:53

Thank you very much.
I watched the video and related document, still have several problems.

1. What is the difference between compound mesh in MESH module and assemble two meshes by ASSE_MAILLAGE in ASTER module?

2. LIAISON_SOLIDE is used to define rigid body.  Why use it to connect structure of element? In my opinion, the option COQ_POU in LIAISON_ELEM make it possible to connect shell with beam.

Looking forward to your reply.

#14 Re: Code_Aster usage » How to export node coordinate? » 2019-11-20 08:10:53

I learned a lot, but didn't find what I want .
Thank you!

#15 Re: Code_Aster usage » How to export node coordinate? » 2019-11-18 11:09:11

Thank you very much.
Actually, I want to create a temperature field. However, the temperature field is based on geometric model. I try to  match node coordinate with the temperature. Thus, how to define or import such a temperature field like this?
x   y    z      temp
*    *     *       *
*    *     *       *
.....

In addition, x,y,z must be the node coordinate?  Or just on the geometric model?

#16 Code_Aster usage » How to export node coordinate? » 2019-11-18 04:47:40

Mahions
Replies: 4

How to export node coordinate in text file form?

#17 Code_Aster usage » How to define a temperature field along the direction of z axis? » 2019-11-13 09:39:40

Mahions
Replies: 0

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.

#18 Code_Aster usage » Il n'y a pas de mailles a projeter. » 2019-11-06 06:59:19

Mahions
Replies: 3

Hi everyone,
I'm pretty new to code_aster and try to simulate a contact problem between two rods caused by thermal thermal. Before that,  I want to test if a rod could work by apply FORCE_FACE. Then, learn to define a temperature field and DEFI_CONTACT. 

When I take the first step ,something is wrong. And this is the error-related part in the message file:
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee. !
   ! Les bases sont fermees.                               !
   ! Type de l'exception : error                           !
   !                                                       !
   !  Il n'y a pas de mailles a projeter.                  !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

In my opinion, LIAISON_MAIL is not defined correctly . But I don't know how to modify.  Follows part of the .comm file.

load = AFFE_CHAR_MECA(
    DDL_IMPO=_F(GROUP_NO=('BC2', ), LIAISON='ENCASTRE'),
    FORCE_FACE=_F(FX=5.0, GROUP_MA=('CONTACT12', )),
    LIAISON_ELEM=_F(
        AXE_POUTRE=(0.0, 0.0, 1.0),
        CARA_ELEM=elemprop,
        GROUP_MA_1=('SHELLBEAM2', ),
        GROUP_NO_2=('BEAMSHELL2', ),
        OPTION='COQ_POU'
    ),
    LIAISON_MAIL=(
        _F(
            GROUP_MA_ESCL=('DOWN2', ),
            GROUP_MA_MAIT=('UP2', ),
            TYPE_RACCORD='MASSIF'
        ), _F(
            GROUP_MA_ESCL=('SHELL2', ),
            GROUP_MA_MAIT=('DOWN2', ),
            TYPE_RACCORD='MASSIF_COQUE'
        )
    ),
    MODELE=model
)

I would really appreciate if someone could help me with questions I have. I'm using Salome-Meca 2017.0.2.

#19 Salome-Meca usage » How to connect different meshes? » 2019-10-30 03:53:29

Mahions
Replies: 0

I have a complicated model which is simplified into beam, shell and solid.  The solid is divided into two parts(solid1、2), for I can’t mesh the solid into hexahedral grids at once.  It turns out that there are four meshes, mesh1(beam),mesh2(shell), mesh3(solid1),mesh4(solid2). When I read all meshes in AsterStudy , I assemble all the meshes by SUPERPOSE in ASSE_MAILLAGE, and use LIAISON_MAIL and LIAISON_ELEM to bind the DOF between mesh3\4 and mesh 1\2,2\3. Is that right?

#20 Salome-Meca usage » How to create a fuel rod wire winding? » 2019-09-23 04:29:37

Mahions
Replies: 0

Actually,I want to create a fuel rod with wire winding,  but I don't know how to create that geometrical object like spring. Also, I can import model from IGES or STEP. However, meshes are wrong. It seems to me that the imported model is a whole. And wire winding's size is much smaller than that of the fuel rod. Could you tell me how to create a fuel rod wire winding or how to mesh a fuel rod with wire winding?

#22 Code_Aster development » What's the meaning of function zi? » 2019-06-24 12:24:25

Mahions
Replies: 1

Hi, I don't understand the meaning of function zi in the source code. And I can not find the definition of the function zi and something like this , zk、zk8、zk32 and so on. Where can I find the definitions or documents of the functions?

#23 Code_Aster development » add a new law of behaviour under irradiation » 2019-06-14 08:58:28

Mahions
Replies: 0

Hello,guys
        Recently I have been tried to add a new non-linear constitutive model. But I don't know where to start. I read the documents (R5.03.23) , validation(ssnv230a),and the code(irrad3m,lc0030).I have several poblems.
    1. Thermoelastic,plasticity, creep under irradiation, swelling and thermal mechanism of creep must be taken into account in the models. R5.03.23 said, "If a thermal mechanism of creep would be added, time would play an explicit part of manner of course." What does this mean?
    2.In order to add such a new material ,what should I do?

  If you have any wonderful experience and suggestions ,please tell me. Thank you very much in advance for your help!

#24 Code_Aster installation » No such files or directory:'as_run':'as_run' » 2019-05-28 03:22:34

Mahions
Replies: 1

$ ./waf test -n sslp114a
checking environment... loading /home/zhlei/dev/codeaster/devtools/etc/env_unstable.sh
executing: ./waf.engine test -n sslp114a --out=build/std --jobs=4
Waf: Entering directory `/home/zhlei/dev/codeaster/src/build/std/release'
Traceback (most recent call last):
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Scripting.py", line 119, in waf_entry_point
    run_commands()
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Scripting.py", line 179, in run_commands
    ctx=run_command(cmd_name)
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Scripting.py", line 170, in run_command
    ctx.execute()
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Scripting.py", line 367, in execute
    return execute_method(self)
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Build.py", line 93, in execute
    self.execute_build()
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Build.py", line 96, in execute_build
    self.recurse([self.run_dir])
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Context.py", line 126, in recurse
    user_function(self)
  File "/home/zhlei/dev/codeaster/src/wscript", line 278, in runtest
    self.load('runtest', tooldir='waftools')
  File "/home/zhlei/dev/codeaster/src/.waf3-2.0.12-842c5bf2f993e0c1b0fa93868c080b69/waflib/Context.py", line 83, in load
fun(self)
  File "waftools/runtest.py", line 74, in runtest
    if not _has_asrun():
  File "waftools/runtest.py", line 128, in _has_asrun
    check_call(['as_run', '--version'], stdout=PIPE, stderr=PIPE)
  File "/usr/lib/python3.6/subprocess.py", line 286, in check_call
    retcode = call(*popenargs, **kwargs)
  File "/usr/lib/python3.6/subprocess.py", line 267, in call
    with Popen(*popenargs, **kwargs) as p:
  File "/usr/lib/python3.6/subprocess.py", line 709, in __init__
    restore_signals, start_new_session)
  File "/usr/lib/python3.6/subprocess.py", line 1344, in _execute_child
    raise child_exception_type(errno_num, err_msg, err_filename)
FileNotFoundError: [Errno 2] No such file or directory: 'as_run': 'as_run'
I compiled CA from source. But when I  ran the Code_Aster validation cases with waf, it didn't work.

#25 Code_Aster installation » could not find the python module 'numpy' » 2019-05-20 08:41:56

Mahions
Replies: 1

The complete log :
Checking for numpy
['/usr/bin/python', '-c', "\nimport numpy as current_module\nversion = getattr(current_module, '__version__', None)\nif version is not None:\n\tprint(str(version))\nelse:\n\tprint('unknown version')\n"]
err: Traceback (most recent call last):
  File "<string>", line 2, in <module>
ModuleNotFoundError: No module named 'numpy'
Error in sys.excepthook:
Traceback (most recent call last):
  File "/usr/lib/python3/dist-packages/apport_python_hook.py", line 63, in apport_excepthook
    from apport.fileutils import likely_packaged, get_recent_crashes
  File "/usr/lib/python3/dist-packages/apport/__init__.py", line 5, in <module>
    from apport.report import Report
  File "/usr/lib/python3/dist-packages/apport/report.py", line 30, in <module>
    import apport.fileutils
  File "/usr/lib/python3/dist-packages/apport/fileutils.py", line 23, in <module>
    from apport.packaging_impl import impl as packaging
  File "/usr/lib/python3/dist-packages/apport/packaging_impl.py", line 24, in <module>
    import apt
  File "/usr/lib/python3/dist-packages/apt/__init__.py", line 23, in <module>
    import apt_pkg
ModuleNotFoundError: No module named 'apt_pkg'

Original exception was:
Traceback (most recent call last):
  File "<string>", line 2, in <module>
ModuleNotFoundError: No module named 'numpy'

   What's wrong with this? I have already installed python-numpy.