1 : C BSIGCO SOURCE CHAT 05/01/12 21:42:54 5004 2 : SUBROUTINE BSIGCO(XXE,DAR,FFE,RPL,FT,TWL,RJ,XR,RC, 3 : * FU,P,FF,XE,YE,ZE) 4 : C 5 : C ***** ELEMENT COQ3 ****** 6 : C CE SUB A COMMME BUT DE CALCULER LES FORCES ( EXPRIMEES DANS LE 7 : C REOERE GLOBALE) QUI EQUILIBRENT UN CHAMP DE CONTRAINTES POUR 8 : C 9 : C EN ENTREE - EP : EPAISSEUR DE LA COQUE 10 : C XXE(I,J) : IEME COORDONNEE DU JEME NOEUD DE L'ELEMENT 11 : C DAR(6) : CONTRAINTES - MEMBRANE PUIS FLEXION 12 : C 13 : C EN SORTIE - FFE(I,J) : IEME COMPOSANTE DU CHAMP DE FORCE A 14 : C APPLIQUER AU JEME NOEUD DE L'ELEMENT 15 : C 16 : C CE SUB VIENT DE TRICO VIA BILBO. 17 : C 18 : C LA LONGUEUR DU SEGMENT DE TRAVAIL NECESSAIRE EST 151 19 : C ATTENTION RPL ET FT SONT EN EQUIVALLENCE AINSI QUE FF(1),XE(1) 20 : C FF(4),YE(1) FF(7),ZE(1) 21 : C 22 : IMPLICIT INTEGER(I-N) 23 : IMPLICIT REAL*8(A-H,O-Z) 24 : DIMENSION FFE(6,2),XXE(3,3),DAR(6) 25 : DIMENSION TWL(6) 26 : DIMENSION RL(9,9),RP(9,9),RPL(9,9),JLIM(2,9),XL(6,3),RJ(2,2), 27 : 1 XR(3,9),FF(9),RC(3,3) 28 : DIMENSION P(3,3),XE(3),YE(3),ZE(3) 29 : DIMENSION FT(18) 30 : DIMENSION FU(6) 31 : DATA JLIM/1,1,2,3,2,3,4,4,5, 6,5,6,7,7,8,9,8,9/ 32 : DATA RL/ 0., 0., 0., 0. ,0. ,0. , 1. ,0. ,0. , 33 : 1 1., 0., 0., 0. ,0. ,0. , -1. ,0. ,0. , 34 : 2 0., 0., 0., 1. ,0. ,0. , -1. ,0. ,0. , 35 : 3 0., -0.5,0.5,0.,0.5,-0.5, 0., 0., 0., 36 : 4 0., 0., 0., 0., 0., -0.5, 0., 0., 0.5, 37 : 6 0.,-0.5, 0., 0., 0., 0., 0., 0.5, 0., 38 : 7 -1., 0.5, -0.5, 1., 0.5, -0.5, 0., 0., 0., 39 : 8 0., 0., 0., -1., 0., 0.5, 1., 0., 0.5, 40 : 9 1., -0.5, 0., 0., 0., 0., -1., -0.5, 0./ 41 : DATA RP/ 1., 0., 0., 0., 0., 0., 0., 0., 0., 42 : 2 0., 0., 0., 0., 0., 0., 0., 0., 0., 43 : 3 0., 0., 0., 0., 0., 0., 0., 0., 0., 44 : 4 0., 0., 0., 1., 0., 0., 0., 0., 0., 45 : 5 0., 0., 0., 0., 0., 0., 0., 0., 0., 46 : 6 0., 0., 0., 0., 0., 0., 0., 0., 0., 47 : 7 0., 0., 0., 0., 0., 0., 1., 0., 0., 48 : 8 0., 0., 0., 0., 0., 0., 0., 0., 0., 49 : 9 0., 0., 0., 0., 0., 0., 0., 0., 0./ 50 : DATA XL/0.,0.,-1., -0.33333333,0.33333333,0., 51 : 2 0., -1., 0., 0.33333333,0.,-0.33333333, 52 : 3 0.5,-0.5,-0.5,0.,0.33333333,-0.33333333/ 53 : DO 185 IP=1,3 54 : XE(IP)=XXE(1,IP) 55 : YE(IP)=XXE(2,IP) 56 : 185 ZE(IP)=XXE(3,IP) 56 : CALL PASSA(XE,YE,ZE,P,X13,X23,Y13,Y23) 57 : 184 CONTINUE 58 : IELE=3 59 : DO 16 IP=1,3 60 : 16 TWL(IP)=DAR(IP+3) 61 : DO 160 IP=1,3 62 : 160 TWL(3+IP)=DAR(IP) 63 : SO= 0.5D0 64 : FU(1)=SO*(TWL(4)*Y23-TWL(6)*X23) 65 : FU(2)=SO*(-TWL(5)*X23+TWL(6)*Y23) 66 : FU(3)=SO*(-TWL(4)*Y13+TWL(6)*X13) 67 : FU(4)=SO*(TWL(5)*X13-TWL(6)*Y13) 68 : FU(5)=SO*TWL(6)*(X23-X13) 69 : FU(6)=SO*(X23-X13)*TWL(5) 70 : RJ(1,1)=X13 71 : RJ(2,1)=X23 72 : RJ(1,2)=Y13 73 : RJ(2,2)=Y23 74 : DT=X13*Y23-Y13*X23 75 : SURF2= ABS(DT) 76 : DT=-1./ABS(DT) 77 : DO 1 K=1,3 78 : IK=3*K-2 79 : DO 1 I=1,2 80 : II=IK+I 81 : DO 1 J=1,2 82 : IJ=IK+J 83 : 1 RP(IJ,II)=RJ(J,I) 84 : DO 2 J=1,9 85 : K1=JLIM(1,J) 86 : K2=JLIM(2,J) 87 : DO 2 I=4,9 88 : S=0.D0 89 : DO 3 K=K1,K2 90 : 3 S=S+RL(K,I)*RP(K,J) 91 : RPL(I,J)=S 92 : 2 CONTINUE 93 : DO 4 J=1,3 94 : DO 4 I=1,9 95 : S=0.D0 96 : DO 5 K=4,9 97 : L=K-3 98 : 5 S=S+XL(L,J)*RPL(K,I) 99 : 4 XR(J,I)=S 100 : DO 8 I=1,9 101 : 8 FF(I)=0.D0 102 : RC(1,1)=Y23*Y23 103 : RC(2,1)=Y13*Y13 104 : RC(3,1)=-2.*Y13*Y23 105 : RC(1,2)=X23*X23 106 : RC(2,2)=X13*X13 107 : RC(3,2)=-2.*X13*X23 108 : RC(1,3)=-2.*X23*Y23 109 : RC(2,3)=-2.*X13*Y13 110 : RC(3,3)=2.*(X13*Y23+X23*Y13) 111 : DO 9 I=1,3 112 : DO 9 J=1,9 113 : S=0.D0 114 : DO 10 K=1,3 115 : 10 S=S+RC(K,I)*XR(K,J) 116 : 9 FF(J)=FF(J)+S*TWL(I) 117 : K=0 118 : DO 11 I=1,IELE 119 : IP=6*(I-1) 120 : FT(IP+1)=FU(2*I-1) 121 : FT(IP+2)=FU(2*I) 122 : FT(IP+3)=DT*FF(3*I-2) 123 : FT(IP+4)=DT*FF(3*I) 124 : FT(IP+5)=-DT*FF(3*I-1) 125 : FT(IP+6)=0. 126 : 11 CONTINUE 127 : DO 12 I=1,IELE 128 : KP=6*(I-1) 129 : MP=KP+3 130 : DO 13 J=1,3 131 : DO 14 JP=1,3 132 : FFE(J,I)=FFE(J,I)+P(JP,J)*FT(JP+KP) 133 : FFE(J+3,I)=FFE(J+3,I)+P(JP,J)*FT(JP+MP) 134 : 14 CONTINUE 135 : 13 CONTINUE 136 : 12 CONTINUE 137 : RETURN 138 : END 139 :
© Cast3M 2003 - All rights reserved.
Disclaimer