1 : C POCOL2 SOURCE CHAT 05/01/13 02:15:51 5004 2 : C POCOL2 SOURCE TC1 92/03/16 21:04:19 501 3 : SUBROUTINE POCOL2(V,CONT,XL,YOUN,CISAIL,SD,SREDY,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 XINZ INERTIE IZ 15 : C YOUN MODULE D YOUNG 16 : C CISAIL MODULE DE CISAILLEMENT 17 : C IDEUX INDICE POUR DEFORMATIONS DU SECOND ORDRE 18 : C SORTIE 19 : C CONT CONTRAINTES AUX 2 NOEUDS 6 PAR NOEUD 20 : C SEULS LES MOMENTS DE FLEXION VARIENT LINEAIREMENT 21 : C 22 : C AOUT 84 23 : C----------------------------------------------------------------------- 24 : IMPLICIT INTEGER(I-N) 25 : IMPLICIT REAL*8(A-H,O-Z) 26 : DIMENSION CONT(*),V(*) 27 : C 28 : RL=1.D00/XL 29 : RL2=RL*RL 30 : RL3= RL2*RL 31 : X2=0. 32 : X3=0. 33 : IF( SREDY. NE.0.) X3= YOUN/CISAIL*RL2*XINZ/SREDY 34 : X13=1./(1.+12.*X3) 35 : X33=(1.+3.*X3)*X13 36 : X63=(1.-6.*X3)*X13 37 : X12=1./(1.+12.*X2) 38 : X32=(1.+3.*X2)*X12 39 : X62=(1.-6.*X2)*X12 40 : C 41 : C DEFORMATIONS GENERALISEES 42 : C 43 : EPS =(V(4)-V(1))*RL 44 : XSIZ1= 45 : .RL2*(-6.*V(2)*X13-4.*XL*X33*V(3)+6.*V(5)*X13-2.*XL*X63*V(6)) 46 : XSIZ2= 47 : .RL2*( 6.*V(2)*X13+2.*XL*X63*V(3)-6.*V(5)*X13+4.*XL*X33*V(6)) 48 : C 49 : C TERMES DU SECOND ORDRE 50 : C 51 : IF(IDEUX.EQ.0) THEN 52 : EPS1=EPS 53 : EPS2=EPS 54 : ELSE IF(IDEUX.NE.0) THEN 55 : EPS1=EPS+0.5D0*(EPS**2+V(3)**2) 56 : EPS2=EPS+0.5D0*(EPS**2+V(6)**2) 57 : ENDIF 58 : C 59 : C CONTRAINTES GENERALISEES 60 : C 61 : SIG1 =YOUN*SD*(EPS1+EPS2)*0.5D0 62 : SIG2 =SIG1 63 : XMOMZ1=YOUN*XINZ*XSIZ1 64 : XMOMZ2=YOUN*XINZ*XSIZ2 65 : TRANY =(XMOMZ1-XMOMZ2)*RL 66 : C 67 : C REMPLISSAGE DU TABLEAU CONT 68 : C 69 : CONT(1)=SIG1 70 : CONT(2)=TRANY 71 : CONT(3)=XMOMZ1 72 : CONT(4)=SIG2 73 : CONT(5)=TRANY 74 : CONT(6)=XMOMZ2 75 : RETURN 76 : END 77 : 78 :
© Cast3M 2003 - All rights reserved.
Disclaimer