Tresca应力 — CAE用语解说
Tresca应力
FEM结果中经常看到von Mises应力,Tresca应力有什么不同?
Tresca应力的理论基础
Tresca应力的定义和物理意义
教科书中出现了"Tresca应力"这个词。书中说它是主应力差最大值的简单表示,但为什么这个量如此重要?
好问题。Tresca应力是基于材料屈服的"最大剪切应力说"的等效应力。从物理角度看,它直接关联到材料内部任意平面上作用的剪切应力的最大值。例如,在拉伸试验中,材料在最大剪切应力面上沿45度方向发生滑动变形。Tresca应力
有三个主应力,为什么中间的
确实违反直觉,但Tresca准则仅关注"最大剪切应力"。剪切应力的大小由主应力对的差决定。从三个主应力可以得到三个最大剪切应力:
与另一种等效应力von Mises的根本区别是什么?我有点搞不清该用哪种。
根本区别在于屈服的物理机制假设。Tresca认为是"最大剪切应力",von Mises认为是"剪切应变能"导致屈服。在实用上,Tresca的屈服曲面是六角柱,von Mises是圆柱形。Tresca通常给出更保守(更安全)的预测,特别是在剪切主导的破坏、延性金属但需考虑剪切破坏的情况,或被规范(如某些压力容器规范)指定时使用。计算成本上Tresca略低一点。
Tresca应力的数值计算手法
FEM中的Tresca应力计算和课题
FEM求解器如何在每个单元的高斯点计算Tresca应力?特别是主应力的排序方法。
内部步骤如下。首先,从Cauchy应力张量
主应力排序中,数值误差会导致顺序变化吗?例如
这是个重要指出。确实会发生。例如,在近乎等向压缩状态(接近静水压)下,三个主应力非常接近。在数值误差范围内,顺序会变得不稳定。商用求解器对此进行了防范,通常在比较运算中引入微小容差(如1.0e-12)。但后处理中看应力云图时,微小棋盘纹理的出现是这个原因之一。
Tresca应力也用于弹塑性分析的屈服判定吗?在这种情况下,求解器如何内部处理?
确实如此,用作屈服函数。例如选择"Tresca屈服条件"的材料模型时,求解器在每个积分点每步检查计算的Tresca应力是否超过单轴屈服应力。如果超过,应力被投影回屈服曲面上(返回映射)。这比用von Mises时计算更复杂一些。因为Tresca的屈服曲面有角(边角),在角的区域流动法则方向不唯一,需要特殊算法(角处理)。Ansys的`TB, TRESCA`命令等可以指定。
Tresca应力的实务应用
分析工作流和应用事例
请告诉我实际设计工作中应该评估Tresca应力的具体情况。
有几个典型情况。第一是压力容器和管道设计。ASME锅炉及压力容器规范第VIII卷第2章沿特定评估路径计算"应力强度",这实质上就是Tresca应力(最大主应力与最小主应力的差)。第二是金属塑性加工(锻造、挤压)分析,剪切变形占主导,Tresca准则更准确表现材料流动。第三是地基和岩体剪切破坏分析。
没想到ASME规范中指定了Tresca。那么从线性静分析结果计算Tresca应力,然后直接与许可应力比较就可以了吗?
没那么简单。ASME第VIII卷第2章的"设计-by-analysis"首先要将应力分为一级、二级、峰值。评估一级一般膜应力
评估Tresca应力时,结果解释需特别注意什么?
最大要点是"应力集中部"的解释。尖锐缺口根部会计算出非常高的Tresca应力,这反映了局部多轴应力状态。但Tresca和von Mises都不能直接用该值评估延性材料的静强度。这时需要采用"结构应力"或"缺口应力"概念。另一点是接触面的结果。接触压和摩擦剪切组合导致复杂的主应力状态。Tresca应力云图在接触边界处出现不连续,可能是网格依赖性或数值微分的影响,需进行网格敏感性研究。
Tresca应力的软件比较
主要CAE软件的输出和设置
请讲Ansys Mechanical、Abaqus/CAE、COMSOL输出Tresca应力的具体操作和变量名的差异。
各有各的特点。首先**Ansys Mechanical**在Solution分支中选择"Stress"→"Stress Intensity"。默认的"Stress Intensity"实际就是Tresca应力。Von Mises是"Equivalent Stress"的单独项。结果浏览器中变量名是`SINT`。**Abaqus/CAE**中在Visualization模块"Result"→"Field Output"打开,选择`TRESC`。或在每个图表中从"Primary Variable"对话框选择`S: TRESC`。**COMSOL Multiphysics**默认有"von Mises stress"分量,但没有Tresca。需在"Definitions"中自己定义"Maximum Shear Stress"变量(即
Ansys中"Stress Intensity"是Tresca这点很容易混淆。为什么这样设定?
这与历史背景及ASME规范术语有关。"Stress Intensity"是ASME规范定义的术语,Ansys特别为压力容器行业的需求采用了这个名称。因此,仔细阅读Ansys文档会看到"Stress Intensity (Tresca)"的明确说明。用户若不知道这个事实,会混淆自己看的是von Mises还是Tresca。养成用变量名(`SINT` 对 `EQV`)确认的习惯很重要。
这些软件中,选择Tresca屈服条件作为材料模型进行弹塑性分析可行吗?
可行,但软件实现状况不同。**Abaqus**最充实,在`*PLASTIC`选项中指定`HARDENING=TRESCA`。**Ansys APDL**可用`TB, TRESCA`命令激活。但**Ansys Workbench** GUI(Engineering Data)的标准材料库没有Tresca模型,通常需要用户插入APDL命令。**COMSOL**的"Plasticity"材料模型界面中可直接选择"Tresca"作为屈服条件。无论如何,与硬化律组合的设置都比von Mises需要更多注意。
Tresca应力的故障排除
计算结果的常见问题
查看Tresca应力云图时,在尖锐角部或接触边界出现比von Mises更极端的"尖峰状"值。这在物理上正确吗,还是数值问题?
多数情况是数值问题的迹象。特别是当尖峰仅出现在1~2个单元时要怀疑。主要原因是主应力排序的不稳定性。在角部,应力张量分量急剧变化,相邻积分点中
用Tresca应力的弹塑性分析中,收敛变差的情况。具体什么情况会发生,如何对策?
收敛恶化的最大原因是前面提到的屈服曲面的"角"。应力点在这个角的区域内(或通过)时,塑性流动方向(塑性势的梯度)不唯一,求解器的迭代计算会振荡,或被迫采用极小增分步。特别在从近乎等向应力状态(静水压支配)加载剪切时易发生。对策:1.进一步细化增分步(Abaqus的`*STATIC`中`INITIAL=1e-5`等),2.适度放松求解器收敛容差(需谨慎),3.实际上,在用户子程序中实现"圆角化Tresca准则"。商用求解器内部也进行微小圆角处理。
虽然是线性静分析,但Tresca应力明显超过屈服应力的部位却实际产品正常运转。这个矛盾怎么理解?
这是CAE初学者必然遇到的壁垒。主要有三个原因。第一是"应力集中的局部性",极小区域的超过被周围弹性域约束,不导致屈服或破坏(塑性区闭合)。第二是"多轴应力状态的影响",Tresca和von Mises等效应力仅是"换算为单轴应力"的值。实际上高静水压时,某些材料的剪切屈服应力本身会上升。第三最重要的是"设计评估思想"。实务中,线性静分析的峰值应力作为疲劳评估或弹塑性分析的输入,不直接与许可值比较。例如汽车车身刚性分析中,即使Tresca超400MPa的部位,若不是持续部位就"不评估"。始终需要将分析目的和评估基准配套考虑。
相关话题
详细信息
错误