Télécharger ared3m.eso

Retour à la liste

Numérotation des lignes :

  1. C ARED3M SOURCE CHAT 05/01/12 21:23:28 5004
  2. SUBROUTINE ARED3M(SIGT,EPST,ENDT,XLIMT,
  3. * SIGD,EPSD,ENDD,XLIMD,
  4. * DEPS,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1,
  5. * VAL2,VAL3,VP2,VP3,EPSTT,FDDD,D,SIFF,VDEII,
  6. * VDEP2,EPSC,PPLUS,U,PROD,ENDOPL,VDEI,KM,AAA)
  7. C======================================================================
  8. C TEST D ENDOMMAGEMENT INITIAL POUR LES ELEMENTS TRIDIMENSIONNELS
  9. C======================================================================
  10. C
  11. C
  12. C MODIF. : F.CORMERY
  13. C E.N.S.M.A - LMPM
  14. C JUIN 1993
  15. C
  16. C======================================================================
  17. C ENTREES : SIGT,EPST,ENDT, VALEURS AU PAS PRECEDENT
  18. C DSIG,DEPS ACCROISSEMENTS INITIAUX
  19. C SORTIES SIGD,EPSD,ENDD, VALEURS A LA FIN DU PAS
  20. C DEPS ACCROISSEMENTS
  21. C VDEP MATRICE TANGENTE
  22. C======================================================================
  23. IMPLICIT INTEGER(I-N)
  24. IMPLICIT REAL*8(A-H,O-Z)
  25. REAL*8 LAMB,MUT
  26. C----------------------------------------------------------------------
  27. C COMMON
  28. C----------------------------------------------------------------------
  29. COMMON /APENTE/ LAMB,MUT,G,ALP1,BETA
  30. C----------------------------------------------------------------------
  31. C DATA ET DIMENSIONS
  32. C----------------------------------------------------------------------
  33. C N9
  34. CC DIMENSION SIGD(6),EPSD(6),ENDD(6),D(3,3),VAL1(3,3),VP1(3),
  35. C N15 N51
  36. CC * SIGT(6),EPST(6),ENDT(6),SIFF(6),VDEII(6,6),
  37. C
  38. CC * DEPS(6),VDEP(6,6),VDEP1(6,6),VAL(3,3),VP(3),
  39. C N87
  40. CC * VAL2(3,3),VAL3(3,3),VP3(3),VP2(3),VDEP2(6,6)
  41. C N93 N174 N255
  42. CC DIMENSION EPSC(6),EPSPLS(6),PPLUS(3,3,3,3),U(3,3,3,3)
  43. C N261 N267 N303 N309
  44. CC DIMENSION PROD(6),ENDOPL(6),VDEI(6,6),KM(6),EPSTT(6)
  45. C
  46. C
  47. DIMENSION AAA(*)
  48. DIMENSION SIGD(*),EPSD(*),ENDD(*),D(3,*),VAL1(3,*),VP1(*),
  49. C N15 N51
  50. * SIGT(*),EPST(*),ENDT(*),SIFF(*),VDEII(6,*),
  51. C
  52. * DEPS(*),VDEP(6,*),VDEP1(6,*),VAL(3,*),VP(*),
  53. C N87
  54. * VAL2(3,*),VAL3(3,*),VP3(*),VP2(*),VDEP2(6,*)
  55. C N93 N174 N255
  56. DIMENSION EPSC(*),EPSPLS(*),PPLUS(3,3,3,*),U(3,3,3,*)
  57. C N261 N267 N303 N309
  58. DIMENSION PROD(*),ENDOPL(*),VDEI(6,*),KM(*),EPSTT(*)
  59. DIMENSION NII(6),MII(6),NIII(6),MIII(6)
  60. DATA NII/1,2,3,3,3,2/
  61. DATA MII/1,2,3,2,1,1/
  62. DATA NIII/1,2,3,2,3,3/
  63. DATA MIII/1,2,3,1,2,1/
  64. DATA ZERO/0.D0/
  65. DATA IMA/6/,DEUX/2.D0/,UD/0.5D0/
  66. C----------------------------------------------------------------------
  67. MT=6
  68. IKK=1
  69. C----------------------------------------------------------------------
  70. C TENSEUR DE DEFORMATION ACTUALISEE
  71. C----------------------------------------------------------------------
  72. MMI=1
  73. GOTO 121
  74. C------------
  75. 3 CONTINUE
  76. C------------
  77. DO 144 I=1,6
  78. M= MII(I)
  79. N=NII(I)
  80. DO 144 J=1,6
  81. K=MII(J)
  82. L=NII(J)
  83. VDEI(I,J)=U(M,N,K,L)
  84. VDEII(I,J)=U(M,N,K,L)
  85. 144 CONTINUE
  86. C---------------------------Calcul de la deformation initiale----------
  87. IF(NPA.EQ.1)THEN
  88. C------------
  89. CALL INVERE(VDEI,6,6,KM,DET)
  90. C------------
  91. DO 4 I=1,6
  92. EPST(I)=0.D0
  93. DO 4 J=1,6
  94. EPST(I)=EPST(I)+VDEI(I,J)*(-G*ENDT(J))
  95. 4 CONTINUE
  96. C------------
  97. MMI=0
  98. ENDIF
  99. C----------------------------------------------------------------------
  100. 234 CONTINUE
  101. DO 10 I=1,IMA
  102. EPSD(I)=EPST(I)+DEPS(I)
  103. EPSC(I)=EPSD(I)
  104. 10 CONTINUE
  105. C----------------------------------------------------------------------
  106. C CALCUL DU TENSEUR DE DEFORMATION POSITIVE
  107. C----------------------------------------------------------------------
  108. 7 CALL EPPLUM(EPSC,PPLUS,1,VAL1,VP1,AAA(442),
  109. . AAA(478),AAA(514),AAA(517),AAA(526),
  110. . AAA(532))
  111. C----------------
  112. EPSD(4)=EPSD(4)/2
  113. EPSD(5)=EPSD(5)/2
  114. EPSD(6)=EPSD(6)/2
  115. C----------------
  116. CALL T4CT2(PPLUS,EPSD,EPSPLS,6)
  117. C----------------
  118. DO 324 I=1,6
  119. IF(ABS(EPSPLS(I)).LE.1E-15)EPSPLS(I)=0.D0
  120. 324 CONTINUE
  121. C----------------
  122. EPSD(4)=EPSD(4)*2
  123. EPSD(5)=EPSD(5)*2
  124. EPSD(6)=EPSD(6)*2
  125. EPSPLS(4)=EPSPLS(4)*2
  126. EPSPLS(5)=EPSPLS(5)*2
  127. EPSPLS(6)=EPSPLS(6)*2
  128. C----------------------------------------------------------------------
  129. C TEST D ENDOMMAGEMENT
  130.  
  131. C----------------------------------------------------------------------
  132. CALL AVENDO(EPSPLS,ENDT,DEQ)
  133. CRIT=DEQ-XLIMT
  134. IF(CRIT.LE.ZERO)GO TO 100
  135. C----------------------------------------------------------------------
  136. C CALCUL DU MULTIPLICATEUR
  137. C----------------------------------------------------------------------
  138. XLMD=ZERO
  139. CALL AMULTM(EPSPLS,ENDT,XLMD,XLIMT)
  140. C----------------------------------------------------------------------
  141. C ACTUALISATION DES GRANDEURS
  142. C----------------------------------------------------------------------
  143. CALL AACTUM(SIGT,EPST,ENDT,XLIMT,SIGD,EPSD,ENDD,XLIMD,
  144. * VDEP,EPSC,EPSPLS,XLMD,PPLUS,
  145. * VDEP1,SIFF,VDEP2,FDDD,AAA(442),AAA(448),
  146. * AAA(454),AAA(460),AAA(466),AAA(472),
  147. * AAA(478),AAA(484),AAA(493),AAA(502),
  148. * AAA(511),AAA(520),AAA(529),AAA(538),
  149. * AAA(547),AAA(628),AAA(709),AAA(790),
  150. * AAA(871))
  151. ISS=2
  152. GO TO 899
  153. C**********************************************************************
  154. C CAS ELASTIQUE
  155. C**********************************************************************
  156. 100 CONTINUE
  157. ISS=1
  158. C----------------------------------------------------------------------
  159. DO 110 I=1,IMA
  160. ENDD(I)=ENDT(I)
  161. 110 CONTINUE
  162. XLIMD=XLIMT
  163. C--------------------------------CALCUL PPLUS:D=D+ ------------
  164. * MODEE1=1
  165. * IF(MODEE1.EQ.1)THEN
  166. DO 123 I=1,6
  167. 123 ENDOPL(I)=ENDD(I)
  168. * GOTO 108
  169. * ENDIF
  170. * CALL T4CT2(PPLUS,ENDD,ENDOPL,6)
  171. C--------------------------------CALCUL DES CONTRAINTES----------------
  172. *108 TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)
  173. TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)
  174. * +EPSD(4)*ENDD(4)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  175. TRAC3=EPSD(1)+EPSD(2)+EPSD(3)
  176. PROD(1)=DEUX*EPSD(1)*ENDD(1)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  177. PROD(2)=DEUX*EPSD(2)*ENDD(2)+EPSD(4)*ENDD(4)+EPSD(6)*ENDD(6)
  178. PROD(3)=DEUX*EPSD(3)*ENDD(3)+EPSD(5)*ENDD(5)+EPSD(4)*ENDD(4)
  179. PROD(4)=UD*EPSD(6)*ENDD(5)+EPSD(2)*ENDD(4)+UD*EPSD(4)*ENDD(3)+
  180. * UD*EPSD(5)*ENDD(6)+UD*EPSD(4)*ENDD(2)+EPSD(3)*ENDD(4)
  181. PROD(5)=UD*EPSD(5)*ENDD(1)+UD*EPSD(4)*ENDD(6)+EPSD(3)*ENDD(5)+
  182. * EPSD(1)*ENDD(5)+UD*EPSD(6)*ENDD(4)+UD*EPSD(5)*ENDD(3)
  183. PROD(6)=UD*EPSD(6)*ENDD(1)+EPSD(2)*ENDD(6)+UD*EPSD(4)*ENDD(5)+
  184. * EPSD(1)*ENDD(6)+UD*EPSD(6)*ENDD(2)+UD*EPSD(5)*ENDD(4)
  185. C---------------------
  186. DO 115 I=1,3
  187. SIGD(I)=G*ENDOPL(I)+LAMB*TRAC3+ALP1*(TRAC2+TRAC3*ENDD(I))
  188. * +DEUX*BETA*PROD(I)+DEUX*MUT*EPSD(I)
  189. 115 CONTINUE
  190. C---------------------
  191. DO 116 I=4,6
  192. SIGD(I)=G*ENDOPL(I)+ALP1*TRAC3*ENDD(I)
  193. * +DEUX*BETA*PROD(I)+MUT*EPSD(I)
  194. 116 CONTINUE
  195. C------------------
  196. DO 17 I=1,6
  197. SIFF(I)=0
  198. DO 17 J=1,6
  199. SIFF(I)=SIFF(I)+VDEII(I,J)*EPSC(J)
  200. 17 CONTINUE
  201. DO 18 I=1,6
  202. SIFF(I)=SIFF(I)+G*ENDOPL(I)
  203. 18 CONTINUE
  204. C----------------------------------------------------------------------
  205. MMI=0
  206. C----------------------------------------------------------------------
  207. C CALCUL MATRICE TANGENTE ELASTIQUE ANISOTROPE
  208. C----------------------------------------------------------------------
  209. C
  210. C---------------------------------Passage en notation indicielle-------
  211. 121 CONTINUE
  212. DO 42 I=1,3
  213. D(I,I)=ENDT(I)
  214. DO 43 J=I+1,3
  215. K = 9 - I -J
  216. * IF ((I.EQ.1).AND.(J.EQ.2)) K=6
  217. * IF ((I.EQ.1).AND.(J.EQ.3)) K=5
  218. * IF ((I.EQ.2).AND.(J.EQ.3)) K=4
  219. D(I,J)=ENDT(K)
  220. D(J,I)=ENDT(K)
  221. 43 CONTINUE
  222. 42 CONTINUE
  223.  
  224. C----------------------------Calcul de la matrice tangente-------------
  225. DO 51 I=1,3
  226. DO 51 J=1,3
  227. DO 51 K=1,3
  228. DO 51 L=1,3
  229. C-------------------------
  230. U(I,J,K,L)=LAMB*DELT(I,J)*DELT(K,L)+MUT*(DELT(I,K)*DELT(J,L)
  231. * +DELT(I,L)*DELT(J,K))+ALP1*(DELT(I,J)*D(K,L)+D(I,J)*DELT(K,L))
  232. * +BETA*(D(I,K)*DELT(J,L)+D(L,J)*DELT(I,K)+DELT(I,L)
  233. * *D(J,K)+D(I,L)*DELT(J,K))
  234. C-------------------------
  235. 51 CONTINUE
  236. C-------------------------------Notation de VOIGT--------------------
  237. IF(MMI.EQ.1)GOTO 3
  238. DO 134 I=1,6
  239. M=MIII(I)
  240. N=NIII(I)
  241. DO 134 J=1,6
  242. K= MIII(J)
  243. L= NIII(J)
  244. VDEP(I,J)=U(M,N,K,L)
  245. 134 CONTINUE
  246. C-------------------------------------------------------------------
  247. DO 142 I=1,6
  248. DO 142 J=1,6
  249. VDEP1(I,J)=VDEP(I,J)
  250. VDEP2(I,J)=VDEP(I,J)
  251. 142 CONTINUE
  252. C********************************************************************
  253. C
  254. C********************************************************************
  255. 899 DO 321 I=1,6
  256. IF(ABS(SIGD(I)).LE.1E-15)SIGD(I)=0.D0
  257. 321 CONTINUE
  258. C--------------------------------------------------------------------
  259. CALL VECT(SIGD,VAL,VP,AAA(442),AAA(451),AAA(457))
  260. DO 332 I=4,6
  261. 332 EPSPLS(I)=EPSPLS(I)/2
  262. CALL VECT(EPSPLS,VAL2,VP2,AAA(442),AAA(451),AAA(457))
  263. DO 333 I=4,6
  264. 333 EPSPLS(I)=EPSPLS(I)*2
  265. CALL VECT(ENDD,VAL3,VP3,AAA(442),AAA(451),AAA(457))
  266. C--------------------------------------------------------------------
  267. CALL INVERE(VDEP2,6,6,KM,DET)
  268. C------------
  269. DO 401 I=1,6
  270. EPSTT(I)=0.D0
  271. DO 401 J=1,6
  272. EPSTT(I)=EPSTT(I)+VDEP2(I,J)*(-G*ENDD(J))
  273. 401 CONTINUE
  274. RETURN
  275. END
  276.  
  277.  
  278.  
  279.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales