Télécharger la plaquette Cast3M


fr en

Télécharger kaxdis.eso

Retour à la liste

Numérotation des lignes : 

kaxdis

  1. C KAXDIS    SOURCE    CB215821  17/11/30    21:16:40     9639           
  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 =MIN(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.  
  95.  

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