Esope source files

   1 : C LISPQ     SOURCE    CHAT      05/01/13    01:23:23     5004           
   2 :       SUBROUTINE LISPQ(XNE,XME,W,SIGY,GA,GB,QSI,Q)                              
   3 : C=======================================================================        
   4 : C    CALCULE LE CRITERE POUR LA PLASTICITE DU LINESPRING                        
   5 : C   ENTREE                                                                      
   6 : C    XNE,XME= EFFORT DE MEMBRANE ET DE TRACTION                                 
   7 : C    W      = EPAISSEUR DE LA COQUE OU DU LINESPRING                            
   8 : C    SIGY   = CONTRAINTES ULTIME                                                
   9 : C    GA GB  = CONSTANTES POUR CALCULER LE CRITERE                               
  10 : C    QSI    = RAPPORT PROFONDEUS SUR EPAISSEUR                                  
  11 : C   SORTIES                                                                     
  12 : C    Q      = VALEUR DE LA FONCTION CRITERE                                     
  13 : C         EBERSOLT AVRIL 85                                                     
  14 : C======================================================================         
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  17 :       PARAMETER(QUATRE=4.D0,UN=1.0D0,DEUX=2.D0)                                 
  18 : C                                                                               
  19 :       XNB= XNE / ( SIGY * W )                                                   
  20 :       XMB= XME *QUATRE / ( SIGY * W *W )                                        
  21 : C                                                                               
  22 :       XXA = XNB + XMB + GA                                                      
  23 :       XXB =-XNB + XMB + GB                                                      
  24 : C                                                                               
  25 :       XXA = XXA * XXA / (UN - GA)                                               
  26 :       XXB = XXB * XXB / (UN + GB)                                               
  27 : C                                                                               
  28 :       Q = XXA + XXB - DEUX +QSI                                                 
  29 :       RETURN                                                                    
  30 :       END                                                                       
  31 :