Esope source files

   1 : C DDRF      SOURCE    CHAT      05/01/12    22:38:25     5004           
   2 : C DDRF     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE DDRF(S2,PHI,V2)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION S2(3),V2(3),TE(3,3)
   8 : C
   9 : C         ************************************
  10 : C         * DETERMINE LE VECTEUR DEFORMATION *
  11 : C         * DANS UN REPERE LIE  A LA FISSURE *
  12 : C         ************************************
  13 : C
  13 :       CALL ZERO(V2,3,1)
  14 :       PI=4.D0*DATAN(1.D0)
  15 :       PHIC=(PHI+90.D0)*PI/180.D0
  16 :       A=DCOS(PHIC)
  17 :       B=DSIN(PHIC)
  18 : *
  19 :       TE(1,1)=A*A
  20 :       TE(1,2)=B*B
  21 :       TE(1,3)=A*B
  22 :       TE(2,1)=B*B
  23 :       TE(2,2)=A*A
  24 :       TE(2,3)=-1.D0*A*B
  25 :       TE(3,1)=-2.D0*A*B
  26 :       TE(3,2)= 2.D0*A*B
  27 :       TE(3,3)= A*A-B*B
  27 :       CALL BST(TE,S2,3,3,V2)
  28 : *
  29 :       RETURN
  30 :       END
  31 :  
  32 :