6
$\begingroup$

I'm using the Finite Element Method to solve an elasticity problem and i need to use a more precise element such as a nine node triangle. For my surprise when I run the code :

Needs["NDSolve`FEM`"]
mesh = 
  ToElementMesh[Disk[], 
    MaxCellMeasure -> {"Length" -> 1}, 
   "MaxBoundaryCellMeasure" -> 2, 
   "MeshOrder" -> 1]
mesh["MeshOrder"]

mesh = 
  ToElementMesh[Disk[], 
    MaxCellMeasure -> {"Length" -> 1}, 
    "MaxBoundaryCellMeasure" -> 2, 
    "MeshOrder" -> 2]
mesh["MeshOrder"]

mesh = 
  ToElementMesh[Disk[], 
    MaxCellMeasure -> {"Length" -> 1}, 
    "MaxBoundaryCellMeasure" -> 2, 
    "MeshOrder" -> 3]
mesh["MeshOrder"]

trying to change the mesh order to 3 I receive the following error:

During evaluation of In[2]:= Throw::nocatch: Uncaught Throw[$Failed, Region`Mesh`RegionException[]] returned to top level.

Out[7]= Hold[Throw[$Failed, Region`Mesh`RegionException[]]]

Does anyone knows how can I build such an element?

$\endgroup$
  • 2
    $\begingroup$ mesh = ToElementMesh[Polygon[CirclePoints[5]], "MeshOrder" -> 3]; shows ToElementMesh::femmonv: 3 is not a valid mesh order. Valid mesh orders are 1 or 2. Using 1 as default order. $\endgroup$ – egwene sedai Nov 24 '17 at 18:54
  • 2
    $\begingroup$ I'm voting to close this question as off-topic because the OP is asking for functionality that is not supported in Mathematica. $\endgroup$ – m_goldberg Nov 25 '17 at 11:19
  • $\begingroup$ @m_goldberg I'm asking for a solution to my problem. In my problem I need this type of element. If anyone has an idea on how to do code an element like this. $\endgroup$ – Diogo Nov 25 '17 at 12:41
  • 1
    $\begingroup$ A solution to this problem will be of great use not only to me, but also to the whole community. $\endgroup$ – Diogo Nov 25 '17 at 12:41
  • $\begingroup$ The implication of @m_goldberg's comment is that if it's not built-in, then someone would have to write such a solver (plus visualization routines?). I think that's out of scope for the site. OTOH, I'm thinking that maybe the Q should be left open because negative answers to some questions are useful. And if higher-order functionality is ever added, then an updated answer can be added. Or maybe someone will point out a 3rd party product. $\endgroup$ – Michael E2 Nov 25 '17 at 13:04
3
$\begingroup$

That is not possible in version 11.2. Second order is the highest order right now.

| improve this answer | |
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.