Esope source files

   1 : C D2VENE    SOURCE    CHAT      05/01/12    22:35:15     5004           
   2 : C DEVENE    SOURCE    IANIS     96/03/15    21:15:48     2089           
   3 :       SUBROUTINE D2VENE (NA1,IND,NPAS,FEXA,Q1,Q2,FTOTA,WEXT,
   4 :      &              WINT,XASM,NPC1)                 
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)        
   7 : 
   8 : **************************************************************
   9 : *
  10 : *  Ce sous-programme calcule le travail des differentes forces
  11 : *  pour chaque mode
  12 : *
  13 : * e   NA1    nombre d'inconnues sur base A
  14 : * e   IND    idice du pas (2 -->1er 1/2 pas, 1 --> 2eme 1/2pas)
  15 : * e   FEXA   forces exterieures sur base A
  16 : * e   NPAS   numero du pas 
  17 : * e   Q1     deplacements generalises
  18 : * e   Q2     vitesses generalisees
  19 : * e   FTOTA  forces totales sur base A
  20 : * e,s WEXT   travail des forces exterieures 
  21 : * e,s WINT   travail des forces interieures (rigidite et
  22 : *            amortissement et forces de liaison )
  23 : * e   XASM   tableau d'amortissements generalises
  24 : * e   NPC1   Nombre de pas de calcul  
  25 : 
  26 : 
  27 : **
  28 : ** ce sous-programme n'est pas, peut - etre, compatible avec 
  29 : ** les liaisons '..._FLUIDE' 
  30 : *
  31 : * I. Politopoulos Mars 1996
  32 : ********************************************************************
  33 : 
  34 :       REAL*8 FEXA(NA1,NPC1,*),FTOTA(NA1,*),Q1(NA1,4),Q2(NA1,4)
  35 :       REAL*8 WEXT(NA1,2),WINT(NA1,2),XASM(NA1)
  36 : 
  37 :         DO 20 I=1, NA1
  38 :           DQ1 = Q1(I,1) - Q1(I,2)
  39 :    
  40 :           F1 = FEXA(I,NPAS,2)
  41 :           F2 = FEXA(I,NPAS,1) 
  42 :           DWEX = (F1 + F2)/2D0 * DQ1
  43 : 
  44 :           F1 = FTOTA(I,2) - F1
  45 :           F2 = FTOTA(I,1) - F2
  46 : 
  47 : *  on ajoute les forces d'amortissement
  48 :           F1 = F1 - XASM(I)*Q2(I,2)
  49 :           F2 = F2 - XASM(I)*Q2(I,1)
  50 :           DWIN = -(F1 + F2)/2D0 * DQ1
  51 : 
  52 :           WEXT(I,1) = WEXT(I,2) + DWEX
  53 :           WINT(I,1) = WINT(I,2) + DWIN
  54 : 
  55 :   20     CONTINUE
  56 :  
  57 :       
  58 :       END
  59 : 
  60 :  
  61 :  
  62 :