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.,)
89 G0=CALC_G( RESULTAT=DEP,
90 THETA=_F(THETA=THETA0),
92 COMP_ELAS=_F( RELATION = 'ELAS',
93 DEFORMATION = 'PETIT')
96 DG0=CALC_G( RESULTAT=DEP,
98 THETA=_F(THETA=THETA0),
100 COMP_ELAS=_F( RELATION = 'ELAS',
101 DEFORMATION = 'PETIT')
104 TEST_TABLE( TABLE=DG0,
106 SENSIBILITE = (FORC),
108 REFERENCE='ANALYTIQUE',
112 THETA01=CALC_THETA( MODELE=MO,
113 THETA_2D=_F( GROUP_NO = ('P0',),
117 DIRECTION=(1., 0., 0.,)
120 G1=CALC_G( RESULTAT=DEP,
121 THETA=_F(THETA=THETA01),
122 COMP_ELAS=_F( RELATION = 'ELAS',
123 DEFORMATION = 'PETIT'),
127 DG1=CALC_G( RESULTAT=DEP,
128 SENSIBILITE = (FORC),
129 THETA=_F(THETA=THETA01),
131 COMP_ELAS=_F( RELATION = 'ELAS',
132 DEFORMATION = 'PETIT')
135 TEST_TABLE( TABLE=DG1,
137 SENSIBILITE = (FORC),
138 REFERENCE='ANALYTIQUE',
143 # DERIVEE PAR RAPPORT A E
146 DEP2=MECA_STATIQUE( MODELE=MO, CHAM_MATER=CHMAT,
147 SENSIBILITE = (YOUNG),
148 EXCIT=_F( CHARGE = CH),
151 DEP2=CALC_ELEM(reuse=DEP2,RESULTAT=DEP2,
152 SENSIBILITE = (YOUNG),
153 OPTION=('EPSI_ELGA_DEPL','SIEF_ELGA_DEPL',),)
156 DG2=CALC_G( RESULTAT=DEP2,
157 SENSIBILITE = (YOUNG),
158 THETA=_F(THETA=THETA0),
160 COMP_ELAS=_F( RELATION = 'ELAS',
161 DEFORMATION = 'PETIT')
164 TEST_TABLE( TABLE=DG2,
166 SENSIBILITE = (YOUNG),
167 REFERENCE='ANALYTIQUE',
172 DG3=CALC_G( RESULTAT=DEP2,
173 SENSIBILITE = (YOUNG),
174 THETA=_F(THETA=THETA01),
176 COMP_ELAS=_F( RELATION = 'ELAS',
177 DEFORMATION = 'PETIT')
180 TEST_TABLE( TABLE=DG3,
183 SENSIBILITE = (YOUNG),
184 REFERENCE='ANALYTIQUE',