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 :
© Cast3M 2003 - All rights reserved.
Disclaimer