C PAMAOT SOURCE CHAT 05/01/13 02:08:54 5004 C 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.12D0 C------------------------------------------------------------------- ALI(1)=(AL(1)+AL(2))/2.D0 ALI(2)=(AL(2)+AL(3))/2.D0 C------------------------------------------------------------------- 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 ENDIF C------------------------------------------------------------------- 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 CONTINUE C------------------------------------------------------------------- RETURN END