Torsion results don't match hand-calcs

I trying to start small and build up to a slightly more complex study. For now, I’m analyzing a simple tube in PreProMax, and I got widely different results from hand-calcs.
PreProMax says: 0.006999728 Rad (~0.4 deg)
Hand Calcs says (13deg)[Torsion of a Thin-walled Open Tube | Online Calculator] EDIT: this is the wrong formula, leaving it here to keep the thread coherent)

8mm OD
6mm ID
300mm length 316 Stainless Steel

I followed this tutorial: https://www.youtube.com/watch?v=2upPxL45OZg

The mesh is 1mm to 0.1mm, with the other settings remaining default.

I applied 500 N⋅mm to a rigid body from a reference point at one end, and made a fixed boundary condition on end-surface of the other end. Then I look at UR3 for the rotational displacement in the Z-axis.

  1. What other information should I include here to ask for help?
  2. If I upload the .pmx file to share, do I also need the .STEP file separately?
  3. Where should I start looking for my problem?

I’d like to apply more force, but I expect I’ll need a non-linear study for more deflection. I need to learn more before I can do that.

Please share the .pmx file. No need for a separate .step file. You can upload it to some hosting website like WeTransfer, Google Drive or Dropbox and paste the link here.

The formulas are for open tube, I assume that your model accounts for that. It would be easier to start from a closed pipe. Formulas from Roark’s:

K is used in θ=(TL)/(KG)

Rotations are nonlinear effects. Default iteration scheme may lead to error. Try to switch to automatic with a minimum of 10 increments for the whole period. Let us know if that fix your discrepancy.

You have to set access to “Anyone with the link” and then copy the link and share it here.

2nd attempt at a link

As I suspected, your part doesn’t have a cutout so you should use the formulas for a closed pipe. I shared the ones from Roark’s above. I would also refine the mesh - it’s very coarse in your model.

(I assume you also replied on reddit?)

You’re right - I’ve checked my math and discovered that I was using the wrong formulas. Thank you for straightening me out. It’s always nice to discover the FEA was right and the hand-math was wrong.

I’m starting with this shape/loading because I’m designing a part that will be a derivative of this. It will get various slots/holes/changes in wall thickness.

As deflections get larger, I’ll need to move to non-linear analysis, no?

No, that wasn’t me.

Right but let’s correct the linear analysis and get a sufficient agreement with the analytical approach first.

I got 0.00693 rad from the formulas above.

I’ve been interupted a few times, but I got there:

this SMath sheet shows that PreProMax and hand-calcs agree.

1 Like

I added PEEQ, and still got good results (0.0698 instead of 0.0699, slightly closer to hand-calcs)

When I increased the torque 10x and re-ran the study, I noticed this in the output:
||b|| getting too small with scaling, boost it statically

what does b mean in this context?

So you added plasticity to the material definition ? Because PEEQ is just an output variable for plasticity.

It’s an array used by PaStiX, you would have to examine the ccx source code to find its meaning.

yes, I added 1 point of plasticity for the material.

My study didn’t converge though. I killed it after 27 min.

It’s better to cover a larger range with plasticity and add some hardening. But of course, the values and load magnitude have to be reasonable.

I’m struggling with the details of how to do this, and what sources to trust.

My understanding is that we’re trying to put in more data points in the relevant region of the stress/strain curve. I assume the software does linear interpolation between them.

Stress up to yield is assumed to be linear, so we’d enter Yield-stress and 0? or 0.2%?
Then we need some additional points along the curve in the plastic region. Should I be grabbing these from a stress/strain curve?

I’m currently working with 316 Stainless Steel (1.4401 in some parts of the world). I found this curve from as figure 9 in this paper.

I tried to grab some values from first part of the graph
205 MPa, 0 mm/mm
310 MPa, 0.01 mm/mm <–QUESTIONable measurement!
675 MPa, 0.2 mm/mm
900 MPa, 0.44 mm/mm

Is there a more accepted method for acquiring / entering this data?

Check this thread: Help Non Linear Aluminium FEM Beginner - #6 by FEAnalyst

Basically, you need true stress vs true plastic strain data points. Here you can find some helpful information, including formulas to calculate true stress/strain from engineering stress/strain and plastic strain from total strain: https://classes.engineering.wustl.edu/2009/spring/mase5513/abaqus/docs/v6.5/books/gsx/default.htm?startat=ch05s02.html

You can use a bilinear model (hardening recommended) like this, for instance:

Re ——- 0
Rm —— A

where Re is yield strength, Rm is tensile strength (UTS) and A is elongation at failure.

Or you can specify more data points to define a while nonlinear stress-strain curve (keeping in mind what I said at the beginning about the required stress and strain measures).

Between them - yes. But outside of the specified data, it assumes plateau (horizontal line). You should avoid this due to potential convergence issues - do the linear extrapolation manually by specifying higher stress-strain values as the last point.

Special attention to this @JKnight ! Some people may be familiar with Nastran in which some non-linear solutions use engineering stress and strain… In CCX we need to use the true ones. The results can be very different.

Thank you for clarifying. I believe I can use the curve in the top-half of the image I posted. I wish I had more of the data they used to create the graph.