Télécharger la plaquette Cast3M


fr en

Télécharger pamaot.eso

Retour à la liste

Numérotation des lignes : 

pamaot

  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. *
  68.       CALL INVER3(D,D1)
  69.       CALL BST(D1,AA,3,3,CC)
  70.       A1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  71.       CALL ZERO(CC,3,1)
  72.       CALL BST(D1,BB,3,3,CC)
  73.       B1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  74.       CALL ZERO(CC,3,1)
  75.       CALL BST(D1,XK1,3,3,CC)
  76.       DK1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  77.       CALL ZERO(CC,3,1)
  78.       CALL BST(D1,XK2,3,3,CC)
  79.       DK2=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)
  80. C-------------------------------------------------------------------
  81.   200 CONTINUE
  82. C-------------------------------------------------------------------
  83.       RETURN
  84.       END
  85.  
  86.  
  87.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales