Esope source files

   1 : C PAMAOT    SOURCE    CHAT      05/01/13    02:08:54     5004           
   2 : C PAMAOT     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE PAMAOT(ALPHA,A1,B1,DK1,DK2)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C             **********************************************
   8 : C             *     CALCUL DES  PARAMETRES DU MATERIAU     *
   9 : C             *           (  MODEL D'OTTOSEN   )           *
  10 : C             **********************************************
  11 : C-------------------------------------------------------------------
  12 :       DIMENSION D(3,3),D1(3,3),AA(3),BB(3),XK1(3),XK2(3),CC(3),AL(3)
  13 :       DIMENSION ALI(2)
  14 : C-------------------------------------------------------------------
  15 :       AA(1) =1.8076D0
  16 :       BB(1) =4.0962D0
  17 :       XK1(1)=14.4863D0
  18 :       XK2(1)=0.9914D0
  19 :       AL(1) =0.08D0
  20 : *
  21 :       AA(2) =1.2759D0
  22 :       BB(2) =3.1962D0
  23 :       XK1(2)=11.7365D0
  24 :       XK2(2)=0.9801D0
  25 :       AL(2) =0.1D0
  26 : *
  27 :       AA(3) =0.9218D0
  28 :       BB(3) =2.5969D0
  29 :       XK1(3)=9.911D0
  30 :       XK2(3)=0.9647D0
  31 :       AL(3) =0.12D0
  32 : C-------------------------------------------------------------------
  33 :       ALI(1)=(AL(1)+AL(2))/2.D0
  34 :       ALI(2)=(AL(2)+AL(3))/2.D0
  35 : C-------------------------------------------------------------------
  36 :       IF(ALPHA.LE.ALI(1)) THEN
  37 :         A1 =AA(1)
  38 :         B1 =BB(1)
  39 :         DK1=XK1(1)
  40 :         DK2=XK2(1)
  41 :         GOTO 200
  42 :        ENDIF
  43 :       IF(ALPHA.LE.ALI(2).AND.ALPHA.GT.ALI(1)) THEN
  44 :         A1 =AA(2)
  45 :         B1 =BB(2)
  46 :         DK1=XK1(2)
  47 :         DK2=XK2(2)
  48 :         GOTO 200
  49 :        ENDIF
  50 :       IF(ALPHA.GT.ALI(2)) THEN
  51 :         A1 =AA(3)
  52 :         B1 =BB(3)
  53 :         DK1=XK1(3)
  54 :         DK2=XK2(3)
  55 :         GOTO 200
  56 :        ENDIF
  57 : C-------------------------------------------------------------------
  58 :       D(1,1)=AL(1)*AL(1)
  59 :       D(1,2)=AL(1)
  60 :       D(1,3)=1.D0
  61 :       D(2,1)=AL(2)*AL(2)
  62 :       D(2,2)=AL(2)
  63 :       D(2,3)=1.D0
  64 :       D(3,1)=AL(3)*AL(3)
  65 :       D(3,2)=AL(3)
  66 :       D(3,3)=1.D0
  67 : *
  67 :       CALL INVER3(D,D1)
  67 :       CALL BST(D1,AA,3,3,CC)
  68 :       A1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  68 :       CALL ZERO(CC,3,1)
  68 :       CALL BST(D1,BB,3,3,CC)
  69 :       B1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  69 :       CALL ZERO(CC,3,1)
  69 :       CALL BST(D1,XK1,3,3,CC)
  70 :       DK1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  70 :       CALL ZERO(CC,3,1)
  70 :       CALL BST(D1,XK2,3,3,CC)
  71 :       DK2=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  72 : C-------------------------------------------------------------------
  73 :   200 CONTINUE
  74 : C-------------------------------------------------------------------
  75 :       RETURN
  76 :       END
  77 :  
  78 :