Télécharger la plaquette Cast3M


fr en

Télécharger cq4ksg.eso

Retour à la liste

Numérotation des lignes : 

cq4ksg

  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