1 : C MATSOU SOURCE CHAT 05/01/13 01:35:34 5004 2 : C MATSOU SOURCE INSL 24/10/96 3 : SUBROUTINE MATSOU(D,EX,PX,NSTRS,IFOUR) 4 : C 5 : C FORMATION DE LA MATRICE DE SOUPLESSE 6 : C----------------------------------------------------------------------- 7 : C VARIABLES PASSEES PAR LES COMMONS COPTIO , ECOU ET NECOU 8 : C 9 : C IFOUR INDICE DU TYPE DE PROBLEME 10 : C -2 CONTRAINTES PLANES 11 : C -1 DEFORMATIONS PLANES 12 : C 0 AXISYMETRIQUE 13 : C 1 SERIE DE FOURIER 14 : C 2 TRIDIMENSIONNEL 15 : C----------------------------------------------------------------------- 16 : IMPLICIT INTEGER(I-N) 17 : IMPLICIT REAL*8(A-H,O-Z) 18 : DIMENSION D5(6,6),D(NSTRS,NSTRS) 19 : C 19 : CALL ZERO(D5,6,6) 19 : CALL ZERO(D,NSTRS,NSTRS) 20 : C 21 : GO TO (10,1000,30,40,50,60),NSTRS 22 : GOTO 1000 23 : C 24 : 60 CONTINUE 25 : C 26 : C COMPORTEMENT TRIDIMENSIONNEL 27 : C 28 : D5(1,1)=1.D0/EX 29 : D5(1,2)=-PX/EX 30 : D5(1,3)=-PX/EX 31 : D5(2,1)=-PX/EX 32 : D5(2,2)=1.D0/EX 33 : D5(2,3)=-PX/EX 34 : D5(3,1)=-PX/EX 35 : D5(3,2)=-PX/EX 36 : D5(3,3)=1.D0/EX 37 : D5(4,4)=2.D0*(1.D0+PX)/EX 38 : D5(5,5)=1.2D0*D5(4,4) 39 : D5(6,6)=D5(5,5) 40 : GO TO 100 41 : C 42 : 50 CONTINUE 43 : IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN 44 : C 45 : C CONTRAINTES PLANES AVEC TOUS LES CISAILLEMENTS 46 : C 47 : D5(1,1)=1.D0/EX 48 : D5(1,2)=-PX/EX 49 : D5(2,1)=-PX/EX 50 : D5(2,2)=1.D0/EX 51 : D5(3,3)=2.D0*(1.D0+PX)/EX 52 : D5(4,4)=D5(3,3)*1.2D0 53 : D5(5,5)=D5(4,4) 54 : C 55 : ELSE IF(IFOUR.EQ.-1) THEN 56 : C 57 : C DEFORMATIONS PLANES AVEC TOUS LES CISAILLEMENTS 58 : C 59 : E1=(1.D0+PX)/EX 60 : D5(1,1)=E1*(1.D0-PX) 61 : D5(1,2)=-PX*E1 62 : D5(2,1)=-PX*E1 63 : D5(2,2)=E1*(1.D0-PX) 64 : D5(3,3)=E1*2.D0 65 : D5(4,4)=D5(3,3)*1.2D0 66 : D5(5,5)=D5(4,4) 67 : ELSE 68 : GO TO 1000 69 : ENDIF 70 : GO TO 100 71 : C 72 : 40 CONTINUE 73 : IF(IFOUR.EQ.0.OR.IFOUR.EQ.-1) THEN 74 : C 75 : C CAS AXISYMETRIQUE OU DEFORMATION PLANES 76 : C 77 : D5(1,1)=1.D0/EX 78 : D5(1,2)=-PX/EX 79 : D5(1,3)=-PX/EX 80 : D5(2,1)=-PX/EX 81 : D5(2,2)=1.D0/EX 82 : D5(2,3)=-PX/EX 83 : D5(3,1)=-PX/EX 84 : D5(3,2)=-PX/EX 85 : D5(3,3)=1.D0/EX 86 : D5(4,4)=2.D0*(1.D0+PX)/EX 87 : C 88 : ELSE IF(IFOUR.EQ.-2) THEN 89 : C 90 : C CONTRAINTES PLANES 91 : C 92 : D5(1,1)=1.D0/EX 93 : D5(1,2)=-PX/EX 94 : D5(2,1)=-PX/EX 95 : D5(2,2)=1.D0/EX 96 : D5(3,3)=D5(1,1)*1.D+6 97 : D5(4,4)=2.D0*(1.D0+PX)/EX 98 : C 99 : ENDIF 100 : GO TO 100 101 : C 102 : 30 CONTINUE 103 : IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN 104 : C 105 : C CONTRAINTES PLANES SANS CISAILLEMENTS TRANSVERSAL 106 : C 107 : D5(1,1)=1.D0/EX 108 : D5(1,2)=-PX/EX 109 : D5(2,1)=-PX/EX 110 : D5(2,2)=1.D0/EX 111 : D5(3,3)=2.D0*(1.D0+PX)/EX 112 : C 113 : ELSE IF(IFOUR.EQ.-1) THEN 114 : C 115 : C DEFORMATIONS PLANES SANS CISAILLEMENTS TRANSVERSAL 116 : C 117 : E1=(1.D0+PX)/EX 118 : D5(1,1)=E1*(1.D0-PX) 119 : D5(1,2)=-PX*E1 120 : D5(2,1)=-PX*E1 121 : D5(2,2)=E1*(1.D0-PX) 122 : D5(3,3)=E1*2.D0 123 : C 124 : ELSE 125 : GO TO 1000 126 : ENDIF 127 : GO TO 100 128 : C 129 : 10 CONTINUE 130 : C 131 : C COMPORTEMENT UNIDIMENSIONNEL 132 : C 133 : D5(1,1)=1.D0/EX 134 : GO TO 100 135 : 1000 WRITE(*,200) IFOUR,NSTRS 136 : STOP 137 : 100 CONTINUE 138 : DO 300 I = 1,NSTRS 139 : DO 400 J = 1,NSTRS 140 : D(I,J) = D5(I,J) 141 : 400 CONTINUE 142 : 300 CONTINUE 143 : 200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS MATSOUP',/, 144 : * 10X,'IFOUR=',I2,' NSTRS=',I2,//) 145 : RETURN 146 : END 147 : 148 : 149 :
© Cast3M 2003 - All rights reserved.
Disclaimer