Télécharger kdom3b.eso

Retour à la liste

Numérotation des lignes :

  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