Télécharger la plaquette Cast3M


fr en

Télécharger kdom4b.eso

Retour à la liste

Numérotation des lignes : 

kdom4b

  1. C KDOM4B    SOURCE    CHAT      05/01/13    00:54:16     5004
  2.       SUBROUTINE KDOM4B(X1,Y1,X2,Y2,X3,Y3,VOLU,SURF,XCEN,YCEN,XC2D,YC2D)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  KDOM4B
  9. C
  10. C DESCRIPTION       :  Subroutine called by KDOM4A
  11. C                      Axial-symmetric case
  12. C                      We compute
  13. C                      the surface of the 2D mesh
  14. C                      the 3D volume of the 3D mesh
  15. C                      its center of gravity
  16. C
  17. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  18. C
  19. C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF
  20. C
  21. C************************************************************************
  22. C
  23. C INPUT             : (X1,Y1), (X2,Y2), (X3,Y3) coordinates of the TRI3
  24. C
  25. C OUTPUT            : VOLU: volume
  26. c                     SURF: surface
  27. C                     (XCEN,YCEN) center of gravity
  28. C                     (XC2D,YC2D) center of gravity of the surface
  29. C
  30. C************************************************************************
  31. C
  32. C Created the 24/02/04
  33. C
  34. C
  35.       IMPLICIT INTEGER(I-N)
  36.       REAL*8 X1,Y1,X2,Y2,X3,Y3,VOLU,SURF,XCEN,YCEN
  37.      &     ,RJACO,XC2D,YC2D
  38. C
  39.       RJACO=(X2-X1)*(Y3-Y1)-(X3-X1)*(Y2-Y1)
  40.       RJACO=ABS(RJACO)
  41.       SURF=0.5D0*RJACO
  42.       XC2D=(X1+X2+X3)/3.0D0
  43.       YC2D=(Y1+Y2+Y3)/3.0D0
  44.       VOLU=RJACO*(X1+X2+X3)/6.0D0
  45. C
  46.       XCEN=(X1+X2)**2
  47.       XCEN=XCEN+((X1+X3)**2)
  48.       XCEN=XCEN+((X2+X3)**2)
  49.       XCEN=XCEN*RJACO/24.0D0
  50.       XCEN=XCEN/VOLU
  51.       YCEN=(X1+X2)*(Y1+Y2)
  52.       YCEN=YCEN+((X1+X3)*(Y1+Y3))
  53.       YCEN=YCEN+((X2+X3)*(Y2+Y3))
  54.       YCEN=YCEN*RJACO/24.0D0
  55.       YCEN=YCEN/VOLU
  56. C
  57.       RETURN
  58.       END
  59.  
  60.  
  61.  
  62.  
  63.  

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