一种基于abaqus的真实粗糙表面建模方法

摘要

本发明公开了一种基于abaqus的真实粗糙表面建模方法，适用于三维真实粗糙表面在abaqus软件中的快速建模。本方法通过三维形貌仪获取真实粗糙表面的三维形貌数据，采用matlab进行数据降噪处理，在abaqus中实现了三维实体模型的建立。采用python语言编写程序，生成点、线、面，生成真实粗糙表面壳体后并进行实体化，最终得到真实粗糙表面的三维实体模型。根据数据直接在有限元分析软件abaqus中进行建模，可以有效避免将数据导入其他软件建好模型之后再导入abaqus中繁琐步骤以及可能会存在的模型失真的情况，并且在很大程度上提升三维模型的准确性，也为后续进一步实现真实粗糙表面接触特性的自动化分析打下良好的基础。

说明书

技术领域

本发明涉及真实粗糙表面的精确建模技术领域，尤其涉及一种基于abaqus二次开发真实粗糙表面的建模方法。

背景技术

实际工程中相互接触的零件表面，从宏观上看是非常光滑的，但是在接触表面的微观尺度上可以看到许多的凸峰和凹谷。当两个表面相互接触时，实际上发生接触的是两个表面上的微凸体，真实的接触面积只是名义接触面积的一部分。因此，真实粗糙表面的建模对粗糙表面间的接触特性分析有着重要意义。

目前对于使用表面测量仪数字化真实粗糙表面得到的数据，通常是使用三维建模软件进行建模，如PREO,SOLIDWORKS,UG，在这些三维软件中建立模型之后然后导入仿真软件中进行分析。但是，首先在三维软件中进行真实粗糙表面建模的精度不是很高，其次将建立好的三维模型导入有限元软件时还会存在模型数据丢失，失真等情况，而这些对于粗糙表面进行仿真分析的准确性是非常不利的。

针对以上不足，有必要基于abaqus二次开发对真实粗糙表面进行精确的建模，实现对真实粗糙表面建模和有限元分析都在abaqus中进行的目的。有效避免了将大量数据倒来倒去的麻烦，同时不用考虑将模型导入abaqus软件时可能会存在数据丢失的问题。

发明内容

本发明的目的在于提供一种基于abaqus二次开发粗糙表面的建模方法，利用abaqus有限元软件建立精确的真实粗糙表面模型。

为了实现上述目的，本发明采用如下技术方案：

步骤一，真实粗糙表面微观形貌数据的获取。

使用三维形貌测量仪获取粗糙表面的三维形貌数据，通过三维表面形貌仪获取的点云数据为matlab软件可直接进行提取的txt格式，其中第一列，第二列，第三列数据分别对应为每个采集点的x，y，z坐标。

步骤二，点云数据去噪处理。

在matlab软件中采用小波降噪的方法对点云数据进行处理，以达到数据滤波去噪的目的。将处理后的数据保存为有限元软件abaqus可以直接进行提取的txt格式。

步骤三，根据去噪后的点云数据，使用python语言编写程序建立真实粗糙表面。具体包含以下步骤：

1)创建粗糙表面的一半模型part-1。

A.创建基准点。

使用python语言进行二次开发将预处理后的三维形貌数据导入，使用while循环在abaqus中生成基准点。

B.由点生线，由线生面。

导入abaqus软件中生成每个点云数据对应坐标的基准点，是一种有序分布的点云数据形式，它的点云形状类似于矩阵，基准点的分布十分规律。首先，从第一个点开始，使得相邻的三个点生成线，使用这三条线生成一个三角面。然后，使用while循环命令，直到生成真实粗糙表面的一半三角化表面。

2)创建粗糙表面的另一半模型part-2。

A.同上面1)中的A步骤一样将含有去噪后点云数据的txt文件导入最终生成基准点。

B.由点生线，由线生面。

根据生成的基准点，首先，从最后一个点开始，使得相邻的三个点生成线，使用这三条线生成一个三角面。然后，使用while循环命令，直到生成真实粗糙表面的另外一半三角化表面。

3)将part-1和part-2进行装配，并生成粗糙表面。

编写程序进入abaqus装配模块，将part-1和part-2两个部件导入完成装配，生成完整的真实粗糙表面part-3，并对装配后的part-3进行光顺处理。

4)创建粗糙表面模型的其余5个面。

通过创建其余五个平面，最终形成真实粗糙表面三维模型的壳体。

5)通过面生成体。

使用通过壳创建实体功能，框选中所有的壳生成实体，完成三维真实粗糙表面的创建。

