1 : C ELLP11 SOURCE CHAT 05/01/12 23:35:11 5004 2 : SUBROUTINE ELLP11(COOR, CORRES, GAMA, CARACT, XCL, FLAG, 3 : * NUMERO, ZA1,ZSM, ZXX, ZSOL, ITERA, MASS, RMAS, NMAS, 4 : * SA, SB, SU, SR, SQ, SDELTA, SDELT1, SP, SP1, SCH, SCH1, 5 : * IPIVO, JPIVO, IAUX, XFRQ, NNPOI, ICHAR, NP, NP24, NP48, NNT, 6 : * NNT12, NFRQ, S0, XPI, METH, IMP) 7 : C 8 : IMPLICIT INTEGER(I-N) 9 : IMPLICIT REAL*8 (A-H,O-Y) 10 : IMPLICIT COMPLEX*16 (Z) 11 : C 12 : C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 13 : C 14 : C OPERATEUR ELFE LAPLACE POUTRE 15 : C 16 : C PROGRAMME PRINCIPAL DE FORTRAN EN SORTIE DE L'INTERFACE ESOPE 17 : C 18 : C L'ALGORITHME EST LE SUIVANT : 19 : C 20 : C 21 : C - ACQUISITION DES DONNEES PAR ARGUMENTS DE ELLP11 22 : C 23 : C - REMPLISSAGE DE LA 2EME PARTIE DE ZA1 (INDEPENDANT DE W) 23 :
Warning: preg_replace(): Compilation failed: unmatched parentheses at offset 6 in C:\xampp\htdocs\cast3m_kit_local_fr_en_2014\sources.php on line 147
24 : C 25 : C - BOUCLE SUR W 26 : C | 27 : C | S = S0 + I W 28 : C | 29 : C | - BOUCLE SUR CHAQUE POUTRE 30 : C | | 30 :
Warning: preg_replace(): Compilation failed: unmatched parentheses at offset 6 in C:\xampp\htdocs\cast3m_kit_local_fr_en_2014\sources.php on line 147
31 : C | | 32 : C | | - ASSEMBLAGE DES MATRICES 33 : C | | 34 : C | | - CALCUL DU SECOND MEMBRE 35 : C | | 35 :
Warning: preg_replace(): Compilation failed: unmatched parentheses at offset 6 in C:\xampp\htdocs\cast3m_kit_local_fr_en_2014\sources.php on line 147
36 : C | 37 : C - RETURN 38 : C 39 : C 40 : C PARAMETRES : 41 : C 42 : C COOR : TABLEAU DES COORDONNEES DES 2*NP NOEUDS FICTIFS 43 : C CORRES : TABLEAU D'ENTIER DES NUMEROS DES NOEUDS REELS 44 : C POUR CHAQUE NOEUD FICTIF 45 : C GAMA : TABLEAU DES VECTEUR DONNANT LE PLAN OXY POUR CHAQUE POUTRE 46 : C CARACT : TABLEAU DES CARACTERISTIQUE DES POUTRES (11 ,NP) 47 : C XCL : TABLEAU DES VALEURS DES CONDITIONS AUX LIMITES 48 : C FLAG : TABLEAU DE POINTEURS SUR XCL 49 : C NUMERO : TABLEAU DE NUMEROTATION GLOBALE <--> LOCALE 50 : C XFRQ : TABLEAU DES POINTS DE CALCUL EN FREQUENCE 51 : C 52 : C NFRQ : NOMBRE DE POINTS DE CALCUL EN FREQUENCE 53 : C NP : NOMBRE DE POUTRES 54 : C NNT : NOMBRE DE NOEUDS REELS 55 : C 56 : C SORTIES : 57 : C 58 : C ZA1 : TABLEAU COMPLEXE REPRESENTANT LA MATRICE DE "RIGIDITE" 59 : C ZSM : VECTEUR COMPLEXE SECOND MEMBRE 60 : C ZXX : VECTEUR COMPLEXE SOLUTION 61 : C 62 : C 63 : C AUTEUR : SAINT-DIZIER 64 : C DATE : 22 DECEMBRE 1989 (VERSION DU 03 AVRIL 1990) 65 : C 66 : C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 67 : C 68 : COMPLEX*16 ZA1(NP24,*),ZSM(*),ZXX(*),ZC1(12,24) 69 : COMPLEX*16 ZSOL(NNT12,*) 70 : REAL*8 COOR(3,*),GAMA(3,*),CARACT(12,*),XCL(12,*),XFRQ(*) 71 : REAL*8 SA(NP48,*),SB(*),SU(*),SR(*),SQ(*) 72 : REAL*8 SDELTA(*),SDELT1(*),SP(*),SP1(*) 73 : REAL*8 SCH(*),SCH1(*) 74 : INTEGER FLAG(*),CORRES(*),NUMERO(*),ITERA(*) 75 : INTEGER MASS(NNT,4) 76 : REAL*8 RMAS(NNT,4) 77 : INTEGER IPIVO(*),JPIVO(*),IAUX(*) 78 : C 79 : C ********************************************************************** 80 : C INITIALISATION 81 : C ********************************************************************** 82 : C 83 : ZI = (0.D0 , 1.D0) 84 : C 85 : C --------------------------------- COEFFICIENT DE NORMALISATION 86 : C (CF ELLP52) 87 : C ********************************************************************** 88 : C BOUCLE SUR LES FREQUENCES 89 : C ********************************************************************** 90 : C 91 : DO 30 IFRQ = 1 , NFRQ 92 : C 93 : C ********************************************************************** 94 : C REMISE A ZERO DES MATRICE POUR CHAQUE FREQUENCE 95 : C ********************************************************************** 96 : C 97 : C 98 : DO 58 I = 1 , 24*NP 99 : DO 59 J = 1 , 24*NP 100 : ZA1(I,J) = (0.D0,0.D0) 101 : 59 CONTINUE 102 : ZSM(I) = (0.D0,0.D0) 103 : 58 CONTINUE 104 : C 105 : C 106 : C ********************************************************************** 107 : C REMPLISSAGE DE LA 2EME PARTIE DE ZA1 (INDEPENDANT DE W) 108 : C + 2EME PARTIE DU SECOND MEMBRE 109 : C ********************************************************************** 110 : C 110 : CALL ELLP12(CORRES,XCL,FLAG,NUMERO,MASS,NMAS,NP,NP24,NNT,ZA1,ZSM) 111 : C 112 : C 113 : ZS = S0 + ZI * 2.D0 * XPI * XFRQ(IFRQ) 114 : C 115 : C ********************************************************************** 116 : C BOUCLE SUR LES POUTRES 117 : C ********************************************************************** 118 : C 119 : DO 40 INP = 1 , NP 120 : C 121 : C ------------------------------ CALCUL DE ZC1 MATRICE ELEMENTAIRE 122 : C 122 : CALL ELLP21 (CARACT,INP,ZS,COOR,GAMA,CORRES,ZC1) 123 : C 124 : C ------------------------------------- ASSEMBLAGE DE ZC1 DANS ZA1 125 : C 126 : DO 41 J = 1 , 24 127 : DO 42 I = 1 , 12 128 : ZA1(12*(INP-1)+I,24*(INP-1)+J) = ZC1(I,J) 129 : 42 CONTINUE 130 : 41 CONTINUE 131 : C 132 : 40 CONTINUE 133 : C 134 : C 135 : C 136 : C -------------------------- PRISE EN COMPTE DES MASSES PONCTUELLES 137 : C 137 : CALL ELLP22 (CORRES,FLAG,COOR,GAMA,ZA1,MASS,RMAS, 138 : * ZS,NP,NP24,NMAS,NNT) 139 : C 140 : C ********************************************************************** 141 : C IMPRESSION DE LA MATRICE ZA1 142 : C ********************************************************************** 143 : C 144 : C 145 : C IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 146 : C 147 : IF ((IMP.NE.0).AND.(IFRQ.EQ.1)) THEN 148 : C 149 : EPS = 1.D-12 150 : C 151 : WRITE(IMP,*)'IMPRESSION DE LA MATRICE ZA1 :' 152 : DO 10 I = 1 , 24*NP 153 : DO 11 J = 1 , 24 * NP 154 : IF (ABS(ZA1(I,J)).GT.EPS) THEN 155 : WRITE(IMP,1001) I,J,ZA1(I,J) 156 : END IF 157 : 11 CONTINUE 158 : 10 CONTINUE 159 : C 160 : C 161 : WRITE(IMP,*) 'IMPRESSION DU SECOND MEMBRE ZSM :' 162 : DO 22 I = 1 , 24*NP 163 : WRITE(IMP,*) 'ZSM(',I,')',ZSM(I) 164 : 22 CONTINUE 165 : C 166 : END IF 167 : C 168 : C IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 169 : C 170 : C ---------------------------------------------------------------------- 171 : C RESOLUTION DU SYSTEME ZA1 * ZXX = ZSM 172 : C ---------------------------------------------------------------------- 173 : C 174 : IF (IFRQ.EQ.1) THEN 174 : CALL ELLP51(ZA1,ZSM,ZXX,NP24,IPIVO,JPIVO,IAUX) 175 : NP2 = NP24 / 12 176 : C 177 : C --------------------------------------------------------------- 178 : C CALCUL DES ORDRES DE GRANDEURS DES DEPLACEMENTS ET 179 : C DES EFFORTS, POUVANT ETRE UTILES POUR UNIFORMISER 180 : C LA MATRICE DE RESOLUTION POUR LES METHODES ITERATIVES 181 : C --------------------------------------------------------------- 182 : C 183 : XUR = 0. 184 : XFM = 0. 185 : DO 111 IN = 1 , NP2 186 : DO 112 JN = 1 , 6 187 : XX = ABS(ZXX((IN-1)*12+JN)) 188 : IF (XX.GT.XUR) THEN 189 : XUR = XX 190 : END IF 191 : 112 CONTINUE 192 : DO 113 JN = 7 , 12 193 : XX = ABS(ZXX((IN-1)*12+JN)) 194 : IF (XX.GT.XFM) THEN 195 : XFM = XX 196 : END IF 197 : 113 CONTINUE 198 : 111 CONTINUE 199 : C 200 : C XUR ET XFM DONNE DES ORDRES DE GRANDEUR SUR LE DEPLACEMENT 201 : C ET SUR LES EFFORTS ; CE RESULTAT PEUT ETRE UTILE POUR LA 202 : C METHODE DE RESOLUTION ITERATIVE 203 : C 204 : ELSE IF (METH.EQ.1) THEN 205 : C 205 : CALL ELLP51(ZA1,ZSM,ZXX,NP24,IPIVO,JPIVO,IAUX) 206 : ELSE IF (METH.EQ.2) THEN 207 : C 207 : CALL ELLP53(ZA1,ZSM,ZXX,NP24,IPIVO,JPIVO,IAUX) 208 : ELSE IF (METH.EQ.3) THEN 209 : C 209 : CALL ELLP54(ZA1,ZSM,ZXX,NP24,NP48,NNPOI,ICHAR, 210 : * XUR,XFM,NIT,SA,SB,SU,SR,SDELTA,SP,SCH,SCH1,EPSY) 211 : ELSE IF (METH.EQ.4) THEN 212 : C 212 : CALL ELLP58(ZA1,ZSM,ZXX,NP24,NP48,NNPOI,ICHAR 213 : * ,NIT,SA,SB,SU,SR,SDELTA,SP,SCH,SCH1,EPSY) 214 : END IF 215 : C 216 : ITERA(IFRQ) = NIT 217 : C 218 : DO 51 I = 1 , 2*NP 219 : KK = CORRES ( I ) 220 : C 221 : DO 52 II = 1 , NNT 222 : IF (NUMERO(II).EQ.KK) THEN 223 : K = II 224 : END IF 225 : 52 CONTINUE 226 : C 227 : DO 53 J = 1 , 12 228 : ZSOL((K-1)*12+J,IFRQ) = ZXX((I-1)*12 + J) 229 : 53 CONTINUE 230 : C 231 : 51 CONTINUE 232 : C 233 : 30 CONTINUE 234 : C 235 : 1001 FORMAT('ZA1(',I3,',',I3,') =',2F12.3) 236 : C 237 : END 238 : 239 :
© Cast3M 2003 - All rights reserved.
Disclaimer