Wheel with moving load in Ansys

Published Categorised as Ansys Tagged , , , , ,
FEM stress results and moving load on a steel wheel displayed in Ansys Mechanical APDL
Steel wheel with MISES stress results and applied force displayed in Ansys Mechanical APDL

About five months ago I started with a steel wheel in PrePoMax. There you can find more information about the design and model setup. Now i reused these wheel geometry and prepared a new mesh in the latest available version 2023 R2 of Ansys Workbench.

Shell FEM mesh on steel wheel in Ansys Workbench
Shell mesh (179k nodes & 60k elements)

For a fatigue calculation I need more load cases and the load must rotate around the wheel to simulate rolling. Even with a tire or a tire model it is a bit tricky to create a moving load. Moving load is typical for load on a bridge or on a floor of a building (Moving load in Wikipedia). I solved that with a Mechanical APDL programming Command.

Create a moving load on a steel wheel with a Mechanical APDL programming Command in Ansys Workbench
Apply the loads and solve the load cases with a Command (APDL) snippet in Ansys Workbench

Now i can create any number of load cases around the 360° of the wheel and i can vary the angle of the load distribution. The Angle and shape of the load distribution are dependent on the tyre and the condition of the road. Even the solution of this large number of load cases is defined inside the APDL command snippet.

force distribution of the moving load on a steel wheel displayed in Ansys Mechanical APDL
First load case with force distribution over 72 degree displayed in Ansys Mechanical APDL

With that approach we are able to simulate the dynamic radial fatigue test according ISO 3006 (Europe) or SAE J328 (US). In this test, a constant radial force Fr is applied. And the wheel is turned with a driven drum.

schematic test setup of dynamic radial fatigue test according ISO 3006
Schematic test setup of the dynamic radial fatigue test according ISO 3006

The results are available inside Ansys Workbench and so we can use all the results and post processing features. So we can also create plots of the maximum absolute principle stress. In Fatlab this is called “Pnmax: Numerically largest principal stress (between P1 and P3)”.

FEM stress results with moving load on a steel wheel displayed in Ansys Workbench
Signed maximum absolute principle stress (Pnmax) in Ansys Workbench

With the integrated “Maximum over time” and “Minimum over time” function we can get the minimum and the maximum stress for every node over all load cases. And from that we can calculate the first simple fatigue result: the stress range with Pnmax_max – Pnmax_min. Here you can already check whether the calculated stress range is greater than the allowed stress range.

FEM stress range results on a steel wheel frontside and backside displayed in Ansys Workbench
Stress range (Pnmax_max – Pnmax_min)

With a S-N curve according DIN EN 1993-1-9 we can direct calculate a fatigue stress

ratio in Ansys Workbench with UR = calculated stress range / allowable stress range at given cycle number. That is comparable with the “Utilization ratio” in Fatlab at the knee point of the SN curve.
400000 test cycles (wheel revolutions) are specified as “minimum cycle requirement for wheels in normal highway service” with material “Ferrous” in SAE J328 standard. For the highest possible fatigue class 160 we can derive an allowable stress range of 274 N/mm² at 400000 cycles.

SN curve according DIN EN 1993-1-9 for R=-1 and fatigue class 160
SN curve according DIN EN 1993-1-9 for fatigue class 160

Now the stress range plot can be converted with the allowable stress range of 274 N/mm² at 400000 cycles in a utilization ratio plot.

fatigue stress utilization ratio on a steel wheel frontside and backside displayed in Ansys Workbench
Utilization ratio = calculated stress range / allowable stress range at 400000 cycles

And with the slope of the SN curve (m1 = 3) we can also calculate a damage plot within Ansys Workbench.

fatigue damage on a steel wheel frontside and backside displayed in Ansys Workbench
Damage frontside and backside

Damage plot and Utilization ratio plot showing values above 1.0. So our calculated design would not have been passed the dynamic radial fatigue test. Now it would be possible to reduce the loads and/or optimize the critical areas.