Télécharger la plaquette Cast3M


fr en

Télécharger timbsg.eso

Retour à la liste

Numérotation des lignes : 

timbsg

  1. C TIMBSG    SOURCE    CHAT      05/01/13    03:40:04     5004
  2.       SUBROUTINE TIMBSG(FORCE,VECT,COORD,CONT,WORK,KERRE)
  3. C-----------------------------------------------------------------------
  4. C  ROUTINE DE CALCUL DES FORCES INTERNES B*SIGMA
  5. C             POUR LA POUTRE DROITE DE TIMOSHENKO SIMPLE ET MULTI-FIBRE
  6. C
  7. C  ENTREE
  8. C          CONT    CONTRAINTES
  9. C          VECT    ORIENTATION DE L'OY LOCAL
  10. C          COORD   COORDONNEES 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  TRAVAIL
  16. C          WORK(300)
  17. C
  18. C-----------------------------------------------------------------------
  19. C     CLONE DE POUBSG PP 4/9/92
  20. C-----------------------------------------------------------------------
  21. C     modifiee par Pierre Pegon (ISPRA) Juillet/Aout 1993
  22. C-----------------------------------------------------------------------
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25.       DIMENSION CONT(*),COORD(3,*),VECT(3),P(3,3),
  26.      .          XE(2),YE(2),ZE(2),FORCE(*),WORK(*)
  27.       DATA P5/0.5D0/
  28. C
  29.       KERRE=0
  30.       DO 2 I=1,2
  31.       XE(I)=COORD(1,I)
  32.       YE(I)=COORD(2,I)
  33.   2   ZE(I)=COORD(3,I)
  34.       XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2
  35.       XL=SQRT(XL)
  36.       IF(XL.EQ.0.) GO TO 999
  37. C
  38. C  VALEURS MOYENNES DANS L ELEMENT
  39. C
  40. C FX
  41.       WORK(1) =-CONT(1)
  42.       WORK(7) = CONT(1)
  43. C FY
  44.       WORK(2) =-CONT(2)
  45.       WORK(8) = CONT(2)
  46. C FZ
  47.       WORK(3) =-CONT(3)
  48.       WORK(9) = CONT(3)
  49. C MX
  50.       WORK(4) =-CONT(4)
  51.       WORK(10)= CONT(4)
  52. C MY
  53.       WORK(5) =-CONT(5)  + CONT(3)*P5*XL
  54.       WORK(11)= CONT(5)  + CONT(3)*P5*XL
  55. C MZ
  56.       WORK(6) =-CONT(6)  - CONT(2)*P5*XL
  57.       WORK(12)= CONT(6)  - CONT(2)*P5*XL
  58. C
  59. C  CALCUL DE LA MATRICE DE PASSAGE
  60. C
  61.       CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)
  62.       IF(KERRE.NE.0) RETURN
  63. C
  64. C  PASSAGE DES FORCES EN AXES GLOBAUX
  65. C
  66.       CALL POUVEC(WORK,FORCE,P,2)
  67.       RETURN
  68. 999   KERRE=1
  69.       RETURN
  70.       END
  71.  
  72.  

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