Nastran to Code_Aster: modal frequency response

In this post we will have a look at a modal frequency response analysis. We will find the frequency response of the structure under a pressure load and a nodal force with a phase lag. A modal damping is also applied.

This case is taken from a publicy available example on the internet, but the original source is probably this one, for all due credit.
To read this article one will need some knowledge both in Nastran and in Code_Aster.

prob6_cover

The translation between Nastran and Code_Aster can be made by using the following correspondances :

Nastran keyword Code_Aster keyword
BEGIN DEBUT
SOL 111 ASSEMBLAGE + CALC_MODES(SOLVEUR_MODAL =_F(METHODE=’TRI_DIAG’),) + RECU_TABLE + EXTR_TABLE + PROJ_BASE + DYNA_VIBRA(BASE_CALCUL=’GENE’,TYPE_CALCUL=’HARM’,AMOR_MODAL) + REST_GENE_PHYS
TIME 600 See “prob24.export” file, line “A tpmax 600
SET, OUTPUT (XYPLOT) RECU_FONCTION + CALC_FONCTION + INFO_FONCTION + CALC_TABLE + IMPR_TABLE
DISPLACEMENT(PHASE, PLOT), XTGRID, YTLOG, … RECU_FONCTION + CALC_FONCTION + IMPR_FONCTION(FORMAT=’XMGRACE’,PILOTE = ‘PNG’)
EIGRL CALC_MODES(OPTION=’BANDE’,CALC_FREQ=_F(FREQ=(10.0,2000.0,),)
TABDMP1 DEFI_FONCTION(NOM_PARA=’FREQ’,NOM_RESU=’AMOR’) + CALC_FONC_INTERP + CREA_TABLE + EXTR_TABLE
TABLED1 DEFI_FONCTION(NOM_PARA=’FREQ’)
PLOAD2 AFFE_CHAR_MECA(FORCE_COQUE)
DPHASE,RLOAD2 DYNA_VIBRA(EXCIT,PHAS_DEG)
DLOAD CALC_VECT_ELEMASSE_VECTEUR
DAREA AFFE_CHAR_MECA(FORCE_NODALE)
FREQ1 DEFI_LIST_FREQ(DEBUT,INTERVALLE)
FREQ4 DEFI_LIST_FREQ(RAFFINEMENT)
PSHELL AFFE_MODELE(MODELISATION=’DKT‘), AFFE_CARA_ELEM(COQUE_F(EPAIS)) + AFFE_MATERIAU
CQUAD4 See “prob24.mail” file: QUAD4 + AFFE_CARA_ELEM(COQUE=_F(GROUPMA))
PARAM,WTMASS, MAT1 DEFI_MATERIAU(ELAS)
GRID See “prob24.mail” file: COOR_3D
SPC1 AFFE_CHAR_MECA(DDL_IMPO)
ENDDATA FIN

After the calculation we can find the frequency response for a given node and compare it with the Nastran results :

n55nastran

Here is the corresponding response from Code_Aster results :

f55zm

First of all, let’s see the Nastran code (file “prob24.dat”) for this example :


ID SEMINAR, PROB24
SOL 111
TIME 30
CEND
TITLE = FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS
SUBTITLE = USING THE MODAL METHOD WITH LANCZOS
ECHO = UNSORTED
SEALL = ALL
SPC = 1
SET 111 = 11, 33, 55
DISPLACEMENT(PHASE, PLOT) = 111
METHOD = 100
FREQUENCY = 100
SDAMPING = 100
SUBCASE 1
DLOAD = 100
LOADSET = 100
$
OUTPUT (XYPLOT)
$
XTGRID= YES
YTGRID= YES
XBGRID= YES
YBGRID= YES
YTLOG= YES
YBLOG= NO
XTITLE= FREQUENCY (HZ)
YTTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER, MAGNITUDE
YBTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER, PHASE
XYPLOT DISP RESPONSE / 11 (T3RM, T3IP)
YTTITLE= DISPLACEMENT RESPONSE AT TIP CENTER, MAGNITUDE
YBTITLE= DISPLACEMENT RESPONSE AT TIP CENTER, PHASE
XYPLOT DISP RESPONSE / 33 (T3RM, T3IP)
YTTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, MAGNITUDE
YBTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, PHASE
XYPLOT DISP RESPONSE / 55 (T3RM, T3IP)
$
BEGIN BULK
$
$ PARAMETERS FOR POST-PROCESSING
PARAM,COUPMASS,1
PARAM,WTMASS,0.00259
$
$ PLATE MODEL DESCRIBED IN NORMAL MODES EXAMPLE
$
INCLUDE 'plate.bdf'
$
$ EIGENVALUE EXTRACTION PARAMETERS
$
EIGRL, 100, 10., 2000.
$
$ SPECIFY MODAL DAMPING
$
TABDMP1, 100, CRIT,
+, 0., .03, 10., .03, ENDT
$ 
$ APPLY UNIT PRESSURE LOAD TO PLATE
$
LSEQ, 100, 300, 400
$
PLOAD2, 400, 1., 1, THRU, 40
$
$ APPLY PRESSURE LOAD
$ 
RLOAD2, 400, 300, , ,310
$
TABLED1, 310,,
, 10., 1., 1000., 1., ENDT
$
$ POINT LOAD
$ 
$ IF 'DAREA' CARDS ARE REFERENCED, THEN 
$ 'DPHASE' AND 'DELAY' CAN BE USED
$
RLOAD2, 500, 600, , 320, 310
$
DPHASE, 320, 11, 3, -45.
$
$
DAREA, 600, 11, 3, 1.0
$
$ COMBINE LOADS
$
DLOAD, 100, 1., .1, 400, 1.0, 500
$
$ SPECIFY FREQUENCY STEPS
$
FREQ1, 100, 20., 20., 49
FREQ4, 100, 20., 1000., .03, 5
$
ENDDATA 

The mesh is in another file “plate.bdf” as indicated by the INCLUDE command.

mesh

Here is the content for this file:


$
$   geometric input file for plate model
$
PSHELL   1       1       .1      1               1
CQUAD4   1       1       1       2       13      12
CQUAD4   2       1       2       3       14      13
CQUAD4   3       1       3       4       15      14
CQUAD4   4       1       4       5       16      15
CQUAD4   5       1       5       6       17      16
CQUAD4   6       1       6       7       18      17
CQUAD4   7       1       7       8       19      18
CQUAD4   8       1       8       9       20      19
CQUAD4   9       1       9       10      21      20
CQUAD4   10      1       10      11      22      21
CQUAD4   11      1       12      13      24      23
CQUAD4   12      1       13      14      25      24
CQUAD4   13      1       14      15      26      25
CQUAD4   14      1       15      16      27      26
CQUAD4   15      1       16      17      28      27
CQUAD4   16      1       17      18      29      28
CQUAD4   17      1       18      19      30      29
CQUAD4   18      1       19      20      31      30
CQUAD4   19      1       20      21      32      31
CQUAD4   20      1       21      22      33      32
CQUAD4   21      1       23      24      35      34
CQUAD4   22      1       24      25      36      35
CQUAD4   23      1       25      26      37      36
CQUAD4   24      1       26      27      38      37
CQUAD4   25      1       27      28      39      38
CQUAD4   26      1       28      29      40      39
CQUAD4   27      1       29      30      41      40
CQUAD4   28      1       30      31      42      41
CQUAD4   29      1       31      32      43      42
CQUAD4   30      1       32      33      44      43
CQUAD4   31      1       34      35      46      45
CQUAD4   32      1       35      36      47      46
CQUAD4   33      1       36      37      48      47
CQUAD4   34      1       37      38      49      48
CQUAD4   35      1       38      39      50      49
CQUAD4   36      1       39      40      51      50
CQUAD4   37      1       40      41      52      51
CQUAD4   38      1       41      42      53      52
CQUAD4   39      1       42      43      54      53
CQUAD4   40      1       43      44      55      54
$ 
MAT1     1       3.+7            .3      .282
$ 
GRID     1               0.      0.      0.
GRID     2               .5      0.      0.
GRID     3               1.      0.      0.
GRID     4               1.5     0.      0.
GRID     5               2.      0.      0.
GRID     6               2.5     0.      0.
GRID     7               3.      0.      0.
GRID     8               3.5     0.      0.
GRID     9               4.      0.      0.
GRID     10              4.5     0.      0.
GRID     11              5.      0.      0.
GRID     12              0.      .5      0.
GRID     13              .5      .5      0.
GRID     14              1.      .5      0.
GRID     15              1.5     .5      0.
GRID     16              2.      .5      0.
GRID     17              2.5     .5      0.
GRID     18              3.      .5      0.
GRID     19              3.5     .5      0.
GRID     20              4.      .5      0.
GRID     21              4.5     .5      0.
GRID     22              5.      .5      0.
GRID     23              0.      1.      0.
GRID     24              .5      1.      0.
GRID     25              1.      1.      0.
GRID     26              1.5     1.      0.
GRID     27              2.      1.      0.
GRID     28              2.5     1.      0.
GRID     29              3.      1.      0.
GRID     30              3.5     1.      0.
GRID     31              4.      1.      0.
GRID     32              4.5     1.      0.
GRID     33              5.      1.      0.
GRID     34              0.      1.5     0.
GRID     35              .5      1.5     0.
GRID     36              1.      1.5     0.
GRID     37              1.5     1.5     0.
GRID     38              2.      1.5     0.
GRID     39              2.5     1.5     0.
GRID     40              3.      1.5     0.
GRID     41              3.5     1.5     0.
GRID     42              4.      1.5     0.
GRID     43              4.5     1.5     0.
GRID     44              5.      1.5     0.
GRID     45              0.      2.      0.
GRID     46              .5      2.      0.
GRID     47              1.      2.      0.
GRID     48              1.5     2.      0.
GRID     49              2.      2.      0.
GRID     50              2.5     2.      0.
GRID     51              3.      2.      0.
GRID     52              3.5     2.      0.
GRID     53              4.      2.      0.
GRID     54              4.5     2.      0.
GRID     55              5.      2.      0.
$
SPC1     1       12345   1       12      23      34      45

In Code_Aster, you will need to create distinct files for the mesh and the commands. Let’s start with the main mesh file (file “prob24.mail”) :



TITRE    
	prob24
	FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS
FINSF

COOR_3D

	N1  0.0000000000e+00 0.0000000000e+00 0.0000000000e+00
	N2  5.0000000000e-01 0.0000000000e+00 0.0000000000e+00
	N13  5.0000000000e-01 5.0000000000e-01 0.0000000000e+00
	N12  0.0000000000e+00 5.0000000000e-01 0.0000000000e+00
	N3  1.0000000000e+00 0.0000000000e+00 0.0000000000e+00
	N14  1.0000000000e+00 5.0000000000e-01 0.0000000000e+00
	N4  1.5000000000e+00 0.0000000000e+00 0.0000000000e+00
	N15  1.5000000000e+00 5.0000000000e-01 0.0000000000e+00
	N5  2.0000000000e+00 0.0000000000e+00 0.0000000000e+00
	N16  2.0000000000e+00 5.0000000000e-01 0.0000000000e+00
	N6  2.5000000000e+00 0.0000000000e+00 0.0000000000e+00
	N17  2.5000000000e+00 5.0000000000e-01 0.0000000000e+00
	N7  3.0000000000e+00 0.0000000000e+00 0.0000000000e+00
	N18  3.0000000000e+00 5.0000000000e-01 0.0000000000e+00
	N8  3.5000000000e+00 0.0000000000e+00 0.0000000000e+00
	N19  3.5000000000e+00 5.0000000000e-01 0.0000000000e+00
	N9  4.0000000000e+00 0.0000000000e+00 0.0000000000e+00
	N20  4.0000000000e+00 5.0000000000e-01 0.0000000000e+00
	N10  4.5000000000e+00 0.0000000000e+00 0.0000000000e+00
	N21  4.5000000000e+00 5.0000000000e-01 0.0000000000e+00
	N11  5.0000000000e+00 0.0000000000e+00 0.0000000000e+00
	N22  5.0000000000e+00 5.0000000000e-01 0.0000000000e+00
	N24  5.0000000000e-01 1.0000000000e+00 0.0000000000e+00
	N23  0.0000000000e+00 1.0000000000e+00 0.0000000000e+00
	N25  1.0000000000e+00 1.0000000000e+00 0.0000000000e+00
	N26  1.5000000000e+00 1.0000000000e+00 0.0000000000e+00
	N27  2.0000000000e+00 1.0000000000e+00 0.0000000000e+00
	N28  2.5000000000e+00 1.0000000000e+00 0.0000000000e+00
	N29  3.0000000000e+00 1.0000000000e+00 0.0000000000e+00
	N30  3.5000000000e+00 1.0000000000e+00 0.0000000000e+00
	N31  4.0000000000e+00 1.0000000000e+00 0.0000000000e+00
	N32  4.5000000000e+00 1.0000000000e+00 0.0000000000e+00
	N33  5.0000000000e+00 1.0000000000e+00 0.0000000000e+00
	N35  5.0000000000e-01 1.5000000000e+00 0.0000000000e+00
	N34  0.0000000000e+00 1.5000000000e+00 0.0000000000e+00
	N36  1.0000000000e+00 1.5000000000e+00 0.0000000000e+00
	N37  1.5000000000e+00 1.5000000000e+00 0.0000000000e+00
	N38  2.0000000000e+00 1.5000000000e+00 0.0000000000e+00
	N39  2.5000000000e+00 1.5000000000e+00 0.0000000000e+00
	N40  3.0000000000e+00 1.5000000000e+00 0.0000000000e+00
	N41  3.5000000000e+00 1.5000000000e+00 0.0000000000e+00
	N42  4.0000000000e+00 1.5000000000e+00 0.0000000000e+00
	N43  4.5000000000e+00 1.5000000000e+00 0.0000000000e+00
	N44  5.0000000000e+00 1.5000000000e+00 0.0000000000e+00
	N46  5.0000000000e-01 2.0000000000e+00 0.0000000000e+00
	N45  0.0000000000e+00 2.0000000000e+00 0.0000000000e+00
	N47  1.0000000000e+00 2.0000000000e+00 0.0000000000e+00
	N48  1.5000000000e+00 2.0000000000e+00 0.0000000000e+00
	N49  2.0000000000e+00 2.0000000000e+00 0.0000000000e+00
	N50  2.5000000000e+00 2.0000000000e+00 0.0000000000e+00
	N51  3.0000000000e+00 2.0000000000e+00 0.0000000000e+00
	N52  3.5000000000e+00 2.0000000000e+00 0.0000000000e+00
	N53  4.0000000000e+00 2.0000000000e+00 0.0000000000e+00
	N54  4.5000000000e+00 2.0000000000e+00 0.0000000000e+00
	N55  5.0000000000e+00 2.0000000000e+00 0.0000000000e+00
FINSF


QUAD4 

 E1  N1 N2 N13 N12
 E2  N2 N3 N14 N13
 E3  N3 N4 N15 N14
 E4  N4 N5 N16 N15
 E5  N5 N6 N17 N16
 E6  N6 N7 N18 N17
 E7  N7 N8 N19 N18
 E8  N8 N9 N20 N19
 E9  N9 N10 N21 N20
 E10  N10 N11 N22 N21
 E11  N12 N13 N24 N23
 E12  N13 N14 N25 N24
 E13  N14 N15 N26 N25
 E14  N15 N16 N27 N26
 E15  N16 N17 N28 N27
 E16  N17 N18 N29 N28
 E17  N18 N19 N30 N29
 E18  N19 N20 N31 N30
 E19  N20 N21 N32 N31
 E20  N21 N22 N33 N32
 E21  N23 N24 N35 N34
 E22  N24 N25 N36 N35
 E23  N25 N26 N37 N36
 E24  N26 N27 N38 N37
 E25  N27 N28 N39 N38
 E26  N28 N29 N40 N39
 E27  N29 N30 N41 N40
 E28  N30 N31 N42 N41
 E29  N31 N32 N43 N42
 E30  N32 N33 N44 N43
 E31  N34 N35 N46 N45
 E32  N35 N36 N47 N46
 E33  N36 N37 N48 N47
 E34  N37 N38 N49 N48
 E35  N38 N39 N50 N49
 E36  N39 N40 N51 N50
 E37  N40 N41 N52 N51
 E38  N41 N42 N53 N52
 E39  N42 N43 N54 N53
 E40  N43 N44 N55 N54 
FINSF




GROUP_MA NOM = G1
    % Original part : 1
  E1 E10 E11 E12 E13 E14 E15 
  E16 E17 E18 E19 E2 E20 E21 
  E22 E23 E24 E25 E26 E27 E28 
  E29 E3 E30 E31 E32 E33 E34 
  E35 E36 E37 E38 E39 E4 E40 
  E5 E6 E7 E8 E9
FINSF

FIN

Then, the command file (file “prob24.comm”) :


# prob24
# FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS

DEBUT(PAR_LOT='NON',
      ERREUR=_F(ERREUR_F='ABORT',),
      IGNORE_ALARM=('SUPERVIS_1',),);

MAIL=LIRE_MAILLAGE(FORMAT='ASTER',
                   VERI_MAIL=_F(VERIF='OUI',),);

MODMECA=AFFE_MODELE(MAILLAGE=MAIL,
                    AFFE=_F(GROUP_MA='G1',
                            PHENOMENE='MECANIQUE',
                            MODELISATION='DKT',),);

FCT00001=DEFI_FONCTION(
                       NOM_PARA='FREQ',
                       NOM_RESU='AMOR',
                       VALE=(0.0 ,0.03 ,
                             10.0 ,0.03),
                       INTERPOL=('LIN','LIN',),
                       PROL_DROITE='LINEAIRE',
                       PROL_GAUCHE='EXCLU',);

FCT00002=DEFI_FONCTION(
                       NOM_PARA='FREQ',
                       VALE=(10.0 ,1.0 ,
                             1000.0 ,1.0 ,),
                       INTERPOL=('LIN','LIN',),
                       PROL_DROITE='LINEAIRE',
                       PROL_GAUCHE='LINEAIRE',);

MAT00001=DEFI_MATERIAU(ELAS=_F(E=30000000.0,
                               NU=0.3,
                               RHO=0.00073038,),);

CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
                    AFFE=_F(GROUP_MA='G1',
                            MATER=MAT00001,),);

