OS-T: 6060 Random Response Fatigue Analysis

The study of fatigue life of a structure subjected to non-deterministic loads.

Power Spectral Density (PSD) from the Random Response Analysis are used to calculate moments that are used to generate the probability density function for the number of cycles versus the stress range.

The following file found in the optistruct.zip file is needed to perform this tutorial. Refer to Access the Model Files.

CAD-box-SN.fem

or

A copy of the model file used in this tutorial is available on <install_directory>/tutorials/hwsolvers/optistruct.

A CAD box which is subjected to random response analysis is utilized to perform the Random Response Fatigue Analysis. The random response analysis setup is already made for this model, an additional loadstep for SN-Fatigue calculation is created in this tutorial. The loading frequency is defined by the FREQ1 card and RANDPS is used to apply the Auto and the Cross PSDs.


Figure 1. CAD Box Model for Fatigue Analysis

Launch HyperMesh and Set the OptiStruct User Profile

The model being used for this exercise is that of a simple CAD box (Figure 1). Four loadsteps have already been defined on this model, each of which represent Normal Modes Analysis, Frequency Response at both the constraint points and the Random Response Analysis.

  1. Launch HyperMesh.
    The User Profile dialog opens.
  2. 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

  1. Click File > Import > Solver Deck.
    An Import tab is added to your tab menu.
  2. For the File type, select OptiStruct.
  3. Select the Files icon files_panel.
    A Select OptiStruct file browser opens.
  4. Select the CAD-box-SN.fem file you saved to your working directory. Refer to Access the Model Files.
  5. Click Open.
  6. Click Import, then click Close to close the Import tab.
    The outline of the Fatigue Analysis setup to be achieved in the following steps.


    Figure 2. Fatigue Setup - Random SN Damage

Set Up the Model

Define FATLOAD Load Collector

The model has a Random Response loadstep defined by the RANDPS card. The RANDPS card contains the two input PSD (AutoPSD_A and AutoPSD_B) applied at the constrained locations and the cross PSD (CrossPSD_AB) between them.

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter FATLOAD_RANDOM.
  3. For Card Image, select FATLOAD.
  4. For LCID(load case ID), select Random Response from the list of load steps.
    Note: TABFAT and scaling parameters are not required for this calculation.

Define FATEVNT Load Collector

Create a random response event for the FATLOAD_RAND created.

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter EVENT_RANDOM.
  3. For Card Image, select FATEVNT.
  4. For FATEVNT_NUM_FLOAD, enter 1.
  5. Click on the Table icon table_pencil next to the Data field and select FATLOAD_RANDOM for FLOAD in the pop-out window.

Define FATSEQ Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter FATSEQ.
  3. For Card Image, select FATSEQ.
  4. For FATSEQ_NUM enter 1, as 1 FATEVENT has been created.
  5. For FID (Fatigue Event Definition), select EVENT_RANDOM and N as 1.


    Figure 3. FATSEQ showing RAND_EVENT created
    Defining the sequence of events for the fatigue analysis is completed. The Fatigue parameters are defined next.

Define Fatigue Parameters

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter Fat-Parameter.
  3. For Card Image, select FATPARM.
  4. Verify TYPE is set to SN.
  5. Set COMBINE to ABSMAXPR.
  6. Set UCORRECT to NONE.
  7. Set STRESSU to MPA (Stress Units).
  8. Set CERTNTY SURVCERT to 0.9.
  9. Set RNDPDF as:

    DM1 = DIRLINK

  10. Check the Random option under RANDOM and edit the following options:

    FACSREND = 8.0 (Calculates the upper limit of the stress range (SREND))

    NBIN = 100 (Calculates the width of the range ( MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqGHciITaaa@3A3D@ S) of stress ranges for which the probability is calculated)



    Figure 4. Probability Density Function. (probability density of number of cycles versus stress range)

Define Fatigue Material Properties

The material curve for the fatigue analysis can be defined on the MAT1 card.

  1. In the Model Browser, click on the MAT1 material.
    The Entity Editor opens.
  2. In the Entity Editor, set MATFAT to SN.
  3. Set UTS (ultimate tensile stress) to 80.0.
  4. For the SN curve set (these values should be obtained from the material's SN curve):
    SRI1
    317.0
    B1
    -0.15
    NC1
    1e6
    SE
    0.279192629
    A/R
    A

Define PFAT Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter Surface_Treatment.
  3. For Card Image, select PFAT.
  4. Set LAYER to WORST.
  5. Set FINISH to NONE.
  6. Set TRTMENT to NONE.
  7. Set Kf to 1.0.

Define FATDEF Load Collector

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter Fat-definitions.
  3. Set the Card Image to FATDEF.
  4. Activate PTYPE and PSOLID in the PTYPE Entity Editor.
  5. Edit FATDEF_PSOLID_NUMIDS to 1.
  6. Click on the Table icon table_pencil next to the Data field and select Box for PID, and Surface Treatment for PFATID in the pop-out window.
  7. Click Close.

Define the Fatigue Load Step

  1. In the Model Browser, right-click and select Create > Load Step.
  2. For Name, enter Vibration_Fatigue.
  3. Set the Analysis type to fatigue.
  4. For FATDEF, select fatdef.
  5. For FATPARM, select fatparam.
  6. For FATSEQ, select fatseq.

Submit the Job

  1. From the Analysis page, enter the OptiStruct panel.
  2. Click save as following the input file field.
    The Save As dialog opens.
  3. For File name, enter the name CAD_Box_Vibration_Fatigue.fem.
  4. Click Save.
  5. Click OptiStruct to submit the analysis.

Review the Results

  1. 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.
  2. Go to the Results tab.
  3. In the Results tab, select Subcase 5 (Vibration_Fatigue) from the subcase field.
  4. On the Results toolbar, click resultsContour-16 to open the Contour panel.
  5. Set Result type to Damage and click on Apply to contour the elements.
    Figure 5. Damage Contour Plot