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

You are not logged in.

#1 Re: Code_Aster usage » [SOLVED] adding nodes using MAIL_PY() » 2021-11-15 14:24:02

Dear Guillermo,

This helps a lot. I really appreciate your help.

Thank you so much.
Cheers,
Zaki

#2 Code_Aster usage » [SOLVED] adding nodes using MAIL_PY() » 2021-10-20 12:18:08

z-1616
Replies: 3

Dear all,

I'm trying to add nodes to my mesh using the following commands in MAIL_PY():

|  # Le décalage de -1 est déjà fait sur les numéros des noeuds et des mailles
|  # Pour ajouter des Noeuds
|     #  Coordonnées
|           LesNoeudsEnPlus = NP.array([[xx,yy,zz]])
|           LesNoeudsEnPlus = NP.concatenate((LesNoeudsEnPlus,NP.array([[xx,yy,zz]])))
|     #  Noms
|           NomNoeudsEnPlus =     ['X%' % in]
|           NomNoeudsEnPlus.append('X%' % in)
|  # Ajout des noeuds au maillage : coordonnées , noms
|     mm.cn = NP.concatenate((mm.cn,LesNoeudsEnPlus))
|     mm.correspondance_noeuds = tuple(linomno + NomNoeudsEnPlus )

but I don't understand what the part ['X%' % in] is can someone give me an example of how should I successfully modify it.
I would really appreciate any help in this matter.

Thank you in advance,
Best Regards,
Zaki

#3 Code_Aster usage » Accessing mesh module within the aster study » 2021-10-19 16:10:31

z-1616
Replies: 0

Dear all,

I'm trying to access the mesh module within the aster study to modify my mesh. It is important for me to do this using the aster study. Does anybody know how to open the module probably using a python command.
please let me know.

thanks in advance,
Zaki

#4 Re: Code_Aster usage » Structure Optimization with Code_Aster » 2021-03-02 12:54:36

Hello Everyone!

I'm really interested in this code but with the new versions of code_aster, I'm facing a lot of problems running it.
Did anyone know how to make this code runs on the new version of code_aster.

Thank you in advance,
Best Regards,
Zaki

#5 Re: Code_Aster development » Parallel simulation code_aster setup » 2020-11-04 19:10:54

Dear Sameer,

Again, thank you so much for your help.
I was able to solve some of the problems I was getting previously but everytime some new errors emerge.
Obviously, It is something hard to solve. Nevertheless, I appreciate your help.

Best Regards,
Zaki Alomar

#6 Re: Code_Aster development » Parallel simulation code_aster setup » 2020-11-03 17:34:03

Dear Sameer,

Thank you for the quick reply.
I'm still having problems when I'm trying to follow the procedures. In the second step, I am getting "no such file or directory" as seen below"
--------------------------------------------------------------------------------------------------------------------
zalomar@ZAKI:~/Downloads/aster-full-src-14.6.0$ sed -i "s:PREFER_COMPILER\ =\ 'GNU':PR                                                                                                                        EFER_COMPILER\ =\'GNU_without_MATH'\nMATHLIB=\ '/opt/OpenBLAS/lib/libopenblas.a':g" se                                                                                                                        tup.cfg
zalomar@ZAKI:~/Downloads/aster-full-src-14.6.0$  python3 setup.py install

--------------------------------------------------------------------------------
Code_Aster Setup version 14.6.0-1
Copyright (c) 2001-2020 EDF R&D - httpwww.code-aster.org
--------------------------------------------------------------------------------


--------------------------------------------------------------------------------
Command line :
  /usr/bin/python3 setup.py install
--------------------------------------------------------------------------------


--------------------------------------------------------------------------------
Reading config file '/home/students/zalomar/Downloads/aster-full-src-14.6.0/setup.cfg'                                                                                                                        ...
      ASTER_ROOT (from cfg) : '/opt/aster'
PREFER_COMPILER (from cfg) : 'GNU_without_MATH'
         MATHLIB (from cfg) : '/opt/OpenBLAS/lib/libopenblas.a'
        MAXDEPTH (from cfg) : 5
      USE_LOCATE (from cfg) : False
