Télécharger la plaquette Cast3M


fr en

Télécharger kgrav3.eso

Retour à la liste

Numérotation des lignes : 

kgrav3

  1. C KGRAV3    SOURCE    CHAT      09/10/09    21:19:24     6519
  2.  
  3.       SUBROUTINE KGRAV3(IPMAIL,XRG,IPOIN1,IPMATR,IDISS)
  4. *---------------------------------------------------------------------*
  5. *      APPELEE PAR KGRAV1
  6. *                                                                     *
  7. * calcule la matrice de rigidité liée à la variation de pression à cause
  8. * d'un mouvement suivant le vecteur IPOIN1.
  9. *
  10. *  éléments COQ3,DKT,DST
  11. *                                                                     *
  12. *---------------------------------------------------------------------*
  13. *                                                                     *
  14. *   ENTREES :                                                         *
  15. *   ________                                                          *
  16. *                                                                     *
  17. *        IPMAIL   Pointeur sur un segment  MELEME                     *
  18. *        XRG      coefficient multiplicateur de la rigidité           *
  19. *        IPOIN1   vecteur(type point), definissant le sens de gravité *
  20. *        IPMATR   pointeur sur la matrice de rigidité de pesanteur
  21. *                 de la sous-zone
  22. *        IDISS : 0 ---> MATRICE SYMETRIQUE
  23. *                1 ---> MATRICE DISSYMETRIQUE
  24. *                                                                     *
  25. *   SORTIES :                                                         *
  26. *   ________
  27. *        on remplit les segments xmatri et imatri pour la sous zone
  28. *
  29. *                                                                     *
  30. *---------------------------------------------------------------------*
  31.       IMPLICIT INTEGER(I-N)
  32.       IMPLICIT REAL*8(A-H,O-Z)
  33.  
  34. -INC PPARAM
  35. -INC CCOPTIO
  36. -INC SMRIGID
  37. -INC SMELEME
  38. -INC SMCOORD
  39.  
  40.       REAL*8 XE(3,3)
  41.  
  42.  
  43.       MELEME=IPMAIL
  44.       NBNN=NUM(/1)
  45.       NBELEM=NUM(/2)
  46. *
  47.       xMATRI=IPMATR
  48.       NLIGRP=18
  49.       NLIGRD=18
  50.  
  51. C
  52. C     BOUCLE DE CALCUL POUR LES DIFFERENTS ELEMENTS
  53. C
  54.       DO 3027 IB=1,NBELEM
  55. C
  56. C     ON CHERCHE LES COORDONNEES DE L ELEMENT IB
  57. C
  58.          CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  59.  
  60. *      SEGINI XMATRI
  61. *      IMATTT(IB)=XMATRI
  62.       IXMATR = XMATRI
  63.       CALL DKTGR1(XE,XRG,IPOIN1,IDISS,IXMATR,ib)
  64. *      SEGDES XMATRI
  65.  3027 CONTINUE
  66.       SEGDES xMATRI
  67.       RETURN
  68.       END
  69.  
  70.  
  71.  

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