1 : C MASBA3 SOURCE CHAT 05/01/13 01:30:56 5004 2 : C MASBA3 SOURCE DJER 93/03/02 22:53:00 3 : SUBROUTINE MASBA3(REL,LRE,ROS,XE,TGL,KERRE) 4 : C----------------------------------------------------------------------- 5 : C ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE DE L ELEMENT DE 6 : C BARRE A 3 NOEUDS 7 : C 8 : C ENTREE 9 : C LRE DIMENSION DU TABLEAU REL 10 : C ROS RHO * SECT 11 : C XE COORDONNEES DES 3 NOEUDS 12 : C 13 : C SORTIE 14 : C REL MATRICE DE MASSE ELEMENTAIRE 15 : C----------------------------------------------------------------------- 16 : IMPLICIT INTEGER(I-N) 17 : IMPLICIT REAL*8(A-H,O-Z) 18 : -INC CCOPTIO 19 : DIMENSION REL(LRE,LRE),TGL(3,LRE),XE(3,3) 20 : dimension xma(3,3),xga(3),pga(3),pn(3,3),dn(3,3),xj(3) 21 : pga(1)=0.5555555555555555D0 *ros 22 : pga(3)=pga(1) 23 : pga(2)=0.8888888888888888d0*ros 24 : xga(1)=-0.774596669241483D0 25 : xga(3)=-xga(1) 26 : xga(2)=0.d0 27 : do 100 i=1,3 28 : xx = xga(i) 29 : pn(i,1)= 0.5d0 * xx * ( xx -1.d0) 30 : pn(i,2)= -(xx -1.d0)*(xx+1.d0) 31 : pn(i,3)= 0.5d0 *xx*(xx+1.d0) 32 : dn(i,1)=xx-0.5d0 33 : dn(i,2)= -2.d0* xx 34 : dn(i,3)=xx +0.5d0 35 : 100 continue 36 : do 101 i=1,3 37 : dx = dn(i,1)*xe(1,1)+dn(i,2)*xe(1,2)+dn(i,3)*xe(1,3) 38 : dy = dn(i,1)*xe(2,1)+dn(i,2)*xe(2,2)+dn(i,3)*xe(2,3) 39 : dz = dn(i,1)*xe(3,1)+dn(i,2)*xe(3,2)+dn(i,3)*xe(3,3) 40 : ds2 = dx*dx + dy * dy + dz*dz 41 : ds = sqrt(ds2) 42 : xj(i)=ds 43 : 101 continue 44 : do 103 i=1,3 45 : do 103 j=1,3 46 : 103 xma(i,j)=0.d0 47 : do 102 ipg=1,3 48 : xjpg=xj(ipg)*pga(ipg) 49 : do 102 i=1,3 50 : do 102 j=1,3 51 : xma(i,j)=xma(i,j)+ xjpg*pn(ipg,i)*pn(ipg,j) 52 : 102 continue 53 : C 54 : KERRE=0 55 : C 56 : do 60 I=1,3 57 : i1=(i-1)*3 58 : do 60 j=1,3 59 : j1=(j-1)*3 60 : do 60 k=1,3 61 : ii= i1 +k 62 : jj = j1 + k 63 : rel(ii,jj)=xma(i,j) 64 : 60 continue 65 : RETURN 66 : 999 KERRE=1 67 : RETURN 68 : END 69 : 70 :
© Cast3M 2003 - All rights reserved.
Disclaimer