PREFER_SHARED_LIBS (from cfg) : True

--------------------------------------------------------------------------------
Installation on :
Ubuntu 18.04.5 LTS

Linux ZAKI 4.15.0-112-generic #113-Ubuntu SMP Thu Jul 9 23:41:39 UTC 2020 x86_64
--------------------------------------------------------------------------------


--------------------------------------------------------------------------------
Checking for max command length...   32768.0
Checking for file... /usr/bin/file
Checking for ar... /usr/bin/ar
Checking for architecture... Linux / posix / x86_64
Checking for number of processors (core)... 10 (will use: make -j 8)
Checking for Code_Aster platform type... LINUX64
Checking for bash... /bin/bash
Checking for Python version... 3.7.5
Checking for python3.7... no
Checking for libpython3.7.so... no
Checking for libpython3.7.a... no
Checking for python3.7m... no
Checking for libpython3.7m.so... no
Checking for libpython3.7m.a... no


Exception raised. See 'setup.log' file for details.
zalomar@ZAKI:~/Downloads/aster-full-src-14.6.0$ echo "$HOSTNAME cpu=$(cat /proc/cpuinf                                                                                                                        o | grep processor | wc -l)" > /opt/aster/etc/codeaster/mpi_hostfile
-bash: /opt/aster/etc/codeaster/mpi_hostfile: No such file or directory

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

Any idea why?

Best Regards,
Zaki.

#7 Code_Aster development » Parallel simulation code_aster setup » 2020-11-03 03:08:30

z-1616
Replies: 5

Hello Everyone,

I am trying to activate the parallel version of Salome. I found the following procedures listed in this link: "sites.google.com/site/codeastersalomemeca/home/code_asterno-heiretuka/parallel-code_aster-12-4-english"
Unfortunately, these steps were provided a long time ago and I think they cannot work on the recent version of Salome (2019).
Anyone can help in this matter, I would greatly appreciate it.

Thank you in advance,

Zaki Alomar.

#8 Re: Code_Aster usage » non-static analysis of cellular structure (contact problem) » 2020-07-21 17:15:23

Dear Sameer,

I have also tried the surfaces

"""contacT = DEFI_CONTACT(ALGO_RESO_GEOM='POINT_FIXE',
                       MODELE=model,
                       ZONE=_F(ALGO_CONT='CONTRAINTE',
                               GROUP_MA_ESCL=('S_face', ),
                               GROUP_MA_MAIT=('T_face', ),
                               SANS_GROUP_NO=('nTop', )))

contactB = DEFI_CONTACT(ALGO_RESO_GEOM='POINT_FIXE',
                        MODELE=model,
                        ZONE=_F(ALGO_CONT='CONTRAINTE',
                                GROUP_MA_ESCL=('S_face', ),
                                GROUP_MA_MAIT=('B_face', ),
                                SANS_GROUP_NO=('nBottom', )))"""

And I got this error:

   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception user raised but not interceptee.                                 !
   ! The bases are fermees.                                                         !
   ! Type of the exception: error                                                   !
   !                                                                                !
   !  La maille de peau : M31161 ne peut pas tre roriente.                       !
   !  Car elle est insre entre 2 mailles "support" places de part et d'autre     !
   ! : M70805 et M92612.                                                            !
   !                                                                                !
   ! Conseils :                                                                     !
   !  Vous pouvez utiliser les mots-cls GROUP_MA_VOL en 3D et GROUP_MA_SURF en 2D  !
   ! pour choisir une des                                                           !
   !                                                                                !
   ! deux mailles supports et ainsi choisir la normale permettant la rorientation. !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Any idea how to solve this error?
Thank you

#9 Code_Aster usage » non-static analysis of cellular structure (contact problem) » 2020-07-20 21:24:22

z-1616
Replies: 2

Dear all,

