Télécharger coq8ri.eso

Retour à la liste

Numérotation des lignes :

  1. C COQ8RI SOURCE CHAT 05/01/12 22:22:13 5004
  2. SUBROUTINE COQ8RI(B,FACT,D,LRE,NBPGAU,IA,NBNN,RIG)
  3. C=======================================================================
  4. C ENTREES
  5. C B(5,LRE) : MATRICE DE DERIVEES DE FONCTIONS DE FORMES
  6. C FACT : DETERMINANT DU JACOBIEN * POIDS DE GAUSS
  7. C D(5,5) : MATRICE DE HOOKE
  8. C LRE : NOMBRE DE DDL
  9. C NBPGAU : NOMBRE DE POINTS DE GAUSS
  10. C IA : NUMERO DE POINT DE GAUSS
  11. C NBNN : NOMBRE DE NOEUDS
  12. C
  13. C SORTIES
  14. C RIG(LRE,LRE): MATRICE DE RIGIDITE ELEMENTAIRE
  15. C
  16. C RECUPERATION BILBO OCTOBRE 85
  17. C=======================================================================
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8 (A-H,O-Z)
  20. PARAMETER (XZER=0.D0)
  21. DIMENSION RIG(LRE,*),D(5,*),B(5,*)
  22. DIMENSION DAM(48,5)
  23. C
  24. C CALCUL DE BTDB
  25. C
  26. DO 60 I =1,LRE
  27. DO 60 J = 1,5
  28. DAM(I,J) = XZER
  29. DO 60 K = 1,5
  30. DAM(I,J) = DAM(I,J) + B(K,I)*D(K,J)
  31. 60 CONTINUE
  32. DO 70 I = 1,LRE
  33. DO 70 J = I,LRE
  34. DO 70 K = 1,5
  35. RIG(I,J) = RIG(I,J) + DAM(I,K)*B(K,J)*FACT
  36. 70 CONTINUE
  37. C
  38. IF (IA.EQ.NBPGAU)THEN
  39. DO 110 I = 1,LRE
  40. DO 110 J = I,LRE
  41. 110 RIG(J,I) = RIG(I,J)
  42. C
  43. C ON AJOUTE LES COEFF POUR LA TORSION NORMALE AU PLAN DE LA COQUE
  44. C
  45. VMAX = XZER
  46. DO 120 IK = 1,NBNN
  47. JK4 = (IK-1)*6 + 4
  48. JK6 = JK4 + 2
  49. DO 120 JK = JK4,JK6
  50. VMIN = RIG(JK,JK)
  51. IF (VMAX.GE.VMIN) GO TO 120
  52. VMAX = RIG(JK,JK)
  53. 120 CONTINUE
  54. DO 130 IK = 1,NBNN
  55. JK4 = (IK-1)*6 + 4
  56. JK6 = JK4 + 2
  57. DO 130 JK = JK4,JK6
  58. 130 RIG(JK,JK)=RIG(JK,JK)+1.D-9*VMAX
  59. ENDIF
  60. C
  61. RETURN
  62. END
  63.  
  64.  

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