Télécharger la plaquette Cast3M


fr en

Télécharger podefo.eso

Retour à la liste

Numérotation des lignes : 

podefo

  1. C PODEFO    SOURCE    CHAT      11/03/22    21:15:13     6910
  2.       SUBROUTINE PODEFO(V,DEFO,XL,IDEUX,carac,you,xnu)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DES DEFORMATIONS DE LA POUTRE DROITE
  6. C
  7. C  ENTREE
  8. C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  9. C    XL      LONGUEUR
  10. C    IDEUX   INDICE POUR DEFORMATIONS DU 2-EME ORDRE
  11. C
  12. C  SORTIE
  13. C    DEFO    DEFORMATIONS AUX 2 NOEUDS    6 PAR NOEUD
  14. C            ( PAS D'EFFET D'EFFORT TRANCHANT || )
  15. C
  16. C-----------------------------------------------------------------------
  17.       IMPLICIT INTEGER(I-N)
  18.       IMPLICIT REAL*8(A-H,O-Z)
  19.       DIMENSION DEFO(*),V(*),carac(*)
  20. C
  21.       RL=1.D00/XL
  22.       RL2=RL*RL
  23.           X2=0.
  24.           X3=0.
  25.           xiny=carac(2)
  26.           xinz=carac(3)
  27.           sredy=carac(5)
  28.           sredz=carac(6)
  29.           cisail= you/(2.*(1+xnu))
  30. * of
  31. *          if( sredy.ne.0.) X2=You/cisail*rl2*xiny/sredy
  32. *          if( sredy.ne.0.) X2=You/cisail*rl2*xinz/sredz
  33.           if( sredy.ne.0.) X3=You/cisail*rl2*xinz/sredy
  34.           if( sredz.ne.0.) X2=You/cisail*rl2*xiny/sredz
  35. * of
  36.              X13=1./(1.+12.*X3)
  37.              X33=(1.+3.*X3)*X13
  38.              X63=(1.-6.*X3)*X13
  39.              X12=1./(1.+12.*X2)
  40.              X32=(1.+3.*X2)*X12
  41.              X62=(1.-6.*X2)*X12
  42. C
  43. C  DEFORMATIONS GENERALISEES
  44. C
  45.       EPS  =(V(7)-V(1))*RL
  46.       XSIX =(V(10)-V(4))*RL
  47.       XSIY1=
  48.      .RL2*( 6.*V(3)*X12-4.*XL*X32*V(5)-6.*V(9)*X12-2.*XL*X62*V(11))
  49.       XSIY2=
  50.      .RL2*(-6.*V(3)*X12+2.*XL*X62*V(5)+6.*V(9)*X12+4.*XL*X32*V(11))
  51.       XSIZ1=
  52.      .RL2*(-6.*V(2)*X13-4.*XL*X33*V(6)+6.*V(8)*X13-2.*XL*X63*V(12))
  53.       XSIZ2=
  54.      .RL2*( 6.*V(2)*X13+2.*XL*X63*V(6)-6.*V(8)*X13+4.*XL*X33*V(12))
  55. C
  56. C  REMPLISSAGE DU TABLEAU DEFO
  57. C
  58.       DEFO(1)=EPS
  59.       DEFO(2)=0.D0
  60.       DEFO(3)=0.D0
  61.       DEFO(4)=XSIX
  62.       DEFO(5)=XSIY1
  63.       DEFO(6)=XSIZ1
  64.       DEFO(7)=EPS
  65.       DEFO(8)=0.D0
  66.       DEFO(9)=0.D0
  67.       DEFO(10)=XSIX
  68.       DEFO(11)=XSIY2
  69.       DEFO(12)=XSIZ2
  70.       IF(IDEUX.EQ.0) RETURN
  71. C
  72. C  TERMES DU SECOND ORDRE
  73. C
  74.       AUX1=0.5D0*(EPS**2+V(5)**2+V(6)**2)
  75.       AUX2=0.5D0*(EPS**2+V(11)**2+V(12)**2)
  76.       DEFO(1)=DEFO(1)+0.5D0*(AUX1+AUX2)
  77.       DEFO(7)=DEFO(7)+0.5D0*(AUX1+AUX2)
  78. C*OF?    DEFO(1)=DEFO(1)+AUX1
  79. C*OF?    DEFO(7)=DEFO(7)+AUX2
  80.       RETURN
  81.       END
  82.  
  83.  
  84.  
  85.  

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