Télécharger devpro.eso

Retour à la liste

Numérotation des lignes :

  1. C DEVPRO SOURCE CHAT 05/01/12 22:47:04 5004
  2. SUBROUTINE DEVPRO(XPHILB,FTOTB,FTOTBA,IBASB,INMSB,IPLSB,IORSB,
  3. & NSB,NPLSB,NA2,IDIMB,NPLB,NA1)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. *--------------------------------------------------------------------*
  7. * *
  8. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  9. * ________________________________________________ *
  10. * *
  11. * Projection des forces base r{elle sur base A. *
  12. * *
  13. * Param}tres: *
  14. * *
  15. * e XPHILB Tableau des vecteurs propres aux points de liaisons. *
  16. * e FTOTB Tableau des forces sur base B. *
  17. * s FTOTBA Tableau des forces base B projet{es sur base A. *
  18. * e IBASB Indique dans quelle sous base appartient le point de *
  19. * liaison. *
  20. * e INMSB Nombre de modes par sous base. *
  21. * e IORSB Donne l'indice du premier mode de la sous base dans *
  22. * l'ensemble des modes. *
  23. * e NSB Nombre total de sous base. *
  24. * e NPLSB Nombre total de points intervenant dans les liaisons *
  25. * d'une sous base. *
  26. * e NPLB Nombre total de points intervenant dans les liaisons. *
  27. * e IDIMB Nombre de ddl retenus. *
  28. * e NA1 Nombre total d'inconnues en base A. *
  29. * *
  30. * *
  31. * Auteur, date de cr{ation: *
  32. * *
  33. * Lionel VIVAN, le 22 Septembre 1989. *
  34. * *
  35. *--------------------------------------------------------------------*
  36. *
  37. INTEGER IBASB(*),INMSB(*),IPLSB(*),IORSB(*)
  38. REAL*8 XPHILB(NSB,NPLSB,NA2,*),FTOTB(NPLB,*),FTOTBA(*)
  39. *
  40. DO 10 IN = 1,NA1
  41. FTOTBA(IN) = 0.D0
  42. 10 CONTINUE
  43. * end do
  44. DO 20 IP = 1,NPLB
  45. ISB = IBASB(IP)
  46. IPLB = IPLSB(IP)
  47. NA3 = INMSB(ISB)
  48. INA2 = IORSB(ISB) - 1
  49. DO 30 IN = 1,NA3
  50. XRET = 0.D0
  51. DO 40 ID = 1,IDIMB
  52. XRET = XRET + XPHILB(ISB,IPLB,IN,ID) * FTOTB(IP,ID)
  53. 40 CONTINUE
  54. * end do
  55. IN2 = INA2 + IN
  56. FTOTBA(IN2) = FTOTBA(IN2) + XRET
  57. 30 CONTINUE
  58. * end do
  59. 20 CONTINUE
  60. * end do
  61. *
  62. END
  63.  
  64.  

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