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

You are not logged in.

Pages: 1

#1 Re: Code_Aster installation » problem installing code Aster version 14.6 » 2022-08-18 10:43:40

ftrillaudp

Dear all,

The problem persists in the testing 15.2 version as well in ubuntu 22.04LTS.

Hopefully, it gets resolved in a near future.

Best regards,

Frederic

#2 Re: Code_Aster installation » installation Ubuntu 20.04LTS - missing libboost-python solved » 2022-07-30 10:09:21

ftrillaudp

Dear all,

I was off the forum for some time now. Indeed, I should have mentioned that the fix depends on the python version installed. You need first to know which library is installed to chose the proper name.

I have not moved to Ubuntu 22.04LTS, hopefully this issue was fixed by the developers.

Best,

Frederic

#3 Re: Code_Aster usage » Chaining code_saturne and code_aster - LIRE_RESU and 'EVOL_CHAR' » 2022-07-30 10:03:58

ftrillaudp

Dear Shayor,

I am very sorry to reply after so many months considering you were under some deadline as far as I understood. I have not used Code_Aster for some time and did not bother to look at this post thinking it was closed. Only by curiosity, I went back to it having discussed this case study with a former student of mine.

I have to get back to it to refresh my mind. If you still need some help, I may be able to provide at least on setting up a case study.

I did this example to understand how to simulate fluid-structure interaction. The case study is random. However, it allowed me to play with the material properties to see the range of values permissible until the coupling fails. If the deformation is relatively small, it works very well.

Best,

Frederic

#4 Code_Aster installation » installation Ubuntu 20.04LTS - missing libboost-python solved » 2020-12-23 02:06:38

ftrillaudp
Replies: 3

Dear all,

Installation on Ubuntu 20.04LTS failed missing libboost_python with default python install version 3.8.x.

Here is the fixed that I used:

in the setup.py file line 758, change:

cfg.get("LIBNAME_BOOST", ["boost_python3", "boost_python3-mt"]),

for

cfg.get("LIBNAME_BOOST", ["boost_python38", "boost_python3-mt"]),

You can also remove "boost_python3-mt".

Best,

Frederic

#5 Re: Code_Aster usage » Chaining code_saturne and code_aster - LIRE_RESU and 'EVOL_CHAR' » 2020-04-12 22:07:16

ftrillaudp

Dear all,

I have a working example of coupling code_saturne v6 and code_aster v14.4. I am joining the full folder including the code_saturne case coupled to the code_aster case.

The meshes for the fluid and for the solid surrounding the fluid are done in Salome platform and exported in med format.

The python main.py must be run to chain the CFD and structural computations.

