CN109492242A - 一种远距离极坐标约束有限元分析方法 - Google Patents

一种远距离极坐标约束有限元分析方法 Download PDF

Info

Publication number
CN109492242A
CN109492242A CN201811004026.8A CN201811004026A CN109492242A CN 109492242 A CN109492242 A CN 109492242A CN 201811004026 A CN201811004026 A CN 201811004026A CN 109492242 A CN109492242 A CN 109492242A
Authority
CN
China
Prior art keywords
rolling element
coordinate system
finite element
remote
constrained
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
CN201811004026.8A
Other languages
English (en)
Other versions
CN109492242B (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.)
715th Research Institute of CSIC
Original Assignee
715th Research Institute of CSIC
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 715th Research Institute of CSIC filed Critical 715th Research Institute of CSIC
Priority to CN201811004026.8A priority Critical patent/CN109492242B/zh
Publication of CN109492242A publication Critical patent/CN109492242A/zh
Application granted granted Critical
Publication of CN109492242B publication Critical patent/CN109492242B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • 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)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Rolling Contact Bearings (AREA)

Abstract

本发明公开了一种远距离极坐标约束有限元分析方法，是基于有限元理论建立球对平面有限元模型，并对目标进行了远距离极坐标约束，同时与Hertz接触理论进行了对比。结果表明随着极坐标系远离被约束模型，所得到的解与理论解越来越接近，当约束坐标系和被约束模型远到一定程度时，结果几乎完全一样。这为该方法推广到轴承有限元求解上提供了一个很好的思路，尤其是为带有游隙的轴承有限元分析提供了一个能够求解且能快速收敛的方法。

Description

