Télécharger vtsimp.eso

Retour à la liste

Numérotation des lignes :

  1. C VTSIMP SOURCE PV 16/11/17 22:01:45 9180
  2. SUBROUTINE VTSIMP(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.  
  54. -INC PPARAM
  55. -INC CCOPTIO
  56. -INC CCGEOME
  57. C***
  58. IKAS=-1
  59. SEGACT IPT1,IZAFM*MOD
  60. IIES=IDIM
  61. NBREF=IPT1.LISREF(/1)-1
  62.  
  63. DO 1 KS=1,NBREF
  64. MELEME=IPT1.LISREF(KS)
  65. SEGACT MELEME
  66. NEL=NUM(/2)
  67. NP=NUM(/1)
  68.  
  69. NOM0=NOMS(ITYPEL)//' '
  70. CALL KALPBG(NOM0,'FONFORM ',IZFFM)
  71. SEGACT IZFFM*MOD
  72. IZHR=KZHR(1)
  73. SEGACT IZHR*MOD
  74. NES=GR(/1)
  75. NPG=GR(/3)
  76.  
  77. NELAM=AM(/1)
  78.  
  79. CALL KAM101(NELAM,NEL,NP,IIES,NUM,AM,XCOOR,
  80. &IZIPAD.LECT,IKAS,IAXI,NES,NPG,FN,GR,HR,PG,PGSQ,RPG)
  81. SEGDES MELEME
  82. SEGSUP IZFFM,IZHR
  83. 1 CONTINUE
  84. IESAM=AM(/3)
  85. SEGACT IZCCH2*MOD
  86. CALL KANOR2(AM,IZCCH2.VPOCHA,NELAM,IESAM)
  87. SEGDES IZCCH2
  88. SEGDES IPT1,IZAFM
  89. C
  90. C
  91. RETURN
  92. 1001 FORMAT(20(1X,I5))
  93. 1002 FORMAT(10(1X,1PE11.4))
  94. END
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales