Esope source files

   1 : C TUFIHO    SOURCE    CHAT      05/01/13    03:55:44     5004           
   2 :       SUBROUTINE TUFIHO(THETA,DTUFI,CTUFI,YOUN,RAYOM,EPAI)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)
   5 : C    Include contenant quelques constantes dont XPI :
   6 : -INC CCREEL
   7 :       DIMENSION DTUFI(2,*),CTUFI(2,*)
   8 : C=======================================================================
   9 : C       CE SOUS-PROGRAMME CALCULE LA MATRICE DE HOOKE
  10 : C           POUR L ELEMENT TUFI
  11 : C        IL EST APPELE PAR TUFIEC
  12 : C     
  13 : C         ENTREE : THETA : ANGLE DE LA FISSURE
  14 : C         SORTIE : DTUFI : MATRICE DE HOOKE
  15 : C                  CTUFI : MATRICE DE COMPLAISANCE
  16 : C=======================================================================
  16 :       CALL TUFIZA(A,RAYOM,EPAI)
  17 :       TSP=THETA/XPI
  18 :       TSP1=TSP**1.5D0
  19 :       TSP2=TSP**4.24D0
  20 :       TSP3=TSP**3.D0
  21 :       TSP4=TSP**5.74D0
  22 :       TSP5=TSP**8.48D0
  23 : C     
  24 : C   CALCUL DES COEF DE LA MATRICE DE COMPLAISANCE
  25 : C     
  26 :       FP=A*(6.09D0*TSP1+12.03D0*TSP2)
  27 :       FP=FP+1.D0
  28 :       FP=FP+A*A*(11.36D0*TSP3+51.70D0*TSP4+67.24D0*TSP5)
  29 :       FM=A*(5.26D0*TSP1+1.69D0*TSP2)
  30 :       FM=FM+1.D0
  31 :       FM=FM+A*A*(8.46D0*TSP3+6.28D0*TSP4+1.33D0*TSP5)
  32 :       FMP=A*(5.67D0*TSP1+6.86D0*TSP2)
  33 :       FMP=FMP+1.D0
  34 :       FMP=FMP+A*A*(9.81D0*TSP3+25.95D0*TSP4+9.46D0*TSP5)
  35 : C     
  36 :       CP=THETA*THETA*FP
  37 :       CP=CP/2.D0/YOUN/XPI/EPAI
  38 :       CM=2.D0*THETA*THETA*FM
  39 :       CM=CM/YOUN/XPI/RAYOM/RAYOM/EPAI
  40 :       CMP=THETA*THETA*FMP
  41 :       CMP=CMP/YOUN/XPI/RAYOM/EPAI
  42 : C     
  43 : C   REMPLISSAGE DE LA MATRICE DE COMPLAISANCE
  44 : C     
  45 :       CTUFI(1,1)=CM
  46 :       CTUFI(1,2)=CMP
  47 :       CTUFI(2,1)=CMP
  48 :       CTUFI(2,2)=CP
  49 : C     
  50 : C   INVERSION ET REMPLISSAGE DE LA MATRICE DE HOOKE
  51 : C     
  52 :       DET=CM*CP-CMP*CMP
  53 :       DTUFI(1,1)=CP/DET
  54 :       DTUFI(1,2)=-CMP/DET
  55 :       DTUFI(2,1)=DTUFI(1,2)
  56 :       DTUFI(2,2)=CM/DET
  57 :       RETURN
  58 :       END
  59 : 
  60 : 
  61 :  
  62 :  
  63 :