Sources Esope | Cast3M

lispst
C LISPST    SOURCE    CHAT      05/01/13    01:23:52     5004
      SUBROUTINE LISPST(XE,CAR,D,XDDL,WORK,NPOINT,MELE,XSTRS,I69,I70,
     1 I195,I157)
C=======================================================================
C
C     CALCULE LES CONTRAINTES DU LINE SPRING
C
C     CODE EBERSOLT NOVEMBRE 84
C
C  ENTREES
C     XE(3,4)  =COORDONNEEES DU LINE SPRING
C     CAR(10)  =STOCKE LES CARACTERISTIQUES DU LINESPRING
C                CAR(1)  =EPAISSEUR NOEUDS 1 4
C                CAR(2)  =PROFONDEUR FISSURE NOEUDS 1 4
C                CAR(3,4,5) =VECTEU ORIENTANT LA FISSSURE NDS 1 4
C                CAR(6)  =EPAISSEUR NOEUDS 2 3
C                CAR(7)  =PROFONDEUR FISSURE NOEUDS 2 3
C                CAR(8,9,10)=VECTEU ORIENTANT LA FISSSURE NDS 2 3
C    D(2,2)    =MATRICE DE HOOKE
C    XDDL(24)  =D.D.L. REPERE GLOBAL
C    NPOINT    = NOMBRE DE POINTS POUR LES CONTRAINTES
C    MELE      = NUMERO  DE L ELEMENT 30 SI LISP   50 LISP
C TRAVAIL
C    WORK(42)  =TABLEAU DE TRAVAIL
C SORTIES
C    XSTRS(3*NPOINT)=  CONTRAINTES  DONT K I
C
C=======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION XE(3,*),CAR(*),D(6,*),XDDL(*),WORK(*),XSTRS(*)
      DIMENSION V1(3),V2(3)
      EPA1 =CAR(1)
      FIS1 =CAR(2)
      V1(1)=CAR(3)
      V1(2)=CAR(4)
      V1(3)=CAR(5)
      EPA2 =CAR(6)
      FIS2 =CAR(7)
      V2(1)=CAR(8)
      V2(2)=CAR(9)
      V2(3)=CAR(10)
C
      IF(MELE.EQ.30) THEN
      CALL LISPS2(XE,EPA1,FIS1,V1,EPA2,FIS2,V2,D,XDDL,
     1   WORK(1),WORK(10),WORK(19),NPOINT,XSTRS,I69,I70,I195,I157)
C
      ELSE IF(MELE.EQ.50) THEN
      CALL LISPS3(XE,EPA1,FIS1,V1,EPA2,FIS2,V2,D,XDDL,
     1   WORK(1),WORK(10),WORK(19),NPOINT,XSTRS,I69,I70,I195,I157)
      ENDIF
C
      RETURN
      END