Télécharger bjo2.eso

Retour à la liste

Numérotation des lignes :

  1. C BJO2 SOURCE CHAT 05/01/12 21:40:22 5004
  2. SUBROUTINE BJO2(IGAU,MFR,IFOU,NIFOU,XEL,BPSS,SHPTOT,SHP,
  3. . BGENE,DJAC,IRRT)
  4. C=======================================================================
  5. C
  6. C CALCUL DE LA MATRICE B = TETA * ( N , -N )
  7. C ET DU JACOBIEN EN IGAU POUR UN JOINT 2D
  8. C ROUTINE FORTRAN PUR
  9. C CODE S. FELIX MAI 92
  10. C=======================================================================
  11. C INPUT
  12. C IGAU = NUMERO DU POINT DE GAUSS
  13. C MFR = NUMERO DE LA FORMULATION
  14. C XEL = COORDONNEES DES NOEUDS DE L'ELEMENT
  15. C BPSS = MATRICE DE PASSAGE
  16. C BPSS(,) = AXE S1
  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 : CE SOUS-PROGRAMME EST VALABLE EN CONTRAINTES PLANES,
  36. C DEFORMATIONS PLANES ET AXISYMETRIQUES. LA MATRICE
  37. C BGENE = TETA ( Ntop , - Nbot )
  38. C EST UNE MATRICE (2x8). EN EFFET, TETA EST UNE MATRICE
  39. C (2x2), ET ( Ntop , - Nbot ) EST UNE MATRICE (2x8).
  40. C LA MATRICE TETA EST ICI APPELEE BPSS.
  41. C
  42. C=======================================================================
  43. IMPLICIT INTEGER(I-N)
  44. IMPLICIT REAL*8(A-H,O-Z)
  45. PARAMETER (LRE=8,NST=2,NBNO=4)
  46. DIMENSION XEL(3,*),BGENE(NST,*),SHP(6,*),SHPTOT(6,NBNO,*)
  47. DIMENSION BPSS(3,3)
  48. C
  49. IRRT = 0
  50. C
  51. C MATRICE JACOBIENNE
  52. C
  53. DO 1 I=1,NBNO
  54. SHP(1,I) = SHPTOT(1,I,IGAU)
  55. SHP(2,I) = SHPTOT(2,I,IGAU)
  56. 1 CONTINUE
  57. C
  58. C !!!!!! ATTENTION : IL FAUT CALCULER LE JACOBIEN AVEC NBNO=2 !!!!!!
  59. C IL NE FAUT SURTOUT PAS METTRE NBNO=4 CAR CA FAUSSE ALORS
  60. C LES RESULTATS ... NBNO DOIT ETRE EGAL A 4 CAR IL YA 4 NOEUDS
  61. C
  62. NBNONN=NBNO/2
  63. CALL DEVOLU(XEL,SHP,MFR,NBNONN,IFOU,NIFOU,1,1.D0,RR,DJAC)
  64.  
  65. IF (DJAC.LT.0.0D0) THEN
  66. IRRT = 1
  67. ELSE IF (DJAC.EQ.0.0D0) THEN
  68. IRRT = 2
  69. ENDIF
  70. C
  71. C MATRICE B
  72. C
  73. C PP ATTENTION: IL FAUT PEUT-ETRE RESPECTER LA NUMEROTATION LOCALE DES
  74. C NOEUDS
  75. C
  76. CALL ZERO(BGENE,NST,LRE)
  77. DO 2 I=1,2
  78. DO 3 J=1,2
  79. DO 4 K=1,2
  80. L=2*(J-1)+K
  81. C PP M=L+4
  82. M=L+2*(2*(2-J)+1)
  83. BGENE(I,L)=BPSS(I,K)*SHP(1,J)
  84. BGENE(I,M)=-BGENE(I,L)
  85. 4 CONTINUE
  86. 3 CONTINUE
  87. 2 CONTINUE
  88. C
  89. RETURN
  90. END
  91.  
  92.  
  93.  

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