1 : C BGRMAT SOURCE CHAT 05/01/12 21:39:51 5004 2 : SUBROUTINE BGRMAT(SX,SY,HXABS,HYABS,HXORD,HYORD,B) 3 : C 4 : C ********************************************************* 5 : C ce sous programme calcule la matrice [BGR] qui 6 : C lie les courbures aux d.d.l. de l element 7 : C ------- 8 : C ********************************************************* 9 : IMPLICIT INTEGER(I-N) 10 : IMPLICIT REAL*8 (A-H,O-Z) 11 : DIMENSION SX(3),SY(3),HXABS(9),HYABS(9),HXORD(9),HYORD(9),B(4,9) 12 : DIMENSION B11(9),B12(9),B21(9),B22(9),B31(9),B32(9),B33(9),B34(9) 13 : DO 10 I=1,9 14 : B11(I)=SY(2)*HXABS(I) 15 : B12(I)=SY(3)*HXORD(I) 16 : B21(I)=-SX(2)*HYABS(I) 17 : B22(I)=-SX(3)*HYORD(I) 18 : B31(I)=-SX(2)*HXABS(I) 19 : B32(I)=-SX(3)*HXORD(I) 20 : B33(I)=SY(2)*HYABS(I) 21 : B34(I)=SY(3)*HYORD(I) 22 : 10 CONTINUE 23 : C ..................................................... 24 : AIR2=DABS(SX(2)*SY(3)-SX(3)*SY(2)) 24 : CALL ZERO(B,4,9) 25 : DO 20 I=1,9 26 : B(1,I)=(B11(I)+B12(I))/AIR2 27 : B(4,I)=(B21(I)+B22(I))/AIR2 28 : B(2,I)=(B31(I)+B32(I))/AIR2 29 : B(3,I)=(B33(I)+B34(I))/AIR2 30 : 20 CONTINUE 31 : RETURN 32 : END 33 : 34 :
© Cast3M 2003 - All rights reserved.
Disclaimer