在有限元仿真当中,应力、应变结果是评判结构强度、安全性的核心依据。
但很多工程师会发现:直接读取单元积分点或节点的原始计算结果,并不能直接拿来做工程判定,往往需要进行应力应变外推处理。
应力外推并不是可选的优化技巧,而是由有限元计算底层原理决定的必要操作,用来解决计算偏差、单元本身缺陷和工程实际要求之间的矛盾。
有限元依靠多项式形函数做插值计算。以4节点四边形单元CPE4为例,单元内部应力被假定为恒定不变。但现实工况中,构件应力大多是连续变化的,比如受弯梁,应力沿着截面线性改变。
这种“以直代曲”的近似会带来计算偏差,单元边界位置误差尤其明显。应力外推利用周边单元数据做高阶插值,还原连续应力分布,弥补低阶形函数带来的精度损失。
有限元应力结果主要来自高斯积分点,例如CPE4单元只有2×2个积分点。积分点处的计算结果虽然精度高,但只代表单元内部局部位置。
像应力集中这类关键位置,积分点无法覆盖,通过外推,把积分点的应力结果外推到节点、构件表面,将离散点数据还原完整连续应力场,捕捉应力渐变趋势,避免关键信息丢失。
部分关键位置,受单元本身特性影响,直接计算的结果会失真,需要通过外推还原真实水平。
孔边、尖角拐角属于典型应力集中区域,应力会在极小范围内急剧升高。即使加密网格,依然会存在单元尺寸效应,出现“应力锁死”现象:仿真直接读出的应力会远低于理论真实值。
举个例子:带圆孔受拉平板,直接输出的结果会比理论峰值低30%-50%。使用线性外推法,把周边单元的应力向集中点外推,就可以还原真实应力峰值,为疲劳评估提供可靠依据。
约束、集中载荷作用的边界位置,很容易出现应力奇异。
固定约束节点因为局部刚度突变,集中力施加点会产生无意义的应力奇异值,直接读取该位置结果会误导判断。
工程上通常选取距离边界2-3个单元位置的有效数据,通过外推推算边界真实应力,避开奇异区域带来的错误数值。
有限元直接输出的离散点结果,并不完全匹配设计规范的校核要求,应力应变外推就是连续仿真计算和工程应用的重要桥梁。
获取构件真实表面应力
绝大多数零部件破坏起源于构件表面,金属混凝土等材料的强度校核标准,也以表面应力作为评判指标。
但高斯积分点大多分布在单元内部,无法直接拿到表面应力。通过外推,将单元内部积分点应力映射到构件表面,得到用于强度校核的表面应力。
满足规范中等效应力评估要求
ASME、GB/T等行业标准,需要基于截面完整应力分布计算von Misese等效应力。
有限元原始输出只是一个个离散点应力,无法直接积分处理。经过外推之后得到连续的截面应力场,才可以完成积分计算,满足规范对于整体应力状态的评估要求。
应力应变外推,不是锦上添花的可选操作,而是修正计算误差、突破单元局限、匹配工程规范的必要步骤。
它本质是依靠数学手段修正离散仿真带来的缺陷,把原始“数值解”转化为可以直接用于设计校核的“工程可用解”。
实际项目中,我们可以根据实体、壳体、梁单元的区别,结合应力集中、疲劳分析等目标,选用线性、二次、子模型外推等不同方法,同时配合网格收敛性校验,保障仿真结果可靠。