Télécharger pamaot.eso

Retour à la liste

Numérotation des lignes :

  1. C PAMAOT SOURCE CHAT 05/01/13 02:08:54 5004
  2. C PAMAOT SOURCE INSL 24/10/96
  3. SUBROUTINE PAMAOT(ALPHA,A1,B1,DK1,DK2)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. C **********************************************
  8. C * CALCUL DES PARAMETRES DU MATERIAU *
  9. C * ( MODEL D'OTTOSEN ) *
  10. C **********************************************
  11. C-------------------------------------------------------------------
  12. DIMENSION D(3,3),D1(3,3),AA(3),BB(3),XK1(3),XK2(3),CC(3),AL(3)
  13. DIMENSION ALI(2)
  14. C-------------------------------------------------------------------
  15. AA(1) =1.8076D0
  16. BB(1) =4.0962D0
  17. XK1(1)=14.4863D0
  18. XK2(1)=0.9914D0
  19. AL(1) =0.08D0
  20. *
  21. AA(2) =1.2759D0
  22. BB(2) =3.1962D0
  23. XK1(2)=11.7365D0
  24. XK2(2)=0.9801D0
  25. AL(2) =0.1D0
  26. *
  27. AA(3) =0.9218D0
  28. BB(3) =2.5969D0
  29. XK1(3)=9.911D0
  30. XK2(3)=0.9647D0
  31. AL(3) =0.12D0
  32. C-------------------------------------------------------------------
  33. ALI(1)=(AL(1)+AL(2))/2.D0
  34. ALI(2)=(AL(2)+AL(3))/2.D0
  35. C-------------------------------------------------------------------
  36. IF(ALPHA.LE.ALI(1)) THEN
  37. A1 =AA(1)
  38. B1 =BB(1)
  39. DK1=XK1(1)
  40. DK2=XK2(1)
  41. GOTO 200
  42. ENDIF
  43. IF(ALPHA.LE.ALI(2).AND.ALPHA.GT.ALI(1)) THEN
  44. A1 =AA(2)
  45. B1 =BB(2)
  46. DK1=XK1(2)
  47. DK2=XK2(2)
  48. GOTO 200
  49. ENDIF
  50. IF(ALPHA.GT.ALI(2)) THEN
  51. A1 =AA(3)
  52. B1 =BB(3)
  53. DK1=XK1(3)
  54. DK2=XK2(3)
  55. GOTO 200
  56. ENDIF
  57. C-------------------------------------------------------------------
  58. D(1,1)=AL(1)*AL(1)
  59. D(1,2)=AL(1)
  60. D(1,3)=1.D0
  61. D(2,1)=AL(2)*AL(2)
  62. D(2,2)=AL(2)
  63. D(2,3)=1.D0
  64. D(3,1)=AL(3)*AL(3)
  65. D(3,2)=AL(3)
  66. D(3,3)=1.D0
  67. *
  68. CALL INVER3(D,D1)
  69. CALL BST(D1,AA,3,3,CC)
  70. A1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  71. CALL ZERO(CC,3,1)
  72. CALL BST(D1,BB,3,3,CC)
  73. B1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  74. CALL ZERO(CC,3,1)
  75. CALL BST(D1,XK1,3,3,CC)
  76. DK1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  77. CALL ZERO(CC,3,1)
  78. CALL BST(D1,XK2,3,3,CC)
  79. DK2=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  80. C-------------------------------------------------------------------
  81. 200 CONTINUE
  82. C-------------------------------------------------------------------
  83. RETURN
  84. END
  85.  
  86.  
  87.  

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