Télécharger bjo4c.eso

Retour à la liste

Numérotation des lignes :

  1. C BJO4C SOURCE CHAT 05/01/12 21:40:33 5004
  2. SUBROUTINE BJO4C(IGAU,XEL,BPSS,SHPTOT,SHP,BGENE,DJAC,IRRT)
  3. C=======================================================================
  4. C
  5. C PPj POUR LE JOINT CISAILLEMENT JCI4 UNIQUEMENT 2D
  6. C CALCUL DE LA MATRICE B = TETA * ( N , -N )
  7. C ET DU JACOBIEN EN IGAU
  8. C ROUTINE FORTRAN PUR
  9. C P. PEGON 2001 d'apres CODE S. FELIX MAI 92
  10. C=======================================================================
  11. C INPUT
  12. C IGAU = NUMERO DU POINT DE GAUSS
  13. C XEL = COORDONNEES DES NOEUDS DE L'ELEMENT
  14. C BPSS = MATRICE DE PASSAGE
  15. C BPSS(,) = AXE S1
  16. C BPSS(,) = AXE S2
  17. C BPSS(,) = AXE SN
  18. C SHPTOT= FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE DE REFERENCE
  19. C SHPTOT(1, ) = FONCTION DE FORME
  20. C SHPTOT(2, ) = DERIVEES PAR RAPPORT A QSI
  21. C SHPTOT(3, ) = DERIVEES PAR RAPPORT A ETA
  22. C OUTPUT
  23. C SHP = FONCTIONS DE FORME ET DERIVEES DANS L'ESPACE GEOMETRIQUE
  24. C SHP(1, ) = FONCTION DE FORME
  25. C SHP(2, ) = DERIVEES PAR RAPPORT A X LOCAL
  26. C SHP(3, ) = DERIVEES PAR RAPPORT A Y LOCAL
  27. C DJAC = JACOBIEN AU POINT D'INTEGRATION IGAU
  28. C BGENE = MATRICE B AU POINT D'INTEGRATION IGAU
  29. C IRRT = DIFFERENT DE ZERO SI ERREUR
  30. C CONSTANTES
  31. C NBNO = NOMBRE DE NOEUDS
  32. C LRE = NOMBRE DE COLONNES DE LA MATRICE B
  33. C NST = NOMBRE DE COMPOSANTES DE CONTRAINTES
  34. C
  35. C REMARQUE : ATTENTION : DANS LES CAS CONTRAINTES PLANES ET
  36. C AXISYMETRIQUE, LA MATRICE BPSS ( RESPECTIVEMENT BGENE )
  37. C N'A PLUS LA DIMENSION (3X3) ( RESPECTIVEMENT (3X24) ).
  38. C TROUVER LEURS DIMENSIONS CORRECTES ET MODIFIER LES
  39. C PARAMETRES LRE, NST, NBNO. CEUX CI CORRESPONDRONT A
  40. C CEUX DU CAS BIDIMENSIONNELS. AINSI PAR EXEMPLE, NBNO
  41. C SERA EGAL A 4, NST SERA EGAL A 2 ETC ...
  42. C
  43. C=======================================================================
  44. IMPLICIT INTEGER(I-N)
  45. IMPLICIT REAL*8(A-H,O-Z)
  46. PARAMETER (LRE=2*8,NST=2,NBNO=8)
  47. DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)
  48. DIMENSION BPSS(3,3)
  49. C
  50. IRRT = 0
  51. C
  52. C MATRICE JACOBIENNE
  53. C
  54. DO 1 I=1,NBNO
  55. SHP(1,I) = SHPTOT(1,I,IGAU)
  56. SHP(2,I) = SHPTOT(2,I,IGAU)
  57. SHP(3,I) = SHPTOT(3,I,IGAU)
  58. 1 CONTINUE
  59. C
  60. C !!!!!! ATTENTION : IL FAUT CALCULER LE JACOBIEN AVEC NBNO=4 !!!!!!
  61. C IL NE FAUT SURTOUT PAS METTRE NBNO=8 CAR CA FAUSSE ALORS
  62. C LES RESULTATS ... NBNO DOIT ETRE EGAL A 8 CAR IL YA 8 NOEUDS
  63. C
  64. NBNONN=NBNO/2
  65. CALL JACOBI(XEL,SHP,2,NBNONN,DJAC)
  66. IF (DJAC.LT.0.0D0) THEN
  67. IRRT = 1
  68. ELSE IF (DJAC.EQ.0.0D0) THEN
  69. IRRT = 2
  70. ENDIF
  71. C
  72. C MATRICE B
  73. C
  74. CALL ZERO(BGENE,NST,LRE)
  75. DO 2 I=1,2
  76. DO 3 J=1,4
  77. DO 4 K=1,2
  78. L=2*(J-1)+K
  79. M=L+8
  80. BGENE(I,L)=BPSS(I,K)*SHP(1,J)
  81. BGENE(I,M)=-BGENE(I,L)
  82. 4 CONTINUE
  83. 3 CONTINUE
  84. 2 CONTINUE
  85. C
  86. RETURN
  87. END
  88.  
  89.  
  90.  

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