Télécharger ottodd.eso

Retour à la liste

Numérotation des lignes :

ottodd
  1. C OTTODD SOURCE CHAT 05/01/13 02:07:21 5004
  2. SUBROUTINE OTTODD(D,YOUN,XNU,G,GS,XINVL,W)
  3. C=========================================================================
  4. C
  5. C ENTREES :
  6. C YOUN,XNU,G,GS,MFR
  7. C XINVL,W
  8. C
  9. C SORTIES :
  10. C D
  11. C
  12. C==========================================================================
  13. C
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16.  
  17. -INC PPARAM
  18. -INC CCOPTIO
  19. C
  20. PARAMETER (XZER=0.D0,UN=1.D0,DEUX=2.D0)
  21. C
  22. DIMENSION W(3),XINVL(3),FK(3),D(18),GS(3)
  23. *
  24. FK(1)=XINVL(1)*W(1)/GS(1)+XINVL(2)*W(2)/GS(2)
  25. FK(2)=XINVL(1)*W(1)/GS(1)+XINVL(3)*W(3)/GS(3)
  26. FK(3)=XINVL(2)*W(2)/GS(2)+XINVL(3)*W(3)/GS(3)
  27. *
  28. DET=(UN+XNU)*(UN+XNU)*(UN-DEUX*XNU)/(YOUN*YOUN*YOUN)
  29. *
  30. D(1)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
  31. D(2)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
  32. D(3)=(UN-XNU*XNU)/(YOUN*YOUN)/DET
  33. D(4)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
  34. D(5)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
  35. D(6)=(XNU*(UN+XNU)/(YOUN*YOUN))/DET
  36. D(7)=G/(UN+G*FK(1))
  37. D(8)=G/(UN+G*FK(2))
  38. D(9)=G/(UN+G*FK(3))
  39. D(10)=XZER
  40. D(11)=XZER
  41. D(12)=XZER
  42. D(13)=XZER
  43. D(14)=XZER
  44. D(15)=XZER
  45. D(16)=XZER
  46. D(17)=XZER
  47. D(18)=XZER
  48. *
  49. * IF(IIMPI.EQ.42) THEN
  50. * WRITE(IOIMP,77100) DET,D(1),D(6),D(5)
  51. *77100 FORMAT(2X,'OTTODD DET=',1PE12.5,2X,'D(1)=',1PE12.5,2X,
  52. * & 'D(6)=',1PE12.5,2X,'D(5)=',1PE12.5/)
  53. * ENDIF
  54. *
  55. RETURN
  56. END
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  

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