1 # MODIF DATE 09/05/2006 AUTEUR REZETTE C.REZETTE
2 # TITRE FISSURE AU CENTRE D'UNE PLAQUE MINCE RECTANGULAIRE FAISANT
3 # CONFIGURATION MANAGEMENT OF EDF VERSION
4 # ======================================================================
5 # COPYRIGHT (C) 1991 - 2004 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.
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.
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 # OBSTACLE A UN FLUC DE CHALEUR UNIFORME EN MILIEU ISOTROPE
22 # TESTE LES DERIVEES DE G DANS LE CAS D'UN CHARGEMENT EN FORCE DE VOLUME
24 DEBUT( CODE=_F( NOM = 'HPLP101B',NIV_PUB_WEB='INTERNET'))
26 YOUNG = DEFI_PARA_SENSI(VALE=1.)
27 FORC = DEFI_PARA_SENSI(VALE=1.)
31 MA=DEFI_GROUP( reuse=MA, MAILLAGE=MA,
33 GROUP_MA = ( 'LP9P10', 'LP9P10B',
38 MO=AFFE_MODELE( MAILLAGE=MA,
39 AFFE=_F( TOUT = 'OUI',
40 PHENOMENE = 'MECANIQUE',
41 MODELISATION = 'C_PLAN')
43 ZERO = DEFI_CONSTANTE(VALE=0.)
44 NU = DEFI_CONSTANTE(VALE=0.3)
45 ALPHA = DEFI_CONSTANTE(VALE=0.)
47 MAT=DEFI_MATERIAU(ELAS_FO=_F( E = YOUNG, NU = NU,
49 ALPHA = ALPHA, RHO=7800.,))
51 CHMAT=AFFE_MATERIAU( MAILLAGE=MA,
52 AFFE=_F( TOUT = 'OUI',
56 CH=AFFE_CHAR_MECA_F( MODELE=MO,
57 FORCE_INTERNE=_F(TOUT='OUI',FY=FORC,),
59 _F( GROUP_NO = 'LP8P9',
61 _F( GROUP_NO = 'LP8P9B',
67 DEP=MECA_STATIQUE( MODELE=MO, CHAM_MATER=CHMAT,
69 EXCIT=_F( CHARGE = CH),
72 DEP=CALC_ELEM(reuse=DEP,RESULTAT=DEP,
74 OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)
76 FOND=DEFI_FOND_FISS( MAILLAGE=MA,
77 FOND_FISS=_F( GROUP_NO = ('P0',)),
81 THETA0=CALC_THETA( MODELE=MO,
82 THETA_2D=_F( GROUP_NO = ('P0',),
86 DIRECTION=(1., 0., 0.,)
91 THETA01=CALC_THETA( MODELE=MO,
92 THETA_2D=_F( GROUP_NO = ('P0',),
96 DIRECTION=(1., 0., 0.,)
100 # DERIVEE PAR RAPPORT A E
103 DEP2=MECA_STATIQUE( MODELE=MO, CHAM_MATER=CHMAT,
104 SENSIBILITE = (YOUNG),
105 EXCIT=_F( CHARGE = CH),
108 DEP2=CALC_ELEM(reuse=DEP2,RESULTAT=DEP2,
109 SENSIBILITE = (YOUNG),
110 OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)