Télécharger bone3d.eso

Retour à la liste

Numérotation des lignes :

  1. C BONE3D SOURCE CHAT 05/01/12 21:42:02 5004
  2. SUBROUTINE BONE3D(SIG0,SIGF,DEPST,NSTRSS,BETINSA)
  3. C
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. DIMENSION SIG0(NSTRSS),SIGF(NSTRSS),DEPST(NSTRSS)
  7. DIMENSION DSIGE(6),D(6,6),SIGR(6)
  8. C
  9. SEGMENT BETINSA
  10. REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR
  11. REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  12. INTEGER IFIS,IPLA,IBB,IGAU
  13. ENDSEGMENT
  14. C
  15. C ***** CALCUL DU TIR ELASTIQUE ***************
  16. C
  17. C
  18. C COMPORTEMENT TRIDIMENSIONNEL
  19. C
  20. CALL ZERO(D,6,6)
  21. E1=YOUN/(1.D0+XNU)/(1.D0-2.D0*XNU)
  22. D(1,1)=E1*(1.D0-XNU)
  23. D(1,2)=E1*XNU
  24. D(1,3)=D(1,2)
  25. D(2,1)=D(1,2)
  26. D(2,2)=D(1,1)
  27. D(2,3)=D(1,3)
  28. D(3,1)=D(1,2)
  29. D(3,2)=D(1,2)
  30. D(3,3)=D(1,1)
  31. D(4,4)=E1*(1.D0-2.D0*XNU)/2.D0
  32. D(5,5)=D(4,4)
  33. D(6,6)=D(5,5)
  34. C
  35. DO 10 I=1,NSTRSS
  36. DSIGE(I)=0.D0
  37. DO 10 J=1,NSTRSS
  38. DSIGE(I)=DSIGE(I)+D(I,J)*DEPST(J)
  39. 10 CONTINUE
  40. C
  41. DO 15 I=1,NSTRSS
  42. SIGR(I)=SIG0(I)/((1.-ENDT)*(1.-ENDC))
  43. 15 CONTINUE
  44. C
  45. DO 30 I=1,NSTRSS
  46. SIGF(I)=SIGR(I)+DSIGE(I)
  47. 30 CONTINUE
  48. C
  49. C ***** TEST DES CRITERES ACTIVES ***************
  50. C
  51. CALL DRUTRA(SIGF,SEQTE,BETINSA)
  52. CALL DRUCOM(SIGF,SEQCE,BETINSA)
  53. C
  54. FCRIT = SEQTE - SEQT
  55. FCRIC = SEQCE - SEQC
  56. C
  57. IF(FCRIT.LT.0.D0.AND.FCRIC.LT.0.D0) IZON=0
  58. IF(FCRIT.GE.0.D0.AND.FCRIC.LT.0.D0) IZON=1
  59. IF(FCRIT.LT.0.D0.AND.FCRIC.GE.0.D0) IZON=2
  60. IF(FCRIT.GE.0.D0.AND.FCRIC.GE.0.D0) IZON=3
  61. C
  62. IF(IZON.EQ.1) THEN
  63. CALL ACTI1(SIGR,SIGF,D,NSTRSS,BETINSA)
  64. CALL ENDAME(1,BETINSA)
  65. GOTO 100
  66. ENDIF
  67. C
  68. IF(IZON.EQ.2) THEN
  69. CALL ACTI2(SIGR,SIGF,D,NSTRSS,BETINSA)
  70. CALL ENDAME(2,BETINSA)
  71. GOTO 100
  72. ENDIF
  73. C
  74. IF(IZON.EQ.3) THEN
  75. CALL ACTI3(SIGR,SIGF,D,NSTRSS,BETINSA)
  76. CALL ENDAME(1,BETINSA)
  77. CALL ENDAME(2,BETINSA)
  78. GOTO 100
  79. ENDIF
  80. C
  81. 100 CONTINUE
  82. C
  83. DO 150 I=1,NSTRSS
  84. SIGF(I)=SIGF(I)*((1.-ENDT)*(1.-ENDC))
  85. 150 CONTINUE
  86.  
  87. RETURN
  88. END
  89.  
  90.  
  91.  
  92.  
  93.  

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