1 : C MATHOO SOURCE CHAT 05/01/13 01:34:21 5004 2 : C MATHOO SOURCE INSL 24/10/96 3 : SUBROUTINE MATHOO(D,EX,PX,NSTRS,IFOUR) 4 : C 5 : C FORMATION DE LA MATRICE D ELASTICITE LINEAIRE 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 : E1=EX/(1.D0+PX)/(1.D0-2.D0*PX) 29 : D5(1,1)=E1*(1.D0-PX) 30 : D5(1,2)=E1*PX 31 : D5(1,3)=D5(1,2) 32 : D5(2,1)=D5(1,2) 33 : D5(2,2)=D5(1,1) 34 : D5(2,3)=D5(1,3) 35 : D5(3,1)=D5(1,2) 36 : D5(3,2)=D5(1,2) 37 : D5(3,3)=D5(1,1) 38 : D5(4,4)=E1*(1.D0-2.D0*PX)/2.D0 39 : D5(5,5)=D5(4,4)/1.2D0 40 : D5(6,6)=D5(5,5) 41 : GO TO 100 42 : C 43 : 50 CONTINUE 44 : IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN 45 : C 46 : C CONTRAINTES PLANES AVEC TOUS LES CISAILLEMENTS 47 : C 48 : E1=EX/(1.D0-PX*PX) 49 : D5(1,1)=E1 50 : D5(1,2)=PX*E1 51 : D5(2,1)=PX*E1 52 : D5(2,2)=E1 53 : D5(3,3)=E1*(1.D0-PX)/2.D0 54 : D5(4,4)=D5(3,3)/1.2D0 55 : D5(5,5)=D5(4,4) 56 : C 57 : ELSE IF(IFOUR.EQ.-1) THEN 58 : C 59 : C DEFORMATIONS PLANES AVEC TOUS LES CISAILLEMENTS 60 : C 61 : E1=EX/(1.D0+PX)/(1.D0-2.D0*PX) 62 : D5(1,1)=E1*(1.D0-PX) 63 : D5(1,2)=PX*E1 64 : D5(2,1)=PX*E1 65 : D5(2,2)=E1*(1.D0-PX) 66 : D5(3,3)=E1*(1.D0-2.D0*PX)/2.D0 67 : D5(4,4)=D5(3,3)/1.2D0 68 : D5(5,5)=D5(4,4) 69 : ELSE 70 : GO TO 1000 71 : ENDIF 72 : GO TO 100 73 : C 74 : 40 CONTINUE 75 : IF(IFOUR.EQ.0.OR.IFOUR.EQ.-1) THEN 76 : C 77 : C CAS AXISYMETRIQUE OU DEFORMATION PLANES 78 : C 79 : E1=EX/(1.D0+PX)/(1.D0-2.D0*PX) 80 : D5(1,1)=E1*(1.D0-PX) 81 : D5(1,2)=E1*PX 82 : D5(1,3)=E1*PX 83 : D5(2,1)=E1*PX 84 : D5(2,2)=E1*(1.D0-PX) 85 : D5(2,3)=E1*PX 86 : D5(3,1)=E1*PX 87 : D5(3,2)=E1*PX 88 : D5(3,3)=E1*(1.D0-PX) 89 : D5(4,4)=E1*(1.D0-2.D0*PX)/2.D0 90 : C 91 : ELSE IF (IFOUR.EQ.-2) THEN 92 : C 93 : C CONTRAINTES PLANES 94 : C 95 : E1=EX/(1.D0-PX*PX) 96 : D5(1,1)=E1 97 : D5(1,2)=PX*E1 98 : D5(2,1)=PX*E1 99 : D5(2,2)=E1 100 : D5(3,3)=D5(1,1)*1.D-6 101 : D5(4,4)=E1*(1.D0-PX)/2.D0 102 : C 103 : ENDIF 104 : GO TO 100 105 : C 106 : 30 CONTINUE 107 : IF(IFOUR.EQ.-2.OR.IFOUR.EQ.2) THEN 108 : C 109 : C CONTRAINTES PLANES SANS CISAILLEMENTS TRANSVERSAL 110 : C 111 : E1=EX/(1.D0-PX*PX) 112 : D5(1,1)=E1 113 : D5(1,2)=PX*E1 114 : D5(2,1)=PX*E1 115 : D5(2,2)=E1 116 : D5(3,3)=E1*(1.D0-PX)/2.D0 117 : C 118 : ELSE IF(IFOUR.EQ.-1) THEN 119 : C 120 : C DEFORMATIONS PLANES SANS CISAILLEMENTS TRANSVERSAL 121 : C 122 : E1=EX/(1.D0+PX)/(1.D0-2.D0*PX) 123 : D5(1,1)=E1*(1.D0-PX) 124 : D5(1,2)=PX*E1 125 : D5(2,1)=PX*E1 126 : D5(2,2)=E1*(1.D0-PX) 127 : D5(3,3)=E1*(1.D0-2.D0*PX)/2.D0 128 : ELSE 129 : GO TO 1000 130 : ENDIF 131 : GO TO 100 132 : C 133 : 10 CONTINUE 134 : C 135 : C COMPORTEMENT UNIDIMENSIONNEL 136 : C 137 : D5(1,1)=EX 138 : GO TO 100 139 : 1000 WRITE(*,200) IFOUR,NSTRS 140 : STOP 141 : 100 CONTINUE 142 : C 143 : DO 300 I = 1,NSTRS 144 : DO 400 J = 1,NSTRS 145 : D(I,J) = D5(I,J) 146 : 400 CONTINUE 147 : 300 CONTINUE 148 : 200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS MATHOOK',/, 149 : * 10X,'IFOUR=',I2,' NSTRS=',I2,//) 150 : RETURN 151 : END 152 : 153 : 154 : 155 :
© Cast3M 2003 - All rights reserved.
Disclaimer