I'm running some preliminary compression simulations on a cellular structure. After a large deformation, the struts penetrate the top and bottom plates as shown in the picture attached. Therefore, I tried to define a contact between the struts and the plates but i got this error "

  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception user raised but not interceptee.                                 !
   ! The bases are fermees.                                                         !
   ! Type of the exception: error                                                   !
   !                                                                                !
   ! At least one of the meshes of contact which you defined is of dimension 3, but !
   ! the dimension of your problem is:  3.                                          !
   ! This mesh is                                                                   !
   ! thus not an edge mesh. There must be an error in your setting in data.         !
   !                                                                                !
   ! Advice:                                                                        !
   ! Check your AFFE_MODELE and the type of your                                    !
   ! meshes in the definition of contact surfaces.                                  !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

The simulation parameters are set as follows:

model = AFFE_MODELE(AFFE=_F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

sscurve = DEFI_FONCTION(INTERPOL=('LIN', ),
                        NOM_PARA='EPSI',
                        NOM_RESU='SIGMA',
                        PROL_DROITE='LINEAIRE',
                        PROL_GAUCHE='CONSTANT',
                        VALE=(0.019, 1270.0, 0.02, 1300.0, 0.021, 1361.0, 0.022, 1398.0, 0.023, 1431.0, 0.025, 1466.0, 0.028, 1534.0, 0.031, 1576.0, 0.037, 1620.0, 0.04, 1635.0, 0.045, 1652.0, 0.05, 1660.0, 0.06, 1670.0, 0.07, 1670.0, 0.08, 1670.0, 0.09, 1670.0))

mater = DEFI_MATERIAU(ELAS=_F(E=190000.0,
                              NU=0.3),
                      TRACTION=_F(SIGM=sscurve))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh)

Loadf = DEFI_FONCTION(INTERPOL=('LIN', ),
                      NOM_PARA='INST',
                      PROL_DROITE='LINEAIRE',
                      PROL_GAUCHE='CONSTANT',
                      VALE=(0.0, 0.0, 0.1, -0.018))

Stepr = DEFI_LIST_REEL(DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=0.1,
                                     PAS=0.001))

times = DEFI_LIST_INST(DEFI_LIST=_F(LIST_INST=Stepr),
                       ECHEC=_F(EVENEMENT='ERREUR',
                                SUBD_METHODE='MANUEL',
                                SUBD_NIVEAU=10,
                                SUBD_PAS=10),
                       METHODE='AUTO')

Fixation = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
                                      DY=0.0,
                                      DZ=0.0,
                                      GROUP_MA=('Bottom_plate', )),
                          MODELE=model)

load0 = AFFE_CHAR_MECA_F(DDL_IMPO=_F(DY=Loadf,
                                     GROUP_MA=('Top_plate', )),
                         MODELE=model)

contacT = DEFI_CONTACT(ALGO_RESO_GEOM='POINT_FIXE',
                       MODELE=model,
                       ZONE=_F(ALGO_CONT='CONTRAINTE',
                               GROUP_MA_ESCL=('Struts', ),
                               GROUP_MA_MAIT=('Top_plate', ),
                               SANS_GROUP_NO=('nTop', )))

contactB = DEFI_CONTACT(ALGO_RESO_GEOM='POINT_FIXE',
                        MODELE=model,
                        ZONE=_F(ALGO_CONT='CONTRAINTE',
                                GROUP_MA_ESCL=('Struts', ),
                                GROUP_MA_MAIT=('Bottom_plate', ),
                                SANS_GROUP_NO=('nBottom', )))

resnonl = STAT_NON_LINE(CHAM_MATER=fieldmat,
                        COMPORTEMENT=_F(DEFORMATION='SIMO_MIEHE',
                                        ITER_INTE_MAXI=20,
                                        RELATION='VMIS_ISOT_TRAC',
                                        RESI_INTE_RELA=1e-06,
                                        TOUT='OUI'),
                        CONVERGENCE=_F(ITER_GLOB_MAXI=100,
                                       RESI_GLOB_MAXI=1e-06,
                                       RESI_GLOB_RELA=1e-06),
                        EXCIT=(_F(CHARGE=Fixation),
                               _F(CHARGE=load0)),
                        INCREMENT=_F(LIST_INST=times),
                        METHODE='NEWTON',
                        MODELE=model,
                        NEWTON=_F(MATRICE='TANGENTE',
                                  REAC_INCR=1,
                                  REAC_ITER=1))

Your help would be greatly appreciated.
Thank you in advance.