src/TestsEltsVisu/Nicolas/ssnl116a.comm

   1 # MODIF  DATE 28/05/2004   AUTEUR LEBOUVIE F.LEBOUVIER
   2 # TITRE CALCUL DU TRONCON BE  DU CIG DES RENARDIERES
   3 #            CONFIGURATION MANAGEMENT OF EDF VERSION
   4 # ======================================================================
   5 # COPYRIGHT (C) 1991 - 2001  EDF R&D                  WWW.CODE-ASTER.ORG
   6 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
   7 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
   8 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
   9 # (AT YOUR OPTION) ANY LATER VERSION.
  10 #
  11 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  12 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  13 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  14 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  15 #
  16 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
  17 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
  18 #    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
  19 # ======================================================================
  20
  21 # REF : ALTHEE STROOBANT,STAGIAIRE MTC AOUT 1999
  22
  23 DEBUT(  CODE=_F(  NOM = 'SSNL116A',NIV_PUB_WEB='INTRANET') )
  24
  25 MA=LIRE_MAILLAGE(  )
  26
  27 MA=DEFI_GROUP( reuse=MA,  MAILLAGE=MA,CREA_GROUP_NO=(
  28                    _F( NOM = 'GONAA',
  29                                   GROUP_MA = 'GOA'),
  30                    _F( NOM = 'GONZZ',
  31                                   GROUP_MA = 'GRBAR_V'),
  32                    _F( NOM = 'GONXX',
  33                                   GROUP_MA = 'GRBAR_T'),
  34                    _F( NOM = 'GON',
  35                                   GROUP_MA = 'GRBAR_H'))
  36                  )
  37
  38 MO=AFFE_MODELE(    MAILLAGE=MA,
  39      VERIF=('MAILLE',  'NOEUD',),AFFE=(
  40      _F(  GROUP_MA = 'GOA',
  41      PHENOMENE = 'MECANIQUE',
  42      MODELISATION = 'POU_D_E'),
  43      _F(  GROUP_MA = ('GRBAR_T','GRBAR_H','GRBAR_V',),
  44      PHENOMENE = 'MECANIQUE',
  45      MODELISATION = 'BARRE'))
  46      )
  47
  48 MATA=DEFI_MATERIAU(  ELAS=_F(  E = 7.2E10,
  49     NU = 0.3,
  50     ALPHA = 22.4E-6),
  51          ECRO_FLEJOU=_F(
  52                EP = 3.7E10,
  53           SY = 75.E6,
  54           SU = 190.E6,
  55           PUISS = 0.29),
  56   VMIS_POUTRE=_F(  NP = 1.2699E6,
  57     MEY = 1.248E5,
  58     MPY = 1.589E5,
  59     CAY = 0.84,
  60     CBY = 0.0012,
  61     MEZ = 1.248E5,
  62     MPZ = 1.589E5,
  63     CAZ = 0.84,
  64     CBZ = 0.0012,
  65      MPX = 1.E10)
  66       )
  67
  68 INCLUDE(UNITE=91)
  69
  70 CARAP=AFFE_CARA_ELEM(      MODELE=MO,
  71       POUTRE=_F( GROUP_MA = 'GOA',
  72      SECTION = 'CERCLE',
  73      CARA = ('R',  'EP',),
  74      VALE = (0.25765,  0.01,)),BARRE=(
  75       _F( GROUP_MA = 'GRBAR_T',
  76      SECTION = 'RECTANGLE',
  77      CARA = 'H',
  78      VALE = 1.),
  79              _F( GROUP_MA = 'GRBAR_H',
  80      SECTION = 'RECTANGLE',
  81      CARA = 'H',
  82      VALE = 1.),
  83                _F( GROUP_MA = 'GRBAR_V',
  84      SECTION = 'RECTANGLE',
  85      CARA = 'H',
  86      VALE = 1.))
  87
  88           )
  89
  90 LITPS=DEFI_LIST_REEL( DEBUT=0.0,INTERVALLE=(
  91     _F( JUSQU_A = 1., NOMBRE = 10),
  92       _F( JUSQU_A = 2., NOMBRE = 2),
  93       _F( JUSQU_A = 3., NOMBRE = 30),
  94       _F( JUSQU_A = 4., NOMBRE = 2))
  95       )
  96
  97 # FONCTION TEMPERATURE POUR LES BARRES : DT = 50
  98 #_______________________________________________________#
  99
 100 DXA=DEFI_FONCTION(   NOM_PARA='INST',
 101      VALE=(  0.,    0.,
 102      1.,   -0.004,
 103      2.,   -0.004,
 104      3.,   0.002,
 105      4.,   0.002,
 106      ),
 107    PROL_DROITE='LINEAIRE',
 108    PROL_GAUCHE='LINEAIRE'
 109    )
 110
 111 # FONCTION TEMPERATURE POUR LES POUTRES : DT = 50
 112 #_______________________________________________________#
 113
 114 NUL=DEFI_FONCTION(   NOM_PARA='INST',
 115      VALE=(  0.,    0.,
 116      1.,   0.,
 117      2.,   0.,
 118      3.,    0.,
 119      4.,    0.,
 120      ),
 121    PROL_DROITE='LINEAIRE',
 122    PROL_GAUCHE='LINEAIRE'
 123    )
 124
 125 CHAMA=AFFE_CHAR_MECA_F(    MODELE=MO,DDL_IMPO=(
 126
 127    _F( GROUP_NO = 'PC01', DX = NUL,  DY = NUL,  DZ = NUL,
 128       DRX = NUL, DRY = NUL, DRZ = NUL),
 129    _F( GROUP_NO = 'PC02', DX = NUL,  DY = NUL,  DZ = DXA,
 130       DRX = NUL, DRY = NUL, DRZ = NUL),
 131
 132    _F( GROUP_NO = 'GNO_EXBA', DX = NUL, DY = NUL, DZ = NUL))
 133
 134         )
 135
 136 RES1=STAT_NON_LINE(
 137                             MODELE=MO,
 138                         CHAM_MATER=CHMATF,
 139                          CARA_ELEM=CARAP,
 140                              EXCIT=_F(
 141                 CHARGE = CHAMA),COMP_INCR=(
 142                          _F(
 143                 RELATION = 'VMIS_POU_FLEJOU',
 144                 GROUP_MA = 'GOA'),
 145                          _F(
 146                 RELATION = 'VMIS_ISOT_LINE',
 147                 GROUP_MA = (
 148                   'GRBAR_T',  'GRBAR_H',)),
 149                          _F(
 150                 RELATION = 'VMIS_ASYM_LINE',
 151                 GROUP_MA = (
 152                   'GRBAR_V',  ))),
 153                          INCREMENT=_F(
 154                 LIST_INST = LITPS,
 155                 NUME_INST_FIN = 44),
 156                             NEWTON=_F(
 157                 MATRICE = 'ELASTIQUE'),
 158                        CONVERGENCE=_F(
 159                 RESI_GLOB_MAXI = 5.E-3,
 160                 ITER_GLOB_MAXI = 20)
 161               )
 162
 163 RES1=CALC_ELEM( reuse=RES1,
 164                             OPTION='SIEF_ELNO_ELGA',
 165                           RESULTAT=RES1,
 166               )
 167
 168 TEST_RESU(RESU=(
 169
 170     _F( RESULTAT = RES1, NOEUD = 'NPC02', NOM_CHAM = 'DEPL', NOM_CMP = 'DZ',
 171      NUME_ORDRE = 1,  VALE = -4.00000E-04,  REFERENCE = 'NON_REGRESSION'),
 172
 173     _F( RESULTAT = RES1, NOEUD = 'NPC02', NOM_CHAM = 'DEPL', NOM_CMP = 'DZ',
 174      NUME_ORDRE = 10, VALE = -4.00000E-03, REFERENCE = 'NON_REGRESSION'),
 175
 176     _F( RESULTAT = RES1, NOEUD = 'NPC02', NOM_CHAM = 'DEPL', NOM_CMP = 'DZ',
 177      NUME_ORDRE = 30, VALE = -4.00000E-04,  REFERENCE = 'NON_REGRESSION'),
 178
 179     _F( RESULTAT = RES1, NOEUD = 'NPC02', NOM_CHAM = 'DEPL', NOM_CMP = 'DZ',
 180      NUME_ORDRE = 44, VALE = 2.00000E-03,  REFERENCE = 'NON_REGRESSION'),
 181
 182 _F( RESULTAT = RES1, MAILLE = 'MV01F010', NOEUD = 'NPC02', NOM_CHAM = 'SIEF_ELNO_ELGA',
 183 NUME_ORDRE = 1,  VALE = -2.00000E+03, REFERENCE = 'NON_REGRESSION', NOM_CMP = 'N'),
 184
 185 _F( RESULTAT = RES1, MAILLE = 'MV01F010', NOEUD = 'NPC02', NOM_CHAM = 'SIEF_ELNO_ELGA',
 186 NUME_ORDRE = 10, VALE = -1.20000E+04, REFERENCE = 'NON_REGRESSION', NOM_CMP = 'N'),
 187
 188 _F( RESULTAT = RES1, MAILLE = 'MV01F010', NOEUD = 'NPC02', NOM_CHAM = 'SIEF_ELNO_ELGA',
 189 NUME_ORDRE = 30, VALE = 5.20000E+03, REFERENCE = 'NON_REGRESSION', NOM_CMP = 'N'),
 190
 191 _F( RESULTAT = RES1, MAILLE = 'MV01F010', NOEUD = 'NPC02', NOM_CHAM = 'SIEF_ELNO_ELGA',
 192 NUME_ORDRE = 44, VALE = 7.60000E+03, REFERENCE = 'NON_REGRESSION', NOM_CMP = 'N'))
 193
 194
 195           )
 196
 197 FIN()
 198 #
 199 #
 200 #
 201 #TE1=POST_RELEVE_T(ACTION:( INTITULE:'T1'
 202 #    GROUP_NO:GONZZ
 203 #    GROUP_MA:GRBAR_V
 204 #    RESULTAT:RES1
 205 #    NOM_CHAM: 'DEPL'
 206 #    TOUT_ORDRE:'OUI'
 207 #    TOUT_CMP:'OUI'
 208 #    OPERATION:'EXTRACTION')
 209 #    );
 210 #
 211 ##MPR_TABLE (  TABLE : TE1
 212 #  NOM_PARA:('COOR_X','DZ')
 213 #  FILTRE : ( NOM_PARA : 'NOEUD' VALE_K : NPC02 )
 214 #  FORMAT:'EXCEL' );
 215 #
 216 #
 217 #
 218 #TB1=POST_RELEVE_T(ACTION:( INTITULE:'TBE1'
 219 #    GROUP_NO:GONZZ
 220 #    GROUP_MA:GRBAR_V
 221 #    RESULTAT:RES1
 222 #    NOM_CHAM: 'SIEF_ELNO_ELGA'
 223 #    TOUT_ORDRE:'OUI'
 224 #    TOUT_CMP:'OUI'
 225 #    OPERATION:'EXTRACTION')
 226 #    );
 227 #
 228 ##MPR_TABLE (  TABLE : TB1
 229 #  NOM_PARA:('COOR_X','N')
 230 #  FILTRE : ( NOM_PARA : 'NOEUD' VALE_K : NV01F010 )
 231 #  FORMAT:'EXCEL' );
 232 #
 233 #
 234 #
 235 #