Télécharger la plaquette Cast3M


fr en

Télécharger bdbst.eso

Retour à la liste

Numérotation des lignes : 

bdbst

  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