1 : C ELLP31 SOURCE CHAT 05/01/12 23:35:39 5004 2 : SUBROUTINE ELLP31 (R,INP,ZS,ZALFAX,ZALFAR,ZALFAY,ZALFAZ, 3 : * ZG01,ZG11,ZG02,ZG12, 4 : * ZG03,ZG13,ZG23,ZG33,ZG43,ZZ3, 5 : * ZG04,ZG14,ZG24,ZG34,ZG44,ZZ4) 6 : C 7 : IMPLICIT INTEGER(I-N) 8 : IMPLICIT REAL*8 (A-H,O-Y) 9 : IMPLICIT COMPLEX*16 (Z) 10 : C 11 : C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 12 : C 13 : C OPERATEUR ELFE LAPLACE POUTRE 14 : C 15 : C CALCUL POUR LA POUTRE N› INP DES FONCTIONS DE GREEN CORRESPONDAN- 16 : C TES EN - TC (N› 1) 17 : C - TORSION (N› 2) 18 : C - FLEXION / OY (N› 3) ) 19 : C - FLEXION / OZ (N› 4) ) --> MODELE DE TIMOSHENKO 20 : C AU POINT DE COORDONNEES LOCALES R. 21 : C 22 : C PARAMETRES : 23 : C 24 : C CARACT : TABLEAU DES CARACTERISTIQUE DES POUTRES (12 ,NP) 25 : C R : POINT DE CALCUL DE G ( 0 < R < L) 26 : C INP : NUMERO DE LA POUTRE 27 : C ZS : VALEUR DE S = S0 + I W 28 : C 29 : C SORTIES : 30 : C 31 : C ZG'N''M' : VALEUR DE LA DERIVEE 'N'IEME DELA FONCTION DE GREEN 32 : C NUMERO M. 33 : C 34 : C 35 : C AUTEUR : SAINT-DIZIER 36 : C DATE : 04 JANVIER 1990 (VERSION DU 03 AVRIL 1990) 37 : C 38 : C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 39 : C 40 : C 41 : C -- CALCUL DE GREEN 1 : TRACTION COMPRESSION 42 : C -------------------- 43 : ZG01 = .5D0*EXP(-ZALFAX*R) 44 : ZG11 =-.5D0*EXP(-ZALFAX*R) 45 : C 46 : C -- CALCUL DE GREEN 2 : TORSION 47 : C ------- 48 : ZG02 = .5D0*EXP(-ZALFAR*R) 49 : ZG12 =-.5D0*EXP(-ZALFAR*R) 50 : C 51 : C -- CALCUL DE GREEN 3 : FLEXION / OY 52 : C ------------ 53 : ZDELTA= ZZ3*ZZ3 - 4.D0 54 : ZDEL = SQRT(ZDELTA) 55 : C 56 : Z1 = (ZZ3+ZDEL)/2.D0 57 : Z2 = (ZZ3-ZDEL)/2.D0 58 : ZL1 = SQRT(Z1) 59 : ZL2 = SQRT(Z2) 60 : PR1 = ZL1*ZALFAY 61 : PR2 = ZL2*ZALFAY 62 : C 63 : C -- SELECTION DES PARTIES REELLES POSITIVES 64 : C 65 : IF (PR1.LT.0.) ZL1 = -ZL1 66 : IF (PR2.LT.0.) ZL2 = -ZL2 67 : C 68 : ZZA = 1.D0/(2.D0*ZL1*(ZL2*ZL2-ZL1*ZL1)) 69 : ZZB = 1.D0/(2.D0*ZL2*(ZL1*ZL1-ZL2*ZL2)) 70 : C 71 : ZE1 = EXP(-ZL1*ZALFAY*R) 72 : ZE2 = EXP(-ZL2*ZALFAY*R) 73 : C 74 : ZG03 = ZZA*ZE1 + ZZB*ZE2 75 : ZG13 = -ZL1*ZZA*ZE1 - ZL2*ZZB*ZE2 76 : ZG23 = ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZZB*ZE2 77 : ZG33 = -ZL1*ZL1*ZL1*ZZA*ZE1 - ZL2*ZL2*ZL2*ZZB*ZE2 78 : ZG43 = ZL1*ZL1*ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZL2*ZL2*ZZB*ZE2 79 : C 80 : C -- CALCUL DE GREEN 4 : FLEXION / OZ 81 : C ------------ 82 : ZDELTA= ZZ4*ZZ4 - 4.D0 83 : ZDEL = SQRT(ZDELTA) 84 : C 85 : Z1 = (ZZ4+ZDEL)/2.D0 86 : Z2 = (ZZ4-ZDEL)/2.D0 87 : ZL1 = SQRT(Z1) 88 : ZL2 = SQRT(Z2) 89 : PR1 = ZL1*ZALFAZ 90 : PR2 = ZL2*ZALFAZ 91 : C 92 : C -- SELECTION DES PARTIES REELLES POSITIVES 93 : C 94 : IF (PR1.LT.0.) ZL1 = -ZL1 95 : IF (PR2.LT.0.) ZL2 = -ZL2 96 : C 97 : ZZA = 1.D0/(2.D0*ZL1*(ZL2*ZL2-ZL1*ZL1)) 98 : ZZB = 1.D0/(2.D0*ZL2*(ZL1*ZL1-ZL2*ZL2)) 99 : C 100 : ZE1 = EXP(-ZL1*ZALFAZ*R) 101 : ZE2 = EXP(-ZL2*ZALFAZ*R) 102 : C 103 : ZG04 = ZZA*ZE1 + ZZB*ZE2 104 : ZG14 = -ZL1*ZZA*ZE1 - ZL2*ZZB*ZE2 105 : ZG24 = ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZZB*ZE2 106 : ZG34 = -ZL1*ZL1*ZL1*ZZA*ZE1 - ZL2*ZL2*ZL2*ZZB*ZE2 107 : ZG44 = ZL1*ZL1*ZL1*ZL1*ZZA*ZE1 + ZL2*ZL2*ZL2*ZL2*ZZB*ZE2 108 : C 109 : END 110 : 111 :
© Cast3M 2003 - All rights reserved.
Disclaimer