B1P0=AFFE_CHAR_MECA(MODELE=MODMECA,
                    DDL_IMPO=_F(NOEUD=('N1','N12','N23','N34','N45',),
                                DX=0.0,
                                DY=0.0,
                                DZ=0.0,
                                DRX=0.0,
                                DRY=0.0,),);

CHARGE1=AFFE_CHAR_MECA(MODELE=MODMECA,
                       FORCE_NODALE=_F(NOEUD='N11',
                                       FZ=1.0,),);

CHARGE2=AFFE_CHAR_MECA(MODELE=MODMECA,
                       FORCE_COQUE=_F(GROUP_MA='G1',
                                      PRES=1.0,),);

CHARGE3=AFFE_CHAR_MECA(MODELE=MODMECA,
                       FORCE_COQUE=_F(GROUP_MA='G1',
                                      PRES=1.0,),);

CAEL=AFFE_CARA_ELEM(MODELE=MODMECA,
                    COQUE=_F(GROUP_MA='G1',
                             EPAIS=0.1,
                             VECTEUR=(0.9,0.1,0.2,),
                             COQUE_NCOU=1,),);

ASSEMBLAGE(MODELE=MODMECA,
           CHAM_MATER=CHMAT,
           CARA_ELEM=CAEL,
           CHARGE=B1P0,
           NUME_DDL=CO('NUMEDDL1'),
           MATR_ASSE=(_F(MATRICE=CO('RIGI1'),
                         OPTION='RIGI_MECA',),
                      _F(MATRICE=CO('MASS1'),
                         OPTION='MASS_MECA',),),);

