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.  
  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