Télécharger la plaquette Cast3M


fr en

Télécharger tconv4.eso

Retour à la liste

Numérotation des lignes : 

tconv4

  1. C TCONV4    SOURCE    CHAT      05/01/13    03:34:29     5004
  2.  
  3. C=======================================================================
  4. C=                            T C O N V 4                              =
  5. C=                            -----------                              =
  6. C=                                                                     =
  7. C=  Fonction :                                                         =
  8. C=  ----------                                                         =
  9. C=   Calcul du jacobien et des derivees des fonctions de forme dans la =
  10. C=   geometrie reelle pour les elements SEG2, SEG3, LIA3 et LIA4.      =
  11. C=                                                                     =
  12. C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  13. C=  ------------                                                       =
  14. C=   XEL    (E)   COORDONNEES des noeuds de l'element (repere GLOBAL)  =
  15. C=   SHP   (E/S)  DERIVEES des fonctions de FORME de l'element         =
  16. C=   IDIM   (E)   DIMENSION du probleme                                =
  17. C=   NBNN   (E)   Nombre de NOEUDS de l'element fini considere         =
  18. C=   DJAC   (S)   JACOBIEN                                             =
  19. C=                                                                     =
  20. C=  Denis ROBERT, le 3 mai 1988.                                       =
  21. C=======================================================================
  22.  
  23.       SUBROUTINE TCONV4 (XEL,SHP,IDIM,NBNN,DJAC)
  24.  
  25.       IMPLICIT INTEGER(I-N)
  26.       IMPLICIT REAL*8 (A-H,O-Z)
  27.  
  28. -INC CCREEL
  29.  
  30.       DIMENSION XEL(3,*),SHP(6,*)
  31.  
  32. C ===========================
  33. C  1 - Cas de la DIMENSION 2
  34. C ===========================
  35.       IF (IDIM.EQ.2) THEN
  36.         A1=XZERO
  37.         A2=XZERO
  38.         DO i=1,NBNN
  39.           A1=A1+SHP(2,i)*XEL(1,i)
  40.           A2=A2+SHP(2,i)*XEL(2,i)
  41.         ENDDO
  42.         DJAC=SQRT(A1*A1+A2*A2)
  43.         IF (NBNN.EQ.2) THEN
  44.           A1=XEL(2,2)-XEL(2,1)
  45.           A2=XEL(1,2)-XEL(1,1)
  46.           DLON=SQRT(A1*A1+A2*A2)
  47.           COSANG=A2/DLON
  48.           SINANG=A1/DLON
  49. C** Ces 4 lignes restent a verifier
  50.           SHP(3,1)=SHP(2,1)*SINANG
  51.           SHP(2,1)=SHP(2,1)*COSANG
  52.           SHP(3,2)=SHP(2,2)*SINANG
  53.           SHP(2,2)=SHP(2,2)*COSANG
  54.         ENDIF
  55. C ===========================
  56. C  2 - Cas de la DIMENSION 3
  57. C ===========================
  58.       ELSE IF (IDIM.EQ.3) THEN
  59.         A1=XZERO
  60.         A2=XZERO
  61.         A3=XZERO
  62.         DO i=1,NBNN
  63.           A1=A1+SHP(2,i)*XEL(1,i)
  64.           A2=A2+SHP(2,i)*XEL(2,i)
  65.           A3=A3+SHP(2,i)*XEL(3,i)
  66.         ENDDO
  67.         DJAC=SQRT(A1*A1+A2*A2+A3*A3)
  68.       ENDIF
  69.  
  70.       RETURN
  71.       END
  72.  
  73.  
  74.  
  75.  

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