Télécharger kdom4b.eso

Retour à la liste

Numérotation des lignes :

  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