|
[Sponsors] |
Order of element vs Degrees of freedom of the element |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 18, 2019, 04:29 |
Order of element vs Degrees of freedom of the element
|
#1 |
Senior Member
Mandeep Shetty
Join Date: Apr 2016
Posts: 188
Rep Power: 10 |
I have been trying to study FEM for past couple of days and this question has been bugging me for a while now.
I have read that the order of the element is the order of the polynomial used to approximate/represent the field variable in that element. If we consider a one-dimensional, 2 degrees of freedom element (with 2 nodes) the polynomial would be a linear polynomial and if this element happens to have an interior node (so totally 3 nodes, and 3 DOF), then the polynomial would be a quadratic polynomial. I have also read that the degree of the polynomial we choose also depends on the degrees of freedom of the element. So, if we consider a 4 DOF element, then the polynomial would be cubic polynomial. But even though we use a cubic polynomial the element is still a linear element. Here is my confusion: Would a quadratic polynomial still be able to represent a one dimensional, 3 noded elements (i.e. with 1 interior node) with 2 DOF at each node (so a 6 DOF element)? I am asking this because usually a 1 dimension 3 noded element would be called a quadratic element...but the textbooks always assume a single degree of freedom at each node. Would this element (with 3 nodes) still be called a 'quadratic-element' if the number of degrees of freedom per node is increased to 2? |
|
June 18, 2019, 13:50 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
The degree g of a lagrangian polynomial in 1D is equal to g=n-1 where n is the number of nodes where the values are prescribed. Given therefore a shape function of a degree g you can use it for any variable you have. |
||
June 18, 2019, 14:42 |
|
#3 | |
Senior Member
Mandeep Shetty
Join Date: Apr 2016
Posts: 188
Rep Power: 10 |
Quote:
Now in the same book when the author is talking about a 1D (2 noded) beam element which has 2 DOF at each node (so totally 4 DOF) again a cubic polynomial is used. Here it is not explicitly stated that this is a cubic element. But since a cubic polynomial is used shouldn't this also be called a cubic element? |
||
June 18, 2019, 15:04 |
|
#4 |
Senior Member
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14 |
Dear granzer,
do you understand the difference between Langrange- and Hermite interpolation? This is the answer to your question. |
|
June 18, 2019, 15:15 |
|
#5 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73 |
Quote:
An element with 2 nodes define g=1 for the esplicit lagrangian polynomial but g=3 for an implicit Hermite (or Padè) interpolation. This is the typical case of the cubic spline. Of course, the resolution requires to solve an implicit equation |
||
July 13, 2019, 05:40 |
|
#6 | ||
Senior Member
Mandeep Shetty
Join Date: Apr 2016
Posts: 188
Rep Power: 10 |
Quote:
Quote:
Now knowing this I would like to refine my previous question with which I think should clear this up for me. Let's take 2 cases. Each case comparing a 1D 2-noded element and a 1D 3-noded element (ie with 1 interior node): Case 1: Each node has 2 DoF for the 2 elements. Let them be x-displacement and rotation. In this case, I know we use hermite shape function. So the approximation function for x-displacement will be given as u(x) = sum(ui * Ni) + sum(u'i * Mi) Where u(x) is the approx function for displacement, ui is the nodal displacement, Ni shape function for nodal displacement, u'i is the derivative of nodal displacement at the node, and Mi shape function for derivative of displacement. a) What about the value for rotation DOF? Will, there be an approximation function for the rotation DOF as well? b) Taking the 2 noded element, I know that this element has a total of 4 DOF. For this will we be using a cubic polynomial? Also, what will the order of the element be? (Is it a cubic element even though it has just 2 nodes?) c) Taking the 3 noded element, I know that this element has a total of 6 DOF. For this will we be using a 5th order polynomial? Also, what will the order of the element be? ( Is it a 5th order element or will it be a quadratic element as it has 3 nodes? Case 2: Lets again consider that each node has 2 DOF. But this time let the DoF be x-displacement and temperature (or x-displacement and y-displacement). a) In this case would the approximation function for x-displacement be : u(x) = sum(ui * Ni) + sum(Ti * Mi) where Ti is the nodal temperature and Mi is the shape function for temprerature DoF. b) What about Temperature? Would there be approximation function for temperature DoF too? (some thing like T(x) = sum(Ti * Si) + sum(ui * Qi) ?? c) What will the order of the element be for the 2 noded element and the 3 noded element? Please help me out with this. Thank you |
|||
Tags |
finite element method |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Gmsh] 3D Mesh conversion from gmsh-2.5.0 to OpenFOAM | Ancioi | OpenFOAM Meshing & Mesh Conversion | 17 | January 9, 2019 00:50 |
Identifying Markers in a CGNS Mesh | tjim | SU2 | 3 | October 12, 2018 02:21 |
How grid resolution for a wave number depends on the scheme order? | marinkobezu | Main CFD Forum | 9 | September 21, 2017 13:11 |
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM | kenn | Main CFD Forum | 2 | July 18, 2004 19:28 |
about the Degrees of Freedom! | KEVIN | FLUENT | 0 | March 23, 2004 09:13 |