Fatigue User’s Guide > Validation Problems > Problem 4: Ten Simple Notched Geometries
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">   
Problem 4: Ten Simple Notched Geometries
Problem Description
The problems here are described in Further Reading (App. 16). They consist of ten simple notched geometries which were tested at various constant amplitude load levels. Four different material types were used, whose properties are given later. In some cases, the measured life was to failure while for others it was to crack initiation. Figure 14‑9 shows the various notched geometries.
Figure 14‑9 Ten Simple Notched Geometries (Dimensions in Inches)
Objectives
1. To predict the initiation life of these ten different specimens using the indicated load levels and material properties.
2. To validate MSC.Fatigue against the obtained measured test data relative to these ten specimens.
Step 1: Geometry and FEA Results
For this series of benchmark problems, the geometries were created in MSC.Patran. Due to symmetry, only half of each specimen was actually modeled. MSC.Patran FEA was used as the analysis code using plane stress elements. The MSC.Patran neutral files and MSC.Patran FEA results files can be found in the examples directory delivered with your MSC.Fatigue system. The neutral files containing the geometries are in files notchxx.out where xx = 01 - 10. The stress-strain results are in files notchxx.txt which need to be converted to binary format using the MSC.Patran utility program RESTXT. See previous examples. Files can be found in
<install_dir>/mscfatigue_files/examples
Step 2:  Material Characterization
Table 14‑11 shows the four materials used in the study of these 10 specimens. They must be input to the MSC.Fatigue materials database manager PFMAT. Follow the method used in Problem 3: Comparison to Another Code, 1237, Step 2: Material Characterization, 1238 to input the four material properties and name them MATA, MATB, MATC, and MATD.
Hint: You can set the units preference to KSI so that you do not have to convert the values in Table 14‑11 to SI units. Use the Preferences option from the PFMAT main menu.
Table 14‑11 Fatigue material properties
 
Material A:
Medium strength steel
Material B:
High strength steel
Material C:
Low strength steel
Material D:
High strength aluminum
 
Material Properties
Material
 
A
B
C
D
Strain Properties
 
 
 
 
Fatigue strength coefficient, Sf ‘(KSI)
169
267
162
241
Fatigue strength exponent, b
-0.081
-0.09
-0.11
-0.15
Fatigue ductility coefficient, Ef’
1.142
0.79
0.338
0.158
Fatigue ductility exponent, c
-0.67
-0.73
-0.48
-0.83
Cyclic strain hardening exponent, n’
0.123
0.112
0.226
0.04
Cyclic strength coefficient, k’ (KSI)
154
235
194
101
Standard error, Se
0
0
0
0
Monotonic Properties
 
 
 
 
Yield strength, Sy, (KSI)
94
158
51
78
Young’s Modulus, E (KSI)
30E3
30E3
30E3
10E3
UTS (KSI)
114
168
78
85
Step 3: Loading Histories
Loading was performed with constant amplitude load levels. The first seven specimens used a stress ratio of R = -1 (fully reversed) while the last three specimens used a R ratio of 0.1. The magnitudes of the loads are given in Table 14‑12. They were created in PTIME by making just one triangular cycle with a max. and min. value of ±20kips/in for R=-1 and entering X-Y points with max=20000, and min=2000 for R=0.1. These time histories were then scaled in FEFAT for the rest of the load levels and specimens. All MSC.Patran FEA loads are in pounds/inch and the scale factors used properly convert the time history loads from ±20kips/in. to the appropriate kip/in depending on the specimen geometry and load level. (Specimen 4 used kip loading as opposed to kip/in.)
The following keystrokes create these two simple constant amplitude load histories. This assumes you are starting in a clean, empty directory:
 
Operation
Comments
ptime
Invoke PTIME either from the system prompt or from the MSC.Fatigue Loading form.
Wave creation
Select the Wave creation option.
Filename: 20KIPR1
Give it a filename called 20KIPR1.
Description 1: Fully reversed
Give it a description.
Load type: Force
Select the load type to be Force.
Units: lbs force
Select the units to be lbs force.
OK
Accept this form.
Triangular
Create a triangular wave. Use the defaults of all values except what is described below.
Amplitude: 20000
Set the amplitude to 20,000.
OK / Finish
Accept this form and select Finish for the next action. The percentage of positive slope should be 50.
Plot an entry
Select the Plot an entry option and accept the default. Close the graphics when done.
Add an entry...
Select the Add option.
X-y time series
Choose the X-Y point creation option to create the other time history.
Filename: 20KIPR01
Give it a filename called 20KIPR01.
Description 1: Not fully reversed
Give it a description.
Load type: Force
Select the load type to be Force.
Units: lbs force
Select the units to be lbs force.
OK
Accept this form.
2000 / 20000 / 2000 / OK / OK
Enter these three values and use the RETURN key each time. Click the OK with a blank entry when done.
End
Press the End button when done defining the X-Y points.
Plot an entry
Select the Plot option and accept the default. Close the graphics when done. Then exit from PTIME.
Step 4: Setup MSC.Fatigue Job
Enter MSC.Patran and use the following operations to setup the MSC.Fatigue job for the first specimen. Accept all defaults unless otherwise specified:
 
