标准的3DEC单元(zones)是4节点的四面体,假设为线性位移插值函数。block zone generate high-order-tetra命令基于二次位移插值函数,将标准单元转换成10节点的高阶四面体单元 (Higher-Order Tetrahedra)。高阶四面体(Higher-Order Tetrahedra)单元可以用于不能用m-d(hex)划分的块体,比标准的四面体单元更精确,因为在单元内有额外的网格点节点。为此,在每个单元边的中点都会创建新的节点。高阶单元(higher-order zones)的建立允许在单元内和单元表面上表示二次位移场。然而,为了接触计算和绘图的目的,单元边界被一个三角形面的网格所近似。一个10节点的四面体的每个面被分成4个平面三角形。该命令只能在使用普通块体生成命令对块进行划分后使用。
model config hotetra
block create tet 0,0,0 1,0,0 0.7,0.7,0 0.5,0.35,0.7
block zone generate edgelength 1
block zone generate high-order-tetra
高阶单元的塑性解比普通的常应变四面体单元更准确。不过高阶单元与自由场边界和节理流体流动不相容。此外,高阶单元是以八个四面体的组合形式来存储和绘制的。默认情况下,每个子四面体中都有一个单独的应力值被绘制出来。然而,高阶单元的应力被投影到网格点上,以得到准确的应力结果(见block gridpoint list stress)。投射的应力也可以通过从block plot的contour-by菜单中选择stress来绘制。
使用高阶四面体可以在可变形块中获得更精确的塑性解。然而由于中间节点的存在,接触力计算的准确性实际上被削弱了。误差随着节理刚度的增加而增加(相对于单元刚度)。对于有刚度的节理模型,建议使用常规的四面体,并开启block zone nodal-mixed-discretization on。