Télécharger barep3.eso

Retour à la liste

Numérotation des lignes :

  1. C BAREP3 SOURCE CHAT 05/01/12 21:33:47 5004
  2. C BAREP3 SOURCE DJER 89/05/26 22:54:40
  3. SUBROUTINE BAREP3(XE,DEPLA,DEFO,QSIGAU,POIGAU,NBPGAU,IB)
  4. C-----------------------------------------------------------------------
  5. C ROUTINE DE CALCUL DES DEFORMATIONS POUR LA BARRE
  6. C
  7. C ENTREE
  8. C DEPLA DEPLACEMENTS
  9. C XE XEONNEES NOEUDS
  10. C
  11. C SORTIE
  12. C DEFO DEFORMATION AUX POINTS D INTEGRATION
  13. C-----------------------------------------------------------------------
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16. -INC CCOPTIO
  17. DIMENSION DEFO(*),XE(3,*),V(3),B3(3),DEPLA(*),U(3)
  18. DIMENSION QSIGAU(*),POIGAU(*)
  19. C
  20. KERRE=0
  21. V(3)=0.D0
  22. DO 1 I=1,IDIM
  23. 1 V(I)=XE(I,3)-XE(I,1)
  24. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  25. XL=SQRT(XL)
  26. IF(XL.EQ.0.) GO TO 999
  27. *
  28. DO 2 I=1,3
  29. U(I)=0.D0
  30. DO 3 J=1,IDIM
  31. IG=(I-1)*IDIM+J
  32. U(I)=U(I)+(DEPLA(IG)*V(J))/XL
  33. 3 CONTINUE
  34. 2 CONTINUE
  35. *
  36. DO 5 IGAU=1,NBPGAU
  37. X1=XL/2.D0*(QSIGAU(IGAU)+1.D0)
  38. B3(1)=(4.D0*X1/XL-3.D0)/XL
  39. B3(2)=(4.D0-8.D0*X1/XL)/XL
  40. B3(3)=(4.D0*X1/XL-1.D0)/XL
  41. DEFO(IGAU)=B3(1)*U(1)+B3(2)*U(2)+B3(3)*U(3)
  42. 5 CONTINUE
  43. C--------------------------------------------
  44. * NRE=LRE
  45. * WRITE(*,*) '** DEPLA GLOBAL BAR3 (dans barep3)** ELEME:',IB
  46. * WRITE(*,1991) (DEPLA(IC),IC=1,NRE)
  47. * WRITE(*,*) '** U LOCAL BAR3 (dans barep3)** '
  48. * WRITE(*,1991) (U(IC),IC=1,3)
  49. * WRITE(*,*) '** QSIGAU NBPGAU=',NBPGAU,'XL=',XL
  50. * WRITE(*,1991) (QSIGAU(IC),IC=1,NBPGAU)
  51. * WRITE(*,*) '** POIGAU '
  52. * WRITE(*,1991) (POIGAU(IC),IC=1,NBPGAU)
  53. * WRITE(*,*) '** B BAR3 (dans barep3)** '
  54. * WRITE(*,1991) (B3(IC),IC=1,3)
  55. * WRITE(*,*) '** DEFO BAR3 (dans barep3)** '
  56. * WRITE(*,1991) (DEFO(IC),IC=1,3)
  57. *1991 FORMAT(18(1X,E12.5))
  58. C--------------------------------------------
  59. 100 CONTINUE
  60. RETURN
  61. 999 KERRE=1
  62. RETURN
  63. END
  64.  
  65.  
  66.  

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