Fatigue User’s Guide > Crack Growth > K Solution Library (PKSOL)
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX''">   
K Solution Library (PKSOL)
The K Solution library module, PKSOL, contains the means to calculate the fracture mechanics stress intensity factor, K, for over 35 geometric cases. The module provides a means to generate a compliance function table for use in the MSC Fatigue fracture mechanics prediction module PCRACK where K is calculated as a function of crack size, stress history, and geometry.
Background Information
The fracture mechanics stress intensity factor, universally known as K (or in fatigue), relates stress and crack size into a single parameter which is the driving force for the behavior of cracks in engineering materials in terms of fracture, stress corrosion or fatigue. Its use relies on the similitude concept which implies that, independent of geometry, if K is the same then the response is the same. The relationship may be viewed as a triangle in Figure 7‑1 in which if two are known the third can be derived. In this way, important decisions can be made regarding, for example, specified material toughness, NDT limits and inspection strategy, and operating stresses in engineering applications.
Figure 7‑1 The Fracture Mechanics Triangle
The three parameters are related by the general expression:
(7‑1)
where Y, the compliance function embodies both the global geometry of the component or structure and the crack size, shape, and location. Consequently, many cases, or K solutions, are available in standard references ((Ref. 3.) through (Ref. 8.)) and more than 35 of the most significant solutions for engineering applications are coded into the PKSOL library module. They may be used to generate a compliance function table for a specific geometry and case, for all crack sizes.
Example
As a specific example of using PKSOL to solve for one corner of the fracture mechanics triangle and its application to a crack growth analysis:
Suppose you want to be able to stop the analysis at a crack length that would cause a failure if the limit load occurred in the next cycle. Imagine you have a time history of random cycles varying between 1 and 3kips that represents a repeatable sequence, say one flight, but does NOT include the once per life 10kips limit load. The analysis would proceed to model the crack growth cycle by cycle including all the retardation, overshoot, and crack closure effects, but there are several ways of ending the process. One of them is to specify the final crack length.
The analyst would use the K solution library module (PKSOL), as a single shot solver, to calculate the crack size when the fracture toughness is exceeded at the maximum stress corresponding to the limit load. (It is also possible to calculate a final crack size when plastic collapse of the remaining section occurs rather than a KIC fracture). This crack size is then used to limit the crack growth process and effectively checks on a cycle by cycle basis "have I reached a crack size where the residual strength compromises the design limit strength?"
Module Operation
PKSOL may be accessed from the MSC Fatigue Solution Parameter form available from the main menu within MSC Fatigue Pre & Post or MSC Patran when the Analysis Type is set to Crack Growth, or it may be run outside of MSC Patran by typing the symbol pksol. A general description of the steps involved in a PKSOL session is shown in Table 7‑1 to calculate a missing corner of the triangle.
 
Table 7‑1
Step
Description
Examples
1
Select units
mm
2
Select calculation type selection
Solve for stress intensity factor, K
3
Select geometry type selection (which general shape)
Cracks at holes
4
Select specific geometry shape selection
Singe crack in hole in tension
5
Input data to define specimen and crack geometry
User enters specimen dimensions (via define and crack size, stress range)
6
Select post-processing options
Plot K against crack sizes for the above case(s)
Once initiated, the PKSOL program will ask the user to select the crack length units to be used as in Figure 7‑3.
When first invoked PKSOL appears displaying two forms.
Figure 7‑2 PKSOL Utility Form
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‑3 PKSOL Units Selections
Having defined the length units, the user will then be asked to define the type of calculation that he wishes to carry out. See Figure 7‑4.
Figure 7‑4 PKSOL Calculation Selections
The user must nominate one unknown on the fracture mechanics triangle (options 1-3), and over the subsequent few screens, supply the known data and define the specimens geometry.
 
