1 : C COQ3ST SOURCE CHAT 05/01/12 22:21:38 5004 2 : SUBROUTINE COQ3ST(XXE,UUE,CONT,DDHOOK) 3 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 4 : C 5 : C ELEMENT COQ3 CALCUL DES EFFORTS AU CENTRE DE GRAVITE. 6 : C 7 : CCCCCCCCCCC ENTREE 8 : C XXE(3,3) COORDONEES DES 3 NOEUDS DE L ELEMENT 9 : C RANGES X(1,2)=COORDONNEE X DU NOEUD 2 10 : C X(2,2)=COORDONNEE Y DU NOEUD 2 11 : C X(3,2)=COORDONNEE Z DU NOEUD 2 12 : C UUE(18) VECTEUR DEPLACEMENT DE L ELEMENT 13 : C 14 : CCCCCCCCCCC SORTIE 15 : C CONT EFFORTS AUX CENTRE DE GRAVITE 16 : C RANGES 3 EFFORTS PUIS 17 : C 3 MOMENTS DE FLEXION 18 : C S1 EST PORTE PAR LE COTE 1 DE L ELEMENT 19 : C S2 EST PORTE PAR LA DIRECTION PERPENDICULAIRE 20 : C A 1 ET DANS LE PLAN DE L ELEMENT 21 : C 22 : CCCCCCCCCCCC TRAVAIL 23 : C 24 : CCCCCCCCCCCCCCCCAUTRES TABLEUX IL FAUT EN TOUT 160 25 : C 26 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 27 : IMPLICIT INTEGER(I-N) 28 : IMPLICIT REAL*8(A-H,O-Z) 29 : DIMENSION XXE(3,*),UUE(*),CONT(*),DDHOOK(6,*) 30 : DIMENSION XGRAV(3),TW(6),CM(3),TM(3),CONT0(6) 31 : DIMENSION UE(18),XX(3),YY(3),P(3,3),XE(3),YE(3),ZE(3), 32 : 1 DT(6),DF(9) 33 : 1,SK(3,6),T(3,3),Q(3,9),SD(3),CI(3),CE(3) 34 : DATA UNTIER,DETIER/.333333333333333D0,.6666666666666667D0/ 35 : DATA QUTIER,CISOIX/1.33333333333333D0,.0833333333333333D0/ 36 : DATA UNSIX/0.1666666666666667D0/ 37 : C 38 : C 39 : C 40 : PS1=0.5 41 : RAI=1. 42 : DO 2 IP=1,18 43 : 2 UE(IP)=UUE(IP) 44 : DO 4 IP=1,3 45 : XE(IP)=XXE(1,IP) 46 : YE(IP)=XXE(2,IP) 47 : ZE(IP)=XXE(3,IP) 48 : 4 CONTINUE 48 : CALL PASSA(XE,YE,ZE, P,X13,X23,Y13,Y23) 49 : XGRAV(1)=(XE(1)+XE(2)+XE(3))*UNTIER 50 : XGRAV(2)=(YE(1)+YE(2)+YE(3))*UNTIER 51 : XGRAV(3)=(ZE(1)+ZE(2)+ZE(3))*UNTIER 52 : DO 10 IP=1,6 53 : IP1=3*(IP-1) 54 : DO 5 IQ=1,3 55 : 5 XX(IQ)=UE(IQ+IP1) 55 : CALL MULMAT(YY,P,XX,3,1,3) 56 : DO 7 IQ=1,3 57 : 7 UE(IQ+IP1)=YY(IQ) 58 : 10 CONTINUE 59 : DO 8 IP=1,3 60 : IP1=IP-1 61 : IP2=2*IP1 62 : IP3=3*IP1 63 : IP6=6*IP1 64 : DO 9 IQ=1,2 65 : 9 DT(IQ+IP2)=UE(IQ+IP6) 66 : DO 11 IQ=3,5 67 : IQ2=IQ-2 68 : 11 DF(IQ2+IP3)=UE(IQ+IP6) 69 : 8 CONTINUE 70 : S=1./(X13*Y23-X23*Y13) 71 : A1=Y13 72 : A2=-X13 73 : A3=Y23 74 : A4=-X23 75 : SK(1,1)=Y23 76 : SK(1,2)=0. 77 : SK(1,3)=-Y13 78 : SK(1,4)=0. 79 : SK(1,5)=Y13-Y23 80 : SK(1,6)=0. 81 : SK(2,1)=0. 82 : SK(2,2)=-X23 83 : SK(2,3)=0. 84 : SK(2,4)=X13 85 : SK(2,5)=0. 86 : SK(2,6)=X23-X13 87 : SK(3,1)=-X23 88 : SK(3,3)=X13 89 : SK(3,5)=X23-X13 90 : SK(3,2)=Y23 91 : SK(3,4)=-Y13 92 : SK(3,6)=Y13-Y23 93 : T(1,1)=Y23*Y23 94 : T(1,2)=Y13*Y13 95 : T(1,3)=-2.*(Y13*Y23) 96 : T(2,1)=X23*X23 97 : T(2,2)=X13*X13 98 : T(2,3)=-2.*(X13*X23) 99 : T(3,1)=-2.*X23*Y23 100 : T(3,2)=-2.*X13*Y13 101 : T(3,3)=2.*(X13*Y23+X23*Y13) 102 : Q(1,1)=DETIER 103 : Q(1,2)=UNTIER*(2.*A1+A3) 104 : Q(1,3)=UNTIER*(2.*A2+A4) 105 : Q(1,4)=-QUTIER 106 : Q(1,5)=UNTIER*(-A1+2.*A3) 107 : Q(1,6)=UNTIER*(-A2+2.*A4) 108 : Q(1,7)=DETIER 109 : Q(1,8)=-A1+UNTIER*A3 110 : Q(1,9)=-A2+UNTIER*A4 111 : Q(2,1)=-QUTIER 112 : Q(2,2)=UNTIER*(2.*A1-A3) 113 : Q(2,3)=UNTIER*(2.*A2-A4) 114 : Q(2,4)=DETIER 115 : Q(2,5)=UNTIER*(A1+2.*A3) 116 : Q(2,6)=UNTIER*(A2+2.*A4) 117 : Q(2,7)=DETIER 118 : Q(2,8)=-A3+UNTIER*A1 119 : Q(2,9)=-A4+UNTIER*A2 120 : Q(3,1)=-DETIER 121 : Q(3,2)=UNTIER*A1+0.5*A3 122 : Q(3,3)=UNTIER*A2+0.5*A4 123 : Q(3,4)=-DETIER 124 : Q(3,5)=0.5*A1+UNTIER*A3 125 : Q(3,6)=0.5*A2+UNTIER*A4 126 : Q(3,7)=QUTIER 127 : Q(3,8)=-UNSIX*(A1+A3) 128 : Q(3,9)=-UNSIX*(A2+A4) 128 : CALL MULMAT(CM,SK,DT,3,1,6) 128 : CALL MULMAT(SD,Q,DF,3,1,9) 128 : CALL MULMAT(TM,T,SD,3,1,3) 129 : CPLA=S*S 130 : CMEM=RAI*S 131 : DO 306 IP=1,3 132 : CONT0(IP+3)=-TM(IP)*CPLA 133 : CONT0(IP )= CM(IP)*CMEM 134 : 306 CONTINUE 134 : CALL MULMAT(CONT,DDHOOK,CONT0,6,1,6) 135 : RETURN 136 : END 137 :
© Cast3M 2003 - All rights reserved.
Disclaimer