Esope source files

   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
  19 :          CALL ZERO(D,6,6)
  20 :          E1=YOUN/(1.D0+XNU)/(1.D0-2.D0*XNU)
  21 :          D(1,1)=E1*(1.D0-XNU)
  22 :          D(1,2)=E1*XNU
  23 :          D(1,3)=D(1,2)
  24 :          D(2,1)=D(1,2)
  25 :          D(2,2)=D(1,1)
  26 :          D(2,3)=D(1,3)
  27 :          D(3,1)=D(1,2)
  28 :          D(3,2)=D(1,2)
  29 :          D(3,3)=D(1,1)
  30 :          D(4,4)=E1*(1.D0-2.D0*XNU)/2.D0
  31 :          D(5,5)=D(4,4)
  32 :          D(6,6)=D(5,5)
  33 : C 
  34 :       DO 10 I=1,NSTRSS
  35 :       DSIGE(I)=0.D0
  36 :       DO 10 J=1,NSTRSS
  37 :       DSIGE(I)=DSIGE(I)+D(I,J)*DEPST(J)
  38 :    10 CONTINUE
  39 : C
  40 :       DO 15 I=1,NSTRSS
  41 :        SIGR(I)=SIG0(I)/((1.-ENDT)*(1.-ENDC))
  42 :    15 CONTINUE
  43 : C    
  44 :       DO 30 I=1,NSTRSS
  45 :          SIGF(I)=SIGR(I)+DSIGE(I)
  46 :    30 CONTINUE
  47 : C
  48 : C ***** TEST DES CRITERES ACTIVES ***************
  49 : C
  49 :       CALL DRUTRA(SIGF,SEQTE,BETINSA)
  49 :       CALL DRUCOM(SIGF,SEQCE,BETINSA)
  50 : C
  51 :       FCRIT = SEQTE - SEQT
  52 :       FCRIC = SEQCE - SEQC
  53 : C
  54 :       IF(FCRIT.LT.0.D0.AND.FCRIC.LT.0.D0) IZON=0
  55 :       IF(FCRIT.GE.0.D0.AND.FCRIC.LT.0.D0) IZON=1
  56 :       IF(FCRIT.LT.0.D0.AND.FCRIC.GE.0.D0) IZON=2
  57 :       IF(FCRIT.GE.0.D0.AND.FCRIC.GE.0.D0) IZON=3
  58 : C
  59 :       IF(IZON.EQ.1) THEN      
  59 :         CALL ACTI1(SIGR,SIGF,D,NSTRSS,BETINSA)
  59 :         CALL ENDAME(1,BETINSA)
  60 :              GOTO 100
  61 :       ENDIF
  62 : C
  63 :       IF(IZON.EQ.2) THEN       
  63 :         CALL ACTI2(SIGR,SIGF,D,NSTRSS,BETINSA)
  63 :         CALL ENDAME(2,BETINSA)
  64 :              GOTO 100
  65 :       ENDIF
  66 : C
  67 :       IF(IZON.EQ.3) THEN      
  67 :         CALL ACTI3(SIGR,SIGF,D,NSTRSS,BETINSA)
  67 :         CALL ENDAME(1,BETINSA)
  67 :         CALL ENDAME(2,BETINSA)
  68 :              GOTO 100
  69 :       ENDIF
  70 : C      
  71 :   100 CONTINUE
  72 : C      
  73 :       DO 150 I=1,NSTRSS
  74 :        SIGF(I)=SIGF(I)*((1.-ENDT)*(1.-ENDC))
  75 :   150 CONTINUE
  76 :       
  77 :       RETURN
  78 :       END
  79 : 
  80 :  
  81 :  
  82 :