Fatigue User’s Guide > Crack Growth > Crack Growth Prediction (PCRACK)
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">   
Crack Growth Prediction (PCRACK)
PCRACK is a cycle-by-cycle fatigue crack growth prediction module which calculates the crack extension, da, in each cycle of a sequence of cycles and adds it to the current crack size, a. This process proceeds until the failure condition is reached or a preset crack size is achieved. The driving force for crack growth in each cycle is the effective range of stress intensity, delta K (). In each cycle, the apparent or applied is calculated from the stress range, the crack size, and the geometry of the component or structure under consideration. This is then modified to effective to account for the possible occurrence of:
crack closure
history effects
notch effects
environmental influence
static fracture mode contribution
Consequently, the crack propagation life is modeled accurately and this can only be achieved on a cycle-by-cycle basis.
There are many applications of this approach including: design analysis, pre‑prediction of test programs, defect tolerance assessment, inspection related strategy, failure investigation and decision support. These are discussed in detail in Fatigue Theory (Ch. 15).
Module Operation
The operation of PCRACK can be in two modes: within the MSC Fatigue Pre & Post and MSC Patran environments or in stand alone mode. The only difference is that in stand alone mode, the user must supply the jobname. (In direct mode from a preprocessor, these are passed to PCRACK automatically.) PCRACK can be accessed directly from the operating system prompt by typing the symbol pcrack. Once PCRACK has been initiated in either of these modes, two windows will be presented. The top, small form is a generic form and allows for general program control. This is discussed in detail in Module Operations (Ch. 17) for the Motif driver.
Figure 7‑28 PCRACK Utility Form
Six items are available from the main menu of PCRACK which are discussed in the table below:
 
Option
Description 
Preprocessing
Perform rainflow cycle counting and then re-orders the cycles in the order in which the component experiences them. The operation of this option is identical to that described in Fatigue Preprocessing, 277.
Analysis
Performs the actual crack growth analysis. This section describes this in detail.
Output a Time History
The combined response time history can be output to a file. See Output Time Histories, 315.
Display a Time History
The output response can be plotted. See Graphical Display of Time Histories, 318.
Utilities
These utilities are identical to those described in Utilities, 324.
eXit
Exits the program.
Results of a crack growth analysis can be viewed using a utility called PCPOST. Detailed operations of PCPOST are in Reviewing Results (PCPOST), 547.
Once an analysis is initiated in interactive mode, five input screens containing various questions appear sequentially in the order shown below.
General Information
Figure 7‑29 shows the first screen that appears. If the user accesses PCRACK from the system level, he will have to supply a jobname. PCRACK already knows the jobname when entered from MSC Fatigue Pre & Post or MSC Patran.
Figure 7‑29 The PCRACK General Info Page
 
The fields to this screen are as follows:
 
Field
Description
Input File Name
The default will be the last jobname used in conjunction with a MSC Fatigue crack growth analysis. The jobname will pick up an existing jobname.tcy file which contains the local stress information. The jobname.tcy file is produced by the program FEFAT, the crack growth fatigue preprocessor. See Fatigue Preprocessing (Ch. 5).
Scaling Factor
This factor allows for scaling the stresses linearly to simulate the effect of a higher level of loading. The default answer is 1 (i.e., leave the stresses unscaled).
Time History Offset
The mean offset of the stresses may be altered here to simulate the effect of a further mean load or residual stresses. Only positive offsets are allowed. The default answer is 0 (i.e., no mean offset will be added). Remember that crack growth is highly sensitive to mean stress since the stress is raised to a power of m in the calculation of crack increment where m is the Paris law slope.
Plane Stress Correction
The materials data describing the crack growth rate relationship are obtained from test coupons which are stresses in a plane strain condition. A correction may be applied to modify these data to simulate a plane stress condition. The default answer to this question is N, or NO (i.e., no plane stress correction will be applied).
Model History Effects
If history effects are to be modeled, answer Yes to this question. Modelling the history effect means that the growth in one cycle is affected by the previous cycle history leading to retardations, accelerations, overshoots and other effects on growth rate and hence propagation life. It may be desirable to “turn off” the history effect to compare life predictions based on linear damage accumulation (no history effect) and those with history effects modeled.
At this point, the screen will clear by clicking OK, and the next page of input questions (Output Parameters) will be presented.
Output Parameters
The output parameters show up on the next screen shown in Figure 7‑30.
Figure 7‑30 The PCRACK Output Parameters Page
The fields to this form are as follows:
 
