1 : C LISPDF SOURCE CHAT 05/01/13 01:22:37 5004 2 : SUBROUTINE LISPDF(XN,XM,GA,GB,QSI,EPAI,SIGY,DFIDQS,DFIDM) 3 : C=================================================================== 4 : C CALCULE LA CONTRIBUTION PLASTIQUE DE J 5 : C ENTREES 6 : C XN,XM EFFORTS MEMBRANE FLEXION 7 : C GA GB PARAMETRES POUR CALCULER CRITERE 8 : C QSI RAPPORT PROFONDEUR FISSURE SUR EPAISSEUR 9 : C EPAI EPAISSEUR 10 : C SIGY CONTRAINTE ULTIME 11 : C SORTIES 12 : C DFIDQS = DERIVEE DU CRITERE PAR RAPPORT A QSI 13 : C EBERSOLT MAI 87 14 : C==================================================================== 15 : IMPLICIT INTEGER(I-N) 16 : IMPLICIT REAL*8(A-H,O-Z) 17 : PARAMETER(UNDEMI=.5D0,TRDEMI=1.5D0,DEUX=2.D0,UN=1.D0) 18 : PARAMETER(QUATRE=4.D0) 19 : C 20 : ST=UN/(SIGY*EPAI) 21 : STT= ST/EPAI 22 : C 23 : DGADQS = TRDEMI - DEUX* QSI 24 : DGBDQS = UNDEMI - DEUX* QSI 25 : XNBAR =XN * ST 26 : XMBAR =XM * STT * QUATRE 27 : AAAA = ( XNBAR+XMBAR+GA)/(UN-GA) 28 : BBBB = (-XNBAR+XMBAR+GB)/(UN+GB) 29 : DFIDQS = (DEUX*AAAA + AAAA*AAAA)*DGADQS 30 : 1 +(DEUX*BBBB - BBBB*BBBB)*DGBDQS+UN 31 : DFIDM = (DEUX*AAAA + DEUX*BBBB)*STT*QUATRE 32 : RETURN 33 : END 34 :
© Cast3M 2003 - All rights reserved.
Disclaimer