The Motivation and Method to Study Cracks in Structures (Part 1 of 3 in a series on Fracture Mechanics)

Before we jump into the topic at hand, I’d like to introduce myself. My name is Alex Austin, and I am the Structural Team Lead at DRD Technology, an Ansys Channel Partner. I studied Mechanical Engineering at the University of Tulsa, OK, from which I graduated with a BS and MS in Mechanical Engineering a little over a decade ago (woo… it’s already been that long!). My graduate work was in the fatigue and fracture space. My primary area of expertise is structural mechanics; as many engineers may know, this is quite a large field of physics when we look at Ansys simulation capabilities. Fracture mechanics is a small part of that overall field, is relatively new in the world of engineering, and is very complex. 

What is Fracture Mechanics? 

When engineers evaluate stress in a structure, the common, simple equation that comes to mind is stress = force/area. This equation carries several assumptions: static equilibrium, uniform cross-sectional area, uniaxial stress, to name a few. With the introduction of a crack to the structure, the state of stress at the crack tip is not uniaxial. Cracks are sharp corners, or notches. In the computer simulation (finite element) world, we call these singularities. In fact, singularities are locations where the theoretical stress is infinite. A strength of materials approach does not account for these singularities. When a crack exists, we need a method to analyze it. Fracture mechanics is that method. Fracture mechanics is the study of crack propagation in materials. 

Image Source: Wikipedia 

Motivation for Fracture Mechanics 

Often, cracks naturally form during the manufacturing process, either through casting or machining methods. Cracks can exist in a product and never cause issues with the working of the structure. In fact, cracks may be invisible to the naked eye. However, when this is not the case, what happens? 

Let’s say we’ve designed and manufactured structure that is currently out in the field, and our customer notices a crack… is this a problem? Let’s say a customer reported cracks popping up on some rotating machinery housings and they simply asked, ‘Is this a problem?’ though, the more likely case is no questions and, ‘Please fix this!!!’. What is the engineer typically tasked with determining? The common questions we ask are: 

  1. Will the crack grow? If so, 
  1. How quickly will the crack grow? And then, 
  1. Will the structure fail catastrophically? 

As stated, the customer absolutely thinks the presence of a crack is a problem. This is commonly the case when the customer is not an engineer, and even if they’re an engineer, they had no insight into the design and manufacture of the product. Answering the above questions will directly determine if the crack is or is not a problem. 

What about the case where the engineer designs a structure and during the design process must consider a structure that has cracks? This is a common practice in regulatory bodies, namely, the FAA (Federal Aviation Administration). In this case, the engineer assumes a crack or cracks exist in the designed components and must design for this potential failure mode. This is referred to as Fatigue and Damage Tolerance. The engineer establishes inspection intervals for components based on this analysis. The maintenance crew knows how many hours the component can be used before it needs to be checked for integrity and possibly replaced. 

 

In the next blog post, we will discuss the methods to evaluate cracked structures. 

 

Customizing the Output in Ansys Mechanical With User Defined Results (UDR)

In many situations, we have seen customers ask for ways to output custom results from ANSYS Mechanical. The usual results like Total Deformation, Equivalent Stress or Equivalent Plastic Strain may not be enough for your needs. Depending on the requirements (say a specification you are designing a part to), you can create a User Defined Result to output the needed result. ANSYS already outputs various quantities via User Defined Results that can be viewed in the Worksheet. Here is a quick look at some of what is available:

These quantities are used to create your own results output. User Defined Results can be operated on in several ways. Here is an excerpt from the ANSYS Help documentation (Mechanical Applications > Mechanical User’s Guide > Using Results > User Defined Results > User Defined Results Expressions):

 

Just as a simple example, say Total Deformation is required and is not output automatically (it is, just an example). If you add an UDR to the results, then type in the expression sqrt(Ux^2+Uy^2+Uz^2), keeping in mind these expressions are case sensitive, you get to resultant deformation from all three component values. Compare this to Total Deformation.

One can also do something more complex, say safety factor calculations. If your specified safety factor is not directly related to the Yield Strength or Ultimate Strength of the material, but some factor of, an UDR can be used; constants can be created and used just like any User Defined Result in the Worksheet. An example is shown here, where a safety factor is calculated based on a value of 6,200 psi. The safety factor looks at the First Principal Stress output, computes the safety factor, then caps the display at 7. Values less than 0 psi (compression when looking at the First Principal Stress) are set to the highest safety factor allowed (7 in this case).

