Télécharger la plaquette Cast3M


fr en

Télécharger tipocq.eso

Retour à la liste

Numérotation des lignes : 

tipocq

  1. C TIPOCQ    SOURCE    CHAT      05/01/13    03:41:06     5004
  2.       SUBROUTINE TIPOCQ(XE, C,Q,MECAN,NCOMP2)
  3. ***********
  4. *
  5. *     ESCLAVE DE MENISM
  6. *
  7. *     P.PEGON (ISPRA) AOUT 1996
  8. *
  9. ***********
  10. *
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8 (A-H,O-Z)
  13.       DIMENSION XE(3,2)
  14.       DIMENSION C(MECAN,NCOMP2),Q(NCOMP2,NCOMP2)
  15. *
  16.       CALL ZERO(C,MECAN,NCOMP2)
  17.       CALL ZERO(Q,NCOMP2,NCOMP2)
  18. *
  19.       NCOMP=NCOMP2/2
  20.       IF(MECAN.EQ.3.AND.NCOMP.EQ.3)THEN
  21.         XL=SQRT((XE(1,2)-XE(1,1))**2+(XE(2,2)-XE(2,1))**2)
  22.         C(1,2)=-1.D0/XL
  23.         C(1,3)= 1.D0
  24.         C(1,5)=+1.D0/XL
  25.         C(2,2)=-1.D0/XL
  26.         C(2,5)=+1.D0/XL
  27.         C(2,6)= 1.D0
  28.         C(3,1)=-1.D0
  29.         C(3,4)= 1.D0
  30. *
  31.         COSDIR=(XE(1,2)-XE(1,1))/XL
  32.         SINDIR=(XE(2,2)-XE(2,1))/XL
  33.         DO I=0,1
  34.           II=I*NCOMP
  35.           Q(II+1,II+1)= COSDIR
  36.           Q(II+1,II+2)= SINDIR
  37.           Q(II+2,II+1)=-SINDIR
  38.           Q(II+2,II+2)= COSDIR
  39.           Q(II+3,II+3)= 1
  40.         ENDDO
  41.       ENDIF
  42. *
  43.       RETURN
  44.       END
  45.  
  46.  
  47.  

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