# MODIF  DATE 04/10/2005   AUTEUR CIBHHPD L.SALMONA 
# TITRE CABLE PESANT AVEC DILATATION THERMIQUE. MODELISATION A : CABLE
#            CONFIGURATION MANAGEMENT OF EDF VERSION
# ======================================================================
# COPYRIGHT (C) 1991 - 2001  EDF R&D                  WWW.CODE-ASTER.ORG
# THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
# THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR   
# (AT YOUR OPTION) ANY LATER VERSION.                                 
#
# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT 
# WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF          
# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU    
# GENERAL PUBLIC LICENSE FOR MORE DETAILS.                            
#
# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE   
# ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,       
#    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.      
# ======================================================================

DEBUT( CODE=_F(  NOM = 'SSNL114A',NIV_PUB_WEB='INTERNET') )

MA2=LIRE_MAILLAGE()

M1=DEFI_MATERIAU( 
         ELAS=_F(  E = 5.700000E+10, ALPHA = 2.300000E-05,
        RHO = 2.844230E+03, NU=0.),
        CABLE=_F(),)

CHM2=AFFE_MATERIAU(    MAILLAGE=MA2,AFFE=(
             _F(  GROUP_MA = (   'O0101010', ),  MATER = M1),
             _F(  GROUP_MA = (   'P0101010', 'P0101011', ),
                   MATER = M1)))

MO2=AFFE_MODELE(   MAILLAGE=MA2,
             AFFE=_F(  GROUP_MA = ( 'P0101010', 'O0101010', 'P0101011', ),
           PHENOMENE = 'MECANIQUE',
         MODELISATION = 'CABLE'))

CARA2=AFFE_CARA_ELEM(    MODELE=MO2, 
         CABLE=_F(  GROUP_MA = ( 'P0101010', 'P0101011','O0101010',  ),
        SECTION = 2.278300E-04)
                        )

CHDDL2=AFFE_CHAR_MECA(    MODELE=MO2,
            DDL_IMPO=_F(  GROUP_NO = ( 'ZFIXES', ),
          DX = 0.0, DY = 0.0, DZ = 0.0))

CHTN2=CREA_CHAMP( OPERATION='AFFE', TYPE_CHAM='NOEU_TEMP_R', MAILLAGE=MA2,
           AFFE=_F(  TOUT = 'OUI', NOM_CMP = 'TEMP', VALE = 39.26))

CHTN0=CREA_CHAMP( OPERATION='AFFE', TYPE_CHAM='NOEU_TEMP_R', MAILLAGE=MA2,
           AFFE=_F(  TOUT = 'OUI', NOM_CMP = 'TEMP', VALE = 0.0))

RESUTHER=CREA_RESU(OPERATION='AFFE', TYPE_RESU='EVOL_THER',    NOM_CHAM='TEMP',AFFE=(
 _F( INST = -1., CHAM_GD = CHTN0),
# MONTEE EN TEMPERATURE PROGRESSIVE
 _F( INST = 0., CHAM_GD = CHTN0),
 _F( INST = 1., CHAM_GD = CHTN2))
                   )

CHTP2=AFFE_CHAR_MECA(    MODELE=MO2,   TEMP_CALCULEE=RESUTHER )

CHPES2=AFFE_CHAR_MECA(    MODELE=MO2,
                   PESANTEUR=(9.81,0.0,0.0,-1.0,) )

L_INST0=DEFI_LIST_REEL(    DEBUT=-1.000000E+00,
         INTERVALLE=_F( JUSQU_A = 1.000000E+00, NOMBRE = 2))

RESU1=STAT_NON_LINE( 
                            MODELE=MO2,
                        CHAM_MATER=CHM2,
                         CARA_ELEM=CARA2,EXCIT=(
                             _F( 
                CHARGE = CHDDL2), _F( 
                CHARGE = CHPES2,
                TYPE_CHARGE = 'SUIV'), _F( 
                CHARGE = CHTP2)),
                         INCREMENT=_F( 
                LIST_INST = L_INST0,
                NUME_INST_FIN = 2),
                         COMP_ELAS=_F( 
                GROUP_MA = (
                  'P0101010',  'P0101011',
                  'O0101010',  ),
                RELATION = 'CABLE',
                DEFORMATION = 'GREEN'),
                       CONVERGENCE=_F( 
                RESI_GLOB_RELA = 1.00E-07,
                ITER_GLOB_MAXI = 50),
                            NEWTON=_F( 
                REAC_ITER = 1)
              )

RESU1=CALC_ELEM( reuse=RESU1,
                            OPTION=(
                  'SIEF_ELNO_ELGA',  'VARI_ELNO_ELGA',),
                          RESULTAT=RESU1,
              )

# SOLUTION ANALYTIQUE : FORTRAN DANS SSNL114.38


TEST_RESU(RESU=(_F( RESULTAT = RESU1,
                 NUME_ORDRE = 1,
                 NOM_CHAM = 'DEPL',
                 NOEUD = 'N36',
                 NOM_CMP = 'DZ',
                 VALE = -6.352,
                 PRECISION = 1.E-3),
           _F( RESULTAT = RESU1,
                 NUME_ORDRE = 2,
                 NOM_CHAM = 'DEPL',
                 NOEUD = 'N36',
                 NOM_CMP = 'DZ',
                 VALE = -8.1955,
                 PRECISION = 1.E-3))
          )

FIN( )
#