A small note on the equation written in the graphic, in order to display a constant value (0 or 7 in this example), it must be multiplied by the identity matrix (matrix of 1’s). If you are just using a constant for equation manipulation, the identity matrix is not required.

User Defined Results can be a powerful tool if the output from Ansys Mechanical isn’t quite tailored to your needs.

Restarts in Ansys Mechanical Can Save Time and Effort!

‘Restarting’ the process in ANSYS Mechanical Products by which a model is solved starting from a previously solved point. The previously solved point contains data for all nodes/elements in the model; there is no results mapping and interpolation when a restart is performed. This is the most accurate method for starting an analysis from a previously solved point.

By default, restarts are not kept when a simulation finishes. The user must modify the Analysis Settings in ANSYS Mechanical to keep the restarts points that will be needed. Those options are shown below:

Restarts do have limitations. For instance, modification of loads in a step before a restart point is to be used invalidates that restart point. Several tables in the ANSYS Help documentation characterize what happens to restart points if objects in the model tree are changed. This section can be found in Mechanical Applications > Understanding Solving > Solution Restarts, as shown.

 

Here’s one of the tables showing where restarts can be used followed by common questions we related to restarts.

What can solution restarts be used for?
For cases where some subset of the loads does not change and only a few loads change/vary. This is commonly referred to as load case modeling. ANSYS Mechanical FEA can do load case modeling.

Can you give an example of load case modeling?
A model that solves bolt pretension in the first several steps and then locks them in place. Commonly, service loads are applied after the bolt tightening is simulated. The service loads can represent multiple load cases, say one load case is a pressure load, another is a force load.

Another example is a ROPS (Roll Over Protection System) analysis, where loads are sequentially applied on a structure until some criteria is met (usually energy dissipated by the structure). Once the criteria are met for one load case, the load is removed and another load is applied. This process can also be done via restarts.

Why would we do this?
This saves time on solving models. Using the example above, without restarts, the bolt pretension steps would need to be solved every time a load case is added/modified. With restarts, loads can just be activated/deactivated in a step following the bolt pretension and the bolt pretension final step used as the restart point.

Things to know before attempting this method:
All of the load cases that will be solved should be known beforehand. The restart analyses require that ALL loads are defined in the initial model. In ANSYS Mechanical, loads cannot be added to the model tree without invalidating the previous results (ignoring the ability to use restarts).

Newly added loads also will not be applied in a restart analysis as the restart method does not create new elements for these new loads.

If bolt pretension is to be used, any loads defined for a load case should be applied in a step after the final bolt pretension, as usual.

Also, if bolt pretension is used, a step after the final pretension occurs IS REQUIRED. If this is not included, ANSYS assumes the last bolt pretension will be Loaded and not Locked. This means that ANSYS would modify the pretension to maintain whatever preload you have assigned, rather than applying the service loads as a working load on the bolt. Physically, it would be like tightening/loosening the bolt as the service loads are applied.

Method:
The general steps to this method are as follows:

  1. Set up a model with all of the loads applied to the structure. This includes bolt pretensions and any loads that need to be simulated for the load cases. Set up the Analysis Settings as required for the analysis (multiple load steps, Large Deflections, etc.).
  2. Modify the Restart Controls to keep all restart points once the model is solved.
  3. Modify the loads for the load case studies to be inactive.
  4. Solve the model.
  5. Duplicate the analysis, sharing the Engineering Data, Geometry and Model cells. This guarantees the model setup remains the same for all models.
  6. In the new system, activate the load for the first load case study.
  7. Use the Tools > Read Results File… and locate the file.rst from the previously run analysis. This imports the results data into the model.
  8. Modify the Restart Analysis options in Analysis Settings to restart from the end of the first simulation (the final bolt pretension loading, for example).
  9. Solve the model.

Example:
The example shown here uses an oil field fluid end model with bolt pretensions applied. There are two load cases: 1) a pressure load on the cylinder bore, 2) a force load applied on one side of the fluid end body.