# Interference fits, radial stress and plotting along a path **URL:** \ **Category:** General Questions\ **Created:** [March 11, 2024, 2:00pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595 "2024-03-11T14:00:26Z")\ **Posts on this page:** 12\ **Page:** 1
**Author:** ![stresser](https://avatars.discourse-cdn.com/v4/letter/s/858c86/32.png) [@stresser](https://prepomax.discourse.group/u/stresser)\ **Post date:** [March 11, 2024, 2:00pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/1 "2024-03-11T14:00:26Z")
Hello, I’m a stress/structures engineer and I regularly use ansys and msc tools for mechanical systems and interactions of those systems with housings and structures. We use a lot of contact and interference fits in our models! Just recently we’ve been base lining contact management and performance in the packages available to us, using the same test - a simple shaft in a socket with interference fit. I decided I would try to solve this using prepomax/calculix over the weekend. The software solution is compared against a classical analytical approach away from stress concentrations/risers. The setup and running was fine, once I understood the steps and generally the solve time was pretty good using 4 cores. The only stumbling block is that gmsh wouldn’t mesh the socket if I split the upper and lower surfaces using the vertexes of the axial edges. I wanted to do this to create a history plot along an edge. After a little searching I see that you can add a keyword to transform the orthogonal axis system to a cylindrical one, and it would be great if that could be defined in the GUI at some point. I need to add that keyword and re-run. We are assessing the contact pressures and the radial stress in the socket, away from the socket edges. In ansys it’s straightforward to create a construction line and plot stresses along that line. How would the forum suggest doing similar with prepomax? I tried to import a line element or axis in the step file (made using freecad) but that doesn’t seem to be permitted. I tried splitting a surface to get a radial edge, but that creates problems with meshing. I can get a history plot along an edge (say the socket axial surface edge) so maybe I could split up the socket in to separate bodies in freecad, import them and make compound meshes from the individual body cells? And maybe I haven’t found the option yet, but it would also be very useful to set fringe plot limits (for example in the attached image, show values between say 50-100MPa All help always appreciated, and maybe I can contribute back in to the group. (btw, my thanks go to the developers, I really like using this software - it’s simple and straightforward, mostly intuitive). The contact pressure results are in line with both analytical and industry standard software results… ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/9/98c55570ed83420afaec21c76d32181808e41557.jpeg) ---
**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\ **Post date:** [March 11, 2024, 2:10pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/2 "2024-03-11T14:10:30Z")
> [@stresser](#): > > I can get a history plot along an edge (say the socket axial surface edge) so maybe I could split up the socket in to separate bodies in freecad, import them and make compound meshes from the individual body cells? Yes, that could work. PrePoMax can create compound parts. Or you could partition only the surface, avoiding the need for compound creation. Plotting along a line is not currently supported but you could use ParaView for this purpose. I showed it in this tutorial: [https://www.youtube.com/watch?v=wRcTD1zSbLU](https://www.youtube.com/watch?v=wRcTD1zSbLU) I also showed there how to convert the results to cylindrical CSYS in ParaView. > [@stresser](#): > > And maybe I haven’t found the option yet, but it would also be very useful to set fringe plot limits (for example in the attached image, show values between say 50-100MPa This can be set by double-clicking on the legend and switching Min/max limit type to Manual. ---
**Author:** ![SergioP1975](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/sergiop1975/32/37_2.png) [@SergioP1975](https://prepomax.discourse.group/u/SergioP1975)\ **Post date:** [March 11, 2024, 2:16pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/3 "2024-03-11T14:16:39Z")
If I remember well you must convert the .frd CalculiX results to some Paraview compatible format such as vtk/vtu/rmed, and then use Paraview to made such kind of plot along a line. For converting the results I use an external utility called Calmed8, I found it years ago in the old Salome forum and is the best by the moment to me after testing several ones. If you need this utility send me a message, because I don´t believe that is available now as it was very old post. I know that sharing directly an executable and libraries is not very GNU/GPL or even secure, but is the best that I can offer 🙂 ---
**Author:** ![stresser](https://avatars.discourse-cdn.com/v4/letter/s/858c86/32.png) [@stresser](https://prepomax.discourse.group/u/stresser)\ **Post date:** [March 11, 2024, 2:19pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/4 "2024-03-11T14:19:48Z")
Thanks for the quick responses. I figured paraview might come in to it. OK, I will look at the tutorial, do the conversion and report back when I’m done. ---
**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\ **Post date:** [March 11, 2024, 2:23pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/5 "2024-03-11T14:23:23Z")
In the tutorial, I used FreeCAD for conversion to VTK but you can find more options here: [Converting results to ParaView format](https://prepomax.discourse.group/t/converting-results-to-paraview-format/453) I mostly use the ccx2paraview tool. ---
**Author:** ![stresser](https://avatars.discourse-cdn.com/v4/letter/s/858c86/32.png) [@stresser](https://prepomax.discourse.group/u/stresser)\ **Post date:** [March 11, 2024, 2:26pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/6 "2024-03-11T14:26:05Z")
I’ll try the ccx2para tool first then ---
**Author:** ![stresser](https://avatars.discourse-cdn.com/v4/letter/s/858c86/32.png) [@stresser](https://prepomax.discourse.group/u/stresser)\ **Post date:** [March 11, 2024, 2:28pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/7 "2024-03-11T14:28:54Z")
Ok, understand the fringe plot limits now… ![image](https://global.discourse-cdn.com/free1/uploads/prepomax/original/2X/0/0f59881deee2c2325fbed818b09c2ae8c79a3f54.jpeg) ---
**Author:** ![lucas\_bueno](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/lucas_bueno/32/1168_2.png) [@lucas\_bueno](https://prepomax.discourse.group/u/lucas_bueno)\ **Post date:** [March 11, 2024, 10:03pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/8 "2024-03-11T22:03:09Z")
> [@stresser](#): > > The contact pressure results are in line with both analytical and industry standard software results… Hi, @stresser ! Can you give more details about this model? It’s always nice to check the performance of CalculiX and PrePoMax as pre/post processor with analytical results 😃 ---
**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\ **Post date:** [March 12, 2024, 9:16am UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/9 "2024-03-12T09:16:40Z")
[https://roymech.org/Useful\_Tables/Mechanics/Cylinders.html](https://roymech.org/Useful_Tables/Mechanics/Cylinders.html) ---
**Author:** ![stresser](https://avatars.discourse-cdn.com/v4/letter/s/858c86/32.png) [@stresser](https://prepomax.discourse.group/u/stresser)\ **Post date:** [March 12, 2024, 5:13pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/10 "2024-03-12T17:13:30Z")
Hey, once I’m done with the paraview side of things I’ll come back and give some more details. I can’t really release some of what I’m doing, but I can give an idea of how well prepomax performs. And, as ANYS posted an analytical methodology we could do the comparison with hand calcs here. ---
**Author:** ![ANYS](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/anys/32/7737_2.png) [@ANYS](https://prepomax.discourse.group/u/ANYS)\ **Post date:** [March 12, 2024, 5:18pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/11 "2024-03-12T17:18:27Z")
Thanks Stresser. At the end of the link there is the tribology calculator to make it easier. [https://www.tribology-abc.com/calculators/e6\_2.htm](https://www.tribology-abc.com/calculators/e6_2.htm) ---
**Author:** ![FEAnalyst](https://yyz2.discourse-cdn.com/free1/user_avatar/prepomax.discourse.group/feanalyst/32/688_2.png) [@FEAnalyst](https://prepomax.discourse.group/u/FEAnalyst)\ **Post date:** [March 12, 2024, 6:29pm UTC](https://prepomax.discourse.group/t/interference-fits-radial-stress-and-plotting-along-a-path/1595/12 "2024-03-12T18:29:31Z")
I covered interference fits in my tutorials once: [https://www.youtube.com/watch?v=sos-6ilPIZc](https://www.youtube.com/watch?v=sos-6ilPIZc) I can recommend _Shigley’s Mechanical Engineering Design_ for hand calcs.