Télécharger la plaquette Cast3M


fr en

Télécharger ddnsth.eso

Retour à la liste

Numérotation des lignes : 

ddnsth

  1. C DDNSTH    SOURCE    CHAT      05/01/12    22:38:08     5004
  2. C DDNSTH    SOURCE    AM        95/12/27    21:26:13     1962
  3.       SUBROUTINE DDNSTH(IR,IC,DP,PN,U,E,E0,EYP,EYN,PCP,PCN,PYP,PYN,
  4.      *                  UYP,UYN,CP,CN,PSP,PSN,BP,BN,UEP,PEP,UEN,
  5.      *                  PEN,UCP,UCN,E1P,E1N,E2P,E2N,USP,USN,UTP,
  6.      *                  PTP,UTN,PTN,EP,EN)
  7. C=======================================================================
  8. C     3-PARAMETER MODEL
  9. C
  10. C     THIS PROGRAM TRACES THE HYSTERETIC BEHAVIOR OF R/C COMPONENTS
  11. C     - "DDNSTH" REPRESENTS DEGRADING, DETERIORATING,
  12. C        NONSYMMETRIC, SLIP-TYPE & TRILINEAR HYSTERETIC MODEL.
  13. C
  14. C    *ESSENTIAL VARIABLES
  15. C          IR= NUMBER OF BRANCH OF HYSTERETIC RULE
  16. C          IC= NON-ZERO VALUE INDICATES CHANGE OF BRANCH
  17. C          DP= INCREMENTAL STRESS
  18. C          PN= NEW STRESS
  19. C          U= NEW DEFORMATION
  20. C          E= INSTANTANIOUS STIFFNESS
  21. C    *PERMANENT VARIABLES FOR ENVELOPE CURVES (LAST LETTER "P" INDICATES
  22. C     POSITIVE SIDE, AND "N" INDICATES NEGATIVE SIDE)
  23. C          E0= INITIAL STIFFNESS
  24. C          EYP,EYN= POST-YIELDING STIFFNESS
  25. C          PCP,PCN= CRACKING STRESS
  26. C          PYP,PYN= YIELDING STRESS
  27. C          UYP,UYN= YIELDING DEFORMATION
  28. C
  29. C    *PERMANENT VARIABLES FOR DEGRADATION
  30. C          CP,CN= COEFFICIENT FOR COMMON POINTS OF UNLOADING BRANCH
  31. C
  32. C    *PERMANENT VARIABLES FOR SLIPPAGE
  33. C          PSP,PSN= SLIPPING YIELD STRESS (WHEN NO SLIP IS EXPECTED,
  34. C                   SET PSP>PYP,PSN<PYN)
  35. C
  36. C    *PERMANENT VARIABLES FOR DETERIORATION
  37. C          BP,BN= RATIO OF PLASTIC DEFORMATION INCREASE PER UNIT
  38. C                 INCREASE OF ABSORBED ENERGY.
  39. C
  40. C    *TEMPORAL VARIABLES FOR DEGRADATION
  41. C          UEP,PEP= EXPERIANCED MAXIMUM (POSITIVE) POINT
  42. C          UEN,PEN= EXPERIANCED MINIMUM (NEGATIVE) POINT
  43. C          UCP,UCN= UNLOADING (RELOADING) TURNING POINT ON X-AXIS
  44. C          E1P,E1N= FIRST UNLOADING (RELOADING) STIFFNESS
  45. C          E2P,E2N= SECOND RELOADING (UNLOADING)STIFFNESS
  46. C
  47. C    *TEMPORAL VARIABLES FOR SLIPPAGE
  48. C          USP,USN= TARGET SLIP YIELD DEF. (CRACK CLOSING POINT)
  49. C
  50. C     *TEMPORAL VARIABLES FOR DETERIORATION
  51. C          UPT,PTP= EXPECTED NEW MAXIMUN POINT
  52. C          UTN,PTN= EXPECTED NEW MINIMUM POINT
  53. C          EP     = INCREASE OF ABSORBED ENERGY TO CACULATE
  54. C                   NEW MAXIMUM POINT
  55. C          EN     = INCREASE OF ABSORBED ENERGY TO CACULATE
  56. C                   NEW MINIMUM POINT
  57. C
  58. C    *NECESSARY INITIAL CONDITION
  59. C          IR= 1 , PN= (LONG TERM LOADING) , PEN=PEP= 0
  60. C
  61. C===================================================================
  62.       IMPLICIT INTEGER(I-N)
  63.       IMPLICIT REAL*8(A-H,O-Z)
  64. C
  65.       IO=IR
  66.       UO=U
  67.       PO=PN
  68.       PN=PO+DP
  69.       GO TO(10,20,30,40,50,60,70,80,90,100,110,120,130,140,150),IR
  70. C
  71. C     -----
  72. C     IR=1
  73. C     -----
  74. C
  75. 10    IF(PN.LT.PCP)GO TO 11
  76.       IR=2
  77.       PO=PCP
  78.       UO=PCP/E0
  79.       E=(PYP-PCP)/(UYP-UO)
  80.       PN=PCP+(PN-PCP)*E/E0
  81.       GO TO 22
  82. 11    IF(PN.LT.PCN)GO TO 12
  83.       E=E0
  84.       U=PN/E
  85.       GO TO 200
  86. 12    IR=3
  87.       PO=PCN
  88.       UO=PCN/E0
  89.       E=(PYN-PCN)/(UYN-UO)
  90.       PN=PCN+(PN-PCN)*E/E0
  91.       GO TO 32
  92. C
  93. C     -----
  94. C     IR=2
  95. C     -----
  96. C
  97. 20    IF(DP.GE.0)GO TO 22
  98. 21    IR=6
  99.       PEP=PO
  100.       UEP=UO
  101.       EP=0.
  102.       E1P=(PO+CP*PYP)/(UO+CP*PYP/E0)
  103.       UCP=U-PO/E1P
  104.       USP=UCP
  105.       PN=PO+DP*E1P/E
  106.       E=E1P
  107.       GO TO 63
  108. 22    IF(PN.GT.PYP)GO TO 23
  109.       U=PCP/E0+(PN-PCP)/E
  110.       EN=EN+0.5*(PN+PO)*(U-UO)
  111.       GO TO 200
  112. 23    IR=4
  113.       PN=PYP+(PN-PYP)*EYP/E
  114.       U=UYP+(PN-PYP)/EYP
  115.       EN=EN+0.5*(PYP+PO)*(UYP-UO)+0.5*(PN+PYP)*(U-UYP)
  116.       E=EYP
  117.       GO TO 200
  118. C
  119. C     -----
  120. C     IR=3
  121. C     -----
  122. C
  123. 30    IF(DP.LE.0.)GO TO 32
  124. 31    IR=7
  125.       PEN=PO
  126.       UEN=UO
  127.       EN=0.
  128.       E1N=(PO+CN*PYN)/(UO+CN*PYN/E0)
  129.       UCN=UO-PO/E1N
  130.       USN=UCN
  131.       PN=PO+DP*E1N/E
  132.       E=E1N
  133.       GO TO 73
  134. 32    IF(PN.LT.PYN)GO TO 33
  135.       U=PCN/E0+(PN-PCN)/E
  136.       EP=EP+0.5*(PN+PO)*(U-UO)
  137.       GO TO 200
  138. 33    IR=5
  139.       PN=PYN+(PN-PYN)*EYN/E
  140.       U=UYN+(PN-PYN)/EYN
  141.       EP=EP+0.5*(PYN+PO)*(UYN-UO)+0.5*(PN+PYN)*(U-UYN)
  142.       E=EYN
  143.       GO TO 200
  144. C
  145. C     -----
  146. C     IR=4
  147. C     -----
  148. C
  149. 40    IF(DP.LT.0.)GO TO 21
  150. 41    U=UYP+(PN-PYP)/EYP
  151.       E=EYP
  152.       EN=EN+0.5*(PO+PN)*(U-UO)
  153.       GO TO 200
  154. C
  155. C     -----
  156. C     IR=5
  157. C     -----
  158. C
  159. 50    IF(DP.GT.0.)GO TO 31
  160. 51    U=UYN+(PN-PYN)/EYN
  161.       E=EYN
  162.       EP=EP+0.5*(PO+PN)*(U-UO)
  163.       GO TO 200
  164. C
  165. C     -----
  166. C     IR=6
  167. C     -----
  168. C
  169. 60    IF(PN.LE.PEP)GO TO 63
  170. 61    IF(PEP.GT.PYP)GO TO 62
  171.       IR=2
  172.       E2=(PYP-PCP)/(UYP-PCP/E0)
  173.       PN=PEP+(PN-PEP)*E2/E
  174.       EN=EN+0.5*(PEP+PO)*(UEP-UO)
  175.       PO=PEP
  176.       UO=UEP
  177.       E=E2
  178.       GO TO 22
  179. 62    IR=4
  180.       PN=PEP+(PN-PEP)*EYP/E
  181.       EN=EN+0.5*(PEP+PO)*(UEP-UO)
  182.       PO=PEP
  183.       UO=UEP
  184.       GO TO 41
  185. 63    IF(PN.LT.0.)GO TO 64
  186.       U=UCP+PN/E1P
  187.       X=0.5*(PO+PN)*(U-UO)
  188.       EP=EP+X
  189.       EN=EN+X
  190.       E=E1P
  191.       GO TO 200
  192. 64    IR=8
  193.       X=0.5*PO**2/E1P
  194.       PO=0.
  195.       UO=UCP
  196.       EP=EP-X
  197.       EN=EN-X
  198.       IF(PEN.LE.PCN)GO TO 65
  199.       E1N=E0
  200.       E2N=-PCN/(UCP-PCN/E0)
  201.       PN=PN*E2N/E1P
  202.       E=E2N
  203.       GO TO 84
  204. CAVP2     65    IF(PSN.GT.PYN)GO TO 66
  205. 65    IF(PSN.GT.PYN.AND.PSN.GT.PEN)GO TO 66
  206.       XE=0.5*PEN*(UEN-UCP)+EN
  207. CAVP1
  208.          IF(XE.LT.0.0)XE=0.0
  209. CAVP1
  210.       UTN=UEN-BN*XE
  211.       E2=(PYN-PCN)/(UYN-PCN/E0)
  212.       PTN=PCN+E2*(UTN-PCN/E0)
  213.       IF(UTN.LT.UYN)PTN=PYN+EYN*(UTN-UYN)
  214.       E1N=(CN*PYN+PTN)/(UTN+CN*PYN/E0)
  215.       E2N=PTN/(UTN-UCP)
  216.       PN=PN*E2N/E
  217.       E=E2N
  218.       GO TO 87
  219. 66    X=UEN-(PEN-PSN)/E1N-UCP
  220.       PX=E0*(PSN*UCP/X)/(PSN/X-E0)
  221.       XE=0.5*(UEN-UCP)*PX+0.5*PEN*(UEN-PX/E0)+EN
  222. CAVP1
  223.          IF(XE.LT.0.0)XE=0.0
  224. CAVP1
  225.       UTN=UEN-BN*XE
  226.       E2=(PYN-PCN)/(UYN-PCN/E0)
  227.       PTN=PCN+E2*(UTN-PCN/E0)
  228.       IF(UTN.LT.UYN)PTN=PYN+EYN*(UTN-UYN)
  229.       E1N=(CN*PYN+PTN)/(UTN+CN*PYN/E0)
  230.       USN=UTN-PTN/E1N
  231.       E2N=PSN/(UTN-(PTN-PSN)/E1N-UCP)
  232.       PN=PN*E2N/E
  233.       E=E2N
  234.       GO TO 89
  235. C
  236. C     -----
  237. C     IR=7
  238. C     -----
  239. C
  240. 70    IF(PN.GE.PEN)GO TO 73
  241. 71    IF(PEN.LT.PYN)GO TO 72
  242.       IR=3
  243.       E2=(PYN-PCN)/(UYN-PCN/E0)
  244.       PN=PEN+(PN-PEN)*E2/E
  245.       EP=EP+0.5*(PEN+PO)*(UEN-UO)
  246.       PO=PEN
  247.       UO=UEN
  248.       E=E2
  249.       GO TO 32
  250. 72    IR=5
  251.       PN=PEN+(PN-PEN)*EYN/E
  252.       EP=EP+0.5*(PEN+PO)*(UEN-UO)
  253.       PO=PEN
  254.       UO=UEN
  255.       GO TO 51
  256. 73    IF(PN.GT.0.)GO TO 74
  257.       U=UCN+PN/E1N
  258.       X=0.5*(PO+PN)*(U-UO)
  259.       EP=EP+X
  260.       EN=EN+X
  261.       E=E1N
  262.       GO TO 200
  263. 74    IR=9
  264.       X=0.5*PO**2/E1N
  265.       PO=0.
  266.       UO=UCN
  267.       EP=EP-X
  268.       EN=EN-X
  269.       IF(PEP.GT.1.001*PCP)GO TO 75
  270.       E1P=E0
  271.       E2P=PCP/(PCP/E0-UCN)
  272.       PN=PN*E2P/E1N
  273.       E=E2P
  274.       GO TO 94
  275. CAVP1   75    IF(PSP.LT.PYP)GO TO 76
  276. 75    IF(PSP.LT.PYP.AND.PSP.LT.PEP)GO TO 76
  277.       XE=0.5*PEP*(UEP-UCN)+EP
  278. CAVP1
  279.          IF(XE.LT.0.0)XE=0.0
  280. CAVP1
  281.       UTP=UEP+BP*XE
  282.       E2=(PYP-PCP)/(UYP-PCP/E0)
  283.       PTP=PCP+E2*(UTP-PCP/E0)
  284.       IF(UTP.GT.UYP)PTP=PYP+EYP*(UTP-UYP)
  285.       E1P=(CP*PYP+PTP)/(UTP+CP*PYP/E0)
  286.       E2P=PTP/(UTP-UCN)
  287.       PN=PN*E2P/E
  288.       E=E2P
  289.       GO TO 97
  290. 76    X=UEP-(PEP-PSP)/E1P-UCN
  291.       PX=E0*(PSP*UCN/X)/(PSP/X-E0)
  292.       XE=0.5*(UEP-UCN)*PX+0.5*PEP*(UEP-PX/E0)+EP
  293. CAVP1
  294.          IF(XE.LT.0.0)XE=0.0
  295. CAVP1
  296.       UTP=UEP+BP*XE
  297.       E2=(PYP-PCP)/(UYP-PCP/E0)
  298.       PTP=PCP+E2*(UTP-PCP/E0)
  299.       IF(UTP.GT.UYP)PTP=PYP+EYP*(UTP-UYP)
  300.       E1P=(CP*PYP+PTP)/(UTP+CP*PYP/E0)
  301.       USP=UTP-PTP/E1P
  302.       E2P=PSP/(UTP-(PTP-PSP)/E1P-UCN)
  303.       PN=PN*E2P/E
  304.       E=E2P
  305.       GO TO 99
  306. C
  307. C     -----
  308. C     IR=8
  309. C     -----
  310. C
  311. 80    IF(DP.LT.0.)GO TO 83
  312.       IR=13
  313. 81    UCN=UO-PO/E1N
  314.       PN=PO+DP*E1N/E2N
  315.       E=E1N
  316.       GO TO 73
  317. 83    IF(PEN.LT.1.001*PCN)GO TO 86
  318. 84    IF(PN.LT.PCN)GO TO 85
  319.       E=E2N
  320.       U=UCP+PN/E
  321.       X=0.5*(PO+PN)*(U-UO)
  322.       EP=EP+X
  323.       GO TO 200
  324. 85    IR=3
  325.       X=0.5*(PCN+PO)*(PCN/E0-UO)
  326.       EP=EP+X
  327.       E2=(PYN-PCN)/(UYN-PCN/E0)
  328.       PN=PCN+E2*(PN-PCN)/E2N
  329.       PO=PCN
  330.       UO=PCN/E0
  331.       E=E2
  332.       GO TO 32
  333. 86    IF(PSN.GT.PYN)GO TO 89
  334. 87    IF(PN.GT.PTN)GO TO 88
  335.       PEN=PTN
  336.       UEN=UTN
  337.       GO TO 71
  338. 88    E=E2N
  339.       U=UO+(PN-PO)/E
  340.       X=0.5*(PO+PN)*(U-UO)
  341.       EP=EP+X
  342.       EN=EN+X
  343.       GO TO 200
  344. CAVP89    PX=E0*E2N*UCP/(E2N-E0)
  345. CPP189    PX=E2N*(USN-UCP)
  346. CAVP89    PX=E2N*UCP*PTN/(E2N*UTN-PTN)
  347. 89    PX=E2N*(USN-UCP)
  348.       IF(PN.GT.PX)GO TO 88
  349.       IR=10
  350.       X=0.5*(USN-UO)*(PX+PO)
  351.       EP=EP+X
  352.       EN=EN+X
  353.       PO=PX
  354. CAVP      UO=PX/E0
  355.       UO=PX/E2N+UCP
  356.       E2N=(PTN-PX)/(UTN-UO)
  357.       PN=PX+(PN-PX)*E2N/E
  358.       E=E2N
  359.       GO TO 87
  360. C
  361. C     -----
  362. C     IR=9
  363. C     -----
  364. C
  365. 90    IF(DP.GT.0.)GO TO 93
  366.       IR=12
  367. 91    UCP=UO-PO/E1P
  368.       PN=PO+DP*E1P/E2P
  369.       E=E1P
  370.       GO TO 63
  371. 93    IF(PEP.GT.1.001*PCP)GO TO 96
  372. 94    IF(PN.GT.PCP)GO TO 95
  373.       E=E2P
  374.       U=UCN+PN/E
  375.       X=0.5*(PO+PN)*(U-UO)
  376.       EN=EN+X
  377.       GO TO 200
  378. 95    IR=2
  379.       X=0.5*(PCP+PO)*(PCP/E0-UO)
  380.       EN=EN+X
  381.       E2=(PYP-PCP)/(UYP-PCP/E0)
  382.       PN=PCP+E2*(PN-PCP)/E2P
  383.       PO=PCP
  384.       UO=PCP/E0
  385.       E=E2
  386.       GO TO 22
  387. 96    IF(PSP.LT.PYP)GO TO 99
  388. 97    IF(PN.LT.PTP)GO TO 98
  389.       PEP=PTP
  390.       UEP=UTP
  391.       GO TO 61
  392. 98    E=E2P
  393.       U=UO+(PN-PO)/E
  394.       X=0.5*(PO+PN)*(U-UO)
  395.       EP=EP+X
  396.       EN=EN+X
  397.       GO TO 200
  398. CAVP99    PX=E0*EP2*UCN/(E2P-E0)
  399. CAVP99    PX=E0*E2P*UCN/(E2P-E0)
  400. CPP199    PX=E2P*(USP-UCN)
  401. CAVP99    PX=E2P*UCN*PTP/(E2P*UTP-PTP)
  402. 99    PX=E2P*(USP-UCN)
  403.       IF(PN.LT.PX)GO TO 98
  404.       IR=11
  405.       X=0.5*(USP-UO)*(PX+PO)
  406.       EP=EP+X
  407.       EN=EN+X
  408.       PO=PX
  409. CAVP      UO=PX/E0
  410.       UO=PX/E2P+UCN
  411.       E2P=(PTP-PX)/(UTP-UO)
  412.       PN=PX+(PN-PX)*E2P/E
  413.       E=E2P
  414.       GO TO 97
  415. C
  416. C     -----
  417. C     IR=10
  418. C     -----
  419. C
  420. 100   IF(DP.LT.0.)GO TO 87
  421.       IR=15
  422.       GO TO 81
  423. C
  424. C     -----
  425. C     IR=11
  426. C     -----
  427. C
  428. 110   IF(DP.GT.0.)GO TO 97
  429.       IR=14
  430.       GO TO 91
  431. C
  432. C     -----
  433. C     IR=12
  434. C     -----
  435. C
  436. 120   PX=(UCP-UCN)/(E1P-E2P)*E1P*E2P
  437.       IF(PN.LT.PX)GO TO 63
  438.       IR=9
  439.       UX=UCP+PX/E1P
  440.       X=0.5*(PX+PO)*(UX-UO)
  441.       PN=PX+(PN-PX)*E2P/E1P
  442.       E=E2P
  443.       PO=PX
  444.       UO=UX
  445.       EP=EP+X
  446.       EN=EN+X
  447.       GO TO 93
  448. C
  449. C     -----
  450. C     IR=13
  451. C     -----
  452. C
  453. 130   PX=(UCN-UCP)/(E1N-E2N)*E1N*E2N
  454.       IF(PN.GT.PX)GO TO 73
  455.       IR=8
  456.       UX=UCN+PX/E1N
  457.       X=0.5*(PX+PO)*(UX-UO)
  458.       PN=PX+(PN-PX)*E2N/E1N
  459.       E=E2N
  460.       PO=PX
  461.       UO=UX
  462.       EP=EP+X
  463.       EN=EN+X
  464.       GO TO 83
  465. C
  466. C     -----
  467. C     IR=14
  468. C     -----
  469. C
  470. 140   PX=(E1P*E2P*(UCP-UTP)+E1P*PTP)/(E1P-E2P)
  471.       IF(PN.LT.PX)GO TO 63
  472.       IR=11
  473.       UX=UCP+PX/E1P
  474.       X=0.5*(PX+PO)*(UX-UO)
  475.       PN=PX+(PN-PX)*E2P/E1P
  476.       E=E2P
  477.       PO=PX
  478.       UO=UX
  479.       EP=EP+X
  480.       EN=EN+X
  481.       GO TO 97
  482. C
  483. C     -----
  484. C     IR=15
  485. C     -----
  486. C
  487. 150   PX=(E1N*E2N*(UCN-UTN)+E1N*PTN)/(E1N-E2N)
  488.       IF(PN.GT.PX)GO TO 73
  489.       IR=10
  490.       UX=UCN+PX/E1N
  491.       X=0.5*(PX+PO)*(UX-UO)
  492.       PN=PX+(PN-PX)*E2N/E1N
  493.       E=E2N
  494.       PO=PX
  495.       UO=UX
  496.       EP=EP+X
  497.       EN=EN+X
  498.       GO TO 87
  499. C     ------------------------------------------------------------------
  500. 200   IC=0
  501.       IF(IO.NE.IR)IC=IO
  502. C=======================================================================
  503.       RETURN
  504.       END
  505.  
  506.  
  507.  

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