近期一直在各种忙, 公 众 号挺久未更新了,借一学员问的问题,写一篇技文:“这种图怎么出呢?”其中的图如下,裂缝宽度在半缝高/半缝长上的变化,通过输出此图可以非常直观的观察:不同工况下,参数对结果的影响。
下面我们就来详细讲解这样曲线的提取过程,以案例分享的案例7:二维单缝压裂Cohesive方法为例,提取不同时刻下的裂缝形态。
步骤如下:
1)公 众 号回复“二维单缝压裂Cohesive方法”,可获得inp文件;
2)按照说明文件提交inp进行计算,获得odb结果文件;
3)在ABAQUS后处理中打开odb文件,变形放大100倍结果如下:
4)创建Path路径,由于我们的这个问题比较特殊,采用的零厚度Cohesive单元模拟水力裂缝,因此Path顶点选取需要格外小心,将裂尖位置放大至足够区分Cohesive单元中间渗流节点和边节点,只有边节点才有位移自由度。
菜单栏Tools -> Path -> Create…,在弹出的对话框中选择Node list,并任意点击Add Befor或Add After按钮。
从图中选择裂缝一侧的所有边节点,如下图所示,随后点击Done,并在弹出对话框中点击OK,Path路径创建完成。
5)云图调整为U1,即X轴(水平)位移云图,Cohesive单元边节点的水平位移量即为裂缝的半缝宽。
6)菜单栏Tools –> XY Data -> Create… 在弹出对话框中选择Path类型。
后续对话框定义如下,选择刚刚定义的路径,并选择输出的Frame时刻和输出的场量为U1,点击Apply可绘制出整个缝长上半缝宽的变化曲线。
7)从软件左侧的结果树中,找到绘制的XYdata数据名称,右击选择Edit,即可复 制并粘贴到Excel中进一步处理。
8)不同时刻的曲线,只需调整第6步中的Step/Frame时刻即可完成。
9)有这些数据,通过EXCEL绘图,添加储隔层分界线和文字注释,这些都是很简单的操作了。
该技巧还可以拓展应用到非常多的情况下,比如提取梁柱挠度曲线、获得指定路径下的场量变化等,灵活应用ABAQUS提供的后处理功能,已经能满足大家绝大多数需求,当然如果觉得这个过程繁琐,还可以写一小段Python程序,从而一劳永逸。