Abaqus这玩意儿,提取应变结果,那可真是头疼啊。先得弄明白积分点计算、节点插值,还得分清实体/壳单元坐标系和平均算法,不然数值准得错成狗。
S11、S22、S33那是正应力分量,实体就简单,全局X/Y/Z轴;壳单元得看面内1/2和法向3,S33默认是0。剪应力分量S12、S13、S23,壳单元一般只看S12,S13、S23都不要。
Mises应力(S)那玩意儿,第四强度理论搞出来的,评估材料屈服风险,公式看着就头疼。
压力(P)简单,就是静水压力,负值表示受压。
最大/最小主应力,判断脆性断裂或拉裂,绝对值最大的主应力,那才是真·关键。
查看结果,先打开.odb文件,Visualization模块走起。选变量,Result > Field Output,S(应力张量)或Mises,壳单元表面默认是中面或平均,想看厚度分布,Options > Contour > Limits或Section Points来指定。
探针查询,Query > Probe Values,节点/单元值一查就知,但别忘了,显示的是插值后节点值。
坐标系确认,壳单元应力用局部坐标系,1-2面内,3法向,Tools > Datum CSYS或显示局部轴,方向不能错。
平均机制,应力从积分点外插到节点再平均,Abaqus两种模式,先算等效再平均,结果偏大偏安全;先平均分量再算等效,云图平滑度受Averaging Threshold控制。
网格太粗,应力梯度太大,云图不连续,那是网格问题,不是平均阈值问题。
显式动力学,接触前的小应力,大多是数值噪声,结合数量级和单位制判断。
壳单元S33,理论上应该是0,显示非零大值,检查载荷方向或单元公式。
主应力排序,Abaqus按代数值排序,拉为正,压为负,不是按绝对值大小。
单位一致性,应力单位看输入模型单位,别忽略了数量级。
提取应变数据,看需求,全局、特定节点/单元集、批量导出。
插件快速提取,推荐用,功能强大,操作简单。
Python脚本定制,内置模块 odbAccess,灵活度高。
可视化模块手动导出,少量数据或临时导出。
应变云图,Visualization模块,Field Output选应变变量,Plot Contours on Deformed Shape,搞定。
注意应变类型,Nlgeom=ON,显示LE(对数应变),否则E(名义应变)。
视图异常,断层或颜色单一,检查未变形形状或数值范围。
数据导出,提取具体数值,Query > Probe values,或File > Write Input。
提取应变数据,前置输出设置E/PE,后处理Report、XY Data插值或Python脚本导出。
确保计算时已输出应变,分析步模块场输出请求里得有应变变量。
CAE图形界面导出,Report > Field Output,设置文件名、小数精度。
沿路径提取分布数据,Tools > Path > Create,生成XY数据,导出为.txt。
Python脚本批量自动化,效率高,可定制。
积分点vs节点,应变默认积分点,节点值是外推平均。
数据量控制,全模型节点导出文件大,筛选关键区域。
变量含义,E为总对数应变,PE为塑性应变,PEEQ为等效塑性应变。