Sources Esope | Cast3M

hs456
C HS456     SOURCE    CHAT      05/01/13    00:26:29     5004
       SUBROUTINE HS456(XE,NSTRS,DDHOMU,HS4,HS5,HS6)
*
*       CALCUL DES FONCTIONS HS4,HS5 ET HS6
*       POUR L'ELEMENT COQUE DST
*
************************************************************************
*
* ENTREES :
*
*     XE : COORDONNEES DES NOEUDS
* NSTRS  :NOMBRE DE COMPOSANTES DE CONTRAINTES
* DDHOMU :MATRICE DE HOOKE
*
* SORTIES :
*
* HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
*
* AUTEUR : X.Z. SUO FEVRIER 96
************************************************************************
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*),HS6(*)
      DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3),HFX(9),HFY(9)
      DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
*
      DO 10 K=4,6
       IF(K.EQ.4)THEN
         IJ=1
         I=2
         J=3
       ELSEIF(K.EQ.5)THEN
         IJ=2
          I=3
          J=1
        ELSE
         IJ=3
          I=1
          J=2
        ENDIF
       SX(IJ)=XE(1,I)-XE(1,J)
       SY(IJ)=XE(2,I)-XE(2,J)
       RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
       EXX(IJ)=-SX(IJ)/RL(IJ)
       EYY(IJ)=-SY(IJ)/RL(IJ)
  10   CONTINUE
       AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
       SOMM=SX(IJ)+SY(IJ)+RL(IJ)+EXX(IJ)+EYY(IJ)
*
       CALL CISAF(XE,EXX,EYY,DDHOMU,AIR,ASS,AXX,
     1            AYY,BSS,BXX,BYY)
*
       CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
     1            HS4,HS5,HS6,HFX,HFY)
*
       RETURN
       END