Sources Esope | Cast3M

pokslo
C POKSLO    SOURCE    CB215821  26/08/24    21:17:44     12622          
      SUBROUTINE POKSLO(RP,XL,SIGM)
C-----------------------------------------------------------------------
C
C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE GEOMETRIQUE
C    ELEMENTAIRE D UNE POUTRE DROITE EN AXES LOCAUX
C    D APRES TEDEL
C
C  ENTREE
C    XL      LONGUEUR
C    SIGM    EFFORT NORMAL MOYEN
C  SORTIE
C    RP      MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE 12X12
C
C     MAI 85 MILLARD
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION RP(12,12)
C
C  INITIALISATIONS
C
      DO 121 J=1,12
      DO 120 I=1,12
      RP(I,J)=0.D00
  120 CONTINUE
  121 CONTINUE
      RL=1.D00/XL
      RN=SIGM*RL
      U1=1.2D00*RN
      U2=0.1D00*SIGM
      U3=0.13333333333333333D00*SIGM*XL
      U4=0.03333333333333333D00*SIGM*XL
      RP(1 ,1)=RN
      RP(1,7)=-RN
      RP(7,7)=RN
      RP(2,2)=U1
      RP(2,6)=U2
      RP(2,8)=-U1
      RP(2,12)=U2
      RP(3,3)=U1
      RP(3,5)=-U2
      RP(3,9)=-U1
      RP(3,11)=-U2
      RP(5,5)=U3
      RP(5,9)=U2
      RP(5,11)=-U4
      RP(6,6)=U3
      RP(6,8)=-U2
      RP(6,12)=-U4
      RP(8,8)=U1
      RP(8,12)=-U2
      RP(9,9)=U1
      RP(9,11)=U2
      RP(11,11)=U3
      RP(12,12)=U3
C
C  REMPLISSAGE DE LA MATRICE
C
      DO 122 I=1,12
      DO 9 J=I,12
      RP(J,I)=RP(I,J)
    9 CONTINUE
  122 CONTINUE
      RETURN
      END