1 : C INTCN2 SOURCE PV 09/03/12 21:25:20 6325 2 : SUBROUTINE INTCN2(IPVANO,IPGEOM,IPINTE,IPVAEQ) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : ************************************************************************ 6 : * 7 : * I N T C N 2 8 : * ----------- 9 : * 10 : * FONCTION: 11 : * --------- 12 : * INTEGRATION NUMERIQUE DANS UN DOMAINE BIDIMENSIONNEL DU PRODUIT: 13 : * COEF. * TRANSPOSEE( N ) 14 : * COEF. : GRANDEUR PHYSIQUE REPRESENTEE PAR UN CHAMELEM 15 : * N : FONCTIONS DE FORME DE L'ELEMENT MAILLANT LE DOMAINE 16 : * CONSIDERE 17 : * 18 : * MODULES UTILISES: 19 : * ----------------- 20 : * 21 : -INC CCOPTIO 22 : -INC CCREEL 23 : -INC SMCHAML 24 : -INC SMELEME 25 : -INC SMINTE 26 : -INC SMCOORD 27 : * 28 : * PARAMETRES: (E)=ENTREE (S)=SORTIE (+ = CONTENU DANS UN COMMUN) 29 : * ----------- 30 : * 31 : * IPVANO (E) POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES 32 : * VALEURS NODALES DE NOTRE COEFFICIENT 33 : * IPGEOM (E) POINTEUR SUR UN OBJET MAILLAGE ELEMENTAIRE 34 : * DU DOMAINE D'INTEGRATION 35 : * IPINTE (E) POINTEUR SUR UN SEGMENT MINTE CONTENANT LES 36 : * CARACTERISTIQUES D'INTEGRATION 37 : * +IDIM (E) VOIR CCOPTIO 38 : * +XPI (E) VOIR CCREEL 39 : * IPVAEQ (S) POINTEUR SUR UN SEGMENT MELVAL CONTENANT LES VALEURS 40 : * NODALES EQUIVALENTES 41 : * 42 : * VARIABLES: 43 : * ---------- 44 : * 45 : * XE(3,NBPTEL) = COORDONNEES DES ELEMENTS DANS LE REPERE GLOBAL 46 : * A ET S = TABLEAUX DE TRAVAIL 47 : * 48 : REAL*8 S(2,3) 49 : SEGMENT,MMAT1 50 : REAL*8 XE(3,NBPTEL),A(NBPTEL,NBPTEL),SHP(6,NBPTEL) 51 : ENDSEGMENT 52 : * 53 : * CONSTANTES: 54 : * ----------- 55 : * 56 : PARAMETER ( O0=0.D0 ) 57 : PARAMETER ( O1=1.D0 ) 58 : PARAMETER ( O2=2.D0 ) 59 : * 60 : * REMARQUES: 61 : * ---------- 62 : * 63 : * L'UTILISATION DE CE S-P PRESUPPOSE UN PRE ET POST-TRAITEMENT 64 : * DES SEGMENTS MELVAL PASSES EN TANT QUE PARAMETRES 65 : * 66 : * AUTEUR, DATE DE CREATION: 67 : * ------------------------- 68 : * 69 : * DENIS ROBERT,LE 15 AVRIL 1988. 70 : * 71 : * LANGAGE: 72 : * -------- 73 : * 74 : * ESOPE + FORTRAN77 75 : * 76 : ************************************************************************ 77 : * 78 : * ON RECUPERE LES VALEURS DU COEFFICIENT 79 : * 80 : MELVA1=IPVANO 81 : SEGACT,MELVA1 82 : NBPTE1=MELVA1.VELCHE(/1) 83 : NEL1=MELVA1.VELCHE(/2) 84 : * 85 : * ON RECUPERE LES CARACTERISTIQUES D'INTEGRATION 86 : * 87 : MINTE=IPINTE 88 : SEGACT,MINTE 89 : NBPGAU=POIGAU(/1) 90 : * 91 : * ON RECUPERE UN DES MAILLAGES ELEMENTAIRES DE L'ENVELOPPE 92 : * 93 : MELEME=IPGEOM 94 : SEGACT,MELEME 95 : NBPTEL=NUM(/1) 96 : NEL=NUM(/2) 97 : * 98 : * INITIALISATION DU MELVAL QUI CONTIENDRA LES VALEURS EQUIVALENTES 99 : * 100 : N1PTEL=NBPTEL 101 : N1EL=NEL 102 : N2PTEL=0 103 : N2EL=0 104 : SEGINI,MELVAL 105 : IPVAEQ=MELVAL 106 : SEGINI,MMAT1 107 : * 108 : * BOUCLE SUR LES ELEMENTS 109 : * 110 : DO 10 IEL=1,NEL 111 : * 112 : * ON CHERCHE LES COORDONNEES DES NOEUDS DANS LE REPERE GLOBAL 113 : * 113 : CALL DOXE(XCOOR,IDIM,NBPTEL,NUM,IEL,XE) 113 : CALL ZERO(A,NBPTEL,NBPTEL) 114 : * 115 : * BOUCLE SUR LES POINTS DE GAUSS 116 : * 117 : DO 20 IGAU=1,NBPGAU 118 : * 119 : * CALCUL DU JACOBIEN AU POINT DE GAUSS CONSIDERE 120 : * 121 : DO 30 IP = 1,2 122 : IP1 = IP + 1 123 : DO 30 IQ = 1,3 124 : S(IP,IQ) = XZERO 125 : DO 40 IR = 1,NBPTEL 126 : S(IP,IQ)=S(IP,IQ)+SHPTOT(IP1,IR,IGAU)*XE(IQ,IR) 127 : 40 CONTINUE 128 : * END DO 129 : 30 CONTINUE 130 : * END DO 131 : S1 = S(1,2)*S(2,3)-S(1,3)*S(2,2) 132 : S2 = S(1,3)*S(2,1)-S(1,1)*S(2,3) 133 : S3 = S(1,1)*S(2,2)-S(1,2)*S(2,1) 134 : DJAC = SQRT ( S1*S1 + S2*S2 + S3*S3 ) 135 : IF (IFOMOD.EQ.0) THEN 136 : * 137 : * CAS DES ELEMENTS AXISYMETRIQUES 138 : * 139 : DO 41 NP=1,NBPTEL 140 : SHP(1,NP)=SHPTOT(1,NP,IGAU) 141 : SHP(2,NP)=SHPTOT(2,NP,IGAU) 142 : SHP(3,NP)=SHPTOT(3,NP,IGAU) 143 : 41 CONTINUE 144 : * END DO 144 : CALL DISTRR(XE,SHP,NBPTEL,RR) 145 : DJAC=O2*XPI*RR*DJAC 146 : ENDIF 147 : DO 50 INOE=1,NBPTEL 148 : DO 60 INO2=1,NBPTEL 149 : A(INOE,INO2)=A(INOE,INO2)+SHPTOT(1,INOE,IGAU)* 150 : & SHPTOT(1,INO2,IGAU)*POIGAU(IGAU)*DJAC 151 : 60 CONTINUE 152 : * END DO 153 : 50 CONTINUE 154 : * END DO 155 : 20 CONTINUE 156 : * END DO 157 : IEMIN=MIN(NEL1,IEL) 158 : DO 70 INOE=1,NBPTEL 159 : DO 80 INO2=1,NBPTEL 160 : INMIN=MIN(NBPTE1,INO2) 161 : VELCHE(INOE,IEL)=VELCHE(INOE,IEL)+MELVA1.VELCHE(INMIN,IEMIN)* 162 : & A(INO2,INOE) 163 : 80 CONTINUE 164 : * END DO 165 : 70 CONTINUE 166 : * END DO 167 : 10 CONTINUE 168 : * END DO 169 : * 170 : SEGSUP,MMAT1 171 : SEGDES,MELEME,MELVA1,MINTE 172 : SEGDES,MELVAL 173 : * 174 : END 175 : 176 : 177 : 178 : 179 :
© Cast3M 2003 - All rights reserved.
Disclaimer