1 : C POCOLO SOURCE CHAT 05/01/13 02:15:54 5004 2 : SUBROUTINE POCOLO(V,CONT,XL,YOUN,CISAIL,SD,SREDY,SREDZ, 3 : . XINX,XINY,XINZ,IDEUX) 4 : C----------------------------------------------------------------------- 5 : C 6 : C ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE 7 : C D APRES TEDEL 8 : C 9 : C ENTREE 10 : C V VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE 11 : C XL LONGUEUR 12 : C SD SECTION DROITE 13 : C SREDY SECTION REDUITE A L EFFORT TRANCHANT TY 14 : C SREDZ SECTION REDUITE A L EFFORT TRANCHANT TZ 15 : C XINX INERTIE IX 16 : C XINY INERTIE IY 17 : C XINZ INERTIE IZ 18 : C YOUN MODULE D YOUNG 19 : C CISAIL MODULE DE CISAILLEMENT 20 : C IDEUX INDICE POUR DEFORMATIONS DU SECOND ORDRE 21 : C SORTIE 22 : C CONT CONTRAINTES AUX 2 NOEUDS 6 PAR NOEUD 23 : C SEULS LES MOMENTS DE FLEXION VARIENT LINEAIREMENT 24 : C 25 : C AOUT 84 26 : C----------------------------------------------------------------------- 27 : IMPLICIT INTEGER(I-N) 28 : IMPLICIT REAL*8(A-H,O-Z) 29 : DIMENSION CONT(*),V(*) 30 : C 31 : RL=1.D00/XL 32 : RL2=RL*RL 33 : RL3= RL2*RL 34 : X2=0. 35 : X3=0. 36 : IF( SREDZ. NE.0.) X2= YOUN/CISAIL*RL2*XINY/SREDZ 37 : IF( SREDY. NE.0.) X3= YOUN/CISAIL*RL2*XINZ/SREDY 38 : X13=1./(1.+12.*X3) 39 : X33=(1.+3.*X3)*X13 40 : X63=(1.-6.*X3)*X13 41 : X12=1./(1.+12.*X2) 42 : X32=(1.+3.*X2)*X12 43 : X62=(1.-6.*X2)*X12 44 : C 45 : C DEFORMATIONS GENERALISEES 46 : C 47 : EPS =(V(7)-V(1))*RL 48 : XSIX =(V(10)-V(4))*RL 49 : XSIY1= 50 : .RL2*( 6.*V(3)*X12-4.*XL*X32*V(5)-6.*V(9)*X12-2.*XL*X62*V(11)) 51 : XSIY2= 52 : .RL2*(-6.*V(3)*X12+2.*XL*X62*V(5)+6.*V(9)*X12+4.*XL*X32*V(11)) 53 : XSIZ1= 54 : .RL2*(-6.*V(2)*X13-4.*XL*X33*V(6)+6.*V(8)*X13-2.*XL*X63*V(12)) 55 : XSIZ2= 56 : .RL2*( 6.*V(2)*X13+2.*XL*X63*V(6)-6.*V(8)*X13+4.*XL*X33*V(12)) 57 : C 58 : C TERMES DU SECOND ORDRE 59 : C 60 : IF(IDEUX.EQ.0) THEN 61 : EPS1=EPS 62 : EPS2=EPS 63 : ELSE IF(IDEUX.NE.0) THEN 64 : EPS1=EPS+0.5D0*(EPS**2+V(5)**2+V(6)**2) 65 : EPS2=EPS+0.5D0*(EPS**2+V(11)**2+V(12)**2) 66 : ENDIF 67 : C 68 : C CONTRAINTES GENERALISEES 69 : C 70 : SIG1 =YOUN*SD*(EPS1+EPS2)*0.5D0 71 : SIG2 =SIG1 72 : TOR =CISAIL*XINX*XSIX 73 : XMOMZ1=YOUN*XINZ*XSIZ1 74 : XMOMZ2=YOUN*XINZ*XSIZ2 75 : XMOMY1=YOUN*XINY*XSIY1 76 : XMOMY2=YOUN*XINY*XSIY2 77 : TRANY =(XMOMZ1-XMOMZ2)*RL 78 : TRANZ =(XMOMY2-XMOMY1)*RL 79 : C 80 : C REMPLISSAGE DU TABLEAU CONT 81 : C 82 : CONT(1)=SIG1 83 : CONT(2)=TRANY 84 : CONT(3)=TRANZ 85 : CONT(4)=TOR 86 : CONT(5)=XMOMY1 87 : CONT(6)=XMOMZ1 88 : CONT(7)=SIG2 89 : CONT(8)=TRANY 90 : CONT(9)=TRANZ 91 : CONT(10)=TOR 92 : CONT(11)=XMOMY2 93 : CONT(12)=XMOMZ2 94 : RETURN 95 : END 96 :
© Cast3M 2003 - All rights reserved.
Disclaimer