1 : C BDKT2 SOURCE CHAT 05/02/25 21:15:02 5035 2 : SUBROUTINE BDKT2(XEL,XDDL,IGAU,BGENE,XSTRS) 3 : C====================================================================== 4 : C 5 : C AJOUTE DANS LES DEFORMATIOND XSTRS LA PARTIE EPS2 6 : C QUE POUR LES ELEMENTS DKT 7 : C ENTREE 8 : C XEL(3,3) =COORDONNEES LOCALES DU DKT 9 : C XDDL(18) =D D L LOCAUX DE LA FORMULATION 10 : C IGAU =NUMERO DU POINT DE GAUSS DE 1 A 3 11 : C XSTRS(6) =DEFORMATIONS (EPS1 ) 12 : C BGENE(6,18) =MATRICE BGENE reliant deformations et XDDL 13 : C SORTIE 14 : C XSTRS(6) =DEFORMATIONS (EPS1 + EPS2 ) 15 : C 16 : C CODE EBERSOLT AOUT 85 17 : C======================================================================= 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : PARAMETER(XZER=0.D0,UNDEMI=.5D0) 21 : DIMENSION XEL(3,*),XSTRS(*),XDDL(*) 22 : DIMENSION SHP(6,9),ROT(9),BGENE(6,18) 23 : C 24 : DUDX=XZER 25 : DUDY=XZER 26 : DVDX=XZER 27 : DVDY=XZER 28 : DWDX=XZER 29 : DWDY=XZER 30 : C 30 : CALL DKTSHP(IGAU,XEL,SHP,DJAC) 31 : C 32 : C ON EXTRAIT LES DEGRES W ROT X ROT Y 33 : C 34 : IACC= 0 35 : IADD= 2 36 : DO 100 IA=1,3 37 : DO 110 IB=1,3 38 : ROT(IACC+IB)=XDDL(IADD+IB) 39 : 110 CONTINUE 40 : IACC=IACC+3 41 : IADD=IADD+6 42 : 100 CONTINUE 43 : C 44 : DO 200 IA=1,9 45 : DWDX =SHP(2,IA)*ROT(IA) +DWDX 46 : DWDY =SHP(3,IA)*ROT(IA) +DWDY 47 : 200 CONTINUE 48 : C 49 : C--Calcul des derivees de u et v 50 : DO 300 IA=0,2 51 : DUDX=DUDX + BGENE(1,1+6*IA)*XDDL(1+6*IA) 52 : DUDY=DUDY + BGENE(3,1+6*IA)*XDDL(1+6*IA) 53 : DVDX=DVDX + BGENE(3,2+6*IA)*XDDL(2+6*IA) 54 : DVDY=DVDY + BGENE(2,2+6*IA)*XDDL(2+6*IA) 55 : 300 CONTINUE 56 : C 57 : C--Mise a jour des deformations au deuxieme ordre 58 : XSTRS(1)=XSTRS(1)+UNDEMI*(DUDX*DUDX+DVDX*DVDX+DWDX*DWDX) 59 : XSTRS(2)=XSTRS(2)+UNDEMI*(DUDY*DUDY+DVDY*DVDY+DWDY*DWDY) 60 : XSTRS(3)=XSTRS(3)+DUDX*DUDY+DVDX*DVDY+DWDX*DWDY 61 : C 62 : RETURN 63 : END 64 :
© Cast3M 2003 - All rights reserved.
Disclaimer