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    CB215821  26/08/24    21:15:17     12622          
  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 151 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.   151 CONTINUE
  41. C
  42.       DO 15 I=1,NSTRSS
  43.        SIGR(I)=SIG0(I)/((1.-ENDT)*(1.-ENDC))
  44.    15 CONTINUE
  45. C
  46.       DO 30 I=1,NSTRSS
  47.          SIGF(I)=SIGR(I)+DSIGE(I)
  48.    30 CONTINUE
  49. C
  50. C ***** TEST DES CRITERES ACTIVES ***************
  51. C
  52.       CALL DRUTRA(SIGF,SEQTE,BETINSA)
  53.       CALL DRUCOM(SIGF,SEQCE,BETINSA)
  54. C
  55.       FCRIT = SEQTE - SEQT
  56.       FCRIC = SEQCE - SEQC
  57. C
  58.       IF(FCRIT.LT.0.D0.AND.FCRIC.LT.0.D0) IZON=0
  59.       IF(FCRIT.GE.0.D0.AND.FCRIC.LT.0.D0) IZON=1
  60.       IF(FCRIT.LT.0.D0.AND.FCRIC.GE.0.D0) IZON=2
  61.       IF(FCRIT.GE.0.D0.AND.FCRIC.GE.0.D0) IZON=3
  62. C
  63.       IF(IZON.EQ.1) THEN
  64.         CALL ACTI1(SIGR,SIGF,D,NSTRSS,BETINSA)
  65.         CALL ENDAME(1,BETINSA)
  66.              GOTO 100
  67.       ENDIF
  68. C
  69.       IF(IZON.EQ.2) THEN
  70.         CALL ACTI2(SIGR,SIGF,D,NSTRSS,BETINSA)
  71.         CALL ENDAME(2,BETINSA)
  72.              GOTO 100
  73.       ENDIF
  74. C
  75.       IF(IZON.EQ.3) THEN
  76.         CALL ACTI3(SIGR,SIGF,D,NSTRSS,BETINSA)
  77.         CALL ENDAME(1,BETINSA)
  78.         CALL ENDAME(2,BETINSA)
  79.              GOTO 100
  80.       ENDIF
  81. C
  82.   100 CONTINUE
  83. C
  84.       DO 150 I=1,NSTRSS
  85.        SIGF(I)=SIGF(I)*((1.-ENDT)*(1.-ENDC))
  86.   150 CONTINUE
  87.  
  88.       RETURN
  89.       END
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  

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