MODES1=CALC_MODES(OPTION='BANDE',
                  MATR_RIGI=RIGI1,
                  MATR_MASS=MASS1,
                  SOLVEUR_MODAL=_F(METHODE='TRI_DIAG'),
                  NORM_MODE=_F(NORME='MASS_GENE'),
                  CALC_FREQ=_F(FREQ=(10.0,2000.0,),),
                  VERI_MODE=_F(STOP_ERREUR='NON',),);

LIMODE1=RECU_TABLE(CO=MODES1,
                   NOM_PARA='FREQ',);

pfreq1 = LIMODE1.EXTR_TABLE().values()['FREQ']

AMOR_I1=CALC_FONC_INTERP(FONCTION=FCT00001,
                         VALE_PARA=pfreq1,);

AMOR_T1=CREA_TABLE(FONCTION=_F(FONCTION=AMOR_I1,),);

AMOR1 = AMOR_T1.EXTR_TABLE().values()['TOUTRESU']


FXEL1_0=CALC_VECT_ELEM(OPTION='CHAR_MECA',
                       CHARGE=(CHARGE2,B1P0,),
                       CARA_ELEM=CAEL,);

FX1_0=ASSE_VECTEUR(VECT_ELEM=FXEL1_0,
                   NUME_DDL=NUMEDDL1,);

