CN114722526A - 一种高转速条件下带式cvt钢带环疲劳寿命计算方法 - Google Patents
一种高转速条件下带式cvt钢带环疲劳寿命计算方法 Download PDFInfo
- Publication number
- CN114722526A CN114722526A CN202210311001.2A CN202210311001A CN114722526A CN 114722526 A CN114722526 A CN 114722526A CN 202210311001 A CN202210311001 A CN 202210311001A CN 114722526 A CN114722526 A CN 114722526A
- Authority
- CN
- China
- Prior art keywords
- speed ratio
- speed
- life
- under
- belt ring
- 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
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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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)
- Geometry (AREA)
- Theoretical Computer Science (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)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种高转速条件下带式CVT钢带环疲劳寿命计算方法,包括:做带式CVT瞬态动力学分析;用等效法得到钢带环微元段运转一周的载荷谱,然后计算带环在此指定条件下的疲劳寿命;分别考虑固定速比下,转速、转矩变化以及速比变化两类情况,使用与前两步相同的寿命计算方法并结合数据插值法得到CVT带环理论寿命图像;考虑一般工况,首先用统计方法得到一般工况速比、输入转矩、输入转速服从的概率分布;再考虑CVT理论寿命图像与输入条件的概率计算定速比时带环寿命;最后利用不同速比带环的寿命与速比服从的概率分布得到带环一般工况的实际寿命。本发明从理论角度给出了钢带环疲劳寿命的计算方法,能够准确的估算出钢带环的疲劳寿命。
Description
技术领域
本发明涉及一种高转速条件下金属带式CVT技术研究领域,具体为一种高转速条件下金属钢带环疲劳寿命计算方法。
背景技术
随着社会的不断发展,汽车已经成为人们出行必备的交通工具,但是随之而来的是大量汽车尾气造成环境污染日益严重的棘手问题,世界各国都极为重视环境问题。电动汽车以电池作为动力源,是一种排放无污染的交通工具。随着国家不断严控汽车排放指标参数,电动汽车有望在未来代替传统汽车成为人们出行的主要交通工具。研究表明,为电动汽车配备变速器能极大改善电动汽车的动力性和经济性,增加其续航能力。在各类变速器中,无级变速器CVT配备电动汽车有着近乎理想的表现,它除了可以极大改善电动汽车的动力性和经济性之外,且相较于传统变速器,还具备结构简单,体积小,换挡平顺性好等独有优点。随着电动汽车的广泛使用,CVT无级变速器有望在未来成为使用最广泛的变速器。
但实验表明,CVT在高转速条件下可靠性会大幅度降低,金属钢带最外层钢带环极易发生断裂,导致CVT失效,失效时间比正常损坏大大提前,极易可能是由于带环的疲劳失效所致。但碍于具体的实验条件和CVT的特殊结构,无法通过实验去获得带环高转速下的实际受载情况,且CVT钢带环的整圈带环上均承受载荷作用,同一时刻载荷作用数量繁多。由于带环的不断旋转,载荷作用点的位置也在不断发生变化,这就使得带环疲劳寿命难以计算。为了能够提高高转速条件下金属带的可靠性,避免疲劳断裂的发生,提出一套可行的方案准确计算高转速条件下金属钢带环的疲劳寿命具有重要工程价值。
发明内容
本发明要解决的技术问题提供一种高转速条件下CVT钢带环疲劳寿命的计算方法。便于企业人员开发设计出性能更加优良的变速器。提高CVT的使用寿命以及扩宽CVT 的应用领域。
为实现上述目的,本发明采用如下技术方案:
一种高转速条件下CVT钢带环疲劳寿命的计算方法,包括如下步骤:
步骤1:做带式CVT瞬态动力学分析;
步骤2:用等效法提取钢带环微元段运转一周的载荷谱,然后求出钢带环在此指定条件下运行的疲劳寿命;
步骤3:分别考虑固定速比下,转速、转矩变化以及速比变化两类情况,使用与前两步相同的寿命计算方法并结合数据插值法得到CVT带环理论寿命图像;
步骤4:考虑一般工况,首先用统计方法得到一般工况的速比、输入转矩、输入转速服从的概率分布;
步骤5:再结合CVT理论寿命图像与输入条件的概率计算定速比时带环寿命;
步骤6:最后利用不同速比带环的寿命与速比服从的概率分布得到带环一般工况的实际寿命;
优选的,步骤1中,利用solid works建模软件建立某一速比时CVT装配体的三维模型,导入ANSYS WORKBENCH做CVT装配体的瞬态动力学分析;
优选的,步骤2中,用等效法得到带环微元段运转一周的载荷谱,并计算带环在此指定条件下的疲劳寿命的步骤包括:
步骤2.1:在瞬态动力学分析结果中提取整个CVT中所有金属片鞍面的接触载荷值,把CVT钢带环看作为与金属片数量相同的带环微元段柔性连接而成,以主动轮入口处第一片金属片对应的带环微元段为研究对象,把与之相邻的、顺时针旋转的下一片金属片等效看作是第一片金属片在下一个时刻的运动状态,循环一周,按顺序排列金属片鞍面载荷值,将之看作是带环微元段运转一周的受力变化情况,得到带环微元段的载荷谱,并记为鞍面载荷矩阵;
步骤2.2:将带环稳定运动一周的总时间划分为与金属片数段相等时间序列,将带环微元段鞍面载荷矩阵转换成时间与应力的关系,得到时序载荷谱;
步骤2.3:使用雨流计数法将时序载荷谱转换成计算所使用的应力载荷谱,得到在不同等级幅值的应力作用下实际作用的循环次数Ai;
步骤2.4:根据材料的S-N疲劳寿命曲线得到不同等级幅值的应力作用下对应的钢带环理论寿命次数Ni;
优选的,步骤3中,分别考虑固定速比时转速、转矩变化以及速比变化两类情况,使用步骤2寿命计算方法并结合数据插值法得到CVT带环理论寿命图像步骤包括:
步骤3.1:首先在固定速比下时,确定这个速比下,车速的范围,反推变速器变速范围,再根据发动机转速和转矩的关系,进而获得变速器输入转矩变化范围;
步骤3.2:在定速比、定转速下,取转矩范围两端点以及在转矩范围内再选几个具有代表性的输入转矩,计算方法不变,分别计算获得寿命对应输入条件下钢带环的疲劳寿命,得到的几个数据点,再运用牛顿插值公式拟合得到转矩与寿命的关系图像;
步骤3.3:依旧保持速比不变,更换一次转速,再选择与步骤3.2中同样的几个输入转矩点,计算方法不变,分别计算获得钢带环在对应输入条件下的疲劳寿命,将所得到的一组数据点,运用牛顿插值公式拟合得到转矩和寿命的关系图;
步骤3.4:重复上述过程几次,取得一些数据点,利用这些个数据点,通过描点作图得到定速比下,输入转速、输入转矩与寿命的关系图像;
步骤3.4:再考虑不同速比下,在速比范围内重新选一个速比,重新建立装配体模型,在ANSYS WORKBENCH里面做瞬态动力学分析,仿真结束后,使用同样的提取方法,提取得到鞍面载荷矩阵,按照上述步骤计算得到对应速比下的转速-转矩-寿命的关系图形,同样通过使用牛顿插值公式可以得到一系列连续点的疲劳寿命值,再改变几次速比,重新建立装配体,重复上述工作,最后描点作图,最终得到速比-转速-转矩 -寿命的理论关系图像。
优选的,步骤4中,考虑一般工况,首先用统计方法得到一般工况速比、输入转矩、输入转速服从的概率分布的步骤包括:
步骤4.1:在一般工况下,采集整个循环过程中,速比随时间的变化情况;不同速比下提取输入转矩随时间的变化情况,输入转速随时间的变化情况;
步骤4.2:依据此数据为前提,把输入转矩时间历程图、输入转速时间历程图、速比时间历程图转换成频度直方图;
步骤4.3:根据统计学里的中心极限定理,当数据量足够大时,任何变量的分布都会近似服从正态分布,对于变速器服役的工况而言,中间档位无疑是使用最多的,低速档和高速档使用均为较少频率使用的档位,同理汽车一般持续保持高速与低速行驶的频率较低,踩油门一般也踩的适中位置,因此同样认定速比、输入转速、输入转矩三个变量服从正态分布,使用正态分布函数拟合直方图,计算出三个正态分布函数变量的均值与标准差,得到正太概率分布函数表达式,进而求得汽车在实际行驶工况时,指定速比使用的概率以及指定速比下,某输入转矩、输入转速使用到的概率,求出的概率分布即是用于表征此时金属带的服役情况的数学函数。
优选的,步骤5中,结合CVT理论寿命图像与输入条件的概率计算定速比时带环寿命的步骤包括:
步骤5.1:步骤4得到三个变量的概率分布后,即明确了一般工况下金属带的服役条件,为了使计算操作简便、实际可操作性强以及结算结果保持较高准确性,使用将连续函数用多段小区间进行逼近的方式进行实际计算,首先将三个变量统计时划分的小区间对应的区间中值求出,然后把此小区间范围代入分布函数求出概率,用这个区间中值代表对应小区间的概率,区间越细,拟合效果越好;
步骤5.2:有了上述概率数据后,再把划分的小区间对应的区间中值的输入转矩、输入转速在步骤3中所建立的钢带环理论寿命图像中查图得到理论寿命值,然后再乘以对应区间中值在这个服役条件下的使用到的概率,计算出金属带在某一定速比下的实际寿命值。
优选的,步骤6中,由步骤5中求出的多个速比下的钢带环的实际寿命值,再结合速比概率分布函数求出对应速比的概率值,则带环在一般工况下的实际寿命就等于对应速比的实际寿命乘以对应速比的概率并相互累加,所用到的CVT钢带环寿命的实际值计算公式为:
优选的,步骤4中,频度直方图的建立方法为:
三个直方图分别以输入转矩M、输入转速n、速比i为横坐标,以转矩频度HMi与输入转矩间隔ΔM之比转速频度Hni与输入转速间隔Δn之比速比频度Hii与速比间隔Δi之比为纵坐标,频度HMi,Hni,Hii反映的是每一级划分的区间内目标值出现次数的百分比,因此直方图反映了对应指标出现的次数比例。
QMi、Qni、Qii为每一级参数出现的次数。
优选的,步骤4中,拟合频度直方图使用的正态分布函数通过以下方式求得,将直方图的纵坐标与横坐标的关系可以用正态分布函数表示:.
式中:x1,x2,x3为采集的数据中某输入转矩、某输入转速、某速比值;P1(x1)、P2(x2)、P3(x3)分别为某一速比下使用某输入转矩的概率值、使用某转速的概率值、使用某速比的概率值;HMi,Hni,Hii:转矩频度、转速频度、速比频度;ΔM,Δn,Δi:扭矩间隔、转速间隔、速比间隔;Mi、nj、ik:每一级的输入转矩、输入转速、速比值;扭矩的均值,转速均值,速比均值,标准差,σx2:标准差,σx3:标准差,e:自然对数的底。
优选的,步骤5中,定速比下钢带环疲劳寿命通过以下公式计算:
本发明具有如下有益效果:
本发明提供的一种高转速条件下CVT带坏疲劳寿命计算方法,该方法能够尽可能准确的计算出CVT在一般工况下的疲劳寿命。对于相关从业人员设计出高可靠性的带式CVT提供研究方法。以前的钢带环疲劳寿命研究中只能通过实验进行,而随着电动汽车的兴起,转速条件变得越来越高时,往往难以达到实验条件,而通过此理论可以准确计算出钢带环的疲劳寿命。
附图说明
图1高转速条件下钢带环疲劳寿命计算方法;
图2为定速比下钢带环理论疲劳寿命计算方法;
图3为不同速比时钢带环理论疲劳寿命计算方法。
具体实施方式
下面结合附图对本发明做进一步的详细说明,本发明采用计算方案如图1所示,具体步骤如下:
步骤1:做带式CVT瞬态动力学分析;
步骤2:用等效法提取钢带环微元段运转一周的载荷谱,然后计算钢带环在此指定条件下运行的疲劳寿命;
步骤3:分别考虑固定速比下,转速、转矩变化以及速比变化两类情况,使用与上步相同的寿命计算方法并结合数据插值法得到CVT带环理论寿命图像;
步骤4:考虑一般工况,首先用统计方法得到一般工况速比、输入转矩、输入转速服从的概率分布;
步骤5:再结合CVT理论寿命图像与输入条件的概率计算定速比时带环寿命;
步骤6:最后利用不同速比带环的寿命与速比服从的概率分布得到带环一般工况的实际寿命。
本发明中,为了精准计算钢带环的疲劳寿命,总体方案可概括如下,首先考虑了钢带环的理论工作寿命值,再考虑了钢带环在一般工况下工作的随机因素,最后联合两者的数据给出了钢带环一般工况下疲劳寿命的计算方式。具体为:第一,在理论疲劳寿命分析阶段,也就是指求解钢带环在固定输入条件下稳定运行直至失效的疲劳寿命时。例如给定一个速比、一个输入转矩和一个输入转速,将此条件看作为钢带环的工作条件,使用本发明提到的计算疲劳寿命方式,得到此固定条件下钢带环的理论寿命值。然后再分别在定速比下,改变输入转速、输入转矩和改变速比两种情况下去选取不同点计算多个钢带环在对应固定条件下的理论寿命,最后结合数据插值方法以得到钢带环的速比、输入转矩、输入转速与寿命的理论寿命图像,以上即理论寿命分析阶段的任务。第二,在一般工况阶段,考虑带环在一般工况下工作的疲劳寿命。本发明将带环在一般工况下的服役的条件用统计学中的分布函数表征出来,提出了考虑更一般的情况下如何去计算钢带环的疲劳寿命值。并且为了提升本发明实际操作可行性,还提出使用划分小区间的方式对表征钢带环服役条件的分布函数进行逼近,以离散的方式进行钢带环疲劳寿命计算,并给出了计算钢带环在定速比下的疲劳寿命公式以及在一般工况下钢带环疲劳寿命公式,总计两个新公式。该方法在计算钢带环疲劳寿命时准确率高,并且可操行强,对钢带环的设计研究具有较大的帮助。
具体实施过程中,步骤1中通过如下方法做带式CVT瞬态动力学分析:使用solidworks三维建模软件,建立速比为iL1的CVT模型(L代表建立带环理论寿命时用到的符号;iL1符号含义:指在建立理论寿命时使用到的某个速比,实际操作中,自己可在速比范围中选取任意速比建模,此符号取法是为了与后文计算实际寿命时使用的速比iS1区分开。后文提到L和S均为理论与实际的含义简称),然后保存为x_t格式,导入 ANSYS WORKBENCH中设置材料、设置接触、设置运动副、划分网格、边界条件设置、分析设置、最后求解做带式CVT的瞬态动力学分析;
具体实施过程中,步骤2中,用等效法得到带环微元段运转一周的载荷谱,并计算带环在此指定条件下的疲劳寿命的步骤包括:
步骤2.1:给定输入转速nLb1(L代表建立带环理论寿命时使用的转速;b代表此时的速比;1代表转速第一次取一个点,下文符号均为此含义),输入转矩MLbk(b代表此时的速比;k代表转速转矩取的某个点位置,下文符号均为此含义),瞬态动力学仿真结束后,提取得到整个金属带所有金属片的鞍面接触应力值为σ1,σ2,σ3....σn。其中n为金属片片数。以主动轮入口处第一片金属片对应的带环微元段为研究对象(以其它带环微元段为研究对象同理),由于整个带环上分布的金属片是几乎完全一样的,因而把与之相邻的、顺时针旋转的下一片金属片等效看作是第一片金属片在下一个时刻的运动状态,进而我们将第二片金属片鞍面的接触载荷看作是第一片金属片在下一个时刻的受载情况,如此顺时针循环下去,再次到回到主动轮第一片时,停止循环。对于带环来说,直接的施力对象只有金属片,因此可以把提取到的金属片鞍面载荷看作是带环微元段运转一周的受力变化情况,将得到的金属片鞍面应力值组建为鞍面载荷矩阵B1,如表1所示:
表1鞍面载荷矩阵B1
金属片鞍面载荷 | σ<sub>1</sub> | σ<sub>2</sub> | σ<sub>3</sub> | ……… | σ<sub>n</sub> |
步骤2.2:将带环稳定运动一周的总时间t划分为金属片数段时间序列t1-tn。将鞍面载荷矩阵B1转换为时间与应力关系图表示出来,得到时序载荷谱,如表2所示:
表2时序载荷谱
时间 | t<sub>1</sub> | t<sub>2</sub> | t<sub>3</sub> | ……… | t<sub>n</sub> |
金属片鞍面载荷 | σ<sub>1</sub> | σ<sub>2</sub> | σ<sub>3</sub> | ……… | σ<sub>n</sub> |
步骤2.3:有了时序载荷谱之后,使用雨流计数法对获得的随机交变应力的幅值进行统计。得到计算使用的应力载荷谱,具体如下:
令Smax1=σmax1,Smax2=σmax2,.......,Smaxn=σmaxn
式中:σmaxi与Smax1为雨流计数法统计过后的不同等级幅值的应力循环的最大应力。雨流计数法没有考虑应力副作用于零件的先后,所以这里没有了时间的序列。表3中Ai为再带环运转一周的过程中,零件在应力Smaxi作用下,实际作用的循环次数:
表3不同等级幅值的应力smaxi作用下实际作用的次数
步骤2.4:然后再根据材料的S-N疲劳寿命曲线得到应力smaxi作用下的理论寿命次数,如表4所示:
表4不同等级幅值的应力smaxi作用下的材料的极限循环次数
对于指标TLbk1含义:用字母L代表建立带环理论寿命时的代号,字母代号b表示指定速比下,第三个数字k取值代表在输入转矩范围内取点(从1到n中随机取)的变化,第四个数字1取值代表转速取点(从1到n中随机取)的变化次数。
具体实施过程中,步骤3中,如图2、图3所示的计算原理,分别考虑在固定速比下,转速、转矩变化以及速比变化两类情况,使用与前两步相同的寿命计算方法并结合数据插值法得到CVT带环理论寿命图像步骤包括:
步骤3.1:首先考虑定速比iL1时,确定这个速比下,车速的范围V1-Vn,反推变速器变速范围nLb1-nLbn,再根据发动机转速和转矩的关系,进而获得输入转矩变化范围MLb1-MLbn;
步骤3.2:在定速比,定转速下,取转矩范围两端点MLb1和MLbn以及在转矩范围内再选1个具有代表性的输入转矩点MLbj,计算方法不变,可以计算获得寿命TLb11、 TLbj1、TLbn1(对于指标TLb11含义,用字母代号b表示指定速比下,第二个数字1取值代表在输入转矩范围内取点(从1-n范围中取点)的变化,第三个数字1取值代表转速取点 (从1-n中范围中取点)的变化次数)。加上最初的取得一个转矩点,总共在此速比下转矩范围内取四个输入转矩点,最后得到四个坐标点(MLb1,TLb11),(MLbk,TLbk1),(MLbj, TLbj1),(MLbn,TLbn1),将四个坐标点再运用牛顿插值公式得到转矩与寿命的关系图像;
步骤3.3:依旧在速比不变,换一个转速nLb2,再选择步骤3.2中同样的四个位置的输入转矩点,计算方法不变,可以计算获得寿命TLb12、TLbk2,TLbj2,TLbn2。将四个坐标点(MLb1,TLb12),(MLbk,TLbk2),(MLbj,TLbj2),(MLbn,TLbn2),运用牛顿插值公式得到转矩和寿命的关系;
步骤3.4:重复上述过程四次,即再取四个转速,总共取五个转速,并全部计算出寿命值,这样可以最终获得一些数据点,所得数据如表5所示(表格中b代表在速比 iL1下用到的参数):
表5速比为iL1时,钢带环在对应输入条件稳定运行直至失效的疲劳寿命
输入转速n | 输入转矩M | 此条件稳定运行下去的疲劳寿命 |
n<sub>Lb1</sub> | M<sub>Lb1</sub> | T<sub>Lb11</sub> |
n<sub>Lb1</sub> | M<sub>Lbk</sub> | T<sub>Lbk1</sub> |
n<sub>Lb1</sub> | M<sub>Lbj</sub> | T<sub>Lbj1</sub> |
n<sub>Lb1</sub> | M<sub>Lbn</sub> | T<sub>Lbn1</sub> |
n<sub>Lb2</sub> | M<sub>Lb1</sub> | T<sub>Lb12</sub> |
n<sub>Lb2</sub> | M<sub>Lbk</sub> | T<sub>Lbk2</sub> |
n<sub>Lb2</sub> | M<sub>Lbj</sub> | T<sub>Lbj2</sub> |
n<sub>Lb2</sub> | M<sub>Lbn</sub> | T<sub>Lbn2</sub> |
n<sub>Lb3</sub> | M<sub>Lb1</sub> | T<sub>Lb13</sub> |
n<sub>Lb3</sub> | M<sub>Lbk</sub> | T<sub>Lbk3</sub> |
……… | ……… | ………… |
n<sub>Lb5</sub> | M<sub>Lbn</sub> | T<sub>Lbn5</sub> |
再通过使用描点作图可以得到定速比下,输入转速、输入转矩与寿命的关系图像;
步骤3.5:如图3所示的计算流程,在速比范围内重新选一个速比iL2,重新建立装配体模型,在ANSYS WORKBENCH里面做瞬态动力学分析,仿真结束后,使用同样的提取方法,提取得到鞍面载荷矩阵B2,按照计算上一速比的寿命步骤可以得到对应速比下的转速,转矩和寿命的关系图形,计算所得数据如表6所示((表格中c代表在速比iL2下用到的参数):
表6速比为iL2时,钢带环在对应输入条件稳定运行的疲劳寿命
输入转速n | 输入转矩M | 此条件稳定运行下去的疲劳寿命 |
n<sub>Lc1</sub> | M<sub>Lc1</sub> | T<sub>Lc11</sub> |
n<sub>Lc1</sub> | M<sub>Lck</sub> | T<sub>Lck1</sub> |
n<sub>Lc1</sub> | M<sub>Lcj</sub> | T<sub>Lcj1</sub> |
n<sub>Lc1</sub> | M<sub>Lcn</sub> | T<sub>Lcn1</sub> |
n<sub>Lc2</sub> | M<sub>Lc1</sub> | T<sub>Lc12</sub> |
……… | ……… | ……… |
n<sub>Lc5</sub> | M<sub>Lcn</sub> | T<sub>Lcn5</sub> |
同样通过使用牛顿插值公式可以得到一系列连续点的疲劳寿命值。类似再改变三次速比iL3、iL4、iL5,重新建立装配体,重复上述工作,然后通过描点作图,最终可以得到速比-转速-转矩-寿命的理论关系图像。上述步骤所运用的牛顿插值公式如下:
P(x)=f(x0)+f[x0,x1](x-x0)+f[x0,…,xn](x-x0)…(x-xn)
式中:P(x)为利用已知数据点得到的插值函数;f[xi,xj]为一阶差商值,xi,xj为已知数据点的横坐标;f[x0,x1,…,xn]为n阶差商值。
本实施例中,当我们掌握了理论寿命图像数据之后,只要能知道一个变速器在某一速比下,其输入转矩和输入转速的值,就可以通过所建立的理论寿命图像查到此变速器在这个工况稳定运行时的理论疲劳寿命。例随意取速比为i0时,转速为n0,转矩为M0,可以通过图像查出此时寿命为T0。
具体实施过程中,步骤4中,考虑一般工况,首先用统计方法得到一般工况速比、输入转矩、输入转速服从的概率分布的步骤包括:
步骤4.1:一般工况下,采集整个循环过程中,速比随时间的变化情况;不同速比下提取输入转矩随时间的变化情况,输入转速随时间的变化情况;
步骤4.2:依据采集的数据为前提,把输入转矩时间历程图、输入转速时间历程图、速比时间历程图转换成频度直方图。三个直方图分别以输入转矩M、输入转速n、速比 i为横坐标,将输入转矩M、输入转速n、速比i划分为相等的份数,例如从0-5划分一份,5-10划分一份等。然后令转矩频度HMi与输入转矩间隔ΔM之比转速频度Hni与输入转速间隔Δn之比速比频度Hii与速比间隔Δi之比为纵坐标,做出图像。频度HMi,Hni,Hii反映的是每一级设定目标范围内值出现次数的百分比,因此直方图反映了对应指标出现的次数比例。
QMi、Qni、Qii为每一级参数出现的次数;
步骤4.3:根据统计学里的中心极限定理,当数据量足够大时,任何变量的分布都会近似服从正态分布,对于变速器服役的工况而言,中间档位无疑是使用最多的,低速档和高速档使用均为较少频率使用的档位,同理汽车一般持续保持高速与低速行驶的频率较低,踩油门一般也踩的适中位置,因此同样认定速比、输入转速、输入转矩三个变量服从正态分布,使用正态分布函数拟合直方图,将纵坐标与横坐标的关系可以用正态分布函数表示:
式中:x1,x2,x3为采集的数据中某输入转矩、某输入转速、某速比的值;P1(x1)、P2(x2)、 P3(x3)分别为某一速比下使用某输入转矩的概率值、使用某转速的概率值、使用某速比的概率值;HMi,Hni,Hii:转矩频度、转速频度、速比频度;ΔM,Δn,Δi:扭矩间隔、转速间隔、速比间隔;Mi、nj、ik:每一级的输入转矩、输入转速、速比值;扭矩的均值,转速均值,速比均值,标准差,σx2:标准差,σx3:标准差,e:自然对数的底。
本实施例中,根据上述变量的均值与标准差的计算公式,计算出均值和标准差,带入分布函数表达式得到概率分布函数,用以代替直方图,进而求得汽车在实际行驶工况时,指定速比使用的概率以及指定速比下,某输入转矩、输入转速使用到的概率。求出的概率分布为表征此时金属带的服役情况的数学函数。
具体实施中,步骤5中,结合CVT理论寿命图像与输入条件的概率计算定速比时带环寿命的步骤包括:
步骤5.1:步骤4得到以上三个变量的概率分布过后,即明确了一般工况下金属带的服役条件,为了使计算操作简便、实际可操性强以及结算结果保持较高准确性,使用将连续函数用多段小区间进行逼近的方式进行实际计算。首先将三个变量统计时划分的小区间对应的区间中值求出,然后把此小区间范围代入分布函数求出概率,用这个区间中值代表对应小区间的概率。区间越细,拟合效果就越好。下面提到的表格中所用的每一个速比值(例如is1)以及输入转速(例如nsb1)、输入转矩(例如Msb1)均为划分的小区间的区间中值,首先分别计算出不同速比下,输入转矩与输入转速的概率值;
具体数据计算结果如表7所示(概率用符号P表示;is1代表计算实际寿命时使用的速比):
表7速比is1(表格中用b代指此速比大小下得到的概率值)
转速小区间中值 | 概率 | 转矩小区间中值 | 概率 |
n<sub>sb1</sub> | P<sub>sbn1</sub> | M<sub>sb1</sub> | P<sub>sbM1</sub> |
n<sub>sb1</sub> | P<sub>sbn1</sub> | M<sub>sb2</sub> | P<sub>sbM2</sub> |
n<sub>sb1</sub> | P<sub>sbn1</sub> | M<sub>sb3</sub> | P<sub>sbM3</sub> |
n<sub>sb1</sub> | P<sub>sbn1</sub> | M<sub>sb4</sub> | P<sub>sbM4</sub> |
n<sub>sb2</sub> | P<sub>sbn2</sub> | M<sub>sb1</sub> | P<sub>sbM1</sub> |
n<sub>sb2</sub> | P<sub>sbn2</sub> | M<sub>sb2</sub> | P<sub>sbM2</sub> |
n<sub>sb2</sub> | P<sub>sbn2</sub> | M<sub>sb3</sub> | P<sub>sbM3</sub> |
…… | …… | …… | …… |
n<sub>sb5</sub> | P<sub>sbn5</sub> | M<sub>sb4</sub> | P<sub>sbM4</sub> |
表7中nsbi,Msbi为is1速比下划分的小区间中值,其中在这个速比下需满足:
Psbn1+Psbn2+Psbn3+Psbn4+Psbn5=1
PsbM1+PsbM2+PsbM3+PsbM4=1
又当速比is2时,具体数据计算结果如表8所示:
表8速比is2(表格中用c代指此速比大小下得到的概率值)
转速小区间中值 | 概率 | 转矩小区间中值 | 概率 |
n<sub>sc1</sub> | P<sub>scn1</sub> | M<sub>sc1</sub> | P<sub>scM1</sub> |
n<sub>sc1</sub> | P<sub>scn1</sub> | M<sub>sc2</sub> | P<sub>scM2</sub> |
n<sub>sc1</sub> | P<sub>scn1</sub> | M<sub>sc3</sub> | P<sub>scM3</sub> |
n<sub>sc1</sub> | P<sub>scn1</sub> | M<sub>sc4</sub> | P<sub>scM4</sub> |
n<sub>sc2</sub> | P<sub>scn2</sub> | M<sub>sc1</sub> | P<sub>scM1</sub> |
n<sub>sc2</sub> | P<sub>scn2</sub> | M<sub>sc2</sub> | P<sub>scM2</sub> |
…… | …… | …… | …… |
n<sub>sc5</sub> | P<sub>scn5</sub> | M<sub>sc4</sub> | P<sub>scM4</sub> |
表8中nsci,Msci为is2速比下划分的小区间中值,其中在这个速比下需满足:
Pscn1+Pscn2+Pscn3+Pscn4+Pscn5=1
PscM1+PscM2+PscM3+PscM4=1
同理再求出is3、is4、is5速比下输入转矩、输入转速被使用到的概率值;
步骤5.2:有了上述概率数据后,再把划分的小区间对应的区间中值的输入转矩、输入转速在所建立的钢带环理论寿命图像中查表得到理论寿命值,再乘以对应区间中值在这个车况下的使用到的概率,即可计算出金属带在某一速比下的实际寿命值,通过使用以下计算公式计算:
式中累加的次数取决于划分小区间的个数;Psbni、Pscni…为某速比下,某一转速的概率; PsbMi、PscMi...为某速比下,某一转矩的概率;Tsb11、Tsb21...为查寻理论图像所得一定输入转矩和转速下的疲劳寿命;定速比下的带式CVT的疲劳寿命;
本实施例中,步骤5.2中,具体操作如下:通过步骤3建立的理论图像,查图得到速比为is1、is2、is3、is4、is5时的寿命值。首先当速比为is1时,用对应输入转速、输入转矩的小区间中值通过查图得到带环的理论疲劳寿命值如表9所示。
表9速比为is1时,输入转速、输入转矩小区间中值通过查图得到的带环疲劳寿命值
再结合步骤5.1中得到的对应小区间中值的输入概率值计算得:
具体实施过程中,步骤6中,由步骤5中求出的多个速比下的钢带环的实际寿命值,再结合速比概率分布函数求出对应速比的概率值,求得带环在一般工况下的实际寿命就等于对应速比的实际寿命乘以对应速比的概率并相互累加,所用到的CVT钢带环寿命的实际值计算公式为:
本实施例中,步骤6中,具体操作如下:根据步骤5求得的取得五个速比点处的疲劳寿命,再结合速比分布函数求得速比小区间中值的概率值,如表10所示:
表10速比本身服从的概率
其中速比概率满足
P1+P2+P3+P4+P5=1
最后带环在所推定的工况下计算寿命为:
本发明是对高转速条件下钢带环疲劳寿命进行计算。由于钢带环的受载形式以及运动方式的特殊性,同一时刻金属钢带受载作用点繁多,且由于带环的不断旋转,力的作用点也在不断变化,因此用理论去计算钢带环的疲劳寿命显得十分困难,这也是首次针对钢带环疲劳寿命提出的理论分析办法,具有一定的新颖性。在以往的研究中,钢带环疲劳寿命通常采用实验研究,但这仅仅限于较低的实验要求,并且计算出的寿命还是在指定的工况下的疲劳寿命,不具有一般性,而本发明提出了一种计算带环在一般工况下工作时疲劳寿命的计算方法,更符合实际条件,能够尽可能准确的估算出钢带环的疲劳寿命。
以上所述的仅是本发明的实施例,方案中公知的具体结构及特性等常识在此未作过多描述,所属领域普通技术人员知晓申请日或者优先权日之前发明所属技术领域所有的普通技术知识,能够获知该领域中所有的现有技术,并且具有应用该日期之前常规实验手段的能力,所属领域普通技术人员可以在本申请给出的启示下,结合自身能力完善并实施本方案,一些典型的公知结构或者公知方法不应当成为所属领域普通技术人员实施本申请的障碍。应当指出,对于本领域的技术人员来说,在不脱离本发明结构的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。
Claims (10)
1.一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,包括以下步骤:
步骤1:做带式CVT瞬态动力学分析;
步骤2:用等效法提取钢带环微元段运转一周的载荷谱,然后求出钢带环在此指定条件下运行的疲劳寿命;
步骤3:分别考虑固定速比下,转速、转矩变化以及速比变化两类情况,使用与前两步相同的寿命计算方法并结合数据插值法得到CVT带环理论寿命图像;
步骤4:考虑一般工况,首先用统计方法得到一般工况速比、输入转矩、输入转速服从的概率分布;
步骤5:再结合CVT理论寿命图像与输入条件的概率计算定速比时带环寿命;
步骤6:最后利用不同速比带环的寿命与速比服从的概率分布得到带环一般工况的实际寿命。
2.根据权力要求1所述的一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,对于步骤1,利用solid works建模软件建立某一速比时CVT装配体的三维模型,导入ANSYS做CVT装配体的瞬态动力学分析。
3.根据权力要求1所述的一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,步骤2中,用等效法得到带环微元段运转一周的载荷谱,并计算带环在此指定条件下的疲劳寿命的步骤包括:
步骤2.1:在瞬态动力学分析结果中提取整个CVT中所有金属片鞍面的接触应力值,把CVT钢带环看作为与金属片数量相同的带环微元段柔性连接而成,以主动轮入口处第一片金属片对应的带环微元段为研究对象,把与之相邻的、顺时针旋转的下一片金属片等效看作是第一片金属片在下一个时刻的运动状态,循环一周,按顺序排列金属片鞍面载荷值,将之看作是带环微元段运转一周的受力变化情况,得到带环微元段的载荷谱,并记为鞍面载荷矩阵;
步骤2.2:将带环稳定运动一周的总时间划分为与金属片数段相等时间序列,将带环微元段鞍面载荷矩阵转换成时间与应力的关系,得到时序载荷谱;
步骤2.3:使用雨流计数法将时序载荷谱转换成计算所使用的应力载荷谱,得到在不同等级幅值的应力作用下实际作用的循环次数Ai;
步骤2.4:根据材料的S-N疲劳寿命曲线得到不同等级幅值的应力作用下对应的钢带环理论寿命次数Ni;
4.根据权力要求1所述的一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,步骤3中,分别考虑固定速比下,转速、转矩变化以及速比变化两类情况,使用上步寿命计算方法并结合数据插值法得到CVT带环理论寿命图像步骤包括:
步骤3.1:首先在固定速比下时,确定这个速比下,车速的范围,反推变速器变速范围,再根据发动机转速和转矩的关系,进而获得变速器输入转矩变化范围;
步骤3.2:在定速比、定转速下,取转矩范围两端点以及在转矩范围内再选几个具有代表性的输入转矩,计算方法不变,分别计算获得对应输入条件下钢带环的疲劳寿命,从而得到的一些数据点,再运用牛顿插值公式拟合得到转矩与寿命的关系图像;
步骤3.3:依旧保持速比不变,更换一次转速,再选择与步骤3.2中同样的几个输入转矩点,计算方法不变,分别计算获得钢带环在对应输入条件下的疲劳寿命,将所得到的一组数据点,运用牛顿插值公式拟合得到转矩和寿命的关系图;
步骤3.4:重复上述过程几次,取得一些数据点,利用这些个数据点,再通过描点作图得到定速比下,输入转速、输入转矩与寿命的关系图像;
步骤3.5:再考虑不同速比下,在速比范围内重新选一个速比,重新建立装配体模型,在ANSYS WORKBENCH里面做瞬态动力学分析,仿真结束后,使用同样的提取方法,提取得到鞍面载荷矩阵,按照上述步骤计算得到对应速比下的转速-转矩-寿命的关系图形,同样通过使用牛顿插值公式可以得到一系列连续点的疲劳寿命值,再改变几次速比,重新建立装配体,重复上述工作,最后通过描点作图,最终得到速比-转速-转矩-寿命的理论关系图像。
5.根据权力要求1所述的一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,步骤4中,考虑一般工况,首先用统计方法得到一般工况速比、输入转矩、输入转速服从的概率分布的步骤包括:
步骤4.1:在一般工况下,采集随机的一个循环过程中,速比随时间的变化情况、不同速比下提取输入转矩随时间的变化情况、输入转速随时间的变化情况;
步骤4.2:依据采集的数据为前提,把输入转矩时间历程图、输入转速时间历程图、速比时间历程图转换成频度直方图;
步骤4.3:根据统计学里的中心极限定理,当数据量足够大时,任何变量的分布都会近似服从正态分布,对于变速器服役的工况而言,中间档位无疑是使用最多的,低速档和高速档使用均为较少频率使用的档位,同理汽车一般持续保持高速与低速行驶的频率较低,踩油门一般也踩的适中位置,因此同样认定速比、输入转速、输入转矩三个变量服从正态分布,使用正态分布函数拟合直方图,先计算出三个正态分布函数变量的均值与标准差,然后得到正态概率分布函数表达式,进而求得汽车在实际行驶工况时,指定的某速比使用的概率以及指定速比下,某输入转矩、输入转速使用到的概率,求出的概率分布即理解为表征此时金属带的服役情况的数学函数。
6.根据权力要求1所述的一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,步骤5中,结合CVT理论寿命图像与输入条件的概率计算定速比时带环寿命的步骤包括:
步骤5.1:步骤4得到三个变量的概率分布后,即明确了一般工况下金属带的服役条件,为了使计算操作简便、实际可操作性强以及结算结果保持较高准确性,使用将连续函数用多段小区间进行逼近的方式进行实际计算,首先将三个变量统计时划分的小区间对应的区间中值求出,然后把此小区间范围代入分布函数求出概率,用这个小区间中值代表对应小区间的概率,区间越细,拟合效果越好;
步骤5.2:有了上述概率数据后,再把划分的小区间对应的区间中值的输入转矩、输入转速在步骤3中所建立的钢带环理论寿命图像中通过查图得到理论寿命值,然后再乘以对应区间中值在这个服役条件下的使用到的概率,计算出金属带在某一定速比下的实际寿命值。
9.根据权利要求5所述的一种高转速条件下带式CVT钢带环疲劳寿命计算方法,其特征在于,步骤4中,拟合频度直方图使用的正态分布函数通过以下方式求得:
将直方图的纵坐标与横坐标的关系可以用正态分布函数表示:.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210311001.2A CN114722526A (zh) | 2022-03-28 | 2022-03-28 | 一种高转速条件下带式cvt钢带环疲劳寿命计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210311001.2A CN114722526A (zh) | 2022-03-28 | 2022-03-28 | 一种高转速条件下带式cvt钢带环疲劳寿命计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114722526A true CN114722526A (zh) | 2022-07-08 |
Family
ID=82240755
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210311001.2A Pending CN114722526A (zh) | 2022-03-28 | 2022-03-28 | 一种高转速条件下带式cvt钢带环疲劳寿命计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114722526A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115389196A (zh) * | 2022-10-26 | 2022-11-25 | 麦格纳动力总成(江西)有限公司 | 一种同步器耐久性试验方法及系统 |
-
2022
- 2022-03-28 CN CN202210311001.2A patent/CN114722526A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115389196A (zh) * | 2022-10-26 | 2022-11-25 | 麦格纳动力总成(江西)有限公司 | 一种同步器耐久性试验方法及系统 |
CN115389196B (zh) * | 2022-10-26 | 2023-03-10 | 麦格纳动力总成(江西)有限公司 | 一种同步器耐久性试验方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2023020628A1 (zh) | 一种基于传动效率的减速器工况及结构参数优选方法 | |
CN103678256B (zh) | 车用发动机载荷谱编制方法 | |
CN106126843A (zh) | 一种基于Matlab的Bladed风机载荷处理系统 | |
CN110550034A (zh) | 一种纯电动汽车两档amt综合换挡方法 | |
CN104102781A (zh) | 一种基于多目标优化的发动机性能优化方法 | |
CN109271702B (zh) | 工程车发动机与液力变矩器逆向匹配的优化方法 | |
CN114722526A (zh) | 一种高转速条件下带式cvt钢带环疲劳寿命计算方法 | |
You et al. | Shift strategy of a new continuously variable transmission based wheel loader | |
CN109000922A (zh) | 基于道路测试的综合传动装置台架耐久性试验方法 | |
CN110243591B (zh) | 一种变速器润滑油加注量的确定方法 | |
CN110843535B (zh) | 车辆动力匹配方法、装置、设备及存储介质 | |
Kęsy et al. | Construction optimization of hydrodynamic torque converter with application of genetic algorithm | |
CN106529161A (zh) | 一种基于火电机组运行数据确定升降负荷速率的方法 | |
Huang et al. | Optimal design for improving operation performance of electric construction machinery collaborative system: Method and application | |
CN111581603A (zh) | 载荷谱处理方法及装置 | |
CN117892459A (zh) | 一种弧齿锥齿轮副磨损故障预测及其啮合特性分析方法 | |
CN102867127A (zh) | 一种用于工程机械液力传动匹配的方法 | |
CN103246759A (zh) | 一种滑动轴承的优化设计方法 | |
CN110469661B (zh) | 一种基于cvt效率的动力性速比优化方法及系统 | |
Li et al. | Shift scheduling strategy development for parallel hybrid construction vehicles | |
CN113700848B (zh) | 一种用于车辆综合传动的功率损失分析方法和装置 | |
Roth et al. | Optimization-based Component Sizing Method for Electrified Heavy-Duty Powertrain Concepts | |
Zou et al. | Establish oriented operating terminals wheel loader duty cycle | |
Thai et al. | Optimization of the gear ratios for a vehicle manual transmission | |
CN108153950A (zh) | 一种发动机与液力变矩器的自动匹配方法 |
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 |