Option
Description
Options 1-3
The first three options solve a different corner of the triangle BUT their operations are essentially the same. The user will be asked for a file name to store the output, however this file name is purely optional. Press the <RETURN> key if it is not necessary to store the information in a file. If a file name is supplied, it will be called filename.lst and will contain the results of the analysis in tabular format for printing.
Option 4
This generates a lookup table of the compliance function, Y, and is used in (7‑1). The data in the Y function table is used in the crack growth calculation module PCRACK. If the user picks this option he will be asked for a file name which is mandatory. The result file will be called filename.ksn. This file if of binary format similar to the DAC files generated for load-time histories.
The K solutions are grouped into nine general areas on the geometric definition menu from which you should make a selection (Figure 7‑5).
Options 5 and 6
Selecting one of these options will display the K Solution Library. Selecting one will graphically display the geometries or list references.
Figure 7‑5 The Ten Basic Geometry Groups
For each group, another menu of specific geometric cases is given and selection from this leads to an input screen incorporating a graphic representation of the geometry highlighting unambiguously the required input dimensions. The various menus of specific cases and their geometries are shown in Figure 7‑6 through Figure 7‑13.
When a graphical image is presented for dimension input, the user must first use the mouse to select Define from the menu bar. He may then input the dimensions. Each parameter is input in turn and at the end of this sequence, alterations can be made if necessary using the syntax:
<parameter>=<value>
the default being <none> indicating no changes or no further changes. When all dimensions have been input, the user must use the mouse to select Calculate from the menu bar. To return to the non-graphic portion of PKSOL, click on Return.
To see graphical images of the different geometry options, press the F1 key while one of the options is highlighted (holding the mouse button down). A graphical display will appear. If all geometries do not show up in the graphical screen, click on Next or Previous to see the rest. Click on Return to return or Hardcopy if you wish to have a hard copy plot.
Figure 7‑6 Standard Specimen Options
Figure 7‑7 Standard Specimen Geometries
Figure 7‑8 Cracks at Holes Options
Figure 7‑9 Cracks at Holes Geometries
Figure 7‑10 Elliptical, Semi-elliptical Cracks in Plates Options
Figure 7‑11 Elliptical, Semi-elliptical Cracks in Plates Geometries
Figure 7‑12 Cracks at Corners Options
Figure 7‑13 Cracks at Corners Geometries
Figure 7‑14 Cracks in Solid Cylinders Options
Figure 7‑15 Cracks in Solid Cylinders Geometries
Figure 7‑16 Cracks in Hollow Cylinders Options
Figure 7‑17 Cracks in Hollow Cylinders Geometries
Figure 7‑18 Cracks in Welded Joints Plate Options
Figure 7‑19 Cracks in Welded Plate Joints Geometries
Figure 7‑20 Cracks in Welded Tubular Joints and
Cracks at Spot Welds in Tension Geometries
Once one of these geometries is selected, a graphical screen is displayed which allows for the definition of the local geometry dimensions. The main form (general options) is explained in Module Operations (Ch. 17). However, the Define, Calculate, and Return options on the main form are specific to this application.
Figure 7‑21 Graphical Definition of a Compliance Function
 
Option
Description
Define
Allows for input of the dimensions through the Command line. Dimensions appear to the right as they are defined. This option must be used first.
Calculate
Calculates the compliance function and puts the user into the post‑analysis options.
File Return
Returns the user to the previous menu without calculating the compliance function. This option is actually found under the File pull-down menu.
If the user picks option A (user parametric definition) from Figure 7‑5, a form will be presented that asks for the constants in the equation:
(7‑2)
Figure 7‑22 User Parametric Definition
This form allows the user to define a compliance function using a 5th order polynomial based on results from test or FEA, or from standard handbooks. The polynomial uses a non-dimensional crack size derived by dividing crack size (a) by total possible crack size (T).
By following the steps described in Table 7‑1 and the on-screen instructions, a results screen is reached. For example, when calculating a critical crack size, a result display, such as in Figure 7‑23, will result.
Figure 7‑23 Execution Results for Critical Crack Size
Post Analysis Options
After successful completion of any of the options from Figure 7‑4, the user will be presented with some post-processing options as shown in Figure 7‑24. If options 1‑ 3 on the calculation menu (Figure 7‑4) are chosen, the postprocessing options in Figure 7‑24 will be presented. If option 4 is chosen, then the post-processing options in Figure 7‑25 will be presented.
Figure 7‑24 Postprocessing Options
Figure 7‑25 More Postprocessing Options
 
Field
Description
Plot (any)
Plots the requested data for the compliance function. It also writes an x-y data file with the name filename.xyd and a filename.tem (template) file in the form of a stress intensity factor vs. crack size when a Y-function is calculated. An example of this type of display is shown in Figure 7‑27. With these files, the user has access to the data for plotting within the pre/postprocessor from the MSC Fatigue solution parameter form when the analysis type is set to crack growth.
Output Y Table
The Y-function table may be displayed to the screen or saved as an ASCII file for use in subsequent report preparation. If an output file is specified when this option is invoked, a file called filename.lst will be saved. If only the display to the screen is required press <RETURN>. An example of part of a Y function table listing is shown in Figure 7‑26. It shows the value of a compliance tabulated as a function of the normalized crack size (a/W). The five columns provide up to 3 decimal places resolution for the a/W parameter.
Calculate a new
compliance function
This option will return the user to the general options screen as shown in Figure 7‑4.
eXit
This option will exit PKSOL and return the user to either the system level or back into the pre-/postprocessor if that is where the user invoked PKSOL.
Figure 7‑26   Function Tabular Display
When a plot is requested within PKSOL, a screen similar to the one below appears. Various utilities are available from the main form. These are explained fully in Module Operations (Ch. 17).
Figure 7‑27 Y Function Plotted as a Function of Crack Size for Constant Stress
PKSOL K Solution References
 
Table 7‑2
Name
Reference
SENT
Murakame
SENB
Murakame
DENT
Murakame
CCP (L>6W)
Murakame
CCP (square)
Murakame
3SENB4
Murakame
3SENB8
Murakame
CTS
Murakame
Round CTS
Murakame
WOL
Walker and May
Single crack at hole
Gray
Double crack at hole
Gray
Surface crack at hole
Newman and Raju
Surface crack in tension
Newman and Raju and PD6493
Surface crack in bending
Newman and Raju and PD6493
Embedded crack in tension
Newman and Raju and PD6493
Embedded crack in bending
Newman and Raju and PD6493
Surface crack in combined tension and bending
Newman and Raju and PD6493
Embedded crack in combined tension and bending
Newman and Raju and PD6493
Quarter elliptical corner crack in tension
Newman and Raju
Quarter elliptical corner crack at a hole
Newman and Raju
Circumferential crack in tension
Gray
Straight crack in tension
James
Semi circular crack in tension
James
Crack at thread in tension
James
Straight crack in tension
James
Semi circular crack in bending
James
Internal surface crack under hoop stress
Newman and Raju
Circumferential crack in thin-walled tube
Murakame
Weld toe surface cracks in tension
PD6493
Weld toe surface cracks in bending
PD6493
Weld toe embedded cracks in tension
PD6493
Weld toe embedded cracks in bending
PD6493
Surface cracks in combined tension and bending
PD6493
Cracks in welded tubular joints
Kam
Cracks at Spot Welds
Smith
 
Murakame
 
“Stress Intensity Factors Handbook” (Ed. Y. Murakame) (Pergamon) 1990 (ISBN 0-08-034809-2)
Walker and May
 
“Compliance Functions for Various Types of Test Specimen Geometry” E. F. Walker and M. J. May, BISRA Open Report MG/E/307/67, 1967
Gray
 
“Convenient Closed Form Stress Intensity Factors for Common Crack Configurations” T.G.F. Gray, Int. Jnl. Fracture 13, 1977, p65
Newman and Raju
 
“Stress Intensity Factor Equations for Cracks in Three‑Dimensional Finite Bodies” J.C. Newman, Jr. and I.S. Raju ASTM STP 791, 1983, pI‑238
“An Empirical Stress Intensity Factor Equation for the Surface Crack” J. C. Newman, Jr. and I. S. Raju, Eng. Fract. Mechs. 15 1981, p185
PD6493
 
“Guidance on Methods for Assessing the Acceptability of Flaws in Fusion Welded Structures” (BSI) 1991 (ISBN 0‑580-19634-8)
James
 
“Review and Synthesis of Stress Intensity Factor Solutions Applicable to Cracks in Bolts” L.A. James and W. J. Mills, Eng. Fract. Mech. 30, 1988, p641
Kam
 
Dr. Jacob Kam, University College London, priv. comm.
Smith
 
Professor R. F. Smith, University of Sheffield priv. comm.