Télécharger la plaquette Cast3M


fr en

Télécharger kdom3b.eso

Retour à la liste

Numérotation des lignes : 

kdom3b

  1. C KDOM3B    SOURCE    CHAT      05/01/13    00:54:03     5004
  2.       SUBROUTINE KDOM3B(X1,Y1,X3,Y3,VOLU,X2,Y2)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  KDOM3B
  9. C
  10. C DESCRIPTION       :  Subroutine called by KDOM3A and KDOM4C
  11. C                      Axial-symmetric case
  12. C                      We compute
  13. C                      the surface of the 2D line
  14. C                      its center of gravity
  15. C
  16. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  17. C
  18. C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF
  19. C
  20. C************************************************************************
  21. C
  22. C INPUT             : (X1,Y1), (X3,Y3) coordinates of the line extrema
  23. C
  24. C OUTPUT            : VOLUME: surface
  25. C                     (X2,Y2) center of gravity
  26. C
  27. C************************************************************************
  28. C
  29. C Created the 24/02/04
  30. C
  31. C
  32.       IMPLICIT INTEGER(I-N)
  33.       REAL*8 X1,Y1,X3,Y3,X2,Y2,RL13,VOLU,W1,W2
  34. C
  35.       W1=1.0D0-(1.0D0/(3.0D0**0.5D0))
  36.       W2=1.0D0+(1.0D0/(3.0D0**0.5D0))
  37. C
  38.       RL13 = ((Y3-Y1)*(Y3-Y1)) + ((X3-X1)*(X3-X1))
  39.       RL13 = RL13**0.5D0
  40.       VOLU=0.5D0*(X1+X3)*RL13
  41.       X2=0.125D0*RL13*(((X1*W1+X3*W2)**2)+((X1*W2+X3*W1)**2))
  42.       X2=X2/VOLU
  43.       Y2=0.125D0*RL13*(((X1*W1+X3*W2)*(Y1*W1+Y3*W2))+
  44.      &     ((X3*W1+X1*W2)*(Y3*W1+Y1*W2)))
  45.       Y2=Y2/VOLU
  46.       RETURN
  47.       END
  48.  
  49.  
  50.  
  51.  
  52.  

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