I am trying to create an element T1 (linear triangle) to solve the steady, incompressible Navier-Stokes equations. Therefore, my SMSTemplate command reads as follows
SMSTemplate[
"SMSTopology"-> "T1",
"SMSNoNodes"->6,
"SMSDOFGlobal"->{2,2,2,1,1,1},
"SMSNodeID"->{"D","D","D","p","p","p"},
"SMSDomainDataNames"->{"nu -kinematic viscosity","rho -density"},
"SMSSymmetricTangent"->False,
"SMSCharSwitch"->{"output"}
];
including a linear interpolation for velocity "D" and pressure "p". In total, this leads to 9 degrees of freedom at 6 nodes per element. Afterwards, I defined subroutines for Gauss points, discretization and constitutive equations to use the SMSStandardModule "Tangent and residual". The residual vector 'p' has then a size of 9 and the stiffness matrix 's' of 9x9. But if I start AceFEM, define a boundary value problem and start the analysis with SMTAnalysis, the following error message occurs:
Does someone know this problem or have an idea? Moreover, I am unsure what the output 'Nodes: {1,4,2}' means.
If you want me to give you more of my code, just let me know.
Thank you in advance