Télécharger la plaquette Cast3M


fr en

Télécharger lonca.eso

Retour à la liste

Numérotation des lignes : 

lonca

  1. C LONCA     SOURCE    CB215821  16/04/21    21:17:42     8920
  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=SQRT((XB-XA)*(XB-XA)+(YB-YA)*(YB-YA))
  77.         XBC=SQRT((XC-XB)*(XC-XB)+(YC-YB)*(YC-YB))
  78. *
  79.         IF(MELE.EQ.8) THEN
  80. *
  81. * ELEMENT LINEAIRE
  82. *
  83.           XLCE=SQRT(2.D0*XAB*XBC)
  84. *
  85.         ELSE
  86. *
  87. * ELEMENT QUADRATIQUE
  88. *
  89.           XLCE=SQRT(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=SQRT(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.  
  160.  
  161.  
  162.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales