Télécharger bdbsts.eso

Retour à la liste

Numérotation des lignes :

  1. C BDBSTS SOURCE CHAT 05/01/12 21:38:00 5004
  2. SUBROUTINE BDBSTS(BGENE,DJAC,DHOOK,LRE,NSTRS,REL)
  3. C=======================================================================
  4. C
  5. C
  6. C INTEGRE LE PRODUIT BDB
  7. C "EXTENSION" A PARTIR DE BDBST AU CAS NON-SYMETRIQUE
  8. C WARNING INB<-->INC
  9. C
  10. C VERSION SEPTEMBRE 92 PIERRE PEGON
  11. C ROUTINE FORTRAN PUR
  12. C
  13. C=======================================================================
  14. C BGENE(NSTRS,LRE)=MATRICE B
  15. C DJAC =JACOBIEN
  16. C DHOOK(NSTRS,NSTRS)=MATRICE DE HOOKE
  17. C LRE=TAILLE DE LA MATRICE DE RIGIDITE
  18. C REL=MATRICE DE RIGIDITE
  19. C LES KIJ SONT AJOUTES SUR REL
  20. C=======================================================================
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8 (A-H,O-Z)
  23. DIMENSION DHOOK(NSTRS,1),REL(LRE,1)
  24. DIMENSION WORK1(10),WORK2(10),BGENE(NSTRS,1)
  25. C
  26. DO 301 IX=1,LRE
  27. C
  28. C BOUCLE SUR LES DDL
  29. C
  30. DO 401 INA=1,NSTRS
  31. C
  32. C CALCUL DE JACOBIAN*POIDS*B POUR LE DDL IX
  33. C
  34. WORK1(INA)=DJAC*BGENE(INA,IX)
  35. 401 CONTINUE
  36. DO 501 INB=1,NSTRS
  37. CC=0.D0
  38. DO 601 INC=1,NSTRS
  39. CC=CC+WORK1(INC)*DHOOK(INC,INB)
  40. 601 CONTINUE
  41. WORK2(INB)=CC
  42. 501 CONTINUE
  43. DO 701 IY=1,LRE
  44. C
  45. C CALCUL DE BDB POUR LES DDL IX ET IY
  46. C
  47. CC=0.D0
  48. DO 801 IND=1,NSTRS
  49. CC=CC+BGENE(IND,IY)*WORK2(IND)
  50. 801 CONTINUE
  51. REL(IX,IY)=REL(IX,IY)+CC
  52. 701 CONTINUE
  53. 301 CONTINUE
  54. C
  55. RETURN
  56. END
  57.  
  58.  
  59.  

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