Télécharger diacl.eso

Retour à la liste

Numérotation des lignes :

  1. C DIACL SOURCE CB215821 16/04/21 21:16:28 8920
  2. SUBROUTINE DIACL(EPS,SIN,COS)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. C EPS = TABLEAU DES DEFORMATIONS (EPS11,EPS22,SQRT(2)*EPS12)
  7. C EPSL1 = DEFORMATIONS PRINCIPALES
  8. C SIN ET COS = SINUS ET COSINUS DE L"ANGLE DE ROTATION
  9. DIMENSION EPS(3)
  10. C
  11. EPS12=EPS(3)/2.D0
  12. X1=2.D0*EPS12
  13. X2=EPS(1)-EPS(2)
  14. X3=SQRT(X1**2+X2**2)
  15. X8=X2*1.D+10
  16. IF (ABS(X1).GT.ABS(X8).OR.X2.EQ.0.D0) THEN
  17. COS=0.5D0*SQRT(2.D0)
  18. SIN=SIGN(COS,X1)
  19. ELSE
  20. X4=X1/X2
  21. X5=1.D0+X4**2
  22. X5=SIGN(1.D0,X2)/SQRT(X5)
  23. X6=(1.D0+X5)*0.5D0
  24. X6=SQRT(X6)
  25. X7=(1.D0-X5)*0.5D0
  26. X7=SIGN(1.D0,X1)*SQRT(X7)
  27. COS=X6
  28. SIN=X7
  29. END IF
  30. C
  31. RETURN
  32. END
  33.  
  34.  
  35.  
  36.  

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