      SUBROUTINE ETC3R1D(COOR,CAR,IOPT,SIGMA,U)
C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  BUT   : CALCUL DE LA CONTRAINTE
C  ---     pentaedre r1 DROIT
c          formule d'integration a 6 points (gauss) 3*2
C
C  PARAMETRES D ENTREE  :
C  -------------------
C   X,Y,Z   : TABLEAUX DES COORDONNEES DES POINTS DE L ELEMENT
c   ijt     : permutation pour oasser de la numerotation par inconnues
c             a celle par noeuds
c   nno     : nombre de noeuds de l'element
c   npo     : nombre de points
c   npi     : nombre de points d'integration
c   dp      : valeur des derivees des polynomes de base aux points
c             d'integration sur l'element de reference
c   vp1     : valeur des polynomes de base aux points
c             d'integration sur l'element de reference
c   poids   : poids de la formule d'integration
c
c  tableaux de travail :
c  -------------------
c   delta   : jacobien aux points d'integration
C
C  PARAMETRE DE SORTIE  :
C  --------------------
C   sigma     : 
C            
C
C +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  PROGRAMMEURS  : Marina Vidrascu 2003
C ...................................................................
C
      parameter (ndim=3, nno=6, npo=6)
      DOUBLE PRECISION SIGMAE(6,NDIM*NNO,NNO),sigma(6,nno),U(ndim*npo)
      DOUBLE PRECISION COOR(npo,3)
      DIMENSION IJT(ndim*nno)
      DOUBLE PRECISION  elas(9),delta(nno), a2(ndim,nno,nno),
     &                  DP(ndim,nno,nno),car(iopt)
c
      DATA IJT/1,4,7,10,13,16,
     +         2,5,8,11,14,17,
     +         3,6,9,12,15,18/
C
      DATA DP/ -1.0 ,-1.0 ,-1.0 ,1.0  ,  0.0  , 0.0  ,0.0 , 1.0  , 0.0 ,
     *     0.0 ,  0.0 ,  1.0 ,   0.0 ,  0.0 ,  0.0 , 0.0 ,  0.0 ,  0.0 ,
     *    -1.0 , -1.0 ,  0.0 ,   1.0 ,  0.0 , -1.0 , 0.0 ,  1.0 ,  0.0 ,
     *     0.0 ,  0.0 ,  0.0 ,   0.0 ,  0.0 ,  1.0 , 0.0 ,  0.0 ,  0.0 ,
     *    -1.0 , -1.0 ,  0.0 ,   1.0 ,  0.0 ,  0.0 , 0.0 ,  1.0 , -1.0 ,
     *     0.0 ,  0.0 ,  0.0 ,   0.0 ,  0.0 ,  0.0 , 0.0 ,  0.0 ,  1.0 ,
     *     0.0 ,  0.0 , -1.0 ,   0.0 ,  0.0 ,  0.0 , 0.0 ,  0.0 ,  0.0 ,
     *    -1.0 , -1.0 ,  1.0 ,   1.0 ,  0.0 ,  0.0 , 0.0 ,  1.0 ,  0.0 ,
     *     0.0 ,  0.0 ,  0.0 ,   0.0 ,  0.0 , -1.0 , 0.0 ,  0.0 ,  0.0 ,
     *    -1.0 , -1.0 ,  0.0 ,   1.0 ,  0.0 ,  1.0 , 0.0 ,  1.0 ,  0.0 ,
     *     0.0 ,  0.0 ,  0.0 ,   0.0 ,  0.0 ,  0.0 , 0.0 ,  0.0 , -1.0 ,
     *    -1.0 , -1.0 ,  0.0 ,   1.0 ,  0.0 ,  0.0 , 0.0 ,  1.0 ,  1.0/


c     si isotrope e=car(1), nu=car(2)
c     si orthotrope elas=car
      if (iopt .eq. 2) then
C
C         CAS ISOTROPE
C         ------------
          CALL EC3C2C (NNO,NNO,COOR(1,1),COOR(1,2),COOR(1,3),IJT,
     +            dp,dp,iopt,car(1),car(2),elas,sigmae,sigma,u,delta,a2)
      else
C
C         CAS ORTHOTROPE
C         --------------
          CALL EC3C2C (NNO,NNO,COOR(1,1),COOR(1,2),COOR(1,3),IJT,
     +            dp,dp,iopt,car(1),car(2),car,sigmae,sigma,u,delta,a2)


      endif


      END
