Professional Documents
Culture Documents
Design Guidance For High Temperature Concentrating Solar Power Components
Design Guidance For High Temperature Concentrating Solar Power Components
DOCUMENT AVAILABILITY
Online Access: U.S. Department of Energy (DOE) reports produced after 1991
and a growing number of pre-1991 documents are available free at OSTI.GOV
(http://www.osti.gov/), a service of the US Dept. of Energy’s Office of Scientific and
Technical Information.
Reports not in digital format are available to DOE and DOE contractors
from the Office of Scientific and Technical Information (OSTI):
U.S. Department of Energy
Office of Scientific and Technical Information
P.O. Box 62
Oak Ridge, TN 37831-0062
www.osti.gov
Phone: (865) 576-8401
Fax: (865) 576-5728
Email: [email protected]
Disclaimer
This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United
States Government nor any agency thereof, nor UChicago Argonne, LLC, nor any of their employees or officers, makes any
warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any
information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights.
Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise,
does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any
agency thereof. The views and opinions of document authors expressed herein do not necessarily state or reflect those of the
United States Government or any agency thereof, Argonne National Laboratory, or UChicago Argonne, LLC.
ANL-20/03
prepared by
Bipul Barua, Argonne National Laboratory
Michael D. McMurtrey, Idaho National Laboratory
Ryann E. Rupp, Idaho National Laboratory
Mark C. Messner, Argonne National Laboratory
January 2020
Abstract
This report provides guidance for the design of components for concentrating power facilities
operating at high temperatures and undergoing high, secondary thermal stresses relative to the
applied primary pressure stress. The design rules were developed for the design of Generation 3
CSP tubular receivers manufactured from Alloy 740H, but are generally applicable to a wide range
of component types undergoing similar loads. Part 1 of the report provides procedural design rules
for components to be used in conjunction with the 2019 edition of the ASME Boiler & Pressure
Vessel Code. Part 2 provides corresponding design data for Alloy 740H. Part 3 of the report is a
commentary describing the rationale behind the design rules and the data underlying the design
material properties. Finally, Part 4 provides an extensive set of worked sample problems detailing
the application of the rules to CSP components.
Table of Contents
Abstract ............................................................................................................................ iii
Table of Contents ............................................................................................................... v
Part 1: Design Rules ........................................................................................................... 1
1-1. General criteria .............................................................................................................. 1
1-2. Design Option A: Design by Elastic Analysis using ASME Section III, Division 5 ........ 2
1-3. Design Option B: Design by Elastic Analysis using ASME Section III, Division 5 with
Reduced Margin and Simplified Creep-Fatigue evaluation .................................................... 4
1-4. Design Option C: Design by Simple Inelastic Analysis .................................................. 6
Part 2: Design Material Data ............................................................................................... 8
2-1. General Criteria ............................................................................................................. 8
2-2. Material Specification .................................................................................................... 8
2-3. Elastic constants ............................................................................................................ 9
2-4. Thermal properties ......................................................................................................... 9
2-5. Mechanical properties .................................................................................................. 10
2-6. Allowable Stresses ....................................................................................................... 12
2-7. Creep, fatigue, and creep-fatigue properties ................................................................. 14
2-8. Weld strength reduction factors.................................................................................... 17
2-9. Isochronous stress-strain relations ................................................................................ 18
2-10. Inelastic constitutive models ...................................................................................... 25
2-11. Temperature limits ..................................................................................................... 26
Part 3: Commentary .......................................................................................................... 27
3-1. Overview ..................................................................................................................... 27
3-2. Commentary on the Design Rules ................................................................................ 29
3-3. Commentary on the 740H Material Data ...................................................................... 42
Part 4: Sample problems ................................................................................................... 58
4-1. Sample problem 1 ........................................................................................................ 58
4-2. Sample problem 2 ........................................................................................................ 80
Acknowledgements ........................................................................................................ 103
Part 1: Design Rules
1-1. General criteria
1-1.1. Applicability
These design criteria apply to components in concentrating solar power facilities at temperatures
above 370° C for ferritic and ferritic-martensitic steels and 425° C for austenitic stainless steels
and nickel-based alloys where creep-fatigue damage in cyclic service or stress relaxation damage
caused by reoccurring application of secondary load is a significant design consideration. The
designer may select any of the following three design options listed in Articles 1-2, 1-3, and 1-4.
A design is required to pass all the checks contained in the selected option to pass the design
criteria.
The design criteria were developed for structures undergoing daily cycling. The rules below are
not suited for structures that see service cycles with holds at constant load longer than 1,000 hours
at temperatures greater than those listed in Part 2, Article 2-10.1.
All references to the ASME Boiler and Pressure Vessel Code (ASME B&PV Code) are to the 2019
Edition.
All three methods require the designer to use a single design composite loading cycle to represent
or bound the service conditions experienced in operation. The definition of this Design Cycle
consists of periodic pressure, thermal, and mechanical force boundary conditions sufficient to
complete a thermal-mechanical analysis of the component along with the number of times this
composite loading cycle will be repeated in service – defined as quantity 𝑁. These boundary
condition histories must include times of application including all relevant hold periods at constant
load. The composite cycle period is defined as 𝑡# .
The Design Specification shall specify a Design Life denoted as 𝑡$%&'() . The composite loading
cycle used for design must cover the entire design life so that 𝑡$%&'() = 𝑁𝑡# .
The following design criteria are to be used in conjunction with the design material data provided
in Part 2.
1-1.4. Limitations
These criteria only cover the design of components and do not provide construction, welding,
examination, or inspection criteria. These additional requirements shall be provided in the Design
Specification, possibly by reference to ASME Section III, Section V, Section VIII, Section IX, or
Section XI or other codes and standards documents, as appropriate.
ANL-20/03 1
1-2. Design Option A: Design by Elastic Analysis using ASME Section
III, Division 5
This option can be applied to any type of component using any material provided in Part 2 up to
the maximum metal temperatures provided in the tabulated design data.
Primary load analysis is to be completed using a steady-state thermal analysis of the Design Cycle
and a linear elastic stress analysis using the material properties provided in Articles 0 and 2-3 of
Part 2.
The component shall meet the criteria of the ASME B&PV Code Section III, Division 5, HBB-
3222.1 with the following modifications:
1. The Design Loading shall be determined using Section III, Division 5, HBB-3113.1 expect
references to “Service Level A Loadings” are replaced by references to the “Design Cycle”
defined in Article 1-1.2.
2. The values of the allowable stress 𝑆- are provided in Article 2-5.1 of Part 2.
3. The provisions of HBB-3222.1(c) do not apply.
The temperatures used in the ratcheting analysis are to be obtained from a transient thermal
analysis of the component subject to the Design Cycle. The stresses and strains are to be obtained
using a small-deformation linear elastic analysis. Material properties are provided in Articles 0
and 2-3 of Part 2.
The component shall meet the criteria of the ASME B&PV Code Section III, Division 5, HBB-T-
1332 Test B-1, including the General Requirements of HBB-T-1331, with the following
modifications:
1. The designer will only analyze a single cycle, the Design Cycle defined in Article 1-1.2.
References to the Section III, Division 5 Service Loadings are replaced by references to
the Design Cycle.
2. The values of 𝑆. are provided in Article 2-4.1 of Part 2.
3. The temperature limit provided by Table HBB-1323 is instead provided in Article 2-10.2
in Part 2.
4. For HBB-T-1332(b) calculate the accumulated strain using a stress of 1.0𝜎3 , rather than
the 1.25𝜎3 in the ASME Code.
5. For HBB-T-1332(b) the isochronous stress strain curves are provided in Article 0 of Part
2.
6. The strain limits in HBB-T-1332(b) are increased to 2% for base metal and 1% for weld
metal.
2 ANL-20/03
1-2.3. Creep-fatigue Criteria
The temperatures used in the creep-fatigue analysis are to be obtained from a transient thermal
analysis of the component subject to the Design Cycle. The stresses and strains are to be obtained
using a small-deformation linear elastic analysis. Material properties are provided in Articles 0
and 2-3 of Part 2.
The component shall meet the criteria contained in the ASME B&PV Code Section III, Division
5, HBB-T-1430 with the following modifications:
1. The designer will only analyze a single cycle, the Design Cycle defined in Article 1-1.2.
References to the Section III, Division 5 Service Loadings are replaced by references to
the Design Cycle.
2. The design fatigue curves are provided in Article 2-6.2 of Part 2.
3. The minimum stress to rupture data are provided in Article 2-6.1 in Part 2.
4. The values of the allowable stress 𝑆6 , required for calculating the quantity 3𝑆6 ̅ , are
provided in Article 0 in Part 2.
5. The values of the relaxation strength 𝑆9: and 𝑆9; , required for calculating the quantity 3𝑆6̅ ,
are provided in Article 0 in Part 2.
6. The isochronous stress-strain curves are provided in Article 0 in Part 2.
7. For HBB-T-1432 the designer shall use HBB-T-1413 for calculating the effective strain
range and not any of the other options allowed in HBB-T-1432(a).
8. The creep-fatigue damage envelope is provided in Article 0 in Part 2.
9. When calculating the creep strain increment in HBB-T-1432(g) use a stress equal to 1.0𝜎3 ,
rather than 1.25𝜎3 as in Section III, Division 5.
10. When using HBB-T-1433(a) Step 5(b) to evaluate creep damage the value of 𝑆:< shall be
1.0𝜎3 , rather than 1.25𝜎3 as in Section III, Division 5.
11. The alternate creep damage calculation procedure in HBB-T-1433(b) shall not be used.
The designer must use the process defined in HBB-T-1433(a).
12. Creep damage shall be calculated from the stress history using a factor 𝐾 > = 0.9.
13. The alternative criteria in HBB-T-1435 shall not be used.
14. The additional factors on weld material properties defined in HBB-T-1710 shall be applied
to weldments. Applicable weld strength reduction factors are provided in Article 2-7 of
Part 2.
For this analysis the Design Cycle loads, defined in Article 1-1.2, shall be increased by a factor of
1.5. This includes thermal stresses which can be increased by factoring the material coefficient of
thermal expansion.
Additionally, the Design Cycle shall be supplemented by intermittent loads, such as wind lateral
loads, that are not periodic. The method for calculating lateral loads shall be specified in the
Design Specification. If this method includes load factors accounting for uncertainty in the applied
loading these load factors shall be used for the lateral loads if they exceed the generic load factor
of 1.5 provided in this Article. If the lateral loads specified in the Design Specification are not
ANL-20/03 3
factored or if the load factor is less than 1.5, the intermittent loads shall be increased by a factor of
1.5.
The factored Design Cycle plus factored intermittent loads shall be used to perform a time-
independent, large deformations, incremental plastic analysis of the component using a von Mises
flow theory and a temperature-dependent flow stress defined by the zero-time, hot-tensile
isochronous stress strain curves given in Part 2, Article 0. If the analysis converges for one
application of the combined load history, i.e. the structure does not undergo plastic collapse or
buckling, the time independent buckling design criteria are satisfied.
For many CSP systems time dependent buckling will not be a significant design issue as primary
stresses and steady, load-controlled lateral loads are typically small. However, if a design has high
primary loads or significant steady, load controlled loadings, such as large lateral self-weight, that
could cause buckling, time dependent buckling may be assessed with the following procedure.
Note that such steady load controlled forces should already be included in the Design Cycle, as
they are primary loads.
For this analysis the Design Cycle loads, defined in Article 1-1.2, shall be increased by a factor of
1.5. This includes thermal stresses, which can be increased by factoring the material coefficient
of thermal expansion.
This factored Design Cycle shall be used to perform a time-independent, large deformations,
incremental plastic analysis of the component using a von Mises flow theory and a temperature-
dependent flow stress defined by isochronous curve for a time equal to the component Design Life
given in Part 2, Article 0. If the analysis converges for one application of the factored Design
Cycle, i.e. the structure does not undergo plastic collapse or buckling, the time dependent buckling
design criteria are satisfied.
This method shall only be applied to Alloy 740H material and to components where the peak stress
is minimal.
4 ANL-20/03
The temperatures used in the creep-fatigue analysis are to be obtained from a transient thermal
analysis of the component subject to the Design Cycle. The stresses and strains are to be obtained
using a small-deformation linear elastic analysis. Material properties are provided in Articles 0
and 2-3 of Part 2.
Consider the total stress intensity (P+Q+F) from the linear elastic analysis. If at any point during
the Design Cycle the total stress intensity exceeds the temperature-dependent value of 𝑆. for the
material, given in Part 2, Article 2-4.1 then Design Option B cannot be used to assess the
component.
To evaluate creep-fatigue damage use the following procedure. The stresses used in these steps
correspond to the maximum total stress (P+Q+F) in the linear elastic analysis.
1. Determine the maximum elastic strain range during the Design Cycle using the elastically
calculated mechanical strains from the analysis and the definition of the strain range given
in HBB-T-1413. Call this strain range Δ𝜀B .
2. Calculate a strain range increase due to creep Δ𝜀C by taking the total creep strain used to
evaluate the ratcheting strain criteria in Article 1-3.2 (using HBB-T-1332 Test B-1),
labeled 𝜀D-DEF and dividing it by the number of repetitions of the design cycle 𝑁, i.e. Δ𝜀C =
𝜀D-DEF /𝑁. As the structure will already have passed the requirements of Article 1-3.2 the
designer may alternatively conservatively use Δ𝜀C = 2%/𝑁.
3. Sum the two strain ranges Δ𝜀 = Δ𝜀B + Δ𝜀C and use this total strain range to determine the
number of cycles to failure using the design fatigue diagrams given in Part 2, Article 2-6.2.
The appropriate diagram corresponds to the maximum metal temperature in the Design
Cycle. Use the design fatigue curve to determine the number of allowable cycles 𝑁J .
Calculate the fatigue damage fraction using the equation 𝐷J = 𝑁/𝑁J .
4. Calculate the von Mises stress history, 𝜎L6 = 𝜎L6 =
C C
MNO𝜎PP − 𝜎.. R + O𝜎.. − 𝜎SS R + (𝜎SS − 𝜎PP )C + 6O𝜎P.
C + 𝜎 C + 𝜎 ^2RX /2
.S PS
corresponding to the total stress history determined in the linear elastic analysis. Use this
von Mises stress history to calculate the creep damage fraction for a single repetition of the
Design Cycle using HBB-T-1411(10) with a factor 𝐾 > = 0.9. That is, the scalar stress used
to determine the time-to-rupture for each time increment Δ𝑡 is the von Mises stress, not the
effective stress given in HBB-T-1411 or the design by elastic analysis procedure given in
HBB-T-1433. The minimum stress to rupture tables are provided in Part 2, Article 2-6.1.
Call this single cycle creep damage fraction 𝐷3B . The total creep damage fraction is 𝐷3 =
𝑁𝐷3B .
5. Use the creep-fatigue damage envelope provided in Part 2, Article 0 with the fatigue and
creep damage fractions 𝐷J and 𝐷3 to determine the acceptability of the component.
In this procedure the modifications to the weld material properties defined in HBB-T-1710 shall
be applied to weldments. Applicable weld strength reduction factors are provided in Article 2-7
of Part 2.
ANL-20/03 5
1-3.4. Time Independent Buckling
This method may be applied to any material for which Part 2, Article 2-9 provides a constitutive
model.
Perform a transient thermal analysis of the component under the Design Cycle using the properties
in Part 2, Article 2-3. Use these temperatures to perform a small deformation inelastic stress
analysis of the component under the Design Cycle loads using the inelastic constitutive model
defined in Part 2, Article 2-9. For each location in the component, 𝒙, this process produces a time
history of stress values 𝝈(𝒙, 𝑡), mechanical strain values 𝜺(𝒙, 𝑡), and temperature values 𝑇(𝒙, 𝑡).
The analysis must repeat the Design Cycle loading until the structure achieves a steady-state
response. This steady-state response is defined by the conditions, for some arbitrary time 𝑡:
1. 𝝈(𝒙, 𝑡) = 𝝈O𝒙, 𝑡 − 𝑡# R
2. 𝑇(𝒙, 𝑡) = 𝑇O𝒙, 𝑡 − 𝑡# R
for all times in the current repetition of the Design Cycle. Both conditions must be met. Once the
analysis reaches this steady-state condition no additional repetitions of the Design Cycle are
required. Extract the stress, strain, and temperature histories from the analysis corresponding to
the final repetitions of the design cycle. Shift these histories in time so that the final stress, strain,
and temperature histories start at time 0 and extend to time 𝑡# . These histories, 𝝈(𝒙, 𝑡), 𝜺(𝒙, 𝑡),
and 𝑇(𝒙, 𝑡), will be used to determine the acceptability of the component.
To assess the component against the ratcheting criteria for each location in the component calculate
the local ratcheting rate
2 2
𝑟(𝒙) = _ 𝜺O𝒙, 𝑡# R: 𝜺O𝒙, 𝑡# R − _ 𝜺(𝒙, 0): 𝜺(𝒙, 0)
3 3
which is the net effective strain accumulated over one cycle in the steady state condition. The ‘:’
in this equation represents tensor contraction. For each point calculate the total ratcheting strain:
𝑅(𝒙) = 𝑁 𝑟(𝒙)
This total ratcheting strain must be less than 10% for base metal and 5% for weldments at all points
in the structure. If this criteria is met the structure passes the ratcheting design criteria.
6 ANL-20/03
1-4.3. Creep-fatigue Criteria
The analysis is identical to the analysis required for the ratcheting check defined in Article 1-4.2.
The design analysis begins with the time-shifted histories 𝝈(𝒙, 𝑡), 𝜺(𝒙, 𝑡), and 𝑇(𝒙, 𝑡).
1. Use HBB-T-1413 with 𝜈 ∗ = 0.5 to calculate the effective strain range Δ𝜖6EP from the
strain history 𝜺(𝒙, 𝑡).
2. Use the design fatigue curves provided in Part 2, Article 2-6.2 to determine the maximum
allowable number of cycles 𝑁J (𝒙). Use the curve corresponding to the maximum metal
temperature at the point under consideration.
e
3. Calculate the fatigue damage fraction using the equation 𝐷J (𝒙) = e (𝒙).
f
4. Use the stress history 𝝈(𝒙, 𝑡) to calculate the von Mises effective stress history at each
point in the structure using the equation 𝜎L6 =
C C
MNO𝜎PP − 𝜎.. R + O𝜎.. − 𝜎SS R + (𝜎SS − 𝜎PP )C + 6O𝜎P.
C + 𝜎 C + 𝜎 ^2RX /2.
.S PS
5. Use the resulting von Mises stress history 𝜎L6 (𝒙, 𝑡)/𝐾′, with 𝐾 > = 0.9, and the local metal
temperature to calculate a time history of allowable time to rupture 𝑇$ (𝒙, 𝑡) using the
minimum stress to rupture provided in Part 2, Article 2-6.1.
D $D
6. Calculate the creep damage fraction using the equation 𝐷3 (𝒙) = 𝑁 ∫l k i (𝒙,D). The integral
j
may be discretized into individual time steps.
7. Use the creep-fatigue damage envelope, provided in Part 2, Article 0, to determine whether
each point in the structure passes or fails the acceptance criteria using the fatigue and creep
damage fractions 𝐷J (𝒙) and 𝐷3 (𝒙).
The structure passes the creep-fatigue design provisions if each point individually passes the
criteria described in 1-7 above.
ANL-20/03 7
Part 2: Design Material Data
2-1. General Criteria
Depending on the property, the data below is provided as tables, charts, or equations or some
combination thereof. If more than one form of the data is presented the designer may use either
the equations or interpolation from table or chart data. For the case where the material property
depends on time or a number of cycle repetitions log-linear interpolation shall be used for that
variable. Because of differences in interpolation versus the mathematical formula and because the
tabulated values are rounded the three sources of data may not precisely agree.
8 ANL-20/03
Elastic constants
The design Poisson’s ratio for Alloy 740H is 0.31 for all temperatures.
ANL-20/03 9
2-4. Mechanical properties
Temperature 𝑆.
(°C) (MPa)
40 621
100 594
150 577
200 562
250 548
300 538
350 531
400 529
450 529
500 529
550 529
600 529
650 529
700 529
750 508
800 463
850 418
Table 2-4.1.1. Design values of yield strength (𝑆. ) for Alloy 740H.
10 ANL-20/03
2-4.2. Tensile strength
Temperature 𝑆m
(°C) (MPa)
40 1034
100 1034
150 1034
200 1030
250 998
300 976
350 967
400 966
450 966
500 966
550 966
600 957
650 921
700 860
750 771
800 651
850 531
Table 2-4.2.1. Design values of tensile strength (𝑆m ) for Alloy 740H.
ANL-20/03 11
2-5. Allowable Stresses
Temperature 𝑆-
(°C) (MPa)
40 295
100 295
150 295
200 279
250 276
300 276
350 276
400 276
450 276
500 276
550 276
600 274
650 226
700 146
750 84.1
800 34.5
850 21.8
Table 2-5.1.1. Allowable stress 𝑆- in MPa for Alloy 740H.
12 ANL-20/03
2-5.2. Allowable Stress 𝑺𝒎
Temperature 𝑆6
(°C) (MPa)
40 345
100 345
150 345
200 343
250 333
300 325
350 322
400 322
450 322
500 322
550 322
600 319
650 307
700 287
750 257
800 217
850 177
Table 2-5.2.1. Allowable stress 𝑆6 in MPa for Alloy 740H.
ANL-20/03 13
2-5.3. Relaxation Strength
Time (hours)
1 10 30 100 300 1000 3000 10000
Temp. 425 483 483 483 483 483 483 483 483
(°C) 450 483 483 483 483 483 483 483 483
475 483 483 483 483 483 483 483 483
500 483 483 483 483 483 483 483 483
525 483 483 483 483 483 483 483 483
550 483 483 483 483 483 483 483 483
575 481 481 481 481 481 481 481 481
600 478 478 478 478 478 476 472 459
625 469 469 469 468 466 460 445 412
650 460 460 459 456 449 428 394 345
675 445 444 441 432 412 374 328 278
700 430 425 417 395 360 310 265 221
725 406 395 377 340 296 248 208 172
750 382 358 326 280 237 195 162 132
775 348 308 270 224 186 151 124 101
800 312 256 217 176 145 116 95 76
825 272 207 171 137 111 88 71 57
850 231 164 133 105 84 66 53 42
Table 2-5.3.1. Relaxation strength as a function of time and temperature for Alloy 740H. Values are in MPa.
The minimum stress to rupture for Alloy 740 is given by the equation
where the temperature 𝑇 is given in units of degrees Celsius, the time 𝑡 is given in units of hours,
and the rupture stress is provided in units of MPa.
14 ANL-20/03
Time (hours)
1 10 30 100 300 1000 3000 10000 30000 100000
Temp. 425 966 966 966 966 966 966 966 966 966 966
(°C) 450 966 966 966 966 966 966 966 966 966 966
475 966 966 966 966 966 966 966 966 966 966
500 966 966 966 966 966 966 966 966 966 966
525 966 966 966 966 966 966 966 966 966 923
550 966 966 966 966 966 966 966 966 864 736
575 962 962 962 962 962 962 950 805 692 586
600 957 957 957 957 957 898 768 647 554 467
625 939 939 939 939 869 729 621 521 444 372
650 921 921 921 836 709 592 502 419 355 296
675 891 891 824 685 578 480 405 337 284 236
700 860 808 679 561 472 390 328 271 228 188
725 816 668 559 460 385 316 265 218 182 150
750 771 553 460 377 314 257 214 175 146 119
775 677 457 379 309 256 208 173 141 117 95
800 565 378 312 253 209 169 140 113 93 76
825 472 313 257 207 170 137 113 91 75 60
850 394 259 212 170 139 111 91 73 60 48
Table 2-6.1.1. Design rupture stresses (𝑆9 ) in MPa for Alloy 740H.
The design fatigue relation for Alloy 740H is defined by the equation:
ANL-20/03 15
Cycles Strain range Strain range
𝑇 ≤ 700℃ 700℃ < 𝑇 ≤ 850℃
(mm/mm) (mm/mm)
10 0.018733 0.014665
20 0.015615 0.012790
40 0.013220 0.011154
100 0.010855 0.009308
200 0.009504 0.008118
400 0.008427 0.007080
1000 0.007236 0.005908
2000 0.006336 0.005153
4000 0.005618 0.004494
10000 0.004875 0.003750
20000 0.004427 0.003271
40000 0.004051 0.002852
100000 0.003639 0.002380
Table 2-7.2.1. Design fatigue relation for Alloy 740H.
16 ANL-20/03
2-6.3. Creep-fatigue damage envelopes
Figure 2-6.3.1 provides the creep-fatigue damage envelope for Alloy 740H. Regions below the
curve are acceptable.
Table 2-7.1 lists weld strength reduction factor corresponding to each allowable Alloy 740H weld
type:
ANL-20/03 17
2-8. Isochronous stress-strain relations
The isochronous stress-strain relations for Alloy 740H are described by the following equations
where 𝜀 is the total strain, 𝜎 the stress, 𝑇 temperature, and 𝑡 time:
𝜀 = 𝜀% + 𝜀# + 𝜀3
𝜎
𝜀% =
𝐸
0 𝜎 ≤ 𝜎l
⎧
⎪ › 𝜎 − 𝜎l ) 600℃ ≤ 𝑇 ≤ 800℃
⎪ 𝐾œ • 𝜎 > 𝜎l
𝜎l
𝜀# = 0 𝜎 ≤ 𝜎B
⎨
⎪›−1 ln Ÿ1 − 𝜎 − 𝜎B
⎪
𝜎 > 𝜎B
𝑇 > 800℃
⎩ 𝛿 𝜎# − 𝜎B
<£¤¥ 𝜎 ‰£¤¥
𝜀3 = 𝜀-̇ 𝑒 ¦§i ( ) ¦§i 𝑡
𝜇
Article 2-2.1.1 provides the Young’s modulus, 𝐸. Tables 2-8.1.1 and 2-8.1.2 lists the parameters
for the equations. For Table 2-8.1.1 parameters shall be interpolated linearly between temperatures
in the table. The units for temperature are ℃, stress is in MPa, time in hours, and strains in
mm/mm. The relations shall not be used for temperatures outside those provided in Table 2-8.1.1,
i.e. the relations are only valid for 600℃ ≤ 𝑇 ≤ 850℃. In these equations the shear modulus 𝜇
shall be calculated with the formula
𝐸
𝜇=
2(1 + 𝜈 )
with 𝐸 the Young’s modulus defined in Article 2-2.1.1 and 𝜈 the Poisson’s ratio defined in Article
2-2.2.1.
18 ANL-20/03
Ramberg-Osgood model parameters Voce hardening model parameters
Temperatures
𝜎l(MPa) K n 𝜎# (MPa) 𝜎B(MPa) 𝛿
600°C - 700°C 400.24 0.0704 6.6480
725°C 374.20 0.0357 7.1315
Parameter Value
𝜀-̇ 1.19x1010 hr-1
k 1.38064x10-20 mJ/K
b 2.53x10-07 mm
A -10.98557
B -0.53098
Table 2-8.1.2. Parameters for the 𝜀3 contribution to the Alloy 740H isochronous stress-strain relation.
Figures 2-8.1.1 to 2-8.1.11 plot the isochronous stress-strain relations for Alloy 740H.
ANL-20/03 19
Figure 2-8.1.1. Design isochronous stress-strain curve for Alloy 740H at 600º C.
Figure 2-8.1.2. Design isochronous stress-strain curve for Alloy 740H at 625º C.
20 ANL-20/03
Figure 2-8.1.3. Design isochronous stress-strain curve for Alloy 740H at 650º C.
Figure 2-8.1.4. Design isochronous stress-strain curve for Alloy 740H at 675º C.
ANL-20/03 21
Figure 2-8.1.5. Design isochronous stress-strain curve for Alloy 740H at 700º C.
Figure 2-8.1.6. Design isochronous stress-strain curve for Alloy 740H at 725º C.
22 ANL-20/03
Figure 2-8.1.7. Design isochronous stress-strain curve for Alloy 740H at 750º C.
Figure 2-8.1.8. Design isochronous stress-strain curve for Alloy 740H at 775º C.
ANL-20/03 23
Figure 2-8.1.9. Design isochronous stress-strain curve for Alloy 740H at 800º C.
Figure 2-8.1.10. Design isochronous stress-strain curve for Alloy 740H at 825º C.
24 ANL-20/03
Figure 2-8.1.11. Design isochronous stress-strain curve for Alloy 740H at 850º C.
The follow inelastic constitutive model may be used for Alloy 740H for 600℃ ≤ 𝑇 ≤ 850℃:
𝝈̇ = 𝑪: (𝜺̇ − 𝜺̇ D − 𝜺3̇ )
with 𝝈̇ the rate of the Cauchy stress tensor, 𝐶 the isotropic elasticity tensor constructed using the
temperature-dependent Young’s modulus provided in Article 2-2.1.1 and Poisson’s ratio provided
in Article 2-2.2.1, 𝜺̇ is the total strain rate, 𝜺̇ D is the thermal strain rate provided by the equation
𝜺̇ D = 𝛼𝑇̇𝑰
with 𝛼 the instantaneous coefficient of thermal expansion provided in Article 2-3.1, 𝑇̇ the
temperature rate, and 𝑰 the identity tensor, and 𝜺̇ 3 is the creep strain rate provided by the equation
<£¤¥ 𝜎L ‰£¤¥ 𝒔
𝜺3̇ = 𝜀-̇ 𝑒 ¦§i ( ) ¦§i
𝜇 𝜎L
ANL-20/03 25
®
where Article 2-8.1 defines the parameters in this equation, 𝜎L = MC 𝒔: 𝒔, and 𝒔 is the deviatoric
stress tensor.
The temperature at which Alloy 740H accumulates 0.1% strain at a stress of 𝑆- over 100,000 hours
life is 595° C.
The cutoff temperature for the B-1 and B-2 tests in ASME Section III, Division 5, HBB-T-1332
is 600° C for Alloy 740.
26 ANL-20/03
Part 3: Commentary
3-1. Overview
The purpose of the design methods and material data contained in this report is to enable the design
of components for solar receivers that operate in conditions where secondary stresses can cause
significant creep or creep-fatigue damage. In the context of solar systems, secondary stresses are
typically thermal stresses and primary stresses are typically caused by pressure. Therefore, these
design criteria are shaped towards structures that see low pressures but high thermal gradients,
operate at high temperatures in the creep regime, and see daily load cycling. These conditions
point towards eventual component failure via creep or creep-fatigue damage accumulated during
daily application and relaxation of the thermal stress, rather than creep failure caused by long times
at constant pressure loads.
Commercial practice, for example Sections I and VIII of the ASME Boiler and Pressure Vessel
Code, guards against creep failure under primary loads by including creep rupture properties in the
high temperature allowable stresses. This approach is reasonable for components that see long
periods of hold time at fixed conditions and limited numbers of load cycles. It is less reasonable
for components which see large numbers of load cycles, like CSP systems, as it does not account
for the damaging effect of secondary load. Figure 3-1.1 illustrates the differences between steady
state load conditions, traditional high temperature nuclear systems assuming a moderate amount
of load cycling interspersed with relatively long hold periods, and CSP systems by illustrating
representative stress/time histories
Figure 3-1.1. Illustration of four different types of loading histories: (1) Steady-state, non-cyclic loading where the
component stresses relax to the primary load. (2) Cyclic load where the secondary stresses do not cause the
component to reset to the initial value of load. (3) Cyclic load where the secondary stresses cause the component to
reset to the initial, high level of stress each cycle. (4) A prototypical CSP load cycle where the loads increase and
decrease gradually throughout the cycle period. Conventional Section I and Section VIII design assumes a load of
type (1), Section III, Division 5 high temperature nuclear design assumes either (2) or (3) with criteria to determine
if the global relaxation mode occurs. This report provides design rules suitable for load type (4).
ANL-20/03 27
The key challenge to address in developing an adequate high temperature CSP design method is
to guard against creep and creep-fatigue damage caused by cyclic secondary stresses. In the
context of a tubular receiver these stresses would be the daily application of the thermal stress
caused by both the through-wall and circumferential thermal gradients. However, the design
criteria also provide:
1. allowable stress criteria for steady-state primary loads to maintain compatibility with
current commercial practice and to provide a minimum section thickness;
2. ratcheting design criteria to prevent service failures caused by excess deformation;
3. buckling criteria, to prevent both time-independent elastic-plastic buckling and time-
dependent buckling caused by creep deformation.
However, for the sample designs considered in Part IV it is creep damage caused by the daily
application of secondary thermal stresses that dictates the maximum design life of the component.
The design rules below lump this damage mechanism in with creep-fatigue design. Indeed, as
Figure 3-1.2 illustrates, creep-relaxation damage caused by secondary loads is one end of the
spectrum of loadings encompassed by the creep-fatigue design procedures. On the other end of
the spectrum is pure fatigue damage and in between lies creep-fatigue interaction. If the
components sampled in Part IV operated at somewhat lower metal temperatures then creep-fatigue
interaction would control the component life. If they operated at much lower temperatures then
pure fatigue damage would control the life. The creep-fatigue interaction diagram therefore
encompasses the full range of failure mechanism for low pressure, high thermal stress components.
Figure 3-1.2. Illustration of how the creep-fatigue checks used here encompass failure by pure fatigue, failure by
pure creep, and failure by creep-fatigue interaction. The region of the design envelope an efficient component
occupies depends on the relative severity of the primary and secondary load, the type of load cycle, and the average
metal temperature.
The design rules in Part I are given procedurally and the component must pass all the required
design checks to be considered adequate. However, for the types of CSP components envisioned
when creating these design processes a more practical process would be to first find the minimum
component thickness using the primary load, allowable stress design procedures. If the primary
stress is low, as for the sample problems in Part IV, this minimum thickness might be too thin to
28 ANL-20/03
feasibly manufacture. If so, this minimum thickness should be increased to the thinnest
manufacturable section.
Once this thinnest section has been determined the designer should evaluate creep-fatigue damage
using one of the options presented in Part I. If the component passes the creep-fatigue criteria
using the minimum section thickness determined above then the designer can proceed to evaluating
the remaining criteria (ratcheting and buckling). These will only rarely constrain the design – most
often the design thickness will be determined by the tradeoff between the primary load allowable
stresses and the secondary creep-fatigue checks.
Figure 3-1.3 illustrates this tradeoff and demonstrates that if the design does not pass the creep-
fatigue check for the minimum section thickness determined by the primary load design then the
component will need to be reconfigured – increasing the section thickness will only increase the
secondary load on the component and make it more difficult to pass the creep-fatigue criteria. At
high metal temperatures CSP components will be constrained by the minimum section thickness
provided by the primary load design criteria and the maximum section thickness provided by the
creep-fatigue design criteria.
Figure 3-1.3. Illustration of the tradeoff between primary load criteria setting a minimum section thickness and
creep-fatigue criteria setting a maximum section thickness. The green shaded region is that feasible design space –
the thicknesses that meet or exceed the target design life. Note that the minimum primary load section thickness
provides the most economic design, hence the design strategy suggested here.
The ratcheting and buckling criteria represent important potential failure modes but rarely
constrain the design of CSP components. The design should check these criteria only after they
achieve a feasible design considering the primary load and creep-fatigue design provisions.
The high temperature/low temperature cutoff is standard ASME practice – 370° C for ferritic and
ferritic-martensitic steels and 425° C for austenitic stainless steels and nickel-based alloys. CSP
ANL-20/03 29
components are likely to spend most or all of their service life in the creep regime and so these
rules focus exclusively on high temperature applications.
The design provisions here are focused on CSP structures undergoing daily cycling. They may
not be adequate for structures with long hold times. The threshold temperature was based on the
temperature at which the material accumulates 0.1% strain at the allowable stress 𝑆- over 100,000
hours. This is essentially a negligible creep criteria.
Several of the Design Options are modifications of the Section III, Division 5 ASME Boiler and
Pressure Vessel Code rules covering the design of high temperature nuclear reactor components.
These modifications were developed starting from the 2019 edition of the ASME Code.
To simplify the design process these rules work with a single, representative Design Cycle. This
is a collection of periodic structural and thermal boundary conditions that adequately represent the
operating cycle of the component. For a CSP system this would likely be a daily cycle that
reasonably represents the standard operation of the facility. An example for a solar receiver could
be the service conditions on the solar equinox. Section III, Division 5 of the ASME Code contains
rules for evaluating several different types of loading cycles with different expected frequencies.
The current rules could be extended to cover multiple cycle types at the expense of a more
complicated design process.
The design methods here use a similar database of material properties to the ASME rules, but alter
the design margin encoded in the ASME data to reflect the relative consequences of failure
between CSP and nuclear systems. As such, the design rules should be used with the provided
material data, not the information in Section III of the ASME Code.
The rules here only cover component design. A complete design and construction code would
require rules for construction, welding, examination, and inspection. The Design Specification
could provide these criteria by reference to the ASME Code. Structures passing these design rules
should also pass the Section VIII, Division 2 design by analysis criteria. These provisions could
then be used to supplement the design of a Code-stamped component to ensure the component will
perform well under cyclic conditions at high temperatures – conditions not currently considered in
detail by either Section I or Section VIII.
This design option is based on the Section III, Division 5 ASME Boiler and Pressure Vessel Code
covering the design and construction of high temperature nuclear reactors. The design methods
contained in Section III, Division 5 are well documented1,2 and so this commentary only provides
a general overview of the approach. The focus of the commentary is on the modifications made
to the base ASME method to account for the unconventional operating cycles and temperatures of
1 R. I. Jetter, “Subsection NH—Class 1 Components in Elevated Temperature Service” in Companion Guide to the
ASME Boiler and Pressure Vessel Code, Volume 1, K. R. Rao ed. ASME Press, New York, NY, 2009.
2 M. H. Jawad, and R. I. Jetter, Design & Analysis of ASME Boiler and Pressure Vessel Components in the Creep
Range. ASME Press, New York, NY, 2009.
30 ANL-20/03
CSP systems as well as modifications made to account for the lower consequences of failure when
comparing nuclear reactors to CSP facilities.
The Section III, Division 5 checks mirror those used in these design rules. The Code mandates a
primary load check against allowable stresses plus checks for ratcheting, creep-fatigue, and
buckling. The methods used in Design Option A are those from the ASME Code simplified,
modified, and with reduced design margin.
Section III, Division 5 actually contains two primary load design checks: one mirroring the Section
I / Section VIII allowable stress checks where the primary stress is compared to some temperature-
dependent but time-independent allowable stress 𝑆- . This check uses a design temperature and
pressure that bound the standard, non-faulted service conditions. The Code also has second
primary stress check against a time- and temperature-dependent allowable stress 𝑆6D . This check
uses the Service Loads provided by the Owner’s Design Specification. The nuclear code accounts
for different expected frequencies and severities of loading by providing A, B, C, and D Service
Load categories. The time-dependent primary load check uses different margins depending on the
load category. Finally, Section III, Division 5 requires not only that each individual Service Load
pass the allowable stress check but also that the time-fraction of all the service loadings combined
remains less than one.
As discussed above, these rules for CSP components operate only with a single Design Cycle.
This is a compromise between Section I and VIII practice, which only typically use Design
Conditions and do not account explicitly for cyclic service, and the Section III, Division 5 practice
of evaluating the component against all expected transients, classifying each particular transient
based on its severity and expected frequency.
The rules retain a check against a time-independent allowable stress using the primary stress from
a set of Design Conditions, now determined as the bounding temperature, pressure, and nozzle
loads from the Design Cycle. This is analogous to current Section I and VIII commercial practice.
The design rules omit the Section III, Division 5 check against time-dependent allowable stress
and the use-fraction check for several reasons:
1. This retains the direct connection to current commercial practice and does not add the
complication of an additional time-dependent allowable stress check.
2. As discussed above, the primary load check for CSP systems in unlikely to be the
controlling factor. The main point of the primary load check in high temperature CSP
systems is to provide a reasonable minimum section thickness, which can be done without
invoking time-dependent allowables.
3. The allowable stress 𝑆- is based on 100,000 hour extrapolated creep properties. This is not
far removed from 10-30 year design lives for actual components.
The determination of the 𝑆- allowable stress values is discussed below in 3-3.7. This process is
essentially consistent across Section III, Division 5 and Section I/Section VIII, Division 2 and so
a reasonable modification of the design rules could be to use Section I or Section VIII rules to first
ANL-20/03 31
design the component and then apply the ratcheting, creep-fatigue, and buckling checks provided
in these rules. This would be consistent with the complete approach described above and yet also
ensure the component passes the Section I or Section VIII criteria.
The Section III, Division 5 allowable stress 𝑆6D has two parts: a time-independent allowable stress
𝑆6 and a time-dependent allowable stress 𝑆D . 𝑆6D is simply the lesser of the two for a given time
and temperature. While the design rules use neither stress as an allowable stress, the Section III,
Division 5 ratcheting and creep-fatigue rules use the values of 𝑆6 to bound the time-independent
initial stresses experienced by a component. As such, Design Option A requires values of 𝑆6 to
complete the ratcheting and creep-fatigue design checks. Section 3-3.8 describes how these values
are obtained.
HBB-3222.1(c) requires a buckling check per the Section III, Division 5 rules. The present design
rules provide an alternative method for assessing buckling failure.
The design rules use a simplified version of the B-1 test from Section III, Division 5, Subsection
HB, Subpart B, Nonmandatory Appendix T. The B-1 test is an implementation of the O’Donnell-
Porowski3 approach. This approach in turn builds on the classical work of Bree4 and so the
ratcheting test is most accurate for cylindrical thin-walled pressure vessels. However, it can be
conservatively used for other types of structures where the peak stress is negligible. If the structure
has significant peak stress Design Option C can be used.
Essentially, the design rules retain the Section III, Division 5 approach and simply remove
additional margin included in the nuclear code. The ASME implementation of the method
increases the core stress calculated by O’Donnell-Porowski by 25%. The design rules here do not
increase the core stress. As such, this method uses a best-estimate of the accumulated strain, rather
than the bounding values used in Section III, Division 5. Additionally, the design rules increase
the allowable accumulated stains from 1% for base metal and ½ % for welds to 2% for base metal
and 1% for welds. This reflects the reduced consequences of failure for CSP components and is
in line with previous CSP design guidance.
3 W. O’Donnell, and J. Porowski, “Upper bounds for accumulated strains due to creep ratcheting” Journal of Pressure
Vessel Technology, 96(3), pp. 150-154, 1974.
4 J. Bree, “Elastic-plastic behavior of thin tubes subjected to internal pressure and intermittent high-heat fluxes with
application to fast-nuclear-reactor fuel elements” Journal of Strain Analysis for Engineering Design, 2(3), pp. 226-
238, 1967.
32 ANL-20/03
Creep-fatigue interaction describes extensive experimental data showing that the combination of
cyclic strain controlled deformation interspersed with holds at constant strain is more damaging
than either holding at constant strain or cycling with the same strain range but no hold periods. A
variety of microstructural mechanisms have been posited to explain the interaction5,6,7, but for the
purposes of developing design rules it is sufficient to note the empirical phenomenon.
The process of going from a full 3D strain state to a scalar strain range requires the definition of
an effective strain range. These rules adopt the ASME definition.
There is an ongoing disagreement as to how to quantify creep damage at constant stress or constant
strain. Several methods have been proposed and are in use in current design and fitness for service
standards. These techniques include time fraction, used in the ASME Boiler and Pressure Vessel
Code, ductility exhaustion8, used in the RCC-MRx and R5 Codes, and the Omega method9 used
in API-579/ASME FFS-1. These different methods will all lead to somewhat different predictions
for rupture given the same stress/strain/time history. However, studies have shown overall they
are all about equally accurate10,11. These design rules adopt the ASME time-fraction approach.
Time fraction is a very straightforward method for calculating a creep damage fraction. It needs a
scalar effective stress measure – a map from a 3D state of stress to a scalar stress – and a database
of the material rupture stress as a function of time and temperature. The rupture stress information
5 S. W. Nam, S. C. Lee, and J. M. Lee, “The effect of creep cavitation on the fatigue life under creep-fatigue
interaction” Nuclear Engineering and Design, 153(2–3), pp. 213–221, 1995.
6 J. Wareing, “Creep-fatigue interaction in austenitic stainless steels” Metallurgical Transactions A, 8(5), pp. 711–
721, 1977
7 M. Sauzay, M. Mottot, L. Allais, M. Noblecourt, I. Monnet, and J. Périnet, “Creep-fatigue behaviour of an AISI
stainless steel at 550°C” Nuclear Engineering and Design, 232(3), pp. 219–236, 2004.
8 R. M. Goldhoff, “Uniaxial Creep-Rupture Behavior of Low-Alloy Steel Under Variable Loading Conditions”
Journal of Basic Engineering, 87(2), pp. 374–378, 1965.
9 M. Prager, “Development of the MPC Omega Method for Life Assessment in the Creep Range” Journal of Pressure
Vessel Technology, 117(2), pp. 95–103, 1995.
10 S.-L. Mannan, F.-Z. Xuan, X.-C. Zhang, Y.-C. Lin, S.-T. Tu, and X.-L. Yan, “Review of creep–fatigue endurance
and life prediction of 316 stainless steels” International Journal of Pressure Vessels and Piping, 126–127, pp. 17–
28, 2014.
11 S. Zhang, and Y. Takahashi, “Evaluation of high temperature strength of a Ni-base alloy 740H for advanced
ultra-supercritical power plant” in Proceedings from the Seventh International Conference on Advances in
Materials Technology for Fossil Power Plants, 2013.
ANL-20/03 33
can be derived from creep test results (see 3-3.10). Given this information creep damage is defined
as 𝐷3 = 𝑡/𝑡Œ where 𝑡 is the time at fixed stress and temperature and 𝑡Œ is the time-to-rupture for
that stress and temperature.
Given a definition of fatigue and creep damage a creep-fatigue interaction diagram can be
constructed from a database of creep-fatigue test results. Creep fatigue tests, sometimes called
dwell-fatigue tests, are strain-controlled cyclic tests with hold periods at constant strain on either
or both of the maximum tension or maximum compression side of the cycle. Often these tests use
fully reversed loading cycles, thought to be the most damaging. The full stress/strain/time
hysteresis information is recorded during the test and this history used to calculate creep and
fatigue damage using the corresponding definitions. Plotting the results of a number of creep-
fatigue tests at different conditions produces a design envelope of the type sketched in Figure 3-
2.2.3.1 (see Section 3-3.12). Note the definition of creep and fatigue damage affects the shape of
the diagram – the same experimental data could produce different interaction diagrams using
different definitions of damage. Therefore it is critical to use an interaction diagram calculated
using consistent definitions of damage.
This process describe the generic creep-fatigue model used to calculate damage in the ASME Code
and in all three design options described here. Design Option A uses the ASME Section III,
Division 5, Subsection HB, Subpart B design by elastic analysis rules to determine the strains,
stresses, temperatures, and times required to calculate damage and check conditions against the
design envelope. While the elastic method is venerable, it has been well-validated by past
experience and has the important advantage over other methods in that it requires only a linear
elastic stress analysis of a single cycle repetition. This means that while the design rules are
somewhat complicated the design analysis is straightforward and easily accomplished.
For full details of the development of the design by elastic analysis method see the reference
material on the Section III, Division 5 rules, cited above. This following provides a brief overview.
This process must be repeated at each point in a component to find the limiting location.
34 ANL-20/03
1. Perform a linear elastic stress analysis of a loading cycle and extract the elastic strain range,
using the ASME definition of effective strain range. The Section III, Division 5 rules
provide methods for combining the effects of multiple types of load cycles. These rules
are not relevant here, as the current design criteria work with a single Design Cycle.
2. Modify the elastically calculated strain range to account for plasticity and creep using
bounding formula provided in the ASME Code.
3. Calculate the fatigue damage using the modified strain range.
4. Use one of the following two methods to construct a stress relaxation history:
a. Use the method of isochronous curves to calculate a stress relaxation history
starting from the modified strain range determined above. The method of
isochronous curves assumes that stress relaxation will follow the monotonic
stress/strain/time relation described by the isochronous curves. This is a
conservative bounding assumption – the real stress relaxation history will generally
be less severe.
b. Use a differential description of the material’s creep deformation to integrate a
stress relaxation history starting from the strain range determined above. This
approach is more accurate than the method of isochronous curves but requires
solving an ordinary differential equation.
5. Use the stress relaxation history, divided by a design factor, to calculate creep damage.
6. Consult the creep-fatigue interaction diagram to determine whether the point passes the
design criteria.
The modified rules for CSP structures given here follow this general process with several
modifications from the ASME rules.
First, the design margin is modified to account for the relative consequences of failure between
CSP and nuclear structures. The fatigue damage in the ASME nuclear code uses fatigue diagrams
with factors of 20 on the number of cycles to failure and 2 on the strain range. These rules, as
described in Section 3-3, adopt factors of 10 and 1.5. These rules divide the stress relaxation
history by a factor 0.9, retaining the ASME design by elastic analysis factor for creep damage.
Additionally, as with the ASME approach, these rules use a statistical minimum stress to rupture
to calculate creep damage, rather than the average property. The O’Donnell Porowski core stress
and a material dependent relaxation strength are used to bound the amount of stress relaxation
calculated with the method of isochronous curves. The ASME rules apply a factor of 1.25 on both
quantities, here the rules do not increase these lower relaxation limits by any factor, instead using
the average material property. Note the current rules provide values of the relaxation strength (see
3-3.9), whereas the current ASME rules require the designer to determine this from material test
data or a separate relaxation analysis.
Second, the ASME method is simplified by disallowing the use of several alternate methods for
calculating strain ranges and relaxation histories. This is done solely to simplify the design
process.
Finally, as noted in Figure 3-1.1 there are two types of cyclic stress relaxation for classical high
temperature structures modeled with elastic, perfectly-plastic, power-law creep material under
cycles consisting of instantaneous loading, hold at constant load, and instantaneous unloading:
ANL-20/03 35
resetting relaxation where each cycle starts out at the same, high value of stress and global
relaxation where the structure follows a global, constant-load relaxation curve with small elastic
transients during the load/unload periods. The global relaxation mode is highly desirable from a
design point of view as these sorts of structures spend more time at lower levels of stress and hence
accumulate less creep damage. The ASME Section III, Division 5 rules contain a test to determine
which mode will occur for a given component, based on the classical condition of reverse plasticity
causing load resets. If the structure falls into the global relaxation mode the ASME rules allow an
alternate, less damaging calculation of creep damage based on a global relaxation history.
The current rules do not allow this alternate method, essentially forcing the designer to treat CSP
components as if they follow a resetting relaxation history. As described in more detail below,
Alloy 740H and other high strength nickel alloys suitable for service at the high metal temperatures
envisioned in future CSP concepts are unlikely to plastically yield during service. This might
suggest that CSP structures will follow a global relaxation mode.
Figure 3-2.2.3.2. Stress relaxation histories for various cycles with the same total strain range but different loading
cycle types.
However, Figure 3-2.2.3.2 demonstrates a problem with the classical division of structures into
global and resetting relaxation. The past work establishing these categories was predominantly
aimed at nuclear applications where the components are loaded and unloaded quickly compared
to long hold periods at fixed conditions. CSP facilities do not operate in this mode, as they
generally load follow the incident solar radiation which increases and decreases continually
through the daily operation cycle. The figure shows three relaxation histories for an elastic-power
law creep material, all starting from the same initial strain: steady relaxation at fixed strain, cyclic
relaxation though a fixed strain range where the strains are applied instantaneously, held fixed,
and instantaneously reversed, and cyclic relaxation through the same range where the strains are
gradually applied and removed in a triangle wave pattern (see inset). Both cyclic histories undergo
the same total strain range and have the same total period, the only difference is in how the strains
and applied and removed.
The instantaneous and gradual cycles produce very different relaxation histories. The
instantaneous cycle validates the classical global relaxation mode where the cyclic relaxation
history follows the steady relaxation curve with short elastic transients. In contrast, the gradual
cycle produces very little global relaxation over the 10 cycles shown the diagram, more closely
resembling a resetting relaxation history.
36 ANL-20/03
Power law creep materials relax faster at high values of strain. This behavior explains the
difference between the instantaneous and global cycles. The instantaneous cycle always relaxes
at the highest value of strain and therefore relaxes comparatively quickly. The gradual cycle, on
average, relaxes at lower values of strain and hence much more slowly. This simple example
demonstrates that creep resistance is not always a good thing – materials that creep faster also shed
load faster. This example also justifies the restriction in the design rules disallowing the use of the
alternate, global relaxation calculation even if the structure would qualify for the alternate using
the ASME rules. CSP cycles are much more akin to the gradual load cycle in the example, as
opposed to nuclear structures which operate more like the instantaneous cycle.
The component is slowly globally relaxing in Figure 3-2.2.3.2 as reverse plasticity is not causing
a full cycle reset. Design Options A and B conservatively neglect this slow relaxation. Option C
accounts for it using inelastic analysis.
Part 4 provides some perspective on the relative importance of time-independent buckling, time-
dependent buckling, and the other design criteria. The fundamental conclusion is that time-
independent buckling under lateral wind loading may be a significant design concern whereas
time-dependent creep buckling is unlikely to constrain designs due to the low steady primary loads
envisioned in future CSP systems. However, the design criteria provide methods for both time
independent and creep buckling to provide guidance for systems that may have non-standard high,
steady primary loads, for example systems with substantial self-weight. Additionally, the creep
buckling rules come at little conceptual cost as they share a similar analysis method to the time
independent buckling method.
For time independent buckling the rules require a direct stability assessment of the structure using
the constitutive response given by the time-independent (hot tensile) isochronous curves. As
described in 3-3.14 these represent a best-estimate, average properties description of the
monotonic tension response of the material. It is therefore reasonable to use this response to
represent the deformation of a structure under time-independent loading conditions, using standard
J2 flow theory to extend the uniaxial tensile response to 3D stress states. The design criteria
requires a large-deformation inelastic analysis of the component and so directly accounts for elastic
and elastic-plastic buckling mechanisms. As described in Part 1 the Design Cycle must be
supplemented in this analysis with any transient, non-cyclic lateral loads – for example wind
loading.
The time-dependent buckling provisions use the method of isochronous curves. This approach has
been successfully used to evaluate buckling in high temperature structures12 and forms the basis
of the buckling exemption charts in ASME Section III, Division 5. The idea of the method is to
replace a time-dependent creep analysis with a time-independent elastic-plastic analysis using the
12 D. S. Griffin, “Design Limits for Elevated-Temperature Buckling” In Welding Research Council Bulletin 443
External Pressure: Effect of Initial Imperfections and Temperature Limits, pp. 11-26, 1999.
ANL-20/03 37
appropriate isochronous curve to represent the stress redistribution caused by creep over the design
life. This concept is empirical, but has been validated by past work13.
The analysis for the time-dependent buckling check is then fundamentally the same as for the time-
dependent check: a large-deformation, inelastic analysis of the component under steady-state
loading. The only differences are that for the time-dependent check the designer uses the
appropriate isochronous curve to represent the material’s constitutive response and the loading for
time-dependent buckling should only include the steady and periodic loads given in the Design
Cycle and not any additional, transient lateral loads. The rational here is the evaluation should
only include the loads causing significant creep deformation – brief transients like wind loading
will not cause creep deformation.
The load factors used in the buckling analysis are 1.5. This provides a 50% safety factor
accounting for uncertainties in the applied loads, the deformation properties of the actual material
(keeping in mind that the isochronous curves represent an average, not a lower-bound, response),
and any geometric imperfections in the fabricated component that might contribute to lower the
buckling load. The factor of 1.5 was extracted from the ASME Boiler & Pressure Vessel Code
Section III, Division 5 for Service Level D (most severe/least frequent) Loads. This means the
current rules adopt the most aggressive/least conservative ASME factor, reflecting the lower
consequences of failure in CSP versus nuclear systems.
High temperature nickel superalloys like 740H are designed to have high strength at elevated
temperatures. Newer alloys like Alloy 617, Alloy 740H, and Alloy 282 are so strong that they are
unlikely to deform plastically in service – note the yield strength for Alloy 740H exceeds 400 MPa
even above 800º C. This does not mean that these materials will remain linear elastic – creep
deformation can and does occur even below the time-independent yield point. However, assuming
the material does not undergo time-independent plastic deformation greatly simplifies the design
process. Options B and C take advantage of this assumption.
CSP components undergo daily cycling and will typically descend to temperatures below the creep
range during the inactive night periods. This means during each daily cycle they have only a
relatively short period of time to undergo stress relaxation. Additionally, high temperature nickel
alloys like 740H are creep-resistant, meaning they undergo relatively little creep deformation
under typical operating conditions.
Design Option B takes advantage of these material characteristics. Using Design Option B
requires that the component remain linear elastic. This criteria is checked by analyzing the Design
Cycle with a linear elastic stress analysis and comparing the resulting stress intensities to the
temperature dependent material yield strength. As the design values of 𝑆. typically lower-bound
the actual material properties such an analysis adequately demonstrates that the component will
not yield, even accounting for uncertainty in material properties.
13 ibid.
38 ANL-20/03
Design Option B uses the same primary load, ratcheting, and buckling methods as Method A.
Simplifications to the O’Donnell-Porowski approach for strong materials would be possible, but
the simplified process is not significantly easier to execute than the base approach required in
Method A.
The simplifications in Design Option B are in the creep-fatigue analysis. The method applies two
basic simplifications to the ASME Section III, Division 5, Subsection HB, Subpart B,
Nonmandatory Appendix T design by elastic analysis method:
1. The elastically-calculated effective strain range used to calculate fatigue damage only
needs to be modified for creep deformation, not plasticity.
2. The elastically-calculated stresses corresponding to the Design Cycle loading can be used
to calculate creep damage.
The first simplification is justified as the material does not undergo time independent plastic
deformation. Figure 3-1.1 describes the second simplification (compare the reset-relaxation to the
CSP-type cycle). It relies on two observations made above:
1. A resetting cyclic stress relaxation history reasonably bounds the response of CSP
components operating at high temperatures. High strength materials like Alloy 740H will
not yield in service and so the elastically-calculated stresses are a reasonable starting point
for a relaxation damage calculation.
2. High strength materials like Alloy 740H will not undergo significant stress relaxation
during the relatively short cycle period typical for CSP components. As such, the elastic
stress history without any additional relaxation caused by creep is reasonably
representative.
Like Option B this design option takes advantage of the assumption that high strength materials
like Alloy 740H will not yield in service. As with Option B this requires a preliminary linear-
elastic analysis to confirm the Design Cycle loads do not cause yielding.
Design by inelastic analysis has long been considered the gold standard for high temperature
design14 – likely to produce the most accurate, least-over conservative design analysis and the
ANL-20/03 39
most efficient final designs. However, there are two practical problems in implementing design
by inelastic analysis in engineering practice:
1. Accurate constitutive models for high temperature cyclic deformation are difficult to
develop. There is comparatively little high temperature deformation data available to
calibrate models and developing a mathematical description of high temperature cyclic
plasticity is notably challenging.
2. A complete inelastic analysis requires the analysis to simulate the full service history of a
component. For CSP components this might be 10,000 daily cycles, each requiring
multiple time steps to capture the details of the daily loading cycle. This results in analysis
with on the order of millions of time steps. These cannot be parallelized using current
methods meaning that the calculation time becomes infeasibly long even for very small
finite element models.
Design Option C attempts to retain some of the advantages of a design by inelastic method while
overcoming these two disadvantages. This option retains the primary load and buckling criteria
from Design Option A. It alters the approach for the ratcheting strain and creep-fatigue checks.
Design by inelastic analysis requires generating a full time/strain/stress history for each point of
the component. Once this history has been generated the accumulated strain over the component’s
service life can be directly calculated and compared to appropriate service limits. Similarly, the
strain/time history can be used to calculate a fatigue damage fraction, the stress/time history used
to calculate a creep damage fraction, and these used in conjunction with the creep-fatigue damage
envelope to assess the component against the creep-fatigue limits.
This design method retains the ASME effective strain range formula for converting a 3D strain
history to a scalar effective strain. The Section III, Division 5 rules use the Huddleston model15
for generating an effective stress from a 3D stress relaxation history for calculating creep damage.
This model is more accurate than simply using the von Mises effective stress, as in the current
design rules, but requires multiaxial creep rupture data to calibrate. This data is not available for
Alloy 740H and so the current rules recommend using the von Mises stress based on experience
with Alloy 617. This material is similar to 740H in that both are γ’ strengthened Ni-alloys and
multiaxial creep testing on Alloy 617 demonstrates that using the von Mises stress adequately
captures multiaxial rupture in that material. The design rules use the von Mises strain as a metric
for assessing strain limits. This choice is somewhat arbitrary, but reasonable given the use of the
von Mises stress in the creep damage calculation.
For ratcheting strain the design method here doubles the pointwise ASME strain limits from 5%
for base metal and 2.5% for welds to 10% for base metal and 5% for welds. Additionally, the
creep damage calculation uses a safety factor of 0.9, instead of the 0.7 recommended by ASME.
Both changes reflect the lower consequences of failure for CSP systems versus nuclear
components.
15 R. L. Huddleston, “An Improved Multiaxial Creep-Rupture Strength Criterion” Journal of Pressure Vessel
Technology, 107(November), pp. 421–429, 1985.
40 ANL-20/03
The key challenge is calculating a representative time/strain/stress history corresponding to many
repetitions of the design cycle. The design method takes advantage of the material strength by
neglecting plastic deformation. This greatly simplifies the task of constructing a material model,
as an elastic-creep is sufficient to describe the material deformation.
The final constitutive model described in 3-3.15 is deliberately further simplified to a temperature
and stress-dependent power law description of creep that does not use any internal variables to
track material state. This greatly abstracts the representation of creep deformation – the model
does not capture primary or tertiary creep – but also allows the designer to use a steady cyclic
analysis of the component to describe deformation and stress relaxation over an arbitrary large
number of load cycles. This means a relatively short cyclic analysis of a few repetitions of the
design cycle can be used to calculate deformation and damage over the full service life of the
component.
A history-independent material model of the type used here will, under repeated applications of a
cyclic load history, eventually reach a cyclic steady state where the stresses and strain rates become
periodic16. Design Option C takes advantage of this concept by requesting the designer first
determine this steady, periodic response and then calculate deformation and damage over the
whole cyclic history of the component as if this steady response was repeated for the entire service
life.
The clear advantage of this approach is that the designer does not need to simulate the full cyclic
history of the component. The steady response can generally be identified by repeating the Design
Cycle for tens of applications, not thousands.
There are however several disadvantages. First, as noted above, the history-independent
description of creep cannot capture real, observed material behavior like primary and tertiary
creep. Tertiary creep can safely be neglected as the creep damage calculation will screen structures
that enter this deformation regime. Neglecting primary creep is conservative for calculating stress
relaxation damage as faster creep rates are actually helpful in relaxing stress, but not necessarily
for calculating ratcheting strain accumulation. However, as creep-fatigue is likely to be the
controlling design criteria for future CSP structures this approximation can be accepted.
Second, applying the steady-cyclic response to the full service history neglects the cyclic-transient
part of the deformation history. Depending on the material response this might mean
approximating the damage and strain accumulated over the first few hundred cycles. However, as
CSP systems are typically designed for thousands of cycles this approximation is reasonable.
Additionally, this design method does not completely neglect damage and strain accumulated
during this period, but rather approximates the response of the component in this regime with the
steady-cyclic response. This approximation is generally reasonable, even if it is inexact.
16 C. O. Frederick, and P.J. Armstrong, “Convergent internal stresses and steady cyclic states of stress” Journal of
Strain Analysis for Engineering Design, 1(2), pp. 154–159, 1966.
ANL-20/03 41
3-3. Commentary on the 740H Material Data
Currently, data in Section III, Division 5 of the ASME Code is provided in either graphs or tables.
A more modern design code would provide equations so that the various material property
correlations could be easily digitized. This set of data endeavors to provide equations, where
possible. However, some of the data is still provided in tables to be compatible with current
Section II properties (for example 𝑆. ).
The intent of this set of design rules is to allow the use of Alloy 740H in any non-cast product
form. Currently welded tube and pipe is not included in ASME Code Case 2702, adding Alloy
740H for Section I use. If the Code Case is amended in the future to include welded tube these
design data could be relatively easily updated with conforming allowable stress values. The weld
reduction factors in 2-8 would provide a way to design seam welded material against the creep-
fatigue and ratcheting limits.
Code Case 2702 includes additional requirements on the material heat treatment, providing two
options (a and b) to satisfy the additional requirements. These design rules retain those additional
restrictions.
This information comes from the Alloy 740H material data sheet.
This information comes from the Alloy 740H material data sheet. The instantaneous coefficients
of thermal expansion were calculated from the provided mean coefficients.
This data comes from ASME Code Case 2702. The values above 800º C are extrapolated.
This data comes from ASME Code Case 2702. The values above 800º C are extrapolated.
ASME Code Case 2702 provides the Section I allowable stresses for Alloy 740H up to 800º C.
This allowable stress is defined by the lesser of (at each temperature);
42 ANL-20/03
2. The minimum specified yield strength of the material multiplied by a factor of 0.67 or 0.9
depending on the material ductility (this criteria does not control for Alloy 740H).
3. The minimum specified tensile strength of the material multiplied by the ratio of the
average temperature-dependent tensile strength to the room temperature tensile strength
multiplied by a factor of 1.1/3.5.
4. The minimum specified yield strength of the material multiplied by the ratio of the average
temperature-dependent yield strength to the room temperature yield strength multiplied by
a factor of 0.67 or 0.9, depending on the material.
5. The average creep rupture strength of the material at 100,000 hours multiplied by a factor
less than 0.67.
6. The minimum creep rupture strength of the material at 100,000 hours multiplied by factor
of 0.8.
7. The full stress to cause a strain of 1% in 100,000 hours.
Full details of this calculation are provided in ASME Section II, Mandatory Appendix 1.
Section III, Division 5 uses a similar allowable stress 𝑆- adopted for use here. The definition of
this allowable stress is identical, except where the value of the time-dependent allowable stress
𝑆6D at 300,000 hours exceeds the value defined above. This extra provision is ignored here, as
these design rules do not use 𝑆6D .
These rules adopt the values provided in ASME Code Case 2702 up to 800º C. The remaining
value at 850 º C is in the regime controlled by creep properties and values were calculated using
the rupture and creep deformation models described below.
This is a time-independent allowable stress used in Section III, Division 5. This quantity is not
used as an allowable stress in the design criteria presented here. However, it is used in the creep-
fatigue and ratcheting rules to determine a typical initial cyclic relaxation stress using the 3𝑆6̅
criteria. As such the designer needs these values to complete a design analysis.
As described in Part 2, the design by elastic analysis creep-fatigue and ratcheting criteria construct
a relaxation analysis based on a starting strain value. Notionally, for a power-law creep material
stress relaxation would continue until the stress in the component reaches the primary stress.
However, in practice material effects, quantified by backstresses in cyclic plasticity models, tend
to provide a material relaxation strength – a stress below which relaxation ceases – that can be
ANL-20/03 43
higher than the component primary stress. The design by elastic analysis process then includes
rules preventing the stress from falling below the relaxation strength.
The material data in Section III, Division 5 and related ASME Code documents are nearly
complete – a designer can complete a design without recourse to additional, outside information.
However, the relaxation strength is an exception – it is not tabulated in the current Code. However,
HBB-T-1324 does allow the designer to use a uniaxial relaxation analysis starting from a stress of
1.5𝑆6 to determine the relaxation strength. Relaxation models for the HBB Class A materials are
not currently provided. For these design rules the tabulated relaxation strengths are based on the
inelastic material model described in 3-3.15.
The minimum stress to rupture describes the stress that will cause rupture in a given time at a given
temperature with reasonable lower-bound material response. The direct experimental
measurement of this value is a creep rupture test. In these tests uniaxial specimens are loaded in
temperature-controlled furnaces with a dead-load stress and left to deform until they rupture. The
minimal experimental measurements are the controlled dead-load stress and temperature and the
corresponding rupture time. More sophisticated tests might measure the deformation as a function
of time, typically using a combination of extensometers and linear variable differential
transformers (LVDTs) to measure deformation directly from the specimen.
However, it is essentially impossible to run sufficient creep tests to fully interpolate between the
two experimental conditions (stress and temperature) for all values needed in a design calculation.
Moreover, for some applications very long-term rupture stresses are required (30+ years), whereas
the test time for the available creep test database is often much shorter, on the order of years at the
most. Therefore, a predictive model for the rupture stress as a function of rupture time and
temperature is first calibrated to the experimental data and then this model is used to generate
design rupture stresses.
For this purpose 54 creep rupture tests on Alloy 740H were identified and collated. Table 3-
3.10.1.1 shows the raw experimental data used to develop the rupture correlation below. The vast
majority of these tests were extracted from the DOE report summarizing the fossil energy
development program for Alloy 740H17, however a few confirmatory tests were completed at INL
(marked with an asterisk in the table).
17 R. Purgert et al. Boiler materials for ultra supercritical coal power plants. DOE technical report,
2015.
44 ANL-20/03
Temperature (ºC) Rupture Stress (MPa) Rupture time (hours)
700 340 5227.4
700 480 512.5
700 300 14871.3
700 480 517.4
700 320 9745.1
700 265 30956.7
700 420 1082
700 340 5688
750 370 142.8
750 320 358.3
750 265 1275.4
750 180 8034
750 180 9787.9
750 370 275.4
750 300 984.6
750 220 7201.6
750 180 22896
750 370 311.9
750 320 658.5
750 300 1020.2
750 265 2185.4
750 220 7355.2
750 180 20789.4
750 370 296.7
750 320 484.8
750 300 914.6
750 220 7382.8
750 180 24061
750 300 723
750 265 2345
750 370 229.7
750 370 242
750 320 553
750 180 18756
800 250 173
800 250 279
800 180 1491
800 130 6883
800 180 1634
800 120 9855
800 110 15864
800 250 259
800 180 1497
800 150 3609
800 250 199
ANL-20/03 45
Temperature (ºC) Rupture Stress (MPa) Rupture time (hours)
800 100 17674
750* 265 3168
750* 320 861
800* 200 1535
850* 100 3053
Table 3-3.10.1.1. Creep rupture database collated to define minimum stress to rupture design data. Tests marked
with an asterisk were conducted at INL as part of the project sponsoring the development of these design rules.
A Larson-Miller18 correlation is used to interpolate and extrapolate the rupture test data to
determine design values of minimum stress-to-rupture. The Larson-Miller parameter
𝐿𝑀𝑃 = 𝑇(log 𝑡9 + 𝐶 ),
with 𝑇 being the absolute temperature, 𝑡9 the rupture time, and 𝐶 a parameter, generally correlates
to the logarithm of the stress to rupture. This correlation can be used to fill in values of rupture
stress for various conditions and extrapolate the short-term data to longer rupture times.
Larson-Miller correlations can be generated by finding the optimal value of the parameter 𝐶 and
the corresponding polynomial relation between the Larson-Miller parameter and the log of the
rupture stress that best fit the experimental data. Figure 3-3.10.2.1 illustrates the process used to
develop the log-linear Larson-Miller relation used here to establish values of the minimum rupture
stress. Current ASME practice is to develop a 95% confidence lower bound prediction interval on
the data and use that lower bound to determine the minimum stress to rupture. That process was
followed here. The figure shows both the mean-property Larson-Miller model (used as the average
stress to rupture in generating 𝑆- and the creep-fatigue interaction diagram) and the 95%
confidence model, used here to generate the minimum stress to rupture table. The figure also
shows the best-fit value of the parameter 𝐶 and the Pearson correlation coefficient of the log-linear
relation.
18 F. R. Larson, J. Miller, “A time-temperature relationship for rupture and creep stresses” Transactions of the ASME,
74, pp. 765–771, 1952.
46 ANL-20/03
Figure 3-3.10.2.1. Larson-Miller correlation used to generate the average and minimum stress to rupture for Alloy
740H.
A high temperature fatigue curve is a relation between a fixed cyclic strain range and temperature
and the resulting number of cycles to failure. Typically these diagrams are generated using fully-
reversed, strain-controlled fatigue tests at fixed temperature. These tests uniaxially cycle a sample
over a fully-reversed strain range, holding temperature constant using a temperature-controlled
furnace or induction heater. The test cycle is repeated until the sample fails. To protect equipment
and specimens, failure may be defined by a set percentage drop in the maximum stress due to the
presence of a crack. Design fatigue curves interpolate these test results, typically using strain
range/cycles to failure relations for fixed temperatures.
Table 3-3.11.1.1 lists the fatigue test data used to generate the design correlations. All these tests
are for fully-reversed strain cycling at a fixed temperature and strain rate. The table lists the
temperature, strain range, number of cycles to failure, and the source of the data. Test results were
collected from the literature and dedicated experiments run at INL as part of this project.
ANL-20/03 47
Temperature (ºC) Strain range (%) Cycles to failure Source
700 1.49 293 Zhang and Takahashi19
700 1.19 467 Zhang and Takahashi
700 0.97 1931 Zhang and Takahashi
700 0.69 12046 Zhang and Takahashi
700 0.66 15104 Zhang and Takahashi
700 0.53 287843 Zhang and Takahashi
750 0.60 14987 INL
750 0.6 19023 INL
750 0.60 7950 INL
750 1.00 1767 INL
750 1.00 1095 INL
760 1.40 147 Jena20
760 1.20 383 Jena
760 1.01 671 Jena
760 0.80 1457 Jena
760 0.60 6566 INL
800 1.00 653 INL
850 0.40 19969 INL
850 0.40 47048 INL
850 1.00 409 INL
850 1.00 402 INL
Table 3-3.11.1.1. Fatigue test data collated to generate the design fatigue relations.
Typically, strain-based fatigue curves show decreasing cycles to failure as the temperature
increases at a fixed strain range. ASME procedure is to provide a set of temperature dependent
curves. The designer uses the curve corresponding to the highest metal temperature at the point
of interest.
The data from 3-3.11.1.1 was grouped into two sets: 700º C and below and 700º C to 850ºC.
Because this is, to our knowledge, the total fatigue dataset for Alloy 740H the design rules use the
curve based on the 700º C for all temperatures less than 700º C. This is almost certainly
conservative. An optimized set of design data would include fatigue tests at lower temperatures,
allowing a less restrictive design fatigue correlation.
Once grouped, two inverse power law curves were fit to the data using least squares regression.
These two curves, corresponding to the two temperature groups, provide the mean property fatigue
correlation.
19 S. Zhang, and Y. Takahashi, “Creep and Creep-Fatigue Deformation and Life Assessment of Ni-Based Alloy 740H
and Alloy 617” In ASME 2018 Pressure Vessels and Piping Conference, pp. V06AT06A060-V06AT06A060.
American Society of Mechanical Engineers, 2018.
20 P. S. Jena et al. “Low cycle fatigue behavior of nickel base superalloy IN 740H at 760° C: Influence of fireside
corrosion atmosphere” International Journal of Fatigue, 116, pp. 623-633, 2018.
48 ANL-20/03
Design fatigue charts include factors of safety. The ASME nuclear code uses a factor of 2 on the
strain range and 20 on the cycles to failure. Based on the lower consequences of failure for CSP
systems, the design fatigue curves here use a factor of 1.5 on strain range and 10 on cycles to
failure. Figure 3-3.11.2.1 shows that these factors bound the uncertainty in the fatigue data and
provide a healthy design margin guarding against detrimental environmental/fatigue interactions,
initial defects, and other uncertainty.
Figure 3-3.11.2.1. Chart comparing the experimental fatigue data, the mean property fatigue correlations, and the
design fatigue correlations.
The use of a creep-fatigue interaction diagram is explained in detail above in 3-1 and 3-2.2.3.
These plots are an allowable damage envelope plotting fatigue damage on one axis and creep
damage on the other axis. Typically, they are generated from a set of creep-fatigue tests. Creep-
fatigue tests are strain-controlled cyclic experiments at fixed temperature and fixed fully-reversed
strain range. They are different from fatigue tests in that a hold period at fixed strain is included
in the cycle either at the maximum tensile, maximum compression, or both ends of the cycle.
ANL-20/03 49
Table 3-3.12.1.1 lists the test used to generate the Alloy 740H creep-fatigue interaction diagram.
All the tests but one were run at INL as part of this project. The single test sourced from the
literature21 is marked with an asterisk. Early in the experimental program two tests with the same
conditions but holds on the opposite ends of the cycle were used to establish that tensile holds are
more damaging for Alloy 740H. This is in agreement with similar results on other high
temperature nickel alloys. As such, tensile holds were used for the remainder of the tests.
The results of a creep-fatigue test are a stress/strain/time history and a number of cycles to failure.
These results can be converted to a fatigue damage fraction by calculating the number of cycles to
failure for a pure fatigue test at the experimental strain range and temperature using the mean
property fatigue curve and dividing the actual number of cycles to failure by this number.
Similarly, a creep damage fraction can be calculated using the time fraction approach via the
integral equation
𝑑𝑡
𝐷3 = ²
𝑡Œ
where 𝑡Œ is the time-to-rupture at a given value of uniaxial stress calculated using the mean
property time to rupture correlation developed above. Notionally this integral would span the full
time test data. In practice the relaxation curve at 𝑁/2, i.e. half the experimentally-observed
number of cycles to failure, is used to represent the whole test.
21 Zhang, Shengde, and Yukio Takahashi. "Creep and Creep-Fatigue Deformation and Life Assessment of Ni-Based
Alloy 740H and Alloy 617." In ASME 2018 Pressure Vessels and Piping Conference, pp. V06AT06A060-
V06AT06A060. American Society of Mechanical Engineers, 2018.
50 ANL-20/03
Figure 3-3.12.2.1 shows the experimental data, converted to creep and fatigue damage fractions
and plotted on an interaction diagram. The figure also shows the interaction point selected for the
design rules: (0.05,0.05). The processes of converting the results of a creep-fatigue test to a point
on the diagram is automatic and requires no judgement beyond selecting the fatigue and creep
damage models. However, going from the collection of experimental data points to the interaction
envelope is a somewhat subjective process. The creep and fatigue damage calculation procedures
contain the design margin in the creep-fatigue design process and so creep-fatigue diagrams are
selected based on average material properties. As such, the final interaction point does not need
to bound the data – it can allow for some outliers falling underneath the curve. However, the final
diagram should not pass through the mean of the data as the points on the diagram represent many
different experimental conditions. That is, the scatter in the plot is not entirely due to variation in
material properties.
Figure 3-3.12.2.1. The creep-fatigue experiments plotted on an interaction diagram along with the final
recommended damage envelope.
Near weldments the creep-fatigue criteria apply a reduction factor on the minimum stress to rupture
in the creep damage calculation. These reduction factors are typically based on cross-weldment
rupture tests and may be temperature and time dependent.
The factor of 0.7 is drawn from ASME Code Case 2702, where it is specified for use in seam
welded construction. The weld rupture factor is based on very limited data – additional weld creep
data would be required to determine a more accurate value.
Note that for seam-welded tube or pipe reduced allowable stress values 𝑆- would be required, in
addition to the reduction factor applied to the rupture strength in the creep-fatigue damage
calculation procedure. As noted above in 3-3.7 the elevated temperature values of these reduced
allowable stresses could be calculated in Code Case 2702 is modified in the future to allow for
seam-welded tube or pipe product forms.
ANL-20/03 51
3-3.14. Isochronous stress-strain relations
ASME Section III, Division 5 design rules use hot tensile and isochronous stress strain curves to
represent creep deformation, both for evaluation of the ratcheting strain limits and in assessing
stress relaxation caused by creep. The isochronous curves can be read as the average stress to
accumulate some amount of total strain over some period of time. The procedure used to create
current Code curves for the Section III, Division 5 materials is to first fit a material model to the
available tensile and creep data and then use that model to generate the hot tensile and isochronous
curves. The difficulty in this analysis stems from the need for strain data during the creep test
rather than just tabulated rupture times. We adopted the same procedure to develop isochronous
stress strain curves for alloy 740H. The tensile test data used to determine the isochronous stress-
strain relations were collected at INL while most of the creep test data were gathered from
literature, except for four creep tests conducted at INL. The tensile test data consist of four elevated
temperature tensile tests spanning from 700°C to 850°C in 50°C increments. Table 3-3.14.1 lists
all the creep tests used. The table lists the temperature, applied constant stress, and the source of
the creep test data.
52 ANL-20/03
𝜀 = 𝜀% + 𝜀# + 𝜀3
The hot tensile curves are the outcome of this model when 𝜀3 = 0, i.e. when 𝑡 = 0, whereas the
isochronous curves are the output of the model for some fixed, non-zero time. The elastic strain is
calculated using the temperature dependent values of Young’s modulus, E for Alloy 740H.
𝜎
𝜀% =
𝐸
The plastic response of Alloy 740H was divided into two regions based on temperature. At
temperatures below and equal to 800° C the composite model uses a Ramberg-Osgood model for
the plastic strain to capture the experimentally-observed smooth transition from elastic to work
hardening plastic behavior. Above this temperature the model uses a Voce hardening model to
capture a quick transition to a nearly perfectly-plastic response. The composite model for the
plastic strain is then
0 𝜎 ≤ 𝜎l
⎧ 𝜎 − 𝜎 )
⎪ › l 600℃ ≤ 𝑇 ≤ 800℃
⎪ 𝐾œ • 𝜎 > 𝜎l
𝜎l
𝜀# = 0 𝜎 ≤ 𝜎B
⎨
⎪›−1 ln Ÿ1 − 𝜎 − 𝜎B
⎪
𝜎 > 𝜎B
𝑇 > 800℃
⎩ 𝛿 𝜎# − 𝜎B
All the model parameters were calibrated based on single hot tensile test data at four different
temperatures – 700°C, 750°C, 800°C, and 850°C. We first fit the experimental data with the
plasticity model (Ramberg-Osgood or Voce hardening) considering experimental elastic modulus
which is then modified to match with the Young’s modulus and the average yield strength. The
values of the average yield strength for Alloy 740 were collected from the background document
of ASME Code Case 2702. These modifications are done only in the elastic regime of the hot
tensile curve and do not affect the shape of the curve in the plastic regime. Figure 3-3.14.2.1 plots
the experimental (red curve) and model (black curve) hot tensile curves. Due to the unavailability
of hot tensile data below 700°C, we considered the same plasticity model constants for all the
temperatures between 600°C and 700°C.
To model the time-dependent strain, 𝜀3 we adopted a simple creep model for alloy 740H.
where 𝜀̇3 is some constant, average creep rate, which is a function of temperature and stress. ASME
Section III, Division 5 design isochronous curves only provide data out to 2.2% total strain which
at most represents about 2% creep strain. The model also assumes that this average rate over the
first 2% of creep strain is approximately equal to the average rate over the first 1% of creep strain.
This allows the model to use the time-to-1% data for calibration. To determine the average creep
rate over the first 1% of creep strain two different methods were used – (1) divide the 1% creep
strain by the time and (2) convert the creep curve data to a plot of creep strain rate versus creep
ANL-20/03 53
strain and then average over the first 1% to produce a similar mean rate. These two rates were then
averaged again to determine the average creep rate over the first 1% of creep strain to the applied
stress and temperature. A model must be developed for this average creep rate in order to
interpolate the data to all the conditions required to generate the design isochronous curves.
Figure 3-3.14.2.1. Experimental and hot tensile model curves for Alloy 740H at different temperatures. Hot tensile
model curves were determined by fitting the experimental data with the plasticity model (Ramberg-Osgood or Voce
hardening) and then modified to match with the Young’s modulus and average yield strength. The average yield
strength data were collected from the background document of ASME Code Case 2702.
54 ANL-20/03
We adopt a form developed by Kocks22 and Mecking23 for creep model of alloy 740H. Their
´
model posits a linear relation between the log-normalized material flow stress log £ and the
§i µ̇
normalized activation energy £¤¥ ln µ̇¶ . If this log-linear relation exists, the Kocks-Mecking model
can be converted into a model for the deformation strain rate as a function of the linear fit slope 𝐴
and intercept 𝐵.
¥ /(¦§i)
<£¤¥ /(¦§i)
𝜎 ‰£¤
𝜀̇ = 𝜀̇l 𝑒 œ •
𝜇
¹
Here 𝜇 is the material shear stress given as 𝜇 = C(Bº»), 𝑘 is the Boltzmann constant, 𝑇 is the
absolute temperature, 𝑏 is a characteristic Burgers vector, and 𝜀̇l is some reference strain rate. This
method was successfully implemented by Messner and Sham24 for modeling the creep
deformation in developing the isochronous stress-strain relationships for a similar nickel based
Alloy 617. Figure 3-3.14.2.2 plots the available Alloy 740H creep data using the average rate to
1% creep strain as the deformation strain rate and the applied values of stress and temperature. As
the Figure 3-3.14.2.2 shows, the Alloy 740 creep data nearly obeys the Kocks-Mecking form.
Based on this diagram, the model for the creep strain adopted for Alloy 740H is
<£¤¥ 𝜎 ‰£¤¥
𝜀3 = 𝜀-̇ 𝑒 ¦§i ( ) ¦§i 𝑡
𝜇
22 U. F. Kocks, “Realistic constitutive relations for metal plasticity” Materials Science and Engineering: A, 317, pp.
181-187, 2001.
23 Y. Estrin, and H. Mecking, “A unified phenomenological description of work hardening and creep based on one-
parameter models” Acta Metallurgica, 32, pp 57-70, 1984.
24 M. C. Messner, and T-L. Sham, "Isochronous Stress-Strain Curves for Alloy 617" In the Proceedings of the ASME
2019 Pressure Vessels and Piping Conference, PVP2019-93642, 2019.
ANL-20/03 55
3-3.15. Inelastic constitutive model
As noted above in the commentary on the design methods, assembling a complete inelastic
constitutive model capable of accurately capturing details of elevated temperature cyclic
plasticity/creep interactions is a very complicated task that has not been completed for the vast
majority of high temperature materials25,26,27.
However, Design Method C detailed above purposely uses a simplified description of the material
response in order to take advantage of the steady cyclic response of the system. In the general
case, the material description could be elastic, perfectly-plastic, with power law creep. However,
as described in 3-2.4, for Alloy 740H plasticity can be neglected, leaving an elastic-creep
constitutive response.
The specific inelastic model was then formulated by taking the elastic and creep strain equations
from 3-3.14 and supplementing them with the standard equation for the thermal strains, using the
coefficients of thermal expansion supplied in Part 2.
This constitutive model is only suitable for use with the simplified inelastic design method detailed
above. It is not suitable for a full inelastic analysis capturing all the details of elevated temperature
plasticity and creep-plasticity interaction. Such a model, for example using the Chaboche form28,
would be much more difficult to develop and require an extensive experimental database of creep,
creep-fatigue, and other cyclic tests including full stress/strain/time hysteresis information.
This temperature limit is the temperature at which the material accumulates 0.1% strain at the
allowable stress 𝑆- over 100,000 hours. It can be determined using the information in 3-3.7 and
3-3.14 above. This particular cutoff was invented for these rules. The rational is that the allowable
stress is a reasonable typical long-term stress level in components design using these rules. The
0.1% over 100,000 hour criteria is a typical ASME negligible creep threshold.
25 M. Yaguchi, and Y. Takahashi, “Unified Inelastic Constitutive Model for Modified 9Cr-1Mo Steel Incorporating
Dynamic Strain Aging Effect” JSME International Journal Series A, 42(1), pp. 1–10, 1999.
26 M. C. Messner, V.-T. Phan, T.-L. Sham, “A Unified Inelastic Constitutive Model for the Average Engineering
Response of Grade 91 Steel” In the Proceedings of the ASME 2018 Pressure Vessels and Piping Conference
PVP2018-84104, 2018.
27 Phan, V.-T., Messner, M. C. and Sham, T.-L. “ A Unified Engineering Inelastic Model for 316H Stainless Steel.”
In the Proceedings of the ASME 2019 Pressure Vessels and Piping Conference, PVP2019-93641, 2019.
28 J. L. Chaboche. “Constitutive equations for cyclic plasticty and cyclic viscoplasticity” International Journal of
Plasticity, 5, pp. 247–302, 1989.
56 ANL-20/03
negligible creep regime – a necessary requirement for applying the O’Donnell-Porowski theorem.
The specific temperature limit is the temperature at which 𝑆6 = 𝑆D at 100,000 hours. Section 3-
3.8 describes the Section III, Division 5 allowable stress 𝑆6 . 𝑆D is a time-dependent allowable
stress used in Section III, Division 5. The details of how it is calculated from creep-data are not
important, as the design document provides the temperature value for Alloy 740H.
ANL-20/03 57
Part 4: Sample problems
This chapter walks through two sample design problems, applying all three design methods to each
problem. The structural material used in all calculations is Alloy 740H, with design properties
from Part 2 above.
Figure 4-1.1 illustrates the first problem considered for evaluating design methods developed for
CSP systems. This problem is an axisymmetric representation of a tube in a cavity receiver. The
tube is 500 mm long and 2 mm thick. The outer diameter is 40 mm. For simplicity, we assumed
uniform heat flux on the outer surface of the tube and that the heat conduction analysis is done in
the steady state, even for the design criteria which require a transient analysis. This results in linear
temperature gradient along the length and thickness of the tube and therefore this problem can be
treated as an axisymmetric problem. This linear temperature distribution can be fully described
by providing the inner and outer tube metal temperatures as a function of time and axial position.
Only two points are required to define the axial gradient. Figure 4-1.1 shows the temperature and
pressure loading considered for this problem. The loading cycle includes warming up of the system
in the morning, steady state operation, five cloud events each with 8 minutes hold, cooling down
in the evening, and no operation during night. The design life of the tube is 30 years.
The purpose of this problem is to provide designers a simple reference problem to check their
understanding of the design methods.
Figure 4-1.1. (Sample problem 1) An axisymmetric representation of a single tube in a cavity receiver. Loading
profile shows the inner and outer wall temperature at the bottom (z = 0 mm) and top (z = 500 mm) ends of the tube,
respectively, and pressure exerted on the inner wall by the salt flowing inside the tube.
58 ANL-20/03
4-1.2. Design calculations based on Method 1
Figure 4-1.2.1 shows the Design Cycle. The daily load cycle can be divided into two service load
types – start-up/shut-down cycle and cloud events. Table provides the details of each event
considered in the composite cycle.
Figure 4-1.2.1 (Sample problem 1) Different load points during the loading cycle considered.
Service load types Associated load points and cycle period Frequency per design cycle
start-up/shut-down cycle A and B ; 12 hours 1
cloud events B and C ; 0.133 hours 5
Table 4-1.2.1 (Sample problem 1) Service load cycles and associated load points (illustrated in Figure 4-1.2.1) in the
daily load cycle and corresponding hold times.
4-1.2.2. Step-2: Transient elastic thermo-mechanical analysis for each service load
case and stress classification
We used MOOSE (Multiphysics Object Oriented Simulation Environment), an open source finite
element solver to perform the elastic thermo-mechanical analyses. We classify stresses due to
pressure as primary load and thermal stresses caused by the temperature gradient as secondary
load. There is no peak load.
ANL-20/03 59
4-1.2.3. Step-3: Primary load design check
Maximum primary load occurs at load point-B. Figure 4-1.2.3 shows the temperature distribution
in the tube and stress components along the thickness of the tube at maximum wall averaged
temperature location. Table 4-1.2.3 reports details of the primary load checks. First, all the stress
components were linearized to divide into membrane and bending components along the stress
classification line. The membrane and bending stress tensors were then used to determine the
stress intensities in Table 4-1.2.3 As indicated in the table the design passes both the criteria in
primary load checks.
Figure 4-1.2.3. (Sample problem 1) Temperature distribution in the tube and through thickness elastic stress
components at maximum wall averaged temperature location under primary load at load point -B.
Design Method 1 uses the O’Donnell-Porowski approach, described in Section III, Division 5,
HBB-1332 Test No. B-1 for ratcheting checks. In this approach, an effective creep stress
parameter, Z is determined from a primary stress parameter, X and a secondary stress parameter,
Y as shown in Figure 4-1.2.4.1 The effective creep stress parameter is used to calculate the
effective creep stress which is then used to determine the ratcheting creep strain using isochronous
stress-strain curves. The definition of X and Y are
60 ANL-20/03
𝑃¤
𝑋 = (𝑃: + ) /𝑆
𝐾D 6EP .:
𝑌 = (𝑄Œ )6EP /𝑆.:
where,
Á
(𝑃: + ÃÂ)6EP = the maximum value of the primary stress intensity, adjusted for bending via 𝐾D ,
Ä
during the cycle being evaluated.
(𝑄Œ )6EP = the maximum range of the secondary stress intensity during the cycle being considered
𝑆.: = is the 𝑆. value corresponds to the lower of the wall averaged temperature for the stress
extremes defining secondary stress range, 𝑄Œ .
𝐾D = (𝐾 + 1)/2
𝐾 is 1.5 for across-the-wall bending of shell structures or rectangular sections, see HBB-3223 (c)
(6) in Section III Division 5.
Figure 4-1.2.4.1. Illustration of determining the effective creep stress parameters from Section III, Division 5, Figure
HBB-T-1332-1.
ANL-20/03 61
It should be noted that, the average wall temperature at one of the stress extremes defining the
secondary stress intensity range must be below the temperature listed in Section III, Division 5,
HBB-T-1323, given as 600℃ for 740H in the description of the design method above.
The creep ratcheting strain increment for a load cycle is evaluated by entering the isochronous
stress strain curves at the maximum wall temperature and effective core, 𝜎3 stress during the load
cycle with the stress held constant for the entire service life. An example of creep ratcheting strain
determination is shown in Figure 4-1.2.4.2.
Since the start-up/shut-down service load includes the extreme temperature profile and the total
time of the day, considering only the start-up/shut-down load should provide conservative
estimation for ratcheting design. Table 4-1.2.4 provides all the calculation details of the ratcheting
design check. Figure 4-1.2.4.3 shows the stress components under secondary loading at load point
B. Stress components are shown at two different locations – maximum wall averaged temperature
and maximum von Mises stress. As indicated in the Table 4-1.2.4, the maximum ratcheting strain
in the structure is less than 2% and therefore the design passes the ratcheting check. Note that, a
design must pass the ratcheting design check before it is checked for creep-fatigue damage.
Figure 4-1.2.4.2. Illustration of determining the creep ratcheting strain increment from isochronous stress strain
curves.
62 ANL-20/03
Figure 4-1.2.4.3. (Sample problem 1) Temperature and von Mises stress distribution in the tube and through
thickness stress components at maximum wall averaged temperature and maximum von Mises stress locations under
secondary load at load point B.
ANL-20/03 63
Stress classification line-1 Stress classification line-2
shown in Figure 4-1.2.4.3 shown in Figure 4-1.2.4.3
6EP
𝑇ÆEFF EL%9E(%$ 755.2 °C 752.6 °C
6') 30 °C 30 °C
𝑇ÆEFF EL%9E(%$
𝑇 6EP 770 °C 767.3 °C
6') 621.0 MPa 621.0 MPa
𝑆.: (at 𝑇ÆEFF EL%9E(%$ )
𝐾 1.5 1.5
𝐾D = (𝐾 + 1)/2 1.25 1.25
Á
(𝑃: + ÃÂ)6EP 20.53 20.53
Ä
(𝑄Œ )6EP 85.06 96.81
Á 0.033 0.033
𝑋 = (𝑃: + ÃÂ)6EP /𝑆.:
Ä
𝑌 = (𝑄Œ )6EP /𝑆.: 0.137 0.156
𝑍 using Section III, Division 0.033 0.033
5, Figure HBB-T-1332-1
´
𝜎3 from 𝑍 = È Ç 20.53 20.53
ÉÊ
Service life 30 years = 131400 hours 30 years = 131400 hours
Ratcheting strain at the end 1.723e-5 % 1.393e-5 %
of service life
Ratcheting design criteria: PASS! PASS!
2% for base metal
Table 4-1.2.4 (Sample problem 1) Ratcheting design check according to Method 1.
According to Section III, Division 5, a design is acceptable if the creep and fatigue damage satisfy
the following relation:
) ∆D
∑Ì ( )Ì + ∑§( )§ ≤ 𝐷
ej i j
where D is the total creep-fatigue damage and the first and second terms on the left side are fatigue
damage, 𝐷J and creep damage, 𝐷3 , respectively. In the fatigue damage term, (𝑛)Ì is the number
of repetitions of cycle type j and (𝑁$ )Ì is the number of design allowable cycles for respective
cycle type; while in the creep damage term, (𝑇$ )§ is the allowable time duration for a given stress
at the maximum temperature occurring in the time interval k and (∆𝑡)§ is the duration of the time
interval k.
The design allowable cycles for fatigue damage is determined by entering fatigue curves at total
strain range, 𝜖D . Total strain range, 𝜖D is calculated using equation HBB-T-1432-16:
𝜖D = 𝐾L ∆𝜖6-$ + 𝐾∆𝜖3
where 𝐾 is the local geometric concentration or equivalent stress concentration factor determined
by dividing effective primary plus secondary plus peak stress divided by the effective primary plus
secondary stress, 𝐾L is the multiaxial plasticity and Poisson ratio adjustment factor, ∆𝜖3 is the
creep strain increment, and ∆𝜖6-$ is the modified maximum equivalent strain range.
64 ANL-20/03
∆𝜖6-$ is calculated using equation Section III, Division 5, HBB-T-1432-12:
È∗
∆𝜖6-$ = ( È̅ )𝐾 C ∆𝜖6EP
where ∆𝜖6EP is the maximum equivalent strain range calculated from the elastic analysis of under
primary and secondary loading together. ∆𝜖6EP is calculated according to Section III, Division 5,
HBB-T-1413 with 𝜈 ∗ = 0.3 for elastic analysis. 𝑆 ∗ and 𝑆̅ are stresses determined by entering the
isochronous stress-strain curves at ∆𝜖6EP and 𝐾∆𝜖6EP , respectively.
where f is the inelastic multiaxial adjustment factor determined using Section III, Division 5,
Figure HBB-T-1432-2 and triaxiality factor, T.F.
where σ’s are principals stresses at the extreme of the stress cycle.
𝐾L> is the adjustment for inelastic biaxial Poisson’s ratio determined from Section III, Division 5,
Figure HBB-T-1432-3 using 𝐾% .
1 ; 𝐾∆𝜖6EP ≤ 3𝑆6̅ /𝐸
𝐾% = × Ã∆ÙÚÛÜ ¹ ̅ /𝐸
̅
®ÈÚ
; 𝐾∆𝜖6EP > 3𝑆6
where
1.5 𝑆6 + 𝑆9; ; when only one of the extreme of the stress difference occurs at a
̅ = Ý temperature above those covered by Division 1, Subsection NB rules
3𝑆6
𝑆9; + 𝑆9: ; when both of the extreme of the stress difference occur at a
temperature above those covered by Division 1, Subsection NB rules
Here 𝑆9; 𝑎𝑛𝑑 𝑆9: are relaxation strengths associated with the temperatures at the hot and cold
extremes of the stress cycle. These values are provided above in the 740H design data. The hot
temperature condition is defined as the maximum operating temperature of the stress cycle. The
hot time is equal to the portion of service life when wall averaged temperatures exceed 425°C. The
cold temperature is defined as the colder of the two temperatures corresponding to the two stress
extremes in the stress cycle. The cold time is again equal to the portion of service life when wall
averaged temperatures exceed 425°C.
The creep strain increment per stress cycle, ∆𝜖3 is determined by entering the isochronous stress-
strain curves at 𝜎3 and maximum metal temperature for the stress cycle time, including hold times
between transients (instead of total service life). Alternatively, the creep accumulated during the
ANL-20/03 65
entire service life divided by the number of stress cycles during the entire service life can also be
used for calculating creep strain increment per stress cycle, ∆𝜖3 . We used the latter option.
The design allowable cycles, 𝑁$ is then calculated from design fatigue curve at maximum metal
temperature and using total strain range, 𝜖D , as illustrated in Figure 4-1.2.5.1 Fatigue damage
fraction, 𝐷J is then determined from the ratio between design cycles and design allowable cycles
for each cycle type and then adding them together. Figure 4-1.2.5.2 shows the equivalent strain
ranges from elastic analysis between load points A and B and between load points B and C along
two stress classification lines. Table 4-1.2.5.1 shows the details of all the relevant calculations to
determine fatigue damage fraction.
66 ANL-20/03
Figure 4-1.2.5.2. (Sample problem 1) Equivalent strain range from elastic analysis.
Creep damage evaluation is done in accordance to HBB-T1433(b)-option (a) but with one
exception as described in the design Method 1. The lower bound stress 𝑆:< is taken as 1.0𝜎3 , rather
than the 1.25𝜎3 specified in the Code. First, stress relaxation profile is determined by entering the
isochronous stress-strain curves at a strain level equal to 𝜖D and at hold-time temperature and
determining the corresponding stress levels at varying times. However, this stress relaxation
process should not be permitted to a stress level less than 𝑆:< . This stress relaxation procedure
results in a stress-time history similar to that illustrated in Figure 4-1.2.5.3 Using the stress-time
history and hold-time temperature during the cycle creep damage fraction can be calculated
according to the illustration in Figure 4-1.2.5.4 For creep damage fraction calculation, we only
considered the start-up/shut-down service load cycle. The time duration of the cloud events is
already included in the start-up/shut-down service load cycle. Creep damage is not expected during
night time. Tables 4-1.2.5.2 and 4-1.2.5.3 show the details of determining creep damage fraction,
𝐷3 from stress relaxation profile.
To determine whether the design passes the creep-fatigue damage check, the fatigue damage
fraction, 𝐷J and creep damage fraction, 𝐷3 are plotted on creep-fatigue interaction diagram as
shown in Figure 4-1.2.5.5 If the (𝐷J , 𝐷3 ) point falls inside the creep-fatigue damage envelop the
design passes. As seen in Figure 4-1.2.5.5, the (𝐷J , 𝐷3 ) points fall inside the creep-fatigue damage
envelop which means the design passes for creep-fatigue damage check.
ANL-20/03 67
At OD on stress classification line-1 At OD on Stress classification line-2
shown in Figure 4-1.2.5.2 shown in Figure 4-1.2.5.2
Start-up/shut- Cloud event Start-up/shut- Cloud event
down cycle down cycle
𝑇 6EP 770°C 770°C 767.3°C 767.3°C
Hot temperature 770°C 770°C 767.3°C 767.3°C
Cold 30°C 550°C 30°C 550°C
temperature
Hot time 12hr*(30*365) 0 12hr*(30*365) 0
=131400 hr =131400 hr
Cold time 12hr*(30*365) 8min*(30*365) 12hr*(30*365) 8min*(30*365)
=131400 hr =1460 hr =131400 hr =1460 hr
𝑆9; 107.2 MPa 354.8 MPa 110.5 MPa 358.9 MPa
𝑆9: Not required 153.3 MPa Not required 158.4 MPa
𝑆6 at 𝑇 6EP 241.0 Not required 243.2 Not required
3𝑆6 ̅ 468.7 MPa 508.1 MPa 475.3 MPa 517.3 MPa
∆𝜖6EP 0.0356 % 0.0360 % 0.0387 % 0.0376 %
𝐾 1 (no peak stress) 1 (no peak stress) 1 (no peak stress) 1 (no peak stress)
𝐾∆𝜖6EP 0.0356 % 0.0360 % 0.0387 % 0.0376 %
𝐸 171700 MPa 171700 MPa 171943 MPa 171943 MPa
3𝑆6 ̅ /𝐸 0.278 % 0.296 % 0.276 % 0.301 %
𝐾% 1 1 1 1
𝐾L> 1 1 1 1
𝐾L 1 1 1 1
È∗ 1 1 1 1
È̅
∆𝜖6-$ 0.0356 % 0.0360 % 0.0387 % 0.0376 %
∆𝜖3 1.57e-9% 0 1.27e-9% 0
𝜖D 0.0356 % 0.0360 % 0.0387 % 0.0376 %
Design 5086519 5047174 4789440 4892815
allowable
cycles, 𝑁$
Design cycles, 𝑛 30*365=10950 30*365*5=54750 30*365=10950 30*365*5=54750
Fatigue damage 0.0130 0.0135
fraction, 𝑫𝒇
Table 4-1.2.5.1. (Sample problem 1) Sample calculation of fatigue damage fraction, 𝐷J according to Method 1.
68 ANL-20/03
Figure 4-1.2.5.3. Illustration of determining stress relaxation profile from isochronous stress-strain curves for creep
damage calculation in Method 1.
Figure 4-1.2.5.4. (Sample problem 1) Illustration of calculating creep damage fraction in Method 1.
ANL-20/03 69
At OD on stress At ID on stress At OD on stress At ID on stress
classification line-1 classification line-1 classification line-2 classification line-2
shown in Figure 4- shown in Figure 4- shown in Figure 4- shown in Figure 4-
1.2.5.2 1.2.5.2 1.2.5.2 1.2.5.2
Start-up/shut-down Start-up/shut-down Start-up/shut-down Start-up/shut-down
cycle cycle cycle cycle
∆𝜖D 0.0356 % 0.0465 % 0.0377 % 0.0537 %
𝑇;i 770°C 740°C 767.3°C 737.3°C
𝑆:< = 𝜎3 20.53 20.53 20.53 20.53
𝐾 > (Table HBB-T- 0.9 0.9 0.9 0.9
1411-1) for elastic
analysis
Creep damage 2.02e-5 2.04e-5 2.53e-5 3.23e-5
fraction per cycle
(from Table 4-
1.2.5.3)
Design cycles, 𝑛 30*365=10950 30*365=10950 30*365=10950 30*365=10950
Creep damage 0.22 0.22 0.28 0.35
fraction, 𝑫𝒄
Table 4-1.2.5.2. (Sample problem 1) Sample creep damage fraction, 𝐷3 calculation according to Method 1.
70 ANL-20/03
At OD on stress classification line-1 shown in Figure 4-1.2.5.2 At OD on stress classification line-2 shown in Figure 4-
1.2.5.2
Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕 Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕
( )𝒌 ( )𝒌
(hr) (MPa) 𝑲> (hr) (hr) 𝑻𝒅 (hr) (MPa) 𝑲> (hr) (hr) 𝑻𝒅
(MPa) (MPa)
0 6.11e1 6.79e1 5.95e5 1 1.68e-6 0 6.65e1 7.39e1 4.73e5 1 2.11e-6
1 6.11e1 6.79e1 5.95e5 1 1.68e-6 1 6.65e1 7.39e1 4.73e5 1 2.11e-6
2 6.11e1 6.79e1 5.95e5 1 1.68e-6 2 6.65e1 7.39e1 4.73e5 1 2.11e-6
3 6.11e1 6.79e1 5.95e5 1 1.68e-6 3 6.65e1 7.39e1 4.73e5 1 2.11e-6
4 6.11e1 6.79e1 5.95e5 1 1.68e-6 4 6.65e1 7.39e1 4.73e5 1 2.11e-6
5 6.11e1 6.79e1 5.95e5 1 1.68e-6 5 6.65e1 7.39e1 4.73e5 1 2.11e-6
6 6.11e1 6.79e1 5.95e5 1 1.68e-6 6 6.65e1 7.39e1 4.73e5 1 2.11e-6
7 6.11e1 6.79e1 5.95e5 1 1.68e-6 7 6.65e1 7.39e1 4.73e5 1 2.11e-6
8 6.11e1 6.79e1 5.95e5 1 1.68e-6 8 6.65e1 7.39e1 4.73e5 1 2.11e-6
9 6.11e1 6.79e1 5.95e5 1 1.68e-6 9 6.65e1 7.39e1 4.73e5 1 2.11e-6
10 6.11e1 6.79e1 5.95e5 1 1.68e-6 10 6.65e1 7.39e1 4.73e5 1 2.11e-6
11 6.11e1 6.79e1 5.95e5 1 1.68e-6 11 6.65e1 7.39e1 4.73e5 1 2.11e-6
12 6.11e1 6.79e1 5.95e5 12 6.65e1 7.39e1 4.73e5
Creep damage fraction per cycle 2.02e-5 Creep damage fraction per cycle 2.53e-5
At ID on stress classification line-1 shown in Figure 4-1.2.5.2 At ID on stress classification line-2 shown in Figure 4-
1.2.5.2
Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕 Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕
( )𝒌 ( )𝒌
(hr) (MPa) 𝑲> (hr) (hr) 𝑻𝒅 (hr) (MPa) 𝑲> (hr) (hr) 𝑻𝒅
(MPa) (MPa)
0 8.11e1 9.01e1 5.88e5 1 1.70e-6 0 9.37e1 1.04e2 3.72e5 1 2.69e-6
1 8.11e1 9.01e1 5.88e5 1 1.70e-6 1 9.37e1 1.04e2 3.72e5 1 2.69e-6
2 8.11e1 9.01e1 5.88e5 1 1.70e-6 2 9.37e1 1.04e2 3.72e5 1 2.69e-6
3 8.11e1 9.01e1 5.88e5 1 1.70e-6 3 9.37e1 1.04e2 3.72e5 1 2.69e-6
4 8.11e1 9.01e1 5.88e5 1 1.70e-6 4 9.37e1 1.04e2 3.72e5 1 2.69e-6
5 8.11e1 9.01e1 5.88e5 1 1.70e-6 5 9.37e1 1.04e2 3.72e5 1 2.69e-6
6 8.11e1 9.01e1 5.88e5 1 1.70e-6 6 9.36e1 1.04e2 3.72e5 1 2.69e-6
7 8.11e1 9.01e1 5.88e5 1 1.70e-6 7 9.36e1 1.04e2 3.72e5 1 2.69e-6
8 8.11e1 9.01e1 5.88e5 1 1.70e-6 8 9.36e1 1.04e2 3.72e5 1 2.69e-6
9 8.11e1 9.01e1 5.88e5 1 1.70e-6 9 9.36e1 1.04e2 3.72e5 1 2.69e-6
10 8.11e1 9.01e1 5.88e5 1 1.70e-6 10 9.36e1 1.04e2 3.72e5 1 2.69e-6
11 8.11e1 9.01e1 5.88e5 1 1.70e-6 11 9.36e1 1.04e2 3.72e5 1 2.69e-6
12 8.11e1 9.01e1 5.88e5 12 9.36e1 1.04e2 3.72e5
Creep damage fraction per cycle 2.04e-5 Creep damage fraction per cycle 3.23e-5
Table 4-1.2.5.3. (Sample problem 1) Sample calculation to determine creep damage fraction, 𝐷3 per cycle from
stress-time history according to Method 1.
ANL-20/03 71
Figure 4-1.2.5.5. (Sample problem 1) Illustration of creep-fatigue design check. Plotted data are results from
analysis according to Method 1.
As buckling is not expected under the salt pressure and thermal loading considered in this problem
and because wind load is not considered, this design check was not performed.
As buckling is not expected under the salt pressure and thermal loading considered in this problem,
this design check was not performed.
4-1.3.2. Step-2: Transient elastic thermo-mechanical analysis for each service load
case and stress classification
72 ANL-20/03
4-1.3.5. Step-5: Creep-fatigue damage check
Method 2 is applicable only if the total stress intensity (𝑃 + 𝑄 + 𝐹) remains less than 𝑆. for all
service loading and if peak stresses are minimal.
where D is the total creep-fatigue damage and the first and second terms on the left side are fatigue
damage, 𝐷J and creep damage, 𝐷3 , respectively. In the fatigue damage term, (𝑛)Ì is the number
of repetitions of cycle type j and (𝑁$ )Ì is the number of design allowable cycles for respective
cycle type; while in the creep damage term, (𝑇$ )§ is the allowable time duration for a given stress
at the maximum temperature occurring in the time interval k and (∆𝑡)§ is the duration of the time
interval k.
The design allowable cycles for fatigue damage is determined by entering fatigue curves at total
strain range, 𝜖D . Total strain range, 𝜖D is calculated using equation HBB-T-1432-16:
∆𝜀 = ∆𝜖B + ∆𝜖C
where ∆𝜖B is the maximum equivalent strain range calculated from the elastic analysis of under
primary and secondary loading together, according to Section III, Division 5, HBB-T-1413. ∆𝜖C
is the creep strain increment per stress cycle. ∆𝜖C can be determined by entering the isochronous
stress-strain curves at the O’Donnell-Porowski core stress, 𝜎3 (determined in Method 1, Step 4)
and maximum metal temperature for the stress cycle time, including hold times between transient
(instead of total service life). Alternatively, ∆𝜖C can be calculated by dividing the creep strain
accumulated during the entire service life by the number of stress cycles during the entire service
life. We used the latter option.
Creep damage for each service load cycle is evaluated from the von Mises stress profile,
determined from elastically calculated stresses, versus time profile for this load cycle. Using the
stress-time profile and the hold time temperature, 𝑇;i during the cycle, creep damage fraction can
be calculated according to the illustration in Figure 4-1.3.5.1 As mentioned before, we only
considered the start-up/shut-down service load cycle for creep damage fraction calculation.
Tables 4-1.3.5.1 and 4-1.3.5.2 show few sample calculations of determining creep damage
fraction, 𝐷3 and fatigue damage fraction, 𝐷J , respectively, according to Method 2. Similar to
Method 1, Method 2 also uses creep-fatigue interaction diagram to determine whether a design
passes creep-fatigue damage check. Comparing (𝐷J , 𝐷3 ) with the damage envelop in creep-fatigue
interaction diagram, as shown in Figure 4-1.3.5.2, the design is found to be passed according to
Method 2.
ANL-20/03 73
Figure 4-1.3.5.1. (Sample problem 1) Illustration of calculating creep damage fraction in Method 2.
At OD on stress classification line-1 shown in Figure 4- At OD on stress classification line-2 shown in Figure 4-
1.2.5.2 1.2.5.2
𝑇;i = 770°C 𝑇;i = 767.3°C
𝑆. 𝑎𝑡 𝑇;i = 490 MPa (Method 2 is applicable!) 𝑆. 𝑎𝑡 𝑇;i = 492.3 MPa (Method 2 is applicable!)
Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕 Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕
( )𝒌 ( )𝒌
(hr) (MPa) 𝑲> (= 𝟎. 𝟗) (hr) (hr) 𝑻𝒅 (hr) (MPa) 𝑲> (= 𝟎. 𝟗) (hr) (hr) 𝑻𝒅
(MPa) (MPa)
0 0 - - - - 0 0 - - - -
1 7.88 8.76 1.47e7 1 6.80e-8 1 9.05 10.06 1.37e7 1 7.23e-8
2 57.17 63.52 7.55e5 10 1.32e-5 2 59.93 66.59 6.95e5 10 1.44e-5
10 57.17 - - - - 10 59.93 - - - -
11 7.88 8.76 1.47e7 1 6.80e-8 11 9.05 10.06 1.37e7 1 7.23e-8
12 0 - - - - 12 0 - - - -
Creep damage fraction per cycle 1.33e-5 Creep damage fraction per cycle 1.45e-5
Creep damage fraction, 𝑫𝒄 0.15 Creep damage fraction, 𝑫𝒄 0.16
At ID on stress classification line-1 shown in Figure 4-1.2.5.2 At ID on stress classification line-2 shown in Figure 4-
1.2.5.2
𝑇;i = 740°C 𝑇;i = 737.3°C
𝑆. 𝑎𝑡 𝑇;i = 512 MPa (Method 2 is applicable!) 𝑆. 𝑎𝑡 𝑇;i = 513 MPa (Method 2 is applicable!)
Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕 Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕
( )𝒌 ( )𝒌
(hr) (MPa) 𝑲> (= 𝟎. 𝟗) (hr) (hr) 𝑻𝒅 (hr) (MPa) 𝑲> (= 𝟎. 𝟗) (hr) (hr) 𝑻𝒅
(MPa) (MPa)
0 0 - - - - 0 0 - - - -
1 9.72 10.80 1.58e7 1 6.33e-8 1 10.55 11.7 1.60e7 1 6.25e-8
2 74.01 82.23 8.16e5 10 1.23e-5 2 86.03 95.59 5.25e5 10 1.90e-5
10 74.01 - - - - 10 86.03 - - - -
11 9.72 10.8 1.58e7 1 6.33e-8 11 10.55 11.7 1.60e7 1 6.25e-7
12 0 - - - - 12 0 - - - -
Creep damage fraction per cycle 1.36e-5 Creep damage fraction per cycle 1.91e-5
Creep damage fraction, 𝑫𝒄 0.15 Creep damage fraction, 𝑫𝒄 0.21
Table 4-1.3.5.1. (Sample problem 1) Sample calculation to determine creep damage fraction, 𝐷3 per cycle from
stress-time history according to Method 2.
74 ANL-20/03
At OD on stress classification line-1 shown At OD on Stress classification line-2 shown
in Figure 4-1.2.5.2 in Figure 4-1.2.5.2
Start-up/shut-down Cloud event Start-up/shut-down Cloud event
cycle cycle
𝑇 6EP 770°C 770°C 767.3°C 767.3°C
∆𝜖B 0.0356 % 0.0360 % 0.0377 % 0.0376 %
∆𝜖C 1.10e-6% 0 1.10e-6% 0
∆𝜀 0.0356 % 0.0360 % 0.0377 % 0.0376 %
Design 5086519 5047174 4883326 4892815
allowable
cycles, 𝑁$
Design cycles, 𝑛 30*365=10950 30*365*5=54750 30*365=10950 30*365*5=54750
Fatigue damage 0.0130 0.0134
fraction, 𝑫𝒇
Table 4-1.3.5.2. (Sample problem 1) Sample calculations of determining fatigue damage fraction, 𝐷J according to
Method 2.
Figure 4-1.3.5.2. (Sample problem 1) Illustration of creep-fatigue design check. Plotted data are results from
analysis according to Method 2.
This method is applicable only if the elastically-calculated stresses remain below the material yield
stress, 𝑆. . In the discussion of design calculation based on Method 2, it is shown that the
elastically-calculated stress for this sample problem is always less than 𝑆. and therefore Method 3
is applicable.
ANL-20/03 75
For primary load design, Method 3 uses the same procedures in Method 1 which is based on elastic
analysis. For ratcheting and creep-fatigue design checks, however, this method uses inelastic
analysis where material’s constitutive response is described by an elastic-creep model. The
description of the elastic-creep material model is provided above. Design calculations related to
ratcheting and creep-fatigue damage are discussed for Method 3.
4-1.4.2.
Step-2a: Transient elastic thermo-mechanical analysis for each service load case
and stress classification (for primary load design check)
We used MOOSE (Multiphysics Object Oriented Simulation Environment), an open source finite
element solver to perform the transient elastic-creep thermo-mechanical analyses under the loading
conditions mentioned in Step 1. The analysis was repeated until a steady state cyclic response was
achieved.
To determine ratcheting strain Method 3 requires to run the analysis using elastic-creep material
C
model, described above, and monitor the maximum effective strain, M® 𝜀: 𝜀 at the beginning and
end of the cycle. The criterion is that the ratcheting strain does not exceed 10% at any point of the
structure for base metal. Figure 4-1.4.4 plots the maximum effective strain at the critical tube
location as a function of cycle count. Extrapolating the maximum effective strain out to design life
of the tube, i.e. 30 years (=30*365 cycles), gives the ratcheting strain of 0.00565% which is less
than 10%.
76 ANL-20/03
Figure 4-1.4.4. (Sample problem 1) Maximum ratcheting strain in the structure versus number of cycles determined
from elastic-creep thermo-mechanical analysis.
Once steady cyclic response was achieved in the analysis, the temperature, stress, strain, time
history for a single cycle of the periodic loading were extracted. To determine fatigue damage
fraction, the effective strain range, ∆ε was first computed from the strain history according to
Section III, Division 5, HBB-T-1413 with 𝜈 ∗ = 0.5 for inelastic analysis. Fatigue damage fraction,
𝐷J was then calculated from ∆ε using rainflow counting and Miner’s rule. Figure 4-1.4.5.1 plots
temperature, von Mises stress, and effective strain range profiles at four critical locations of the
tube after a steady cyclic response was achieved. Table 4-1.4.5.1 shows details of the fatigue
damage fraction calculation according to Method 3. The von Mises effective stress, 𝜎%JJ (𝑡) was
used to determine the creep damage fraction. Figure 4-1.4.5.2 illustrates the method of creep
damage fraction, 𝐷3 calculation and Table 4-1.4.5.2 reports details of the calculation for four
critical locations in the structure. All four sets of (𝐷J , 𝐷3 ) fall inside the damage envelop in the
creep-fatigue interaction diagram, as shown in Figure 4-1.4.5.3, which means according to Method
3 the design passes creep-fatigue damage check.
Note, according to Table HBB-1411-1 in Section III, Division 5 of ASME Code, a design margin
(i.e, 𝐾 > ) is applied to the effective stress while determining the allowable rupture time from the
design rupture table. ASME Code recommends to use 𝐾 > = 0.9 for elastic analysis and 𝐾 > = 0.67
for inelastic analysis. However, the history behind the ASME stress factors is somewhat murky
and given the lower consequences of failure for CSP systems, the design method developed here
uses 𝐾 > = 0.9 for all analysis methods, including the inelastic analysis.
ANL-20/03 77
Figure 4-1.4.5.1. (Sample problem 1) Temperature, effective strain range, and von Mises stress profiles at four
critical location of the tube after a steady cyclic response is achieved in the elastic-creep thermo-mechanical
analysis.
78 ANL-20/03
Figure 4-1.4.5.2. Illustration of calculating creep damage fraction in Method 3.
ANL-20/03 79
Figure 4-1.4.5.3. (Sample problem 1) Illustration of creep-fatigue design check. Plotted data are results from
analysis according to Method 3.
For sample problem 2 we considered a tube in an external tubular receiver. The receiver has
parabolic reflectors at the back of the tubes which help reduce the variation in the circumferential
heat flux distribution on the tube. Figure 4-2.1.1 shows the schematic of the tubular receiver. The
tube is 10.5 m long, 42.2 mm diameter, and 1 mm thick. Heat flux on the tube is non uniform along
both the length and circumference. We considered only one type of cycle for this problem which
represent heat flux on the day of spring equinox. Figure 4-2.1.2 plots the loading profiles of
maximum flux incident, salt inlet and outlet temperature, and salt pressure during day (10 hrs).
Table 4-2.1.1 lists values of flux incident loading at noon (i.e. time = 5hr in Figure 4-2.1.2) for
different locations on the tube outer surface. To determine the wind load on the receiver for time-
independent buckling check, we considered the external receiver located in an open terrain in the
U.S. and 100 m above the ground.
The salt considered for this problem is MgCl2/KCl (mole: 32/68%) binary molten salt. The mass
flow rate of the salt is 44.5 kg/s. Salt inlet and outlet temperatures are 700°C and 720°C,
80 ANL-20/03
respectively. To determine the convective heat transfer coefficient between salt and metal
Gnielinski correlation29 for turbulent flow (forced convection) in tubes was used.
f
ú ü(Œ%ý ‰Blll)Á9
𝑁𝑢ù = û
Ò
Ó
; for 3000 < 𝑅𝑒ù < 5 × 10þ and 0.5 < 𝑃𝑟 < 2000
f Ó
BºBC.“ ú ü (Á9 ¥‰B)
û
where 𝑁𝑢ù is the Nusselt number, 𝑅𝑒ù is the Reynolds number, 𝑃𝑟 is the Prandtl number, and 𝑓
is the Darcy friction factor which can be obtained from the following equation.
All the required properties of salt, as listed in Table 4-2.1.2, to determine convective heat transfer
coefficient were taken from30. Note the pressure loading profile shown in Figure 4-2.1.2 was
selected based on the pressure loss due to friction and elevation. We assumed a surface roughness
of 0.049 mm to determine the frictional pressure loss. In the heat transfer calculation for this
sample problem, we also considered heat loss due to radiation and natural convection of air. The
ÿ
convective heat transfer coefficient of air is assumed to be 14 6Ó Ã. For structural boundary
condition, we considered the tube can freely expand both in radial and axial direction but no
warping in the axial direction. The design life of the tube is considered to be 4.4 years. Design
calculations according all three methods are provided below.
29 V. Gnielinski, “Neue Gleichungen für den Wärme-und den Stoffübergang in turbulent durchströmten Rohren und
Kanälen” Forschung im Ingenieurwesen A, 41(1), pp. 8-16, 1975.
30 Xu, X, X. Wang, P. Li, Y. Li, Q. Hao, B. Xiao, H. Elsentriecy, and D. Gervasio, “Experimental
Test of Properties of KCl–MgCl2 Eutectic Molten Salt for Heat Transfer and Thermal Storage
Fluid in Concentrated Solar Power Systems” Journal of Solar Energy Engineering, 140(5), pp.
051011, 2018.
ANL-20/03 81
Figure 4-2.1.2. (Sample problem 2) Loading profiles of maximum flux incident, salt inlet and outlet temperature,
and salt pressure during day. Only one type of cycle is considered. Receiver operation time per day is 10 hours.
Loading profiles shown are only for the tube considered for design study.
82 ANL-20/03
Angular location, θ
along the -0.438 0.988 0.955 0.922 0.889 0.856 0.823 0.790
axis of 0.000 0.976 0.943 0.911 0.878 0.846 0.813 0.781
the tube, 0.437 0.972 0.940 0.907 0.875 0.842 0.810 0.778
m 0.875 0.980 0.947 0.915 0.882 0.849 0.817 0.784
1.312 0.993 0.960 0.927 0.894 0.861 0.828 0.794
1.750 1.000 0.967 0.933 0.900 0.867 0.833 0.800
2.188 0.994 0.961 0.928 0.895 0.861 0.828 0.795
2.625 0.969 0.937 0.904 0.872 0.840 0.808 0.775
3.062 0.918 0.887 0.857 0.826 0.796 0.765 0.734
3.500 0.834 0.806 0.778 0.751 0.723 0.695 0.667
3.937 0.714 0.690 0.666 0.643 0.619 0.595 0.571
4.375 0.568 0.549 0.530 0.511 0.492 0.473 0.454
4.812 0.414 0.400 0.386 0.373 0.359 0.345 0.331
5.250 0.272 0.263 0.254 0.245 0.236 0.227 0.218
Table 4-2.1.1. (Sample problem 2) Flux incident loading at noon.
Table 4-2.1.2. (Sample problem 2) Temperature dependent properties of MgCl2/KCl (mole: 32/68%) binary molten
salt taken from30.
ANL-20/03 83
4-2.2. Design calculations based on Method 1
4-2.2.2. Step-2: Transient elastic thermo-mechanical analysis for each service load
case and stress classification
We used MOOSE (Multiphysics Object Oriented Simulation Environment), an open source finite
element solver to perform the elastic thermo-mechanical analyses. We classify pressure as primary
load and temperature gradient as secondary load. There is no peak load.
Figure 4-2.2.3.1 shows the heat flux and tube outer wall temperature distribution at noon. Figure
4-2.2.3.2 shows through thickness elastic stress components at noon at a critical location of the
tube.
Figure 4-2.2.3.1. (Sample problem-2) Contour plot of flux incident on tube and tube outer wall temperature at noon.
84 ANL-20/03
Figure 4-2.2.3.2. (Sample problem-2) Through thickness temperature and elastic stresses at a critical location of the
tube at noon.
Maximum primary load occurs at noon. Table 4-2.2.3 reports details of the primary load checks.
First, all the stress components were linearized to divide into membrane and bending components
along the stress classification line. The membrane and bending stress components were then used
to determine the stress intensities in Table 4-2.2.3 As indicated in the table the design passes both
the criteria in primary load checks.
Design method 1 uses the O’Donnell-Porowski approach, described in Section III, Division 5,
HBB-1332 for ratcheting checks. In this approach, an effective creep stress parameter, Z is
determined from a primary stress parameter, X and a secondary stress parameter, Y as shown in
Figure 4-1.2.4.1 The effective creep stress parameter is used to calculate the effective creep stress
which is then used to determine the ratcheting creep strain using isochronous stress-strain curves.
The definition of X and Y are
𝑃¤
𝑋 = (𝑃: + ) /𝑆
𝐾D 6EP .:
𝑌 = (𝑄Œ )6EP /𝑆.:
ANL-20/03 85
where,
Á
(𝑃: + ÃÂ)6EP = the maximum value of the primary stress intensity, adjusted for bending via 𝐾D ,
Ä
during the cycle being evaluated.
(𝑄Œ )6EP = the maximum range of the secondary stress intensity during the cycle being considered
𝑆.: = is the 𝑆. value corresponds to the lower of the wall averaged temperature for the stress
extremes defining secondary stress range, 𝑄Œ .
𝐾D = (𝐾 + 1)/2
𝐾 is 1.5 for across-the-wall bending of shell structures or rectangular sections, see HBB-3223 (c)
(6) in Section III Division 5. Once Z is found, effective core, 𝜎3 stress is determined from
𝜎3
𝑍=
𝑆.:
It should be noted that, the average wall temperature at one of the stress extremes defining the
secondary stress intensity range must be below the temperature listed in Section III, Division 5,
HBB-T-1323, given as 600℃ for 740H in the description of the design Method 1. The creep
ratcheting strain increment for a load cycle is evaluated by entering the isochronous stress strain
curves at the maximum wall temperature and effective core, 𝜎3 stress during the load cycle with
the stress held constant for the entire service life. An example of creep ratcheting strain
determination is shown in Figure 4-1.2.4.2 Table 4-2.2.4 provides all the calculation details of the
ratcheting design check. As indicated in the table, the design passes the ratcheting check.
86 ANL-20/03
According to Section III, Division 5, a design is acceptable if the creep and fatigue damage satisfy
the following relation:
𝑛 ∆𝑡
&( )Ì + &( )§ ≤ 𝐷
𝑁$ 𝑇$
Ì §
where D is the total creep-fatigue damage and the first and second terms on the left side are fatigue
damage, 𝐷J and creep damage, 𝐷3 , respectively. In the fatigue damage term, (𝑛)Ì is the number
of repetitions of cycle type j and (𝑁$ )Ì is the number of design allowable cycles for respective
cycle type; while in the creep damage term, (𝑇$ )§ is the allowable time duration for a given stress
at the maximum temperature occurring in the time interval k and (∆𝑡)§ is the duration of the time
interval k.
The design allowable cycles for fatigue damage is determined by entering fatigue curves at total
strain range, 𝜖D . Total strain range, 𝜖D is calculated using equation HBB-T-1432-16:
𝜖D = 𝐾L ∆𝜖6-$ + 𝐾∆𝜖3
where 𝐾 is the local geometric concentration or equivalent stress concentration factor determined
by dividing effective primary plus secondary plus peak stress divided by the effective primary plus
secondary stress, 𝐾L is the multiaxial plasticity and Poisson ratio adjustment factor, ∆𝜖3 is the
creep strain increment, and ∆𝜖6-$ is the modified maximum equivalent strain range.
𝑆∗ C
∆𝜖6-$ = ( ̅ )𝐾 ∆𝜖6EP
𝑆
where ∆𝜖6EP is the maximum equivalent strain range calculated from the elastic analysis of under
primary and secondary loading together. ∆𝜖6EP is calculated according to Section III, Division 5,
HBB-T-1413 with 𝜈 ∗ = 0.3 for elastic analysis. 𝑆 ∗ and 𝑆̅ are stresses determined by entering the
isochronous stress-strain curves at ∆𝜖6EP and 𝐾∆𝜖6EP , respectively.
where f is the inelastic multiaxial adjustment factor determined using Section III, Division 5,
Figure HBB-T-1432-2 and triaxiality factor, T.F.
|𝜎B + 𝜎C + 𝜎® |
𝑇. 𝐹. =
1
[(𝜎B − 𝜎C )C + (𝜎C − 𝜎® )C + (𝜎® − 𝜎B )C ]
√2
ANL-20/03 87
where σ’s are principals stresses at the extreme of the stress cycle.
𝐾L> is the adjustment for inelastic biaxial Poisson’s ratio determined from Section III, Division 5,
Figure HBB-T-1432-3 using 𝐾% .
1 ; 𝐾∆𝜖6EP ≤ 3𝑆6̅ /𝐸
𝐾% = › 𝐾∆𝜖6EP 𝐸 ̅ /𝐸
̅ ; 𝐾∆𝜖6EP > 3𝑆6
3𝑆6
where
1.5 𝑆6 + 𝑆9; ; when only one of the extreme of the stress difference occurs at a
̅ = Ý temperature above those covered by Division 1, Subsection NB rules
3𝑆6
𝑆9; + 𝑆9: ; when both of the extreme of the stress difference occur at a
temperature above those covered by Division 1, Subsection NB rules
Here 𝑆9; 𝑎𝑛𝑑 𝑆9: are relaxation strengths associated with the temperatures at the hot and cold
extremes of the stress cycle. These values are provided above in the 740H design data. The hot
temperature condition is defined as the maximum operating temperature of the stress cycle. The
hot time is equal to the portion of service life when wall averaged temperatures exceed 425°C. The
cold temperature is defined as the colder of the two temperatures corresponding to the two stress
extremes in the stress cycle. The cold time is again equal to the portion of service life when wall
averaged temperatures exceed 425°C.
The creep strain increment per stress cycle, ∆𝜖3 is determined by entering the isochronous stress-
strain curves at 𝜎3 and maximum metal temperature for the stress cycle time, including hold times
between transients (instead of total service life). Alternatively, the creep accumulated during the
entire service life divided by the number of stress cycles during the entire service life can also be
used for calculating creep strain increment per stress cycle, ∆𝜖3 . We used the latter option.
The design allowable cycles, 𝑁$ is then calculated from design fatigue curve at maximum metal
temperature and using total strain range, 𝜖D . Fatigue damage fraction, 𝐷J is then determined from
the ratio between design cycles and design allowable cycles for each cycle type and then adding
them together. Note, we considered only one type of cycle in the sample problem 2. Table 4-2.3.5.1
shows the details of all the relevant calculations to determine fatigue damage fraction.
Creep damage evaluation is done in accordance to HBB-T1433(b)-option (a) but with one
exception as described in the design Method 1. The lower bound stress 𝑆:< is taken as 1.0𝜎3 , rather
than the 1.25𝜎3 specified in the Code. First, stress relaxation profile is determined by entering the
isochronous stress-strain curves at a strain level equal to 𝜖D and at hold-time temperature and
determining the corresponding stress levels at varying times. However, this stress relaxation
process should not be permitted to a stress level less than 𝑆:< . This stress relaxation procedure
results in a stress-time history. Using the stress-time history and hold-time temperature during the
cycle creep damage fraction can be calculated according to the illustration in Figure 4-1.2.5.4.
Tables 4-2.3.5.2 and 4-2.3.5.3 show the details of determining creep damage fraction, 𝐷3 from
stress relaxation profile.
88 ANL-20/03
To determine whether the design passes the creep-fatigue damage check, the fatigue damage
fraction, 𝐷J and creep damage fraction, 𝐷3 are plotted on creep-fatigue interaction diagram as
shown in Figure 4-2.3.5.1 As indicated in the figure, the (𝐷J , 𝐷3 ) points fall inside the creep-fatigue
damage envelop which means the design passes for creep-fatigue damage check.
ANL-20/03 89
At OD on stress classification line shown in Figure 4- At ID on stress classification line shown in Figure 4-
2.2.3.2 2.2.3.2
Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕 Time 𝑺𝒌 𝑺𝒌 (𝑻𝒅 )𝒌 ∆𝒕𝒌 ∆𝒕
( )𝒌 ( )𝒌
(hr) (MPa) 𝑲> (hr) (hr) 𝑻𝒅 (hr) (MPa) 𝑲> (hr) (hr) 𝑻𝒅
(MPa) (MPa)
0 96.28 106.98 16341 1 6.12e-05 0 111.48 123.87 68686 1 1.46e-05
1 96.25 106.94 16378 1 6.11e-05 1 111.47 123.86 68718 1 1.46e-05
2 96.20 106.89 16424 1 6.09e-05 2 111.46 123.84 68782 1 1.45e-05
3 96.17 106.86 16452 1 6.08e-05 3 111.46 123.84 68782 1 1.45e-05
4 96.13 106.81 16498 1 6.06e-05 4 111.45 123.83 68814 1 1.45e-05
5 96.09 106.77 16535 1 6.05e-05 5 111.44 123.82 68846 1 1.45e-05
6 96.06 106.73 16573 1 6.03e-05 6 111.43 123.81 68878 1 1.45e-05
7 96.02 106.69 16610 1 6.02e-05 7 111.43 123.81 68878 1 1.45e-05
8 95.98 106.64 16657 1 6.00e-05 8 111.42 123.80 68910 1 1.45e-05
9 95.94 106.60 16695 1 5.99e-05 9 111.41 123.79 68942 1 1.45e-05
10 95.91 106.57 16723 10 111.41 123.79 68942
Creep damage fraction per cycle 6.05e-4 Creep damage fraction per cycle 1.45e-4
Table 4-2.3.5.3. (Sample problem-2) Sample calculation of creep damage fraction, 𝐷3 per cycle from stress-time
history according to Method 1.
Figure 4-2.3.5.1. (Sample problem-2) Illustration of creep-fatigue design check. Plotted data are results from
analysis according to Method 1.
Since this is a tube in an external receiver, wind load should be considered in the time-independent
buckling check. To determine the design wind load on the receiver tube, we first determine the
wind pressure along the wind direction. We adopted the method provided in ASME STS-1-2016
Steel Stacks31 for calculating the wind pressure. According to Equation 4-4 in STS-1-2016, the
deign wind pressure along the wind direction is,
31 American Society of Mechanical Engineers Code for Steel Stacks, ASME STS-1-2016.
90 ANL-20/03
𝑞S = 0.00256𝑉 C 𝐼𝐾SD 𝐾S 𝑝𝑠𝑓
where 𝑉 is the wind speed in mph, 𝐼 is the importance factor, 𝐾SD is the topographic factor, and 𝐾S
is the velocity pressure exposure coefficient evaluated at height, 𝑧.
Importance factor, 𝐼 can be determined according to Tables I-2 and I-3 in STS-1-2016 Mandatory
Appendix I. As a part of a power generation facility, the CSP receiver falls under category IV in
Table I-2. From Table I-3 for category IV structure, the importance factor is found to be 1.15.
Considering the CSP tower located in an open terrain, 𝐾SD is 1 according to Fig. I-2 in STS-1-2016
Mandatory Appendix I. 𝐾S is determined from Table I-4 in STS-1-2016 Mandatory Appendix I.
This table provides the values of 𝐾S as a function of height, 𝑧 and the exposure category of the
structure. According to Paragraph 4.3.3.4 of STS-1-2016, the exposure category for a structure
located in an open terrain is C. As mentioned in the problem definition, the CSP receiver stands
on a 100 m tall tower which makes the value of 𝐾S to vary from 1.61 at the bottom of the receiver
to 1.66 at the top of the receiver.
A structure under wind load is analyzed by computing the total force exerted on the structure from
the drag coefficient for the whole structure and the design wind pressure. Then, use the computed
force to determine stresses on the structure. However, this method cannot be applied for analyzing
individual tube of a CSP receiver. We, therefore, used a pressure coefficient, 𝐶# to determine the
distribution of wind pressure around the receiver and use the maximum positive and negative
values of the distributed wind pressure as wind load on the tube. The value of 𝐶# depends on the
geometry and dimensions of the structure. The overall structure of the external receiver considered
for this design problem can be considered as a cylinder. The variation in 𝐶# around a cylindrical
structure can be determined from the following equation32.
𝐷 𝐷 𝐷
𝐶# = −0.54 + 0.16 ∗ + œ0.28 + 0.04 ∗ • ∗ 𝑐𝑜𝑠𝜃 + œ1.04 − 0.2 ∗ • ∗ 𝑐𝑜𝑠2𝜃
𝐿 𝐿 𝐿
𝐷 𝐷
+ œ0.36 − 0.05 ∗ • ∗ 𝑐𝑜𝑠3𝜃 − œ0.14 − 0.05 ∗ • 𝑐𝑜𝑠4𝜃
𝐿 𝐿
where 𝐷 is the diameter of the whole receiver which is 8.5 m for our design problem, 𝐿 is the
height of the receiver which should be approximately equal to the length of the tube, and 𝜃 is the
circumferential coordinate from the direction of wind. The distribution of 𝐶# along the
circumference of the external receiver considered in this sample problem is shown in Figure 4-
2.2.6.1.
32 C. Lei, and J. M. Rotter, “Buckling of anchored cylindrical shells of uniform thickness under
wind load” Engineering Structures, 41, pp. 199-208, 2012.
ANL-20/03 91
Figure 4-2.2.6.1. Distribution of wind pressure coefficient, 𝐶# around the circumference of the external cylindrical
receiver. Receiver diameter, 𝐷 = 8.5 𝑚 and height, 𝐿 = 10.5 𝑚. θ = 0 is along the direction of wind velocity.
By multiplying 𝑞S with 𝐶# , the distribution of design wind pressure around the cylindrical external
receiver can be found. Depending on the location of the tube with respect to the wind direction the
value of 𝐶# and the corresponding design wind pressure acting on the tube will vary. Considering
wind can blow from any direction, the maximum acting positive and negative design wind pressure
on the tube, according to Figure 4-2.2.6.1, are 𝑞S and -1.43𝑞S , respectively. Thus, the buckling
analysis is performed for two cases of design wind load – the maximum positive wind pressure
which pushes the tube toward the cylindrical axis of the overall receiver and the maximum negative
wind pressure which radially pulls the tube away from the overall cylindrical shape of the receiver.
For structural analysis, the design wind pressure load is applied as a pressure field on the outer
half (that faces heliostat) of the outer surface of the tube. The direction of the positive pressure
field is toward the cylindrical axis of the overall receiver.
Design wind V
𝐶# 𝐾S (𝑧)
𝐼 𝐾SD
load (MPa) mph
Case-1 1.0 Varies linearly from
1.228 ∗ 10‰l“
90.0 1.15 1.0 1.62 (bottom of the tube)
Case-2 ∗ 𝑉 C 𝐼𝐾SD 𝐾S 𝐶# -1.43 to 1.66 (top of the tube)
Table 4-2.2.6.1. Applied design wind load on the tube for time-independent buckling analysis.
According to the design methods developed for the high temperature components, time-
independent buckling can be checked by analyzing the structure assuming a constitutive response
given by the hot tensile curves and a load factor of 1.5. Note that, the load factor 1.5 is a slight
reduction from the ASME Section III Division 5 HBB factor of 1.67, reflecting the lower
consequences of failure for CSP systems. Thus, we used a load factor of 1.5 on the primary internal
pressure load due to the salt pressure and on the thermal loading due to the temperature gradient.
92 ANL-20/03
To apply the load factor on thermal loading in finite element calculations we multiplied the
coefficient of thermal expansion, CTE by 1.5. In the case of wind load, the design wind load
expression, provided in Table 4-2.2.6.1, already includes several safety factors and therefore we
did not use any additional load factor. Table 4-2.2.6.2 lists all the factors used for different types
of loads.
The constitutive response of alloy 740H described by the hot tensile or isochronous stress strain
curves can be implemented in commercial finite element software such as Abaqus by providing
the tabulated values of flow stress and corresponding plastic strain. However, this option is not yet
available in the open source finite element package we used for this project. We, therefore, used
an elastic perfectly-plastic material model already implemented into our finite element solver for
analyzing the structure for buckling checks. Figure 4-2.2.6.2 plots in red the hot tensile curves at
different temperatures for alloy 740H. Curves in black represent the elastic perfectly-plastic
material models used in the finite element calculations. As seen in Figure 4-2.2.6.2, the curves for
elastic perfectly-plastic material models are the original stress-strain curves without the hardening
portion which means use of these models would be a more conservative estimation of the buckling.
ANL-20/03 93
Figure 4-2.2.6.2. Plots in red are hot tensile curves and plots in black represent the elastic perfectly-plastic material
models used in finite element simulation for buckling analysis.
94 ANL-20/03
Figure 4-2.2.6.3. Distribution of temperature, von Mises stress, and radial (with respect to the cylindrical axis of the
whole receiver) displacement at noon under the load applied for time-independent buckling check.
As time-dependent buckling is not expected under the salt pressure and thermal loading considered
in this problem, this design check was not performed.
4-2.3.2. Step-2: Transient elastic thermo-mechanical analysis for each service load
case and stress classification
Same as in Method 1 Step 2.
ANL-20/03 95
4-2.3.4. Step-4: Ratcheting check
Method 2 is applicable only if the primary plus secondary stress intensity (𝑃 + 𝑄) remains less
than 𝑆. for all service loading and if peak stresses are minimal.
𝑛 ∆𝑡
&( )Ì + &( )§ ≤ 𝐷
𝑁$ 𝑇$
Ì §
where D is the total creep-fatigue damage and the first and second terms on the left side are fatigue
damage, 𝐷J and creep damage, 𝐷3 , respectively. In the fatigue damage term, (𝑛)Ì is the number
of repetitions of cycle type j and (𝑁$ )Ì is the number of design allowable cycles for respective
cycle type; while in the creep damage term, (𝑇$ )§ is the allowable time duration for a given stress
at the maximum temperature occurring in the time interval k and (∆𝑡)§ is the duration of the time
interval k.
The design allowable cycles for fatigue damage is determined by entering fatigue curves at total
strain range, 𝜖D . Total strain range, 𝜖D is calculated using equation HBB-T-1432-16:
∆𝜀 = ∆𝜖B + ∆𝜖C
where ∆𝜖B is the maximum equivalent strain range calculated from the elastic analysis of under
primary and secondary loading together, according to Section III, Division 5, HBB-T-1413. ∆𝜖C
is the creep strain increment per stress cycle. ∆𝜖C can be determined by entering the isochronous
stress-strain curves at the O’Donnell-Porowski core stress, 𝜎3 (determined in Method 1, Step 4)
and maximum metal temperature for the stress cycle time, including hold times between transient
(instead of total service life). Alternatively, ∆𝜖C can be calculated by dividing the creep strain
accumulated during the entire service life by the number of stress cycles during the entire service
life. We used the latter option.
Creep damage for each service load cycle is evaluated from the von Mises stress profile,
determined from elastically calculated stresses, versus time profile for this load cycle. Using the
stress-time profile and the hold time temperature, 𝑇;i during the cycle, creep damage fraction can
be calculated according to the illustration in Figure 4-2.3.5.1.
Tables 4-2.3.5.1 and 4-2.3.5.2 show example calculations of determining creep damage fraction,
𝐷3 and fatigue damage fraction, 𝐷J , respectively, according to Method 2. Comparing (𝐷J , 𝐷3 ) with
the damage envelop in creep-fatigue interaction diagram, as shown in Figure 4-2.3.5.2, the design
is found to be passed according to Method 2.
96 ANL-20/03
Figure 4-2.3.5.1. (Sample problem-2) Illustration of calculating creep damage fraction in Method 2.
ANL-20/03 97
Figure 4-2.3.5.2. (Sample problem-2) Illustration of creep-fatigue design check. Plotted data are results from
analysis according to Method 2.
4-2.4.2.
Step-2a: Transient elastic thermo-mechanical analysis for each service load case
and stress classification (for primary load design check)
We used MOOSE (Multiphysics Object Oriented Simulation Environment), an open source finite
element solver to perform the transient elastic-creep thermo-mechanical analyses under the loading
conditions mentioned in Step 1. The analysis was repeated until a steady state cyclic response was
achieved.
98 ANL-20/03
4-2.4.4. Step-4: Ratcheting check
To determine ratcheting strain Method 3 requires to run the analysis using elastic-creep material
C
model, described above, and monitor the maximum effective strain, M® 𝜀: 𝜀 at the beginning and
end of the cycle. The criterion is that the ratcheting strain does not exceed 10% at any point of the
structure for base metal. Figure 4-2.4.4 plots the maximum effective strain at the critical tube
location as a function of cycle count. Extrapolating the maximum effective strain out to design life
of the tube, i.e. 4.4 years (=4.4*365 cycles), gives the ratcheting strain of 0.745% which is less
than 10%.
Figure 4-2.4.4. (Sample problem-2) Maximum ratcheting strain in the structure versus number of cycles determined
from elastic-creep thermo-mechanical analysis.
Once steady cyclic response was achieved in the analysis, the temperature, stress, and strain-time
history for a single cycle of the periodic loading were extracted. Figure 4-2.4.5.1 plots the
temperature, and the steady cyclic effective strain range and von Mises effective stress at the
critical location of the tube. Details of the fatigue damage fraction calculation is provided in Table
4-2.4.5.1. The tube experiences negligible fatigue damage. Figure 4-2.4.5.2 illustrates creep
damage fraction evaluation from the steady cyclic von Mises effective stress profile. Detailed
calculation of creep damage fraction is provided in Table 4-2.4.5.2. Comparing (𝐷J , 𝐷3 ) with the
damage envelop in creep-fatigue interaction diagram, as shown in Figure 4-2.4.5.3, the design is
found to be passed according to Method 3.
ANL-20/03 99
Figure 4-2.4.5.1. (Sample problem-2) Temperature, effective strain range, and von Mises stress profiles at the
critical location of the tube after a steady cyclic response is achieved in the elastic-creep thermo-mechanical
analysis.
100 ANL-20/03
Figure 4-2.4.5.2. (Sample problem-2) Illustration of calculating creep damage fraction in Method 3.
Figure 4-2.4.5.3. (Sample problem-2) Illustration of creep-fatigue design check. Plotted data are results from
analysis according to Method 3.
ANL-20/03 101
4-2.4.6. Step-6: Time-independent buckling check
102 ANL-20/03
Acknowledgements
This work was funded by the U.S. Department of Energy through the Office of Energy Efficiency
and Renewable Energy, Solar Energy Technologies Office, CSP Program. We gratefully
acknowledge the guidance of Mark Lausten, the DOE project manager for this project as well as
the contributions of Robert Jetter who served as a high temperature design consultant in the
development of the design rules and corresponding material data.
ANL-20/03 103
August 2018
www.anl.gov