Télécharger rigli2.eso

Retour à la liste

Numérotation des lignes :

  1. C RIGLI2 SOURCE CHAT 05/01/13 03:00:02 5004
  2. C RIGLI2 SOURCE DJER 94/12/19 21:15:07 654
  3. SUBROUTINE RIGLI2(REL,LRE,BPSS,WORK)
  4. C====================================================================
  5. C
  6. C MODIFIE REL (LRE,LRE) MATRICE DE RIGIDITE EN INTRODUISANT DES
  7. C RAIDEUR LOCALES EQUIVALENTES : Klux Kluy Kluz Klrx Klry Klrz
  8. C====================================================================
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8(A-H,O-Z)
  11. DIMENSION REL(LRE,LRE),BPSS(3,3),WORK(9),C(6),RE(6,6),C1(3)
  12. C
  13. CALL ZERO(REL,LRE,LRE)
  14. CALL ZERO(RE,6,6)
  15. XZER = 0.D0
  16. C
  17. DO I =1,6
  18. RE(I,I)=WORK(I)
  19. END DO
  20. C
  21. DO 5 I =1,3
  22. DO 6 I1=1,3
  23. C1(I1)=BPSS(I1,I)
  24. 6 CONTINUE
  25. DO 1 J=1,3
  26. S1=XZER
  27. S2=XZER
  28. DO 2 K=1,3
  29. S1=S1+C1(K)*RE(K,J)
  30. S2=S2+C1(K)*RE(K+3,J+3)
  31. 2 CONTINUE
  32. C(J) =S1
  33. C(J+3)=S2
  34. 1 CONTINUE
  35. DO 3 J=1,3
  36. S1=XZER
  37. S2=XZER
  38. DO 4 K=1,3
  39. S1=S1+C(K)*BPSS(K,J)
  40. S2=S2+C(K+3)*BPSS(K,J)
  41. 4 CONTINUE
  42. REL(I,J) =S1
  43. REL(I+3,J+3)=S2
  44. REL(I+6,J+6)=S1
  45. REL(I+9,J+9)=S2
  46. 3 CONTINUE
  47. 5 CONTINUE
  48.  
  49. *-------------------------------------------------------------
  50. * WRITE(*,*) '** REL APRES ( dans RIGLI2) ** '
  51. * DO 112 IA=1,LRE
  52. * WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
  53. * 112 CONTINUE
  54. *1991 FORMAT(18(1X,E12.5))
  55. *-------------------------------------------------------------
  56. RETURN
  57. END
  58.  
  59.  
  60.  

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