Aster/Cata/cataSTA8/materiau/Z2CN1810_REF_A.NOMI

   1 # MODIF  DATE 24/10/2006   AUTEUR DURAND C.DURAND
   2 #            CONFIGURATION MANAGEMENT OF EDF VERSION
   3 # ======================================================================
   4 # COPYRIGHT (C) 1991 - 2001  EDF R&D                  WWW.CODE-ASTER.ORG
   5 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
   6 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
   7 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
   8 # (AT YOUR OPTION) ANY LATER VERSION.
   9 #
  10 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  11 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  12 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  13 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  14 #
  15 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
  16 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
  17 #    1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
  18 # ======================================================================
  19 #
  20 # LAMBDA EN FONCTION DE LA TEMPéRATURE
  21 #
  22 coef0=1.E ## UNIT -3
  23 _A0=DEFI_FONCTION(      NOM_PARA='TEMP',
  24                             PROL_DROITE='EXCLU',
  25                             PROL_GAUCHE='EXCLU',
  26                                    VALE=(
  27           20.            ,14.7            *coef0,
  28           50.            ,15.2            *coef0,
  29           100.           ,15.8            *coef0,
  30           150.           ,16.7            *coef0,
  31           200.           ,17.2            *coef0,
  32           250.           ,18.             *coef0,
  33           300.           ,18.6            *coef0,
  34           350.           ,19.3            *coef0,
  35           400.           ,20.             *coef0,
  36           450.           ,20.5            *coef0,
  37           500.           ,21.1            *coef0,
  38           550.           ,21.7            *coef0,
  39           600.           ,22.2            *coef0,
  40           650.           ,22.7            *coef0,
  41           700.           ,23.2            *coef0,
  42           750.           ,23.7            *coef0,
  43           800.           ,24.1            *coef0,
  44           950.           ,26.67           *coef0,
  45           1150.          ,29.24           *coef0,
  46           1370.          ,32.06           *coef0,
  47           ))
  48
  49 #
  50 # RHO_CP EN FONCTION DE LA TEMPéRATURE
  51 #
  52 coef1=1.E ## UNIT -9
  53 _A1=DEFI_FONCTION(      NOM_PARA='TEMP',
  54                             PROL_DROITE='EXCLU',
  55                             PROL_GAUCHE='EXCLU',
  56                               VALE=(
  57           20.            ,3.60E+6         *coef1,
  58           50.            ,3.74E+6         *coef1,
  59           100.           ,3.90E+6         *coef1,
  60           150.           ,4.10E+6         *coef1,
  61           200.           ,4.16E+6         *coef1,
  62           250.           ,4.27E+6         *coef1,
  63           300.           ,4.30E+6         *coef1,
  64           350.           ,4.35E+6         *coef1,
  65           400.           ,4.39E+6         *coef1,
  66           450.           ,4.39E+6         *coef1,
  67           500.           ,4.44E+6         *coef1,
  68           550.           ,4.47E+6         *coef1,
  69           600.           ,4.49E+6         *coef1,
  70           650.           ,4.53E+6         *coef1,
  71           700.           ,4.58E+6         *coef1,
  72           750.           ,4.61E+6         *coef1,
  73           800.           ,4.72E+6         *coef1,
  74           1000.          ,4.99E+6         *coef1,
  75           1200.          ,5.04E+6         *coef1,
  76           1500.          ,5.04E+6         *coef1,
  77           ))
  78
  79 #
  80 # E EN FONCTION DE LA TEMPéRATURE
  81 #
  82 coef2=1.E ## UNIT -6
  83 _A2=DEFI_FONCTION(      NOM_PARA='TEMP',
  84                             PROL_DROITE='EXCLU',
  85                             PROL_GAUCHE='EXCLU',
  86                               VALE=(
  87           0.             ,198.5E+9        *coef2,
  88           20.            ,197.E+9         *coef2,
  89           50.            ,195.E+9         *coef2,
  90           100.           ,191.5E+9        *coef2,
  91           150.           ,187.5E+9        *coef2,
  92           200.           ,184.E+9         *coef2,
  93           250.           ,180.E+9         *coef2,
  94           300.           ,176.5E+9        *coef2,
  95           350.           ,172.E+9         *coef2,
  96           400.           ,168.E+9         *coef2,
  97           450.           ,164.E+9         *coef2,
  98           500.           ,160.E+9         *coef2,
  99           550.           ,155.5E+9        *coef2,
 100           600.           ,151.5E+9        *coef2,
 101           700.           ,142.5E+9        *coef2,
 102           800.           ,130.E+9         *coef2,
 103           1000.          ,81.5E+9         *coef2,
 104           1200.          ,7.4E+9          *coef2,
 105           1400.          ,0.              *coef2,
 106           ))
 107
 108 #
 109 # NU EN FONCTION DE LA TEMPéRATURE
 110 #
 111 _A3=DEFI_FONCTION(      NOM_PARA='TEMP',
 112                             PROL_DROITE='CONSTANT',
 113                             PROL_GAUCHE='CONSTANT',
 114                                    VALE=(   0.,    0.3,
 115                                         ))
 116
 117 #
 118 # ALPHA EN FONCTION DE LA TEMPéRATURE
 119 #
 120
 121 _A4=DEFI_FONCTION(      NOM_PARA='TEMP',
 122                             PROL_DROITE='EXCLU',
 123                             PROL_GAUCHE='EXCLU',
 124                               VALE=(   20., 16.4E-6,
 125                                        50., 16.54E-6,    100.,  16.8E-6,
 126                                       150., 17.04E-6,    200.,  17.2E-6,
 127                                       250., 17.5E-6,     300.,  17.7E-6,
 128                                       350., 17.9E-6,     400.,  18.1E-6,
 129                                       450., 18.24E-6,    500.,  18.4E-6,
 130                                       600., 18.7E-6,     700.,  18.9E-6,
 131                                       800., 19.1E-6,     900.,  19.3E-6,
 132                                      1000., 19.5E-6,    1400.,  20.E-6,
 133                                      1600., 20.2E-6,
 134                                    ))
 135
 136 #
 137 # COURBE DE TRACTION à LA TEMPéRATURE 20°C
 138 #
 139
 140 coef5=1.E ## UNIT -6
 141 _A5=DEFI_FONCTION(      NOM_PARA='EPSI',
 142                            PROL_DROITE='EXCLU',
 143                            PROL_GAUCHE='EXCLU',
 144                                    VALE=(
 145           9.39E-4,        185.E+6         *coef5,
 146           1.07E-3,        190.53E+6       *coef5,
 147           1.64E-3,        225.E+6         *coef5,
 148           3.32E-3,        261.E+6         *coef5,
 149           7.45E-3,        287.E+6         *coef5,
 150           0.011,          300.E+6         *coef5,
 151           0.032,          360.E+6         *coef5,
 152           0.05,           404.E+6         *coef5,
 153           0.1,            491.E+6         *coef5,
 154                                         ))
 155
 156 #
 157 # COURBE DE TRACTION à LA TEMPéRATURE 200°C
 158 #
 159 coef6=1.E ## UNIT -6
 160 _A6=DEFI_FONCTION(      NOM_PARA='EPSI',
 161                            PROL_DROITE='EXCLU',
 162                            PROL_GAUCHE='EXCLU',
 163                                    VALE=(
 164           7.61E-4        ,140.E+6         *coef6,
 165           1.E-3          ,150.E+6         *coef6,
 166           2.9E-3         ,168.E+6         *coef6,
 167           4.E-3          ,181.E+6         *coef6,
 168           7.E-3          ,195.E+6         *coef6,
 169           1.E-2          ,205.E+6         *coef6,
 170           0.034          ,276.E+6         *coef6,
 171           0.05           ,303.E+6         *coef6,
 172           0.15           ,450.E+6         *coef6,
 173           ))
 174
 175 #
 176 # COURBE DE TRACTION à LA TEMPéRATURE 350°C
 177 #
 178 coef7=1.E ## UNIT -6
 179 _A7=DEFI_FONCTION(      NOM_PARA='EPSI',
 180                            PROL_DROITE='EXCLU',
 181                            PROL_GAUCHE='EXCLU',
 182                                    VALE=(
 183           7.56E-4        ,130.E+6         *coef7,
 184           8.9E-4         ,136.E+6         *coef7,
 185           1.84E-3        ,145.E+6         *coef7,
 186           2.9E-3         ,151.E+6         *coef7,
 187           4.9E-3         ,160.E+6         *coef7,
 188           8.9E-3         ,174.E+6         *coef7,
 189           0.011          ,180.E+6         *coef7,
 190           0.051          ,261.E+6         *coef7,
 191           #0.01          ,353.E+6         *coef7,
 192           0.1            ,353.E+6         *coef7,
 193           ))
 194
 195 #
 196 # COURBE DE TRACTION à LA TEMPéRATURE 400°C
 197 #
 198 coef8=1.E ## UNIT -6
 199 _A8=DEFI_FONCTION(      NOM_PARA='EPSI',
 200                            PROL_DROITE='EXCLU',
 201                            PROL_GAUCHE='EXCLU',
 202                                    VALE=(
 203           7.14E-4        ,120.E+6         *coef8,
 204           1.E-3          ,125.E+6         *coef8,
 205           2.E-3          ,134.E+6         *coef8,
 206           3.E-3          ,141.E+6         *coef8,
 207           8.E-3          ,157.E+6         *coef8,
 208           19.E-3         ,185.E+6         *coef8,
 209           25.E-3         ,200.E+6         *coef8,
 210           0.05           ,240.E+6         *coef8,
 211           0.1            ,320.E+6         *coef8,
 212           ))
 213
 214 #
 215 # COURBE DE TRACTION à LA TEMPéRATURE 600°C
 216 #
 217 coef9=1.E ## UNIT -6
 218 _A9=DEFI_FONCTION(      NOM_PARA='EPSI',
 219                            PROL_DROITE='EXCLU',
 220                            PROL_GAUCHE='EXCLU',
 221                                    VALE=(
 222           5.94E-4        ,90.E+6          *coef9,
 223           1.E-3          ,104.E+6         *coef9,
 224           2.E-3          ,112.E+6         *coef9,
 225           3.E-3          ,119.E+6         *coef9,
 226           6.5E-3         ,131.E+6         *coef9,
 227           1.E-2          ,141.E+6         *coef9,
 228           2.14E-2        ,174.E+6         *coef9,
 229           0.05           ,224.E+6         *coef9,
 230           0.15           ,350.E+6         *coef9,
 231           ))
 232
 233 #
 234 # COURBE DE TRACTION à LA TEMPéRATURE 800°C
 235 #
 236 coefA=1.E ## UNIT -6
 237 _AA=DEFI_FONCTION(      NOM_PARA='EPSI',
 238                            PROL_DROITE='EXCLU',
 239                            PROL_GAUCHE='EXCLU',
 240                                    VALE=(
 241           5.77E-4        ,75.E+6          *coefA,
 242           1.5E-3         ,81.E+6          *coefA,
 243           3.E-3          ,92.5E+6         *coefA,
 244           4.E-3          ,95.5E+6         *coefA,
 245           6.E-3          ,104.E+6         *coefA,
 246           0.01           ,115.E+6         *coefA,
 247           0.0278         ,141.E+6         *coefA,
 248           0.05           ,159.E+6         *coefA,
 249           0.1            ,170.E+6         *coefA,
 250           ))
 251
 252 #
 253 # COURBE DE TRACTION
 254 #
 255
 256 _AB=DEFI_NAPPE(
 257                    NOM_PARA='TEMP',
 258                PROL_DROITE='EXCLU',
 259                PROL_GAUCHE='EXCLU',
 260                  PARA=(20.,200.,350.,400.,600.,800., ),
 261                  FONCTION=(_A5,_A6,_A7,_A8,_A9,_AA, ))
 262
 263 # COURBE DE FATIGUE DE WOHLER RCCM Z1.4.2
 264
 265 _AZ=DEFI_FONCTION(NOM_PARA='SIGM',
 266                      INTERPOL='LOG',
 267                      VALE=(
 268                      180.0E6,   1000000.0  ,
 269                      200.0E6,   500000.0  ,
 270                      230.0E6,   200000.0  ,
 271                      260.0E6,   100000.0  ,
 272                      295.0E6,   50000.0  ,
 273                      350.0E6,   20000.0  ,
 274                      405.0E6,   10000.0  ,
 275                      485.0E6,   5000.0  ,
 276                      615.0E6,   2000.0  ,
 277                      750.0E6,   1000.0  ,
 278                      940.0E6,   500.0  ,
 279                      1275.0E6,   200.0  ,
 280                      1655.0E6,   100.0  ,
 281                      2190.0E6,   50.0  ,
 282                      3240.0E6,   20.0  ,
 283                      4480.0E6,   10.0,
 284                           ),);
 285 _BZ=CALC_FONCTION(INVERSE=_F(FONCTION=_AZ,),);
 286 _CZ=CALC_FONCTION(COMB=_F(FONCTION = _BZ,
 287                          COEF     =  1.E ## UNIT -6
 288                          ),);
 289 DETRUIRE( CONCEPT =_F( NOM=_AZ), )
 290 _AZ=CALC_FONCTION(INVERSE=_F(FONCTION=_CZ,),);
 291
 292
 293
 294 # SM EN FONCTION DE LA TEMPERATURE
 295
 296 _BD=DEFI_FONCTION(
 297             NOM_PARA='TEMP',
 298          PROL_DROITE='EXCLU',
 299          PROL_GAUCHE='CONSTANT',
 300                 VALE=(
 301                      20.0,115.0E6,
 302                      50.0,115.0E6,
 303                      100.0,115.0E6,
 304                      150.0,115.0E6,
 305                      200.0,109.0E6,
 306                      250.0,103.0E6,
 307                      300.0,96.0E6,
 308                      340.0,94.0E6,
 309                      350.0,94.0E6,
 310                      ))
 311 _AD=CALC_FONCTION(COMB=_F(FONCTION = _BD,
 312                          COEF     =  1.E ## UNIT -6
 313                          ),);
 314
 315 # N_KE CONSTANTE EN FAIT
 316 _AE=DEFI_FONCTION(
 317             NOM_PARA='TEMP',
 318          PROL_DROITE='CONSTANT',
 319          PROL_GAUCHE='CONSTANT',VALE=(20.,0.3,100.,0.3),)
 320
 321 # M_KE CONSTANTE EN FAIT
 322 _AF=DEFI_FONCTION(
 323             NOM_PARA='TEMP',
 324          PROL_DROITE='CONSTANT',
 325          PROL_GAUCHE='CONSTANT',VALE=(20.,1.7,100.,1.7),)
 326
 327
 328
 329 MAT=DEFI_MATERIAU(
 330               THER_NL=_F(                ## SUBST THER
 331                RHO_CP = _A1,             ## EVAL _A1
 332                LAMBDA = _A0,             ## EVAL _A0
 333                ),
 334               ELAS_FO=_F(                ## SUBST ELAS
 335                E = _A2,                  ## EVAL _A2
 336                NU = _A3,                 ## EVAL _A3
 337                ALPHA = _A4,              ## EVAL _A4
 338                TEMP_DEF_ALPHA = 20.,     ## SUPPR
 339                ),
 340               TRACTION=_F( SIGM = _AB,    ),
 341
 342               FATIGUE=_F(WOHLER=_AZ,
 343                          E_REFE=1.79E11*1.E ## UNIT -6
 344                         ),
 345
 346               RCCM_FO=_F(                ## SUBST RCCM
 347                SM=_AD,                   ## EVAL _AD
 348                N_KE=_AE,                 ## EVAL _AE
 349                M_KE=_AF,                 ## EVAL _AF
 350                )
 351           )
 352
 353 #
 354