The case study is setup through fictitious material properties so that there is enough deformation to impact the fluid but not enough to lead to large deformations. Large deformations may require a different mechanical solver as well as recompute the meshes (regularization as elements with negative volume may arise due to gherkin's large deformation).

Best regards,

Frederic

#6 Re: Code_Aster usage » Chaining code_saturne and code_aster - LIRE_RESU and 'EVOL_CHAR' » 2020-03-19 03:59:53

ftrillaudp

Dear AsterO'dactyle and Johannes_ACKVA,

Thanks to take the time to look at my issue. My knowledge of Code_aster is quite limited and the explaining is much appreciated.

My goal is to couple code_saturne and code_Aster so that one can evaluate the deformation of the walls on the fluid behavior. Some way to go still. If I manage, I will post the model on the forum.

Thanks again, best regards

Frederic

#7 Re: Code_Aster usage » Chaining code_saturne and code_aster - LIRE_RESU and 'EVOL_CHAR' » 2020-03-16 21:54:09

ftrillaudp

Hi,

It took some time but I figured out how to save only the pressure field on the surface of the fluid mesh. However, it did not solve the issue. The same error keeps popping up:

Erreur lors de la fabrication d'un champ par éléments :                    !
   !  Il manque la composante : PRES  sur la maille : M349                      !
   !  Champ : ch_pres .002.000000     !

I have no idea how to get the name of the component from the *.med file.

I have joined the new surface mesh with the pressure field.

Best,

Frederic

#8 Code_Aster usage » Chaining code_saturne and code_aster - LIRE_RESU and 'EVOL_CHAR' » 2020-03-16 03:28:30

ftrillaudp
Replies: 9

Dear all,

I am stuck trying too read a result from code_saturne in med format. Code_Aster does not recognizes the component 'PRES':

<A> <CALCULEL_58>                                                          !
   !                                                                            !
   ! Erreur lors de la fabrication d'un champ par éléments :                    !
   !  Il manque la composante : PRES  sur la maille : M349                      !
   !  Champ : ch_pres .002.000000                                               !
   !                                                                            !
   ! Risques et conseils :                                                      !
   !  Si cette erreur se produit lors de l'exécution de la commande PROJ_CHAMP, !
   !  il est possible de poursuivre le calcul en choisissant PROL_ZERO='OUI'

Here is the *.comm and attached the case study for someone to run it a give it a try:

DEBUT(PAR_LOT='OUI',);

import math
import numpy
                         
############
### MESH ###
############

mail_mec = LIRE_MAILLAGE(UNITE=20, FORMAT='MED',);
mail_cfd = LIRE_MAILLAGE(UNITE=21, FORMAT='MED',);


################
### MATERIAL ###
################

acier = DEFI_MATERIAU(ELAS = _F(E=190E9, NU=0.3,),);

mat1 = AFFE_MATERIAU(
                     MAILLAGE = mail_mec,
                     AFFE = (
                              _F(GROUP_MA='SSOL001', MATER=acier,),
                            ),
                    );
                         
##############
### MODELE ###
##############

mod_cfd = AFFE_MODELE(
                      MAILLAGE = mail_cfd,
                      AFFE = _F(TOUT='OUI', PHENOMENE='MECANIQUE', MODELISATION='3D',),
                      );

mod_mec = AFFE_MODELE(
                    MAILLAGE = mail_mec,
                    AFFE = _F(TOUT='OUI', PHENOMENE='MECANIQUE', MODELISATION='3D',),
                   );
                     
##################
### CHARGEMENT ###
##################

pres_cfd = LIRE_RESU(
                          TYPE_RESU = 'EVOL_CHAR',
                          FORMAT = 'MED',
                          MODELE = mod_cfd,
                          UNITE = 21,
                          FORMAT_MED = _F(
                                          NOM_CHAM = 'PRES',
                                          NOM_CMP = 'PRES',
                                          NOM_CHAM_MED = 'Pressure',
                                          NOM_CMP_MED = 'Scalar',
                                        ),
                          TOUT_ORDRE = 'OUI',
                         );

ch_pres = PROJ_CHAMP(
                      METHODE = 'COLLOCATION',
                      RESULTAT = pres_cfd,
                      MODELE_1 = mod_cfd,
                      MODELE_2 = mod_mec,
                      TOUT_ORDRE = 'OUI',
                      VIS_A_VIS = _F(
                                     TOUT_1 = 'OUI',
                                     GROUP_MA_2 = 'SFAC001',
                                     ),
                    );

chg_mec = AFFE_CHAR_MECA(
                            MODELE = mod_mec,
                            EVOL_CHAR = ch_pres,
                            DDL_IMPO = _F(
                                          GROUP_MA = ('SFAC002','SFAC003'),
                                          DX=0, DY=0, DZ=0
                                          ),
                            );
                       
meca = MECA_STATIQUE(
                     MODELE = mod_mec,
                     CHAM_MATER = mat1,
                     EXCIT = (
                              _F(CHARGE = chg_mec),
                             ),
                    );

meca = CALC_CHAMP(
                  reuse = meca,
                  RESULTAT = meca,
                  CONTRAINTE = 'SIGM_ELNO',
                  FORCE = 'REAC_NODA',
                  CRITERES = ('SIEQ_ELNO', 'SIEQ_NOEU',),
                );

## Ne fonctionne pas avec AFFE_CHAR_MECA_F
IMPR_CONCEPT(
              UNITE = 81,
              FORMAT = 'MED',
              VERSION_MED = '4.0.0',
              CONCEPT = (_F(CHARGE = chg_mec),),
            );
           
IMPR_RESU(
          UNITE = 80,
          FORMAT = 'MED',
          VERSION_MED = '4.0.0',
          RESU= (
                  _F(
                    GROUP_MA = 'SSOL001',
                    RESULTAT = meca,
                    NOM_CHAM = 'DEPL',
                  ),
                  _F(
                      GROUP_MA = 'SSOL001',
                      RESULTAT = meca,
                      NOM_CHAM = 'SIEQ_NOEU',
                      NOM_CMP = ('VMIS',),
                    ),
                ),
        );

FIN();


Best,

Frederic

#9 Re: Code_Aster usage » import volumetric force from text files for loading » 2020-03-11 16:28:37

ftrillaudp

Dear all,

I figure out how to import a volumetric force from an external file. Here is the process:

Format of external file (space delimited):

X1 Y1 Z1 FX1 FY1 FZ1
X2 Y2 Z2 FX2 FY2 FZ2
.
.
.
Xn Yn Zn FXn FYn FZn

n nodes in the mesh.

Here is the annoted *.comm to read in the file and perform mechanical computation. We tested on a case suited with homogeneous and inhomogeneous force field.


### Code_Aster 14.4 ###
### Ubuntu 18.04LTS ###

### Two bodies: SSOL001 and SSOL002. One face: SFAC003 belonging to SSOL002 on which nodal fixed displacements are imposed to 0.0

DEBUT(PAR_LOT='OUI',);

import math
import numpy
                         
############
### MESH ###
############

mail = LIRE_MAILLAGE(UNITE=20, FORMAT='MED',);

mail = DEFI_GROUP(
                  reuse = mail,
                  MAILLAGE = mail,
                  CREA_GROUP_NO = (
                                                     _F(NOM='SFAC003_NO', GROUP_MA = 'SFAC003', CRIT_NOEUD = 'TOUS'),
                                                   ),
             );


################
### MATERIAL ###
################

### Arbitrary values:
acier = DEFI_MATERIAU(ELAS = _F(E=128E9, NU=0.3,),);
cuivre = DEFI_MATERIAU(ELAS = _F(E=128E9, NU=0.3,),);

mat1 = AFFE_MATERIAU(
                     MAILLAGE = mail,
                     AFFE = (
                              _F(GROUP_MA='SSOL001', MATER=cuivre,),
                              _F(GROUP_MA='SSOL002', MATER=acier,),
                            ),
                    );
                         
##############
### MODELE ###
##############

mod1 = AFFE_MODELE(
                    MAILLAGE = mail,
                    AFFE = _F(TOUT='OUI', PHENOMENE='MECANIQUE', MODELISATION='3D',),
                   );
                     
##################
### CHARGEMENT ###
##################

### Read in external forces, definition of function per component of vector force. We get the closest node in the mesh corresponding to the (X,Y,Z) position in the external file. The force will be imposed on the body SSOL001.
### replace fullpathToFile by the absolute path to the external fiel, here named "nodalForces.txt"
### We need numpy to use "genfromtxt" and deal with arrays.

def func_FX(x, y, z):
  distance, FX = list(), list()
  array = numpy.genfromtxt("/fullpathTopFile/nodalForces.txt", delimiter=None)
  for ii in range(0,array.shape[0]):
    lineArray = array[ii][:]
    XD = lineArray[0]; YD = lineArray[1]; ZD = lineArray[2]
    distance.append((x-XD)**2+(y-YD)**2+(z-ZD)**2)
    FX.append(lineArray[3])
  return FX[distance.index(min(distance))]
 
def func_FY(x, y, z):
  distance, FY = list(), list()
  array = numpy.genfromtxt("/fullpathTopFile/nodalForces.txt", delimiter=None)
  for ii in range(0,array.shape[0]):
    lineArray = array[ii][:]
    XD = lineArray[0]; YD = lineArray[1]; ZD = lineArray[2]
    distance.append((x-XD)**2+(y-YD)**2+(z-ZD)**2)
    FY.append(lineArray[4])
  return FY[distance.index(min(distance))]

def func_FZ(x, y, z):
  distance, FZ = list(), list()
  array = numpy.genfromtxt("/fullpathTopFile/nodalForces.txt", delimiter=None)
  for ii in range(0,array.shape[0]):
    lineArray = array[ii][:]
    XD = lineArray[0]; YD = lineArray[1]; ZD = lineArray[2]
    distance.append((x-XD)**2+(y-YD)**2+(z-ZD)**2)
    FZ.append(lineArray[5])
  return FZ[distance.index(min(distance))]

### Construction of the function:
forceX = FORMULE(VALE='func_FX(X,Y,Z)', NOM_PARA=('X','Y','Z',), func_FX=func_FX,);
forceY = FORMULE(VALE='func_FY(X,Y,Z)', NOM_PARA=('X','Y','Z',), func_FY=func_FY,);
forceZ = FORMULE(VALE='func_FZ(X,Y,Z)', NOM_PARA=('X','Y','Z',), func_FZ=func_FZ,);

forceVol = AFFE_CHAR_MECA_F(
                            MODELE=mod1,
                            FORCE_INTERNE = _F(GROUP_MA='SSOL001', FX=forceX, FY=forceY, FZ=forceZ,),
                            INFO=2,
                          );
                       

deplImpo = AFFE_CHAR_MECA(
                          MODELE = mod1,
                          DDL_IMPO = _F(GROUP_MA = 'SFAC003', DX=0, DY=0, DZ=0),
                         );
                       
meca = MECA_STATIQUE(
                     MODELE = mod1,
                     CHAM_MATER = mat1,
                     EXCIT = (
                              _F(CHARGE = forceVol),
                              _F(CHARGE = deplImpo),
                             ),
                    );

meca = CALC_CHAMP(
                  reuse = meca,
                  RESULTAT = meca,
                  CONTRAINTE = 'SIGM_ELNO',
                  FORCE = 'REAC_NODA',
                  CRITERES = ('SIEQ_ELNO', 'SIEQ_NOEU',),
                );

## It does not work with  AFFE_CHAR_MECA_F
IMPR_CONCEPT(
              UNITE = 81,
              FORMAT = 'MED',
              VERSION_MED = '4.0.0',
              CONCEPT = (_F(CHARGE = forceVol),),
            );
           
IMPR_RESU(
          UNITE = 80,
          FORMAT = 'MED',
          VERSION_MED = '4.0.0',
          RESU= (
                  _F(
                    GROUP_MA = ('SSOL001','SSOL002',),
                    RESULTAT = meca,
                    NOM_CHAM = 'DEPL',
                  ),
                  _F(
                      GROUP_MA = ('SSOL001','SSOL002',),
                      RESULTAT = meca,
                      NOM_CHAM = 'SIEQ_NOEU',
                      NOM_CMP = ('VMIS',),
                    ),
                ),
        );

FIN();

In the attached file, we provide a results showing the displacement field for a pure vertical Fz force (Fx = Fy = 0).

Best regards,

Frederic

#10 Code_Aster usage » import volumetric force from text files for loading » 2020-01-23 20:34:25

ftrillaudp
Replies: 2

Dear all,

I am stuck trying to import Lorentz force (volumetric force) from an external electromagnetic solver as a load for a mechanical computation in Code_aster. I have individual files for the x-, y- and z-component of the force (force-FX.text, force-FY.txt, force-FZ.txt). The files have the following structure (force-FX.txt):

X1 Y1 Z1 FX1
X2 Y2 Z2 FX2
X3 Y3 Z3 FX1
.....

I would like to import the forces in the loading:

AFFE_CHAR_MECA_F(MODELE=Model, FORCE_NODAL=_F(GROUP_MA='Nodes', FX=FX, FY=FY, FZ=FZ,),);

How could it be done?

Best regards,

Frederic

#11 Re: Salome-Meca installation » [SOLVED] Salome-meca 2019.0.0-1 » 2019-10-01 16:54:24

ftrillaudp

I did all the recommended suggestions. Here is what I got as error:

runSalome running on debian-gnu-linux-vm
Searching for a free port for naming service: 2819 - OK
Searching Naming Service  + found in 0.1 seconds
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by SALOME_Session_Server)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSalomeSession.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSalomeApp.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libLightApp.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libEvent.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libCAM.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libstd.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libLogWindow.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libGraphicsView.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libImageComposer.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libPyViewer.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libPyEditor.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libQxScene.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSUITApp.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSalomeStyle.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libToolsGUI.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSVTK.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libVTKViewer.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libPyConsole.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libGLViewer.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSPlot2d.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libPlot2d.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Qwt-612/lib/libqwt.so.6)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libSOCC.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libOCCViewer.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libViewerTools.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libPVViewer.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libsuit.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libObjBrowser.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libqtx.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /home/parallels/../../opt/salome_meca2019/appli_V2019_univ/lib/salome/libPyInterp.so)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Xml.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5OpenGL.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5PrintSupport.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Test.so.5)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libvtkpqApplicationComponents-pv5.6.so.1)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libvtkpqComponents-pv5.6.so.1)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libvtkpqPython-pv5.6.so.1)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libvtkpqCore-pv5.6.so.1)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libvtkpqWidgets-pv5.6.so.1)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Network.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Help.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Sql.so.5)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libQtTesting-pv5.6.so)
SALOME_Session_Server: /usr/lib/x86_64-linux-gnu/libQt5Core.so.5: version `Qt_5.9' not found (required by /opt/salome_meca2019/V2019_univ/prerequisites/Paraview-v560p1/lib/paraview-5.6/libvtkGUISupportQt-pv5.6.so.1)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Widgets.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Gui.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Core.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Svg.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5Concurrent.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `GLIBCXX_3.4.21' not found (required by /usr/lib/x86_64-linux-gnu/libQt5CLucene.so.5)
SALOME_Session_Server: /opt/salome_meca2019/V2019_univ/prerequisites/debianForSalome/lib/libstdc++.so.6: version `CXXABI_1.3.9' not found (required by /usr/lib/x86_64-linux-gnu/libQt5CLucene.so.5)
Traceback (most recent call last):
  File "/opt/salome_meca2019/V2019_univ/modules/KERNEL_V9_3_0/bin/salome/orbmodule.py", line 172, in waitNSPID
    os.kill(thePID,0)
ProcessLookupError: [Errno 3] No such process

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/opt/salome_meca2019/appli_V2019_univ/bin/salome/runSalome.py", line 633, in useSalome
    clt = startSalome(args, modules_list, modules_root_dir)
  File "/opt/salome_meca2019/appli_V2019_univ/bin/salome/runSalome.py", line 578, in startSalome
    session=clt.waitNSPID("/Kernel/Session",mySessionServ.PID,SALOME.Session)
  File "/opt/salome_meca2019/V2019_univ/modules/KERNEL_V9_3_0/bin/salome/orbmodule.py", line 174, in waitNSPID
    raise RuntimeError("Process %d for %s not found" % (thePID,theName))
RuntimeError: Process 18951 for /Kernel/Session not found
Searching /Kernel/Session in Naming Service 

--- Error during Salome launch ---

#12 Re: Salome-Meca installation » [SOLVED] Salome-meca 2019.0.0-1 » 2019-10-01 15:53:51

ftrillaudp

IN my case, it did not work. Here is the output:

runSalome running on debian-gnu-linux-vm
Searching for a free port for naming service: 2815 - OK
Searching Naming Service  + found in 0.1 seconds
Searching /Kernel/Session in Naming Service  ++Traceback (most recent call last):
  File "/opt/salome_meca/V2019_univ/modules/KERNEL_V9_3_0/bin/salome/orbmodule.py", line 172, in waitNSPID
    os.kill(thePID,0)
ProcessLookupError: [Errno 3] No such process

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/opt/salome_meca/appli_V2019_univ/bin/salome/runSalome.py", line 633, in useSalome
    clt = startSalome(args, modules_list, modules_root_dir)
  File "/opt/salome_meca/appli_V2019_univ/bin/salome/runSalome.py", line 578, in startSalome
    session=clt.waitNSPID("/Kernel/Session",mySessionServ.PID,SALOME.Session)
  File "/opt/salome_meca/V2019_univ/modules/KERNEL_V9_3_0/bin/salome/orbmodule.py", line 174, in waitNSPID
    raise RuntimeError("Process %d for %s not found" % (thePID,theName))
RuntimeError: Process 5864 for /Kernel/Session not found


--- Error during Salome launch ---

I have added:
export LD_LIBRARY_PATH=/opt/salome_meca/V2019_univ/prerequisites/Qt-591/lib:$LD_LIBRARY_PATH

I have installed net-tools and libnlopt0

Best,

Frederic

#13 Re: Salome-Meca installation » [SOLVED] Salome-meca 2019.0.0-1 » 2019-09-27 19:41:27

ftrillaudp

Same error here on Debian 9.

runSalome running on debian-gnu-linux-vm
Error when trying to access active network connections.
Traceback (most recent call last):
  File "/opt/salome_meca/appli_V2019_univ/bin/salome/PortManager.py", line 111, in __isNetworkConnectionActiveOnPort
    p = Popen(cmd, stdout=PIPE, stderr=STDOUT)
  File "/opt/salome_meca/V2019_univ/prerequisites/Python-365/lib/python3.6/subprocess.py", line 709, in __init__
    restore_signals, start_new_session)
  File "/opt/salome_meca/V2019_univ/prerequisites/Python-365/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: 'netstat': 'netstat'
Searching for a free port for naming service: 2811 - OK
Searching Naming Service  + found in 0.1 seconds
SALOME_Session_Server: error while loading shared libraries: libQt5Xml.so.5: cannot open shared object file: No such file or directory
Traceback (most recent call last):
  File "/opt/salome_meca/V2019_univ/modules/KERNEL_V9_3_0/bin/salome/orbmodule.py", line 172, in waitNSPID
    os.kill(thePID,0)
ProcessLookupError: [Errno 3] No such process

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/opt/salome_meca/appli_V2019_univ/bin/salome/runSalome.py", line 633, in useSalome
    clt = startSalome(args, modules_list, modules_root_dir)
  File "/opt/salome_meca/appli_V2019_univ/bin/salome/runSalome.py", line 578, in startSalome
    session=clt.waitNSPID("/Kernel/Session",mySessionServ.PID,SALOME.Session)
  File "/opt/salome_meca/V2019_univ/modules/KERNEL_V9_3_0/bin/salome/orbmodule.py", line 174, in waitNSPID
    raise RuntimeError("Process %d for %s not found" % (thePID,theName))
RuntimeError: Process 8757 for /Kernel/Session not found
Searching /Kernel/Session in Naming Service 

--- Error during Salome launch ---

#14 Re: Code_Aster usage » warning "type K16" during thermal resolution - code_aster » 2016-05-06 22:48:56

ftrillaudp

Dear all,

Thanks a lot for the explanations.

Best,

Frederic

#15 Code_Aster usage » warning "type K16" during thermal resolution - code_aster » 2016-04-27 22:12:14

ftrillaudp
Replies: 3

Dear all,

During the resolution of a transient thermal problem, a warning appears which I do not understand. Here is the *.hdf file generated in Salome-Meca as well as the solution and the code_aster *.comm files.

The warning is the following:


   !-----------------------------------------------------------------!
   ! <A> <PREPOST_91>                                                !
   !                                                                 !
   !  On ne sait pas imprimer au format 'MED' les cartes de type K16 !
   !                                                                 !
   !  On ne sait pas imprimer au format 'MED' les cartes de type K16 !
   !                                                                 !
   !  On ne sait pas imprimer au format 'MED' les cartes de type K16 !
   !                                                                 !
   !  On ne sait pas imprimer au format 'MED' les cartes de type K16 !
   !                                                                 !
   !  On ne sait pas imprimer au format 'MED' les cartes de type K16 !
   !                                                                 !
   !                                                                 !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette  !
   ! alarme, vous pouvez obtenir des résultats inattendus !          !
   !-----------------------------------------------------------------!

If someone understands the warning, please do share. It is not the first time that I am experiencing this warning "type K16".

Best,

Frederic

#16 Code_Aster usage » Thermal interface between two solids-howto in Code_aster » 2013-10-18 19:17:29

ftrillaudp
Replies: 1

Dear all,

I am wondering how it is possible to define the quality of the thermal interface between two solids. In some codes such as FEMAP TMG and Thermal Desktop, a transfer coefficient in W/m²K can be defined to account for the thermal quality of the interface without having to create any mesh. What would be the best implementation in Code_Aster to account for the thermal interface?

Best,

Frederic

#17 Introduce yourself / Présentez vous » From Mexico with regards » 2013-07-10 19:58:23

ftrillaudp
Replies: 0

Never late than ever. Hi, I am a research professor in the Institute of Engineering of the National Autonomous University of Mexico in Mexico city.

I came from an experimental background even though I am doing a lot more numerical work those days. My fields of interest are electromagnetism, thermodynamics and mechanics applied to superconducting or/and cryogenic systems for power and scientific applications. lately, I have been wandering in space technologies and electromagnetic compatibility. I am into open source software and use various known codes in addition to Code_Aster such as Elmer-CSC, Gmsh/Getdp, SYRTHES, Code_Saturne and Cast3M (under academic license). I would not say that I am a proficient user but I get around and the supporting community is very often responded n a timely manner.

Best of luck with Code_Aster,

Frederic Trillaud, PhD

#18 Re: Code_Aster usage » How to use FORMULE and print out results » 2013-04-20 23:29:22

ftrillaudp

Thanks, I am waiting for the next realease of code_aster then.

Best,

Frederic

#19 Re: Code_Aster usage » How to use FORMULE and print out results » 2013-04-18 16:12:16

ftrillaudp

HI,

I do have any error message (see the *.mess file below). Here is the mesh.
I am using Salome-Meca 2013-1 64 bits. I get in the postprocessing the flux and the temperature but I do not see the user result.


Thanks,


Frederic


Message:

# ------------------------------------------------------------------------------------------
# Impression du contenu du fichier de commandes à exécuter :
# ------------------------------------------------------------------------------------------
#############
### ASTER ###
#############

### NOM DE FICHIER : solenoid2D.comm
### AUTEUR : Frederic Trillaud <ftrillaudp@gmail.com>
### INSTITUTION : IINGEN-UNAM
### DATE (Anglais) : 04/13/2013
### UNITE : SI


### INTRODUCTION : calculs magnetique lineaire se basant sur une analogie thermique.


DEBUT();


### Parametres geometriques ###
### Ruban et isolation :
largRub = 4e-3;
epaisRub = 100e-6;
largIso = largRub;
epaisIso = 50e-6;
### Bobine :
nTours = 10;
nBob = 1;
epSepBob = 50e-6;
epaisBob = nBob*largRub+(nBob-1)*epSepBob;
largBob = (epaisRub+epaisIso)*nTours;
rIntBob = 0.1;
sTranBob = largBob*epaisBob;
factEnc = (nBob*nTours*epaisRub*largRub)/sTranBob;

### Parametres electriques ###
courElec = 100.;
je = courElec/(factEnc*sTranBob);

### Reluctivite dans le vide (equivalent a la conductivite thermique) ###
k0 = 1e7/(4*pi);


magMESH = LIRE_MAILLAGE(
                        UNITE = 20,
                        FORMAT = 'MED',
                       );

magMATER = DEFI_MATERIAU(
                         THER = _F(
                                   LAMBDA = k0,
                                  ),
                        );

magMODE = AFFE_MODELE(
                     MAILLAGE = magMESH,
                     AFFE=_F(
                             TOUT = 'OUI',
                             PHENOMENE = 'THERMIQUE',
                             MODELISATION = 'AXIS',
                            ),
                    );

matFIELD = AFFE_MATERIAU(
                         MAILLAGE = magMESH,
                         AFFE = _F(
                                   TOUT = 'OUI',
                                   MATER = magMATER,
                                  ),
                        );

magLOAD = AFFE_CHAR_THER(
                         MODELE = magMODE,
                         TEMP_IMPO = (
                                      _F(
                                         GROUP_MA = 'WIRE001',
                                         TEMP = 0.0,
                                        ),
                                     ),
                         SOURCE = (
                                   _F(
                                      GROUP_MA = 'COND001',
                                      SOUR = je,
                                     ),
                                  ),
                        );

potVec = THER_LINEAIRE(
                       MODELE = magMODE,
                       CHAM_MATER = matFIELD,
                       EXCIT = _F(
                                  CHARGE = magLOAD,
                                 ),
                      );


BMod = FORMULE(
               NOM_PARA = (
                           "FLUX",
                           "FLUY",
                          ),
               VALE = """ sqrt(FLUX**2+FLUY**2)/k0 """
              );


potVec = CALC_CHAMP(
                    reuse = potVec,
                    RESULTAT = potVec,
                    MODELE = magMODE,
                    CHAM_MATER = matFIELD,
                    THERMIQUE = 'FLUX_ELGA',
                    CHAM_UTIL = _F(
                                   NOM_CHAM = 'FLUX_ELGA',
                                   FORMULE = BMod,
                                   NUME_CHAM_RESU = 1,
                                  ),
                 );

                                 
IMPR_RESU(
          FORMAT = 'MED',
          RESU = (
                  _F(
                     RESULTAT = potVec,
                     NOM_CHAM = 'TEMP',
                    ),
                  _F(
                     RESULTAT = potVec,
                     NOM_CHAM = 'FLUX_ELGA',
                    ),
                  _F(
                     RESULTAT = potVec,
                     NOM_CHAM = 'UT01_ELGA',
                    ),
                 ),
          );
         

FIN();

# ------------------------------------------------------------------------------------------
                                                                               
                                                                               
        -- CODE_ASTER -- VERSION : DÉVELOPPEMENT STABILISÉE (testing) --       
                                                                               
                          Version 11.3.0 du 12/12/2012                         
                         Copyright EDF R&D 1991 - 2013                         
                                                                               
                    Exécution du : Thu Apr 18 10:16:27 2013                     
                      Nom de la machine : ELC121P000A095M                       
                              Architecture : 64bit                             
                          Type de processeur : x86_64                           
                Système d'exploitation : Linux 3.2.0-41-generic                 
                        Langue des messages : en (UTF-8)                       
                                                                               
                           Parallélisme MPI : inactif                           
                          Parallélisme OpenMP : actif                           
                        Nombre de processus utilisés : 2                       
                      Version de la librairie HDF5 : 1.8.8                     
                      Version de la librairie MED : 3.0.6                       
                          Librairie MUMPS : installée                           
                    Version de la librairie SCOTCH : 5.1.10                     

                  Mémoire limite pour l'exécution : 512.00 Mo                   
                   consommée par l'initialisation : 209.95 Mo                   
                  par les objets du jeu de commandes : 0.00 Mo                 
                 reste pour l'allocation dynamique : 302.05 Mo                 
                Taille limite des fichiers d'échange : 48.00 Go                 
                                                                               


  # ------------------------------------------------------------------------------------------
  # Commande No :  0001            Concept de type : -
  # ------------------------------------------------------------------------------------------
  DEBUT(INFO=1,
        PAR_LOT='OUI',
        IMPR_MACRO='NON',
        MEMOIRE=_F(TAILLE_BLOC=800.0,
                   TAILLE_GROUP_ELEM=1000,),
        RESERVE_CPU=_F(BORNE=900,),
        MESURE_TEMPS=_F(MOYENNE='OUI',
                        NIVE_DETAIL=1,),
        )


  Valeur initiale du temps CPU maximum =   120 secondes
  Valeur du temps CPU maximum passé aux commandes =   108 secondes
  Réserve CPU prévue = 12 secondes

  Début de lecture...


     Réouverture de la base

     Nom de la base                          :  ELEMBASE
     Créée avec la version                   :  11.03.00
     Nombre d'enregistrements utilisés       :  24
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre d'identificateurs utilisés       :  123
     Taille maximum du répertoire            :  300
     Pourcentage d'utilisation du répertoire :  41 %


     Fermeture de la base

     Nom de la base                          :  ELEMBASE
     Nombre d'enregistrements utilisés       :  24
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  49
     Volume des accès en lecture             :         38.28 Mo.
     Nombre total d'accès en écriture        :  3
     Volume des accès en écriture            :          2.34 Mo.
     Nombre d'identificateurs utilisés       :  123
     Taille maximum du répertoire            :  300
     Pourcentage d'utilisation du répertoire :  41 %

  Relecture des catalogues des éléments faite.

  Fin de lecture (durée  0.020000  s.)
  # Mémoire (Mo) :   212.52 /   209.96 /    27.51 /    14.95 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0001   user+syst:        0.06s (syst:        0.02s, elaps:        0.07s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0002            Concept de type : maillage_sdaster
  # ------------------------------------------------------------------------------------------
  magMESH = LIRE_MAILLAGE(FORMAT='MED',
                          UNITE=20,
                          INFO_MED=1,
                          INFO=1,
                          VERI_MAIL=_F(APLAT=1.E-3,
                                       VERIF='OUI',),
                          )

  ====== VERIFICATION DU MAILLAGE ======


------------ MAILLAGE magMESH  - IMPRESSIONS NIVEAU  1 ------------

JE-18-AVRI-2013 10:16:27                                                       

NOMBRE DE NOEUDS                        7763

NOMBRE DE MAILLES                      15712
                              SEG2                  515
                              TRIA3               15197

NOMBRE DE GROUPES DE MAILLES               5
                              WIRE002             22
                              COND001             50
                              COND002          10992
                              COND003           4155
                              WIRE001            327

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

  # Mémoire (Mo) :   212.52 /   209.96 /    27.51 /    14.95 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0002   user+syst:        0.05s (syst:        0.00s, elaps:        0.14s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0003            Concept de type : mater_sdaster
  # ------------------------------------------------------------------------------------------
  magMATER = DEFI_MATERIAU(THER=_F(LAMBDA=7.9577471545947669E5,),
                           INFO=1,)

  La mémoire consommée actuellement hors JEVEUX est de 185.33 Mo.
  La limite de l'allocation dynamique JEVEUX est fixée à 325.08 Mo.

MATERIAU : magMATER
RELATION DE COMPORTEMENT: THER           


  # Mémoire (Mo) :   212.52 /   209.96 /    27.51 /    14.95 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0003   user+syst:        0.01s (syst:        0.00s, elaps:        0.01s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0004            Concept de type : modele_sdaster
  # ------------------------------------------------------------------------------------------
  magMODE = AFFE_MODELE(MAILLAGE=magMESH,
                        AFFE=_F(PHENOMENE='THERMIQUE',
                                TOUT='OUI',
                                MODELISATION=('AXIS', ),),
                        INFO=1,
                        VERI_JACOBIEN='OUI',
                        PARTITION=_F(PARALLELISME='GROUP_ELEM',),
                        )

  La mémoire consommée actuellement hors JEVEUX est de 185.32 Mo.
  La limite de l'allocation dynamique JEVEUX est fixée à 326.67 Mo.

SUR LES        15712 MAILLES DU MAILLAGE magMESH
    ON A DEMANDE L'AFFECTATION DE        15712
    ON A PU EN AFFECTER                  15712

    MODELISATION      ELEMENT FINI      TYPE MAILLE          NOMBRE
    AXIS              THAXSE2           SEG2                    515
    AXIS              THAXTR3           TRIA3                 15197

#2      Calculs elementaires et assemblages               CPU (USER+SYST/SYST/ELAPS):      0.02      0.00      0.05
  # Mémoire (Mo) :   212.52 /   210.74 /    27.51 /    23.62 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0004   user+syst:        0.09s (syst:        0.01s, elaps:        0.21s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0005            Concept de type : cham_mater
  # ------------------------------------------------------------------------------------------
  matFIELD = AFFE_MATERIAU(MAILLAGE=magMESH,
                           AFFE=_F(MATER=(magMATER, ),
                                   TOUT='OUI',),
                           INFO=1,
                           )

  La mémoire consommée actuellement hors JEVEUX est de 185.32 Mo.
  La limite de l'allocation dynamique JEVEUX est fixée à 326.68 Mo.
  # Mémoire (Mo) :   212.52 /   210.74 /    27.51 /    23.62 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0005   user+syst:        0.00s (syst:        0.00s, elaps:        0.01s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0006            Concept de type : char_ther
  # ------------------------------------------------------------------------------------------
  magLOAD = AFFE_CHAR_THER(MODELE=magMODE,
                           SOURCE=_F(GROUP_MA=('COND001', ),
                                     SOUR=2.5000000E7,),
                           TEMP_IMPO=_F(GROUP_MA=('WIRE001', ),
                                        TEMP=0.0,),
                           INFO=1,
                           )

  # Mémoire (Mo) :   212.52 /   210.74 /    27.51 /    23.62 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0006   user+syst:        0.03s (syst:        0.00s, elaps:        0.08s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0007            Concept de type : evol_ther
  # ------------------------------------------------------------------------------------------
  potVec = THER_LINEAIRE(MODELE=magMODE,
                         CHAM_MATER=matFIELD,
                         EXCIT=_F(CHARGE=magLOAD,),
                         PARM_THETA=0.56999999999999995,
                         ARCHIVAGE=_F(PRECISION=9.9999999999999995E-07,
                                      CRITERE='RELATIF',),
                         SOLVEUR=_F(RENUM='METIS',
                                    STOP_SINGULIER='OUI',
                                    NPREC=8,
                                    METHODE='MULT_FRONT',),
                         )


Le système linéaire à résoudre contient 8417 noeuds dont:
   - 7763 noeuds portant des degrés de liberté physiques
   - 654 noeuds portant des degrés de liberté de Lagrange
Pour un total de 8417 équations.

La matrice est de taille 8417 équations.
Elle contient 32357 termes non nuls si elle est symétrique et 56297 termes non nuls si elle n'est pas symétrique (le nombre de
termes non nuls est susceptible de varier si l'on utilise le contact en formulation continue ou la méthode XFEM avec contact).
Soit un taux de remplissage de   0.079 %.

Il n'y a pas d'état initial défini. On prend un état initial nul.

  Archivage des champs
    Champ stocké <TEMP> à l'instant  0.000000000000e+00 pour le numéro d'ordre 0
NUMERO D'ORDRE:      0 INSTANT:   0.000E+00 DUREE MOYENNE:   2.100E-01
#1      Resolution des systemes lineaires                 CPU (USER+SYST/SYST/ELAPS):      0.10      0.00      0.12
#2      Calculs elementaires et assemblages               CPU (USER+SYST/SYST/ELAPS):      0.10      0.00      0.13
  # Mémoire (Mo) :   471.01 /   343.01 /    32.18 /    23.62 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0007   user+syst:        0.24s (syst:        0.00s, elaps:        0.28s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0008            Concept de type : formule
  # ------------------------------------------------------------------------------------------
  BMod = FORMULE(VALE=' sqrt(FLUX**2+FLUY**2)/k0 ',
                 NOM_PARA=('FLUX', 'FLUY'),)

  # Mémoire (Mo) :   471.01 /   343.01 /    32.18 /    23.62 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0008   user+syst:        0.00s (syst:        0.00s, elaps:        0.00s)
  # ------------------------------------------------------------------------------------------

  # ------------------------------------------------------------------------------------------
  # Commande No :  0009            Concept de type : evol_ther
  # ------------------------------------------------------------------------------------------
  potVec = CALC_CHAMP(reuse=potVec,
                      CHAM_MATER=matFIELD,
                      RESULTAT=potVec,
                      MODELE=magMODE,
                      CHAM_UTIL=_F(FORMULE=(BMod, ),
                                   NOM_CHAM='FLUX_ELGA',
                                   NUME_CHAM_RESU=1,),
                      PRECISION=9.9999999999999995E-07,
                      THERMIQUE=('FLUX_ELGA', ),
                      INFO=1,
                      CRITERE='RELATIF',)


    Numéro d'ordre    0 : 15197 mailles ont été affectées (15712 mailles dans le maillage).
   
   !-------------------------------------------------------!
   ! <EXCEPTION> <DVP_1>                                   !
   !                                                       !
   !    Erreur de programmation : condition non respectée. !
   !                                                       !
   !                                                       !
   !                                                       !
   ! Il y a probablement une erreur dans la programmation. !
   ! Veuillez contacter votre assistance technique.        !
   !-------------------------------------------------------!
   

  Validation du concept 'potVec'.

  # ------------------------------------------------------------------------------------------
  # Commande No :  0010            Concept de type : -
  # ------------------------------------------------------------------------------------------
  FIN(RETASSAGE='NON',
      INFO_RESU='OUI',
      FORMAT_HDF='NON',
      UNITE=6,)

   
   !-------------------------------------------------------------------------!
   ! <I> <CATAMESS_89>                                                       !
   !                                                                         !
   !     Liste des alarmes émises lors de l'exécution du calcul.             !
   !                                                                         !
   !     Les alarmes que vous avez choisies d'ignorer sont précédées de (*). !
   !     Nombre d'occurrences pour chacune des alarmes :                     !
   !            aucune alarme                                                !
   !-------------------------------------------------------------------------!
   

======>

STRUCTURE DU CONCEPT potVec   CALCULE POUR 1 NUMERO D'ORDRE

LISTE DES NOMS SYMBOLIQUES:
! ---------- !----------------!----------------!
! NUME_ORDRE !      TEMP      !   FLUX_ELGA    !
! ---------- !----------------!----------------!
!          0 !     TEMP_R     !     FLUX_R     !
! ---------- !----------------!----------------!

LISTE DES NOMS DE VARIABLES D'ACCES:
                                      INST             DE TYPE  R

LISTE DES NOMS DE PARAMETRES:
! ---------- !----------------!----------------!----------------!----------------!----------------!----------------!----------------!----------------!
! NUME_ORDRE !     MODELE     !    CHAMPMAT    !    CARAELEM    !     EXCIT      !   ITER_GLOB    !   ITER_LINE    !   PARM_THETA   !     DELTAT     !
! ---------- !----------------!----------------!----------------!----------------!----------------!----------------!----------------!----------------!
!          0 !       K8       !       K8       !       K8       !      K24       !       I        !       I        !       R        !       R        !
! ---------- !----------------!----------------!----------------!----------------!----------------!----------------!----------------!----------------!
---------------------------------------------------------------------------
Concepts de la base: G
    Nom       Type                 Taille (Mo)         Nombre      Nombre de
                                                     d'objets       segments
    TOTAL                                19.89            150            278
 
    magMESH    MAILLAGE_SDASTER           1.59             34             52
    magMATER   MATER_SDASTER              0.00              4              4
    magMODE    MODELE_SDASTER             0.60              8             13
    matFIELD   CHAM_MATER                 0.00              5             10
    magLOAD    CHAR_THER                  0.16             24             49
    potVec     EVOL_THER                  2.09             28             42
    BMod       FORMULE                    0.00              2              2
    &FOZERO                               0.00              2              2
    &&_NUM_C                              0.00              1              1
    &CATA.AC                              0.00              2              4
    &CATA.CL                              0.23              1              3
    &CATA.GD                              0.04              4             11
    &CATA.ME                              0.13              2              4
    &CATA.OP                              0.27              4             19
    &CATA.PH                              0.00              1              1
    &CATA.PR                              0.00              2              4
    &CATA.TE                             14.74             17             42
    &CATA.TH                              0.01              2              4
    &CATA.TM                              0.02              7             11
---------------------------------------------------------------------------


     Fermeture de la base

     Nom de la base                          :  GLOBALE
     Nombre d'enregistrements utilisés       :  34
     Nombre d'enregistrements maximum        :  62914
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  7
     Volume des accès en lecture             :          5.47 Mo.
     Nombre total d'accès en écriture        :  48
     Volume des accès en écriture            :         37.50 Mo.
     Nombre d'identificateurs utilisés       :  293
     Taille maximum du répertoire            :  2000
     Pourcentage d'utilisation du répertoire :  14 %


     Fermeture de la base

     Nom de la base                          :  VOLATILE
     Nombre d'enregistrements utilisés       :  3
     Nombre d'enregistrements maximum        :  62914
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  0
     Volume des accès en lecture             :          0.00 Mo.
     Nombre total d'accès en écriture        :  4
     Volume des accès en écriture            :          3.12 Mo.
     Nombre d'identificateurs utilisés       :  572
     Taille maximum du répertoire            :  2000
     Pourcentage d'utilisation du répertoire :  28 %

<FIN> Arrêt normal dans "FIN".
 
<I>       FERMETURE DES BASES EFFECTUEE
 
   STATISTIQUES CONCERNANT L'ALLOCATION DYNAMIQUE :
     TAILLE CUMULEE MAXIMUM            :                    32  Mo.
     TAILLE CUMULEE LIBEREE            :                    12  Mo.
     NOMBRE TOTAL D'ALLOCATIONS        :                 46179
     NOMBRE TOTAL DE LIBERATIONS       :                 46179
     APPELS AU MECANISME DE LIBERATION :                     0
     TAILLE MEMOIRE CUMULEE RECUPEREE  :                     0  Mo.
     VOLUME DES LECTURES               :                     0  Mo.
     VOLUME DES ECRITURES              :                     0  Mo.
 
   MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION :      24.09 Mo
     - IMPOSE DE NOMBREUX ACCES DISQUE
     - RALENTIT LA VITESSE D'EXECUTION
   MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION :      32.18 Mo
     - LIMITE LES ACCES DISQUE
     - AMELIORE LA VITESSE D'EXECUTION
   MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS     :     471.01 Mo
     - COMPREND LA MEMOIRE CONSOMMEE PAR  JEVEUX,
       LE SUPERVISEUR PYTHON, LES LIBRAIRIES EXTERNES
 
<I>       FIN D'EXECUTION LE : JE-18-AVRI-2013 10:16:29


  <I> Informations sur les temps d'exécution
      Temps cpu total                         2.03 s
      Temps cpu user total                    1.94 s
      Temps cpu systeme total                 0.09 s
      Temps cpu restant                     105.97 s


  La signature de la base sauvegardée est (à l'adresse 0) :
    738c07e67d5e30bebdaade547a232edfdd6771f4

********************************************************************************
* COMMAND                  :       USER :     SYSTEM :   USER+SYS :    ELAPSED *
********************************************************************************
* init (jdc)               :       0.54 :       0.04 :       0.58 :       0.60 *
*  . compile               :       0.00 :       0.00 :       0.00 :       0.00 *
*  . exec_compile          :       0.08 :       0.01 :       0.09 :       0.08 *
*  . report                :       0.01 :       0.00 :       0.01 :       0.00 *
*  . build                 :       0.00 :       0.00 :       0.00 :       0.00 *
* DEBUT                    :       0.04 :       0.02 :       0.06 :       0.07 *
* LIRE_MAILLAGE            :       0.05 :       0.00 :       0.05 :       0.14 *
* DEFI_MATERIAU            :       0.01 :       0.00 :       0.01 :       0.01 *
* AFFE_MODELE              :       0.08 :       0.01 :       0.09 :       0.21 *
* AFFE_MATERIAU            :       0.00 :       0.00 :       0.00 :       0.01 *
* AFFE_CHAR_THER           :       0.03 :       0.00 :       0.03 :       0.08 *
* THER_LINEAIRE            :       0.24 :       0.00 :       0.24 :       0.28 *
* FORMULE                  :       0.00 :       0.00 :       0.00 :       0.00 *
* CALC_CHAMP               :       0.95 :       0.02 :       0.97 :       0.83 *
* FIN                      :       0.02 :       0.02 :       0.04 :       0.05 *
*  . part Superviseur      :       0.58 :       0.06 :       0.64 :       0.79 *
*  . part Fortran          :       1.36 :       0.03 :       1.39 :       1.55 *
********************************************************************************
* TOTAL_JOB                :       1.94 :       0.09 :       2.03 :       2.34 *
********************************************************************************

JDC.py : ERREUR A L'EXECUTION - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR d'execution de JDC en MIXTE
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee. !
   ! Les bases sont fermees.                               !
   ! Type de l'exception : error                           !
   !                                                       !
   !    Erreur de programmation : condition non respectée. !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE

>> JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_0001-6793-ELC121P000A095M=1

#20 Code_Aster usage » How to use FORMULE and print out results » 2013-04-15 17:06:24

ftrillaudp
Replies: 5

Hi,

I am trying out FORMULE. The idea is to do a magnetostatic computation relying on an analogy with thermal problem. I have computed the potential vector as the field temperature, now I am trying to get the module of the magnetic field using the computed thermal field (FLUX_ELGA). However, I cannot recover the user field (UT01_ELGA), here is the code:

Soeone has an idea of what I am doing wrong?

DEBUT();


### Parametres geometriques ###
### Ruban et isolation :
largRub = 4e-3;
epaisRub = 100e-6;
largIso = largRub;
epaisIso = 50e-6;
### Bobine :
nTours = 10;
nBob = 1;
epSepBob = 50e-6;
epaisBob = nBob*largRub+(nBob-1)*epSepBob;
largBob = (epaisRub+epaisIso)*nTours;
rIntBob = 0.1;
sTranBob = largBob*epaisBob;
factEnc = (nBob*nTours*epaisRub*largRub)/sTranBob;

### Parametres electriques ###
courElec = 100.;
je = courElec/(factEnc*sTranBob);

### Reluctivite dans le vide (equivalent a la conductivite thermique) ###
k0 = 1e7/(4*pi);


magMESH = LIRE_MAILLAGE(
                        UNITE = 20,
                        FORMAT = 'MED',
                       );

magMATER = DEFI_MATERIAU(
                         THER = _F(
                                   LAMBDA = k0,
                                  ),
                        );

magMODE = AFFE_MODELE(
                     MAILLAGE = magMESH,
                     AFFE=_F(
                             TOUT = 'OUI',
                             PHENOMENE = 'THERMIQUE',
                             MODELISATION = 'AXIS',
                            ),
                    );

matFIELD = AFFE_MATERIAU(
                         MAILLAGE = magMESH,
                         AFFE = _F(
                                   TOUT = 'OUI',
                                   MATER = magMATER,
                                  ),
                        );

magLOAD = AFFE_CHAR_THER(
                         MODELE = magMODE,
                         TEMP_IMPO = (
                                      _F(
                                         GROUP_MA = 'WIRE001',
                                         TEMP = 0.0,
                                        ),
                                     ),
                         SOURCE = (
                                   _F(
                                      GROUP_MA = 'COND001',
                                      SOUR = je,
                                     ),
                                  ),
                        );

TEMP = THER_LINEAIRE(
                       MODELE = magMODE,
                       CHAM_MATER = matFIELD,
                       EXCIT = _F(
                                  CHARGE = magLOAD, 
                                 ),
                      );


BMod = FORMULE(
               NOM_PARA = (
                           "FLUX",
                           "FLUY",
                          ),
               VALE = """ sqrt(FLUX**2+FLUY**2)/k0 """
              );


TEMP = CALC_CHAMP(
                 reuse = TEMP,
                 RESULTAT = TEMP,
                 MODELE = magMODE,
                 CHAM_MATER = matFIELD,
                 THERMIQUE = 'FLUX_ELGA',
                 CHAM_UTIL = _F(
                                NOM_CHAM = 'FLUX_ELGA',
                                FORMULE = BMod,
                                NUME_CHAM_RESU = 1,
                               ),
                 );

                                 
IMPR_RESU(
          FORMAT = 'MED',
          RESU = (
                  _F(
                     RESULTAT = TEMP,
                     NOM_CHAM = 'TEMP',
                    ),
                  _F(
                     RESULTAT = TEMP,
                     NOM_CHAM = 'FLUX_ELGA',
                    ),
                  _F(
                     RESULTAT = TEMP,
                     NOM_CHAM = 'UT01_ELGA',
                    ),
                 ),
          );
         

FIN();

#21 Code_Aster usage » heat flux calculation - error message CALC_ELEM » 2013-04-13 22:34:50

ftrillaudp
Replies: 21

Hi,

I am trying a simple heat model. I can compute the temperature, but I am stuck while computing derived entities such as the heat flux. An erro message pops up:

erreur de syntaxe,  Erreur de nom: name 'CALC_ELEM' is not defined ligne 55

Best,

Here is the code:

DEBUT();


magMESH = LIRE_MAILLAGE(
                        UNITE = 20,
                        FORMAT = 'MED',
                       );

magMATER = DEFI_MATERIAU(
                         THER = _F(
                                   LAMBDA = 100,
                                  ),
                        );

magMODE = AFFE_MODELE(
                     MAILLAGE = magMESH,
                     AFFE=_F(
                             TOUT = 'OUI',
                             PHENOMENE = 'THERMIQUE',
                             MODELISATION = 'AXIS',
                            ),
                    );

matFIELD = AFFE_MATERIAU(
                         MAILLAGE = magMESH,
                         AFFE = _F(
                                   TOUT = 'OUI',
                                   MATER = magMATER,
                                  ),
                        );

magLOAD = AFFE_CHAR_THER(
                         MODELE = magMODE,
                         TEMP_IMPO = (
                                      _F(
                                         GROUP_MA = 'WIRE001',
                                         TEMP = 0.0,
                                        ),
                                     ),
                         SOURCE = (
                                   _F(
                                      GROUP_MA = 'COND001',
                                      SOUR = 2500000,
                                     ),
                                  ),
                        );

TEMP = THER_LINEAIRE(
                       MODELE = magMODE,
                       CHAM_MATER = matFIELD,
                       EXCIT = _F(
                                  CHARGE = magLOAD,
                                 ),
                      );

TEMP = CALC_ELEM(
                 REUSE = TEMP,
                 RESULTAT = TEMP,
                 OPTION = (
                           'FLUX_ELGA',
                           'FLUX_ELNO_TEMP',
                          ),
                );

TEMP = CALC_NO(
               reuse = TEMP,
               RESULTAT = TEMP,
               OPTION = 'FLUX_NOEU_TEMP',
              );
                                 
IMPR_RESU(
          FORMAT = 'MED',
          RESU = _F(
                    RESULTAT = TEMP,
                    NOM_CHAM = (
                                'FLUX_ELGA_TEMP',
                                'FLUX_ELNO_TEMP',
                               ),
                   ),
          );
         

FIN();

#22 Code_Aster usage » calcule du volume d'un solide-Measure of voluminal solid » 2012-05-11 21:59:00

ftrillaudp
Replies: 1

Bonjour,

Existe-t'il une fonction pour reccuperer le volume de solid pour definir les termes sources volumiques?
Is there a function in Code_Aster to calculate the volume of a solid to be used in the definition of a volumic sourca?

Cordialement,

Frederic

#23 Re: Code_Aster usage » issues with radiant heat implementation (thermique non lineaire) » 2012-04-26 02:20:41

ftrillaudp

Merci pour la reponse, le code tourne!

Frederic

#24 Code_Aster usage » issues with radiant heat implementation (thermique non lineaire) » 2012-04-24 22:21:23

ftrillaudp
Replies: 2

Bonjour,

J'essaie de resoudre un probleme de thermo-mecanique avec rayonnement. Le code genere une erreur. Voici le fichier *.comm and *.mess.
Je souhaiteria recupere le champ de temperature au dernier pas de temps pour le calcul mecanique, comment?

Cordialement,

Frederic Trillaud
UNAM

#25 Introduce yourself / Présentez vous » Introduction » 2012-04-24 03:30:57

ftrillaudp
Replies: 0

Good evening,

I am a French research professor at the National Autonomous University of Mexico (UNAM) working on thermal an mechanical analysis of superconducting systems and recently on space technologies.
I have been using Cast3M for quite sometime and Elmer. I am looking for new free tools.

Regards,

Frederic Trillaud, PhD
Institute of Engineering, UNAM

Pages: 1