Télécharger bgfcq4.eso

Retour à la liste

Numérotation des lignes :

  1. C BGFCQ4 SOURCE CHAT 06/03/29 21:15:48 5360
  2. SUBROUTINE BGFCQ4(IGAU,XEL,SHPTOT,SHP,BGF,DJAC,IARR)
  3. C |====================================================================|
  4. C | |
  5. C | CALCULE LA MATRICE DE GRADIAN EN FLEXION BGF POUR COQ4
  6. C | CODE SUO X.Z. MARS 87
  7. C |INPUT |
  8. C | IGAU=NUMERO DU POINT DE GAUSS |
  9. C | NBNO=NOMBRE DE NOEUDS |
  10. C | LRE =NOMBRE DE COLONNES DE LA MATRICE BGR |
  11. C | XEL =COORDONNEES LOCALES DE L ELEMENT |
  12. C | SHPTOT(6,NBNO,NBGAU)=FONCTIONS DE FORMES ET DERIVEES |
  13. C |ZONE DE TRAVAIL |
  14. C | SHP(6,NBNO)=TABLEAU DE TRAVAIL |
  15. C |OUTPUT |
  16. C | DJAC=JACOBIEN |
  17. C | BGF(9,LRE)=MATRICE BGF |
  18. C |====================================================================|
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21. PARAMETER (LRE=24,NBNO=4)
  22. DIMENSION XEL(3,*),BGF(9,*),SHP(6,*),SHPTOT(6,NBNO,*)
  23. DATA ZER/0.D0/
  24. IARR=0
  25. C
  26. CALL ZERO(BGF,9,LRE)
  27. C
  28. C
  29. DO 101 NP=1,NBNO
  30. SHP(1,NP)=SHPTOT(1,NP,IGAU)
  31. SHP(2,NP)=SHPTOT(2,NP,IGAU)
  32. SHP(3,NP)=SHPTOT(3,NP,IGAU)
  33. 101 CONTINUE
  34. CALL JACOBI(XEL,SHP,2,NBNO,DJAC)
  35. IF(DJAC.LE.ZER) IARR=1
  36. IF(DJAC.LE.ZER) RETURN
  37. C
  38. C ====================== CALCUL
  39. C
  40. K=1
  41. DO 102 NP=1,NBNO
  42. BGF(1,K+4)=SHP(2,NP)
  43. BGF(2,K+4)=SHP(3,NP)
  44. BGF(3,K+4)=SHP(1,NP)
  45. BGF(4,K+3)=-SHP(2,NP)
  46. BGF(5,K+3)=-SHP(3,NP)
  47. BGF(6,K+3)=-SHP(1,NP)
  48. K=K+6
  49. 102 CONTINUE
  50. RETURN
  51. END
  52.  
  53.  
  54.  

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