Esope source files

   1 : C VNSIMP    SOURCE    MAGN      10/05/18    21:16:38     6675           
   2 :       SUBROUTINE VNSIMP(IPT1,IZAFM,IZCCH2,IZIPAD)
   3 : C************************************************************************       
   4 : C   CALCUL DE LA MATRICE M ( P*DIV(U) )  --> AM(NP,IES,NEL)                     
   5 : C                                                                               
   6 : C   IKAS=1  V NORMALE IMPOSEE
   7 : C   IKAS=-1 V TANGENTE IMPOSEE
   8 : C
   9 : C************************************************************************       
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  12 :                                                                                 
  13 : C***                                                                            
  14 :       CHARACTER*8 NOM0
  15 :       DIMENSION XYZ(3,9)
  16 : C***                                                                            
  17 :                                                                                 
  18 : -INC SMELEME
  19 : C-INC SMMATRAKANC
  20 : C*************************************************************************
  21 : C                                                                               
  22 : C    REPERAGE ET STOKAGE DES MATRICES ELEMENTAIRES  puis assemblees
  23 : C                                                                               
  24 : 
  25 : * LGEOC SPG de la pression et/ou des multiplicateurs de Lagrange
  26 : * (points CENTRE ) pour chaque operateur de contrainte
  27 : * KGEOC SPG pour la totalite des points CENTRE.
  28 : * KGEOS SPG pour la totalite des points SOMMET (Diagonale vitesse)
  29 : * KLEMC Connectivites de l'ensemble des contraintes
  30 : * LIZAFM(NBSOUS) contient les pointeurs IZAFM des sous-zones
  31 : 
  32 :       SEGMENT MATRAK
  33 :       INTEGER LGEOC(NBOP),IDEBS(NBOP),IFINS(NBOP)
  34 :       INTEGER LIZAFM(NBSOUS)
  35 :       INTEGER IKAM0 (NBSOUS)
  36 :       INTEGER IMEM  (NBELC)
  37 :       INTEGER KLEMC,KGEOS,KGEOC,KDIAG,KCAC,KIZCL,KIZGC
  38 :       ENDSEGMENT
  39 : 
  40 :       SEGMENT IZAFM                                                             
  41 :       REAL*8  AM(NNELP,NP,IESP),RPGI(NELAX)
  42 :       ENDSEGMENT                                                                
  43 : 
  44 :       POINTEUR IPMJ.IZAFM,IPMK.IZAFM
  45 : 
  46 : C*******************************************************************            
  47 : -INC SMLENTI
  48 :       POINTEUR IZIPAD.MLENTI
  49 : -INC SIZFFB
  50 : -INC SMCOORD
  51 : -INC SMCHPOI
  52 :       POINTEUR IZCCH2.MPOVAL
  53 : -INC CCOPTIO
  54 : -INC CCGEOME
  55 : C***                                                                            
  56 :       IKAS=1
  57 :       SEGACT IPT1,IZAFM*MOD
  58 :       IIES=IDIM
  59 :       NBREF=IPT1.LISREF(/1)-1
  60 : C     write(6,*)' nbref=',nbref
  61 : 
  62 :       DO 1 KS=1,NBREF
  63 :       MELEME=IPT1.LISREF(KS)
  64 : C     write(6,*)' meleme=',meleme
  65 :       SEGACT MELEME
  66 :       NEL=NUM(/2)
  67 :       NP=NUM(/1)
  68 : 
  69 :       NOM0=NOMS(ITYPEL)//'    '
  69 :       CALL KALPBG(NOM0,'FONFORM ',IZFFM)
  70 :       SEGACT IZFFM*MOD
  71 :       IZHR=KZHR(1)
  72 :       SEGACT IZHR*MOD
  73 :       NES=GR(/1)
  74 :       NPG=GR(/3)
  75 : 
  76 :       NELAM=AM(/1)
  77 : 
  77 :       CALL KAM101(NELAM,NEL,NP,IIES,NUM,AM,XCOOR,
  78 :      &IZIPAD.LECT,IKAS,IAXI,NES,NPG,FN,GR,HR,PG,PGSQ,RPG)
  79 :       SEGDES MELEME
  80 :       SEGSUP IZFFM,IZHR
  81 :  1    CONTINUE
  82 :       IESAM=AM(/3)
  83 :       SEGACT IZCCH2*MOD
  83 :       CALL KANOR2(AM,IZCCH2.VPOCHA,NELAM,IESAM)
  84 :       SEGDES IZCCH2
  85 :       SEGDES IPT1,IZAFM
  86 : C                                                                               
  87 : C                                                                               
  88 :       RETURN                                                                    
  89 :  1001 FORMAT(20(1X,I5))                                                         
  90 :  1002 FORMAT(10(1X,1PE11.4))                                                    
  91 :       END                                                                       
  92 :  
  93 :  
  94 :  
  95 :  
  96 :  
  97 :  
  98 :  
  99 :  
 100 :  
 101 :  
 102 :