CN109191582A - 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法 - Google Patents
多尺度耦合双正弦波形面弹塑性接触修正模型建立方法 Download PDFInfo
- Publication number
- CN109191582A CN109191582A CN201811064901.1A CN201811064901A CN109191582A CN 109191582 A CN109191582 A CN 109191582A CN 201811064901 A CN201811064901 A CN 201811064901A CN 109191582 A CN109191582 A CN 109191582A
- Authority
- CN
- China
- Prior art keywords
- matrix
- grid
- sinusoidal waveform
- large scale
- scale grid
- 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.)
- Withdrawn
Links
Classifications
-
- 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
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种多尺度耦合双正弦波形面弹塑性接触修正模型建立方法,可用于生成更加高效、精确的弹塑性接触分析模型;通过采用不同的步长将接触区域投影平面划分为不同尺度的网格,确定网格交叉点坐标,结合双正弦波形面的光滑连续性与周期性特点,用大尺度网格控制宏观基体高度与双正弦波形面的幅值,小尺度网格控制双正弦面的光滑连续性,生成粗糙表面的高度矩阵,采用数字编程技术和云点处理技术,在有限元软件中生成粗糙表面三维实体模型与接触平面,并定义非线性材料属性、划分网格、定义接触参数和边界条件,生成多尺度耦合双正弦波形面弹塑性接触修正模型,为多尺度耦合接触问题的研究提供有效的途径。
Description
技术领域
本发明涉及一种结合面接触模型建模方法,尤其涉及一种多尺度耦合双正弦波形面弹塑性接触修正模型建立方法,可用于粗糙表面接触问题的研究。
背景技术
真实的接触表面在微观尺度下都是粗糙的,表面在接触过程中的相互作用对齿轮、轴承、离合器、轮轨等典型机械零部件的接触与摩擦学性能至关重要。接触应力、接触面积与法向变形量等接触特性对机械设备与零件的摩擦、磨损、润滑密封和导热等性能起着关键影响。
接触问题是高度非线性的,其分析模型在理论和算法上都存在着许多有待解决的难点。国内外学者对于粗糙表面的表征方法与接触分析开展了大量的研究,先后基于以下假设条件,如微凸体接触相互独立、微凸体峰顶具有相同的曲率半径、忽略微凸体相互作用和宏观基体变形的影响等,先后提出GW模型、CEB模型、KE模型等经典接触模型。上述模型的假设条件,在一定程度上简化了复杂的接触问题,降低了分析与求解难度,但是与真实接触问题存在较大的差异,影响分析结果的真实性和准确性。
早期研究中建立的经典模型受客观条件的限制,除以上假设前提的局限性,往往忽略微观与宏观的尺度耦合问题。随着计算机仿真技术和数值求解方法的发展,有限元数值模拟成为研究粗糙表面摩擦学特性的重要方法。目前研究大多基于经典的接触模型与粗糙表面数字表征方法,从微观尺度对单一微凸体以及粗糙表面的接触特性进行分析研究。
发明内容
本发明的目的是为了重点解决圆弧面宏观结构参数与微观接触表面形貌参数耦合的建模问题,提供多尺度耦合双正弦波形面弹塑性接触修正模型建立方法。
为了实现上述目的,本发明提供的技术方案为:
多尺度耦合双正弦波形面弹塑性接触修正模型建立方法,该方法包括以下步骤:
(1)根据研究的接触表面几何特征,确定接触表面在水平面xoy的投影区域,记为区域D;定义投影区域D内任一点P坐标为(x,y),其中x∈[a,b],y∈[c,d],Lx,Ly分别为投影区域D在x轴方向和y轴方向的边长,则有Lx=b-a,Ly=c-d;D=x×y={(x,y):a≤x≤b,c≤y≤d};
(2)以ΔxR,ΔyR为步长,分别将区域D的边长Lx,Ly等分得到M和N个点,将区域D等分为(M-1)(N-1)个大尺度网格,投影区域D内大尺度网格各交叉点的坐标记为其坐标值满足 且有M=Lx/ΔxR+1,N=Ly/ΔyR+1,定义由交叉点围成的区域为大尺度网格,记为
(3)以Δxs,Δys为步长,将各大尺度网格再进一步等分为小尺度网格,则分别将区域边长Lx,Ly等分为m和n个点,则投影区域D等分为(m-1)(n-1)个小尺度网格,投影区域D内小尺度网格各交叉点坐标记为其坐标值满足其中i=1,2,···,m,j=1,2,···,n,且有m=Lx/Δxs+1,n=Ly/Δys+1,定义由交叉点 围成的区域为小尺度网格,记为
(4)定义(M-1)×(N-1)阶矩阵R,矩阵不同位置的元素用于存储相应位置的大尺度网格宏观基体法向高度,即在每一个大尺度网格范围内具有确定的统一高度,不同的大尺度网格的宏观基体高度不同;采用表示大尺度网格的宏观基体法向高度,其取值大小分别为:
其中,RI,J表示矩阵R的第I行、第J列元素,1≤I≤M-1,1≤J≤N-1;
(5)确定投影区域D中小尺度网格所在的大尺度网格定义(M-1)×(N-1)阶矩阵A存储每个大尺度网格中的双正弦波形面幅值大小,即在每一个大尺度网格范围内具有确定的统一的幅值,不同大尺度网格内幅值大小不同;定义表示大尺度网格中双正弦波形面的幅值,其取值记为:其中,AI,J表示矩阵A的第I行第J列元素,1≤I≤M-1,1≤J≤N-1;则小尺度网格投影点对应的双正弦波形面上点的法向高度的坐标值为:
其中,沿x轴方向的波长为λx=ΔxR,沿y轴方向的波长为λy=ΔyR;
(6)定义为小尺度网格投影点对应粗糙表面上点相对于水平面xoy的法向高度坐标值为:
(7)以上步骤将粗糙表面转化为一系列空间直角坐标系的坐标点,并进行插值运算令粗糙表面更为连续,最终生成粗糙表面上各点的空间直角坐标系坐标值,获得点云矩阵;将点云矩阵转换为一个三列矩阵,每列分别对应着粗糙表面上各点的x、y、z坐标值,修改上述三列矩阵为五列矩阵,创建空间三维实体的另外五个表面,配合曲面生成三维实体模型;
(8)选择粗糙表面点云中的最高点,并过该点创建与粗糙表面相切的平面,与粗糙表面一起组成接触模型的三维实体模型;
(9)对接触目标体,依次定义非线性材料属性、划分网格、定义接触参数和边界条件,即生成一种多尺度耦合双正弦波形面弹塑性接触修正模型,进行多尺度耦合问题研究。
其中,步骤(4)矩阵R中的对应位置元素服从概率密度函数为的高斯分布,其中,zR为概率密度函数的随机变量,μR和σR为离散后的宏观基体法向高度的均值和标准差。
其中,步骤(5)矩阵A中的对应位置元素服从概率密度函数为高斯分布,其中zA为概率密度函数的随机变量,μA和σA双正弦波形面幅值的均值和标准差。
本发明的有益效果:
本发明采用数字编程技术、云点处理技术、计算机协同仿真技术,在双正弦波形面的基础上,提出一种多尺度耦合的弹塑性接触修正模型。通过对不同参数的选择与控制,可有效解决微宏观尺度耦合问题,对于接触问题的进一步研究,尤其是多尺度耦合问题研究提供了一个有效的途径;方法中大尺度网格和小尺度网格的设置与耦合,宏观基体法向高度与双正弦波形面幅值变量的控制,保证接触表面的表征具有随机特性,符合实际情况;大尺度网格与小尺度网格耦合的方法,解决了工程分析中研究对象的宏观尺度与接触的微观尺度不统一的建模难题,可更加方便的开展接触问题中微宏观耦合问题的研究;双正弦波幅值变量的设置,使得不同区域内接触微凸体峰顶具有不同的曲率半径,更符合实际,避免了现有方法中机械地假设曲率近似相同的缺点,如认为是曲率一致的球体;接触界面的双正弦波具有光滑连续性和周期性特点,即方便了粗糙面的数字建模与表征,也避免了分析模型不成熟的收敛,从而提高了建模、分析的速度与准确性;采用了计算机协同仿真技术,通过数字编程和有限元方法建立的双正弦波形面弹塑性接触修正模型,计及了粗糙表面微凸体接触过程中的相互影响,相比现有方法更符合实际。
附图说明
图1为建模流程示意图;
图2为定义大尺度网格示意图;
图3为定义小尺度网格示意图;
图4为宏观基体高度示意图;
图5为图4的侧视图;
图6为幅值为变量的波形面示意图;
图7为图6的侧视图;
图8为基体高度与波形面耦合示意图;
图9为图8的侧视图;
图10为宏观基体与双正弦波形面叠加示意图;
图11为有限元三维实体模型示意图;
图12为有限元接触模型示意图。
具体实施方式
以下对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅仅用于说明和解释本发明,并不用于限制本发明。
参见图1至图12。多尺度耦合的圆弧面粗糙表面建模方法,具体包括以下几个步骤:
(1)参见图2,根据研究的接触表面几何特征,确定接触表面在水平面xoy的投影区域,记为区域D;定义投影区域D内任一点P坐标为(x,y),其中x∈[a,b],y∈[c,d],Lx,Ly分别为投影区域D在x轴方向和y轴方向的边长,则有Lx=b-a,Ly=c-d;D=x×y={(x,y):a≤x≤b,c≤y≤d}。
(2)参见图2,以ΔxR,ΔyR为步长,分别将区域D的边长Lx,Ly等分得到M和N个点,将区域D等分为(M-1)(N-1)个大尺度网格,投影区域D内大尺度网格各交叉点的坐标记为其坐标值满足
且有M=Lx/ΔxR+1,N=Ly/ΔyR+1,定义由交叉点围成的区域为大尺度网格,记为
(3)参见图3,以Δxs,Δys为步长,将各大尺度网格再进一步等分为小尺度网格,则分别将区域边长Lx,Ly等分为m和n个点,则投影区域D等分为(m-1)(n-1)个小尺度网格,投影区域D内小尺度网格各交叉点坐标记为
其坐标值满足其中i=1,2,···,m,j=1,2,···,n,且有m=Lx/Δxs+1,n=Ly/Δys+1,定义由交叉点 围成的区域为小尺度网格,记为
(4)参见图4和图5,定义(M-1)×(N-1)阶矩阵R,矩阵不同位置的元素用于存储相应位置的大尺度网格宏观基体法向高度,即在每一个大尺度网格范围内具有确定的统一高度,不同的大尺度网格的宏观基体高度不同;采用表示大尺度网格的宏观基体法向高度,其取值大小分别为:其中,RI,J表示矩阵R的第I行第J列元素,1≤I≤M-1,1≤J≤N-1;这里矩阵R中对应位置元素确定,服从概率密度函数为的高斯分布,其中,zR为概率密度函数的随机变量,μR和σR为离散后的宏观基体高度的均值和标准差。
(5)参见图6和图7,确定投影区域D中小尺度网格所在的大尺度网格定义(M-1)×(N-1)阶矩阵A存储每个大尺度网格中的双正弦波幅值大小,即在每一个大尺度网格范围内具有确定的统一的幅值,不同大尺度网格内幅值大小不同;定义表示大尺度网格中双正弦波形面的幅值,其取值记为:
其中,AI,J表示矩阵A的第I行第J列元素,1≤I≤M-1,1≤J≤N-1;这里矩阵A对应位置元素确定,服从概率密度函数为高斯分布,其中,zA为概率密度函数的随机变量,μA和σA双正弦波形面幅值的均值和标准差,则小尺度网格投影点对应的双正弦波形面上点的法向高度的坐标值为:
其中,沿x轴方向的波长为λx=ΔxR,沿y轴方向的波长为λy=ΔyR。
(6)参见图8至图10,定义为小尺度网格投影点对应粗糙表面上点相对于水平面xoy的法向高度坐标值为:
用于表征是宏观基体法向高度和双正弦波形面的法向高度的耦合。
(7)参见图11,以上步骤将粗糙表面转化为一系列空间直角坐标系的坐标点,并进行插值运算令粗糙表面更为连续,最终生成粗糙表面上各点的空间直角坐标系坐标值,获得点云矩阵,将点云矩阵转换为一个三列的矩阵,每列分别对应着粗糙表面上各点的x、y、z坐标值,编写程序修改上述三列矩阵为五列矩阵,将第一列和第二列依次添加字母k和从一开始的自然数编号,后三列存放x、y、z坐标值,生成属性为xls格式文件,并保存为csv格式文件,提取csv文件信息嵌入ANSYS APDL命令流文件,在ANSYS中生成一系列关键点,根据各关键点的位置坐标依次选择并连接任意四个相邻点,运用Coons pathces方法由上述四点间线段构造曲面,使用循环语句生成整个粗糙表面;创建空间三维实体的另外五个表面,配合曲面生成三维实体模型。
(8)参见图12,采用APDL命令,选择粗糙表面点云中的最高点,并过该点创建与粗糙表面相切的平面,与粗糙表面一起组成接触模型的三维实体模型。
(9)对三维实体模型,依次定义非线性材料属性、划分网格、定义接触参数和边界条件,即可生成一种多尺度耦合双正弦波形面弹塑性接触修正模型。
本发明采用数字编程技术、云点处理技术、计算机协同仿真技术,在双正弦波形面的基础上,提出一种多尺度耦合的弹塑性接触修正模型。通过对不同参数的选择与控制,可有效解决微宏观尺度耦合问题,对于接触问题的进一步研究,尤其是多尺度耦合问题研究提供了一个有效的途径;方法中大尺度网格和小尺度网格的设置与耦合,宏观基体法向高度与双正弦波形面幅值变量的控制,保证接触表面的表征具有随机特性,符合实际情况;大尺度网格与小尺度网格耦合的方法,解决了工程分析中研究对象的宏观尺度与接触的微观尺度不统一的建模难题,可更加方便的开展接触问题中微宏观耦合问题的研究;双正弦波幅值变量的设置,使得不同区域内接触微凸体峰顶具有不同的曲率半径,更符合实际,避免了现有方法中机械地假设曲率近似相同的缺点,如认为是曲率一致的球体;接触界面的双正弦波具有光滑连续性和周期性特点,即方便了粗糙面的数字建模与表征,也避免了分析模型不成熟的收敛,从而提高了建模、分析的速度与准确性;采用了计算机协同仿真技术,通过数字编程和有限元方法建立的双正弦波形面弹塑性接触修正模型,计及了粗糙表面微凸体接触过程中的相互影响,相比现有方法更符合实际。
Claims (3)
1.多尺度耦合双正弦波形面弹塑性接触修正模型建立方法,其特征在于,该方法包括以下步骤:
(1)根据研究的接触表面几何特征,确定接触表面在水平面xoy的投影区域,记为区域D;定义投影区域D内任一点P坐标为(x,y),其中x∈[a,b],y∈[c,d],Lx,Ly分别为投影区域D在x轴方向和y轴方向的边长,则有Lx=b-a,Ly=c-d;D=x×y={(x,y):a≤x≤b,c≤y≤d};
(2)以ΔxR,ΔyR为步长,分别将区域D的边长Lx,Ly等分得到M和N个点,将区域D等分为(M-1)(N-1)个大尺度网格,投影区域D内大尺度网格各交叉点的坐标记为其坐标值满足 I=1,2,···,M,J=1,2,···,N,且有M=Lx/ΔxR+1,N=Ly/ΔyR+1,定义由交叉点围成的区域为大尺度网格,记为
(3)以Δxs,Δys为步长,将各大尺度网格再进一步等分为小尺度网格,则分别将区域边长Lx,Ly等分为m和n个点,则投影区域D等分为(m-1)(n-1)个小尺度网格,投影区域D内小尺度网格各交叉点坐标记为其坐标值满足其中i=1,2,···,m,j=1,2,···,n,且有m=Lx/Δxs+1,n=Ly/Δys+1,定义由交叉点 围成的区域为小尺度网格,记为
(4)定义(M-1)×(N-1)阶矩阵R,矩阵不同位置的元素用于存储相应位置的大尺度网格宏观基体法向高度,即在每一个大尺度网格范围内具有确定的统一高度,不同的大尺度网格的宏观基体高度不同;采用表示大尺度网格的宏观基体法向高度,其取值大小分别为:其中,RI,J表示矩阵R的第I行、第J列元素,1≤I≤M-1,1≤J≤N-1;
(5)确定投影区域D中小尺度网格所在的大尺度网格定义(M-1)×(N-1)阶矩阵A存储每个大尺度网格中的双正弦波形面幅值大小,即在每一个大尺度网格范围内具有确定的统一的幅值,不同大尺度网格内幅值大小不同;定义表示大尺度网格中双正弦波形面的幅值,其取值记为:其中,AI,J表示矩阵A的第I行第J列元素,1≤I≤M-1,1≤J≤N-1;则小尺度网格投影点对应的双正弦波形面上点的法向高度的坐标值为:
其中,沿x轴方向的波长为λx=ΔxR,沿y轴方向的波长为λy=ΔyR;
(6)定义为小尺度网格投影点对应粗糙表面上点相对于水平面xoy的法向高度坐标值为:
(7)以上步骤将粗糙表面转化为一系列空间直角坐标系的坐标点,并进行插值运算令粗糙表面更为连续,最终生成粗糙表面上各点的空间直角坐标系坐标值,获得点云矩阵;将点云矩阵转换为一个三列矩阵,每列分别对应着粗糙表面上各点的x、y、z坐标值,修改上述三列矩阵为五列矩阵,创建空间三维实体的另外五个表面,配合曲面生成三维实体模型;
(8)选择粗糙表面点云中的最高点,并过该点创建与粗糙表面相切的平面,与粗糙表面一起组成接触模型的三维实体模型;
(9)对接触目标体,依次定义非线性材料属性、划分网格、定义接触参数和边界条件,即生成一种多尺度耦合双正弦波形面弹塑性接触修正模型,进行多尺度耦合问题研究。
2.根据权利要求1所述的多尺度耦合双正弦波形面弹塑性接触修正模型建立方法,其特征在于:步骤(4)矩阵R中的对应位置元素服从概率密度函数为的高斯分布,其中,zR为概率密度函数的随机变量,μR和σR为离散后的宏观基体法向高度的均值和标准差。
3.根据权利要求1所述的多尺度耦合双正弦波形面弹塑性接触修正模型建立方法,其特征在于:步骤(5)矩阵A中的对应位置元素服从概率密度函数为高斯分布,其中zA为概率密度函数的随机变量,μA和σA双正弦波形面幅值的均值和标准差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811064901.1A CN109191582A (zh) | 2018-09-12 | 2018-09-12 | 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811064901.1A CN109191582A (zh) | 2018-09-12 | 2018-09-12 | 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109191582A true CN109191582A (zh) | 2019-01-11 |
Family
ID=64910745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811064901.1A Withdrawn CN109191582A (zh) | 2018-09-12 | 2018-09-12 | 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109191582A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114757047A (zh) * | 2022-04-28 | 2022-07-15 | 西安交通大学 | 一种轴承钢材料m50合金的多尺度建模计算方法 |
-
2018
- 2018-09-12 CN CN201811064901.1A patent/CN109191582A/zh not_active Withdrawn
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114757047A (zh) * | 2022-04-28 | 2022-07-15 | 西安交通大学 | 一种轴承钢材料m50合金的多尺度建模计算方法 |
CN114757047B (zh) * | 2022-04-28 | 2024-04-02 | 西安交通大学 | 一种轴承钢材料m50合金的多尺度建模计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102040988B1 (ko) | 복수의 3d 곡선들로부터의 표면 생성 | |
CN102779358B (zh) | 用于设计几何三维建模对象的方法及设备 | |
US20080275677A1 (en) | System, methods, and computer readable media, for product design using coupled computer aided engineering models | |
US20080162090A1 (en) | System, methods, and computer readable media, for product design using t-spline deformation | |
US7643026B2 (en) | NURBS surface deformation apparatus and the method using 3D target curve | |
EP0195958A2 (en) | Method of generating finite elements using the symmetric axis transform | |
Montenegro et al. | An automatic strategy for adaptive tetrahedral mesh generation | |
CN109002662A (zh) | 一种多尺度参数耦合的圆弧粗糙表面建模方法 | |
CN105160700B (zh) | 一种用于三维模型重建的截面曲线重构方法 | |
Samulyak et al. | Lagrangian particle method for compressible fluid dynamics | |
CN107038308B (zh) | 一种基于线性内插的规则格网地形建模方法 | |
CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
CN104821005A (zh) | 可模拟笔势风格的汉字书法作品书写过程动态复现的方法 | |
CN115758938A (zh) | 面向粘性边界流场数值模拟的附面层网格生成方法 | |
Whitehurst | A free Lagrange method for gas dynamics | |
Frade et al. | Evolution of artificial terrains for video games based on obstacles edge length | |
Sun et al. | Adaptive generation and local refinement methods of three-dimensional hexahedral element mesh | |
CN101694726A (zh) | 一种基于多源地形数据的融合与绘制方法 | |
CN109191582A (zh) | 多尺度耦合双正弦波形面弹塑性接触修正模型建立方法 | |
CN117113772A (zh) | 一种漂浮式风机的混合动网格模拟方法及装置 | |
CN105302951A (zh) | 一种有限元网格曲面剖分方法 | |
JP2603902B2 (ja) | 自動メッシュ生成方法及びシステム | |
CN105427362A (zh) | 一种快速aiap形状插值算法 | |
Hung et al. | Automatic clustering method for real-time construction simulation | |
Zhang et al. | Conformal adaptive hexahedral-dominant mesh generation for CFD simulation in architectural design applications |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20190111 |