FXEL1_1=CALC_VECT_ELEM(OPTION='CHAR_MECA',
                       CHARGE=(CHARGE1,B1P0,),
                       CARA_ELEM=CAEL,);

FX1_1=ASSE_VECTEUR(VECT_ELEM=FXEL1_1,
                   NUME_DDL=NUMEDDL1,);
PROJ_BASE(BASE=MODES1,
          MATR_ASSE_GENE=(_F(MATRICE=CO('MASSGEN1'),
                             MATR_ASSE=MASS1,),
                          _F(MATRICE=CO('RIGIGEN1'),
                             MATR_ASSE=RIGI1,),),
          VECT_ASSE_GENE=(_F(VECTEUR=CO('VG1_0'),
                             VECT_ASSE=FX1_0,),
                          _F(VECTEUR=CO('VG1_1'),
                             VECT_ASSE=FX1_1,),
                          ),);

LSF00001 = DEFI_LIST_FREQ(DEBUT=20.0,
                         INTERVALLE=_F(JUSQU_A=1000.0,
                                       NOMBRE=49),
                         RAFFINEMENT=_F(LIST_RAFFINE=pfreq1,
                                        CRITERE='RELATIF',
                                        DISPERSION=0.03),);

GENE1=DYNA_VIBRA(BASE_CALCUL='GENE',
                 TYPE_CALCUL='HARM',
                 MATR_MASS=MASSGEN1,
                 MATR_RIGI=RIGIGEN1,
                 AMOR_MODAL=_F(AMOR_REDUIT=AMOR1,),
                 LIST_FREQ=LSF00001,
                 EXCIT=(
                 _F(VECT_ASSE_GENE=VG1_0,
                    FONC_MULT=FCT00002,),
                 _F(VECT_ASSE_GENE=VG1_1,
                    FONC_MULT=FCT00002,
                    PHAS_DEG=-45.0,),
                 ),);

