CN106844912A - 一种基于网格变形技术的机构运动副磨损仿真方法 - Google Patents
一种基于网格变形技术的机构运动副磨损仿真方法 Download PDFInfo
- Publication number
- CN106844912A CN106844912A CN201710010871.5A CN201710010871A CN106844912A CN 106844912 A CN106844912 A CN 106844912A CN 201710010871 A CN201710010871 A CN 201710010871A CN 106844912 A CN106844912 A CN 106844912A
- Authority
- CN
- China
- Prior art keywords
- model
- mechanism kinematic
- grid node
- bead
- target area
- 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
Links
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于网格变形技术的机构运动副磨损仿真方法,包括以下步骤:创立机构运动副模型,判断机构运动副模型的接触区域,记为目标区域,采用网格划分目标区域,并输出网格节点的编号和原始坐标;在网格节点处创建无质量的小球模型,并用小球模型代替目标区域与机构运动副模型接触;执行动力学仿真,获取小球模型与机构运动副模型间的相互作用力及相对滑动速度;根据阿查德磨损计算公式计算网格节点处的磨损深度;根据磨损深度计算出网格节点的新坐标。通过划分网格的方式将目标区域离散,建立起接触压力、运动副元素间相对滑动速度与磨损量之间的定量关系,并通过可视化展现目标区域磨损后特征,为结构的设计提供了参考。
Description
技术领域
本发明涉及机械结构仿真设计技术领域,具体来说,涉及一种基于网格变形技术的机构运动副磨损仿真方法。
背景技术
机械产品的一个主要失效形式就是零部件表面磨损。绝大多数零件都与别的零件发生静或动的接触和配合关系,载荷作用于表面,摩擦和磨损就发生了,磨损积累到一定程度,零件的尺寸变化超出允许范围造成零件报废。零件的使用寿命在很大程度上收到表面磨损的限制。
零件表面的磨损机理十分复杂,不仅与零件本身的材料性质有关,而且与零件之间的相互作用有着密切联系,目前对于零部件表面磨损的理论研究和试验研究都有着较高的技术难度。网格变形法是一种从应变量换算应力分布的构造模拟法。先在模型表面划出规则的小网格(网格可为方形、矩形或者三角形),然后对模型施加载荷,用读数显微镜度量加载前后网格坐标之差求出线应变或角应变,根据应力应变关系求出应力分布。网格变形技术是数字几何处理的一个核心技术,近年来得到了国内外研究者的极大关注,且在逆向工程、工业品创新设计及计算机动画等领域得到了广泛的应用。
1953年阿查德提出了磨损体积与载荷、滑动距离成正比、接触面中较软金属的硬度成反比的经验磨损表达式。阿查德磨损方程为:
式中:V为磨损掉的材料体积,S接触面相对滑动距离,F是接触面间法向载荷,H是较软材料布氏硬度值,K是产生磨粒的概率即磨损系数。阿查德磨损方程左右两边同除以面积A可得如下表达式:
又dS=vdt,用表示磨损深度,那么式(2)可以表示为:
δ=∫kpvdt (3)
式中,是与材料有关的磨损系数,v表示接触面相对滑动速度,是单位面积接触压力。
发明内容
本发明的目的在于提出一种基于网格变形技术的机构运动副磨损仿真方法,能够实现对零件表面磨损情况的仿真。
为实现上述发明目的,本发明的技术方案是这样实现的:
一种基于网格变形技术的机构运动副磨损仿真方法,包括以下步骤:
1)创立机构运动副模型,判断机构运动副模型的接触区域,记为目标区域,采用网格划分所述目标区域,并输出网格节点的编号和原始坐标;
2)在所述网格节点处创建无质量的小球模型,并用小球模型代替目标区域与机构运动副模型接触,建立起小球模型与机构运动副模型的动力学仿真模型;
3)执行动力学仿真,获取小球模型与机构运动副模型间的相互作用力及相对滑动速度;
4)根据阿查德磨损计算公式计算网格节点处的磨损深度,所述阿查德磨损计算公式如下所示:
δ=∫kpvdt
其中:δ表示磨损深度,k为磨损系数,p为压强,v是小球模型与机构运动副接触面间的相对滑动速度;
5)根据所述磨损深度计算出网格节点的新坐标。
进一步的,步骤1)中采用网格划分软件将目标区域划分为由若干矩形单元和三角形单元组成的混合网格,并将网格节点的编号和原始坐标输出并保存于如下的矩阵A中,
其中:矩阵A的第一列为网格节点的编号,其余三列为网格节点的原始坐标。
进一步的,所述机构运动副模型包括支架模型、电磁铁模型、两个抱爪模型、两个连接臂模型和导向块模型,两个所述抱爪模型通过销轴模型对称设于支架模型的横向两端,所述导向块模型包括中间的导架模型,所述导架模型的横向尺寸沿纵向向下逐渐增大,所述导架模型上沿纵向设有滑道模型,所述滑道模型内设有挡块模型,所述挡块模型固定设于支架模型上,两个所述连接臂模型之间连接有拉簧模型,两个所述连接臂模型的两端在拉簧模型的作用下分别与导架模型的横向两端和抱爪模型靠近销轴模型的一端相抵接,所述电磁铁模型设于导架模型的下面,所述导架模型在电磁铁模型的作用下能够在所述滑道模型的行程内沿纵向向下滑动,所述导架模型还连接有用于使导架模型向上滑动的复位弹簧模型。
进一步的,步骤5)中按照如下公式计算出网格节点的新坐标
其中:Δh=δ表示磨损深度,表示磨损方向,A′表示磨损后网格节点的编号坐标矩阵。
进一步的,步骤2)中所述小球模型的半径小于0.1mm,且其球心与网格节点相重合。
进一步的,步骤2)中采用matlab编写并输出ADAMS的命令流CMD文件的方式来设置小球模型的分布,其中CMD文件的作用为在DAMS工作环境中,通过命令在网格节点处创建无质量小球模型,并在小球模型和导向块模型之间创建接触。
本发明的有益效果:本发明通过将网格变形技术应用于机构运动副的仿真,具体来说,通过划分网格的方式将目标区域离散,建立起接触压力、运动副元素间相对滑动速度与磨损量之间的定量关系,并通过可视化展现目标区域磨损后特征;本发明基于网格变形技术的运动副磨损仿真法为机构的结构设计、寿命设计提供了参考。
附图说明
图1是本发明所述的抱爪机构模型的仿真结构示意图;
图2是本发明所述的抱爪机构模型的局部仿真结构示意图;
图3是本发明所述的磨损仿真的流程示意图;
图4是本发明所述的抱爪模型与连接臂模型作用区域的小球模型分布示意图;
图5是本发明所述的抱爪模型与连接臂模型作用区域的磨损结构示意图;
图6是本发明所述的抱爪模型与连接臂模型作用区域的磨损量投影图。
图中所示:
1-支架模型;2-电磁铁模型;3-抱爪模型;4-连接臂模型;5-导向块模型;6-导架模型;7-滑道模型;8-挡块模型;9-拉簧模型;10-销轴模型。
具体实施方式
下面结合本发明的附图,对本发明的技术方案进行清楚、完整地描述。
如图3所示,本发明实施例所述的一种基于网格变形技术的机构运动副磨损仿真方法,包括以下步骤:
1)创立机构运动副模型,判断机构运动副模型的接触区域,记为目标区域,采用网格划分所述目标区域,并输出网格节点的编号和原始坐标;
在本实施例中,如图1-2所示,所述机构运动副模型包括支架模型1、电磁铁模型2、两个抱爪模型3、两个连接臂模型4和导向块模型5,两个所述抱爪模型3通过销轴模型10对称设于支架模型1的横向两端,所述导向块模型5包括中间的导架模型6,所述导架模型6的横向尺寸沿纵向向下逐渐增大,所述导架模型6上沿纵向设有滑道模型7,所述滑道模型7内设有挡块模型8,所述挡块模型8固定设于支架模型1上,两个所述连接臂模型4之间连接有拉簧模型9,两个所述连接臂模型4的两端在拉簧模型9的作用下分别与导架模型6和抱爪模型3相抵接,所述电磁铁模型2设于导架模型6的下面,所述导架模型6在电磁铁模型2的作用下能够在所述滑道模型7的行程内沿纵向向下滑动,所述导架模型6还连接有用于使导架模型6向上滑动的复位弹簧模型(图中未示)。
在本实施例中,步骤1)中采用网格划分软件将目标区域划分为由若干矩形单元和三角形单元组成的混合网格,并将网格节点的编号和原始坐标输出并保存于如下的矩阵A中,
其中:矩阵A的第一列为网格节点的编号,其余三列为网格节点的原始坐标。
上述矩形单元或三角形单元是平直不可弯曲的区域面,其数量和尺寸决定了网格表面的几何精度和磨损表面的平滑程度。
2)在所述网格节点处创建无质量的小球模型,并用小球模型代替目标区域与机构运动副接触,建立起小球模型与机构运动副的动力学仿真模型;
在本实施例中,采用matlab编写并输出ADAMS的命令流CMD文件的方式来设置小球模型的分布,其中CMD文件的作用为在DAMS工作环境中,通过命令在网格节点处创建无质量小球模型,并在小球模型和连接臂模型之间创建接触。在小球模型的创建过程中,要求所述小球模型的半径小于0.1mm,且其球心与网格节点相重合。在网格节点处创建无质量小球模型目的是用离散的小球模型代替目标表面与导向块模型接触,理论上小球模型数目越多越接近实际情况,本实施例共创建了964个小球,参考图4显示了创建后小球的分布。
3)执行动力学仿真,获取小球模型与机构运动副间的相互作用力及相对滑动速度;
抱爪模型3与连接臂模型4属于点面的接触,抱爪模型3与连接臂模型4接触部位(即目标区域)容易磨损,当磨损量累积到一定程度,抱爪模型3就难以抱紧物料而失效。考虑到两个抱爪模型3的对称性,取左抱爪模型作为目标进行磨损仿真分析,右抱爪模型将具有与其相似的磨损特征。
具体操作时,在ADAMS工作环境中,导架模型6在电磁铁模型2和复位弹簧模型的作用下上下移动,从而控制抱爪模型3的张合。抱爪模型3在张开和闭合过程中目标区域表面受到挤压和摩擦,记录目标区域表面的接触信息。执行动力学仿真时,运动副元素间的碰撞摩擦过程划分为若干时间步长,记录每个时间步长下目标区域表面的接触信息,接触信息包括运动副元素间接触点位置、接触力、接触过程中的速度位移。
4)根据阿查德磨损计算公式计算网格节点处的磨损深度,所述阿查德磨损计算公式如下所示:
δ=∫kpvdt
其中:δ表示磨损深度,k为磨损系数,p为压强,v是小球模型与机构运动副接触面间的相对滑动速度;
接触面相对滑动速度v通过步骤3)的动力学仿真得来;单位面积接触压力其中FN是法向接触压力,式中Si表示节点所在单元面积,三角形单元时a=3,四边形单元时a=4;磨损系数k是根据接触面间材料属性和润滑条件查表得到。
5)根据所述磨损深度计算出网格节点的新坐标。如此得到了磨损后目标表面的磨损特征,如下循环实现了机构运动副磨损的动态仿真。图5和图6显示了目标表面的磨损特征.
在本实施例中,步骤5)中按照如下公式计算出网格节点的新坐标
其中:Δh=δ表示磨损深度,表示磨损方向,即为运动副接触面上网格节点位置的法向向量,A′表示磨损后网格节点的编号坐标矩阵。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于网格变形技术的机构运动副磨损仿真方法,其特征在于,包括以下步骤:
1)创立机构运动副模型,判断机构运动副模型的接触区域,记为目标区域,采用网格划分所述目标区域,并输出网格节点的编号和原始坐标;
2)在所述网格节点处创建无质量的小球模型,并用小球模型代替目标区域与机构运动副模型接触,建立起小球模型与机构运动副模型的动力学仿真模型;
3)执行动力学仿真,获取小球模型与机构运动副模型间的相互作用力及相对滑动速度;
4)根据阿查德磨损计算公式计算网格节点处的磨损深度,所述阿查德磨损计算公式如下所示:
δ=∫kpvdt
其中:δ表示磨损深度,k为磨损系数,p为压强,v是小球模型与机构运动副接触面间的相对滑动速度;
5)根据所述磨损深度计算出网格节点的新坐标。
2.根据权利要求1所述的仿真方法,其特征在于,步骤1)中采用网格划分软件将目标区域划分为由若干矩形单元和三角形单元组成的混合网格,并将网格节点的编号和原始坐标输出并保存于如下的矩阵A中,
其中:矩阵A的第一列为网格节点的编号,其余三列为网格节点的原始坐标。
3.根据权利要求1所述的仿真方法,其特征在于,所述机构运动副模型包括支架模型、电磁铁模型、两个抱爪模型、两个连接臂模型和导向块模型,两个所述抱爪模型通过销轴模型对称设于支架模型的横向两端,所述导向块模型包括中间的导架模型,所述导架模型的横向尺寸沿纵向向下逐渐增大,所述导架模型上沿纵向设有滑道模型,所述滑道模型内设有挡块模型,所述挡块模型固定设于支架模型上,两个所述连接臂模型之间连接有拉簧模型,两个所述连接臂模型的两端在拉簧模型的作用下分别与导架模型的横向两端和抱爪模型靠近销轴模型的一端相抵接,所述电磁铁模型设于导架模型的下面,所述导架模型在电磁铁模型的作用下能够在所述滑道模型的行程内沿纵向向下滑动,所述导架模型还连接有用于使导架模型向上滑动的复位弹簧模型。
4.根据权利要求2所述的仿真方法,其特征在于,步骤5)中按照如下公式计算出网格节点的新坐标
其中:Δh=δ表示磨损深度,表示磨损方向,A′表示磨损后网格节点的编号坐标矩阵。
5.根据权利要求1所述的仿真方法,其特征在于,步骤2)中所述小球模型的半径小于0.1mm,且其球心与网格节点相重合。
6.根据权利要求1或5所述的仿真方法,其特征在于,步骤2)中采用matlab编写并输出ADAMS的命令流CMD文件的方式来设置小球模型的分布,其中CMD文件的作用为在DAMS工作环境中,通过命令在网格节点处创建无质量小球模型,并在小球模型和导向块模型之间创建接触。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710010871.5A CN106844912B (zh) | 2017-01-06 | 2017-01-06 | 一种基于网格变形技术的机构运动副磨损仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710010871.5A CN106844912B (zh) | 2017-01-06 | 2017-01-06 | 一种基于网格变形技术的机构运动副磨损仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106844912A true CN106844912A (zh) | 2017-06-13 |
CN106844912B CN106844912B (zh) | 2020-07-24 |
Family
ID=59117901
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710010871.5A Active CN106844912B (zh) | 2017-01-06 | 2017-01-06 | 一种基于网格变形技术的机构运动副磨损仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106844912B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107908811A (zh) * | 2017-09-28 | 2018-04-13 | 武汉开目信息技术股份有限公司 | 一种添加弹性挡圈的仿真方法及装置 |
CN108595887A (zh) * | 2018-05-10 | 2018-09-28 | 南京航空航天大学 | 一种聚晶cbn磨粒磨损演变的仿真计算与预测方法 |
CN110287627A (zh) * | 2019-06-28 | 2019-09-27 | 浙江大学 | 一种基于包络的大规模串联传动机构生成方法 |
CN112347590A (zh) * | 2020-11-30 | 2021-02-09 | 北海职业学院 | 一种基于逆向技术的失效零件优化设计方法 |
CN113092299A (zh) * | 2021-04-12 | 2021-07-09 | 湖州中科晟合信息科技有限公司 | 一种接触物体间的磨损量预测方法和装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1534099A (zh) * | 2003-03-28 | 2004-10-06 | 中国科学院力学研究所 | 强化轨迹及分布可控的材料表面激光强化方法和系统 |
US20070168081A1 (en) * | 2005-05-06 | 2007-07-19 | Sung-Ho Shin | Analytic integration of tolerances in designing precision interfaces for modular robotics |
CN203852451U (zh) * | 2014-05-23 | 2014-10-01 | 北京威高亚华人工关节开发有限公司 | 旋转平台型膝关节假体 |
CN104268305A (zh) * | 2014-07-31 | 2015-01-07 | 西安电子科技大学 | 一种基于位置动力学的角度弯曲模型的仿真方法 |
CN104390012A (zh) * | 2014-11-04 | 2015-03-04 | 南京工业大学 | 磁流体润滑的机械密封装置及自适应密封控制方法 |
CN104794266A (zh) * | 2015-04-03 | 2015-07-22 | 上海交通大学 | 基于摩擦功和数值模拟的成形界面损伤预测方法 |
-
2017
- 2017-01-06 CN CN201710010871.5A patent/CN106844912B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1534099A (zh) * | 2003-03-28 | 2004-10-06 | 中国科学院力学研究所 | 强化轨迹及分布可控的材料表面激光强化方法和系统 |
US20070168081A1 (en) * | 2005-05-06 | 2007-07-19 | Sung-Ho Shin | Analytic integration of tolerances in designing precision interfaces for modular robotics |
CN203852451U (zh) * | 2014-05-23 | 2014-10-01 | 北京威高亚华人工关节开发有限公司 | 旋转平台型膝关节假体 |
CN104268305A (zh) * | 2014-07-31 | 2015-01-07 | 西安电子科技大学 | 一种基于位置动力学的角度弯曲模型的仿真方法 |
CN104268305B (zh) * | 2014-07-31 | 2017-03-01 | 西安电子科技大学 | 一种基于位置动力学的角度弯曲模型的仿真方法 |
CN104390012A (zh) * | 2014-11-04 | 2015-03-04 | 南京工业大学 | 磁流体润滑的机械密封装置及自适应密封控制方法 |
CN104794266A (zh) * | 2015-04-03 | 2015-07-22 | 上海交通大学 | 基于摩擦功和数值模拟的成形界面损伤预测方法 |
Non-Patent Citations (1)
Title |
---|
唐琴: "基于摩擦功原理的凸轮-挺住磨损仿真", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107908811A (zh) * | 2017-09-28 | 2018-04-13 | 武汉开目信息技术股份有限公司 | 一种添加弹性挡圈的仿真方法及装置 |
CN107908811B (zh) * | 2017-09-28 | 2021-06-15 | 武汉开目信息技术股份有限公司 | 一种添加弹性挡圈的仿真方法及装置 |
CN108595887A (zh) * | 2018-05-10 | 2018-09-28 | 南京航空航天大学 | 一种聚晶cbn磨粒磨损演变的仿真计算与预测方法 |
CN110287627A (zh) * | 2019-06-28 | 2019-09-27 | 浙江大学 | 一种基于包络的大规模串联传动机构生成方法 |
CN110287627B (zh) * | 2019-06-28 | 2022-09-27 | 浙江大学 | 一种基于包络的大规模串联传动机构生成方法 |
CN112347590A (zh) * | 2020-11-30 | 2021-02-09 | 北海职业学院 | 一种基于逆向技术的失效零件优化设计方法 |
CN113092299A (zh) * | 2021-04-12 | 2021-07-09 | 湖州中科晟合信息科技有限公司 | 一种接触物体间的磨损量预测方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106844912B (zh) | 2020-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106844912A (zh) | 一种基于网格变形技术的机构运动副磨损仿真方法 | |
Vollebregt | Numerical modeling of measured railway creep versus creep-force curves with CONTACT | |
Oh et al. | Capabilities and applications of FEM code DEFORM: the perspective of the developer | |
Podra et al. | Wear simulation with the Winkler surface model | |
CN105569641A (zh) | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 | |
CN105243189B (zh) | 一种车窗密封系统的等效建模及快速评价方法 | |
CN103366377B (zh) | 一种基于形貌特征的微接触特性与图像阈值处理方法 | |
CN107169158A (zh) | 一种基于流固耦合作用的静压滑座工作性能计算方法 | |
CN106886685A (zh) | 一种计及三维分形的界面加载力预测方法 | |
CN102521434A (zh) | 一种轨道交通车辆门系统锁闭装置动力学仿真分析方法 | |
Falconnet et al. | Numerical and experimental analyses of punch wear in the blanking of copper alloy thin sheet | |
CN115470682A (zh) | 基于数字孪生的泵站机组水力部件磨损及寿命预测方法 | |
Zheng et al. | Analysis and prediction of surface wear resistance of ball-end milling topography | |
Thompson et al. | A Proposal for the Calculation of Wear | |
Obermayr et al. | A discrete element model for cohesive soil | |
Hartley et al. | Experimental and theoretical studies of workpiece deformation, stress, and strain during flat rolling | |
Imamovic et al. | Examples of kinematic analysis of complex mechanism using modern software applications | |
Rao | Analysis of aerostatic porous rectangular thrust bearings with offset loads | |
Desai et al. | Review on tribological modeling of worm gear | |
van der Wekken et al. | Local plasticity modelling and its influence on wheel-rail friction | |
Prabhu et al. | Experimental Study and Verification of Wear for Glass Reinforced Polymer using ANSYS | |
Wiklund et al. | Phenomenological friction model in deep drawing of aluminum sheet metals | |
Zhang et al. | Numerical simulation and experimental research on mechanical behaviour of hydraulic disc brakes based on multi-body dynamics | |
Wang et al. | Wear analysis of slideway in emulsion pumps based on finite element method | |
Wang et al. | Thermal elastohydrodynamic lubrication of non-Newtonian fluids with rough surfaces |
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 |