Télécharger kgrav3.eso

Retour à la liste

Numérotation des lignes :

  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. -INC CCOPTIO
  34. -INC SMRIGID
  35. -INC SMELEME
  36. -INC SMCOORD
  37.  
  38. REAL*8 XE(3,3)
  39.  
  40.  
  41. MELEME=IPMAIL
  42. NBNN=NUM(/1)
  43. NBELEM=NUM(/2)
  44. *
  45. xMATRI=IPMATR
  46. NLIGRP=18
  47. NLIGRD=18
  48.  
  49. C
  50. C BOUCLE DE CALCUL POUR LES DIFFERENTS ELEMENTS
  51. C
  52. DO 3027 IB=1,NBELEM
  53. C
  54. C ON CHERCHE LES COORDONNEES DE L ELEMENT IB
  55. C
  56. CALL DOXE(XCOOR,IDIM,NBNN,NUM,IB,XE)
  57.  
  58. * SEGINI XMATRI
  59. * IMATTT(IB)=XMATRI
  60. IXMATR = XMATRI
  61. CALL DKTGR1(XE,XRG,IPOIN1,IDISS,IXMATR,ib)
  62. * SEGDES XMATRI
  63. 3027 CONTINUE
  64. SEGDES xMATRI
  65. RETURN
  66. END
  67.  
  68.  
  69.  

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