一种远距离极坐标约束有限元分析方法
技术领域
本发明涉及有限元分析方法，主要是一种远距离极坐标约束有限元分析方法。
背景技术
在静态分析中，开始不连接在一起的两个(或多个)物体(如图1)，在建立接触前，会产生刚体运动。如果在求解过程中的任一时刻，两物体没有联系，刚度矩阵就会奇异，有限元分析将产生出负主元，最终导致求解无法进行。要想克服以上问题所带来的刚体位移，以往常常利用以下几种方法：
a.在恰好碰上的位置上建立几何模型。但对于解决位置不确定或表面是曲面及不规则的面，这样求解仍然很难进行下去；另外由于有限元网格修整，仍不可避免存在间隙和小侵入，这样同样会带来求解难以进行下去。
b.动力学，在动力学分析中，惯性作用可以阻止刚体运动，克服刚体运动的一种选择是动态求解。这需要加上质量和阻尼使求解从静态转化为动态，但要缓冲掉不想要的动态影响不总是那么容易，另外这种方法求解效率也比较慢。
c.位移控制，该方法利用三个载荷步：载荷步1，用一个小位移使物体进入初始接触状态；载荷步2，从位移控制转换到力控制，删除强加的位移，加上反作用力，用一个子步求解；载荷步3，继续进行求解加载。但对于复杂加载，要施加哪个位移约束并不是那么明显。因此给求解带来了困难。
d.软弹簧，该方法用接地软弹簧来求解刚体运动，与系统刚度相比弹簧刚度应该小到可以被忽略，以确保对求解不会产生影响，但这需要反复的实验，对于复杂问题求解仍然有很大的困难。
发明内容
本发明的目的在于克服现有技术存在的不足，而提供一种远距离极坐标约束有限元分析方法,能够在不影响精度的情况下，实现球对平面(或类似结构)非接触的快速求解。尤其对于求解带有游隙的轴承保证求解精度的情况下，快速求解极其有帮助。
本发明的目的是通过如下技术方案来完成的。这种远距离极坐标约束有限元分析方法，该方法是基于有限元原理，建立球对平面或轴承模型，同时建立远距离极坐标系，并将球或轴承滚动体约束在极坐标系上(1)，通过对球或者轴承滚动体的约束，提高求解的收敛性；在外部边界条件作用下，不断调整远距离坐标系的距离(2)，而求解精度随着坐标系远离球或轴承滚动体的过程中求解精度将会不断地提高，使求解精度变的可控。
更进一步的，远距离坐标约束的方法：为每个滚动体建立一个局部坐标系坐标系位置为沿滚动体中心到坐标原点的射线上，调节局部坐标位置并逐渐远离滚动体中心，随着局部坐标系距对应滚动体距离的增大，内圈与滚动体接触力逐渐趋近外圈与滚动体接触力，最终内外圈与滚动体接触力几乎相等；通过合理调节局部坐标距滚动体中心的距离能够提高求解精度，并能够保证求解的收敛性。
如图1所示，R为球体半径，λ为球体与接触平台的间隙，D为球心O1距局部坐标系原点O的距离，α为点N1与原点O的连线与极坐标轴X的夹角，β为N2与原点O的连线与极坐标轴X的夹角，F为施加到球体上的外力。边界条件：(1)将整个球体(或球体上任一两点)周向约束于局部极坐标系，则α和β值将保持不变；(2)将平台两端施加全约束；(3)球体上施加如图1所示的力F。通过改变距离D(可以通过改变局部坐标系的位置)，则接触平台的法向接触的合力F1将逐渐变大，最终和力F趋于相同，而球体和平面的接触力也逐渐和理论解相同。
结果表明该方法不仅能够有效的求得精确的解，而且能够使收敛速率大大的提高。
由于滚动轴承主要由内、外圈及滚动体组成，未受力前难免部分滚动体处于不完全约束状态。受力过程中部分滚动体逐渐与内外圈接触并最终产生接触力，最后仍不可避免一部分滚动体仍处于不完全约束状态，这为有限元求解带来了一定的困难。许多相关文献的作者往往引入一些附加的约束以提高收敛性及收敛速度，但这却牺牲了求解的精度甚至正确性。针对这个难题本文提出远距离极坐标约束方法。首先基于有限元软件建立远距离坐标系及球对平面有限元模型，并对目标进行了远距离极坐标约束，同时与Hertz接触理论进行了对比。通过这种方法来证明远距离坐标系约束求解的精度和准确性。然后建立球轴承有限元模型运用远距离坐标约束方法与理论求解接触力进行对比。最后得出远距离极坐标约束法的有效性。
本发明的有益效果为：该方法基于有限元软件建立球对平面有限元模型，并对目标进行了远距离极坐标约束，同时与Hertz接触理论进行了对比。结果表明随着约束坐标系远离被约束模型，与理论解越来越接近，当约束坐标系和被约束模型远到一定程度时，结果几乎完全一样。这为该方法推广到轴承有限元求解上提供了一个很好的思路，尤其是为带有游隙的轴承有限元分析提供了一个能够且能快速收敛的方法。
(1)提高了求解精度，并使求解精度可控。
(2)提高了求解速度。
(3)提供了带有间隙物体接触的静应力的求解方法。
附图说明
图1为带有间隙的球对平面静态分析模型；
图2为极坐标约束示意图；
图3为轴承有限元分析爆炸图；
图4为有游隙时轴承受力示意图；
图5为轴承的受力分析图；
图6为单个滚动体的受力分析图；
图7为轴承滚动体与内外圈接触力；
图8局部坐标系建立的示意图；
图9为滚动体与内圈环向接触应力随距离的变化；
图10为滚动体与外圈环向接触应力随距离的变化；
图11为滚动体分别与内圈和外圈接触力随距离的变化。
具体实施方式
下面将结合附图对本发明做详细的介绍：
采用远距离坐标约束在轴承上的应用原理及求解结果
1有限元模型的建立
所求解的模型为NJP3226轴承，其中该圆柱滚子轴承的结构参数为：
表3-1单个轴承尺寸示意图
图1为有限元模型，轴承材料采用GCr15，额定工作载荷为30kN，其弹性模量为200GPa，松比为0.3。约束方法：
1)内圈内表面全约束；
2)滚动体周向约束于局部坐标系上；
3)外圈外表面受余弦分布载荷p(θ)＝2.63cos[1.087(θ-90)](7.2°≤θ≤172.8°)
Fr＝∫∫p(θ)cos(θ)rdθdx＝B∫p(θ)cos(θ)rdθ (3-1)
2滚子无周向约束
2.1有游隙时
用有限元进行圆柱滚子轴承建模时，由于保持架作用力比较小，而且保持架建模比较复杂。且往往忽略保持架的作用，实际情况是，内圈固定，外圈加外部载荷，但由于轴承受力，接触区会产生一定的变形，进而导致非承载区滚动体出现一定的间隙，使部分滚动体处于无约束状态，这样往往会导致计算无法收敛。
特别是当滚动体与内外圈间有游隙时，假设内圈固定外圈受外部载荷，滚动体与外圈间有游隙时，则外圈处于不静定状态，计算往往不能够收敛；假设滚动体与内圈间有游隙时，则外载荷作用于外圈，外圈将力传递到滚动体上，此时滚动体处于不静定状态，计算时往往也不能够收敛。
2.2无游隙时
如图所示为滚子受力示意图假设计算能够收敛，即滚动体能够满足平衡条件。
由力的平衡得滚子与外圈接触力：
Qeq＝Qiqcosθ+Ff2sinθ (3-2)
式中Qiq为滚子与内圈的接触力；Ff2为滚子与内圈的摩擦力。
滚子与外圈的摩擦力：
Ff1＝Qiqsinθ+Ff2cosθ (3-3)
Ff1＝fQeq (3-4)
滚子与内圈的摩擦力：
Ff2＝fQiq (3-5)
式中f为摩擦系数
由上可得
Qeq＝Qiqcosθ+fQiqsinθ (3-6)
Qeq＝Qiqcosθ-Qiqsinθ/f (3-7)
两式相减得到
Qiqsinθ(f2+1)＝0 (3-8)
其中Qiq≠0，f2+1≠0，sinθ≠0因此假设不成立，即滚动体不能够满足平衡条件。
这是在很多有限元数值计算中，没有保持架的圆柱滚子轴承很难收敛的原因。而在很多计算中习惯将圆柱滚子约束到轴承几何中心上，下一节将进行讨论。
3滚子有周向约束
3.1将滚子约束到几何中心上
在很多计算中习惯将滚动体周向约束于几何中心上，之所以这样约束的原因是很多作者认为：1，将滚子约束到几何中心上这样可以减少计算资源，加快求解速度。2，省略了轴承的保持架，而保持架的作用就是保持滚动体之间的距离，防止滚动体产生滑动，为达到同样的效果，可以在圆柱坐标系下，约束滚动体节点的切向UY。但这种约束方法本身存在一定的缺陷，因为外力作用于外圈时，力将会通过外圈将力传递到滚动体上，而滚动体在力的作用下将沿着节点初始位置向坐标原点直线运动，但由于滚子，而这样与滚子真实的结果会有很大不同。对于真实结果，根据力的平衡可知，滚子与内圈接触力应等于滚子与外圈接触力，即Qiq＝Qeq。
但经有限元计算发现，滚子与内圈接触力并不等于滚子与外圈接触力，如图3-5所示。滚子与内圈接触力小于与外圈与滚动体约束力。因此这种约束明显带着一定的缺陷。
3.2采用远距离坐标约束
方法：为每个滚动体建立一个局部坐标系坐标系位置为沿滚动体中心到坐标原点的射线上，调节局部坐标位置并逐渐远离滚动体中心，可得到这样的结果：由图8可知，当局部坐标系在滚子中心时，内圈环向没有应力的作用，而当局部坐标距离逐渐远离整体坐标系中心时，内圈所受Von mises应力逐渐增大，以致于最后该应力趋近于确定值。同时图9可知外圈与滚动体的Von mises应力不随局部坐标系距离的远离而变化。分析原因可能是由于力直接作用于外圈，外圈并没有受到一定的约束干扰，因此外圈受力不会因为局部坐标距几何中心距离的变化而变化。而图10可知内圈与滚动体的接触力随局部坐标系距离的远离而增大以至于最后达到和外圈与滚动体接触力相等，这一点和理论公认的内圈与滚动体接触力和外圈与滚动体接触力相等一致。因此合理调节局部坐标距滚动体中心的距离能够提高求解精度，并能够保证求解的收敛性。
图11从分析结果上可以看出，随着局部坐标系距对应滚动体距离的增大，内圈与滚动体接触力逐渐趋近外圈与滚动体接触力，最终内外圈与滚动体接触力几乎相等，而过去一些学者将滚动体周向约束带来内外圈接触力不相等的结果有很大程度的改善。说明这种方法具有一定的优越性。
可以理解的是，对本领域技术人员来说，对本发明的技术方案及发明构思加以等同替换或改变都应属于本发明所附的权利要求的保护范围。

Claims (2)

1.一种远距离极坐标约束有限元分析方法，其特征在于：该方法是基于有限元原理，建立球对平面或轴承模型，同时建立远距离极坐标系，并将球或轴承滚动体约束在极坐标系上，通过对球或者轴承滚动体的约束，提高求解的收敛性；在外部边界条件作用下，不断调整远距离坐标系的距离，而求解精度随着坐标系远离球或轴承滚动体的过程中求解精度将会不断地提高。
2.根据权利要求1所述的远距离极坐标约束有限元分析方法，其特征在于：远距离坐标约束的方法：为每个滚动体建立一个局部坐标系坐标系位置为沿滚动体中心到坐标原点的射线上，调节局部坐标位置并逐渐远离滚动体中心，随着局部坐标系距对应滚动体距离的增大，内圈与滚动体接触力逐渐趋近外圈与滚动体接触力，最终内外圈与滚动体接触力几乎相等；通过合理调节局部坐标距滚动体中心的距离能够提高求解精度，并能够保证求解的收敛性。
CN201811004026.8A 2018-08-30 2018-08-30 一种远距离极坐标约束有限元分析方法 Active CN109492242B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811004026.8A CN109492242B (zh) 2018-08-30 2018-08-30 一种远距离极坐标约束有限元分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811004026.8A CN109492242B (zh) 2018-08-30 2018-08-30 一种远距离极坐标约束有限元分析方法

