Télécharger la plaquette Cast3M


fr en

Télécharger kgrav2.eso

Retour à la liste

Numérotation des lignes : 

kgrav2

  1. C KGRAV2    SOURCE    CHAT      09/10/09    21:19:23     6519
  2.  
  3.       SUBROUTINE KGRAV2(IPMAIL,XRG,IPOIN1,IPMATR,IDISS,IPMINT)
  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 COQ4
  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. *        ipmint:  pointeur sur le segment d'integration
  25. *                                                                     *
  26. *   SORTIES :                                                         *
  27. *   ________
  28. *        on remplit les segments xmatri et imatri pour la sous zone
  29. *
  30. *                                                                     *
  31. *---------------------------------------------------------------------*
  32.       IMPLICIT INTEGER(I-N)
  33.       IMPLICIT REAL*8(A-H,O-Z)
  34.  
  35. -INC PPARAM
  36. -INC CCOPTIO
  37. -INC SMRIGID
  38. -INC SMELEME
  39. -INC SMCOORD
  40.  
  41.  
  42.       REAL*8 XE(3,4)
  43.  
  44.  
  45.       MELEME=IPMAIL
  46.       NBNN=NUM(/1)
  47.       NBELEM=NUM(/2)
  48. *
  49.       xMATRI=IPMATR
  50.       NLIGRP=24
  51.       NLIGRD=24
  52.  
  53. C
  54. C     BOUCLE DE CALCUL POUR LES DIFFERENTS ELEMENTS
  55. C
  56.       DO 3027 IB=1,NBELEM
  57. C
  58. C     ON CHERCHE LES COORDONNEES DE L ELEMENT IB
  59. C
  60.          CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  61.  
  62. *      SEGINI XMATRI
  63. *      IMATTT(IB)=XMATRI
  64. *      IXMATR = XMATRI
  65.       CALL CQ4GR1(XE,XRG,IPOIN1,IDISS,xmatri,IPMINT,ib)
  66. *      SEGDES XMATRI
  67.  3027 CONTINUE
  68.       SEGDES xMATRI
  69.       RETURN
  70.       END
  71.  
  72.  
  73.  

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