Field
Description
Create Result FIle
If a results file is desired answer YES.
Results File Name
The PCRACK analysis generates a results file which will have the default name jobname.crg. An alternative name may be supplied in response to this question. If an output file already exists with the same name as the one supplied, overwrite permission will be requested.
Results File Output Interval
A snapshot of the calculation parameters will be taken at regular intervals and stored in the output results file (jobname.crg). The user may select how often the snapshot is taken by defining the number of life units between snapshots. The default will cause the gap to be between 100 and 1000 cycles, converted into user life units. If the snapshots are too frequent, the crack growth results postprocessor (PCPOST) allows the user to decimate or compact the results file to save disk space.
Type of Display
The graphical status display shows the rate of crack advance using an a-N graph and various computational parameters. If the user requests not to have a graphical display, a tabular display will be produced indicating crack length, elapsed life in repeats of the loading history and elapsed life in user defined units.
Screen Update Interval
The PCRACK graphical screen display will be updated at regular intervals. The update frequency should be chosen to reflect the speed of the computer on which the analysis is being carried out. An update interval of between 100 and 1000 cycles will be chosen using the default. Faster computers will require longer update intervals to prevent the computations becoming bound by the screen update rate.
Target Design Life
The crack growth analysis may be terminated after a set number of user life units as specified here. If the number of life units is left undefined as for the default answer, the analysis will proceed until another failure criterion is satisfied.
At this point, the screen will clear by clicking OK, and the next page of input questions (Local Geometry) will be presented.
Local Geometry
The next screen displayed is shown in Figure 7‑31. The user will first be prompted for a file name containing the compliance function for the crack geometry of interest. This file was generated by PKSOL (see the previous section). Once the file name has been selected by clicking the OK button, the screen updates and shows the following fields which can be modified.
Figure 7‑31 The Local Geometry Definition Screen
The fields of this form are described below:
 
Field
Description
K Solution File Name
The K solution is generated by the program PKSOL and is stored in a file with file type filename.ksn. The K solution relates crack size and compliance. For more detail, refer to K Solution Library (PKSOL), 508. The default K solution file name offered will be the last K solution used or generated.
The local geometry will now be defined. Note that the local geometry fields do not appear until the K solution file name is supplied and accepted. The length units will be reported to you as already defined in the K solution file. For example:
Units of length: inches
The user will then be asked to define the physical size of any additional notch not modeled using FEA and not accounted for in the compliance function.
Notch Depth
If the crack under consideration is a surface crack, then a surface notch can be introduced into the global geometry from which the crack can grow. The crack will be influenced by the presence of the notch and its associated stress field. To correct the K solution for this, the local geometry of the notch has to be defined. The notch depth must be defined in the same length units as used for the crack definition when generating the K solution. Valid units will be either mm, m, milliinches or inches. If a notch is not present, (i.e., notch depth is zero), then the next two questions will not be asked and their defaults will be used.
Notch Root Radius
This input is the notch root radius which defaults to 0.000039in if the Notch Depth above is zero. It must be defined in the same length units as used for the crack definition when generating the K solution.
Sharp Crack Root Radius
This is the user-specified value of the root radius of a “sharp notch” (i.e., a real fatigue crack, which defaults to 0.000039in if the Notch Depth above is zero). The sharp crack root radius must be defined in the same length units as used for the crack definition when generating the K solution. Valid units are as for the previous question.
Maximum Crack Size
The default maximum crack size (n) is the component dimension in the cracking direction. This field enables you to specify a smaller value to take into account plastic collapse, yielding, etc. Note that if fracture toughness is exceeded at a crack size less than the minimum, then this condition will define the end of fatigue crack growth. All crack lengths MUST be defined in the same units as used for crack geometry when creating the K solution and you will be reminded of these units at the top of the input screen.
Initial Crack Size
This is a critical value. It should be either:
1. The size of a discovered defect.
2. A postulated defect.
3. The NDT limit (i.e., the largest defect that cannot reliably be detected by nondestructive testing).
The initial crack size may be chosen in the range 0 to n where n is the maximum possible physical length of the crack based on the component geometry. A zero initial crack size may only be specified if a notch has been defined in the local geometry definition (see above). If no notch has been defined and you enter zero for the initial crack size, then PCRACK will compute the short crack parameter, lo, and substitute this for the initial crack size. See Cycle-by-Cycle Approach, 579, for details on lo.
Material and Environment
The next form enables the user to specify the material type and the test environment (Figure 7‑32).
Figure 7‑32 The Material and Environment Selection Screen
 
