endame
C ENDAME SOURCE CHAT 05/01/12 23:38:48 5004 C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) SEGMENT BETINSA REAL*8 RT,RC,YOUN,XNU,GFT,GFC,HC REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC INTEGER IFIS,IPLA,IBB,IGAU ENDSEGMENT C *--------------------------------------------------------------------* * * * CALCUL DE L'ENDOMMAGEMENT * * _________________________ * * * * ENTREES: * * * * DKT PARAMETRE D'ECROUISSAGE EN TRACTION * * DKC PARAMETRE D'ECROUISSAGE EN COMPRESSION * * IC =1 : COURBE de TRACTION * * =2 : COURBE de COMPRESSION * * * * SORTIES: * * * * ENDC ENDOMMAGEMENT DE COMPRESSION * * ENDT ENDOMMAGEMENT DE TRACTION * * * *--------------------------------------------------------------------* C C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C Initialisation C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C BB=3.D0 DT = 0.5D0 DC = 0.25D0 ENDL = 0.99D0 C C C C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C Direction de traction C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C IF (IC.EQ.1) THEN FT=RT AT= -0.5D0 BT=FT*(HC/GFT)*(1.D0+AT/2.D0) CT=BT*LOG(1.D0-DT)/(LOG(((1+AT)-SQRT(1+AT*AT))/(2.D0*AT))) C ENDT=1- EXP(-CT*DKT) C IF (ENDT.GE.ENDL) THEN IFIS=2 ENDIF IF (ENDT.GT.0.D0) THEN IFIS=1 ENDIF ENDIF C C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C Direction de compression C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C IF (IC.EQ.2) THEN FC=RC AC=(2.D0*BB-1.D0)+2.D0* SQRT(BB**2-BB) BC=FC*(HC/GFC)*(1.D0+AC/2.D0) CC=BC*LOG(1.D0-DC)/(LOG((1+AC)/(2.D0*AC))) C ENDC=1- EXP(-CC*DKC) C IF (ENDC.GE.ENDL) THEN IPLA=3 ENDIF IF (ENDC.GE.DC) THEN IPLA=2 ENDIF IF (ENDC.GT.0.D0) THEN IPLA=1 ENDIF ENDIF C C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C CAS OU ENDO=ENDO_LIMITE C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C IF(ENDT.GE.ENDL) ENDT= ENDL IF(ENDC.GE.ENDL) ENDC= ENDL C IF ((ENDT.LT.0.D0).OR.(ENDC.LT.0.D0)) THEN WRITE(*,*)'PROBLEME DANS DAMAGE' WRITE(*,*)'ENDT=',ENDT,'ENDC=',ENDC C STOP ENDIF C C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ C RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales