Sources Esope | Cast3M

lispdf
C LISPDF    SOURCE    CHAT      05/01/13    01:22:37     5004
      SUBROUTINE LISPDF(XN,XM,GA,GB,QSI,EPAI,SIGY,DFIDQS,DFIDM)
C===================================================================
C      CALCULE LA CONTRIBUTION PLASTIQUE DE J
C   ENTREES
C    XN,XM EFFORTS MEMBRANE FLEXION
C      GA GB  PARAMETRES POUR CALCULER CRITERE
C      QSI    RAPPORT PROFONDEUR FISSURE SUR EPAISSEUR
C      EPAI EPAISSEUR
C      SIGY CONTRAINTE ULTIME
C   SORTIES
C     DFIDQS = DERIVEE DU CRITERE PAR RAPPORT A QSI
C           EBERSOLT MAI  87
C====================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      PARAMETER(UNDEMI=.5D0,TRDEMI=1.5D0,DEUX=2.D0,UN=1.D0)
      PARAMETER(QUATRE=4.D0)
C
      ST=UN/(SIGY*EPAI)
      STT= ST/EPAI
C
      DGADQS = TRDEMI - DEUX* QSI
      DGBDQS = UNDEMI - DEUX* QSI
      XNBAR   =XN * ST
      XMBAR   =XM * STT * QUATRE
      AAAA   = ( XNBAR+XMBAR+GA)/(UN-GA)
      BBBB   = (-XNBAR+XMBAR+GB)/(UN+GB)
      DFIDQS = (DEUX*AAAA + AAAA*AAAA)*DGADQS
     1        +(DEUX*BBBB - BBBB*BBBB)*DGBDQS+UN
      DFIDM  = (DEUX*AAAA + DEUX*BBBB)*STT*QUATRE
      RETURN
      END