Esope source files

   1 : C LIRANG    SOURCE    CHAT      05/01/13    01:20:21     5004           
   2 :       SUBROUTINE LIRANG(XTOT,XABSCI,XORDON,I,NLIAB,NRG,XPENTE,NIP)             
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : *--------------------------------------------------------------------*          
   6 : *                                                                    *          
   7 : *     Operateur DYNE : algorithme de Fu - de Vogelaere               *          
   8 : *     ________________________________________________               *          
   9 : *                                                                    *          
  10 : *     Repere l'indice NRG du tableau XABSCI representant la liste de *
  11 : *     reels des abscisses de la loi de comportement elasto-plastique.*
  12 : *     Pour une valeur XTOT de l'abscisse, cet indice sera tel que :  *
  13 : *              XABSCI(I,NRG)<= XTOT <= XABSCI(I,NRG+1)               *
  14 : *                                                                    *          
  15 : *     Parametres:                                                    *          
  16 : *                                                                    *          
  17 : * e   XTOT    Valeur donnee de l'abscisse                            *
  18 : * e   XABSCI  Tableau contenant les abscisses de la loi plastique    *
  19 : *             pour la liaison point-point-plastique                  *
  20 : * e   I       numero de la liaison.                                  *
  21 : *             pour la liaison point-point-plastique                  *
  22 : * s   NRG     Rang du point repere                                   *          
  23 : *                                                                    *                                                                          *          
  24 : *     Auteur, date de creation:                                      *          
  25 : *                                                                    *          
  26 : *     Nicolas WECXSTEEN modifications 04/96 point-point-..-plastique *
  27 : *                                                                    *
  28 : *--------------------------------------------------------------------*
  29 : *
  30 :       REAL*8 XABSCI(NLIAB,*),XORDON(NLIAB,*)
  31 : *
  32 : 
  33 : 
  34 :       NRG = 1
  35 :       XVAL = XABSCI(I,2)
  36 :  10   IF (XTOT.GE.XVAL) THEN
  37 :         NRG = NRG + 1
  38 :         if (nrg.eq.NIP) then 
  39 :           XPENTE = (XORDON(I,NRG)-XORDON(I,NRG-1))/
  40 :      &               (XABSCI(I,NRG)-XABSCI(I,NRG-1))
  41 :           go to 20
  42 :         endif
  43 : 
  44 :         XVAL = XABSCI(I,NRG+1)
  45 :         GOTO 10
  46 :       endif
  47 : 
  48 :         XPENTE = (XORDON(I,NRG+1)-XORDON(I,NRG))/
  49 :      &               (XABSCI(I,NRG+1)-XABSCI(I,NRG))
  50 : 
  51 :  20    continue
  52 : 
  53 :       END
  54 :        
  55 :  
  56 :