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

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

Info

Publication number
CN112613218A
CN112613218A CN202011569705.7A CN202011569705A CN112613218A CN 112613218 A CN112613218 A CN 112613218A CN 202011569705 A CN202011569705 A CN 202011569705A CN 112613218 A CN112613218 A CN 112613218A
Authority
CN
China
Prior art keywords
rough surface
data
abaqus
point
dimensional
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011569705.7A
Other languages
English (en)
Other versions
CN112613218B (zh
Inventor
赵永胜
刘盼
刘志峰
李迎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202011569705.7A priority Critical patent/CN112613218B/zh
Publication of CN112613218A publication Critical patent/CN112613218A/zh
Application granted granted Critical
Publication of CN112613218B publication Critical patent/CN112613218B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Processing Or Creating Images (AREA)

Abstract

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

Description

一种基于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所示。

Claims (2)

1.一种基于abaqus真实粗糙表面的建模方法，其特征在于，包含下述步骤：
1)使用三维形貌测量仪获取粗糙表面的三维形貌数据，通过三维表面形貌仪获取的点云数据为matlab软件直接进行提取的txt格式，其中第一列，第二列，第三列数据分别对应为每个采集点的x，y，z坐标；
2)根据去噪后的点云数据，在abaqus软件中基于python语言进行二次开发建立真实粗糙表面模型。
2.根据权利要求1所述的一种基于python语言进行二次开发的方法，其特征在于，具体包含下述步骤：
1)创建真实粗糙表面的一半三角化表面part-1，具体如下：
将含有去噪后点云数据的txt文件，使用open命令以只读的方式打开；
txt文件中一行数据即为一个点的坐标，使用readlines命令按每一行读取txt文件中的数据保存在命名为content的变量中，content变量以列表的形式保存了每一个坐标点，其中每一个坐标点都是字符串的形式；
字符串中含有x，y，z坐标，且含有空格和回车符；接下来使用for循环遍历content中每一个含有坐标数据的字符串，对获得的每一个字符串使用split命令进行分割去除空格和回车符，得到的每一个坐标点都以列表的形式存储，x，y，z坐标都是字符串形式；
最后，将每一个列表中字符串形式的x，y，z点转换为对应的浮点型数据，并生成基准点；
从第一个点开始，将点的标号为i与i+1，i与i+a，i+1与i+a连接成线，然后生成一个以这三条线为边界的面，使用while循环一直到生成粗糙表面的一半模型part-1；i是除了每一行最后一个点及最后一行点的标号；a为x坐标方向每一行点的个数；
2)创建真实粗糙表面的另一半三角化表面part-2，具体如下：
生成基准点之后，从最后一点开始，将点的标号为j与j-1，j与j-a，j-1与j-a连接成线，然后生成一个以这三条线为边界的面，使用while循环一直到生成粗糙表面的另一半模型part-2；j是除了每一行第一个点及第一行点的标号；
3)将part-1和part-2进行装配，并生成粗糙表面，具体如下：
进入abaqus装配模块中，将part-1和part-2导入，自动完成装配；
使用装配模块中的合并功能，选中两个几何实例进行合并，合并后生成part-3；
然后进入到部件模块，使用部件模块中的修复面功能，选中part-3中所有的面，得到一个去除三角边界的光顺的三维真实粗糙表面；
4)创建真实粗糙表面三维模型，具体如下：
创建四个基准点，这四个基准点的坐标分别是第一行和最后一行两个端点的坐标并将它们的z坐标改为c值；c为粗糙面模型底面在z方向上的坐标；
将这四个基准点与粗糙面上的基准点相互之间连接成线；
对要生成其它5个面的线分别创建五个线的集合set-1,set-2,set-3,set-4,set-5；设置一个集合生成一个面之后再设置下一个集合，最终生成五个面；
使用通过壳创建实体功能，选中生成的六个面，最终生成真实粗糙表面三维模型。
CN202011569705.7A 2020-12-26 2020-12-26 一种基于abaqus的真实粗糙表面建模方法 Active CN112613218B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011569705.7A CN112613218B (zh) 2020-12-26 2020-12-26 一种基于abaqus的真实粗糙表面建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011569705.7A CN112613218B (zh) 2020-12-26 2020-12-26 一种基于abaqus的真实粗糙表面建模方法

