积分点 — CAE术语解说
积分点
老师,听说积分点在FEM中很重要,但我不清楚它与节点有什么区别?
积分点的理论基础
积分点的基本概念
经常听到"积分点"这个词,但在有限元法中它具体做什么呢?与节点不同吗?
完全正确,两者是不同的。积分点是各要素内部设置的"采样点"。在这些位置计算应变和应力,然后通过数值积分求得要素刚度矩阵。例如,在二维四节点等参数要素中,标准使用
为什么需要在要素内部设置这些点?不能直接使用节点应力吗?
这是个好问题。节点位移是连续的,但节点应力在要素之间往往不连续。而积分点处的应力是直接评估要素本构关系的"一次"信息,精度最高。例如在弹性分析中,在积分点处计算
"高斯积分点",听说还有其他种类?位置是如何确定的?
高斯-勒让德积分最常见,因为它具有最高的积分精度。还有对应梯形公式的"罗伯特积分点",以及在要素角或边上放置点的"鲁棒积分"等种类。位置和权重用要素的自然坐标系(ξ, η, ζ)定义。一维两点高斯积分的位置是
积分点的数值计算方法
积分规则与要素刚度
积分点的数量应该怎么确定?是否越多精度越高?
并不一定。存在"完全积分"和"低减积分"两个概念。例如二维四节点要素的完全积分就是前面提到的2×2=4点。而低减积分只使用1×1=1点。完全积分会过于刚硬(锁定),低减积分过于柔软而有零能量变形(沙漏)的风险。软件通常采用"选择性低减积分"。
沙漏现象具体是什么,怎样防止?
在低减积分要素中,出现不贡献于剪切变形或体积变形的零能量模式,导致要素"软软的"没有刚度。主要有两种防止方法。一是添加"沙漏控制"。二是使用非线性分析中常见的"混合要素",它单独处理体积应变,引入压力场作为额外自由度。Abaqus的CPE4H或C3D8H就是这方面的例子。
要素刚度矩阵在积分点处具体如何计算?
通过以下数值积分计算。这是FEM的核心部分。
积分点的实务应用
网格与积分点的选择
实际分析中选择要素类型时,应如何考虑积分点?
首先,线性要素和二次要素差异很大。在Ansys中,线性固体要素(SOLID185)默认完全积分(8点),但可通过关键选项切换到低减积分(1点)。二次要素(SOLID186)由于高阶,默认是低减积分(14点)。一般来说,弯曲为主的情况或接触问题中,低减积分要素更有利;但当考虑体积锁定时,应选择完全积分或混合要素。
后处理时是否需要确认积分点结果?需要注意什么?
非常重要。特别是非线性分析(塑性、超弹性)中,材料的屈服或损伤是在积分点评估的。在Abaqus/Viewer中使用"探针值"功能可直接查看特定要素的积分点应力。需要注意的是,积分点值在要素内变化,同一要素有多个值。另外,极端变形会导致积分点"反转",雅可比矩阵行列式变为负,分析停止。这是"负体积错误"的原因之一。
网格细化后,积分点数简单地增加吗?对计算成本有多大影响?
完全正确,总积分点数=要素数×(每个要素的积分点数)。计算成本影响很大。要素刚度矩阵计算、材料状态更新(非线性)、内部变量保存都在积分点执行。例如,百万个要素的模型若每个要素使用8个积分点,就要处理800万个数据点,每步都如此。这会压缩内存和CPU时间。因此,在确保必要精度的同时选择合适的积分规则是优化的关键。
积分点的软件比较
各求解器的实现差异
Ansys、Abaqus、COMSOL在积分点的设置和术语上有区别吗?
区别很大。首先术语。Abaqus称"积分点",Ansys Mechanical APDL也称"积分点",但Workbench也叫"解内点"。COMSOL称"积分点"或"高斯点"。设置方法也不同。Abaqus在要素类型定义时用后缀字符指定(R:低减积分,无后缀:完全积分),如"C3D8R"。Ansys APDL用"KEYOPT"设置。COMSOL的特点是在"物理场"的"离散化"设置中可独立选择多项式次数和积分阶次。
非线性材料的评估在所有软件中都在积分点进行吗?
是的,这是标准实现。例如弹塑性材料中,每个积分点评估屈服函数
结果输出时,积分点数据直接写入文件的功能各软件有差异吗?
有差异。Abaqus可通过"*EL FILE"输出要求将积分点数据(S、E、PE、PEEQ等)保存到.fil或.odb文件。Ansys用"ETABLE"命令将积分点结果读入表格,或结合"NLIST"和"PRNSOL"。COMSOL选择"积分点"作为"数据集",从中衍生数据进行评估和导出。第三方后处理工具如EnSight和FieldView也能读取并可视化这些原始积分点数据。
积分点的故障应对
与积分点相关的错误和对策
分析中出现"积分点塑性应变发散"的错误,原因是什么?
主要有三个原因。首先,局部应变集中导致现实不可能的大变形,超过材料模型的适用范围。其次,材料参数(如硬化规则)设置错误。第三,时间增量(Δt)过大导致数值不稳定。对策包括细化网格、检查材料参数、更严格的步长"自动增量"设置(减小最大增量尺寸)等。
用低减积分要素出现"检测到零能量变形模式"的警告,能忽略吗?
绝对不能忽略。这个警告表示模型可能存在物理上无意义的变形模式,导致刚度不足。放任不管会导致应力被低估或收敛异常快。对策包括增强沙漏控制(如Abaqus的"Hourglass stiffness"),改变网格(如单层低减积分要素很危险),或在变形较大的区域使用完全积分要素或混合要素。
后处理中应力等高线图在要素内出现"斑纹",相邻要素间有间隙。这与积分点有关吗?
正是如此。这是积分点的不连续值外推和平均化到节点过程中产生的。特别是低减积分要素,积分点少,外推影响更大。有两个对策。一是改变后处理设置中"平均化"的方法(Abaqus可调整"平均化阈值")。二是提高分析本身精度,即细化网格或使用二次要素。二次要素有更多积分点,能向节点外推更光滑的应力分布。
苛刻接触条件下出现"负雅可比"错误,与积分点有关吗?
直接相关。要素极度变形导致积分点处的雅可比矩阵行列式
相关主题
价值
详细
错误