Télécharger mapalu.eso

Retour à la liste

Numérotation des lignes :

  1. C MAPALU SOURCE CB215821 16/04/21 21:17:43 8920
  2. SUBROUTINE MAPALU(NMATT,WORK,B,IDIM)
  3. C=======================================================================
  4. C RECUPERATION DU REPERE LOCAL DE JOI1 ( element de liaison a 2 noeuds)
  5. C=======================================================================
  6. IMPLICIT INTEGER (I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. DIMENSION B(3,3),WORK(NMATT)
  9. * ---------------------------------------------------
  10. C
  11. IF(IDIM.EQ.3) THEN
  12. VX1=WORK(1)
  13. VY1=WORK(2)
  14. VZ1=WORK(3)
  15. VL1=SQRT(VX1*VX1+VY1*VY1+VZ1*VZ1)
  16. C
  17. VX2=WORK(4)
  18. VY2=WORK(5)
  19. VZ2=WORK(6)
  20. VL2=SQRT(VX2*VX2+VY2*VY2+VZ2*VZ2)
  21. *
  22. B(1,1)=VX1/VL1
  23. B(1,2)=VY1/VL1
  24. B(1,3)=VZ1/VL1
  25. B(2,1)=VX2/VL2
  26. B(2,2)=VY2/VL2
  27. B(2,3)=VZ2/VL2
  28. B(3,1)=B(1,2)*B(2,3)-B(1,3)*B(2,2)
  29. B(3,2)=B(1,3)*B(2,1)-B(1,1)*B(2,3)
  30. B(3,3)=B(1,1)*B(2,2)-B(1,2)*B(2,1)
  31. *
  32. ELSE IF (IDIM.EQ.2) THEN
  33. VX1=WORK(1)
  34. VY1=WORK(2)
  35. VL1=SQRT(VX1*VX1+VY1*VY1)
  36. *
  37. B(1,1)=VX1/VL1
  38. B(1,2)=VY1/VL1
  39. B(1,3)=0
  40. B(2,1)=B(1,2)*(-1.D0)
  41. B(2,2)=B(1,1)
  42. B(2,3)=0
  43. B(3,1)=0
  44. B(3,2)=0
  45. B(3,3)=1
  46. ENDIF
  47. *
  48. C
  49. * WRITE(*,*) ' matrice BPSS '
  50. * WRITE(*,1190) ((B(IA,IC),IA=1,3),IC=1,3)
  51. *1199 FORMAT(12(1X,E8.2))
  52. RETURN
  53. END
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  

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