Publications (2)

Publication Number Publication Date
CN112613218A true CN112613218A (zh) 2021-04-06
CN112613218B CN112613218B (zh) 2024-03-15

Family

ID=75248148

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011569705.7A Active CN112613218B (zh) 2020-12-26 2020-12-26 一种基于abaqus的真实粗糙表面建模方法

Country Status (1)

Country Link
CN (1) CN112613218B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102254354A (zh) * 2011-06-27 2011-11-23 中国科学院武汉岩土力学研究所 一种实现土的微观结构表面形态三维可视化的方法
CN106096149A (zh) * 2016-06-15 2016-11-09 东南大学 一种微观路表形貌的数字化生成方法
CN108168492A (zh) * 2017-12-17 2018-06-15 北京工业大学 一种基于三维微观表面形貌的接触特性分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102254354A (zh) * 2011-06-27 2011-11-23 中国科学院武汉岩土力学研究所 一种实现土的微观结构表面形态三维可视化的方法
CN106096149A (zh) * 2016-06-15 2016-11-09 东南大学 一种微观路表形貌的数字化生成方法
CN108168492A (zh) * 2017-12-17 2018-06-15 北京工业大学 一种基于三维微观表面形貌的接触特性分析方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
G YI: "Application of python-based Abaqus preprocess and postprocess technique in analysis of gearbox vibration and noise reduction", FRONTIERS OF MECHANICAL ENGINEERING *
WEIXIN_39743423: "python显示点云_python处理点云数据并生成三维点云模型", Retrieved from the Internet <URL:https://blog.csdn.net/weixin_39743423/article/details/110322162> *
毕舒心等: "基于粗糙齿面的渐开线齿轮齿面接触刚度建模", 现代制造技术与装备 *
陈庆发等: "大型复杂地质体三维数值模型构建方法比较研究", 岩土力学 *
魏兵辉等: "高速干滑动摩擦系数的有限元仿真计算研究", 西安交通大学学报, pages 2 *

Also Published As

Publication number Publication date
CN112613218B (zh) 2024-03-15

Similar Documents

Publication Publication Date Title
CN101067872B (zh) 从三维逆向建模的源模型数据检测二维草图数据的方法
CN111077844B (zh) 一种基于实测数据特征引导的零件精确加工方法
CN107045526A (zh) 一种电子建筑施工图的图形识别方法
CN109146711A (zh) 一种综合利用bim技术和3d激光扫描系统优化施工的方法
CN103488838B (zh) 一种余热锅炉管道计算机辅助出图方法
CN109859317B (zh) 一种基于cass和catia的3dgis地形模型快速建模方法
CN103678737B (zh) 一种杠杆系统虚拟装配方法
CN116402866A (zh) 基于点云的零件数字孪生几何建模与误差评定方法及系统
CN101281403A (zh) 异型玻璃加工图形定位方法及系统
CN109359333B (zh) 一种包含多尺度形貌特征的体模型构建方法
CN110795835A (zh) 一种基于自动同步建模的三维工序模型逆向生成方法
CN112506474A (zh) 一种基于mbd模型的3d测头在机测量自动编程方法
CN114608461A (zh) 一种非均匀壁厚零件的激光扫描测量方法
CN117057064A (zh) 一种机械设备零部件参数化快速出图的方法
CN105045958A (zh) 基于bp神经网络的gps高程拟合的实现系统及方法
CN109859301A (zh) 一种岩石结构面粗糙度系数精细化表征方法
CN112613218A (zh) 一种基于abaqus的真实粗糙表面建模方法
CN117392200A (zh) 一种基于点云融合的车轮尺寸测量系统及方法
CN115641405A (zh) 一种三维建模方法及系统
CN116255930A (zh) 一种基于点云切片的截面提取及测量方法及系统
CN112990373B (zh) 基于多尺度特征融合的卷积孪生点网络叶片轮廓拼接系统
CN115147471A (zh) 一种基于曲率密度特征的激光点云自动配准方法
CN110555245B (zh) 一种三维尺度精准关注部位提取应力梯度的方法
CN111590265A (zh) 基于激光扫描技术的模具修复方法
CN111260783B (zh) 一种基于K近邻和Poisson曲面的矿体三维自动建模方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant