PDF

Bracket — Fatigue Evaluation
Introduction
The S-N curve, also called the Wöhler curve, is one of the most popular methods for fatigue evaluation. The curve relates the stress amplitude to the limiting fatigue life and can be obtained directly from a set of standard fatigue test. Often, a structure is however subjected to conditions different from the experimental conditions of the fatigue tests. The fatigue data must then be appropriately modified in order to take the actual operating conditions into account.
In this example it is shown how to perform a fatigue evaluation when the material data needs to account for harsh environmental conditions and poor manufacturing process.
Model Definition
The bracket geometry can be seen in Figure 1. The component is subjected to external loads that lift one arm up and pull the other arm down. The load magnitude is cycled between a zero load, a peak load and back to zero load. Additional information regarding the model set up can be found in the documentation of the application Bracket-Spring Foundation Analysis, found in the Structural Mechanics Module.
Figure 1: Bracket geometry.
The structural steel in the bracket has poor cleanliness and contains fairly large inclusions. It is well known that with decreasing purity the endurance limit decreases and thus large inclusions shorten the lifetime. The fatigue properties of the material have been tested in a material laboratory and are summarized in Table 1. The material has a distinct endurance limit at 110 MPa. Tests were performed in nominal testing conditions and on specimens with a very good surface finish.
As opposite to the testing conditions, the bracket was machined and has a rough surface as the result. Moreover the bracket operates in salt water that has a strong influence on the fatigue resistance. Since there is lack of test data for the operating conditions, the S-N curve data must be modified. One approach is to use
(1)
where σa is the stress amplitude, k is the modification factor, N is the fatigue lifetime, and fSN is the S-N curve. Based on the experience, the modification factor for the corrosion in salt water can be set to 0.4. The modification factor for the surface finish, due to machining manufacturing technique, can be set to 0.7. The modification factor in Equation 1 for the combined effect is simply the product of the modification factors for all operating conditions and thus 0.28.
In the fatigue analysis the modification due to machining, k = 0.7, and due to the combined effect is evaluated, k = 0.28.
Results and Discussion
The stress distribution in the bracket a the peak load is shown in Figure 2.
Figure 2: Stress in the bracket at the peak operating load.
A close-up of the part with stress concentrations reveals that also on the inner side of the bracket there are significant stresses, see Figure 3.
Figure 3: A close-up of the part of the bracket that experiences highest stresses.
The fatigue analysis of the bracket material modified for the manufacturing conditions predicts infinite life. When submerged in saltwater, the stresses in most of the bracket are not high enough to cause fatigue. Only in a few local points is there is a risk of fatigue, see Figure 4. The bracket most probably fails in the connection between the two flat parts since the stresses on both the upper and the lower sides are high enough to initiate fatigue.
Figure 4: Fatigue lifetime when bracket is submerged in salt water.
In Figure 4 a finite life is predicted close to the holes that fasten the bracket. Those results are highly dependent on the boundary conditions that were applied there. In this model a spring foundations was used to fasten the bracket, but the constraint could have also been applied using a fixed connection of with a full contact analysis. This would change the stress field in the vicinity of the holes and result in a different fatigue lifetime.
Notes About the COMSOL Implementation
The S-N curve can be defined using different function types in COMSOL Multiphysics. When using an interpolation function it is important to use many points to specify the relation between fatigue life and stress amplitude, since the range of the fatigue life is very large and a small change in the stress amplitude results in a large change in the fatigue life. The function value between two data points can be evaluated in different ways but linear interpolation is the most common one. In Figure 5 the S-N curve is specified with eleven and with four measurement points respectively. Note that the fatigue life is displayed in a logarithmic scale and therefore the interpolation is not a straight line. The two curves display large differences. At 200 MPa one curve predicts 3.23·104 cycles while the other one predicts 5.80·104 cycles; a difference by almost a factor two.
Figure 5: S-N curve based different number of measurement points.
Application Library path: Fatigue_Module/Stress_Life/bracket_fatigue
Modeling Instructions
Root
In this example you will start from an existing model from the Structural Mechanics Module.
1
From the File menu, choose Open.
2
From the Application Libraries root, browse to the folder Structural_Mechanics_Module/Tutorials and double-click the file bracket_spring.mph.
Root
Fatigue is calculated based on a load cycle. Create a parameter to scale the magnitude of the external load.
Global Definitions
Parameters 1
1
In the Model Builder window, under Global Definitions click Parameters 1.
2
In the Settings window for Parameters, locate the Parameters section.
3
Component 1 (comp1)
In the Model Builder window, expand the Component 1 (comp1) node.
Definitions
Analytic 1 (load)
1
In the Model Builder window, expand the Component 1 (comp1)>Definitions node, then click Analytic 1 (load).
2
In the Settings window for Analytic, locate the Definition section.
3
In the Expression text field, type para*F*cos(atan2(py,abs(px))).
Start by improving the existing mesh in order to get a better resolution of stress gradients at the critical fillets.
Mesh 1
In the Model Builder window, expand the Component 1 (comp1)>Mesh 1 node.
Size 2
1
In the Model Builder window, expand the Component 1 (comp1)>Mesh 1>Free Tetrahedral 1 node.
2
Right-click Free Tetrahedral 1 and choose Size.
3
In the Settings window for Size, locate the Geometric Entity Selection section.
4
From the Geometric entity level list, choose Boundary.
5
6
Locate the Element Size section. Click the Custom button.
7
Locate the Element Size Parameters section. Select the Maximum element size check box.
8
9
Click  Build Selected.
Study 1
Since it is an elastic case and the external load is proportional, only two load cases are necessary to capture the stress amplitude evaluated in the S-N curve.
Step 1: Stationary
1
In the Model Builder window, expand the Study 1 node, then click Step 1: Stationary.
2
In the Settings window for Stationary, click to expand the Study Extensions section.
3
Select the Auxiliary sweep check box.
4
5
6
In the Home toolbar, click  Compute.
Global Definitions
Interpolation 1 (int1)
1
In the Home toolbar, click  Functions and choose Global>Interpolation.
2
In the Settings window for Interpolation, locate the Definition section.
3
From the Data source list, choose File.
4
Find the Functions subsection. In the table, enter the following settings:
5
Click Browse.
6
7
Click Import.
Add Physics
1
In the Home toolbar, click  Add Physics to open the Add Physics window.
2
Go to the Add Physics window.
3
Find the Physics interfaces in study subsection. In the table, clear the Solve check box for Study 1.
4
In the tree, select Structural Mechanics>Fatigue (ftg).
5
Click Add to Component 1 in the window toolbar.
6
In the Home toolbar, click  Add Physics to close the Add Physics window.
Fatigue (ftg)
Stress-Life 1
1
Right-click Component 1 (comp1)>Fatigue (ftg) and choose the boundary evaluation Stress-Life.
2
In the Settings window for Stress-Life, locate the Boundary Selection section.
3
From the Selection list, choose All boundaries.
4
Locate the Fatigue Model Selection section. From the σ list, choose Signed von Mises (principal).
5
From the Modification list, choose Stress factor.
6
Locate the Solution Field section. From the Physics interface list, choose Solid Mechanics (solid).
7
Locate the Fatigue Model Parameters section. From the fSN(N) list, choose wohler.
8
In the k text field, type k.
9
Locate the Evaluation Settings section. In the Ncut text field, type 1e6.
Add Study
1
In the Home toolbar, click  Add Study to open the Add Study window.
2
Go to the Add Study window.
3
Find the Physics interfaces in study subsection. In the table, clear the Solve check box for Solid Mechanics (solid).
4
Find the Studies subsection. In the Select Study tree, select Preset Studies for Selected Physics Interfaces>Fatigue.
5
Click Add Study in the window toolbar.
6
In the Home toolbar, click  Add Study to close the Add Study window.
Study 2
Step 1: Fatigue
1
In the Settings window for Fatigue, locate the Values of Dependent Variables section.
2
Find the Values of variables not solved for subsection. From the Settings list, choose User controlled.
3
From the Method list, choose Solution.
4
From the Study list, choose Study 1, Stationary.
Parametric Sweep
1
In the Study toolbar, click  Parametric Sweep.
2
In the Settings window for Parametric Sweep, locate the Study Settings section.
3
4
5
In the Study toolbar, click  Compute.
Results
Surface 1
1
In the Model Builder window, expand the Cycles to Failure (ftg) node, then click Surface 1.
2
In the Settings window for Surface, click to expand the Quality section.
3
In the Cycles to Failure (ftg) toolbar, click  Plot.
To get a better view of the critical points, you can rotate and zoom the bracket. For this purpose, it is useful to have a dedicated View node:
4
Click the  Show More Options button in the Model Builder toolbar.
5
In the Show More Options dialog box, in the tree, select the check box for the node Results>Views.
6
View 3D 2
1
In the Model Builder window, right-click Views and choose View 3D.
2
3
In the Settings window for View 3D, locate the View section.
4
Select the Lock camera check box.
Now apply this view to the plots.
Stress (solid)
1
In the Model Builder window, click Stress (solid).
2
In the Settings window for 3D Plot Group, locate the Plot Settings section.
3
From the View list, choose View 3D 2.
4
In the Stress (solid) toolbar, click  Plot.
Compare the resulting plot with that in Figure 3.
Cycles to Failure (ftg)
1
In the Model Builder window, click Cycles to Failure (ftg).
2
In the Settings window for 3D Plot Group, locate the Plot Settings section.
3
From the View list, choose View 3D 2.
4
In the Cycles to Failure (ftg) toolbar, click  Plot.
Compare with Figure 4.
Max/Min Surface 1
1
In the Model Builder window, right-click Max/Min Surface 1 and choose Disable.
2
Right-click Max/Min Surface 1 and choose Enable.