CN116186933A - 一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 - Google Patents
一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 Download PDFInfo
- Publication number
- CN116186933A CN116186933A CN202310176497.1A CN202310176497A CN116186933A CN 116186933 A CN116186933 A CN 116186933A CN 202310176497 A CN202310176497 A CN 202310176497A CN 116186933 A CN116186933 A CN 116186933A
- Authority
- CN
- China
- Prior art keywords
- piston
- moment
- pin
- oil film
- skirt
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005461 lubrication Methods 0.000 title claims abstract description 28
- 238000004088 simulation Methods 0.000 title claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 23
- 230000010354 integration Effects 0.000 claims abstract description 7
- 238000006073 displacement reaction Methods 0.000 claims description 11
- 239000007789 gas Substances 0.000 claims description 9
- 230000005484 gravity Effects 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 5
- 238000004880 explosion Methods 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 239000012809 cooling fluid Substances 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000002737 fuel gas Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Images
Classifications
-
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Shafts, Cranks, Connecting Bars, And Related Bearings (AREA)
Abstract
本发明的目的在于提供一种考虑活塞与缸套热变形的裙部润滑特性仿真方法,包括以下步骤:设定活塞二阶运动变量初值,开始曲轴转角循环;结合活塞型线计算油膜厚度;计算油膜压力和微凸体接触压力,利用高斯‑勒让德积分得到油膜承载力和微凸体接触力;通过动力学方程计算得到当前时刻连杆对销的力和力矩以及其它力和力矩,对动力学方程进行求解;判段曲轴转角是否满足循环终止条件;判断本次工作循环是否满足收敛条件,若不满足更新二阶运动变量初值进行下一计算循环,若满足计算完成输出结果。本发明克服了爆发压力附近和最小油膜数值较小等特殊条件下计算结果失真和求解不稳定的问题。为对实际柴油机活塞裙部润滑预测及优化设计提供一定支撑。
Description
技术领域
本发明涉及的是一种内燃机仿真方法,具体地说是内燃机摩擦副润滑特性仿真方法。
背景技术
在各类动力机械中,柴油机以其热效率高、扭矩大、经济性好、功率范围广、启动迅速、维修简单、运行安全、使用寿命长等优点而得到广泛应用,主要用于重型汽车、大型客车、工程机械、坦克、船舶、发电机组等,在国民经济和国防建设中处于重要地位。活塞作为柴油机最重要的零部件之一,其与缸套之间的润滑状态对柴油机动力性和可靠性的保证十分关键。
活塞在柴油机工作过程中不但要承受气缸压力、惯性力、侧推力和摩擦力等机械负荷的作用,还要承受来自高温燃气的热负荷,进而使得活塞和缸套发生热变形、活塞裙-缸套摩擦副润滑状态改变。出于研究对象爆发压力较低的原因,过去相当一部份的研究不考虑活塞与缸套的热变形或者仅考虑活塞的热变形,但近年来随着柴油机向着高强化、高功率密度方向发展,柴油机的机械负荷和热负荷得到了进一步提高,活塞与缸套发生了严重的热变形,对摩擦副的配合间隙与动力油膜产生了极大的影响。在这种更为严苛的工作条件下裙部型面微小的不合理设计会显著增加活塞与缸套的磨损,甚至会出现拉缸等故障。因此考虑活塞与缸套热变形,对活塞裙部润滑特性更准确的计算是十分必要的。
在考虑活塞和缸套热变形对裙部润滑计算时,若将每个节点对应的变形计算结果均作为输入,会使得数据处理量巨大且计算时间较长,为降低计算成本、减少设计校核时间有必要对热变形数据做适当的简化,根据李闯等在文章《发某柴油机活塞裙部变形分析》和毕玉华等在文章《冷却液流动均匀性对缸套热变形的影响研究》的研究结果,活塞轴线截面上的径向变形为椭圆形、缸套径向变形不均匀呈豌豆形,需对活塞和缸套采用不同的简化方法。
综上所述,为准确的计算活塞裙部润滑特性,考虑活塞与缸套的热变形是十分必要的,同时为兼顾准确性和计算成本,对计算方法进行简化和优化又十分重要。
发明内容
本发明的目的在于提供能够更准确、高效的对活塞裙部接触和润滑状态进行预测,进而大幅度提高活塞裙部设计精度和效率的一种考虑活塞与缸套热变形的裙部润滑特性仿真方法。
本发明的目的是这样实现的:
本发明一种考虑活塞与缸套热变形的裙部润滑特性仿真方法,其特征是:
(1)设定活塞二阶运动变量初值,开始曲轴转角循环;
(2)对活塞和缸套热变形量进行处理,结合活塞型线计算油膜厚度;
(3)通过混合润滑模型计算油膜压力和微凸体接触压力,利用高斯-勒让德积分得到油膜承载力和微凸体接触力及力矩;
(4)根据已知的当前时刻缸内气体力和活塞环摩擦力,通过动力学方程计算得到当前时刻连杆对销的力和力矩以及其它力和力矩,采用变步长四阶龙格库塔法对动力学方程进行求解,由当前时刻所受的力及力矩得到下一时刻的二阶运动变量值;
(5)令当前时间和曲轴转角均增加一步,判段曲轴转角是否满足循环终止条件,若不满足更新当前时刻变量值继续计算,若满足进入下一步;
(6)判断本次工作循环是否满足收敛条件,若不满足更新二阶运动变量初值进行下一工作循环计算,若满足计算完成输出结果。
本发明还可以包括:
2、步骤二具体为:采用沿活塞中心轴方向长、短轴不断变化的椭圆面对考虑活塞热变形后的活塞型面进行简化,采用沿缸套中心轴方向长、短轴不断变化的4个1/4椭圆面对考虑缸套热变形的缸套内表面进行简化,具体的油膜厚度计算方法为:
确定活塞型线凸点处的曲率Rp、缸套半径R、活塞推力侧和非推力侧热变形后的轴向型线f1,0°(y)和f1,90°(y),缸套周向角度为0°、90°、180°和270°处的缸套轴向变形曲线d1,0°(y)、d1,90°(y)、d1,180°(y)和d1,270°(y),计算得到裙部相对于凸点处的缩减量和缸套某处相对于冷态的变形量为:
其中yp为活塞裙部顶端距上止点的轴向位移;
计算由活塞二阶运动导致的推力面油膜厚度变化量:
其中,a为活塞重心相对于活塞裙顶部的轴向位置,b为活塞销中心相对于活塞重心是轴向位置;
得到裙部与缸套之间的油膜厚度为:
3、步骤三具体为:
利用经过Vogenpohl变换的平均雷诺方程计算油膜压力,方程为:
计算微凸体接触压力,方程为:
其中,E′为活塞裙和缸套综合弹性模量,η为峰元密度,β为峰元曲率半径,H为膜厚比(h/σ);
采用高斯-勒让德积分得到油膜承载力和微凸体接触力及力矩。
4、步骤四具体为:
根据已知的当前时刻缸内气体力和活塞环摩擦力以及与油膜和微凸体相关的力及力矩,得到连杆作用于活塞销的力,再由活塞销与连杆之间的摩擦系数得到连杆与活塞销之间的摩擦力矩,利用建立的活塞裙部动力学方程计算当前时刻的活塞销横向运动加速度和活塞摆动加速度动力学方程为:
其中,Mpin=μpinFL,μ为摩擦系数,下角标f代表与摩擦力相关的力和力矩,下角标c代表与微凸体接触相关的力和力矩,下角标j代表与角点接触相关的力和力矩,下角标h代表与油膜动压润滑相关的力和力矩,G代表与缸内气体相关的力和力矩,L代表连杆作用于活塞销相关的力和力矩,pin代表和活塞销相关的力和力矩,Ip为活塞绕重心的转动惯量;
根据当前时刻的销横向运动加速度和摆动加速度采用变步长四阶龙格库塔法计算下一时刻的销横向运动速度和摆动速度同时根据上一时刻计算出的当前时刻销横向运动速度和摆动速度采用变步长四阶龙格库塔法计算下一时刻的销横向位移xp和摆动角度β。
本发明的优势在于:本发明考虑了实际柴油机稳定运行时发生的活塞裙部和缸套热变形对裙部润滑的影响,增加了活塞销摩擦力矩、活塞环对活塞的作用等受力因素,提出了较为全面、科学的裙部润滑特性计算方法;运用了经过Vogenpohl变换的平均雷诺方程和变步长四阶龙格库塔法等优化算法,利用相同的计算资源,克服了爆发压力附近和最小油膜数值较小等特殊条件下计算结果失真和求解不稳定的问题。为对实际柴油机活塞裙部润滑预测及优化设计提供一定支撑。
附图说明
图1为本发明的流程图;
图2为油膜厚度沿轴向变化示意图;
图3为油膜厚度沿周向变化示意图;
图4为活塞受力示意图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1-4,本发明包括以下步骤:
步骤一、设定活塞二阶运动变量初值,开始曲轴转角循环。
步骤二、对活塞和缸套热变形量进行处理,结合活塞型线计算油膜厚度。
采用沿活塞中心轴方向长、短轴不断变化的椭圆面对考虑活塞热变形后的活塞型面进行简化,采用沿缸套中心轴方向长、短轴不断变化的4个1/4椭圆面对考虑缸套热变形的缸套内表面进行简化,具体的油膜厚度计算方法为:
首先,确定活塞型线凸点处的曲率Rp,缸套半径R,活塞推力侧和非推力侧热变形后的轴向型线f1,0°(y)和f1,90°(y),缸套周向角度为0°、90°、180°和270°处的缸套轴向变形曲线d1,0°(y)、d1,90°(y)、d1,180°(y)和d1,270°(y),计算得到裙部某处相对于凸点处的缩减量和缸套某处相对于冷态的变形量为:
其中yp为活塞裙部顶端距上止点的轴向位移。
其次,计算由活塞二阶运动导致的推力面油膜厚度变化量:
其中,a为活塞重心相对于活塞裙顶部的轴向位置,b为活塞销中心相对于活塞重心是轴向位置。
最后,得到裙部某处与缸套之间的油膜厚度为:
步骤三、通过混合润滑模型计算油膜压力和微凸体接触压力,后利用高斯-勒让德积分得到油膜承载力和微凸体接触力及力矩。
首先,利用经过Vogenpohl变换的平均雷诺方程计算油膜压力,方程为:
其次,计算微凸体接触压力,方程为:
其中,E′为活塞裙和缸套综合弹性模量,η为峰元密度,β为峰元曲率半径,H为膜厚比(h/σ)。
最后,采用高斯-勒让德积分得到油膜承载力和微凸体接触力及力矩。
在步骤三中,将经过Vogenpohl变换的平均雷诺方程应用于混合润滑模型计算油膜压力,修正了网格不是非常细密情况下经典平均雷诺方程在最小油膜厚度很小时计算结果严重失真的问题。
步骤四、根据已知的当前时刻缸内气体力和活塞环摩擦力,通过动力学方程计算得到当前时刻连杆对销的力和力矩以及其它力和力矩,采用变步长四阶龙格库塔法对动力学方程进行求解,由当前时刻所受的力及力矩得到下一时刻的二阶运动变量值。
首先,根据已知的当前时刻缸内气体力和活塞环摩擦力以及上一步计算出来的与油膜和微凸体相关的力及力矩,得到连杆作用于活塞销的力,再由活塞销与连杆之间的摩擦系数得到连杆与活塞销之间的摩擦力矩,后利用建立的活塞裙部动力学方程计算当前时刻的活塞销横向运动加速度和活塞摆动加速度动力学方程为:
其中,Mpin=μpinFL,μ为摩擦系数,下角标f代表与摩擦力相关的力和力矩,下角标c代表与微凸体接触相关的力和力矩,下角标j代表与角点接触相关的力和力矩,下角标h代表与油膜动压润滑相关的力和力矩,G代表与缸内气体相关的力和力矩,L代表连杆作用于活塞销相关的力和力矩,pin代表和活塞销相关的力和力矩,Ip为活塞绕重心的转动惯量。
其次,根据当前时刻的销横向运动加速度和摆动加速度采用变步长四阶龙格库塔法计算下一时刻的销横向运动速度和摆动速度同时根据上一时刻计算出的当前时刻销横向运动速度和摆动速度采用变步长四阶龙格库塔法计算下一时刻的销横向位移xp和摆动角度β。
在步骤四中,动力学方程中考虑了活塞销与活塞之间的摩擦力矩、活塞环对活塞的横向摩擦力等力及力矩,并采用变步长四阶龙格库塔法对动力学方程进行求解,对爆发压力附近的计算步长进行自动加密。
步骤五、令当前时间和曲轴转角均增加一步,判段曲轴转角是否满足循环终止条件,若不满足更新当前时刻变量值继续计算,若满足进入下一步;
步骤六、判断本次工作循环是否满足收敛条件,若不满足更新二阶运动变量初值进行下一计算循环,若满足计算完成输出结果。
以活塞销横向位移和活塞摆动角度作为收敛判定参数,当不满足收敛条件时,将计算出最大曲轴转角处的销横向位移、摆动角度、销横向速度、摆动速度作为下一计算循环初值。
Claims (5)
1.一种考虑活塞与缸套热变形的裙部润滑特性仿真方法,其特征是:
(1)设定活塞二阶运动变量初值,开始曲轴转角循环;
(2)对活塞和缸套热变形量进行处理,结合活塞型线计算油膜厚度;
(3)通过混合润滑模型计算油膜压力和微凸体接触压力,利用高斯-勒让德积分得到油膜承载力和微凸体接触力及力矩;
(4)根据已知的当前时刻缸内气体力和活塞环摩擦力,通过动力学方程计算得到当前时刻连杆对销的力和力矩以及其它力和力矩,采用变步长四阶龙格库塔法对动力学方程进行求解,由当前时刻所受的力及力矩得到下一时刻的二阶运动变量值;
(5)令当前时间和曲轴转角均增加一步,判段曲轴转角是否满足循环终止条件,若不满足更新当前时刻变量值继续计算,若满足进入下一步;
(6)判断本次工作循环是否满足收敛条件,若不满足更新二阶运动变量初值进行下一工作循环计算,若满足计算完成输出结果。
3.根据权利要求1所述的一种考虑活塞与缸套热变形的裙部润滑特性仿真方法,其特征是:步骤二具体为:采用沿活塞中心轴方向长、短轴不断变化的椭圆面对考虑活塞热变形后的活塞型面进行简化,采用沿缸套中心轴方向长、短轴不断变化的4个1/4椭圆面对考虑缸套热变形的缸套内表面进行简化,具体的油膜厚度计算方法为:
确定活塞型线凸点处的曲率Rp、缸套半径R、活塞推力侧和非推力侧热变形后的轴向型线f1,0°(y)和f1,90°(y),缸套周向角度为0°、90°、180°和270°处的缸套轴向变形曲线d1,0°(y)、d1,90°(y)、d1,180°(y)和d1,270°(y),计算得到裙部相对于凸点处的缩减量和缸套某处相对于冷态的变形量为:
其中yp为活塞裙部顶端距上止点的轴向位移;
计算由活塞二阶运动导致的推力面油膜厚度变化量:
其中,a为活塞重心相对于活塞裙顶部的轴向位置,b为活塞销中心相对于活塞重心是轴向位置;
得到裙部与缸套之间的油膜厚度为:
5.根据权利要求1所述的一种考虑活塞与缸套热变形的裙部润滑特性仿真方法,其特征是:步骤四具体为:
根据已知的当前时刻缸内气体力和活塞环摩擦力以及与油膜和微凸体相关的力及力矩,得到连杆作用于活塞销的力,再由活塞销与连杆之间的摩擦系数得到连杆与活塞销之间的摩擦力矩,利用建立的活塞裙部动力学方程计算当前时刻的活塞销横向运动加速度和活塞摆动加速度动力学方程为:
其中,Mpin=μpinFL,μ为摩擦系数,下角标f代表与摩擦力相关的力和力矩,下角标c代表与微凸体接触相关的力和力矩,下角标j代表与角点接触相关的力和力矩,下角标h代表与油膜动压润滑相关的力和力矩,G代表与缸内气体相关的力和力矩,L代表连杆作用于活塞销相关的力和力矩,pin代表和活塞销相关的力和力矩,Ip为活塞绕重心的转动惯量;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310176497.1A CN116186933A (zh) | 2023-02-28 | 2023-02-28 | 一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310176497.1A CN116186933A (zh) | 2023-02-28 | 2023-02-28 | 一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116186933A true CN116186933A (zh) | 2023-05-30 |
Family
ID=86436313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310176497.1A Pending CN116186933A (zh) | 2023-02-28 | 2023-02-28 | 一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116186933A (zh) |
-
2023
- 2023-02-28 CN CN202310176497.1A patent/CN116186933A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Singh et al. | Design, analysis and optimization of three aluminium piston alloys using FEA | |
CN116186933A (zh) | 一种考虑活塞与缸套热变形的裙部润滑特性仿真方法 | |
Kumar et al. | Dynamic analysis of Bajaj Pulsar 150cc connecting rod using ANSYS 14.0 | |
Gai et al. | Optimization of Forged 42CrMo4 Steel Piston Pin Hole Profile Using Finite Element Method | |
MAHROOGI et al. | Dynamic-analysis-of-piston-secondary-motion-using FEM method | |
Bifeng et al. | Friction and wear performance of double-bump design of piston skirt main thrust side | |
Qin et al. | Tribological characteristics of piston rings in a single-piston hydraulic free-piston engine | |
CN112836445A (zh) | 一种用于求解滚动轴承的动态接触混合润滑参数的方法 | |
Abramchuk et al. | Prospects of using steel pistons in transport diesel engines | |
Borysevych et al. | Сurrent state of development of advanced forgings processes | |
Boon et al. | Predicting frictional losses generated by piston connecting-rod big end journal bearing for an internal combustion engine | |
Ye et al. | Study on the Lubrication Performances of Crankshaft Main Bearings in a 16V Marine Diesel Engine | |
CN116415462B (zh) | 基于浮动衬套的双层油膜润滑分析方法及系统 | |
Wable et al. | A Review on Design and Optimization of connecting rod | |
HariPriya et al. | Materialized Optimization of Connecting Rod for Four Stroke Single Cylinder Engine | |
Tripathi et al. | A review on fatigue behaviour of connecting rod made of Austempered Ductile Iron | |
UOI et al. | Elastohydrodynamic Lubrication Analysis of Journal Bearing | |
Petrovic | Fracture Conditions of Piston Rod Beam in Internal Combustion Engine Due to Damage | |
Kumari et al. | Research Article Failure Modes and Empirical Relations to Design Piston Pins for IC Engine: Analysis of piston pin | |
Joshi | Design, Analysis and Performance Testing of a Diesel Engine as a Portable Electrical Generator | |
Kayode et al. | Modelling and Optimization of Torsional Stiffness of a Crankshaft | |
Guo | Fatigue Strength Detection Method for Crankshaft of Vehicle Internal Combustion Engine Based on Dynamics | |
Tiwari et al. | A Review on Static structural behaviour of Piston | |
Kaya et al. | Optimum connecting rod design for diesel engines | |
Prakash et al. | Optimization of Connecting Rod for Weight Reduction by Finite Element Analysis |
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 |