Aster/Tests/torsion.comm

   1               DEBUT();
   2
   3
   4               INCLUDE(UNITE: 80);
   5
   6 TPS =  1.     ;
   7
   8
   9
  10 SY = 90. ;
  11
  12
  13
  14 A   =  1.     ;
  15 H   =  0.01   ;
  16 N   =  50     ;
  17 T   =  0.01   ;
  18 TD = DEFI_VALEUR(R8: EVAL(T*180./PI)) ;
  19 CT = DEFI_VALEUR(R8: EVAL(COS(T))) ;
  20 ST = DEFI_VALEUR(R8: EVAL(SIN(T))) ;
  21 MST= DEFI_VALEUR(R8: EVAL(-ST));
  22 RAT= DEFI_VALEUR(R8: EVAL(H/A));
  23
  24
  25 K_REF = DEFI_VALEUR( R8: EVAL(0.9845615*Q1) );
  26 EPS_Y = DEFI_VALEUR( R8: EVAL( (1./K_REF)**5. ));
  27 YOUNG = DEFI_VALEUR( R8: EVAL(SY/EPS_Y) );
  28
  29 !FORMULE( REEL: (
  30   RP(REEL:EPSI) = K_REF * SY * EPSI**0.2
  31   ) );
  32
  33 ECH_EPSI = DEFI_LIST_REEL (
  34               DEBUT: EPS_Y
  35               INTERVALLE: (
  36                 JUSQU_A: 3.
  37                 NOMBRE : 300
  38                 )
  39               ) ;
  40
  41 TRAC      = CALC_FONC_INTERP (
  42               FONCTION: RP
  43               LIST_PARA  :  ECH_EPSI
  44               PROL_GAUCHE: 'LINEAIRE'
  45               PROL_DROIT : 'LINEAIRE'
  46               ) ;
  47
  48
  49             PRE_GIBI();
  50  MA =       LIRE_MAILLAGE();
  51
  52 &MA =       DEFI_GROUP
  53               (
  54               MAILLAGE: MA
  55               CREA_GROUP_NO:
  56                 (
  57                 TOUT_GROUP_MA: 'OUI'
  58                 )
  59               ) ;
  60
  61
  62  CUIVRE   = DEFI_MATERIAU
  63               (
  64               ELAS :
  65                 (
  66                 E : YOUNG
  67                 NU: 0.3
  68                 )
  69               TRACTION : (
  70                 SIGM: TRAC
  71                 )
  72               NON_LOCAL: (
  73                 LONG_CARA: LC
  74                 )
  75               ) ;
  76
  77  CHMAT =    AFFE_MATERIAU
  78               (
  79               MAILLAGE: MA
  80               AFFE :
  81                 (
  82                 GROUP_MA: TRANCHE
  83                 MATER   : CUIVRE
  84                 )
  85               ) ;
  86
  87
  88  MO_MECA =   AFFE_MODELE
  89               (
  90               MAILLAGE: MA
  91               AFFE :
  92                 (
  93                 GROUP_MA     : TRANCHE
  94                 PHENOMENE    : 'MECANIQUE'
  95                 MODELISATION : '3D_SI'
  96                 )
  97               ) ;
  98
  99  MO_DELO = AFFE_MODELE (
 100              MAILLAGE: MA
 101              AFFE:
 102                (
 103                GROUP_MA: TRANCHE
 104                PHENOMENE: 'NON_LOCAL'
 105                MODELISATION: '3D'
 106                )
 107              ) ;
 108
 109
 110 !FORMULE ( REEL : ( UX_IMP(REEL: INST, REEL:X, REEL:Y) =
 111                    -INST*RAT*Y
 112           ));
 113 !FORMULE ( REEL : ( UY_IMP(REEL: INST, REEL:X, REEL:Y) =
 114                     INST*RAT*X
 115           ));
 116 !FORMULE ( REEL : ( UZ_IMP(REEL:X, REEL:Y) =  0. ));
 117
 118
 119  ANCRAGE   = AFFE_CHAR_MECA   (
 120               MODELE  : MO_MECA
 121               FACE_IMPO:
 122                 (
 123                 GROUP_MA: BAS
 124                 DX: 0.
 125                 DY: 0.
 126                 DZ: 0.
 127                 )
 128               ) ;
 129
 130  TORSION  = AFFE_CHAR_MECA_F (
 131               MODELE  : MO_MECA
 132               FACE_IMPO:
 133                 (
 134                 GROUP_MA: HAUT
 135                 DX      : UX_IMP
 136                 DY      : UY_IMP
 137                 DZ      : UZ_IMP
 138                 )
 139               ) ;
 140
 141
 142  SYMETRIE = AFFE_CHAR_MECA (
 143               MODELE: MO_MECA
 144               LIAISON_GROUP: (
 145                 GROUP_MA_1   : SYM_CY_0
 146                 GROUP_MA_2   : SYM_CY_1
 147                 SANS_GROUP_NO: (HAUT, BAS)
 148                 DDL_1        : ('DX' 'DY' )
 149                 DDL_2        : ('DX' )
 150                 COEF_MULT_1  : (CT MST)
 151                 COEF_MULT_2  : (-1.)
 152                 COEF_IMPO    : 0.
 153                 CENTRE       : (0. 0. 0.)
 154                 ANGL_NAUT    : (TD 0. 0.)
 155                 SOMMET       : 'OUI'
 156                 )
 157               LIAISON_GROUP: (
 158                 GROUP_MA_1   : SYM_CY_0
 159                 GROUP_MA_2   : SYM_CY_1
 160                 SANS_GROUP_NO: (HAUT, BAS)
 161                 DDL_1        : ('DX' 'DY' )
 162                 DDL_2        : ('DY' )
 163                 COEF_MULT_1  : (ST CT)
 164                 COEF_MULT_2  : (-1.)
 165                 COEF_IMPO    : 0.
 166                 CENTRE       : (0. 0. 0.)
 167                 ANGL_NAUT    : (TD 0. 0.)
 168                 SOMMET       : 'OUI'
 169                 )
 170               LIAISON_GROUP: (
 171                 GROUP_MA_1   : SYM_CY_0
 172                 GROUP_MA_2   : SYM_CY_1
 173                 SANS_GROUP_NO: (HAUT, BAS)
 174                 DDL_1        : ('DZ')
 175                 DDL_2        : ('DZ')
 176                 COEF_MULT_1  : (1.)
 177                 COEF_MULT_2  : (-1.)
 178                 COEF_IMPO    : 0.
 179                 CENTRE       : (0. 0. 0.)
 180                 ANGL_NAUT    : (TD 0. 0.)
 181                 SOMMET       : 'OUI'
 182                 )
 183               ) ;
 184
 185
 186  INSTANTS = DEFI_LIST_REEL
 187               (
 188               DEBUT: 0.
 189               INTERVALLE:
 190                 (JUSQU_A: TPS  NOMBRE: 50)
 191               ) ;
 192
 193
 194
 195  EVOL =     STAT_NON_LINE
 196               (
 197               MODELE: MO_MECA
 198               CHAM_MATER: CHMAT
 199               EXCIT :
 200                (CHARGE: TORSION)
 201                (CHARGE: ANCRAGE)
 202                (CHARGE: SYMETRIE)
 203               INCREMENT:
 204                 (
 205                 LIST_INST: INSTANTS
 206                 SUBD_PAS : 5
 207                 SUBD_PAS_MINI: 0.0001
 208                 )
 209               COMP_INCR:
 210                 (
 211                 RELATION: 'PLAS_GRAD_TRAC'
 212                 DEFORMATION: 'PETIT'
 213                 )
 214               NEWTON:
 215                 (
 216                 MATRICE:'ELASTIQUE'
 217                 )
 218               CONVERGENCE:
 219                 (
 220                 RESI_GLOB_RELA: 1.E-4
 221                 ITER_GLOB_MAXI: 500
 222                 ITER_INTE_MAXI: 250
 223                 RESI_INTE_RELA: 1.E-9
 224                 )
 225               ARCHIVAGE: (
 226                 ARCH_ETAT_INIT: 'OUI'
 227                 LIST_INST     : INSTANTS
 228                 )
 229
 230               MODELE_NON_LOCAL: MO_DELO
 231               LAGR_NON_LOCAL: (
 232                 RESI_PRIM_ABSO: 1.E-9
 233                 RESI_DUAL_ABSO: 1.E-4
 234                 RHO: 1000.
 235                 ITER_PRIM_MAXI: 100
 236                 ITER_DUAL_MAXI: 10
 237                 )
 238               ) ;
 239
 240
 241
 242 &EVOL       = CALC_ELEM (
 243                 MODELE: MO_MECA
 244                 CHAM_MATER: CHMAT
 245                 RESULTAT  : EVOL
 246                 OPTION: ('SIEF_ELNO_ELGA' 'VARI_ELNO_ELGA')
 247                 ) ;
 248
 249 &EVOL       = CALC_NO (
 250                 RESULTAT: EVOL
 251                 OPTION: ('FORC_NODA'
 252                          'SIEF_NOEU_ELGA' 'VARI_NOEU_ELGA')
 253                 ) ;
 254
 255
 256 ECHAN       = INTE_MAIL_3D (
 257                 MAILLAGE    : MA
 258                 DEFI_SEGMENT: (
 259                   GROUP_NO_ORIG: O1
 260                   GROUP_NO_EXTR: A1
 261                   )
 262                 ) ;
 263
 264  PLAS_CUM = POST_RELEVE_T (
 265               ACTION: (
 266                 INTITULE : 'P'
 267                 CHEMIN   :  ECHAN
 268                 RESULTAT :  EVOL
 269                 NOM_CHAM : 'VARI_NOEU_ELGA'
 270                 NOM_CMP  : 'V1'
 271                 INST     :  TPS
 272                 OPERATION: 'EXTRACTION'
 273                 )
 274               ) ;
 275
 276
 277  SIGMA    = POST_RELEVE_T (
 278               ACTION: (
 279                 INTITULE : 'SIYZ'
 280                 CHEMIN   :  ECHAN
 281                 RESULTAT :  EVOL
 282                 NOM_CHAM : 'SIEF_NOEU_ELGA'
 283                 NOM_CMP  : 'SIYZ'
 284                 INST     :  TPS
 285                 OPERATION: 'EXTRACTION'
 286                 )
 287               ) ;
 288
 289
 290
 291  MOM_A   = POST_RELEVE_T (
 292               ACTION: (
 293                 INTITULE  : 'MOMENT'
 294                 GROUP_NO  :  HAUT
 295                 RESULTAT  :  EVOL
 296                 NOM_CHAM  : 'FORC_NODA'
 297                 RESULTANTE: ('DX' 'DY' 'DZ')
 298                 MOMENT    : ('DRX' 'DRY' 'DRZ')
 299                 POINT     : (0. 0. H)
 300                 TOUT_ORDRE: 'OUI'
 301                 OPERATION : 'EXTRACTION'
 302                 )
 303               ) ;
 304
 305  MOM_KA_A   = RECU_FONCTION (
 306               TABLE : MOM_A
 307               PARA_X: 'INST'
 308               PARA_Y: 'MOMENT_Z'
 309               ) ;
 310
 311 !FORMULE ( REEL:( MA3_KA(REEL: INST) = MOM_KA_A(INST)*360/TD/A**3/SY));
 312
 313
 314
 315               DEFUFI(
 316                 IMPRESSION: (
 317                   NOM: 'POST_PY'
 318                   UNITE: 81
 319                 )
 320               ) ;
 321
 322
 323
 324  IMPR_COURBE(
 325    FORMAT: 'AGRAF'
 326    FICHIER: 'POST_PY'
 327    COURBE: (
 328      FONCTION : MA3_KA
 329      LIST_PARA: INSTANTS
 330      )
 331    ) ;
 332
 333
 334               FIN();
 335
 336
 337 IMPR_TABLE(TABLE:SIGMA);
 338 IMPR_TABLE(TABLE:PLAS_CUM);
 339
 340             FIN();