abaqus失效单元删除_情人节？小同人教你如何应景的使用Abaqus/Explicit

    每次写文案时总是第一段话最难，想以客套的寒暄开场，却怕显得生疏；想以熟悉的问候致意，却又彷徨言辞稚嫩。所以不同人的生活有着不同的2月14日，不同人的眼中也有不同的小同人。或许最好的代入，莫过开门见山。

   那么，在还春暖花开的将来未来之际，不妨用Abaqus/Explict看一看如何的“怦然心动“与“一箭钟情”。

怦然心动

心动过度

一箭钟情的外表

一箭钟情的内在

怦然心动

    总的来说，怦然心动的实现，需要彼此认识(接触设置)、缘分使然(初始速度)以及缘分到位(显式分析)。

    心形模型的建模其实不难，先从CAD中绘制好样条曲线，导入abaqus中成为sketch，再拉伸成实体，进而对实体edge倒角，最后通过抽离形成shell即可。

心形建模的基本设置

    比较关键的一步在于设置材性参数。照顾广大单身狗今日份的玻璃心，便选用了适用于陶瓷的脆性材料，材料设置基本如下：

脆性材料的基本设置

    为了营造遇见对的人之后刹那间的惊鸿一瞥，便定义了一个又冷又硬的shell part作为碰撞的边界。由于shell会被约束成刚体，材性可以随便定义，注意网格划分的时候，尺寸因略大于心形体的网格尺寸，并需要选用explict单元。刚性边界的实现通过constraint中的rigid body模拟，将shell绑在预先定义好的reference point上即可，如下。

刚性边界的定义

    剩下的就是显式分析荷载步与初始速度的施加，较为容易，就不多做赘述乐。但是这里需要注意的是，为了模拟材料失效，必须在mesh中赋予element delete。同时在field output中，必须勾选status。如果不勾选status，odb的显示效果将相当劲爆，有兴趣的朋友可以自行尝试。

单元删除的定义

输出设置

    这些都弄好之后，怦然心动就可以提交分析了。如果两个人缘分正好(心形part初始速度设置为10m/s)，或许正如下图中小鹿乱撞：

    可是如果缘分没到位(心形part初始速度设置为60m/s)，那或许就是火星撞地球刹不住车，只剩满地残骸了：

一箭钟情

    “一箭钟情”其本质上同源于军工领域的侵彻分析。需要分别定义一个靶体和一个侵彻体。故而这个模型里不再使用shell单元而改用solid单元。

    心形体建模同前，只是不再抽离成壳。箭头(侵彻体)的建模采用3D solid revolve。如下图所示。

箭头几何模型

   同理，箭头通过和预先定义的reference point进行rigid body约束，形成刚性体，与心形呈45度角贯入。如下图。

箭头与心形的assembly

    材料仍然采用适用于玻璃心的陶瓷脆性材料，对心形part定义了element delete，并且在field output中勾选了status选项。同上，此处不再赘述。

    但是，这个分析中最为关键的一个细节在于接触参数的设置。因为箭头贯入心形solid，导致表面单元失效后，其内部单元将重新与箭头刚性体接触。如果直接采用general contact中的all with，势必造成单元失效之后的接触失效，虽然不会导致求解失败，但此时箭头已经无法继续向心形体的内部单元传递压力。

    为了解决这个bug，在定义接触之前，预先定义两个基于mesh的surface。一个surface选中箭头所有内部单元的表面，另一个surface选中心形体中所有内部单元的表面。继而，接触设置如下：

    可以说这是保障侵彻分析最为关键、最容易忽视、也最重要的一步，完成之后，就是简单的explicit荷载步定义与初始条件的施加，本文不多做赘述。

    此后，提交job分析，即获得了缘分来临时的沦陷与沉浸(右图是定义了view cut后的内部应力云图)：

一箭钟情

    本文谨以Abaqus/Explict的两个非线性算例，介绍了单元删除、脆性破坏、材料失效、通用接触与显式分析的用法与关键细节(打赏15元并留言邮箱可获得CAE分析模型与inp文件哦)。在这个普天同庆(！！！/？？？)的日子里，小同人衷心祝愿缘分已经来临的朋友天长地久，祝愿缘分还在路上的朋友早日收获幸福。

南乡子.小同人

    掩卷且翻篇。一任繁华艳阳天，新绿渐稠花尽染。开颜。扶摇万里是炊烟。

    往来尽酒酣。极目中流共江山，秦淮烟波为谁散。樽前，应许策马五百年。

——小同人随笔附赠