Numérotation des lignes :

pamaot
C PAMAOT    SOURCE    CHAT      05/01/13    02:08:54     5004C PAMAOT     SOURCE    INSL       24/10/96      SUBROUTINE PAMAOT(ALPHA,A1,B1,DK1,DK2)C      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)C             **********************************************C             *     CALCUL DES  PARAMETRES DU MATERIAU     *C             *           (  MODEL D'OTTOSEN   )           *C             **********************************************C-------------------------------------------------------------------      DIMENSION D(3,3),D1(3,3),AA(3),BB(3),XK1(3),XK2(3),CC(3),AL(3)      DIMENSION ALI(2)C-------------------------------------------------------------------      AA(1) =1.8076D0      BB(1) =4.0962D0      XK1(1)=14.4863D0      XK2(1)=0.9914D0      AL(1) =0.08D0*      AA(2) =1.2759D0      BB(2) =3.1962D0      XK1(2)=11.7365D0      XK2(2)=0.9801D0      AL(2) =0.1D0*      AA(3) =0.9218D0      BB(3) =2.5969D0      XK1(3)=9.911D0      XK2(3)=0.9647D0      AL(3) =0.12D0C-------------------------------------------------------------------      ALI(1)=(AL(1)+AL(2))/2.D0      ALI(2)=(AL(2)+AL(3))/2.D0C-------------------------------------------------------------------      IF(ALPHA.LE.ALI(1)) THEN        A1 =AA(1)        B1 =BB(1)        DK1=XK1(1)        DK2=XK2(1)        GOTO 200       ENDIF      IF(ALPHA.LE.ALI(2).AND.ALPHA.GT.ALI(1)) THEN        A1 =AA(2)        B1 =BB(2)        DK1=XK1(2)        DK2=XK2(2)        GOTO 200       ENDIF      IF(ALPHA.GT.ALI(2)) THEN        A1 =AA(3)        B1 =BB(3)        DK1=XK1(3)        DK2=XK2(3)        GOTO 200       ENDIFC-------------------------------------------------------------------      D(1,1)=AL(1)*AL(1)      D(1,2)=AL(1)      D(1,3)=1.D0      D(2,1)=AL(2)*AL(2)      D(2,2)=AL(2)      D(2,3)=1.D0      D(3,1)=AL(3)*AL(3)      D(3,2)=AL(3)      D(3,3)=1.D0*      CALL INVER3(D,D1)      CALL BST(D1,AA,3,3,CC)      A1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)      CALL ZERO(CC,3,1)      CALL BST(D1,BB,3,3,CC)      B1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)      CALL ZERO(CC,3,1)      CALL BST(D1,XK1,3,3,CC)      DK1=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)      CALL ZERO(CC,3,1)      CALL BST(D1,XK2,3,3,CC)      DK2=CC(1)*ALPHA*ALPHA+CC(2)*ALPHA+CC(3)C-------------------------------------------------------------------  200 CONTINUEC-------------------------------------------------------------------      RETURN      END   

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