DEBUT();


              INCLUDE(UNITE: 80);
 
TPS =  1.     ;



SY = 90. ;



A   =  1.     ;
H   =  0.01   ;
N   =  50     ;
T   =  0.01   ;
TD = DEFI_VALEUR(R8: EVAL(T*180./PI)) ;
CT = DEFI_VALEUR(R8: EVAL(COS(T))) ;
ST = DEFI_VALEUR(R8: EVAL(SIN(T))) ;     
MST= DEFI_VALEUR(R8: EVAL(-ST));
RAT= DEFI_VALEUR(R8: EVAL(H/A));


K_REF = DEFI_VALEUR( R8: EVAL(0.9845615*Q1) );
EPS_Y = DEFI_VALEUR( R8: EVAL( (1./K_REF)**5. ));
YOUNG = DEFI_VALEUR( R8: EVAL(SY/EPS_Y) );

!FORMULE( REEL: (
  RP(REEL:EPSI) = K_REF * SY * EPSI**0.2
  ) );
  
ECH_EPSI = DEFI_LIST_REEL (
              DEBUT: EPS_Y
              INTERVALLE: (
                JUSQU_A: 3.
                NOMBRE : 300
                )
              ) ;
  
TRAC      = CALC_FONC_INTERP (
              FONCTION: RP
              LIST_PARA  :  ECH_EPSI
              PROL_GAUCHE: 'LINEAIRE'
              PROL_DROIT : 'LINEAIRE'
              ) ;
              
              
            PRE_GIBI();
 MA =       LIRE_MAILLAGE();

&MA =       DEFI_GROUP
              (
              MAILLAGE: MA
              CREA_GROUP_NO:
                (
                TOUT_GROUP_MA: 'OUI'
                )
              ) ;
              

 CUIVRE   = DEFI_MATERIAU
              (
              ELAS :
                (
                E : YOUNG               
                NU: 0.3             
                )
              TRACTION : (
                SIGM: TRAC
                )
              NON_LOCAL: (
                LONG_CARA: LC
                )
              ) ;

 CHMAT =    AFFE_MATERIAU
              (
              MAILLAGE: MA
              AFFE :
                (
                GROUP_MA: TRANCHE
                MATER   : CUIVRE
                )
              ) ;
              

 MO_MECA =   AFFE_MODELE
              (
              MAILLAGE: MA
              AFFE :
                (
                GROUP_MA     : TRANCHE
                PHENOMENE    : 'MECANIQUE'
                MODELISATION : '3D_SI'
                )
              ) ;

 MO_DELO = AFFE_MODELE (
             MAILLAGE: MA
             AFFE: 
               (
               GROUP_MA: TRANCHE
               PHENOMENE: 'NON_LOCAL'
               MODELISATION: '3D'
               )
             ) ;
             

!FORMULE ( REEL : ( UX_IMP(REEL: INST, REEL:X, REEL:Y) =
                   -INST*RAT*Y
          ));
!FORMULE ( REEL : ( UY_IMP(REEL: INST, REEL:X, REEL:Y) = 
                    INST*RAT*X
          ));
!FORMULE ( REEL : ( UZ_IMP(REEL:X, REEL:Y) =  0. ));


 ANCRAGE   = AFFE_CHAR_MECA   (
              MODELE  : MO_MECA
              FACE_IMPO:
                (
                GROUP_MA: BAS
                DX: 0.
                DY: 0.
                DZ: 0.
                )
              ) ;

 TORSION  = AFFE_CHAR_MECA_F (
              MODELE  : MO_MECA
              FACE_IMPO:
                (
                GROUP_MA: HAUT
                DX      : UX_IMP
                DY      : UY_IMP
                DZ      : UZ_IMP
                )
              ) ;
                

 SYMETRIE = AFFE_CHAR_MECA (
              MODELE: MO_MECA
              LIAISON_GROUP: (
                GROUP_MA_1   : SYM_CY_0
                GROUP_MA_2   : SYM_CY_1
                SANS_GROUP_NO: (HAUT, BAS)
                DDL_1        : ('DX' 'DY' )
                DDL_2        : ('DX' )
                COEF_MULT_1  : (CT MST)
                COEF_MULT_2  : (-1.)
                COEF_IMPO    : 0.
                CENTRE       : (0. 0. 0.)
                ANGL_NAUT    : (TD 0. 0.)
                SOMMET       : 'OUI'
                )
              LIAISON_GROUP: (
                GROUP_MA_1   : SYM_CY_0
                GROUP_MA_2   : SYM_CY_1
                SANS_GROUP_NO: (HAUT, BAS)
                DDL_1        : ('DX' 'DY' )
                DDL_2        : ('DY' )
                COEF_MULT_1  : (ST CT)
                COEF_MULT_2  : (-1.)
                COEF_IMPO    : 0.
                CENTRE       : (0. 0. 0.)
                ANGL_NAUT    : (TD 0. 0.)
                SOMMET       : 'OUI'
                )
              LIAISON_GROUP: (
                GROUP_MA_1   : SYM_CY_0
                GROUP_MA_2   : SYM_CY_1
                SANS_GROUP_NO: (HAUT, BAS)
                DDL_1        : ('DZ')
                DDL_2        : ('DZ')
                COEF_MULT_1  : (1.)
                COEF_MULT_2  : (-1.)
                COEF_IMPO    : 0.
                CENTRE       : (0. 0. 0.)
                ANGL_NAUT    : (TD 0. 0.)
                SOMMET       : 'OUI'
                )
              ) ;
              
            
 INSTANTS = DEFI_LIST_REEL
              (
              DEBUT: 0.
              INTERVALLE:
                (JUSQU_A: TPS  NOMBRE: 50)
              ) ;



 EVOL =     STAT_NON_LINE
              (
              MODELE: MO_MECA
              CHAM_MATER: CHMAT
              EXCIT :
               (CHARGE: TORSION)
               (CHARGE: ANCRAGE)
               (CHARGE: SYMETRIE)                
              INCREMENT:
                (
                LIST_INST: INSTANTS
                SUBD_PAS : 5
                SUBD_PAS_MINI: 0.0001
                )
              COMP_INCR:
                (
                RELATION: 'PLAS_GRAD_TRAC'
                DEFORMATION: 'PETIT'
                )
              NEWTON:
                (
                MATRICE:'ELASTIQUE'
                )
              CONVERGENCE:
                (
                RESI_GLOB_RELA: 1.E-4
                ITER_GLOB_MAXI: 500
                ITER_INTE_MAXI: 250
                RESI_INTE_RELA: 1.E-9
                )
              ARCHIVAGE: (
                ARCH_ETAT_INIT: 'OUI'
                LIST_INST     : INSTANTS
                )

              MODELE_NON_LOCAL: MO_DELO
              LAGR_NON_LOCAL: (
                RESI_PRIM_ABSO: 1.E-9
                RESI_DUAL_ABSO: 1.E-4
                RHO: 1000.
                ITER_PRIM_MAXI: 100
                ITER_DUAL_MAXI: 10
                )
              ) ;

          

&EVOL       = CALC_ELEM (
                MODELE: MO_MECA
                CHAM_MATER: CHMAT
                RESULTAT  : EVOL
                OPTION: ('SIEF_ELNO_ELGA' 'VARI_ELNO_ELGA')
                ) ;

&EVOL       = CALC_NO (
                RESULTAT: EVOL
                OPTION: ('FORC_NODA' 
                         'SIEF_NOEU_ELGA' 'VARI_NOEU_ELGA')
                ) ;

  
ECHAN       = INTE_MAIL_3D (
                MAILLAGE    : MA
                DEFI_SEGMENT: (
                  GROUP_NO_ORIG: O1
                  GROUP_NO_EXTR: A1
                  )
                ) ;
                
 PLAS_CUM = POST_RELEVE_T (
              ACTION: (
                INTITULE : 'P'
                CHEMIN   :  ECHAN
                RESULTAT :  EVOL
                NOM_CHAM : 'VARI_NOEU_ELGA'
                NOM_CMP  : 'V1'
                INST     :  TPS
                OPERATION: 'EXTRACTION'
                )
              ) ;              


 SIGMA    = POST_RELEVE_T (
              ACTION: (
                INTITULE : 'SIYZ'
                CHEMIN   :  ECHAN
                RESULTAT :  EVOL
                NOM_CHAM : 'SIEF_NOEU_ELGA'
                NOM_CMP  : 'SIYZ'
                INST     :  TPS
                OPERATION: 'EXTRACTION'
                )
              ) ;              



 MOM_A   = POST_RELEVE_T (
              ACTION: (
                INTITULE  : 'MOMENT'
                GROUP_NO  :  HAUT
                RESULTAT  :  EVOL
                NOM_CHAM  : 'FORC_NODA'
                RESULTANTE: ('DX' 'DY' 'DZ')
                MOMENT    : ('DRX' 'DRY' 'DRZ')
                POINT     : (0. 0. H)
                TOUT_ORDRE: 'OUI'
                OPERATION : 'EXTRACTION'
                )
              ) ;

 MOM_KA_A   = RECU_FONCTION (
              TABLE : MOM_A
              PARA_X: 'INST'
              PARA_Y: 'MOMENT_Z'
              ) ;

!FORMULE ( REEL:( MA3_KA(REEL: INST) = MOM_KA_A(INST)*360/TD/A**3/SY));



              DEFUFI(
                IMPRESSION: (
                  NOM: 'POST_PY'
                  UNITE: 81
                )
              ) ;
              
              
              
 IMPR_COURBE( 
   FORMAT: 'AGRAF'
   FICHIER: 'POST_PY'
   COURBE: (
     FONCTION : MA3_KA
     LIST_PARA: INSTANTS
     )
   ) ;


              FIN();
              
              
IMPR_TABLE(TABLE:SIGMA);
IMPR_TABLE(TABLE:PLAS_CUM);
 
            FIN();