Angle of twist from in-plane displacements in FEA

AI Thread Summary
Solid elements in finite element analysis (FEA) typically lack rotational degrees of freedom, making it challenging to directly measure the angle of twist from torsion simulations. However, in-plane displacement values can be utilized to calculate the angle of twist by dividing the displacement magnitude by its distance from the center of the cross-section. This method has been validated for elliptical bars, yielding results that closely match analytical values. While this approach works well for symmetric cross-sections like rectangles and hexagons, it may not be suitable for asymmetric shapes such as C or L sections, which require alternative methods for accurate calculations. Understanding the geometry and ensuring proper discretization are crucial for minimizing errors in these analyses.
FEAnalyst
Messages
348
Reaction score
149
TL;DR Summary
How to calculate the angle of twist from in-plane displacements?
Hi,

solid elements used in FEA don't have rotational DOFs so normally it's not possible to read the angle of twist from simulation involving torsion. Some programs allow the transformation of results to cylindrical coordinate system but it's not always the case. From what I've heard, it's possible to use in-plane displacement values to determine the angle of twist. However, I'm not sure how to actually do it. Here's an example of an elliptical bar subjected to torsion:
- X displacements:
U1.PNG

- Y displacements:
U2.PNG


How can I calculate the angle of twist from these displacement values? Or maybe it's not enough and nodal coordinates are needed as well?

Thanks in advance for your help.
 
Engineering news on Phys.org
The maximum and minimum numbers in the color bars are the maximum X or Y displacement at the node with the maximum displacement. This is shown in the sketch below:
Ellipse.jpg

Note that the maximum Y displacement is exactly twice the maximum X displacement, and that the X dimension is twice the Y dimension of the ellipse. Also, the magnitudes of the positive X & Y displacements are equal to the magnitudes of the negative X & Y displacements. This shows that the total displacement is a rotation about the center of the ellipse. Any one of the displacements divided by its distance from the center gives the rotation angle in radians.

All of this is helped by the discretization which puts nodes at the extreme X and Y coordinates. Larger elements with nodes away from those points would add some error due to the lack of nodes at the extreme X and Y points on the ellipse.

It's easier to see this if the displacement plots are shown with deformed geometry.
 
jrmichler said:
Any one of the displacements divided by its distance from the center gives the rotation angle in radians.
Thank you very much. Let's take the displacement of the node on the left. The magnitude is ##0.03882 \ mm## and the semi-major axis is ##100 \ mm## so the angle should be ##0.0003882 \ rad##. And the analytical value is ##0.000394 \ rad## so the results agree really well.

However, what about other types of cross-sections? Could I use the same method (nodal displacement divided by its distance from the center) for any other cross-section like rectangular for example ?
 
Last edited:
You should be able to if it has appropriate symmetry. Rectangular, hexagon, or hollow round tube should be good, while a C-channel would require more effort.
 
jrmichler said:
You should be able to if it has appropriate symmetry. Rectangular, hexagon, or hollow round tube should be good, while a C-channel would require more effort.
Right, asymmetric sections like C or L (probably also T) will be problematic here. Do you know how to calculate the angle of twist from in-plane displacements in such cases? Maybe different method should be used.
For example:
- cross-section:
section.JPG

- X displacement:
L section 1.JPG

- Y displacement:
L section 2.JPG

And the analytical result is ##0.0257936508 \ rad## in this case.
 
X displacement divided by distance from zero X displacement should be fairly close if the cross section maintains its shape. Similar for Y. Look at the deformed geometry plot, try some calculations, and find out.
 
Thread 'What type of toilet do I have?'
I was enrolled in an online plumbing course at Stratford University. My plumbing textbook lists four types of residential toilets: 1# upflush toilets 2# pressure assisted toilets 3# gravity-fed, rim jet toilets and 4# gravity-fed, siphon-jet toilets. I know my toilet is not an upflush toilet because my toilet is not below the sewage line, and my toilet does not have a grinder and a pump next to it to propel waste upwards. I am about 99% sure that my toilet is not a pressure assisted...
After over 25 years of engineering, designing and analyzing bolted joints, I just learned this little fact. According to ASME B1.2, Gages and Gaging for Unified Inch Screw Threads: "The no-go gage should not pass over more than three complete turns when inserted into the internal thread of the product. " 3 turns seems like way to much. I have some really critical nuts that are of standard geometry (5/8"-11 UNC 3B) and have about 4.5 threads when you account for the chamfers on either...
Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'
Scenario 1 (figure 1) A continuous loop of elastic material is stretched around two metal bars. The top bar is attached to a load cell that reads force. The lower bar can be moved downwards to stretch the elastic material. The lower bar is moved downwards until the two bars are 1190mm apart, stretching the elastic material. The bars are 5mm thick, so the total internal loop length is 1200mm (1190mm + 5mm + 5mm). At this level of stretch, the load cell reads 45N tensile force. Key numbers...
Back
Top