Télécharger hs456.eso

Retour à la liste

Numérotation des lignes :

  1. C HS456 SOURCE CHAT 05/01/13 00:26:29 5004
  2. SUBROUTINE HS456(XE,NSTRS,DDHOMU,HS4,HS5,HS6)
  3. *
  4. * CALCUL DES FONCTIONS HS4,HS5 ET HS6
  5. * POUR L'ELEMENT COQUE DST
  6. *
  7. ************************************************************************
  8. *
  9. * ENTREES :
  10. *
  11. * XE : COORDONNEES DES NOEUDS
  12. * NSTRS :NOMBRE DE COMPOSANTES DE CONTRAINTES
  13. * DDHOMU :MATRICE DE HOOKE
  14. *
  15. * SORTIES :
  16. *
  17. * HS4,HS5,HS6 : LES FONCTIONS QUI RELIE LES TERMES (BETA)SK A {U}
  18. *
  19. * AUTEUR : X.Z. SUO FEVRIER 96
  20. ************************************************************************
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. DIMENSION XE(3,*),DDHOMU(NSTRS,*),HS4(*),HS5(*),HS6(*)
  24. DIMENSION SX(3),SY(3),RL(3),EXX(3),EYY(3),HFX(9),HFY(9)
  25. DIMENSION ASS(6),BSS(6),AXX(3),BXX(3),AYY(3),BYY(3)
  26. *
  27. DO 10 K=4,6
  28. IF(K.EQ.4)THEN
  29. IJ=1
  30. I=2
  31. J=3
  32. ELSEIF(K.EQ.5)THEN
  33. IJ=2
  34. I=3
  35. J=1
  36. ELSE
  37. IJ=3
  38. I=1
  39. J=2
  40. ENDIF
  41. SX(IJ)=XE(1,I)-XE(1,J)
  42. SY(IJ)=XE(2,I)-XE(2,J)
  43. RL(IJ)=SQRT(SX(IJ)*SX(IJ)+SY(IJ)*SY(IJ))
  44. EXX(IJ)=-SX(IJ)/RL(IJ)
  45. EYY(IJ)=-SY(IJ)/RL(IJ)
  46. 10 CONTINUE
  47. AIR=ABS(0.5D0*(SX(1)*SY(2)-SX(2)*SY(1)))
  48. SOMM=SX(IJ)+SY(IJ)+RL(IJ)+EXX(IJ)+EYY(IJ)
  49. *
  50. CALL CISAF(XE,EXX,EYY,DDHOMU,AIR,ASS,AXX,
  51. 1 AYY,BSS,BXX,BYY)
  52. *
  53. CALL CISAR(ASS,BSS,AXX,AYY,BXX,BYY,RL,EXX,EYY,
  54. 1 HS4,HS5,HS6,HFX,HFY)
  55. *
  56. RETURN
  57. END
  58.  
  59.  
  60.  

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