Télécharger dyne.eso

Retour à la liste

Numérotation des lignes :

  1. C DYNE SOURCE BP208322 18/01/25 21:15:39 9710
  2. SUBROUTINE DYNE
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Sous-programme de l'operateur DYNE *
  8. * __________________________________ *
  9. * *
  10. * Calcul d'une reponse dynamique a l'aide d'algorithmes *
  11. * explicites pour les forces non-lineaire: *
  12. * *
  13. * 1. algorithme de Fu - de Vogelaere *
  14. * 2. algorithme des differences centrees *
  15. * 3. algorithme de Newmark (implicite sauf pour les forces NL *
  16. * *
  17. * auteurs: Denis ROBERT-MOUGIN & Lionel VIVAN ,le 19 mai 1989. *
  18. * derniere modif: BP, 2018-01-16 *
  19. * *
  20. *--------------------------------------------------------------------*
  21. -INC CCOPTIO
  22. *
  23. CHARACTER*6 MOPLEX
  24. CHARACTER*20 MODEVO,MOALG,MODIFF,MORK4
  25. CHARACTER*20 MONEWM(2)
  26. DATA MOPLEX /'PLEXUS'/
  27. DATA MODEVO /'DE_VOGELAERE'/
  28. DATA MODIFF /'DIFFERENCES_CENTREES'/
  29. DATA MONEWM /'ACCELERATION_MOYENNE','FOX_GOODWIN'/
  30. c 12345678901234567890 12345678901
  31. DATA MORK4 /'RUNGE_KUTTA_4'/
  32. *
  33. * Quel algorithme a-t-on choisi ?
  34. *
  35. CALL MESLIR(-272)
  36. CALL LIRCHA(MOALG,1,IRET)
  37. IF (IRET.EQ.0) RETURN
  38. *
  39. * Algorithme de Plexus
  40. * on s'en sert plus. Il ne figure pas dans la notice
  41. IF (MOALG(1:6).EQ.MOPLEX) THEN
  42. c CALL DYNPLX
  43. CALL ERREUR(477)
  44. RETURN
  45.  
  46. * Algorithme de Fu - de Vogelaere
  47. ELSE IF (MOALG(1:12).EQ.MODEVO) THEN
  48. CALL DYNDEV
  49. RETURN
  50.  
  51. * Algorithme des differences centrees (le meme que plexus)
  52. ELSE IF (MOALG(1:20).EQ.MODIFF) THEN
  53. CALL DYNDIF(2)
  54. RETURN
  55.  
  56. * Algorithmes implicites de Newmark
  57. ELSE IF (MOALG(1:20).EQ.MONEWM(1)) THEN
  58. CALL DYNDIF(3)
  59. RETURN
  60. ELSE IF (MOALG(1:11).EQ.MONEWM(2)) THEN
  61. CALL DYNDIF(4)
  62. RETURN
  63.  
  64. * Erreur dans le choix de l'algorithme
  65. ELSE
  66. CALL ERREUR(477)
  67. RETURN
  68. ENDIF
  69. *
  70. END
  71.  
  72.  
  73.  
  74.  

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