Télécharger epsig0.eso

Retour à la liste

Numérotation des lignes :

  1. C EPSIG0 SOURCE CHAT 05/01/12 23:42:30 5004
  2. SUBROUTINE EPSIG0(S,EPSI,MFR,IFOUR,YUNG,XNU,CARAC,NSTRS)
  3.  
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. DIMENSION S(*),EPSI(*),CARAC(*)
  7.  
  8. DO I=1,NSTRS
  9. EPSI(I)=0.D0
  10. ENDDO
  11.  
  12. AUX1=1.D0/YUNG
  13. AUX2=2.D0*(1.D0+XNU)
  14. AUX3=AUX1*AUX2
  15. IF(MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  16. EP1=CARAC(1)
  17. IF(MFR.NE.5) ALFAH=CARAC(2)**2
  18. AUX4=AUX1/EP1
  19. AUX5=12.D0*AUX4/(EP1*EP1)
  20. ENDIF
  21.  
  22.  
  23. C---------ELEMENTS MASSIFS MFR=1----------------------------------------
  24.  
  25. IF(MFR.EQ.1) THEN
  26.  
  27. EPSI(1)=(S(1)-XNU*(S(2)+S(3)))*AUX1
  28. EPSI(2)=(S(2)-XNU*(S(1)+S(3)))*AUX1
  29. EPSI(3)=(S(3)-XNU*(S(1)+S(2)))*AUX1
  30. DO I=4,NSTRS
  31. EPSI(I)=AUX3*S(I)
  32. ENDDO
  33.  
  34. C---------ELEMENTS COQUES SANS ET AVEC C.T. MFR=3 ET MFR=9--------------
  35.  
  36. ELSE IF(MFR.EQ.3.OR.MFR.EQ.9) THEN
  37.  
  38. IF(IFOUR.GE.1) THEN
  39. EPSI(1)=(S(1)-XNU*S(2))*AUX4
  40. EPSI(2)=(S(2)-XNU*S(1))*AUX4
  41. EPSI(3)=AUX2*AUX4*S(3)
  42. EPSI(4)=(S(4)-XNU*S(5))*AUX5
  43. EPSI(5)=(S(5)-XNU*S(4))*AUX5
  44. EPSI(6)=AUX2*AUX5*S(6)
  45. ELSE
  46. EPSI(1)=(S(1)-XNU*S(2))*AUX4
  47. EPSI(2)=(S(2)-XNU*S(1))*AUX4
  48. EPSI(3)=(S(3)-XNU*S(4))*AUX5
  49. EPSI(4)=(S(4)-XNU*S(3))*AUX5
  50. ENDIF
  51.  
  52. C---------ELEMENTS COQUES EPAISSES MFR=5--------------------------------
  53.  
  54. ELSE IF(MFR.EQ.5) THEN
  55.  
  56. EPSI(1)=(S(1)-XNU*S(2))*AUX1
  57. EPSI(2)=(S(2)-XNU*S(1))*AUX1
  58. EPSI(3)=AUX3*S(3)
  59. EPSI(4)=AUX3*S(4)
  60. EPSI(5)=AUX3*S(5)
  61.  
  62. C---------ELEMENTS POUTRES MFR=7----------------------------------------
  63.  
  64. ELSE IF(MFR.EQ.7) THEN
  65.  
  66. IF(CARAC(4).NE.0.D0) EPSI(1)=AUX1*S(1)/CARAC(4)
  67. IF(CARAC(5).NE.0.D0) EPSI(2)=AUX3*S(2)/CARAC(5)
  68. IF(CARAC(6).NE.0.D0) EPSI(3)=AUX3*S(3)/CARAC(6)
  69. IF(CARAC(1).NE.0.D0) EPSI(4)=AUX3*S(4)/CARAC(1)
  70. IF(CARAC(2).NE.0.D0) EPSI(5)=AUX1*S(5)/CARAC(2)
  71. IF(CARAC(3).NE.0.D0) EPSI(6)=AUX1*S(6)/CARAC(3)
  72.  
  73. ENDIF
  74. END
  75.  
  76.  
  77.  
  78.  

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