FEM Torque
Forum rules
and Helpful information
and Helpful information
IMPORTANT: Please click here and read this first, before asking for help
Also, be nice to others! Read the FreeCAD code of conduct!
Also, be nice to others! Read the FreeCAD code of conduct!
-
- Posts: 8
- Joined: Tue Dec 06, 2016 9:12 pm
FEM Torque
I am looking to apply a rotational force to a cylindrical object without introducing orthogonal forces in FEM analysis. Is there a good way to do this?
OS: Windows 7
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.16.6706 (Git)
Build type: Release
Branch: releases/FreeCAD-0-16
Hash: f86a4e411ff7848dea98d7242f43b7774bee8fa0
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 6.8.0.oce-0.17
OS: Windows 7
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.16.6706 (Git)
Build type: Release
Branch: releases/FreeCAD-0-16
Hash: f86a4e411ff7848dea98d7242f43b7774bee8fa0
Python version: 2.7.8
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 6.8.0.oce-0.17
Re: FEM Torque
Hopefully I understand what you would like to do. The pressure constraint goes in normal direction but you would like to apply the force in tangential direction!? I'm afraid no. You can fix faces along rotated faces but not apply forces. I do even not know of a work around by edition the input file. May be some CalculiX crack knows any better if there is something in CalculiX. BTW. We had this problem a few month (may be a year) ago already but I can not find the topic any more.timgross1986 wrote:I am looking to apply a rotational force to a cylindrical object without introducing orthogonal forces in FEM analysis. Is there a good way to do this?
...
Version: 0.16.6706 (Git)
BTW: see post viewtopic.php?f=3&t=18991#p148047 in the regard of FEM and 0.16.6707
-
- Posts: 8
- Joined: Tue Dec 06, 2016 9:12 pm
Re: FEM Torque
That's correct. Perhaps I can simulate it by making a bunch of spokes coming out of the shape and applying a bunch of smaller forces that add up to the total torque? This would mess up my stress values right? I guess in a practical application the force will be acting on a bolt of some sort. Maybe I can start with that. Thank you for replying .
Re: FEM Torque
Use a octagon (or more sides) after meshing this is more or less what you have any way. The more sides the better the approximation.
"fight the good fight"
Re: FEM Torque
Calculix has the keyword *TRANSFORM for this. You can change the coordinate system for a node set from orthogonal to cylindrical. The z-coordinate is now in axial direction. The x-coordinate points in radial direction and the y-coordinate is in tangential direction.timgross1986 wrote:I am looking to apply a rotational force to a cylindrical object without introducing orthogonal forces in FEM analysis. Is there a good way to do this?
I set up a model with an cylindrical face. In order to get a node set of this face, I used a fixed constraint on it. I had to use another constraint with a zero force(pressure), to overcome the FreeCAD-restriction to resist to write an input file without forces defined. The input file has to be edited, to create a torque:
Code: Select all
*TRANSFORM,NSET=FemConstraintFixed001,TYPE=C
0.,0.,0.,0.,1.,0.
***BOUNDARY
**FemConstraintFixed001,1
*CLOAD
FemConstraintFixed001,2,0.1
For not known reason, I could not get it to work with a Gmsh-mesh. So I used a mesh from netgen.
Here are my FreeCAD testcase (actual FreeCAD 0.17 Branch from Bernd) and the modified input file.
Ulrich
- Attachments
-
- Torque_test.png (18.54 KiB) Viewed 4842 times
-
- Pad002_Mesh.inp.zip
- (229.78 KiB) Downloaded 97 times
-
- Torque_test.fcstd
- Without mesh, use netgen for it.
- (23.19 KiB) Downloaded 102 times
-
- Posts: 8
- Joined: Tue Dec 06, 2016 9:12 pm
Re: FEM Torque
Wow, thank you.
I am having some trouble implementing the code. As you can see below. I am not as familiar with modifying the input in calculix as I could be. Still a beginner really. I will need to do some more playing with python to understand this really I'm afraid.
"TRANSFORM,NSET=FemConstraintFixed001,TYPE=C
Traceback (most recent call last):
File "<input>", line 1, in <module>
NameError: name 'C' is not defined
>>> 0.,0.,0.,0.,1.,0.
...
... ***BOUNDARY
... **FemConstraintFixed001,1
... *CLOAD
(0.0, 0.0, 0.0, 0.0, 1.0, 0.0)
>>> FemConstraintFixed001,2,0.1"
I appreciate the thorough answer,
I am having some trouble implementing the code. As you can see below. I am not as familiar with modifying the input in calculix as I could be. Still a beginner really. I will need to do some more playing with python to understand this really I'm afraid.
"TRANSFORM,NSET=FemConstraintFixed001,TYPE=C
Traceback (most recent call last):
File "<input>", line 1, in <module>
NameError: name 'C' is not defined
>>> 0.,0.,0.,0.,1.,0.
...
... ***BOUNDARY
... **FemConstraintFixed001,1
... *CLOAD
(0.0, 0.0, 0.0, 0.0, 1.0, 0.0)
>>> FemConstraintFixed001,2,0.1"
I appreciate the thorough answer,
Re: FEM Torque
Maybe the Calculix version? The version bundled with Win 0.16.6706 is Calculix 2.8p2 from bConverged, maybe too old?timgross1986 wrote:I appreciate the thorough answer,
"fight the good fight"
Re: FEM Torque
@ulirch
It is what we use in FreeCADs plane rotation constraint for rotation of a constraint displacement. But it was implemented by makkemal so I do not know if it could be used for the force constraint too.ulrich1a wrote:Calculix has the keyword *TRANSFORM for this. You can change the coordinate system for a node set from orthogonal to cylindrical. The z-coordinate is now in axial direction. The x-coordinate points in radial direction and the y-coordinate is in tangential direction.timgross1986 wrote:I am looking to apply a rotational force to a cylindrical object without introducing orthogonal forces in FEM analysis. Is there a good way to do this?
Re: FEM Torque
The *Transform should work on the force constraint
To test this build a model with transform on a displacement on the surface of the cylinder
Then add force to the surface of the cylinder, and edit the input file.
Replace the displacement node selection name with the force node selection name in transform card
Remove displacement from boundary and check
If it works we can add it to transformable constraints permanently.
1 or 2 lines of code @kgoao
To test this build a model with transform on a displacement on the surface of the cylinder
Then add force to the surface of the cylinder, and edit the input file.
Replace the displacement node selection name with the force node selection name in transform card
Remove displacement from boundary and check
If it works we can add it to transformable constraints permanently.
1 or 2 lines of code @kgoao
Re: FEM Torque
There is no python involved here.timgross1986 wrote:I will need to do some more playing with python to understand this really I'm afraid.
What you need to do is to use the Button. There you click on write inp file and then edit inp file.
After editing close the edit window and click on run calculix.
Ulrich