Télécharger devfxa.eso

Retour à la liste

Numérotation des lignes :

devfxa
  1. C DEVFXA SOURCE BP208322 18/12/20 21:15:31 10048
  2. * SUBROUTINE DEVFXA(FEXA,IFEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS,
  3. SUBROUTINE DEVFXA(FEXA,FTOTA,NPFEXA,NA1,NPC1,NPAS,
  4. & FTEXB,FEXB,NPLB,IDIMB,RIGIDE)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. *--------------------------------------------------------------------*
  8. * *
  9. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  10. * ________________________________________________ *
  11. * *
  12. * Ajout des termes de forces ext{rieures libres. *
  13. * *
  14. * Param}tres: *
  15. * *
  16. * e FEXA Evolution des forces extérieures, en base A *
  17. * e IFEXA Numéro du mode correspondant au point de chargement *
  18. * (supprime le 2018-12-14 par bp)
  19. * es FTOTA Forces extérieures libres totalisées, en base A *
  20. * e NPFEXA Nombre de points de chargement *
  21. * e NA1 Nombre total d'inconnues dans la base A *
  22. * e NPC1 Nombre de pas de calcul *
  23. * e NPAS Num{ro du pas courant *
  24. * e FTEXB Evolution des forces extérieures sur base B, pour le *
  25. * calcul des moments *
  26. * es FEXB Forces extérieures sur base B, pour le calcul des *
  27. * moments * *
  28. * *
  29. * Auteur, date de création: *
  30. * *
  31. * Denis ROBERT-MOUGIN, le 22 mai 1989. *
  32. * *
  33. *--------------------------------------------------------------------*
  34. REAL*8 FEXA(NPFEXA,NPC1,*),FTOTA(NA1,*)
  35. REAL*8 FEXB(NPLB,2,*),FTEXB(NPLB,NPC1,2,*)
  36. c INTEGER IFEXA(*)
  37. *
  38. LOGICAL RIGIDE
  39. *
  40. * Boucle sur les points de chargement:
  41. *
  42. NPA = NPAS + 1
  43. DO 100 I=1,NPFEXA
  44. FTOTA(I,1) = FEXA(I,NPA,1)
  45. FTOTA(I,2) = FEXA(I,NPA,2)
  46. 100 CONTINUE
  47. * END DO
  48. *
  49. * Dans le cas des corps rigides, boucle sur les chargements Base B
  50. IF (RIGIDE) THEN
  51. IF (IDIMB.EQ.6) THEN
  52. IDIM = 3
  53. ELSE
  54. IDIM = 2
  55. ENDIF
  56. DO 200 IP=1,NPLB
  57. DO 210 ID=1,IDIM
  58. FEXB(IP,1,ID) = FTEXB(IP,NPA,1,ID)
  59. FEXB(IP,2,ID) = FTEXB(IP,NPA,2,ID)
  60. 210 CONTINUE
  61. 200 CONTINUE
  62. ENDIF
  63. RETURN
  64. END
  65.  
  66.  
  67.  
  68.  

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