Esope source files

   1 : C LONCA     SOURCE    CHAT      05/01/13    01:25:06     5004           
   2 : C LONCA     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE LONCA(EX,RB,ALPHA,EPSU,XE,NBNN,ENGF,MELE)
   4 : *
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : *
   8 :       DIMENSION XE(3,NBNN)
   9 : *
  10 : ***********************************************************************
  11 : *  CALCUL DE LA LONGUEUR CARACTERISTIQUE POUR LES TRI3 TRI6 QUA4 QUA8
  12 : *  SOIT LA RACINE CARREE PONDERRE DE LA SURFACE DE L ELEMENT.
  13 : *  COEFFICIENT DE PONDERATION :   1  ELEMENTS QUADRATIQUES  TRI6 QUA8
  14 : *                                ___
  15 : *                              \/ 2  ELEMENTS LINEAIRES     TRI3 QUA4
  16 : *
  17 : * ENTREES   EX    : MODULE DE YOUNG
  18 : *           RB    : RESISTANCE LIMITE EN COMPRESSION SIMPLE
  19 : *           ALPHA : RAPPORT ENTRE LIMITE EN TRACTION ET LIMITE EN
  20 : *                   COMPRESSION
  21 : *           EPSU  : DEFORMATION LIMITE INITIALE EN TRACTION
  22 : *           ENGF  : ENERGIE DE FISSURATION ( N/mm ou daN/cm )
  23 : *
  24 : * SORTIES   EPSU  : DEFORMATION ULTIME EN TRACTION EN FONCTION
  25 : *                   DE LA LONGUEUR CARACTERISTIQUE
  26 : **********************************************************************
  27 : *       
  28 : *      POIDS DANS LE CAS OU LA SURFACE EST CALCULEE EN UTILISANT
  29 : *      L INTEGRATION NUMERIQUE
  30 : *
  31 :        P555=0.5555555555555D0
  32 :        P888=0.8888888888888D0
  33 :        P125=0.1259391805448D0
  34 :        P132=0.1323941527885D0
  35 : *
  36 : *-----------------------------------------------------------------------
  37 : *     RECHERCHE DES COORDONNESS DANS XE
  38 : *-----------------------------------------------------------------------
  39 : *     ELEMENTS TRI3 ET QUA4
  40 : *
  41 :       IF(MELE.EQ.4.OR.MELE.EQ.8) THEN
  42 :          I1=1
  43 :          J1=2
  44 :          K1=3
  45 :       ELSE
  46 : *     ELEMENTS TRI6 ET QUA8
  47 :         IF(MELE.EQ.6.OR.MELE.EQ.10) THEN
  48 :            I1=1
  49 :            J1=3
  50 :            K1=5
  51 :         ELSE
  52 :           WRITE(*,*)'**ERREUR DANS LONCA!!! **ELEMENT NON DISPONIBLE**'
  53 :           STOP
  54 :         ENDIF
  55 :       ENDIF
  56 :       XA=XE(1,I1)
  57 :       YA=XE(2,I1)
  58 :       ZA=XE(3,I1)
  59 :       XB=XE(1,J1)
  60 :       YB=XE(2,J1)
  61 :       ZB=XE(3,J1)
  62 :       XC=XE(1,K1)
  63 :       YC=XE(2,K1)
  64 :       ZC=XE(3,K1)
  65 : *
  66 : *-----------------------------------------------------------------------
  67 : *     CALCUL DE LA LONGUEUR CARACTERISTIQUE : XLCE
  68 : *-----------------------------------------------------------------------
  69 : *
  70 :       IF(MELE.EQ.8.OR.MELE.EQ.10) THEN
  71 : *
  72 : * CAS DU QUA4 ET QUA8 : ATTENTION MAILLAGE PARFAITEMENT RECTANGULAIRE
  73 : *                       UTILISATION DE L OPERATEUR DALLE
  74 : *                       LE RECTANGLE EST SUPPOSE ETRE DANS LA PLAN (Z=0)
  75 : *
  76 :         XAB=DSQRT((XB-XA)*(XB-XA)+(YB-YA)*(YB-YA))
  77 :         XBC=DSQRT((XC-XB)*(XC-XB)+(YC-YB)*(YC-YB))
  78 : *
  79 :         IF(MELE.EQ.8) THEN
  80 : *
  81 : * ELEMENT LINEAIRE
  82 : *
  83 :           XLCE=DSQRT(2.D0*XAB*XBC)
  84 : *
  85 :         ELSE
  86 : *
  87 : * ELEMENT QUADRATIQUE
  88 : *
  89 :           XLCE=DSQRT(XAB*XBC)
  90 : *
  91 :         ENDIF
  92 :       ELSE
  93 : *
  94 : * CAS DU TRI3 ET TRI6 : LA SURFACE EST CALCULEE A PARTIR DU PRODUIT
  95 : *                       VECTORIEL DE DEUX COTES DU TRIANGLE
  96 : *                       LE TRIANGLE EST SUPPOSE ETRE DANS LA PLAN (Z=0)
  97 : *
  98 : * DEFINITION DES DEUX VECTEURS
  99 : *
 100 :         XAB=XB-XA
 101 :         YAB=YB-YA
 102 :         XAC=XC-XA
 103 :         YAC=YC-YA
 104 : *
 105 : * SURFACE = NORME DU PRODUIT VECTORIEL / 2.
 106 : *
 107 :           SURF=(XAB*YAC)-(YAB*XAC)
 108 : *
 109 :         IF(MELE.EQ.4) THEN
 110 : *
 111 : * ELEMENT LINEAIRE
 112 : *
 113 :           XLCE=DSQRT(2.D0)*SURF
 114 : *
 115 :         ELSE
 116 : *
 117 : * ELEMENT QUADRATIQUE
 118 : *
 119 :           XLCE=SURF
 120 : *
 121 :         ENDIF
 122 :       ENDIF
 123 : *
 124 : *
 125 :       GF0=ENGF
 126 :       FTG=ALPHA*RB
 127 : *
 128 : * LONGUEUR CARACTERISTIQUE CRITIQUE : softening lineaire.
 129 : *
 130 :       XLCC=2.D0*EX*GF0/(FTG**2.D0)
 131 : *
 132 : * CALCUL DE LA DEFORMATION ULTIME
 133 : *
 134 :       IF(XLCE.LE.XLCC) THEN
 135 :         EPGF=2.D0*GF0/(FTG*XLCE)
 136 :       ELSE
 137 : *
 138 : * LE CRITERE LIMITANT LA TAILLE DE L ELEMENT EST VIOLE
 139 : * DANS CE CAS IL FAUT REDUIRE LA RESISTANCE EN TRACTION
 140 : * POUR NE PAS AVOIR DE DEFORMATION ULTIME DONNANT UN
 141 : * SNAP BACK : voir these Feenstra pp. 26-27
 142 : * DANS NOTRE CAS - critere d Ottosen - ON NE PEUT TROP
 143 : * MODIFIER LA LIMITE EN TRACTION C EST A DIRE LE RAPPORT
 144 : * ALPHA. AINSI ON ADOPTE LE CALAGE SUIVANT :
 145 : *
 146 :         EPGF=2.D0*GF0/(FTG*XLCC)
 147 : *
 148 : * DANS CE CAS ON NE DESSIPE PAS L ENERGIE CONSTANTE ENGF
 149 : *
 150 :       ENDIF
 151 : *
 152 :       EPSU=EPGF  
 153 : *
 154 :  1000 CONTINUE
 155 : *
 156 :       RETURN
 157 :       END
 158 :  
 159 :