Télécharger kaxdis.eso

Retour à la liste

Numérotation des lignes :

  1. C KAXDIS SOURCE CB215821 16/04/21 21:17:28 8920
  2. SUBROUTINE KAXDIS(A1,A2,NS,KIMP,NG0,NP0,NG,NP)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C Include contenant quelques constantes dont XPI :
  6. -INC CCREEL
  7. C*********************************************************************
  8. C SP appele par KAXC ou KAXK en axisymetrique
  9. C determine si les elements 1 et 2 sont proches ou non
  10. C si oui: methode des Trapezes
  11. C si non: methode de Gauss
  12. C entree: "
  13. C A1 :COORDONNEES DES SOMMETS de L'ELEMENT 1 (reduites)
  14. C points I1,J1
  15. C A2 :id element 2
  16. C NS :DIMENSION (2)
  17. C NG0,NP0 :nombre de points d'integration (Gauss et Trapezes) lu
  18. C dans FFOR
  19. C resultat:
  20. C NG,NP :nombre de points d'integration (Gauss et Trapezes)
  21. C
  22. C*********************************************************************
  23. DIMENSION A1(NS,2),A2(NS,2),U(3)
  24.  
  25. C on recalcule les surfaces dans le plan r-z
  26.  
  27. CALL KNORM2(NS,A1,NS,U,S1)
  28. CALL KNORM2(NS,A2,NS,U,S2)
  29. IF (KIMP.GE.3) write(6,*) ' S1 S2 ',S1,S2
  30.  
  31. C on calcule le minimum des cordes croisees
  32.  
  33. C..I1J1
  34.  
  35. DR = A2(1,1)-A1(1,1)
  36. DZ = A2(2,1)-A1(2,1)
  37. XI1J1 = SQRT(DR*DR+DZ*DZ)
  38.  
  39. C..I1J2
  40.  
  41. DR = A2(1,1)-A1(1,2)
  42. DZ = A2(2,1)-A1(2,2)
  43. XI1J2=SQRT(DR*DR+DZ*DZ)
  44.  
  45. C..I2J1
  46.  
  47. DR = A2(1,2)-A1(1,1)
  48. DZ = A2(2,2)-A1(2,1)
  49. XI2J1= SQRT(DR*DR+DZ*DZ)
  50.  
  51. C..I2J2
  52.  
  53. DR = A2(1,2)-A1(1,2)
  54. DZ = A2(2,2)-A1(2,2)
  55. XI2J2 = SQRT(DR*DR+DZ*DZ)
  56.  
  57. DIST = DMIN1(XI1J1,XI2J1,XI1J2,XI2J2)
  58.  
  59. IF (KIMP.GE.3) write(6,*) ' DIST',DIST
  60.  
  61. D1=DIST/S1
  62. D2=DIST/S2
  63.  
  64. IF (KIMP.GE.3) write(6,*) ' D1 D2 ',D1,D2
  65.  
  66. C cas standard: au maximum 10 points de Gauss
  67. C
  68. C CRIT=0.05
  69. C 10/2015
  70. C augmentation du critere pour les faces proches qui etait
  71. C insuffisant pour les configurations de type face-a-face
  72. CRIT = 0.2
  73. C
  74. IF ((D1.GE.CRIT).AND.(D2.GE.CRIT)) THEN
  75. NG=NG0
  76. NP=0
  77. ELSE
  78. C cas des faces jointives ou tres proches:
  79. C integration par la methode des Trapezes
  80. C
  81. NG=0
  82. NP=NP0
  83. ENDIF
  84.  
  85. IF(KIMP.GE.3) write(6,*) ' kaxdis NG NP ',NG,NP
  86.  
  87. RETURN
  88. END
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  

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