# MODIF  DATE 09/05/2006   AUTEUR REZETTE C.REZETTE 
# TITRE FISSURE AU CENTRE D'UNE PLAQUE MINCE RECTANGULAIRE FAISANT
#            CONFIGURATION MANAGEMENT OF EDF VERSION
# ======================================================================
# COPYRIGHT (C) 1991 - 2004  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.        
# ======================================================================
#       OBSTACLE A UN FLUC DE CHALEUR UNIFORME EN MILIEU ISOTROPE
#
# TESTE LES DERIVEES DE G DANS LE CAS D'UN CHARGEMENT EN FORCE DE VOLUME

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

YOUNG = DEFI_PARA_SENSI(VALE=1.)
FORC = DEFI_PARA_SENSI(VALE=1.)

MA=LIRE_MAILLAGE( )

MA=DEFI_GROUP( reuse=MA,   MAILLAGE=MA,
                           CREA_GROUP_NO=_F( 
                         GROUP_MA = ( 'LP9P10',  'LP9P10B', 
                                      'LP8P9',   'LP8P9B',
                                         ))
                       )

MO=AFFE_MODELE( MAILLAGE=MA, 
                             AFFE=_F(  TOUT = 'OUI',
                                    PHENOMENE = 'MECANIQUE',
                                    MODELISATION = 'C_PLAN')
                       )
ZERO =    DEFI_CONSTANTE(VALE=0.)
NU =    DEFI_CONSTANTE(VALE=0.3)
ALPHA =    DEFI_CONSTANTE(VALE=0.)

MAT=DEFI_MATERIAU(ELAS_FO=_F(  E = YOUNG, NU = NU,
                                TEMP_DEF_ALPHA =20.,
                               ALPHA = ALPHA, RHO=7800.,))

CHMAT=AFFE_MATERIAU(    MAILLAGE=MA,
                                  AFFE=_F(  TOUT = 'OUI',
                                         MATER = MAT)
                          )

CH=AFFE_CHAR_MECA_F(           MODELE=MO,
                       FORCE_INTERNE=_F(TOUT='OUI',FY=FORC,),
                               DDL_IMPO=(
                                    _F(  GROUP_NO = 'LP8P9',
                                               DX = ZERO),
                                             _F(  GROUP_NO = 'LP8P9B',
                                               DX = ZERO),
                                             _F(  GROUP_NO = 'P11',
                                               DY = ZERO))
                            )

DEP=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT, 
                      SENSIBILITE = (FORC), 
                              EXCIT=_F( CHARGE = CH),
                     )

DEP=CALC_ELEM(reuse=DEP,RESULTAT=DEP,   
                          SENSIBILITE = (FORC),
              OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)

FOND=DEFI_FOND_FISS(    MAILLAGE=MA,
                              FOND_FISS=_F( GROUP_NO = ('P0',)),
                           NORMALE=(0., 1., 0.,)
                      )

THETA0=CALC_THETA(       MODELE=MO,
                         THETA_2D=_F(  GROUP_NO = ('P0',),
                                    MODULE = 1.,
                                    R_INF = 3.75E-5,
                                    R_SUP = 7.50E-5),
                        DIRECTION=(1., 0., 0.,)
                     )

G0=CALC_G(      RESULTAT=DEP,
                        THETA=_F(THETA=THETA0),
                        SYME_CHAR='SANS',
                        COMP_ELAS=_F( RELATION = 'ELAS',
                                      DEFORMATION = 'PETIT')
                      )

DG0=CALC_G(     RESULTAT=DEP,
                        SENSIBILITE = (FORC), 
                        THETA=_F(THETA=THETA0),
                        SYME_CHAR='SANS',
                        COMP_ELAS=_F( RELATION = 'ELAS',
                                      DEFORMATION = 'PETIT')
                      )

TEST_TABLE(       TABLE=DG0,
                       NOM_PARA='DG/DF',
                          SENSIBILITE = (FORC), 
                           VALE=1.06E-2,
                      REFERENCE='ANALYTIQUE',
                      PRECISION=0.0005    )
FIN()

THETA01=CALC_THETA(       MODELE=MO,
                          THETA_2D=_F(  GROUP_NO = ('P0',),
                                     MODULE = 1.,
                                     R_INF = 7.50E-5,
                                     R_SUP = 1.125E-4),
                         DIRECTION=(1., 0., 0.,)
                     )

G1=CALC_G(        RESULTAT=DEP,
                          THETA=_F(THETA=THETA01),
                          COMP_ELAS=_F( RELATION = 'ELAS',
                                     DEFORMATION = 'PETIT'),
                          SYME_CHAR='SANS'
                     )

DG1=CALC_G(      RESULTAT=DEP,
                         SENSIBILITE = (FORC), 
                         THETA=_F(THETA=THETA01),
                         SYME_CHAR='SANS',
                          COMP_ELAS=_F( RELATION = 'ELAS',
                                      DEFORMATION = 'PETIT')
                      )

TEST_TABLE(       TABLE=DG1,
                       NOM_PARA='G',
                          SENSIBILITE = (FORC), 
                      REFERENCE='ANALYTIQUE',
                           VALE=1.06E-2,
                      PRECISION=0.0005    )

#
#  DERIVEE PAR RAPPORT A E
#

DEP2=MECA_STATIQUE(   MODELE=MO,   CHAM_MATER=CHMAT, 
                      SENSIBILITE = (YOUNG), 
                              EXCIT=_F( CHARGE = CH),
                   )

DEP2=CALC_ELEM(reuse=DEP2,RESULTAT=DEP2,   
                          SENSIBILITE = (YOUNG), 
              OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)


DG2=CALC_G(     RESULTAT=DEP2,
                        SENSIBILITE = (YOUNG), 
                        THETA=_F(THETA=THETA0),
                        SYME_CHAR='SANS',
                         COMP_ELAS=_F( RELATION = 'ELAS',
                                      DEFORMATION = 'PETIT')
                      )

TEST_TABLE(       TABLE=DG2,
                       NOM_PARA='G',
                          SENSIBILITE = (YOUNG), 
                      REFERENCE='ANALYTIQUE',
                           VALE=-5.3E-3,
                      PRECISION=0.0005    )


DG3=CALC_G(     RESULTAT=DEP2,
                        SENSIBILITE = (YOUNG), 
                        THETA=_F(THETA=THETA01),
                        SYME_CHAR='SANS',
                        COMP_ELAS=_F( RELATION = 'ELAS',
                                      DEFORMATION = 'PETIT')
                      )

TEST_TABLE(       TABLE=DG3,
                       NOM_PARA='G',
                           VALE=-5.3E-3,
                          SENSIBILITE = (YOUNG), 
                      REFERENCE='ANALYTIQUE',
                      PRECISION=0.0005    )


FIN()