Télécharger la plaquette Cast3M


fr en

Télécharger mapaex.eso

Retour à la liste

Numérotation des lignes : 

mapaex

  1. C MAPAEX    SOURCE    BL174031  17/08/01    21:15:01     9523           
  2. C MAPAEX     SOURCE    DJ        93/01/12    21:20:55     781
  3.       SUBROUTINE MAPAEX(XE,NBNN,WORK,AL,TGL,LRE,KERRE)
  4. C=======================================================================
  5. C     RECHERCHE DU REPERE LOCAL DE BAEX ( barre 3D excentre)
  6. C=======================================================================
  7.       IMPLICIT INTEGER(I-N)
  8.       IMPLICIT REAL*8(A-H,O-Z)
  9.       DIMENSION XE(3,NBNN),B(3,3),TGL(2,LRE),WORK(6)
  10. * ---------------------------------------------------
  11.       EXCZ=WORK(2)
  12.       EXCY=WORK(3)
  13.       VX3P=WORK(4)
  14.       VY3P=WORK(5)
  15.       VZ3P=WORK(6)
  16. C
  17.       DO 13 I=1,2
  18.       DO 13 J=1,LRE
  19.        TGL(I,J)=0.D0
  20.    13 CONTINUE
  21. C
  22.       VX1=XE(1,2)-XE(1,1)
  23.       VY1=XE(2,2)-XE(2,1)
  24.       VZ1=XE(3,2)-XE(3,1)
  25.       VL1=SQRT(VX1*VX1+VY1*VY1+VZ1*VZ1)
  26.       AL=VL1
  27. C
  28.       KERRE=0
  29.       IF(AL.EQ.0.D0) THEN
  30.          KERRE=1
  31.          GO TO 999
  32.       ENDIF
  33. C
  34.       VX2=VY3P*VZ1-VY1*VZ3P
  35.       VY2=VX1*VZ3P-VX3P*VZ1
  36.       VZ2=VX3P*VY1-VX1*VY3P
  37.       VL2=SQRT(VX2*VX2+VY2*VY2+VZ2*VZ2)
  38. C
  39. *     correction 1 BL : passage au repere local seg2
  40.       B(1,1)=VX1/VL1
  41.       B(1,2)=VY1/VL1
  42.       B(1,3)=VZ1/VL1
  43.       B(3,1)=-1*VX2/VL2
  44.       B(3,2)=-1*VY2/VL2
  45.       B(3,3)=-1*VZ2/VL2
  46.       B(2,1)=B(3,2)*B(1,3)-B(1,2)*B(3,3)
  47.       B(2,2)=B(1,1)*B(3,3)-B(3,1)*B(1,3)
  48.       B(2,3)=B(3,1)*B(1,2)-B(1,1)*B(3,2)
  49. *     fin correction 1 BL
  50. C
  51.       LIB=6
  52. C
  53.       DO 20 I=1,2
  54.       DO 10 J=1,3
  55.       L1=(I-1)*LIB+J
  56.       TGL(I,L1)=B(1,J)
  57.       IF(LIB.EQ.6) THEN
  58.         L2=I*LIB-2
  59.         L3=I*LIB-1
  60.         L4=I*LIB
  61. *     correction 2 BL : prise en compte des 2 excentrements
  62.         TGL(I,L2)= B(2,1)*EXCZ-B(3,1)*EXCY
  63.         TGL(I,L3)= B(2,2)*EXCZ-B(3,2)*EXCY
  64.         TGL(I,L4)= B(2,3)*EXCZ-B(3,3)*EXCY
  65. *     fin correction 2 BL 
  66.       ELSE
  67.         L2=I*LIB-1
  68.         L3=I*LIB
  69.         TGL(I,L2)= B(2,1)*EXCZ
  70.         TGL(I,L3)= B(2,2)*EXCZ
  71.       ENDIF
  72.    10 CONTINUE
  73.    20 CONTINUE
  74.  999  CONTINUE
  75. C
  76. *      WRITE(*,*) ' ======================== calcul TGL ds mapaex '
  77. *      WRITE(*,*) ' matrice  B  elem =',IB
  78. *      WRITE(*,1190) ((B(IA,IC),IA=1,3),IC=1,3)
  79. *      WRITE(*,*) ' XGENE  passage'
  80. *      DO 110 I =1,2
  81. *      WRITE(*,1199) (TGL(I,IC),IC=1,LRE)
  82. * 110  CONTINUE
  83. *1190   FORMAT(3(3(2X,E10.4),/))
  84. *1199   FORMAT(12(1X,E8.2))
  85.       RETURN
  86.       END
  87.  
  88.  
  89.  
  90.  
  91.  

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