Télécharger bdbst.eso

Retour à la liste

Numérotation des lignes :

  1. C BDBST SOURCE CHAT 05/01/12 21:37:58 5004
  2.  
  3. C=======================================================================
  4. C= B D B S T =
  5. C= --------- =
  6. C= =
  7. C= Fonction : =
  8. C= ---------- =
  9. C= Calcul de la contribution d'un point de Gauss a la matrice de =
  10. C= RIGIDITE elastique REL d'un element fini. =
  11. C= Contribution : Jacobien.transpose(BGENE)*DHOOK*BGENE =
  12. C= =
  13. C= Parametres : (E)=Entree (S)=Sortie =
  14. C= ------------ =
  15. C= BGENE (E) Matrice B =
  16. C= DJac (E) "Volume" associe au point de Gauss =
  17. C= DHOOK (E) Matrice de Hooke au point de Gauss =
  18. C= LRE (E) Nombre de DDL de l'element fini =
  19. C= NSTRS (E) Nombre de composantes de contraintes/deformations =
  20. C= REL (E/S) Matrice de rigidite de l'element fini =
  21. C= =
  22. C= Remarque : Actuellement, NSTRS doit etre inferieur ou egal a 10 ! =
  23. C=======================================================================
  24.  
  25. SUBROUTINE BDBST(BGENE,DJAC,DHOOK,LRE,NSTRS,REL)
  26.  
  27. IMPLICIT INTEGER(I-N)
  28. IMPLICIT REAL*8 (A-H,O-Z)
  29.  
  30. -INC CCREEL
  31.  
  32. DIMENSION DHOOK(NSTRS,*),REL(LRE,*),BGENE(NSTRS,*)
  33. DIMENSION WORK1(10),WORK2(10)
  34.  
  35. C= Boucle sur les DDL
  36. DO IX=1,LRE
  37. C== Calcul de JACobien*Poids*B pour le DDL IX
  38. DO i=1,NSTRS
  39. WORK1(i)=DJAC*BGENE(i,IX)
  40. ENDDO
  41. C== Calcul de Vol.*B*D pour le DDL IX
  42. DO i=1,NSTRS
  43. CC=XZero
  44. DO j=1,NSTRS
  45. CC=CC+WORK1(j)*DHOOK(i,j)
  46. ENDDO
  47. WORK2(i)=CC
  48. ENDDO
  49. C== Calcul de Vol.Bt.D.B pour les DDL IX et IY
  50. DO IY=1,IX
  51. CC=XZero
  52. DO i=1,NSTRS
  53. CC=CC+BGENE(i,IY)*WORK2(i)
  54. ENDDO
  55. REL(IX,IY)=REL(IX,IY)+CC
  56. REL(IY,IX)=REL(IX,IY)
  57. ENDDO
  58. ENDDO
  59.  
  60. RETURN
  61. END
  62.  
  63.  
  64.  
  65.  

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