Esope source files

   1 : C EPSIG0    SOURCE    CHAT      05/01/12    23:42:30     5004           
   2 :       SUBROUTINE EPSIG0(S,EPSI,MFR,IFOUR,YUNG,XNU,CARAC,NSTRS)
   3 : 
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 :       DIMENSION S(*),EPSI(*),CARAC(*)
   7 : 
   8 :       DO I=1,NSTRS
   9 :         EPSI(I)=0.D0
  10 :       ENDDO
  11 : 
  12 :       AUX1=1.D0/YUNG
  13 :       AUX2=2.D0*(1.D0+XNU)
  14 :       AUX3=AUX1*AUX2
  15 :       IF(MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  16 :         EP1=CARAC(1)
  17 :         IF(MFR.NE.5) ALFAH=CARAC(2)**2
  18 :         AUX4=AUX1/EP1
  19 :         AUX5=12.D0*AUX4/(EP1*EP1)
  20 :       ENDIF
  21 : 
  22 : 
  23 : C---------ELEMENTS MASSIFS MFR=1----------------------------------------
  24 : 
  25 :       IF(MFR.EQ.1) THEN
  26 : 
  27 :         EPSI(1)=(S(1)-XNU*(S(2)+S(3)))*AUX1
  28 :         EPSI(2)=(S(2)-XNU*(S(1)+S(3)))*AUX1
  29 :         EPSI(3)=(S(3)-XNU*(S(1)+S(2)))*AUX1
  30 :         DO I=4,NSTRS
  31 :           EPSI(I)=AUX3*S(I)
  32 :         ENDDO
  33 : 
  34 : C---------ELEMENTS COQUES SANS ET AVEC C.T. MFR=3 ET MFR=9--------------
  35 : 
  36 :       ELSE IF(MFR.EQ.3.OR.MFR.EQ.9) THEN
  37 : 
  38 :         IF(IFOUR.GE.1) THEN
  39 :           EPSI(1)=(S(1)-XNU*S(2))*AUX4
  40 :           EPSI(2)=(S(2)-XNU*S(1))*AUX4
  41 :           EPSI(3)=AUX2*AUX4*S(3)
  42 :           EPSI(4)=(S(4)-XNU*S(5))*AUX5
  43 :           EPSI(5)=(S(5)-XNU*S(4))*AUX5
  44 :           EPSI(6)=AUX2*AUX5*S(6)
  45 :         ELSE
  46 :           EPSI(1)=(S(1)-XNU*S(2))*AUX4        
  47 :           EPSI(2)=(S(2)-XNU*S(1))*AUX4
  48 :           EPSI(3)=(S(3)-XNU*S(4))*AUX5
  49 :           EPSI(4)=(S(4)-XNU*S(3))*AUX5    
  50 :         ENDIF
  51 :     
  52 : C---------ELEMENTS COQUES EPAISSES MFR=5--------------------------------
  53 : 
  54 :       ELSE IF(MFR.EQ.5) THEN
  55 : 
  56 :         EPSI(1)=(S(1)-XNU*S(2))*AUX1
  57 :         EPSI(2)=(S(2)-XNU*S(1))*AUX1 
  58 :         EPSI(3)=AUX3*S(3)
  59 :         EPSI(4)=AUX3*S(4)
  60 :         EPSI(5)=AUX3*S(5)
  61 : 
  62 : C---------ELEMENTS POUTRES MFR=7----------------------------------------
  63 : 
  64 :       ELSE IF(MFR.EQ.7) THEN
  65 : 
  66 :         IF(CARAC(4).NE.0.D0) EPSI(1)=AUX1*S(1)/CARAC(4)
  67 :         IF(CARAC(5).NE.0.D0) EPSI(2)=AUX3*S(2)/CARAC(5)
  68 :         IF(CARAC(6).NE.0.D0) EPSI(3)=AUX3*S(3)/CARAC(6)
  69 :         IF(CARAC(1).NE.0.D0) EPSI(4)=AUX3*S(4)/CARAC(1)
  70 :         IF(CARAC(2).NE.0.D0) EPSI(5)=AUX1*S(5)/CARAC(2)
  71 :         IF(CARAC(3).NE.0.D0) EPSI(6)=AUX1*S(6)/CARAC(3)
  72 : 
  73 :       ENDIF
  74 :       END
  75 :  
  76 :  
  77 :