1 : C CDCY SOURCE CHAT 05/01/12 21:51:18 5004 2 : C CDCY SOURCE INSL 24/10/96 3 : SUBROUTINE CDCY(EPEQ,SEQ,IPL,PENT,EX,TU,EMAX,EPO,ICAL) 4 : C 5 : IMPLICIT INTEGER(I-N) 6 : IMPLICIT REAL*8(A-H,O-Z) 7 : C 8 : C ******************************* 9 : C *** COURBE DE COMPRESSION *** 10 : C *** UNIAXIALE *** 11 : C ******************************* 12 : C------------------------------------------------------------------- 13 : IF(EPEQ .GE.EMAX) THEN 14 : SEQ =0.D0 15 : IPL = 4 16 : PENT=0.D0 17 : GOTO 10 18 : ENDIF 19 : C------------------------------------------------------------------- 20 : IF(ICAL.EQ.1) THEN 21 : C ****************************** 22 : C ** PRE-PIC LINEAIRE ** 23 : C ****************************** 24 : IF(EPEQ.LT.EPO) THEN 25 : SEQ =EX*EPEQ 26 : PENT=EX 27 : IPL = 0 28 : ELSE 29 : SEQ =TU*(EMAX-EPEQ)/(EMAX-EPO) 30 : PENT=-TU/(EMAX-EPO) 31 : IPL =2 32 : ENDIF 33 : GOTO 10 34 : ENDIF 35 : C------------------------------------------------------------------- 36 : C ****************************** 37 : C ** COURBE DE COMPRESSION *** 38 : C ** INSA-LBS1 *** 39 : C ****************************** 40 : IF(EPEQ.LT.EPO) THEN 41 : B3=EX/TU-2.D0/EPO 42 : C3=1.D0/EPO/EPO 43 : DIN3=1.D0+B3*EPEQ+C3*EPEQ*EPEQ 44 : SEQ=EX*EPEQ/DIN3 45 : PENT=EX*(1.D0-C3*EPEQ*EPEQ)/DIN3/DIN3 46 : IPL = 0 47 : ELSE 48 : SEQ=TU*(EMAX-EPEQ)/(EMAX-EPO) 49 : PENT=-TU/(EMAX-EPO) 50 : IPL=2 51 : ENDIF 52 : C------------------------------------------------------------------- 53 : 10 CONTINUE 54 : RETURN 55 : END 56 : 57 : 58 :
© Cast3M 2003 - All rights reserved.
Disclaimer