Operation
Comments
patran
Invoke MSC.Patran (or MSC.Fatigue Pre & Post) if you have not already done so.
File / New...
Open a new database from the File pull -down menu. Call it “notch.” Set the analysis preference to MSC.Patran FEA if asked. Ignore any warning messages.
File / Import...
Import the neutral file notch01.out into the database. At this time you may wish to manipulate the model, turn labels off, or other MSC.Patran operations. When ready, go on to the next step.
Tools / FATIGUE...
(Analysis)
Invoke MSC.Patran’s FATIGUE interface by selecting it from under the Tools pull-down menu (or the Analysis application switch in MSC.Fatigue Pre & Post).
General Setup Parameters:
Analysis: Initiation
Set the analysis type to Crack Initiation on the main form.
Units: PSI
Set the Stress Units to PSI.
Jobname: notch1a
Give the job a name. Use notch1a.
Title: Notch specimen one with material A
Give the job a title.
Solution Parameter Form: Accept all defaults.
Materials Information Form:
Material: mata
Place the cursor in the cell under the word Material and click on the mouse. A listbox will appear. Select the material MATA from this listbox.
Finish: Polished
Select Polished from the option menu that appears. The word polished appears in the Finish cell. All ten specimens were polished specimen with no surface treatment.
Treatment: No Treatment
Select No Treatment from the option menu that appears.
Region: default_group
Select the default group as the region. It contains all the nodes of the entire model for which the fatigue analysis will be applied to. The Materials Information form can be closed down now by clicking the OK button.
Loading Information Form:
Results From: MSC.Patran FEA
The results are from an external MSC.Patran FEA results file.
Select a MSC.Patran FEA Job: notch01
Press the Select File button and select the notch01.res file from the listbox.
Load Case ID: 1
Place the cursor in the Load Case ID cell and click the mouse button. A databox appears in which the load case ID from the FE analysis is to be entered. Press RETURN to accept the default (1).
Time History: 20KIPR1
Select a time history from the list that appears. Our time history is the one that was previously created and scaled, SAETRN15.
Load Magnitude: 20000
Enter 20,000 lbs as the load magnitude. This is used to normalize the stresses from the FE analysis. Press RETURN to enter this value into the cell. This form can now be closed down by clicking the OK button.
Job Control Form:
Partial Analysis
Set the action to Partial Analysis and click the Apply button. Open the main FATIGUE and Job Control forms again when the translation is complete.
Monitor Job
At this point the job has been submitted and can be monitored as to its progress if desired. Set the Action to Monitor and click Apply each time you wish to see the progress state of the job.
At this point, the job has been submitted and can be monitored as to its progress if desired. Once the message Preprocessing completed successfully appears as the status message, the results can be determined by FEFAT. Note that we did not run the analysis to completion and only went through the preprocessor stage.
Step 5: Evaluate Results
Table 14‑12 and Figure 14‑10 show results for measured vs. MSC.Fatigue predicted lives. Node 1 for all notch specimens has the most damage. To duplicate the results presented in Table 14‑12 for the first specimen using FEFAT, use the following keystrokes.
 
Operation
Comments
Results Form:
Optimize / Node 1
Set the Action to Optimize and click the Apply button in invoke FEFAT. Put “Node 1” in the databox as the node to be analyzed.
Design Life / OK
Press the OK button when FEFAT appears after entering a design life of 1000.
End
Close the summary page down.
Sensitivity Analysis...
Select the sensitivity analysis option.
Scaling Factors
Enter some scaling factors separated by commas from Table 14‑12.
Recalculate
Recalculate
End
Clear the summary page.
Results Display...
Select the results display option.
Sensitivity Plot
Plot the scaled fatigue lives just calculated. Select the Close option from the File menu when done.
Change Parameters...
Select the Change parameters option.
Mean Stress Correction: Morrow
Change the mean stress correction method to Morrow.
Recalculate
Recalculate. You can repeat these operations for no mean stress correction. Compare the values to those in Table 14‑12.
End / eXit
Quit from FEFAT when done.
You may repeat this entire exercise for all other notch configurations if you wish to reproduce all the results in Table 14‑12. If you do, do not forget to use the correct material (A,B,C, or D) and the correct load time history (R=1, or R=0.1).
It might be pointed out that the lives at high load levels do not appear to give good results compared to those at lower loading levels. This is due to the fact that at these higher loading levels the specimens are being literally torn apart where wide areas are experiencing stresses exceeding the yield stress. When this happens the fatigue theory begins to break down since it is based on a local strain approach. There are other correction methods in MSC.Fatigue, namely the Mertens‑Dittmann and the Seeger‑Beste, that can better compensate for high plasticity.
 