RESU1=REST_GENE_PHYS(RESU_GENE=GENE1,
                     TOUT_ORDRE='OUI',
                     NOM_CHAM='DEPL');

FODZ11=RECU_FONCTION(RESULTAT = RESU1,
                      NOM_CMP = 'DZ',
                      NOEUD = 'N11',
                      NOM_CHAM = 'DEPL');

FODZ11M=CALC_FONCTION(EXTRACTION=_F(FONCTION=FODZ11,PARTIE='MODULE'),);                      

TB11M = INFO_FONCTION(MAX=_F(FONCTION=FODZ11M),);  

FODZ33=RECU_FONCTION(RESULTAT = RESU1,
                      NOM_CMP = 'DZ',
                      NOEUD = 'N33',
                      NOM_CHAM = 'DEPL');

FODZ33M=CALC_FONCTION(EXTRACTION=_F(FONCTION=FODZ33,PARTIE='MODULE'),);                                     

TB33M = INFO_FONCTION(MAX=_F(FONCTION=FODZ33M),);

FODZ55=RECU_FONCTION(RESULTAT = RESU1,
                      NOM_CMP = 'DZ',
                      NOEUD = 'N55',
                      NOM_CHAM = 'DEPL');                    
                      

FODZ55M=CALC_FONCTION(EXTRACTION=_F(FONCTION=FODZ55,PARTIE='MODULE'),);

