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    CHAT      05/01/12    21:33:47     5004
  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.   1   V(I)=XE(I,3)-XE(I,1)
  26.       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  27.       XL=SQRT(XL)
  28.       IF(XL.EQ.0.) GO TO 999
  29. *
  30.       DO 2 I=1,3
  31.       U(I)=0.D0
  32.       DO 3 J=1,IDIM
  33.        IG=(I-1)*IDIM+J
  34.        U(I)=U(I)+(DEPLA(IG)*V(J))/XL
  35.   3   CONTINUE
  36.   2   CONTINUE
  37. *
  38.       DO 5 IGAU=1,NBPGAU
  39.          X1=XL/2.D0*(QSIGAU(IGAU)+1.D0)
  40.          B3(1)=(4.D0*X1/XL-3.D0)/XL
  41.          B3(2)=(4.D0-8.D0*X1/XL)/XL
  42.          B3(3)=(4.D0*X1/XL-1.D0)/XL
  43.          DEFO(IGAU)=B3(1)*U(1)+B3(2)*U(2)+B3(3)*U(3)
  44.   5   CONTINUE
  45. C--------------------------------------------
  46. *      NRE=LRE
  47. *      WRITE(*,*) '**  DEPLA GLOBAL BAR3 (dans barep3)** ELEME:',IB
  48. *      WRITE(*,1991) (DEPLA(IC),IC=1,NRE)
  49. *      WRITE(*,*) '**  U LOCAL BAR3 (dans barep3)** '
  50. *      WRITE(*,1991) (U(IC),IC=1,3)
  51. *      WRITE(*,*) '**  QSIGAU   NBPGAU=',NBPGAU,'XL=',XL
  52. *      WRITE(*,1991) (QSIGAU(IC),IC=1,NBPGAU)
  53. *      WRITE(*,*) '**  POIGAU   '
  54. *      WRITE(*,1991) (POIGAU(IC),IC=1,NBPGAU)
  55. *      WRITE(*,*) '**  B BAR3 (dans barep3)** '
  56. *      WRITE(*,1991) (B3(IC),IC=1,3)
  57. *      WRITE(*,*) '**  DEFO BAR3 (dans barep3)** '
  58. *      WRITE(*,1991) (DEFO(IC),IC=1,3)
  59. *1991  FORMAT(18(1X,E12.5))
  60. C--------------------------------------------
  61.  100  CONTINUE
  62.       RETURN
  63.  999  KERRE=1
  64.       RETURN
  65.       END
  66.  
  67.  
  68.  

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