Télécharger la plaquette Cast3M


fr en

Télécharger bone3d.eso

Retour à la liste

Numérotation des lignes : 

bone3d

  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