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    CB215821  26/10/26    21:16:17     12682          
  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.       ZE(I)=COORD(3,I)
  34.   2   CONTINUE
  35.       XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2
  36.       XL=SQRT(XL)
  37.       IF(XL.EQ.0.) GO TO 999
  38. C
  39. C  VALEURS MOYENNES DANS L ELEMENT
  40. C
  41. C FX
  42.       WORK(1) =-CONT(1)
  43.       WORK(7) = CONT(1)
  44. C FY
  45.       WORK(2) =-CONT(2)
  46.       WORK(8) = CONT(2)
  47. C FZ
  48.       WORK(3) =-CONT(3)
  49.       WORK(9) = CONT(3)
  50. C MX
  51.       WORK(4) =-CONT(4)
  52.       WORK(10)= CONT(4)
  53. C MY
  54.       WORK(5) =-CONT(5)  + CONT(3)*P5*XL
  55.       WORK(11)= CONT(5)  + CONT(3)*P5*XL
  56. C MZ
  57.       WORK(6) =-CONT(6)  - CONT(2)*P5*XL
  58.       WORK(12)= CONT(6)  - CONT(2)*P5*XL
  59. C
  60. C  CALCUL DE LA MATRICE DE PASSAGE
  61. C
  62.       CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)
  63.       IF(KERRE.NE.0) RETURN
  64. C
  65. C  PASSAGE DES FORCES EN AXES GLOBAUX
  66. C
  67.       CALL POUVEC(WORK,FORCE,P,2)
  68.       RETURN
  69. 999   KERRE=1
  70.       RETURN
  71.       END
  72.  
  73.  
  74.  

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