OS-T: 6020 Multiaxial Fatigue Analysis (Strain - Life) using E-N Approach
The E-N (Strain - Life) approach should be chosen to predict the fatigue life when plastic strain occurs under the given cyclic loading.
![rd2070_SN_curve](../../../images/solvers/rd2070_sn_curve.png)
![rd2080a_strainlife_curve](../../../images/solvers/rd2080a_strainlife_curve.png)
In OptiStruct various strain combination types are available with the default being "Absolute maximum principle strain". In general "Absolute maximum principle stain" is recommended for brittle materials, while "Signed von Mises strain" is recommended for ductile material. The sign on the signed parameters is taken from the sign of the Maximum Absolute Principal value.
![](../../../images/solvers/multiaxial_en_fatigue_setup_flowchart.png)
- FATDEF
- Defines the elements and associated fatigue properties that will be used for the fatigue analysis.
- PFAT
- Defines the finish, treatment, layer and the fatigue strength reduction factors for the elements.
- MATFAT
- Defines the material properties for the fatigue analysis. These properties should be obtained from the material's E-N curve (Figure 2). The E-N curve, typically, is obtained from completely reversed bending on mirror polished specimen.
- Fatigue Parameters
Figure 4. Mean Stress Correction
- FATPARM
- Defines the parameters for the fatigue analysis. These include stress combination method, mean stress correction method (Figure 4), Rainflow parameters, and Stress Units.
- Fatigue Sequence and Event Definition
Figure 5. Load Time History
- FATSEQ
- Defines the loading sequence for the fatigue analysis. This card can refer to another FATSEQ card or a FATEVNT card.
- FATEVNT
- Defines loading events for the fatigue analysis.
- FATLOAD
- Defines fatigue loading parameters.
- TABLEFAT
- Defines the y values for each point on the time loading history (Figure 5).
Copy the following model files to your working directory. Refer to Access the Model Files.
ctrlarm_EN.fem, load1.csv and load2.csv
or
A copy of the model files used in this tutorial are available on <install_directory>/tutorials/hwsolvers/optistruct.
![rd2070a_control_arm](../../../images/solvers/rd2070a_control_arm.png)
![rd2070a_load_time](../../../images/solvers/rd2070a_load_time.png)
![rd2070a_vertical_force](../../../images/solvers/rd2070a_vertical_force.png)
![rd2080a_en_curve](../../../images/solvers/rd2080a_en_curve.png)
Launch HyperMesh and Set the OptiStruct User Profile
The model being used for this exercise is that of a control arm (Figure 6). Loads and boundary conditions and two static loadcases have already been defined on this model.
-
Launch HyperMesh.
The User Profile dialog opens.
-
Select OptiStruct and click
OK.
This loads the user profile. It includes the appropriate template, macro menu, and import reader, paring down the functionality of HyperMesh to what is relevant for generating models for OptiStruct.
Import the Model
-
Click
.An Import tab is added to your tab menu.
- For the File type, select OptiStruct.
-
Select the Files icon
.
A Select OptiStruct file browser opens. - Select the ctrlarm_EN.fem file you saved to your working directory.
- Click Open.
- Click Import, then click Close to close the Import tab.
Set Up the Model
Define TABFAT Load Collector
- Make sure the Utility menu is selected in the View menu. Click .
- Click on the Utility menu beside the Model tab in the browser. In the Tools section, click on TABLE Create.
- Set Options to Import table.
- Set Tables to TABFAT.
- Click Next.
- Browse for the loading file.
- In the Open the XY Data File dialog box, set the Files of type filter to CSV (*.csv).
- Open the load1.csv file you saved to your working directory.
- Create New Table with Name table1.
-
Click Apply to save the table.
The curve table1 with TABFAT card image is created.
- Browse for a second loading file load2.csv.
- Create New Table with Name table2.
-
Click Apply to save the table.
The curve table2 with TABFAT card image is created.
-
Exit from the Import TABFAT window.
Tables appear under Curve in the Model Browser.Note: A file in DAC format can very easily be imported in HyperGraph and converted to CSV format to be read in HyperMesh.
Define FATLOAD Load Collector
- In the Model Browser, right-click and select .
- For Name, enter FATLOAD1.
- Click Color and select a color from the color palette.
- For Card Image, select FATLOAD.
- For TID(table ID), select table1 from the list of curves.
- For LCID (load case ID), select SUBCASE1 from the list of load steps.
- Set LDM (load magnitude) to 1.
- Set Scale to 5.0.
- Repeat the process to create another load collector named FATLOAD2 with FATLOAD Card Image and pointing to table2 and SUBCASE2.
- Set LDM to 1 and Scale to 5.0.
Define FATEVNT Load Collector
- In the Model Browser, right-click and select .
- For Name, enter FATEVENT.
- For Card Image, select FATEVNT.
- For FATEVNT_NUM_FLOAD, enter 2.
-
Click on the Table icon
next to the Data field and select FATLOAD1 for FLOAD(1) and FATLOAD2 for FLOAD(2) in the pop-out window.
Define FATSEQ Load Collector
- In the Model Browser, right-click and select .
- For Name, enter FATSEQ.
- For Card Image, select FATSEQ.
-
For FID (Fatigue Event Definition), select FATEVENT
.
Defining the sequence of events for the fatigue analysis is completed. The Fatigue parameters are defined next.
Define Fatigue Parameters
- In the Model Browser, right-click and select .
- For Name, enter fatparam.
- For Card Image, select FATPARM.
- Verify TYPE is set to EN.
- Set MAXLFAT to Yes for the multiaxial method.
- Set STRESSU to MPA (Stress Units).
- Set RAINFLOW RTYPE to LOAD.
- Set CERTNTY SURVCERT to 0.5.
Define Fatigue Material Properties
The material curve for the fatigue analysis can be defined on the MAT1 card.
-
In the Model Browser, click on the Aluminum material.
The Entity Editor opens.
- In the Entity Editor, set MATFAT to EN.
- Set UTS (ultimate tensile stress) to 600.
-
For the EN curve set (these values should be
obtained from the material's EN curve):
- SF
- 1002.000
- B
- -0.095
- C
- -0.690
- EF
- 0.350
- NP
- 0.110
- KP
- 966.000
- NC
- 2E+08
- SEE
- 0.100
- SEP
- 0.100
Define PFAT Load Collector
- In the Model Browser, right-click and select .
- For Name, enter pfat.
- For Card Image, select PFAT.
- Set LAYER to TOP.
- Set FINISH to NONE.
- Set TRTMENT to NONE.
- Set Kf to 1.0.
Define FATDEF Load Collector
- In the Model Browser, right-click and select .
- For Name, enter fatdef.
- Set the Card Image to FATDEF.
- Select the PTYPE check box and activate PSOLID.
- Click <Unspecified> in PID: and select Property collector and then select shell.
- Click <Unspecified> in PFATID: and select Loadcol and then select pfat.
- Click Close.
Define the Fatigue Load Step
- In the Model Browser, right-click and select .
- For Name, enter Fatigue.
- Set the Analysis type to fatigue.
- For FATDEF, select fatdef.
- For FATPARM, select fatparam.
- For FATSEQ, select fatseq.
Submit the Job
- From the Analysis page, enter the OptiStruct panel.
-
Click save as following the input file field.
The Save As dialog opens.
- For File name, enter the name ctrlarm_EN_fatigue.fem.
- Click Save.
- Click OptiStruct to submit the analysis.
Review the Results
-
From the OptiStruct panel, click HyperView.
HyperView is launched and the results are loaded. A message window appears to inform of the successful model and result files loading into HyperView.
- Go to the Results tab.
- Change the Load Case to Subcase 3 - fatigue.
- From the Results Browser (left window pane), toggle Components PSOLID_2 and PSOLID_5 off.
-
On the Results toolbar, click
to open the Contour panel.
- Set Result type to Life and click on Apply to contour the elements.
- Right-click on top of the Contour Plot numeric values (in the graphics window) and select Edit Legend.
-
Select Interpolation: Log and click
OK.
Figure 10. Elemental Life results indicating ~4500 cycles before the first element fails