Table 14‑12
Specimen No. (mat.)
Load Amplitude (KIP)
Scale Factor
Initiation Life (cycles)
MSC.Fatigue
STW Morrow None
1(A)
R=-1
20
1.0
62
~300
~280
~280
16
0.8
635
~630
~560
~560
14
0.7
1300
~1,000
~880
~880
12
0.6
2400
~1,750
~1,500
~1,500
10
0.5
6000
~3,600
~3100
~3,100
8
0.4
14000
~10,400
~8,700
~8,800
7
0.35
19000
~22,400
~18,350
~18,600
6
0.3
100000
~6.5E4
~5.27E4
~5.38E4
6
0.3
8850
"
"
5.4
0.27
1.8e6
~1.55E5
~1.27E5
~1.31E5
5.2
0.26
500000
~2.19E5
~1.81E5
~1.87E5
2(A)
R=-1
13.99
0.6995
68
~490
~440
~440
12.65
0.6325
190
~690
~600
~600
12.11
0.6055
265
~770
~680
~680
10.65
0.5325
1250
~1,230
~1,080
~1,080
9.02
0.4510
3600
~2,300
~2,000
~2,010
9.03
0.4515
2400
~2,290
~1,990
~2,000
7.00
0.3500
11500
~6,940
~5,880
~5,940
5.68
0.2840
5540
~21,700
~17,940
~18,250
4.95
0.2475
160780
~5.53E4
~44,550
~45,740
4.70
0.2350
188000
~8.26E4
~6.76E4
~6,98E4
3(A)
R=-1
18
0.90
4
~30
~25
~25
14
0.70
47
~50
~50
~50
11.4
0.57
133
~100
~90
~90
9
0.45
1000
~190
~170
~170
6.6
0.33
2000
~500
~450
~450
5
0.25
1000
~1,290
~1,130
~1,140
3
0.15
11150
~11,840
~9,900
~10,020
2
0.10
62180
~2.16E5
~1.81E5
~1.89E5
1.8
0.09
77000
~6.1E5
~5.34E5
~5.65E5
4(B)
R=-1
40
2
17
~35
~30
~30
27.5
1.375
150
~100
~90
~90
18
0.9
742
~440
~370
~370
12.5
0.625
3600
~2,280
~1,920
~1940
9.5
0.475
32500
~13,880
~12,120
~12,480
8
0.4
600000
~6.17E4
~5.82E4
~6.11E4
8
0.4
97000
7.5
0.375
120000
~1.15E5
~1.12E5
~1.19E5
7
0.35
190000
~2.31E5
~2.31E5
~2.46E5
5(B)
R=-1
4.2
0.3
15000
~530
~450
~450
3.5
0.25
35000
~1,160
~980
~990
2.8
0.2
92000
~3,800
~3,240
~3,300
2.3
0.164
250000
~15,420
~13,620
~14,140
6(C)
R=-1
17.64
0.441
3603
~2,380
~2,230
~2,240
17.64
0.441
2800
14.70
0.3675
13340
~5,073
~4,760
~4,790
11.76
0.294
51960
~13,540
~12,750
~12,850
10.59
0.265
129000
~21,940
~20,700
~20,890
9.41
0.2353
230460
~39,030
~36,930
~37,330
8.82
0.2205
424750
~5.41E4
~5.13E4
~5.19E4
8.24
0.206
589510
~7.7E4
~7.32E4
~7.41E4
7.75
0.1938
1180700
~1.07E5
~1.02E5
~1.03E5
7(C)
R=-1
20.25
0.5
924
~490
~450
~450
18.00
0.45
2400
~730
~670
~680
15.75
0.3938
5400
~1,220
~1,130
~1,140
12.80
0.32
14000
~2,800
~2,610
~2,620
11.25
0.2813
25000
~4,800
~4,480
~4,500
9.00
0.225
87315
~12,800
~11,960
~12,040
7.43
0.1858
170000
~31,500
~29,500
~29,760
8(D)
R=0.1
4.8
0.12
4850
~5,560
14,110
~41,040
4.0
0.1
8850
~11,930
~41,330
~1.37E5
3.2
0.08
22600
~41,640
~2.04E5
~6.02E5
2.4
0.06
27500
~2.79E5
~1.83E6
~4.09E6
1.6
0.04
2650000
~4.16E6
~3.61E7
~6.1E7
9(D)
R=0.1
4.8
0.12
2190
~1,950
~3,240
~6,570
4.0
0.1
5900
~3,790
~8,120
~21,340
3.2
0.08
10800
~9,100
~28,430
~9.26E4
2.4
0.06
27500
~41,700
~2.08E5
~6.25E5
1.6
0.04
2100000
~6.13E5
~4.53E6
~9.31E6
10(D)
R=0.1
4.8
0.12
945
~620
~650
~840
4.0
0.1
1940
~1,320
~1,900
~3,180
3.2
0.08
4380
~2,960
~5,670
~13,290
2.4
0.06
10000
~9,280
~28,260
~8.77E4
1.6
0.04
70500
~1.02E5
~5.35E5
~1.3E6
Figure 14‑10 Measured and Predicted Lives of Notched Specimens