TB55M = INFO_FONCTION(MAX=_F(FONCTION=FODZ55M),);

TBCOMB=CALC_TABLE(TABLE=TB11M,
                  ACTION =(
                    _F(OPERATION='COMB',TABLE = TB33M),
                    _F(OPERATION='COMB',TABLE = TB55M)
                    )
                  )

IMPR_TABLE(TABLE=TBCOMB);

IMPR_FONCTION(FORMAT='XMGRACE',
              UNITE=25,
              ECHELLE_Y='LOG',
              PILOTE = 'PNG',
              COURBE=(_F(FONCTION=FODZ11M,COULEUR=2,MARQUEUR=0,)))              

IMPR_FONCTION(FORMAT='XMGRACE',
              UNITE=26,
              ECHELLE_Y='LOG',
              PILOTE = 'PNG',
              COURBE=(_F(FONCTION=FODZ33M,COULEUR=2,MARQUEUR=0,)))

IMPR_FONCTION(FORMAT='XMGRACE',
              UNITE=27,
              ECHELLE_Y='LOG',
              PILOTE = 'PNG',
              COURBE=(_F(FONCTION=FODZ55M,COULEUR=3,MARQUEUR=0,)))                          

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=(_F(RESULTAT=RESU1,PARTIE='REEL',NOM_CHAM="DEPL",NOM_CHAM_MED="DEPLM"),
                _F(RESULTAT=RESU1,PARTIE='IMAG',NOM_CHAM="DEPL",NOM_CHAM_MED="DEPLP")));

FIN(FORMAT_HDF='OUI',);

You will also need a file to keep everything together ( file “prob24.export”). This file is the one needed to start the calculation on command line (“as_run prob24.export”) :


P actions make_etude
P memory_limit 256.0
P mode interactif
P mpi_nbcpu 1
P mpi_nbnoeud 1
P ncpus 1
P testlist verification sequential
P time_limit 60.0
P version testing
A memjeveux 375.0
A tpmax 172800
P memjob 3072000
P tpsjob 2880
F comm prob24.comm D 1
F mail prob24.mail D 20
F mess prob24.mess R 6
F resu prob24.resu R 8
F rmed prob24.rmed R 80
F libr f11zm.png R 25
F libr f33zm.png R 26
F libr f55zm.png R 27

 

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s