1 : C COQ8RI SOURCE CHAT 05/01/12 22:22:13 5004 2 : SUBROUTINE COQ8RI(B,FACT,D,LRE,NBPGAU,IA,NBNN,RIG) 3 : C======================================================================= 4 : C ENTREES 5 : C B(5,LRE) : MATRICE DE DERIVEES DE FONCTIONS DE FORMES 6 : C FACT : DETERMINANT DU JACOBIEN * POIDS DE GAUSS 7 : C D(5,5) : MATRICE DE HOOKE 8 : C LRE : NOMBRE DE DDL 9 : C NBPGAU : NOMBRE DE POINTS DE GAUSS 10 : C IA : NUMERO DE POINT DE GAUSS 11 : C NBNN : NOMBRE DE NOEUDS 12 : C 13 : C SORTIES 14 : C RIG(LRE,LRE): MATRICE DE RIGIDITE ELEMENTAIRE 15 : C 16 : C RECUPERATION BILBO OCTOBRE 85 17 : C======================================================================= 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8 (A-H,O-Z) 20 : PARAMETER (XZER=0.D0) 21 : DIMENSION RIG(LRE,*),D(5,*),B(5,*) 22 : DIMENSION DAM(48,5) 23 : C 24 : C CALCUL DE BTDB 25 : C 26 : DO 60 I =1,LRE 27 : DO 60 J = 1,5 28 : DAM(I,J) = XZER 29 : DO 60 K = 1,5 30 : DAM(I,J) = DAM(I,J) + B(K,I)*D(K,J) 31 : 60 CONTINUE 32 : DO 70 I = 1,LRE 33 : DO 70 J = I,LRE 34 : DO 70 K = 1,5 35 : RIG(I,J) = RIG(I,J) + DAM(I,K)*B(K,J)*FACT 36 : 70 CONTINUE 37 : C 38 : IF (IA.EQ.NBPGAU)THEN 39 : DO 110 I = 1,LRE 40 : DO 110 J = I,LRE 41 : 110 RIG(J,I) = RIG(I,J) 42 : C 43 : C ON AJOUTE LES COEFF POUR LA TORSION NORMALE AU PLAN DE LA COQUE 44 : C 45 : VMAX = XZER 46 : DO 120 IK = 1,NBNN 47 : JK4 = (IK-1)*6 + 4 48 : JK6 = JK4 + 2 49 : DO 120 JK = JK4,JK6 50 : VMIN = RIG(JK,JK) 51 : IF (VMAX.GE.VMIN) GO TO 120 52 : VMAX = RIG(JK,JK) 53 : 120 CONTINUE 54 : DO 130 IK = 1,NBNN 55 : JK4 = (IK-1)*6 + 4 56 : JK6 = JK4 + 2 57 : DO 130 JK = JK4,JK6 58 : 130 RIG(JK,JK)=RIG(JK,JK)+1.D-9*VMAX 59 : ENDIF 60 : C 61 : RETURN 62 : END 63 :
© Cast3M 2003 - All rights reserved.
Disclaimer