Publications (2)

Publication Number Publication Date
CN109492242A true CN109492242A (zh) 2019-03-19
CN109492242B CN109492242B (zh) 2023-04-07

Family

ID=65690364

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811004026.8A Active CN109492242B (zh) 2018-08-30 2018-08-30 一种远距离极坐标约束有限元分析方法

Country Status (1)

Country Link
CN (1) CN109492242B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114510797A (zh) * 2022-01-25 2022-05-17 哈尔滨理工大学 一种无保持架轴承控制滚动体离散的方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009150434A (ja) * 2007-12-19 2009-07-09 Ntn Corp 転がり軸受の保持器およびその設計方法
CN102909643A (zh) * 2012-10-22 2013-02-06 哈尔滨工业大学 磁场分布均匀化的小直径永磁体球形抛光头及该永磁体球形抛光头结构参数优化设计方法
CN103116672A (zh) * 2013-01-30 2013-05-22 河南科技大学 利用有限元建模判断转盘轴承是否满足要求的方法
CN103530468A (zh) * 2013-10-23 2014-01-22 清华大学 一种考虑轴承刚度耦合非线性的多支撑轴系有限元方法
CN103927428A (zh) * 2014-05-09 2014-07-16 清华大学 一种考虑多因素影响的锥齿轮错位量有限元计算方法
JP2015032097A (ja) * 2013-08-01 2015-02-16 Ntn株式会社 転がり軸受の動力学解析方法および解析装置
CN105608297A (zh) * 2016-02-17 2016-05-25 中国科学院沈阳自动化研究所 一种角接触球轴承动刚度的数值模拟方法
CN106560815A (zh) * 2016-02-02 2017-04-12 梁明轩 球轴承可靠性设计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009150434A (ja) * 2007-12-19 2009-07-09 Ntn Corp 転がり軸受の保持器およびその設計方法
CN102909643A (zh) * 2012-10-22 2013-02-06 哈尔滨工业大学 磁场分布均匀化的小直径永磁体球形抛光头及该永磁体球形抛光头结构参数优化设计方法
CN103116672A (zh) * 2013-01-30 2013-05-22 河南科技大学 利用有限元建模判断转盘轴承是否满足要求的方法
JP2015032097A (ja) * 2013-08-01 2015-02-16 Ntn株式会社 転がり軸受の動力学解析方法および解析装置
CN103530468A (zh) * 2013-10-23 2014-01-22 清华大学 一种考虑轴承刚度耦合非线性的多支撑轴系有限元方法
CN103927428A (zh) * 2014-05-09 2014-07-16 清华大学 一种考虑多因素影响的锥齿轮错位量有限元计算方法
CN106560815A (zh) * 2016-02-02 2017-04-12 梁明轩 球轴承可靠性设计方法
CN105608297A (zh) * 2016-02-17 2016-05-25 中国科学院沈阳自动化研究所 一种角接触球轴承动刚度的数值模拟方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
SARABJEETSINGH: "Analysesofcontactforcesandvibrationresponse", 《ELSEVIER》 *
吴云鹏,: "滚动轴承力学模型的研究及其进展", 《煤矿机械》 *
甘伟等: "基于ABAQUS的深沟球轴承接触应力有限元分析", 《科技通报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114510797A (zh) * 2022-01-25 2022-05-17 哈尔滨理工大学 一种无保持架轴承控制滚动体离散的方法

