基础-ANSYS非线性分析.pdf

ANSYS非线性分析 3D-e教育 非线性结构分析简介 非线性结构的定义 非线性结构分析的类型 非线性结构分析注意事项 非线性结构分析的分析过程 分析过程简介 建模 加载求解 结果后处理 几何非线性 几何非线性简介 大应变效应 对几何非线性情况的处理方法 几何非线性实例 材料非线性 材料非线性简介 塑性分析中的一些基本原则 查看结果 材料非线性实例 状态非线性 状态非线性简介 点点接触分析 点面接触分析 面面的接触分析 状态非线性实例 3D-e教育 非线性分析简介非线性分析简介 非线性结构分析简介 非线性结构的定义 非线性结构分析的分类 非线性结构分析注意事项 章目录 小节目录 | 3D-e教育 非线性结构分析定义非线性结构分析定义 非线性结构分析定义 非线性分析的基本特征变化的结构刚度 章目录 小节目录 | 3D-e教育 非线性结构分析分类非线性结构分析分类 非线性结构分析分类 非线性结构分析 几何非线性 材料非线性 状态非线性 结构变形的大位移造成的刚度变化 材料的物理定律是非线性的 系统状态的改变造成系统刚度的变化 章目录 小节目录 | 3D-e教育 非线性结构分析注意事项非线性结构分析注意事项 非线性结构分析注意事项 1了解程序的运作方式和结构的表现行为 2简化模型 3采用足够的网格密度 4逐步加载 5合理地使用平衡迭代 章目录 小节目录 | 3D-e教育 非线性结构的分析过程 非线性结构的分析过程 非线性结构分析的分析过程 分析过程概述 建模 加载求解 结果后处理 章目录 小节目录 | 3D-e教育 分析过程概述 分析过程概述 非线性分析的基本步骤与线性分析一样, 主要包括三步: 前处理 后处理 加载求解 建立有限元模型，定义 材料非线性参数。 施加约束载荷，定义求 解模式和载荷步设置。 观察计算结果，读取应 力、应变等数据。 章目录 小节目录 | 3D-e教育 建模 建模 塑性 塑性是一种材料行为，在施加载荷的作用下，材料发生永久的变形。 在比例极限下，材料表现为线性行为。 在屈服点下，材料表现为弹性行为。 一般对其进行简化，认为这两点重合。 章目录 小节目录 | 通过材料的材性实验得到材料的应力应变关系。 3D-e教育 建模 建模 塑性理论 材料发生塑性变形的两类理论：增量理论、全量理论。 增量理论描述的是应力增量和应变增量的数学关系。 增量理论的三个基本组成部分： 增量理论 屈服准则 流动准则 硬化准则 解决材料受三向应力的是否发生塑 形变形的应力等效准则。 塑性应变沿垂直于屈服面的 方向发展。 描述初始屈服准则如何随不断发展的 塑性应变的变化而发生变化 章目录 小节目录 | 3D-e教育 建模建模 屈服准则 在单向受力状态，通过比较其轴向应力与材料屈服应力，可以确定其是否屈服。 在三向应力状态，则必须确定一个屈服准则。屈服准则是应力的单值度量， 如果结构的应力状态和屈服准则已知，则可确定结构是否出现了塑性应变。 常用的屈服准则：Von Misses屈服准则，其等效应力计算公式： 屈服准则在主应力空间的表示为： 当等效应力超出材料 的屈服极限时，发生 屈服。 章目录 小节目录 | 3D-e教育 建模建模 流动、硬化准则 流动准则从屈服准则导出，表示为塑性应变沿垂直于屈服面的方向发展。 硬化准则则描述初始屈服准则如何随不断发展的塑性应变的变化而发生变化， 即如何修正屈服面。 章目录 小节目录 | 3D-e教育 建模建模 材料性质 定义材料性质时： 首先给出弹性材料的性质（弹性模量、泊松比） 然后给出非线性材料性质 GUI：Main MenuPreprocessorMaterial 定义线性材 料的弹性模 量和泊松比 定义材料 的非线性 性质 章目录 小节目录 | 3D-e教育 建模建模 单元类型 在ANSYS中并不是所有的单元类型都支持塑性。如：BEAM3、BEAM4、SHELL63 单元都为纯弹性单元。 对于相对小的应变情况，可以考虑采用附加位移模式的不协调单元，如平 面42、实体45单元。 对于一般的大变形情况，考虑用有中间节点的单元类型，平面183、实体 186、187单元。 另外，对于塑性分析，建议使用SHELL181、BEAM188、BEAM189单元。 章目录 小节目录 | 3D-e教育 加载求解 加载求解 基本概念-载荷步、子步、平衡迭代 载荷步：最高级，在一定“时间”范围内你明确定义的。假定载荷在载荷步内 是线性变化的。 子步：载荷增量，即在每一个载荷步内，为了逐步加载可以控制程序来执行多 次求解（子步或时间步）。 平衡迭代：为保证收敛，在每一个子步内程序所进行的一系列求解（平衡迭 代）。 章目录 小节目录 | 3D-e教育 加载求解加载求解 基本概念-时间 每一个载荷步和子步都与一个明确、唯一的“时间”值相对应。因此，子步也 称为时间步。 载荷步是通过“时间”施加的，如：time 1；time 2； ，即“时间”值是作 为每一个载荷步的重点。 每一个子步都有一个唯一的“时间”值对应，也可以激活自动时间步长，让 ANSYS自动计算并控制载荷增量。 对于速率无关的静态分析，可以采用任意的单位来定义“时间”。 如果在载荷终点时间不指定任何值，则在每一个载荷步终点的“时间”认为是 载荷步数。 章目录 小节目录 | 3D-e教育 加载求解加载求解 基本概念-时间 GUI：Main MenuSolutionLoad Step OptsTime/FrequencTime and Substps 载荷步时间 子步数 加载方式 时间步长控制 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性分析的求解控制措施 GUI：Main MenuSolutionSoln Control 1、时间、时间步选项 2、方程求解器选择 3、输出数据的管理 4、设定重启动控制 5、定义收敛误差 6、控制平衡迭代次数 7、增强求解收敛 8、在不收敛情形下控制程序行为 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-自动求解控制 关于非线性分析的求解控制很多，如何进行选择，这是一个值得仔细分析的问 题。但是在大多数情况下，我们只需要比较简单的步骤就可以进行非线性分析。 在默认状态下，自动求解控制处于自动激活状态。 提供全面、自动及智能的非线性工具设置，一般能获得有效的收敛解。 一般可以解决大部分非线性问题。 推荐首先选用该控制选项进行非线性分析，如果收敛，则OK；若收敛速度 较慢、或者不收敛，通过调整求解选项去保证结果收敛。 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-自动求解控制 在求解过程中，将会出现迭代历史的图形显示。 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-结果文件选项 缺省时，程序只将最后一个子步的计算结果写入到结果文件中去。 输出结果类 型选择 为动画、结 构历史写出 更多的子步 结果 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-求解器 ANSYS提供了5个求解器： ANSYS求 解器 波前求解器（默认） （Frontal solver） 稀疏矩阵直接求解器 （Sparse solver） 雅可比共轭梯度求 解器 预置条件共轭梯度求 解器（PCG solver） 不完全乔列斯基共 轭梯度求解器 直接求解器，适用 于自由度在5万或50 万以内的求解器； 迭代求解器，能够处 理自由度在5万100 万以上的大模型 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-非线性选项 非线性选项主要作用是帮助控制收敛。 收敛增加工具 平衡迭代最大数 蠕变选项 收敛容差 二分（后 退）准则 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-收敛准则 程序将连续进行平衡迭代直到满足收敛准则、或者直到达到允许的平衡迭代的 最大数NEQIT为止。 收敛准则：通常基于力、位移中的一种或者两种。分析时可以采用缺省的收敛 准则，也可以自己定义收敛准则。 以力为基础的收敛准则提供了收敛的绝对量度； 以位移为基础的收敛准则仅提供了表现收敛的相对量度。 应当几乎总是使用力（力矩）收敛检查，同时可以使用位移（平动、转动）收 敛进行检查。对于位移，程序将收敛检查建立在当前（i）和前面（i-1）次迭 代之间的位移改变上。 章目录 小节目录 | 3D-e教育 加载求解 加载求解 求解策略-收敛准则 对于多数工程问题，采用缺省收敛准则进行分析，一般均能得到较好的结果。 对于特殊受力情形，需要修改收敛准则: 加强、或放松收敛准则 增加力、力矩、平动、转动的收敛准则 改变收敛项的范数（无限、L1、L2） 注：如果明确地定义了任何收敛准则（CNVTOL，缺省准则将“失效”。因 此，一旦定义了位移收敛检查，必须再定义力收敛检查（使用多个CNVTOL命 令来定义多个收敛准则）。 章目录 小节目录 | 3D-e教育 加载求解 加载求解 求解策略-收敛准则 GUI：Main MenuSolutionSoln Control 章目录 小节目录 | 3D-e教育 加载求解加载求解 求解策略-最大迭代次数 非线性分析的求解若不能在一个合理的平衡迭代数内收敛，程序将荷载进行二 分并再次求解。 可以指定平衡迭代次数。 程序缺省值为15到26次迭代。 设置最大迭代次数 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性求解两个收敛工具 两个加强的收敛工具 线性搜索（Line Search） 预报器（Predictor） 线性搜索对克服振荡收敛非常有效，使求解趋于稳定。但需要花费额外的计算时间。对 解决超弹性、接触问题比较有效。 缺省时，预报器打开，适用 于大多数的情形。 预报器通过对每一个子步的 第一次平衡迭代结果预报自 由度解而加快收敛。 当非线性响应比较光滑时， 预报器非常有用，并且可以 加速收敛。如塑性分析。 非线性响应不光滑时，打开 预报器会导致发散，此时须 关闭预报器 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性求解-后退选项 如果一个子步收敛失败，程序将目前荷载进行二分并再次进行求解。 缺省的后退控制适用于大多数情况。 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性求解-弧长法 对于比例加载，若采用荷载增量法，一般只能得到结构的承受的最大承载能力 （最大响应点），无法得到荷载下降点（变形最大值）。 如果预料到结构在它的载荷历史内在某些点将变得物理意义上不稳定（ 结构 的载荷位移曲线的斜度将为0或负值），这是可以使用弧长法来帮助稳定数 值求解。 和其它分析过程基本一致，不同的有如下几点： 弧长法只对力控制的比例加载分析有效，不允许使用面荷载。 荷载因子将作用于所有的荷载。 施加的荷载应略高于（20）特征值分析计算得到的临界荷载。 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性求解-弧长法 弧长法只对静力分析有效，同时激活大位移效应（NLGEOM，ON） 推荐使用缺省的求解控制。 不要使用线性搜索（LNSRCH）、预报器（PRED）、自适应下降、自动时间步 （AUTOTS，DELTIM）、时间（TIME）、或带弧长的时间积分（TIMINT）。 结果文件中应写入足够多的子步结果（OUTRES），以便在时间历程处理器中 检查位移结果数据。 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性求解弧长法 弧长的计算： Reference Arc-Length Radius=Total Load/Nsubstep 上式中，Nsubstep为指定的子步数，采用较多的子步数可使计算结果更接近真 实值，但也导致较长的计算时间。 最大弧长半径因子(10) 最小弧长半径因子(0.001) 章目录 小节目录 | 3D-e教育 加载求解加载求解 非线性求解弧长法 弧长法中的“时间”不一定是单调上升的，在一些特殊受力情形（如后屈曲） 中，可能出现负值。 此时的结果后处理须采用荷载步、子步来处理。 章目录 小节目录 | 3D-e教育 结果后处理结果后处理 后处理所要做的工作 1、收敛检查 2、通用后处理 3、时间历史后处理器 4、验证 章目录 小节目录 | 3D-e教育 结果后处理结果后处理 后处理收敛检查 在进行后处理之前，必须核实结果有一个收敛解。 验证收敛的方法很多： 1. 检查三个文件：输出、错误和监控文件 2. 通用后处理器中检查结果 3. 绘制时间历史曲线图 章目录 小节目录 | 3D-e教育 结果后处理结果后处理 在通用后处理器中，结果列表： 子步999999写不收敛的结果；相应的结果不要进行处理； 在时间历史结果曲线图中； 收敛的情况下均能得到结构在荷载作用下的光滑曲线； 而在不收敛的情况下，响应曲线“跳跃”。 章目录 小节目录 | 3D-e教育 结果后处理结果后处理 后处理通用后处理器 允许在某个时间点上查看结果。 1. 尽量采用结果概况去选择想要的结果 GUI：Main MenuGeneral PostprocResults Summary 2. 尽量少采用通过时间去选取结果避免出现如图所示的内插误差。 GUI：Main MenuGeneral PostprocRead Results-By Time/Freq 注意：非线性响应不成比例，此时 叠加原理无效, 荷载工况无效 结果比例缩放无效。 章目录 小节目录 | 3D-e教育 结果后处理结果后处理 结果后处理时间-历程后处理器 在时间历程后处理器中，通过定义临时变量，可以得到结构： 荷载挠度曲线 应力应变曲线 章目录 小节目录 | 3D-e教育 结果后处理结果后处理 后处理验证 一个没有错误信息的收敛分析并不一定就能给出可靠的结果。 对于结果的验证主要从以下几个方面考虑： 1、经验（定性） 2、标准考题（经典考题） 3、建模会导致求解不正确 网格太粗糙（查看单元等值线：不应该出现中断、丢失） 网格扭曲（通过查看变形形状图：显示高度扭曲的单元） 材料错误 边界条件不正确 后处理（大部分位移响应曲线为光滑的，若出现锯齿状，说明 时间步太大）。 章目录 小节目录 | 3D-e教育 几何非线性几何非线性 几何非线性 几何非线性简介 大应变效应 对几何非线性情况的处理方法 几何非线性实例 章目录 小节目录 | 3D-e教育 几何非线性简介几何非线性简介 几何非线性简介 如果结构经受大变形，它变化的几何形状可能会引起结构的非线性地响应。 如图 所示的钓鱼杆，在轻微的垂向作用下，会产生很大的变形。随着垂向载 荷的增加，杆不断弯曲以致于动力臂明显地减少，导致杆端显示出在较高载荷 下不断增长的刚性。一般来说，随着位移增长，一个有限单元已移动的坐标可 以以多种方式改变结构的刚度。这类问题总是是非线性的，需要进行迭代获得 一个有效的解。 章目录 小节目录 | 3D-e教育 大应变效应大应变效应 一个结构的总刚度依赖于它的组成部件（单元）的方向和单刚。当一个单元的结 点经历位移后，那个单元对总体结构刚度的贡献可以以两种方式改变变。 首先，如果这个单元的形状改变，它的单元刚度将改变。 其次，如果这个单元的取向改变，它的局部刚度转化到全局部件的变换也将改 变。 章目录 小节目录 | 3D-e教育 大应变效应大应变效应 大应变效应 大应变分析说明由单元的形状和取向改变导致的刚度改变。因为刚度受位移影 响，且反之亦然，所以在大应变分析中需要迭代求解来得到正确的位移。通过 GUI路径 Main Menu | Solution | Analysis Options 来激活大应变效应。这效应改 变单元的形状和取向，且还随单元转动表面载荷。在大多数实体单元（包括所 有的大应变和超弹性单元），以及部分的壳单元中大应变特性是可用的。在 ANSYS/Linear Plus程序中大应变效应是不可用的。 大应变处理对一个单元经历的总旋度或应变没有理论限制。然而，应限制应变 增量以保持精度。 因此， 总载荷应当被分成几个较小的步， 这可以 通过 GUI 路径 Main Menu | Solution | Time/Prequent。无论何时当系统是非保守系统， 来自动实现如在模型中有塑性或摩擦，或者有多个大位移解存在，如具有突然 转换现象，使用小的载荷增量具有双重重要性。 章目录 小节目录 | 3D-e教育 对几何非线性情况的处理方法对几何非线性情况的处理方法 对几何非线性情况的处理方法 1应力应变 在大应变求解中，所有应力应变输入和结果将依据真实应力和真实（或对数 ）应变。 2单元的形状 应该认识到在大应变分析的任何迭代中低劣的单元形状（也就是，大的纵横比 ，过度的顶角以及具有负面积的已扭曲单元）将是有害的。 3应力刚化 结构的面外刚度可能严重地受那个结构中面内应力的状态的影响。面内应力和 横向刚度之间的耦合，通称为应力刚化。 4旋转软化 旋转软化为动态质量效应调整（软化）旋转物体的刚度矩阵。在小位移分析中 这种调整近似于由于大的环形运动而导致几何形状改变的效应。 章目录 小节目录 | 3D-e教育 几何非线性实例几何非线性实例 几何非线性实例详解-细长杆屈曲分析 下图为一横截面为正方形的细长杆，一端固定，另一端受集中力F的作用，求 细长杆变形后的形状及A点的位移。 长杆几何参数：L=1200mm；B=6mm；H=6mm 长杆材料参数：弹性模量E=85GPa；泊松比=0.3。 载荷：F=20N。 章目录 小节目录 | 3D-e教育 材料非线性材料非线性 材料非线性 材料非线性简介 塑性分析中的一些基本原则 查看结果 材料非线性实例 章目录 小节目录 | 3D-e教育 材料非线性简介材料非线性简介 材料非线性 非线性的应力应变关系是结构非线性的常见原因。许多因素可以影响材料的 应力应变性质，包括加载历史（如在弹塑性响应状况下），环境状况（如温 度），加载的时间总量（如在蠕变响应状况下）。 ANSYS的材料非线性分析能力包括： 弹塑性分析 超弹分析 蠕变分析等。 章目录 小节目录 | 3D-e教育 材料非线性简介材料非线性简介 材料非线性 塑性是一种在某种给定载荷下，材料产生永久变形的材料特性。对大多的工程 材料来说，当其应力低于比例极限时，应力应变关系是线性的。另外，大多 数材料在其应力低于屈服点时，表现为弹性行为。也就是说，当移走载荷时， 其应变也完全消失。 由于屈服点和比例极限相差很小，因此在 ANSYS 程序中，假定它们相同。在应 力应变的曲线中，低于屈服点的叫作弹性部分，超过屈服点的叫作塑性部分 ，也叫作应变强化部分。 塑性分析中考虑了塑性区域的材料特性。当材料中的应力超过屈服点时，塑性 被激活,也就是说，有塑性应变发生。而屈服应力本身可能是下列某个参数的 函数：温度、应变率、以前的应变历史、侧限压力和其它参数。 章目录 小节目录 | 3D-e教育 塑性分析中的一些基本原则塑性分析中的一些基本原则 下面的这些原则应该有助于可执行一个精确的塑性分析。 1、所需要的塑性材料常数必须能够足以描述所经历的应力或应变范围内的材 料特性。 2、缓慢加载，应该保证在一个时间步内，最大的塑性应变增量小于 5%，一般 来说，如果 Fy是系统刚开始屈服时的载荷，那么在塑性范围内的载荷增量应 近似为： 0.05*Fy对用面力或集中力加载的情况； Fy对用位移加载的情况。 3、当模拟类似梁或壳的几何体时，必须有足够的网格密度，为了能够足够的 模拟弯曲反应，在厚度方向必须至少有二个单元。 章目录 小节目录 | 3D-e教育 塑性分析中的一些基本原则塑性分析中的一些基本原则 下面的这些原则应该有助于可执行一个精确的塑性分析。 4、除非那个区域的单元足够大，应该避免应力奇异，由于建模而导致的应力 奇异有： 单点加载或单点约束； 凹角； 模型之间采用单点连接； 单点耦合或接触条件。 5、如果模型的大部分区域都保持在弹性区内，那么可以采用下列方法来降低 计算时间： 在弹性区内仅仅使用线性材料特性，不使用TB命令； 在线性部分使用子结构。 章目录 小节目录 | 3D-e教育 查看结果查看结果 查看结果 1感兴趣的输出项（例如应力，变形，支反力等等）对加载历史的响应应该 是光滑的，一个不光滑的曲线可能表明使用了太大的时间步长或太粗的网格。 2 每个时间步长内的塑性应变增量应该小于5，这个值在输出文件中以 “Max plastic Strain Step”输出，也可以使用POST26 来显示这个值 （Main Menu | Time Hist Postpro Define Variables）。 3塑性应变等值线应该是光滑的，通过任一单元的梯度不应该太大。 4画出某点的应力应变图，应力是指输出量 SEQV（Mises 等效应力），总 应变由累加的塑性应变EPEQ和弹性应变得来。 章目录 小节目录 | 3D-e教育 材料非线性实例材料非线性实例 材料非线性实例圆盘大应变分析 有一圆盘，其纵截面形状如图所示。圆盘由弹塑性材料构成，其顶面承受均布 压力载荷P的作用，求圆盘的应力和位移响应。 几何参数：半径R=10mm；高度H=5mm。 载荷：P=3MPa 材料参数：弹性模量E=2500MPa；泊松比 ；屈服强度 ；剪切模量 。 35. 0 MPa s 5 . 1 MPaET6 . 3 章目录 小节目录 | 3D-e教育 状态非线性状态非线性 状态非线性 状态非线性分析简介 点点接触分析 点面接触分析 面面的接触分析 状态非线性实例 章目录 小节目录 | 3D-e教育 状态非线性简介状态非线性简介 系统的刚度由于系统状态的改变在不同的值之间突然变化。 状态改变也许和载荷直接有关，也可能由某种外部原因引起（如在冻土中的紊 乱热力学条件） 。ANSYS程序中单元的激活与杀死选项用来给这种状态的变化 建模。 接触是一种很普遍的非线性行为。接触是状态变化非线性类型形中一个特殊而 重要的子集。接触问题是一种高度非线性行为，需要较大的计算资源，为了进 行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点： 其一，在求解问题之前，不知道接触区域，表面之间是接触或分开是未知 的，突然变化的。这随载荷、材料、边界条件和其它因素而定； 其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们 都是非线性的，摩擦使问题的收敛性变得困难。 接触问题分为两种基本类型：刚体柔体的接触、柔体柔体的接触。 ANSYS 支持三种接触方式：点点，点面，平面面，每种接触方式使用的接 触单元适用于某类问题，并且每种类型的问题都有不同的特点。 章目录 小节目录 | 3D-e教育 点点点接触分析点接触分析 点点接触分析 点点接触单元主要用于模拟点点的接触行为，为了使用点点的接触单元， 需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动 的情况（即使在几何非线性情况下）。 如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动） 保持小量，那么可以用点点的接触单元来求解面面的接触问题，过盈装配问 题是一个用点点的接触单元来模拟面与的接触问题的典型例子。 在 ANSYS程序中提供了如下三种点点的接触单元： CONTAC12：2D点点的接触单元，这个单元是通过总体坐标系XY平 面内的二个结点来定义的，可以用于2D平面应力，平面应变和轴对称分 析中。 CONTAC52：3D在点点的接触单元， COMBIN40 章目录 小节目录 | 3D-e教育 点点面接触分析面接触分析 点面接触分析 点面接触单元主要用于给点面的接触行为建模，如果通过一组结点来定义接 触面，生成多个单元，那么可以通过点面的接触单元来模拟面面的接触问题 ，面即可以是刚性体也可以是柔性体。 使用这类接触单元，不需要预先知道确切的接触位置，接触面之间也不需要保 持一致的网格，并且允许有大的变形和大的相对滑动。 ANSYS程序的点面接触单元允许下列非线性行为： 有大变形的面面接触分析 接触和分开 库仑摩擦滑动 热传递 章目录 小节目录 | 3D-e教育 点点面接触分析面接触分析 点面接触分析 在 ANSYS 程序中点面的接触是通过跟踪一个表面（接触面）上的点相对于另 一表面（目标面）上的线或面的位置来表示的，程序使用接触单元来跟踪两个 面的相对位置。 接触单元的形状为三角形，四面体或锥形，其底面由目标面上的