1 : C TCOQ4A SOURCE CHAT 05/01/13 03:34:56 5004 2 : SUBROUTINE TCOQ4A(POINTS,QSIGAU,ETAGAU,DZEGAU,POIGAU,SHPTOT) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : ************************************************************************ 6 : * 7 : * T C O Q 4 A 8 : * ----------- 9 : * 10 : * FONCTION: ON ACQUIERT: 11 : * --------- - LES COORDONNEES INTRINSEQUES 12 : * - LES FONCTIONS DE FORME ET LEURS DERIVEES 13 : * POUR L'ELEMENT DE THERMIQUE COQ4 14 : * 15 : * MODULES UTILISES: 16 : * ----------------- 17 : * 18 : -INC CCOPTIO 19 : * 20 : * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN) 21 : * ----------- 22 : * 23 : * POINTS (E) CHAINE DE CARACTERES: 24 : * 'NOEUDS' SI ON VEUT LES VALEURS AUX NOEUDS 25 : * 'GAUSS' POUR LES VALEURS AUX POINTS DE GAUSS. 26 : * QSIGAU (S) | 27 : * ETAGAU (S) | COORDONNEES INTRINSEQUES. 28 : * POIGAU (S) POIDS D'INTEGRATION,SI ON A 'GAUSS'. 29 : * DZEGAU (S) | 30 : * SHPTOT (S) FONCTIONS DE FORME ET LEURS DERIVEES. 31 : * 32 : REAL*8 QSIGAU(8),ETAGAU(8),POIGAU(8) 33 : REAL*8 SHPTOT(6,4,8),DZEGAU(8) 34 : CHARACTER*(*) POINTS 35 : * 36 : * CONSTANTES: 37 : * ----------- 38 : * 39 : * 40 : PARAMETER ( X=0.577350269189626D0 ) 41 : PARAMETER ( O1=1.D0 ) 42 : PARAMETER ( O25=0.25D0 ) 43 : * 44 : * 45 : * AUTEUR, DATE DE CREATION: 46 : * ------------------------- 47 : * 48 : * J-P.SAUSSAIS OCT.77,REPRIS PAR DENIS ROBERT,LE 23 OCTOBRE 1987. 49 : * ADAPTE POUR LES COQUES P DOWLATYARI MAI 92 50 : * 51 : * LANGAGE: 52 : * -------- 53 : * 54 : * FORTRAN77 55 : * 56 : ************************************************************************ 57 : * 58 : IF (POINTS(1:5).EQ.'NOEUD') GOTO 66 59 : IF (POINTS(1:5).NE.'GAUSS') GOTO 100 60 : * 61 : * COORDONNEES DES 4 POINTS DE GAUSS 62 : * 63 : NBPGAU=8 64 : NBNO=4 65 : QSIGAU(3)=X 66 : QSIGAU(2)=X 67 : ETAGAU(3)=X 68 : ETAGAU(4)=X 69 : QSIGAU(4)=-X 70 : QSIGAU(1)=-X 71 : ETAGAU(2)=-X 72 : ETAGAU(1)=-X 73 : * 74 : DO 500 I=1,4 75 : QSIGAU(I+4)=QSIGAU(I) 76 : ETAGAU(I+4)=ETAGAU(I) 77 : DZEGAU(I) =-X 78 : DZEGAU(I+4)= X 79 : 500 CONTINUE 80 : POIGAU(1)=O1 81 : POIGAU(2)=O1 82 : POIGAU(3)=O1 83 : POIGAU(4)=O1 84 : DO 600 I=1,4 85 : POIGAU(I+4)=POIGAU(I) 86 : 600 CONTINUE 87 : GOTO 2 88 : * 89 : * COORDONNEES DES NOEUDS 90 : * 91 : 66 CONTINUE 92 : NBNO=4 93 : NBPGAU=NBNO 94 : QSIGAU(3)=O1 95 : ETAGAU(3)=O1 96 : QSIGAU(4)=-O1 97 : ETAGAU(4)=O1 98 : QSIGAU(1)=-O1 99 : ETAGAU(1)=-O1 100 : QSIGAU(2)=O1 101 : ETAGAU(2)=-O1 102 : * 103 : * FONCTIONS DE FORME ET LEURS DERIVEES 104 : * 105 : 2 DO 1 I=1,NBPGAU 106 : XI=QSIGAU(I) 107 : YI=ETAGAU(I) 108 : SHPTOT(1,3,I)=O25*(O1+XI)*(O1+YI) 109 : SHPTOT(1,4,I)=O25*(O1-XI)*(O1+YI) 110 : SHPTOT(1,1,I)=O25*(O1-XI)*(O1-YI) 111 : SHPTOT(1,2,I)=O25*(O1+XI)*(O1-YI) 112 : C 113 : SHPTOT(2,3,I)=O25*(O1+YI) 114 : SHPTOT(2,4,I)=-O25*(O1+YI) 115 : SHPTOT(2,1,I)=-O25*(O1-YI) 116 : SHPTOT(2,2,I)=O25*(O1-YI) 117 : C 118 : SHPTOT(3,3,I)=O25*(O1+XI) 119 : SHPTOT(3,4,I)=O25*(O1-XI) 120 : SHPTOT(3,1,I)=-O25*(O1-XI) 121 : SHPTOT(3,2,I)=-O25*(O1+XI) 122 : 1 CONTINUE 123 : RETURN 124 : 100 CONTINUE 124 : CALL ERREUR(19) 125 : * 126 : END 127 :
© Cast3M 2003 - All rights reserved.
Disclaimer