棱柱单元 — CAE术语解说
棱柱单元的理论基础
基本概念和控制方程
棱柱单元在教科书中有时也被称为"楔形单元"，但具体是什么形状呢？
正如其名，底面为三角形，沿高度方向拉伸而成的六面体。准确地说，它由2个三角形面和3个四边形面组成。例如，一次棱柱单元由底面三角形的3个节点和沿高度移动的3个节点（共6个节点）定义。二次单元在边的中点处添加节点，共有15个节点。
为什么要使用这种三角形和四边形混合的复杂单元形状？使用六面体单元不就可以了吗？
这是个很好的问题。主要有两个原因。第一，复杂几何结构（特别是曲面）的网格划分时，三角形单元对曲面的追随性强。但三角形单元存在剪切锁定问题，特别不适合板壳结构的弯曲分析。通过结合三角形网格的面内灵活性和棱柱单元厚度方向的直线性，可以精确捕捉薄壁结构的弯曲行为。
厚度方向的直线性是什么意思？
在板壳结构中，厚度方向的应变分布被假定为线性的（如Mindlin板理论）。棱柱单元的高度方向具有直线边，因此能够自然地表现这种厚度方向的线性变形模式。相比之下，用四面体单元建立板的模型时，厚度方向的节点配置不规则，难以表现这种线性分布，容易引起过度刚性（锁定）。棱柱单元的形状函数可以分离为面内方向（ξ, η）和厚度方向（ζ），这是一个重要优势。
形状函数可以分离有什么具体好处？
计算效率和积分精度得到提高。例如，一次棱柱单元的形状函数可以表示为 $$ N_i(\xi, \eta, \zeta) = L_i(\xi, \eta) \cdot M_i(\zeta) $$ 的形式，其中 $$ L_i $$ 是三角形部分的面积坐标，$$ M_i $$ 是高度方向的线性函数。这样，单元刚性矩阵的数值积分可以分离为面内方向的三角形积分和厚度方向的高斯积分，从而降低计算成本。特别是对于厚度方向材料特性变化的层积复合材料，可以高效地在厚度方向配置积分点。
棱柱单元的数值计算方法
离散化和求解器设置
棱柱单元的积分点数应该怎样确定？与四面体和六面体有区别吗？
标准配置中，一次单元通常采用面内1个积分点、厚度方向2个高斯积分点。这被称为"减积分"，具有缓解剪切锁定的效果。Ansys Solid186（20节点二次单元）的默认设置是完全积分3x3x3=27个点，但对于壳状结构，推荐选择"URI（均匀减积分）"选项，将面内2x2、厚度方向2个积分点的共8个点作为积分点。过度的减积分会产生零能量变形（沙漏）的风险，需要谨慎。
除了减积分，还有其他缓解剪切锁定的方法吗？
有的。代表性的方法是"假设自然应变法（ANS）"和"增强假设应变法（EAS）"。Abaqus的C3D6（一次棱柱）和C3D15（二次棱柱）单元默认已经施加了剪切锁定对策。特别是对于用实体单元建立复合壳体模型的情况，Abaqus通过选择启用了"非协调模式"的单元类型，可以大幅提高弯曲变形的精度。在这种情况下，单元内部引入了额外的变形自由度。
与四面体网格相比，棱柱网格对整体刚性矩阵的带宽和条件数有什么影响？
一般来说，在相同节点数的情况下，棱柱网格的连接性更好（每个节点连接的单元数更少），结果是刚性矩阵的带宽比四面体网格更小。这有助于直接法求解器的计算时间缩短。但是，棱柱网格的单元宽高比（厚度方向长度与面内边长的比）容易变得非常大。当宽高比超过1000时，单元为极薄的形状，矩阵的条件数恶化，迭代法求解器的收敛性下降。MSC Nastran为缓解这个问题提供了缩放选项。
棱柱单元的实务应用
工作流程和检查列表
请教我在实际分析中使用棱柱单元的典型工作流程。
以飞机复合材料面板或汽车点焊周边的应力集中分析为例。首先，从CAD导入中面几何，生成壳网格（主要为三角形）。接下来，使用"拉伸"或"实体网格生成"功能，将壳网格沿厚度方向拉伸为若干层棱柱单元（如复合材料，每层纤维方向为一层）。此时，厚度方向的单元数至少要2层，最好是3层以上。如果只有1层，剪切应变被假定为常数，精度会下降。
厚度方向的单元尺寸与面内单元尺寸的比例应该是多少？
这是一个非常重要的问题。根据一般经验法则，宽高比（面内边长/单元厚度）应控制在10:1到20:1以内。航空航天领域的薄壁结构有时可能达到100:1以上，但这种情况下，单元公式和积分方法的选择对结果有很大影响。Ansys手册建议，静分析可以达到1000:1，但特征值分析和动分析应该控制在100:1以下。需要通过网格敏感性分析来确认应力和位移对宽高比不敏感。
棱柱网格与其他单元（四面体或六面体）混合的模型会有问题吗？
在接口部分需要小心。例如，用棱柱单元建立薄壁部分的模型，用四面体建立肋或支架等厚壁部分的模型时，如果两者的界面节点重合，可以直接连接。但当棱柱单元的三角形面与四面体的三角形面连接时没问题，而棱柱的四边形面与四面体连接时，软件会在内部对四边形面进行"三角形连接"处理。此时局部刚性可能变化，在应力集中的连接部分，最好采用统一的单元类型（如全部使用二次单元）作为最佳实践。
棱柱单元的软件对比
Ansys / Abaqus / COMSOL等
主要CAE软件中棱柱单元的名称和特点有什么不同？
各公司的命名和实现方式不同。Ansys Mechanical中，一次单元是"SOLID185"，二次单元是"SOLID186"。Abaqus/Standard中，一次完全积分是"C3D6"，一次减积分是"C3D6R"，二次单元是"C3D15"。COMSOL Multiphysics中有"棱柱"单元分类，可以选择从一次到三次的拉格朗日单元。特别是COMSOL在电磁场和流体分析中也广泛使用棱柱边界层网格。
哪个软件在复合材料分析中比较强？
在按层定义材料和纤维配向角的"复合实体"功能方面，各公司都在努力。Ansys Composite PrepPost (ACP) 专注于前处理，可以直观地进行基于棱柱单元的分层网格生成和破坏准则（Tsai-Wu、Puck等）评估。Abaqus/CAE也有"复合铺层"功能，支持使用C3D6R和C3D15单元的分析。MSC Nastran的"PCOMP"属性卡也是业界标准之一，可与SOL 400等非线性求解器联动，追踪复合材料的损伤扩展。
在CFD（流体分析）中如何使用棱柱单元？
在CFD中"棱柱边界层网格"是必要技术。在壁面附近流速梯度急峻的区域（边界层），将多层棱柱单元垂直于壁面堆积排列。例如，Star-CCM+的"棱柱层网格生成器"和Fluent Meshing的"充气"功能都属于这一类。第一层的厚度由无量纲壁面距离y+的值（如y+≤1）反推确定。这样可以大幅提高乱流模型（SST k-ω等）的精度。流域内部通常用四面体或多面体单元填充，形成与边界层网格混合的混合网格。
棱柱单元的故障对应
常见错误和对策
分析运行时出现"负雅可比"或"过度变形"错误。棱柱单元特有的原因是什么？
棱柱单元常见的原因是宽高比极端过大，以及拉伸方向与曲面法向偏差大。后者是指，在复杂曲面上进行网格拉伸生成时，相邻单元的拉伸方向可能不连续，导致单元出现"扭曲"形状。对策包括：1) 细化面内网格以提高曲面追随性；2) 将拉伸方向统一为"平均曲面法向"；3) 对问题局部使用二次单元（C3D15等）提高形状表现力。
在复合材料层间裂开分析中，结果对网格分层数过于敏感。怎样才能得到稳定的结果？
这是经典问题。用一次棱柱单元（C3D6）的1层评估层间剪切应力时，剪切应变被假定为常数，导致应力被低估。厚度方向至少需要3层以上的单元，或者使用二次单元（C3D15）的2层也能改善。此外，可以使用Abaqus的"连续体壳单元"（SC6R、SC8R）等特殊单元，它们在厚度方向有更多积分点。对于层间裂开本身的建模，目前的最佳实践是在层间插入"内聚力单元"（COH3D6等）。在这种情况下，内聚力单元的厚度设为接近零的虚拟值，并定义其牵引-分离关系。
使用棱柱边界层网格的CFD分析收敛性变差了。可能的原因和需要检查的要点有什么？
首先确认"网格品质"。棱柱层的增长率（相邻层的厚度比）是否超过1.2？增长过快会导致数值扩散加大。其次，边界层网格与内部核心网格（四面体等）的接口处，单元尺寸是否急剧变化。尺寸比超过5倍会产生数值"反射"，干扰收敛。最后，y+值是否适合所用的乱流壁面模型（壁函数还是低雷诺数模型）。例如SST k-ω模型在低雷诺数条件下求解时，理想的y+应该≤1，如果达不到这个条件，物理结果会错误。
相关主题
了
细节
错误