Télécharger cq4ksg.eso

Retour à la liste

Numérotation des lignes :

  1. C CQ4KSG SOURCE CHAT 05/01/12 22:27:01 5004
  2. SUBROUTINE CQ4KSG(DJAC,XSTRS,SHP,REL)
  3. C=====================================================================
  4. C MATICE KSIGMA POUR L ELEMENT COQ4
  5. C
  6. C ENTREES
  7. C DJAC = JACOBIEN
  8. C XSTRS(6) = CONTRAINTES AU POINT DE GAUSS
  9. * SHP = Fonctions de forme et d{riv{es dans l'espace
  10. * g{om{trique, au centre de l'{l{ment.
  11. * SHP(2,... = d{riv{e par rapport @ "X" local,
  12. * SHP(3,... = d{riv{e par rapport @ "Y" local.
  13. *
  14. C SORTIES
  15. C REL (24,24) = MATRICE K SIGMA
  16. C
  17. C DELHOMMEAU
  18. * P.M. 17/06/91: reprise pour nouveau COQ4, avec utilisation de
  19. * "SHP" au lieu de "BGENE".
  20. *
  21. C=======================================================================
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8(A-H,O-Z)
  24. DIMENSION REL(24,*),XSTRS(*),SHP(6,*)
  25. C
  26. C
  27. SIGXX=XSTRS(1)
  28. SIGYY=XSTRS(2)
  29. SIGXY=XSTRS(3)
  30. C
  31. DO 100 IA=1,4
  32. JJ=(IA-1)*6+1
  33. DO 200 IB=1,4
  34. II=(IB-1)*6+1
  35. CC= SHP(2,IA)*SHP(2,IB)*SIGXX + SHP(3,IA)*SHP(3,IB)*SIGYY
  36. DD= (SHP(2,IA)*SHP(3,IB) + SHP(2,IB)*SHP(3,IA)) * SIGXY
  37. S =(CC+DD)*DJAC*4.D0
  38. REL(II,JJ) = S
  39. REL(II+1,JJ+1)= S
  40. REL(II+2,JJ+2)= S
  41. 200 CONTINUE
  42. 100 CONTINUE
  43. C
  44. RETURN
  45. END
  46.  
  47.  
  48.  

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