F.D.& E. Experimental Fatigue life Prediction Program: saefcalc2

Note! This calculator uses "Digital Curves" as material property inputs.

How to run simulation:

• pdrain.f Program to perform a Rainflow cycle count on your peak/valley load history
  (at github: Click on "Raw" button to download)
  
  
• saefcalc2.f Program to compute crack initiation life
  
  
• makeRepSaefcalc2 Script to process output from saefcalc2.f
  (needs gnuplot and htmldoc on your Linux system)
  
  
• delete1arg.f (used by makeRepSaefcalc2 to create report)
  
  
• merged_a36_fitted.html Sample "fitted" fatigue curve for A36 Steel.
  
  
• loads4rain.out Sample loads file for input to pdrain.f

You can perform a crack initiation prediction with the commands (example):

Edit the bash script makeRepSaefcalc2 and add user name and affiliation etc.

If you have not already done so, you will need to make  the script "makeRepSaefcalc2" 
   executable with command:
      chmod 744  makeRepSaefcalc2


Rainflow Cycle Count your load history:
      ./pdrain  2 < loads4rain.out   > sample.rain 

Run the simulation:
      ./saefcalc2   merged_a36_fitted.html  5.4 < sample.rain   > a36sample_5.4

Create the report:
      ./makeRepSaefcalc2  a36sample_5.4 

View the *.pdf report:
      evince  a36sample_5.4.pdf 

Example Reports :

• a36sample_5.4.pdf Example output report in pdf format
• a36sample_5.4.html Example output report in html format

Local Files for download (mostly same as on github.com):

• Calculation Programs:
  saefcalc2.f.txt Vers 2.2, Feb2016 Computes life using standard rainflow file input.
  pdrain.f.txt Rainflow cycle counting program. vers.1.2 Feb 2016
  
  
• Material fitted curve fatigue file examples:
  merged_a36_fitted.html.txt A36 steel merged file
  AA7075_fitted.txt sample fitted digital fatigue curve
  
  
• History and Rainflow count file examples:
  loads4rain.out.txt short history file
  sample.rain.txt Rainflow count file for above file loads4rain.out.
  
  grpltorweld.rain.txt grapple log skidder sample rainflow file
  
  
• Sample Output file from saefcalc2 run:
  a36sample_5.4.txt
  
  
• Auxillary programs to create report:
  delte1arg.f.txt Prog. to delete 1st string in each line of input
  "htmldoc download" Download if not available on your Linux sys.
  "gnuplot download" Download if not available on your Linux sys.
  
  
• Script to create *.pdf and *.html reports:
  makeRepSaefcalc2.txt Script to create the *.pdf and *.html reports
  
  

After download a fortran program can be compiled with the statements:   

      mv saefcalc2.f.txt saefcalc2.f

In newer Linux versions the compile command to use is probably:

      gfortran -g -w -fbounds-check saefcalc2.f  -o saefcalc2

Also note that the new compiler version do not preserve variables between calls to 
subroutines, and that one must enter a "      SAVE" line at the top of each subroutine. 
This changed "feature" can cause a lot of trouble when compiling older programs. 



For details of output file data extraction see (and change if you wish) the
script makeRepSaefcalc2


The output from saefcalc2 goes into a large results file that contains 
all sorts of intermediate information plus the final results.  The final 
result lines each have a "#name " at the begining of each line, and can 
be extracted from the general file using grep, as below.

(Assuming "temp" is the saefcalc2  output file name)
View all output:
       more temp

View life results:
       grep xcalc2 temp
       grep xcalc3 temp

View stress-strain loops for plotting:

       grep plotloops temp | delete1arg >temploop

  Then use 
           gnuplot

  In gnuplot type the commands:
           set grid
           plot "temploop" with lines 

  to get the hysteresis loop shapes.

View damage per input block:
       grep plotdam temp