Télécharger cmazzz.eso

Retour à la liste

Numérotation des lignes :

  1. C CMAZZZ SOURCE BP208322 17/03/01 21:15:49 9325
  2. SUBROUTINE CMAZZZ(WRK52,WRK53,WRK54,NVARI,Iecou)
  3. C MAZZZ SOURCE AM 98/12/23 21:39:16 3409
  4. C
  5. C calcule la deformation initiale et l'increment de deformation
  6. C a partir de la contrainte initiale et l'increment de contrainte
  7. C elastique puis appelle la subroutine MAZARS
  8. C
  9. C
  10. C variables en entree
  11. C
  12. C WRK0,KRK1,WRK5 pointeurs sur des segments de travail
  13. C
  14. C NSTRS1 nombre de composantes dans les vecteurs des contraintes
  15. C et les vecteurs des deformations
  16. C
  17. C NVARI nombre de variables internes (doit etre egal a 2)
  18. C
  19. C NMATT nombre de constantes du materiau
  20. C
  21. C ISTEP flag utilise pour separer les etapes dans un calcul non local
  22. C ISTEP=0 -----> calcul local
  23. C ISTEP=1 -----> calcul non local etape 1 on calcule les seuils
  24. C ISTEP=2 -----> calcul non local etape 2 on continue le calcul
  25. C a partir des seuils moyennes
  26. C
  27. C
  28. C variables en sortie
  29. C
  30. C VARF variables internes finales dans WRK0
  31. C
  32. C SIGF contraintes finales dans WRK0
  33. C
  34. C Modif L.Bode - 09/10/92 - Traitement particulier des coques
  35. C Modif L.Bode - 14/10/92 - Modifications complementaires
  36. C
  37. c KICH XCARB<- XCAR , CMAZAR <- MAZARS
  38.  
  39. IMPLICIT INTEGER(I-N)
  40. IMPLICIT REAL*8(A-H,O-Z)
  41. -INC CCOPTIO
  42. -INC DECHE
  43. *
  44. SEGMENT IECOU
  45. * COMMON/IECOU/NYOG,NYNU,NYALFA,NYSMAX,NYN,NYM,NYKK,
  46. INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7,
  47. C INTEGER NYOG, NYNU, NYALFA,NYSMAX,NYN, NYM, NYKK,
  48. 1 icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16,
  49. C . NYALF1,NYBET1,NYR, NYA, NYRHO,NSIGY, NNKX, NYKX, IND,
  50. 2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  51. C . NSOM, NINV, NINCMA,NCOMP, JELEM, LEGAUS,INAT, NCXMAT,
  52. 3 icow25,icow26,icow27,icow28,icow29,icow30,ICARA,
  53. C . LTRAC, MFR, IELE, NHRM, NBNN, NBELEM,ICARA,
  54. 4 icow32,icow33,NSTRS1,MFR1,icow36,icow37,icow38,
  55. C . LW2, NDEF, NSTRSS,MFR1, NBGMAT,NELMAT,MSOUPA,
  56. 5 icow39,icow40,icow41,icow42,icow43,icow44
  57. C . NUMAT1,LENDO, NBBB, NNVARI,KERR1, MELEME
  58. INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  59. . icow51,icow52,icow53,icow54,icow55,icow56
  60. . icow57,icow58
  61. ENDSEGMENT
  62.  
  63.  
  64. SEGMENT WRKK2
  65. REAL*8 EPSILI(NSTRSV)
  66. END SEGMENT
  67. *
  68. SEGMENT WRK6
  69. REAL*8 SIG0S(NSTRS),DEPSTS(NSTRS)
  70. END SEGMENT
  71. *
  72. INTEGER NVARI
  73. INTEGER KCAS,IRTD,ISTRS
  74. REAL*8 PREC,EPAI,FAC,AUX,AUX1,AUX2,YOUN,XNU
  75. REAL*8 UN
  76. PARAMETER (UN=1.D0)
  77. KERRE=0
  78. YOUN = XMAT(1)
  79. XNU = XMAT(2)
  80. C
  81. C ON CALCULE LES CONTRAINTES VRAIES DANS LE CAS DES COQUES
  82. C
  83. C print*,'dans mazzzz MFR1=', mfr1
  84. IF (MFR1 .EQ. 9) THEN
  85. EPAI=XCARB(1)
  86. SEGINI WRK6
  87. DO 101 ISTRS=1,NSTRS1
  88. SIG0S(ISTRS)=SIG0(ISTRS)
  89. DEPSTS(ISTRS)=DEPST(ISTRS)
  90. 101 CONTINUE
  91. NSTRSV=4
  92. IFOUR2=-2
  93. C Modif L.Bode - 14/10/92
  94. C Dans le cas des coques, on force la dimension a 2 pour MAZARS
  95. C ie on travaille en contraintes planes
  96. JDIM =2
  97. C Fin modif L.Bode
  98. DO 102 ISTRS=1,2
  99. SIG0(ISTRS)=SIG0S(ISTRS)/EPAI
  100. 102 CONTINUE
  101. DEPST(3)=0.D0
  102. DEPST(4)=DEPSTS(3)
  103. SIG0(3)=0.D0
  104. SIG0(4)=SIG0S(3)/EPAI
  105. ELSE IF (MFR1 .EQ. 1) THEN
  106. NSTRSV=NSTRS1
  107. IFOUR2=IFOUR
  108. C Modif L.Bode - 14/10/92
  109. C Pour les elts massifs, on utilise la vraie dimension
  110. JDIM = IDIM
  111. C Fin modif L.Bode
  112. ELSE
  113. PRINT*,'MFR1=',MFR1
  114. KERRE=57
  115. RETURN
  116. END IF
  117. SEGINI WRKK2
  118.  
  119.  
  120. C
  121. C calcul de la matrice elastique
  122. C
  123. CMATE = 'ISOTROPE'
  124. KCAS=1
  125. * print*,'increment de deformation elastique'
  126. CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,KCAS,DDHOOK,IRTD)
  127. * DO ISTRS=1,NSTRS1
  128. * print*,(DDHOOK(ISTRS,J),j=1,nstrs)
  129. * print*,DEPST(ISTRS)
  130. * END DO
  131. IF ( IRTD .EQ. 1) THEN
  132. C
  133. C calcul de l'increment de contrainte
  134. C
  135. CALL MATVE1 (DDHOOK,DEPST,NSTRSV,NSTRSV,DSIGT,1)
  136. C
  137. C
  138. C inversion de cette matrice
  139. C
  140. PREC=1.D-08
  141. CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,2,DDHOOK,IRTD)
  142. * DO ISTRS=1,NSTRSV
  143. * print*,(DDHOOK(ISTRS,J),j=1,nstrsv)
  144. * END DO
  145. * print*,'appel a invalm'
  146. CALL INVALM(DDHOOK,NSTRSV,NSTRSV,IRTD,PREC)
  147. IF (IRTD.EQ.0)THEN
  148. C
  149. C calcul des deformations du materiau elastique lineaire
  150. C
  151. * print*,'appel a matve1'
  152. CALL MATVE1 (DDHOOK,SIG0,NSTRSV,NSTRSV,EPSILI,1)
  153. C
  154. C modification pour tenir compte de l'endommagement
  155. C
  156. DO 100 ISTRS=1,NSTRSV
  157. EPSILI(ISTRS)=EPSILI(ISTRS)+EPIN0(ISTRS)
  158. 100 CONTINUE
  159.  
  160. C
  161. C appel a la routine MAZARS
  162. C
  163. C Modif L.Bode - 14/10/92
  164. C On envoie la dimension et le numero de la formulation
  165. C ( Elts Coques JDIM =2 , IFOUR2 = -2 => contraintes planes
  166. C Elts Massifs JDIM = IDIM ,IFOUR2 = IFOUR)
  167. * print*,'appel a mazars'
  168. icarbi=icara
  169. CALL CMAZAR (WRK52,WRK53,WRK54,WRKK2,NSTRSV,NVARI,
  170. 1 ICARbi,JDIM,IFOUR2)
  171. icara=icarbi
  172. C Fin modif L.Bode
  173. C
  174. C ON RECALCULE LES CONTRAINTES EFFECTIVES POUR LES COQUES
  175. C
  176. * print*,'apres mazars'
  177. IF (MFR1 .EQ. 9) THEN
  178.  
  179.  
  180. C
  181. C ON MET A JOUR DE FACON NON LINEAIRE LA PARTIE MEMBRANE
  182. C ET LES PARTIES FLEXION ET EFFORTS
  183. C TRANCHANTS LE CAS ECHEANT
  184. C
  185. FAC=(EPAI**3)/12.D0
  186. AUX =FAC*YOUN/(1.D0-XNU*XNU)
  187. AUX1=FAC*YOUN*.5D0/(1.D0+XNU)
  188. AUX2=EPAI*YOUN*.5D0/(1.d0+XNU)/1.2d0
  189. DO 103 ISTRS=1,NSTRS1
  190. SIG0(ISTRS)=SIG0S(ISTRS)
  191. DEPST(ISTRS)=DEPSTS(ISTRS)
  192. 103 CONTINUE
  193.  
  194. DO 104 ISTRS=1,2
  195. SIGF (ISTRS)=SIGF(ISTRS)*EPAI
  196. 104 CONTINUE
  197. SIGF(3)=SIGF(4)*EPAI
  198. SIGF(4)=SIG0(4)+AUX*(DEPST(4)+XNU*DEPST(5))
  199. SIGF(5)=SIG0(5)+AUX*(DEPST(5)+XNU*DEPST(4))
  200. SIGF(6)=SIG0(6)+AUX1*DEPST(6)
  201. SIGF(7)=SIG0(7)+AUX2*DEPST(7)
  202. SIGF(8)=SIG0(8)+AUX2*DEPST(8)
  203. SEGSUP WRK6
  204. END IF
  205. ELSE
  206. print*,'erreur dans invalm'
  207. KERRE=56
  208. END IF
  209. ELSE
  210. print*,'erreur dans dohmas'
  211. KERRE=56
  212. END IF
  213. SEGSUP WRKK2
  214. RETURN
  215. END
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  

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