Télécharger kgrav2.eso

Retour à la liste

Numérotation des lignes :

  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. -INC CCOPTIO
  35. -INC SMRIGID
  36. -INC SMELEME
  37. -INC SMCOORD
  38.  
  39.  
  40. REAL*8 XE(3,4)
  41.  
  42.  
  43. MELEME=IPMAIL
  44. NBNN=NUM(/1)
  45. NBELEM=NUM(/2)
  46. *
  47. xMATRI=IPMATR
  48. NLIGRP=24
  49. NLIGRD=24
  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 CQ4GR1(XE,XRG,IPOIN1,IDISS,xmatri,IPMINT,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