Télécharger la plaquette Cast3M


fr en

Télécharger barep3.eso

Retour à la liste

Numérotation des lignes : 

barep3

  1. C BAREP3    SOURCE    CB215821  26/10/26    21:15:11     12682          
  2. C BAREP3    SOURCE    DJER      89/05/26    22:54:40
  3.       SUBROUTINE BAREP3(XE,DEPLA,DEFO,QSIGAU,POIGAU,NBPGAU,IB)
  4. C-----------------------------------------------------------------------
  5. C  ROUTINE DE CALCUL DES DEFORMATIONS POUR LA BARRE
  6. C
  7. C  ENTREE
  8. C          DEPLA   DEPLACEMENTS
  9. C          XE   XEONNEES NOEUDS
  10. C
  11. C  SORTIE
  12. C          DEFO    DEFORMATION AUX POINTS D INTEGRATION
  13. C-----------------------------------------------------------------------
  14.       IMPLICIT INTEGER(I-N)
  15.       IMPLICIT REAL*8(A-H,O-Z)
  16.  
  17. -INC PPARAM
  18. -INC CCOPTIO
  19.       DIMENSION DEFO(*),XE(3,*),V(3),B3(3),DEPLA(*),U(3)
  20.       DIMENSION QSIGAU(*),POIGAU(*)
  21. C
  22.       KERRE=0
  23.       V(3)=0.D0
  24.       DO 1 I=1,IDIM
  25.       V(I)=XE(I,3)-XE(I,1)
  26.   1   CONTINUE
  27.       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  28.       XL=SQRT(XL)
  29.       IF(XL.EQ.0.) GO TO 999
  30. *
  31.       DO 2 I=1,3
  32.       U(I)=0.D0
  33.       DO 3 J=1,IDIM
  34.        IG=(I-1)*IDIM+J
  35.        U(I)=U(I)+(DEPLA(IG)*V(J))/XL
  36.   3   CONTINUE
  37.   2   CONTINUE
  38. *
  39.       DO 5 IGAU=1,NBPGAU
  40.          X1=XL/2.D0*(QSIGAU(IGAU)+1.D0)
  41.          B3(1)=(4.D0*X1/XL-3.D0)/XL
  42.          B3(2)=(4.D0-8.D0*X1/XL)/XL
  43.          B3(3)=(4.D0*X1/XL-1.D0)/XL
  44.          DEFO(IGAU)=B3(1)*U(1)+B3(2)*U(2)+B3(3)*U(3)
  45.   5   CONTINUE
  46. C--------------------------------------------
  47. *      NRE=LRE
  48. *      WRITE(*,*) '**  DEPLA GLOBAL BAR3 (dans barep3)** ELEME:',IB
  49. *      WRITE(*,1991) (DEPLA(IC),IC=1,NRE)
  50. *      WRITE(*,*) '**  U LOCAL BAR3 (dans barep3)** '
  51. *      WRITE(*,1991) (U(IC),IC=1,3)
  52. *      WRITE(*,*) '**  QSIGAU   NBPGAU=',NBPGAU,'XL=',XL
  53. *      WRITE(*,1991) (QSIGAU(IC),IC=1,NBPGAU)
  54. *      WRITE(*,*) '**  POIGAU   '
  55. *      WRITE(*,1991) (POIGAU(IC),IC=1,NBPGAU)
  56. *      WRITE(*,*) '**  B BAR3 (dans barep3)** '
  57. *      WRITE(*,1991) (B3(IC),IC=1,3)
  58. *      WRITE(*,*) '**  DEFO BAR3 (dans barep3)** '
  59. *      WRITE(*,1991) (DEFO(IC),IC=1,3)
  60. *1991  FORMAT(18(1X,E12.5))
  61. C--------------------------------------------
  62.       CONTINUE
  63.       RETURN
  64.  999  KERRE=1
  65.       RETURN
  66.       END
  67.  
  68.  
  69.  
  70.  

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