接触面节点力提取及再施加方法心得分享

这几天做了个，需要将接触面上的力提取做被接触部件的进一步分析，而ABAQUS中没有自动提取施加接触力的功能，所以自己动手用matlab+inp的方式，实现了一个很原始的接触力提取和施加。可能用户自定义载荷（VDLoad）能够实现，但是时间不够没有尝试这方面的东西。

首先说一下，abaqus中的一个不是很完整的功能。

就是Analytical Field - Mapped field，Mapped field 是个很有意思的功能，它包含了两种一个是Point cloud field，即用一组对应空间点的载荷，来表达载荷的空间分布，然后通过插值得到节点上载荷。

第二个叫做.odb mapped field（看到odb你懂的吧），这个就是自动提取odb文件中的载荷分布，这个的用法很简单，在viewport菜单中新建一个viewport，叫做viewport2然后在这个viewport2里打开载荷存储的odb文件，并打开你要提取的载荷的Field output，然后在viewport1中进入到.odb mapped field，有按钮可选提取哪个viewport的分布。很简单，很智能。

但是，还是要但是，在纯力学的范畴内，这两个功能都只能用在dload就是，pressure里，其余的载荷类型均无法识别。（除非只要接触正压强，切向的力完全弄不出来）

所以另辟道路，简单的说吧，只提个方法好了！就是用编程语言改动后处理输出的接触面节点力的prt文件成为inp的格式，然后放入对应代码模块即可！！！

很简单的思路，但是有一些注意的地方：

1. 节点对应，即网格划分得相同；所以最好copy一个part，然后换成mesh part；

2. 注意提取接触力时别弄错了，接触力有方向的问题，自己注意选择；

3. 不要忘记了，节点还有contact shear force，如果照此方法，应该是 press force + shear force。