本发明提出一种基于abaqus二次开发粗糙表面的建模方法。利用三维形貌仪采集到包含表面波纹度、表面形状误差、表面粗糙度等成分的三维表面形貌数据，将数据导入matlab进行数据预处理，最终基于abaqus二次开发建立三维真实粗糙表面模型。本发明减少了大量的数据倒来倒去的步骤，避免了三维模型导入abaqus中可能会出现失真的情况，有效保证了真实粗糙表面接触分析模型的准确性。并且后续要批量建立三维真实粗糙表面模型，只需要修改程序中的某些数值即可，避免了繁琐的操作步骤。

附图说明

图1三维形貌仪采集的0.38mmx0.35mm试件表面形貌。

图2滤波后将数据导入abaqus中生成基准点

图3生成三维粗糙表面的半模型part-1

图4生成三维粗糙表面的另外一半模型part-2

图5将part-1和part-2装配后生成的表面part-3

图6对part-3做光顺处理

图7生成真实粗糙表面模型壳

图8实体化后生成的真实粗糙表面

具体实施方式

以下结合图1至图8所示，及具体实例对本发明的具体实施步骤进行说明

1、真实粗糙表面三维形貌数据的获取

使用三维表面形貌仪获取三维形貌数据具体步骤如下：

首先，将待测工件放在酒精溶液中使用超声清洗10min，清洗完成后进行吹干，保证表面无杂质。其次，选择合适量程的光学测量探头，进行采样区域大小和采样步长的设定，并对待测加工试样进行数据测量。本例中，X轴方向采样距离为0.38mm；Y轴方向采样距离为0.35mm；X方向采样间隔为20um，Y方向采样间隔为50um。所获得的数据如图1所示。

2、将采集到的三维形貌数据导入matlab中进行降噪处理。对数据进行预处理之后，将数据保存为abaqus软件可以直接进行提取的txt格式。

3、创建真实粗糙表面的一半三角化表面part-1具体步骤如下：

A.将含有去噪后点云数据的txt文件，使用open命令以只读的方式打开。

B.txt文件中一行数据即为一个点的坐标，使用readlines命令按每一行读取txt文件中的数据保存在content变量中，content变量以列表的形式保存了每一个坐标点，其中每一个坐标点都是字符串的形式。

C.字符串中含有x，y，z坐标，且含有空格和回车符。接下来使用for循环遍历content中每一个含有坐标数据的字符串，对获得的每一个字符串使用split命令进行分割去除空格和回车符，得到的每一个坐标点都以列表的形式存储，x，y，z坐标都是字符串形式。

D.最后，将每一个列表中字符串形式的x，y，z点转换为对应的浮点型数据，并生成基准点，如图2所示。

E.从第一个点开始，将点的标号为i(i是除了每一行最后一个点及最后一行点的标号)与i+1，i与i+a(a为x坐标方向每一行点的个数)，i+1与i+a连接成线，然后生成一个以这三条线为边界的面，使用while循环一直到生成粗糙表面的一半模型part-1，如图3所示。

4、创建真实粗糙表面的另一半三角化表面part-2具体步骤如下：

A.前几个步骤与步骤3中的A、B、C、D相同。

B.生成基准点之后，从最后一点开始，将点的标号为j(j是除了每一行第一个点及第一行点的标号)与j-1，j与j-a，j-1与j-a连接成线，然后生成一个以这三条线为边界的面，使用while循环一直到生成粗糙表面的另一半模型part-2，如图4所示。

5、将part-1和part-2进行装配，并生成粗糙表面，具体步骤如下：

A.进入abaqus装配模块中，将part-1和part-2导入，自动完成装配，如图5所示。

B.使用装配模块中的合并功能，选中两个几何实例进行合并，合并后生成part-3。

C.然后进入到部件模块，使用部件模块中的修复面功能，选中part-3中所有的面，得到一个去除三角边界的光顺的三维真实粗糙表面，如图6所示。

6、创建真实粗糙表面三维模型，具体步骤如下：

A.创建四个基准点，这四个基准点的坐标分别是第一行和最后一行点的坐标并将它们的z坐标改为c值(c为粗糙面模型底面在z方向上的坐标)。

B.将这四个基准点与粗糙面上的基准点相互之间连接成线。

C.对生成其它5个面的线别创建五个线的集合，set-1,set-2,set-3,set-4,set-5。设置一个集合生成一个面后再设置下一个集合，最终生成五个面，如图7所示。

D.使用通过壳创建实体功能，选中生成的六个面，最终生成真实粗糙表面三维模型，如图8所示。