Télécharger kelvin.eso

Retour à la liste

Numérotation des lignes :

  1. C KELVIN SOURCE BP208322 17/03/01 21:17:48 9325
  2. C KELVIN SOURCE MOUSTACAS 21/07/15 01:38:12 5004
  3.  
  4. SUBROUTINE KELVIN (wrk52,wrk53,wrk54,IB,IGAU,NBPGAU)
  5. C
  6. C==============================================================
  7. C chaine de KELVIN
  8. C==============================================================
  9. C----------------------------------------------------------------------
  10. C-ENTREES:
  11. C-------
  12. C-6 = NBR. DE COMPOSANTES DES CONTR. OU DES DEFORM.
  13. C-SIG0(6) = CONTR. AU DEBUT DU PAS D'INTEGRATION
  14. C-DEPST(6) = INCREMENT DES DEFORM. TOTALES
  15. C-NVARI = NBR. DE VARIABLES INTERNES
  16. C-VAR0(NVARI) = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  17. C
  18. C-NCOMAT = NBR. DE CARACTERISTIQUES MECANIQUES DU MATERIAU
  19. C-IVAL(NCOMAT) = INDICE DES COMPOSANTES DE MATERIAU
  20. C-XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  21. C-MFR = INDICE DE LA FORMULATION MECANIQUE
  22. C-ICARA = NBR. DE CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  23. C-XCAR(ICARA) = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
  24. C-IFOUR = DIMENSION DU MODÈLE 3D = 2
  25. C
  26.  
  27.  
  28. IMPLICIT INTEGER(I-N)
  29. IMPLICIT REAL*8 (A-H,O-Z)
  30.  
  31. -INC CCOPTIO
  32. -INC CCREEL
  33. -INC DECHE
  34.  
  35. PARAMETER (UNDEMI=.5D0,UN=1.D0,DEUX=2.D0,TROIS=3.D0)
  36. PARAMETER (QUATRE=4.D0)
  37. PARAMETER (HEURE= 24.D0,SECONDE=3600.D0)
  38. REAL*8 K,K1,K2,K3,TR1,TR2,TR3,ACN
  39. REAL*8 RMF1,RMF2,RMF3
  40. REAL*8 JPF1,JPF2,JPF3
  41. DIMENSION AN(6,6),BN(6,6)
  42. DIMENSION SR0(6,6),SR1(6,6),SR2(6,6),SR3(6,6)
  43. DIMENSION EF0(6,6),EF1(6,6),EF2(6,6),EF3(6,6)
  44. DIMENSION DEFTOT(6)
  45. DIMENSION SIGF(6)
  46. NVARI = VAR0(/1)
  47. c WRITE(*,*) 'NVARI', NVARI
  48.  
  49. C TESTER IFOUR POUR VÉRIFIER SI LE MODELE EST BIEN 3D
  50. IF (INPLAS.EQ.174) THEN
  51. c WRITE(*,*) 'INPLAS', INPLAS
  52. c WRITE(*,*) 'CALCUL FLUAGE CHAINE DE KELVIN'
  53. c WRITE(*,*) 'ENTRER TAU EN JOURS'
  54. ENDIF
  55.  
  56. c WRITE(*,*) 'MATEPL kelvin.eso ', MATEPL
  57. c write(*,*) ' inplas jnppla ',inplas,jnpla
  58. IF (IFOUR.NE.2) THEN
  59. c WRITE(*,*) 'IFOUR', IFOUR
  60. c WRITE(*,*) 'L OPTION DE CALCUL N EST PAS IMPLEMENTE'
  61. c WRITE(*,*) 'SIMULATION MASSIF 3D ISOTROPE SEULEMENT'
  62. KERRE = 99
  63. RETURN
  64. ENDIF
  65.  
  66. C----------------------------------------------------------------------
  67. C CARACTERISTIQUES
  68. C----------------------------------------------------------------------
  69.  
  70. KERRE=0
  71. C list xmax pour verifier la positon de rho et alpha
  72. epf1=epstf(1)
  73. epf2=epstf(2)
  74. epf3=epstf(3)
  75. epf4=epstf(4)
  76. epf5=epstf(5)
  77. epf6=epstf(6)
  78.  
  79. YOUN=XMAT(1)
  80. XNU=XMAT(2)
  81. RHO=XMAT(3)
  82. ALP=XMAT(4)
  83. YF1=XMAT(5)
  84. TF1=XMAT(6)
  85. YF2=XMAT(7)
  86. TF2=XMAT(8)
  87. YF3=XMAT(9)
  88. TF3=XMAT(10)
  89.  
  90. SSXX = SIG0(1)
  91. SSYY = SIG0(2)
  92. SSZZ = SIG0(3)
  93. SSXY = SIG0(4)
  94. SSXZ = SIG0(5)
  95. SSYZ = SIG0(6)
  96.  
  97. TINI = TEMP0
  98. TIFI = TEMPF
  99. DELT = (TEMPF-TEMP0)
  100.  
  101.  
  102. EPSFL01 = VAR0(1)
  103. EPSFL02 = VAR0(2)
  104. EPSFL03 = VAR0(3)
  105. EPSFL04 = VAR0(4)
  106. EPSFL05 = VAR0(5)
  107. EPSFL06 = VAR0(6)
  108.  
  109. SIGRE01 = VAR0(7)
  110. SIGRE02 = VAR0(8)
  111. SIGRE03 = VAR0(9)
  112. SIGRE04 = VAR0(10)
  113. SIGRE05 = VAR0(11)
  114. SIGRE06 = VAR0(12)
  115.  
  116. ELAS01 = VAR0(13)
  117. ELAS02 = VAR0(14)
  118. ELAS03 = VAR0(15)
  119. ELAS04 = VAR0(16)
  120. ELAS05 = VAR0(17)
  121. ELAS06 = VAR0(18)
  122.  
  123.  
  124. C-MODULE D YOUNG A LA FIN DE LA RELAXATION =
  125. C-SOMME DE TOUS LES RESSORT EN PARALLELE
  126. K = UN/((UN/YOUN)+(UN/YF1)+(UN/YF2)+(UN/YF3))
  127. C-MODULE D YOUNG RESSORT 1
  128. K1 = (YOUN*YF1)/(YOUN+YF1)
  129. K2 = (YOUN*YF2)/(YOUN+YF2)
  130. K3 = (YOUN*YF3)/(YOUN+YF3)
  131. C -(TEMPF)AU STEME PARALLELE 1
  132. TR1 = (YF1*TF1)/(YOUN+YF1)
  133. TR2 = (YF2*TF2)/(YOUN+YF2)
  134. TR3 = (YF3*TF3)/(YOUN+YF3)
  135.  
  136. ACN = UN/((UN+XNU)*(UN-(DEUX*XNU)))
  137.  
  138. C-CALCUL TRIDIMENSIONEL
  139. C-CREATION DE LA MATRICE DE RIGIDITE ELASTIQUE
  140.  
  141. DO I=1,6
  142. DO J=1,6
  143. AN(I,J)=XZERO
  144. BN(I,J)=XZERO
  145. SR0(I,J)=XZERO
  146. SR1(I,J)=XZERO
  147. SR2(I,J)=XZERO
  148. SR3(I,J)=XZERO
  149. EF0(I,J)=XZERO
  150. EF1(I,J)=XZERO
  151. EF2(I,J)=XZERO
  152. EF3(I,J)=XZERO
  153. ENDDO
  154. ENDDO
  155.  
  156. IF(IFOUR.EQ.2) THEN
  157. AN(1,1)=1.D0
  158. AN(1,2)=-XNU
  159. AN(1,3)=-XNU
  160. AN(2,1)=-XNU
  161. AN(2,2)=1.D0
  162. AN(2,3)=-XNU
  163. AN(3,1)=-XNU
  164. AN(3,2)=-XNU
  165. AN(3,3)=1.D0
  166. AN(4,4)=1.D0+XNU
  167. AN(5,5)=1.D0+XNU
  168. AN(6,6)=1.D0+XNU
  169.  
  170. BN(1,1)=1.D0-XNU
  171. BN(1,2)=XNU
  172. BN(1,3)=XNU
  173. BN(2,1)=XNU
  174. BN(2,2)=1.D0-XNU
  175. BN(2,3)=XNU
  176. BN(3,1)=XNU
  177. BN(3,2)=XNU
  178. BN(3,3)=1.D0-XNU
  179. BN(4,4)=(1.D0-(2.D0*XNU))/2.D0
  180. BN(5,5)=(1.D0-(2.D0*XNU))/2.D0
  181. BN(6,6)=(1.D0-(2.D0*XNU))/2.D0
  182. C MATRICE RIGIDITÉ SANS YOUNG
  183. ENDIF
  184.  
  185. DO I=1,6
  186. DO J=1,6
  187. AN(I,J) = AN(I,J)
  188. BN(I,J) = ACN*BN(I,J)
  189. ENDDO
  190. ENDDO
  191.  
  192. C===========RELAXATION : CALCUL DES CONTRAINTES=======================
  193. C=====================================================================
  194. RMF1 = (YOUN-K1)*EXP(-(TEMPF/TR1))
  195. RMF2 = (YOUN-K2)*EXP(-(TEMPF/TR2))
  196. RMF3 = (YOUN-K3)*EXP(-(TEMPF/TR3))
  197.  
  198. DO I=1,6
  199. DO J=1,6
  200. SR0(I,J) = BN(I,J)*K
  201. SR1(I,J) = BN(I,J)*RMF1
  202. SR2(I,J) = BN(I,J)*RMF2
  203. SR3(I,J) = BN(I,J)*RMF3
  204. ENDDO
  205. ENDDO
  206.  
  207. C======CAS ISOTROPE 3D
  208. C======RESSORT ELASTICITE INSTANTANNÉE
  209. S0XX = SR0(1,1)+SR0(2,1)+SR0(3,1)
  210. C WRITE(*,*) 'SOXX', S0XX
  211. S0YY = SR0(1,2)+SR0(2,2)+SR0(3,2)
  212. S0ZZ = SR0(1,3)+SR0(2,3)+SR0(3,3)
  213. S0XY = SR0(4,4)
  214. S0XZ = SR0(5,5)
  215. S0YZ = SR0(6,6)
  216. C======SYSTEME KV1
  217. S1XX = SR1(1,1)+SR1(2,1)+SR1(3,1)
  218. S1YY = SR1(1,2)+SR1(2,2)+SR1(3,2)
  219. S1ZZ = SR1(1,3)+SR1(2,3)+SR1(3,3)
  220. S1XY = SR1(4,4)
  221. S1XZ = SR1(5,5)
  222. S1YZ = SR1(6,6)
  223.  
  224. C======SYSTEME KV2
  225. S2XX = SR2(1,1)+SR2(2,1)+SR2(3,1)
  226. S2YY = SR2(1,2)+SR2(2,2)+SR2(3,2)
  227. S2ZZ = SR2(1,3)+SR2(2,3)+SR2(3,3)
  228. S2XY = SR2(4,4)
  229. S2XZ = SR2(5,5)
  230. S2YZ = SR2(6,6)
  231.  
  232. C======SYSTEME KV3
  233. S3XX = SR3(1,1)+SR3(2,1)+SR3(3,1)
  234. S3YY = SR3(1,2)+SR3(2,2)+SR3(3,2)
  235. S3ZZ = SR3(1,3)+SR3(2,3)+SR3(3,3)
  236. S3XY = SR3(4,4)
  237. S3XZ = SR3(5,5)
  238. S3YZ = SR3(6,6)
  239.  
  240. C======TOTAL DES CONTRAINTES
  241. SIGREF1 = (S0XX+S1XX+S2XX+S3XX)*epf1
  242. SIGREF2 = (S0YY+S1YY+S2YY+S3YY)*epf2
  243. SIGREF3 = (S0ZZ+S1ZZ+S2ZZ+S3ZZ)*epf3
  244. SIGREF4 = (S0XY+S1XY+S2XY+S3XY)*epf4
  245. SIGREF5 = (S0XZ+S1XZ+S2XZ+S3XZ)*epf5
  246. SIGREF6 = (S0YZ+S1YZ+S2YZ+S3YZ)*epf6
  247.  
  248. SIGF(1) = SIGREF1
  249. SIGF(2) = SIGREF2
  250. SIGF(3) = SIGREF3
  251. SIGF(4) = SIGREF4
  252. SIGF(5) = SIGREF5
  253. SIGF(6) = SIGREF6
  254.  
  255. c DO I=1,6
  256. c RELL=SIGF(I)
  257. c WRITE(6,*) ' I= ',I, 'SIGF = ',SIGF(I), 'RELL', RELL
  258. c ENDDO
  259.  
  260. C=============FLUAGE : CALCUL DES DÉFORMATIONS=========================
  261. C======================================================================
  262.  
  263. JPF1 = (UN/YF1)*(UN-EXP(-(TEMPF/TF1)))
  264. JPF2 = (UN/YF2)*(UN-EXP(-(TEMPF/TF2)))
  265. JPF3 = (UN/YF3)*(UN-EXP(-(TEMPF/TF3)))
  266.  
  267. DO I=1,6
  268. DO J=1,6
  269. EF0(I,J)=AN(I,J)/YOUN
  270. EF1(I,J)=AN(I,J)*JPF1
  271. EF2(I,J)=AN(I,J)*JPF2
  272. EF3(I,J)=AN(I,J)*JPF3
  273. ENDDO
  274. ENDDO
  275.  
  276. C======RESSORT ELASTIQUE
  277. E0XX = EF0(1,1)+EF0(2,1)+EF0(3,1)
  278. E0YY = EF0(1,2)+EF0(2,2)+EF0(3,2)
  279. E0ZZ = EF0(1,3)+EF0(2,3)+EF0(3,3)
  280. E0XY = EF0(4,4)
  281. E0XZ = EF0(5,5)
  282. E0YZ = EF0(6,6)
  283.  
  284. C======SYSTEME KV1
  285. E1XX = EF1(1,1)+EF1(2,1)+EF1(3,1)
  286. E1YY = EF1(1,2)+EF1(2,2)+EF1(3,2)
  287. E1ZZ = EF1(1,3)+EF1(2,3)+EF1(3,3)
  288. E1XY = EF1(4,4)
  289. E1XZ = EF1(5,5)
  290. E1YZ = EF1(6,6)
  291.  
  292. C======SYSTEME KV2
  293. E2XX = EF2(1,1)+EF2(2,1)+EF2(3,1)
  294. E2YY = EF2(1,2)+EF2(2,2)+EF2(3,2)
  295. E2ZZ = EF2(1,3)+EF2(2,3)+EF2(3,3)
  296. E2XY = EF2(4,4)
  297. E2XZ = EF2(5,5)
  298. E2YZ = EF2(6,6)
  299.  
  300. C======SYSTEME KV3
  301. E3XX = EF3(1,1)+EF3(2,1)+EF3(3,1)
  302. E3YY = EF3(1,2)+EF3(2,2)+EF3(3,2)
  303. E3ZZ = EF3(1,3)+EF3(2,3)+EF3(3,3)
  304. E3XY = EF3(4,4)
  305. E3XZ = EF3(5,5)
  306. E3YZ = EF3(6,6)
  307.  
  308. C======DÉFORMATIONS DUES AU FLUAGE
  309. EPSFLF1 = (E1XX+E2XX+E3XX)*SIGREF1
  310. EPSFLF2 = (E1YY+E2YY+E3YY)*SIGREF2
  311. EPSFLF3 = (E1ZZ+E2ZZ+E3ZZ)*SIGREF3
  312. EPSFLF4 = (E1XY+E2XY+E3XY)*SIGREF4
  313. EPSFLF5 = (E1XZ+E2XZ+E3XZ)*SIGREF5
  314. EPSFLF6 = (E1YZ+E2YZ+E3YZ)*SIGREF6
  315.  
  316. c======DEFORMATIONS ELASTIQUES
  317. ELAS1 = E0XX*SIGREF1
  318. ELAS2 = E0YY*SIGREF2
  319. ELAS3 = E0ZZ*SIGREF3
  320. ELAS4 = E0XY*SIGREF4
  321. ELAS5 = E0XZ*SIGREF5
  322. ELAS6 = E0YZ*SIGREF6
  323.  
  324. C======TOTAL DES DÉFORMATION : DEFORMATIONS DUES AU FLUAGE +DEFORMATIONS
  325. C======DUES A L'ELASTICITÉ INSTANTANNÉE DU RESSORT
  326. DEF1 = EPSFLF1+ ELAS1
  327. DEF2 = EPSFLF2+ ELAS2
  328. DEF3 = EPSFLF3+ ELAS3
  329. DEF4 = EPSFLF4+ ELAS4
  330. DEF5 = EPSFLF5+ ELAS5
  331. DEF6 = EPSFLF6+ ELAS6
  332.  
  333. DEFTOT(1) = DEF1
  334. DEFTOT(2) = DEF2
  335. DEFTOT(3) = DEF3
  336. DEFTOT(4) = DEF4
  337. DEFTOT(5) = DEF5
  338. DEFTOT(6) = DEF6
  339.  
  340. c DO I=1,6
  341. c IF(EPSTF(I).NE.0.D0) THEN
  342. c RELL=(DEFTOT(I)-EPSTF(I))/EPSTF(I)
  343. c ELSE
  344. c RELL=(DEFTOT(I)-EPSTF(I))
  345. c ENDIF
  346.  
  347. c WRITE(6,*) ' I= ',I, 'EPSTF = ',EPSTF(I),
  348. c & 'DEFTOT = ', DEFTOT(I), 'RELL', RELL
  349. c ENDDO
  350.  
  351. C===== VARIABLE INTERNE AU LA FIN DU PAS
  352. VARF(1)= EPSFLF1
  353. VARF(2)= EPSFLF2
  354. VARF(3)= EPSFLF3
  355. VARF(4)= EPSFLF4
  356. VARF(5)= EPSFLF5
  357. VARF(6)= EPSFLF6
  358. VARF(7)= SIGREF1
  359. VARF(8)= SIGREF2
  360. VARF(9)= SIGREF3
  361. VARF(10)= SIGREF4
  362. VARF(11)= SIGREF5
  363. VARF(12)= SIGREF6
  364. VARF(13)=ELAS1
  365. VARF(14)=ELAS2
  366. VARF(15)=ELAS3
  367. VARF(16)=ELAS4
  368. VARF(17)=ELAS5
  369. VARF(18)=ELAS6
  370.  
  371. END
  372.  
  373.  
  374.  
  375.  

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