Télécharger lonca.eso

Retour à la liste

Numérotation des lignes :

  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