Télécharger dtak.eso

Retour à la liste

Numérotation des lignes :

  1. C DTAK SOURCE CHAT 05/01/12 22:57:00 5004
  2. C SUBRUTINE DTAK : APPLICATION D'UNE LOI DE COMPORTEMENT POUR DES
  3. C DEFORMATIONS DUES A LA FLEXION
  4. C======================================================================
  5. SUBROUTINE DTAK(IFLAG,IDEP,ELEF1,ELDEP1,ELEF2,ELDEP2,
  6. 1 YLEF1,YLDEP1,YLEF2,YLDEP2,ULEF1,ULDEP1,
  7. 2 ULEF2,ULDEP2,EF,VARDEP,DEP,DEPE,
  8. 3 DEPMAX,DEPMIN,EFMAX,EFMIN,X0,RR)
  9. C
  10. C =====================================================================
  11. C ELEF1 - EFFORT MAX. ELASTIQUE DANS LE SENS +
  12. C ELDEP1 - DEFORMATION MAX. ELASTIQUE DANS LE SENS +
  13. C ELEF1 - EFFORT MAX. ELASTIQUE DANS LE SENS -
  14. C ELDEP1 - DEFORMATION MAX. ELASTIQUE DANS LE SENS -
  15. C YLEF1 - EFFORT MAX. LINEAIRE DANS LE SENS +
  16. C YLDEP1 - DEFORMATION MAX. LINEAIRE DANS LE SENS +
  17. C YLEF1 - EFFORT MAX. LINEAIRE DANS LE SENS -
  18. C YLDEP1 - DEFORMATION MAX.LINEAIRE DANS LE SENS -
  19. C ULEF1 - EFFORT MAX. PLASTIQUE DANS LE SENS +
  20. C ULDEP1 - DEFORMATION MAX. PLASTIQUE DANS LE SENS +
  21. C ULEF2 - EFFORT MAX. PLASTIQUE DANS LE SENS -
  22. C ULDEP2 - DEFORMATION MAX. PLASTIQUE DANS LE SENS -
  23. C EF - EFFORT DU PAS PRECEDENT
  24. C VARDEP - L'INCREMENT DE DEFORMATION
  25. C DEP - DEFORMATION DU PAS PRECEDENT
  26. C DEPMAX - DEFORMATION MAX. +
  27. C DEPMIN - DEFORMATION MAX. -
  28. C EFMAX - EFFORT MAX. +
  29. C EFMIN - EFFORT MAX. -
  30. C X0 - DERNIERE INTERSECTION AVEC L'AXE OX
  31. C IFLAG - FLAG POUR LE STADE DE COMPORTEMENT
  32. C REL - RIGIDITE ELASTIQUE COHERENTE AVEC LA MATRICE DE
  33. C RIGIDITE
  34. C REL1 - RIGIDITE ELASTIQUE DANS LE SENS +
  35. C REL2 - RIGIDITE ELASTIQUE DANS LE SENS -
  36. C RR1 - RIGIDITE SECANTE DANS LE SENS +
  37. C RR2 - RIGIDITE SECANTE DANS LE SENS -
  38. C RIN1 - RIGIDITE INELASTIQUE DANS LE SENS +
  39. C RIN2 - RIGIDITE INELASTIQUE DANS LE SENS -
  40. C DEPT - DEFORMATION TOTALE
  41. C RR - RIGIDITATE POUR LE STADE DE COMPORTEMENT ACTUEL
  42. C
  43. C ----------------------------------------------------------------------
  44. C
  45. IMPLICIT INTEGER(I-N)
  46. IMPLICIT REAL*8(A-H,O-Z)
  47. C
  48. C INITIALISATION
  49. C
  50. REL = (ELEF1 - ELEF2)/(ELDEP1 - ELDEP2)
  51. C
  52. IF (ELDEP1.LT.YLDEP1) THEN
  53. REL1 = (YLEF1 - ELEF1)/(YLDEP1 - ELDEP1)
  54. ELSE
  55. REL1 = REL
  56. END IF
  57. C
  58. IF (ELDEP2.GT.YLDEP2) THEN
  59. REL2 = (YLEF2 - ELEF2)/(YLDEP2 - ELDEP2)
  60. ELSE
  61. REL2 = REL
  62. END IF
  63. C
  64. RIN1 = (ULEF1 - YLEF1) / (ULDEP1 - YLDEP1)
  65. RIN2 = (ULEF2 - YLEF2) / (ULDEP2 - YLDEP2)
  66. RR1 = (YLEF1 / YLDEP1)
  67. RR2 = (YLEF2 / YLDEP2)
  68. C
  69. IF(IFLAG.EQ.0) THEN
  70. IFLAG = 1
  71. RR = REL
  72. DEPMAX = ELDEP1
  73. EFMAX = ELEF1
  74. DEPMIN = ELDEP2
  75. EFMIN = ELEF2
  76. DEP = 0.0
  77. X0 = 0.0
  78. END IF
  79. C
  80. DEPT = DEP + VARDEP
  81. C
  82. C CAS EF = 0
  83. C
  84. IF(EF.EQ.0.)THEN
  85. IF(IFLAG.EQ.1)THEN
  86. IF(DEPT.GT.0.)GO TO 11
  87. IF(DEPT.LT.0.)GO TO 12
  88. RETURN
  89. END IF
  90. IF((VARDEP.GT.0.).AND.(RR.EQ.(EFMIN / (DEPMIN - X0))))THEN
  91. RR = EFMAX / (DEPMAX - X0)
  92. GO TO 3
  93. ELSE IF((VARDEP.LT.0.).AND.(RR.EQ.(EFMAX/(DEPMAX-X0))))THEN
  94. RR = EFMIN / (DEPMIN - X0)
  95. GO TO 3
  96. END IF
  97. IF(IFLAG.EQ.3)GO TO 31
  98. END IF
  99. C
  100. GO TO(1,2,3,4),IFLAG
  101. C
  102. C STADE ELASTIQUE
  103. C
  104. 1 IFLAG = 1
  105. 11 IF(EF.GE.0)THEN
  106. IF(DEPT.GT.YLDEP1)THEN
  107. DEPMAX = YLDEP1
  108. EFMAX = YLEF1
  109. VARDEP = DEPT - YLDEP1
  110. DEP = YLDEP1
  111. EF = YLEF1
  112. RR = RIN1
  113. GO TO 2
  114. ELSE
  115. IF ((ELDEP1.LT.YLDEP1).AND.(DEPT.GE.0.))THEN
  116. IF (VARDEP.GE.0.)THEN
  117. IF ((DEP.GE.DEPMAX).OR.(DEPMAX.EQ.ELDEP1)) THEN
  118. IF ((DEP.LE.ELDEP1).AND.(DEPT.LE.ELDEP1))THEN
  119. RR=REL
  120. ELSE IF ((DEP.GE.ELDEP1).AND.(DEPT.GE.ELDEP1))THEN
  121. RR=REL1
  122. ELSE IF ((DEP.LE.ELDEP1).AND.(DEPT.GE.ELDEP1))THEN
  123. VARDEP=DEPT-ELDEP1
  124. EF=ELEF1
  125. DEP=ELDEP1
  126. RR=REL1
  127. ENDIF
  128. ELSE
  129. IF (DEPT.GE.DEPMAX) THEN
  130. VARDEP=DEPT-DEPMAX
  131. EF=EFMAX
  132. DEP=DEPMAX
  133. RR=REL1
  134. ELSE
  135. RR=(EFMAX/DEPMAX)
  136. ENDIF
  137. ENDIF
  138. ELSE
  139. IF (DEP.GT.DEPMAX) THEN
  140. DEPMAX=DEP
  141. EFMAX=EF
  142. ENDIF
  143. IF (DEPMAX.EQ.ELDEP1)THEN
  144. RR=REL
  145. ELSE
  146. RR=(EFMAX/DEPMAX)
  147. ENDIF
  148. ENDIF
  149. EF = EF + VARDEP * RR
  150. DEP = DEPT
  151. RETURN
  152. END IF
  153. END IF
  154. IF(DEPT.LT.0.) THEN
  155. DEP = 0.
  156. VARDEP = DEPT
  157. EF = 0.
  158. RR = REL
  159. GO TO 12
  160. END IF
  161. EF = EF + VARDEP * RR
  162. DEP = DEPT
  163. RETURN
  164. END IF
  165. 12 IF(EF.LE.0.)THEN
  166. IF(DEPT.LT.YLDEP2) THEN
  167. DEPMIN = YLDEP2
  168. EFMIN = YLEF2
  169. VARDEP = DEPT - YLDEP2
  170. DEP = YLDEP2
  171. EF = YLEF2
  172. RR = RIN2
  173. GO TO 2
  174. ELSE
  175. IF ((ELDEP2.GT.YLDEP2).AND.(DEPT.LE.0.))THEN
  176. IF (VARDEP.LE.0.)THEN
  177. IF ((DEP.LE.DEPMIN).OR.(DEPMIN.EQ.ELDEP2)) THEN
  178. IF ((DEP.GE.ELDEP2).AND.(DEPT.GE.ELDEP2))THEN
  179. RR=REL
  180. ELSE IF ((DEP.LE.ELDEP2).AND.(DEPT.LE.ELDEP2))THEN
  181. RR=REL2
  182. ELSE IF ((DEP.GE.ELDEP2).AND.(DEPT.LE.ELDEP2))THEN
  183. VARDEP=DEPT-ELDEP2
  184. EF=ELEF2
  185. DEP=ELDEP2
  186. RR=REL2
  187. ENDIF
  188. ELSE
  189. IF (DEPT.LE.DEPMIN) THEN
  190. VARDEP=DEPT-DEPMIN
  191. EF=EFMIN
  192. DEP=DEPMIN
  193. RR=REL2
  194. ELSE
  195. RR=(EFMIN/DEPMIN)
  196. ENDIF
  197. ENDIF
  198. ELSE
  199. IF (DEP.LT.DEPMIN)THEN
  200. DEPMIN=DEP
  201. EFMIN=EF
  202. ENDIF
  203. IF (DEPMIN.EQ.ELDEP2)THEN
  204. RR=REL
  205. ELSE
  206. RR=(EFMIN/DEPMIN)
  207. ENDIF
  208. ENDIF
  209. EF = EF + VARDEP * RR
  210. DEP = DEPT
  211. RETURN
  212. END IF
  213. END IF
  214. IF (DEPT.GT.0.) THEN
  215. DEP = 0.
  216. VARDEP = DEPT
  217. EF = 0.
  218. RR = REL
  219. GO TO 11
  220. END IF
  221. EF = EF + VARDEP * RR
  222. DEP = DEPT
  223. RETURN
  224. END IF
  225. C
  226. C STADE PLASTIQUE
  227. C
  228. 2 IFLAG = 2
  229. IF((EF * VARDEP).GE.0.) THEN
  230. EF = EF + VARDEP * RR
  231. DEP = DEPT
  232. IF(EF.GT.ULEF1) IDEP = 1
  233. IF(EF.LT.ULEF2) IDEP = -1
  234. RETURN
  235. END IF
  236. IF(EF.GE.0.)THEN
  237. DEPMAX = DEP
  238. EFMAX = EF
  239. RR = SQRT(YLDEP1 / DEP) * RR1
  240. ELSE
  241. DEPMIN = DEP
  242. EFMIN = EF
  243. RR = SQRT(YLDEP2 / DEPMIN) * RR2
  244. END IF
  245. C
  246. C DECHARGE SANS CHANGEMENT DE SENS
  247. C
  248. 3 IFLAG = 3
  249. IF(EF.EQ.0.)GO TO 4
  250. IF((EF * VARDEP).GT.0.)THEN
  251. IF(EF.GT.0.)RR = (EFMAX - EF) / (DEPMAX - DEP)
  252. IF(EF.LT.0.)RR = (EFMIN - EF) / (DEPMIN - DEP)
  253. GO TO 4
  254. END IF
  255. EFT1 = EF + VARDEP * RR
  256. IF((EFT1 * EF).GT.0.0)THEN
  257. EF = EFT1
  258. DEP = DEPT
  259. RETURN
  260. END IF
  261. 31 X0 = DEP - EF / RR
  262. IF(VARDEP.GT.0.)RR = EFMAX / (DEPMAX - X0)
  263. IF(VARDEP.LT.0.)RR = EFMIN / (DEPMIN - X0)
  264. DEP = X0
  265. VARDEP = DEPT - DEP
  266. EF = 0.
  267. IF(VARDEP.EQ.0.)RETURN
  268. C
  269. C DECHARGE APRES LE CHANGEMENT DE SENS
  270. C
  271. 4 IFLAG = 4
  272. IF((EF * VARDEP).GE.0.)THEN
  273. IF(DEPT.GT.DEPMAX)THEN
  274. VARDEP = DEPT - DEPMAX
  275. DEP = DEPMAX
  276. EF = EFMAX
  277. IF(DEPMAX.GT.YLDEP1)THEN
  278. RR = RIN1
  279. GO TO 2
  280. ELSE
  281. RR = (EFMAX/DEPMAX)
  282. GO TO 11
  283. END IF
  284. END IF
  285. IF(DEPT.LT.DEPMIN)THEN
  286. VARDEP = DEPT - DEPMIN
  287. DEP = DEPMIN
  288. EF = EFMIN
  289. IF(DEPMIN.LT.YLDEP2)THEN
  290. RR = RIN2
  291. GO TO 2
  292. ELSE
  293. RR = (EFMIN/DEPMIN)
  294. GO TO 11
  295. END IF
  296. END IF
  297. EF = EF + VARDEP * RR
  298. DEP = DEPT
  299. RETURN
  300. END IF
  301. IF(EF.GE.0.)RR = SQRT(YLDEP1 / DEPMAX) * RR1
  302. IF(EF.LT.0.)RR = SQRT(YLDEP2 / DEPMIN) * RR2
  303. GO TO 3
  304. END
  305.  
  306.  
  307.  

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