src/TestsEltsVisu/Nicolas/ssnl114a.comm

   1 # MODIF  DATE 04/10/2005   AUTEUR CIBHHPD L.SALMONA
   2 # TITRE CABLE PESANT AVEC DILATATION THERMIQUE. MODELISATION A : CABLE
   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 DEBUT( CODE=_F(  NOM = 'SSNL114A',NIV_PUB_WEB='INTERNET') )
  22
  23 MA2=LIRE_MAILLAGE()
  24
  25 M1=DEFI_MATERIAU(
  26          ELAS=_F(  E = 5.700000E+10, ALPHA = 2.300000E-05,
  27         RHO = 2.844230E+03, NU=0.),
  28         CABLE=_F(),)
  29
  30 CHM2=AFFE_MATERIAU(    MAILLAGE=MA2,AFFE=(
  31              _F(  GROUP_MA = (   'O0101010', ),  MATER = M1),
  32              _F(  GROUP_MA = (   'P0101010', 'P0101011', ),
  33                    MATER = M1)))
  34
  35 MO2=AFFE_MODELE(   MAILLAGE=MA2,
  36              AFFE=_F(  GROUP_MA = ( 'P0101010', 'O0101010', 'P0101011', ),
  37            PHENOMENE = 'MECANIQUE',
  38          MODELISATION = 'CABLE'))
  39
  40 CARA2=AFFE_CARA_ELEM(    MODELE=MO2,
  41          CABLE=_F(  GROUP_MA = ( 'P0101010', 'P0101011','O0101010',  ),
  42         SECTION = 2.278300E-04)
  43                         )
  44
  45 CHDDL2=AFFE_CHAR_MECA(    MODELE=MO2,
  46             DDL_IMPO=_F(  GROUP_NO = ( 'ZFIXES', ),
  47           DX = 0.0, DY = 0.0, DZ = 0.0))
  48
  49 CHTN2=CREA_CHAMP( OPERATION='AFFE', TYPE_CHAM='NOEU_TEMP_R', MAILLAGE=MA2,
  50            AFFE=_F(  TOUT = 'OUI', NOM_CMP = 'TEMP', VALE = 39.26))
  51
  52 CHTN0=CREA_CHAMP( OPERATION='AFFE', TYPE_CHAM='NOEU_TEMP_R', MAILLAGE=MA2,
  53            AFFE=_F(  TOUT = 'OUI', NOM_CMP = 'TEMP', VALE = 0.0))
  54
  55 RESUTHER=CREA_RESU(OPERATION='AFFE', TYPE_RESU='EVOL_THER',    NOM_CHAM='TEMP',AFFE=(
  56  _F( INST = -1., CHAM_GD = CHTN0),
  57 # MONTEE EN TEMPERATURE PROGRESSIVE
  58  _F( INST = 0., CHAM_GD = CHTN0),
  59  _F( INST = 1., CHAM_GD = CHTN2))
  60                    )
  61
  62 CHTP2=AFFE_CHAR_MECA(    MODELE=MO2,   TEMP_CALCULEE=RESUTHER )
  63
  64 CHPES2=AFFE_CHAR_MECA(    MODELE=MO2,
  65                    PESANTEUR=(9.81,0.0,0.0,-1.0,) )
  66
  67 L_INST0=DEFI_LIST_REEL(    DEBUT=-1.000000E+00,
  68          INTERVALLE=_F( JUSQU_A = 1.000000E+00, NOMBRE = 2))
  69
  70 RESU1=STAT_NON_LINE(
  71                             MODELE=MO2,
  72                         CHAM_MATER=CHM2,
  73                          CARA_ELEM=CARA2,EXCIT=(
  74                              _F(
  75                 CHARGE = CHDDL2), _F(
  76                 CHARGE = CHPES2,
  77                 TYPE_CHARGE = 'SUIV'), _F(
  78                 CHARGE = CHTP2)),
  79                          INCREMENT=_F(
  80                 LIST_INST = L_INST0,
  81                 NUME_INST_FIN = 2),
  82                          COMP_ELAS=_F(
  83                 GROUP_MA = (
  84                   'P0101010',  'P0101011',
  85                   'O0101010',  ),
  86                 RELATION = 'CABLE',
  87                 DEFORMATION = 'GREEN'),
  88                        CONVERGENCE=_F(
  89                 RESI_GLOB_RELA = 1.00E-07,
  90                 ITER_GLOB_MAXI = 50),
  91                             NEWTON=_F(
  92                 REAC_ITER = 1)
  93               )
  94
  95 RESU1=CALC_ELEM( reuse=RESU1,
  96                             OPTION=(
  97                   'SIEF_ELNO_ELGA',  'VARI_ELNO_ELGA',),
  98                           RESULTAT=RESU1,
  99               )
 100
 101 # SOLUTION ANALYTIQUE : FORTRAN DANS SSNL114.38
 102
 103
 104 TEST_RESU(RESU=(_F( RESULTAT = RESU1,
 105                  NUME_ORDRE = 1,
 106                  NOM_CHAM = 'DEPL',
 107                  NOEUD = 'N36',
 108                  NOM_CMP = 'DZ',
 109                  VALE = -6.352,
 110                  PRECISION = 1.E-3),
 111            _F( RESULTAT = RESU1,
 112                  NUME_ORDRE = 2,
 113                  NOM_CHAM = 'DEPL',
 114                  NOEUD = 'N36',
 115                  NOM_CMP = 'DZ',
 116                  VALE = -8.1955,
 117                  PRECISION = 1.E-3))
 118           )
 119
 120 FIN( )
 121 #