Esope source files

   1 : C LISPML    SOURCE    CHAT      05/01/13    01:23:01     5004           
   2 :       SUBROUTINE LISPML(QSI,A)                                                  
   3 : C=======================================================================        
   4 : C          PLASTICITE LINESPRING                                                
   5 : C     CALCULE LA FONCTION   A(QSI) POUR LE MOMENT LIMITE                        
   6 : C   ENTREE                                                                      
   7 : C    QSI  = RAPPORT PROFONDEUR DE FISSURE SUR EPAISSEUR                         
   8 : C   SORTIE                                                                      
   9 : C     A   = FONCTION DE QSI     M=(T**2)*SIGF*A(QSI)/4.                         
  10 : C=======================================================================        
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  13 :       PARAMETER(UNDEMI= .5D0,TRDEMI=1.5D0,DEUX=2.D0,UN=1.D0)                    
  14 :       PARAMETER(QUATRE=4.D0)                                                    
  15 : C                                                                               
  16 :       GA= TRDEMI*QSI - QSI*QSI                                                  
  17 :       GB= UNDEMI*QSI - QSI*QSI                                                  
  18 : C                                                                               
  19 :       EA=(UN+GB) + (UN-GA)                                                      
  20 :       EB=DEUX*GA*(UN+GB) + DEUX*GB*(UN-GA)                                      
  21 :       EC=GA*GA*(UN+GB) + GB*GB*(UN-GA) -                                        
  22 :      S   (DEUX-QSI)*(UN-GA)*(UN+GB)                                             
  23 : C                                                                               
  24 :       DELTA=DSQRT(EB*EB - QUATRE*EA*EC)                                         
  25 :       A = (DELTA -EB) /DEUX/EA                                                  
  26 : C                                                                               
  27 :       RETURN                                                                    
  28 :       END                                                                       
  29 :