Field
Description
Data Source
This option allows users to select materials data from a user-defined database. The default is to search the standard database, which will exist in either the local or central directory.
Database Name
The name of the user-defined database to use. This database may exist anywhere on the system, but the default is to look in the local directory.
Material Name
The material name is the name used in the materials database under which are stored the parameters defining the crack growth properties of the material. The default name offered will be the last materials data set used. If the user wishes to specify another material, he may either type in a name that he knows exists or may use the List capability to select a material from a list.
The screen that appears when a materials listing is shown allows the user to view not only all the available materials but also to view the material parameters associated with each material. This eliminates the need to run PFMAT to view these parameters during an interactive run.
Environment
There may be crack growth materials data sets stored in the database for a range of different environments. The user may chose which environment to use in the analysis. If an environment is specified for which there are no data in the database, the parameters for AIR will be assumed and an environment definition will be requested (see Environment Definition below). If the user wishes to specify another environment, he/she may either type in a name that he knows exists or may use the List capability to select an environment from a list.
Edit Data
This option allows the user to edit the materials data for the selected environment.
Material Parameter Editing
At this point, the form will clear and the materials parameters retrieved from the database presented either for air or for the chosen environment (Figure 7‑33). The user may change any of these parameters for the duration of the PCRACK session. See Rules for Changing Young’s Modulus, 133.
Figure 7‑33 The Material and Environment Edit Screen
The analysis will now commence. If the status has been requested to appear graphically, then the screen will clear, and a graphics screen display will be updated regularly. The parameters on the top of this screen display are: elapsed life, accumulated crack size, current apparent , current effective and current growth rate da/dN.
Also, highlighted on the right-hand side, are the complexing effects currently modifying the apparent to effective (i.e., closure, history, notch, environment and static fracture modes).
Within the central rectangular area is a graphical display which updates and rescales automatically at regular intervals as specified. The graph shows crack size plotted against life, as shown in Figure 7‑34.
Figure 7‑34 On-Line Graphics Display of Crack Size — Cycles
If the crack growth is negligibly small, then the process will cease after one repeat of the history or the design life, whichever is greater; this is to avoid infinite or exceptionally long runs. The crack growth analysis may be aborted at any time by clicking on the Abort option on the main form. The other options, Alpha and Replot, will convert the graphical display to an alphanumeric display and replot the screen respectively. If an alphanumeric display is requested it may be converted back to graphics by selecting the Pause button. A window will appear giving the options of Resuming, Quitting, or returning to Graphical display.
When the process stops, whether it is due to the final crack size or fracture toughness being exceeded or due to negligible growth or to a user abort, the final situation will be presented to you as shown in Figure 7‑35.
Figure 7‑35 Final Situation Display
The start time and end time can be used to calculate the run time and processing speed in user-defined events per second. The list of parameters gives the current values at the point of termination. The final situation results page is also written to the pfatigue.prt file for later reference.
From this screen, the user may then proceed to the post analysis menu which follows on the next screen display.
Post Analysis Menu
After an initial crack growth job has been run, the Post Analysis menu appears as in Figure 7‑36. The options are described below.
Figure 7‑36 The Post Analysis Menu
 
Field
Description
Final Situation
This option gives a screen of information as shown in Figure 7‑35. The results are saved in a jobname.crg type file.
Final a-N Graph
This gives the autoscaled final graphics plot ( Figure 7‑34 ).
Fatigue Life
Interpolation
This option gives the user the opportunity to interpolate the crack size and elapsed life data for different values of either the initial crack size or the final crack size or both. The defaults offered in the following example (Figure 7‑37) are the original initial and final crack sizes. It is possible to interpolate life beyond the crack size limits of the initial calculation. However, such interpolation is only approximate as the crack size-life relationship is both nonlinear and bounded by certain physical constraints.
Note: The initial crack size used when the user enters a value of zero is in fact the short crack parameter, lo. The value of this parameter for the current job is the default value shown in the initial crack size field. See Cycle-by-Cycle Approach, 579.
Start Another Job
As the title suggests, this loops the user back to screen 1 of this program and allows him to go through the whole process again.
Edit Analysis Parameters
This will bring up a submenu ( Figure 7‑38 ) which will allow for editing parameters associated with the initial job setup each of which has been discussed in detail previously.
Preferences
This option allows for turning material checking on or off, setting a failure criterion, or doing an auto-overwrite.
The material parameters generally have certain valid ranges which the program checks for. If the user enters a value out of the valid range he will be told and requested to enter a valid number. He may see what the valid range is for a given parameter by using the Help on Question option from the help pull-down menu on the top bar when the appropriate field is active. If the user does not want to be bothered with this, he must set the material checking preference off ( Figure 7‑39 ).
The failure criterion option determines whether residual k is included in the calculation of k to exceed klc.
The auto-overwrite option sets a flag for this job. If the answer is Yes, then all output files are automatically deleted without conformation. If No, the program will request confirmation prior to deleting any output files.
Recalculate
After editing parameters, changing materials, etc., the user will need to select this option to see the effect of any of these changes. This reruns the analysis.
Exit
Allows for exiting the program and puts the user back either into the MSC Fatigue forms or back to the system prompt from which PCRACK was invoked.
Figure 7‑37 The Interpolate Fatigue Life Screen
Figure 7‑38 The Edit Parameters Submenu
Figure 7‑39 The Material Checking Preference