Also Published As

Publication number Publication date
CN109492242B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
US20180153444A1 (en) Body movement tracking
US11573086B2 (en) Electronic apparatus of estimation of movement direction and method thereof
JP2012502331A (ja) ロール補正を備えたハンドヘルドポインティングデバイス
Chen et al. Effects of geometric parameters on static load capacity of a double-row four-point contact ball bearing
CN109492242A (zh) 一种远距离极坐标约束有限元分析方法
Li et al. Dynamic characteristics of ball bearing with flexible cage lintel and wear
Halpin et al. An analytical model of four-point contact rolling element ball bearings
JPWO2014156385A1 (ja) 測位装置とその制御方法、および、プログラム
Du et al. A self-adaptive alignment strategy for large components based on dynamic compliance center
Liu et al. Mass and mass center identification of target satellite after rendezvous and docking
Wang et al. Three-dimensional wear prediction of four-degrees-of-freedom parallel mechanism with clearance spherical joint and flexible moving platform
JP2008191025A (ja) 共形接触下の接触面圧および表面下応力計算方法
Harsha et al. Effects of preload and number of balls on nonlinear dynamic behavior of ball bearing system
Bons et al. Compensating for soft-tissue artifact using the orientation of distal limb segments during electromagnetic motion capture of the upper limb
You et al. Load distribution calculation of flexible ball bearing with elliptical cam wave generator
Szumiński Determination of the stiffness of rolling kinematic pairs of manipulators
CN108592835B (zh) 一种基于应变的天线副面位姿实时测量方法
Qiu et al. Static analysis on slewing bearing with negative clearance based on ANSYS
汪坤 et al. IMU-based attitude optimization of the robotic arm end effector
Zhang et al. Research on adaptive trajectory following of lower limb rehabilitation exoskeleton based on sliding mode control
Yuan et al. Optimization design of the web press’s fold mechanism based on robustness
Zhang et al. Research on Dynamic Equation Construction Method with Clearance
Damian et al. The influence of the misalignment on load distribution in angular contact ball bearings
Zhao et al. Dynamic analysis of the rotating thin-walled structure with absolute nodal coordinate formulation
Liu et al. Application of Flores model considering variable recovery coefficient in dynamics analysis of ball bearings

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