Télécharger barbs3.eso

Retour à la liste

Numérotation des lignes :

  1. C BARBS3 SOURCE CB215821 16/04/21 21:15:11 8920
  2. C BARBS3 SOURCE CHAT 89/05/26 22:56:13
  3. SUBROUTINE BARBS3(FORCE,XE,CONTG,KERRE,QSIGAU,POIGAU,NBPGAU,IB)
  4. C-----------------------------------------------------------------------
  5. C ROUTINE DE CALCUL DES FORCES INTERNES B*SIGMA
  6. C POUR LA BARRE
  7. C
  8. C ENTREE
  9. C EFFORT L'EFFORT DANS LA BARRE
  10. C XE XEONNEES DES 2 NOEUDS
  11. C
  12. C SORTIE
  13. C FORCE FORCES INTERNES
  14. C KERRE INDICE D'ERREUR ( 0 SI PAS DE PB )
  15. C
  16. C-----------------------------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. -INC CCOPTIO
  20. DIMENSION XE(3,*),V(3),FORCE(*),CONTG(*)
  21. DIMENSION QSIGAU(*),POIGAU(*),B3(3),FL(3)
  22. C
  23. KERRE=0
  24. V(3)=0.D0
  25. DO 1 I=1,IDIM
  26. 1 V(I)=XE(I,3)-XE(I,1)
  27. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  28. XL=SQRT(XL)
  29. IF(XL.LE.0.D0) GO TO 999
  30. *
  31. SECT=1.D0
  32. CALL ZERO(FL,3,1)
  33. DO 5 IGAU=1,NBPGAU
  34. ATDH=SECT*POIGAU(IGAU)*XL/2.D0
  35. X1=XL/2.D0*(QSIGAU(IGAU)+1.D0)
  36. B3(1)=(4.D0*X1/XL-3.D0)/XL
  37. B3(2)=(4.D0-8.D0*X1/XL)/XL
  38. B3(3)=(4.D0*X1/XL-1.D0)/XL
  39. DO 3 I=1,3
  40. FL(I)=FL(I)+B3(I)*CONTG(IGAU)*ATDH
  41. 3 CONTINUE
  42. 5 CONTINUE
  43. C
  44. C REMPLISSAGE DU VECTEUR FORCE
  45. C
  46. DO 2 IN=1,3
  47. DO 2 I=1,IDIM
  48. IG=(IN-1)*IDIM+I
  49. FORCE(IG)=V(I)*FL(IN)/XL
  50. 2 CONTINUE
  51. C--------------------------------------------
  52. * NRE=LRE
  53. * WRITE(*,*) '** B BAR3 (dans barbs3)** ELEME:',IB
  54. * WRITE(*,1991) (B3(IC),IC=1,3)
  55. * WRITE(*,*) '** CONTG BAR3 (dans barbs3)** '
  56. * WRITE(*,1991) (CONTG(IC),IC=1,3)
  57. * WRITE(*,*) '** FL LOCAL BAR3 (dans barbs3)** '
  58. * WRITE(*,1991) (FL(IC),IC=1,3)
  59. * WRITE(*,*) '** FORCE GLOBAL BAR3 (dans barbs3)** '
  60. * WRITE(*,1991) (FORCE(IC),IC=1,NRE)
  61. *1991 FORMAT(18(1X,E12.5))
  62. C--------------------------------------------
  63. RETURN
  64. 999 KERRE=1
  65. RETURN
  66. END
  67.  
  68.  
  69.  
  70.  
  71.  

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