CN115098975A - 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 - Google Patents
一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 Download PDFInfo
- Publication number
- CN115098975A CN115098975A CN202210857752.4A CN202210857752A CN115098975A CN 115098975 A CN115098975 A CN 115098975A CN 202210857752 A CN202210857752 A CN 202210857752A CN 115098975 A CN115098975 A CN 115098975A
- Authority
- CN
- China
- Prior art keywords
- friction plate
- inner hub
- tooth
- representing
- gear teeth
- 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
- 238000004364 calculation method Methods 0.000 title claims description 22
- 238000000034 method Methods 0.000 claims abstract description 29
- 238000013016 damping Methods 0.000 claims abstract description 28
- 230000005484 gravity Effects 0.000 claims abstract description 25
- 238000004519 manufacturing process Methods 0.000 claims abstract description 13
- 238000007667 floating Methods 0.000 claims abstract description 10
- 230000001133 acceleration Effects 0.000 claims description 30
- 238000006073 displacement reaction Methods 0.000 claims description 29
- 238000005452 bending Methods 0.000 claims description 18
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 238000005381 potential energy Methods 0.000 claims description 9
- 238000000926 separation method Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000006008 O'Donnell synthesis reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000009415 formwork Methods 0.000 claims description 3
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 229910052731 fluorine Inorganic materials 0.000 description 2
- 125000001153 fluoro group Chemical group F* 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002360 preparation method Methods 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
- 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/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)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mechanical Operated Clutches (AREA)
Abstract
本发明公开一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,具体包括以下步骤:S1:获取待计算的浮动支撑摩擦片组件参数,所述摩擦片组件参数包括基本参数和随机齿形制造误差量,所述摩擦片组件包括摩擦片和内毂;S2、计算摩擦片与内毂单齿接触的综合刚度;S3、计算摩擦片在重力作用下引起的偏心距,以及每时刻摩擦片的时变偏心距和时变偏心相位;S4、判断摩擦片与内毂是否发生碰撞,确定内毂与摩擦片碰撞齿数;S5、计算摩擦片与内毂轮齿间碰撞的冲击反弹系数和迟滞阻尼系数;S6、计算摩擦片与内毂每个接触齿的碰撞力;S7、构建内毂与摩擦片平动‑扭转多自由度冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载。
Description
技术领域
本发明涉及摩擦片技术领域,特别涉及一种考虑随机齿向误差的离合器摩擦片冲击力计算方法。
背景技术
摩擦片与内毂轮齿在加工过程中不可避免的会使各轮齿间的齿侧间隙呈随机性分布,受摩擦片重力作用的影响摩擦片与内毂会存在一定的偏心距。摩擦片组件运转过程中,由于变速机构的扭转振动等引起内毂转速表现为高频特征的非平稳性,造成摩擦片与内毂发生冲击碰撞,使摩擦片部分齿的齿面受到较大的冲击载荷,引起摩擦片轮齿的齿面与齿根发生疲劳损坏。
现有技术对具有齿侧间隙的传动件的冲击动载的研究很少考虑齿形制造误差引起的随机齿向误差的影响,尚未考虑随机齿向误差在重力作用下引起的偏心作用,使最终计算结果与实际值发生较大偏差。同时,考虑随机齿向误差的影响,研究浮动齿圈摩擦片组件的冲击载荷也未发现相关研究。
发明内容
针对现有技术中未考虑随机齿向误差在重力作用下引起的偏心作用,使获取的摩擦片冲击力与实际值发生较大偏差的问题,本发明提出一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,充分考虑了摩擦片组件轮齿的制造误差引起的随机齿向误差对浮动齿圈摩擦片组件冲击载荷的影响,保证较高计算效率的同时使获取的摩擦片冲击力与实际值误差较小,提高精准性。
为了实现上述目的,本发明提供以下技术方案:
一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,具体包括以下步骤:
S1:获取待计算的浮动支撑摩擦片组件参数,所述摩擦片组件参数包括基本参数和随机齿形制造误差量,所述摩擦片组件包括摩擦片和内毂;
S2、基于摩擦片组件的基本参数,利用势能法计算摩擦片与内毂单齿接触的综合刚度;
S3、基于摩擦片组件的基本参数和随机齿形制造误差量,计算摩擦片在重力作用下引起的偏心距,以及每时刻摩擦片的时变偏心距和时变偏心相位;
S4、根据所述摩擦片与内毂的随机齿侧间隙和相对位置关系,判断摩擦片与内毂是否发生碰撞,确定内毂与摩擦片碰撞齿数和随机碰撞齿位置;
S5、基于能量平衡理论和内毂与摩擦片的相对运动关系,计算摩擦片与内毂轮齿间碰撞的冲击反弹系数和迟滞阻尼系数;
S6、根据S1中确定的随机齿形制造误差量、S2中确定的摩擦片单齿综合刚度、S3中确定的摩擦片与内毂的时变偏心距和时变偏心相位、S4中确定的内毂与摩擦片碰撞齿数和S5中确定的摩擦片与内毂轮齿间的冲击反弹系数和迟滞阻尼系数,计算摩擦片与内毂每个接触齿的碰撞力;
S7、根据S6中确定的内毂与摩擦片每个轮齿的碰撞力,构建内毂与摩擦片平动-扭转多自由度冲击动载动力学模型,输出每时刻内毂与摩擦片的冲击动载。
优选地,所述S1中,摩擦片组件的基本参数包括:摩擦片与内毂的齿数、模数、齿宽、分度圆压力角、弹性模量、泊松比、摩擦片与内毂质量、转动惯量、齿侧间隙、内毂的转速、转速波动频率和幅值。
优选地,所述S2包括以下步骤:
S201、根据摩擦片组件的参数,利用势能法分别计算内毂或摩擦片的轮齿单齿对的弯曲刚度kbj、剪切刚度ksj和轴向压缩刚度kaj,对应的关系式如下:
公式(1)中,j=1时代表内毂;j=2时代表摩擦片;Uaj、Ubj和Usj分别为内毂轮齿或摩擦片轮齿的轴向压缩、弯曲和剪切变形存储的弹性势能;Fj表示内毂轮齿或摩擦片轮齿的啮合点接触载荷;Faj为内毂轮齿或摩擦片轮齿的啮合力在水平方向的分量;Ej表示内毂轮齿或摩擦片轮齿的弹性模量;Axj表示距离内毂轮齿或摩擦片轮齿的固定端x处截面的面积;dj表示内毂轮齿或摩擦片轮齿啮合位置距齿根圆距离;Mj为距离内毂轮齿或摩擦片轮齿的固定端x处的弯矩;Ixj表示距离内毂轮齿或摩擦片轮齿的固定端x处截面的惯性矩;Fbj为内毂轮齿或摩擦片轮齿的啮合力在竖直方向的分量;Gj表示内毂轮齿或摩擦片轮齿的剪切模量;
S202、根据摩擦片组件的参数,利用O’Donnell变形理论计算轮齿单齿对的基体刚度kfj,表达式如下:
公式(2)中,j=1时代表内毂;j=2时代表摩擦片;θfj为内毂或摩擦片的轮齿基体的变形;Mfj为距离内毂轮齿或摩擦片轮齿的固定端x处的单位弯矩;Vj表示距离内毂轮齿或摩擦片轮齿的固定端x处的单位剪切力;Ej表示内毂轮齿或摩擦片轮齿的弹性模量;υj表示内毂轮齿或摩擦片轮齿的泊松比;hj’表示距离内毂轮齿或摩擦片轮齿固定端x处截面;Fj表示内毂轮齿或摩擦片轮齿的啮合点接触载荷;
S203、根据S201中内毂和摩擦片轮齿的弯曲刚度kbj、剪切刚度ksj、轴向压缩刚度kaj和S202中基体刚度kfj,获取摩擦片组件轮齿单齿对接触的综合刚度K,计算如下:
公式(3)中,K表示摩擦片组件轮齿单齿对接触的综合刚度。
优选地,所述S3包括以下步骤:
S301、内毂与摩擦片之间存在随机齿侧间隙且服从高斯分布,形式如下:
ci=normrnd(μ,σ2) (4)
公式(4)中,ci表示不偏心状态第i个齿的齿侧间隙;μ表示齿侧间隙理论值;σ2表示齿侧间隙的离散程度,齿侧间隙的范围在(μ-3σ,μ+3σ)之间分布;
S302、重力作用下摩擦片工作齿面先接触产生的偏心距计算如下:
公式(5)中,aR表示重力作用下摩擦片工作齿面接触引起的偏心距;abs(…)表示绝对值;min(…)表示最小值;ci表示不偏心状态第i个齿的齿侧间隙;α0表示压力角;Z为摩擦片和内毂的齿数,i为第i个齿的齿序号;
重力作用下摩擦片非工作齿面先接触产生的偏心距计算如下:
公式(6)中,aL表示重力作用下摩擦片非工作齿面接触引起的偏心距;abs(…)表示绝对值;min(…)表示最小值;ci表示不偏心状态第i个齿的齿侧间隙;α0表示压力角;Z为摩擦片和内毂的齿数,i为第i个齿的齿序号;
则重力作用下摩擦片产生的初始自偏心距表示为:
ab=(aL+aR)/2;
公式(7)中,a表示任意时刻摩擦片与内毂的偏心距;表示偏心相位;x为每时刻摩擦片的质心相对于内毂质心在水平方向(x方向)的位移;y为每时刻摩擦片的质心相对于内毂质心在竖直方向(y方向)的位移;sign(x)表示为符号函数。
优选地,所述S4包括以下步骤:
S401、利用S1所述的内毂与摩擦片的基本参数,计算任一时刻第i个齿的工作齿面与非工作齿面的齿侧间隙,计算形式如下:
公式(8)中,cRi表示自偏心状态第i个齿的工作齿面的齿侧间隙;cLi表示自偏心状态第i个齿的非工作齿面的齿侧间隙;ci表示不偏心状态第i个齿的齿侧间隙;a表示任意时刻摩擦片与内毂的偏心距;α0表示压力角;Z为摩擦片和内毂的齿数;i为第i个齿的齿序号;
S402、考虑摩擦副组件在扭矩作用下的相对位移,计算内毂与摩擦片的工作齿面与非工作齿面的变形量,计算形式如下:
δRi=R(θ1-θ2)cosα0-cRi,δLi=R(θ2-θ1)cosα0-cLi (9)
公式(9)中,δRi表示内毂与摩擦片第i个齿的工作齿面的变形量;δLi表示内毂与摩擦片第i个齿的非工作齿面的变形量;R为分度圆半径;θ1为内毂角位移;θ2为摩擦片角位移;α0为压力角;cRi表示自偏心状态第i个齿的工作齿面的齿侧间隙;cLi表示自偏心状态第i个齿的非工作齿面的齿侧间隙;
S403、根据S402的工作齿面变形量δRi和非工作齿面变形量δLi的计算结果进行碰撞判断:
若δRi>0则摩擦片与内毂的第i个工作齿面发生碰撞(正碰),如图4所示,为正碰接触状态示意图,记CRi=1;若δLi>0则摩擦片与内毂的第i个非工作齿面发生碰撞(反碰),如图5所示,为反碰接触状态示意图,记CLi=1;若δRi≤0或δLi≤0,则摩擦片与内毂的第i个齿不发生碰撞,记CRi=0或CLi=0;
S404、根据S403的计算结果,分别获得内毂与摩擦片正碰和反碰的碰撞齿数,表示为:
公式(10)中,Z为摩擦片和内毂的齿数;i为第i个齿的齿序号;
则内毂与摩擦片轮齿总碰撞齿数为:
N=NR+NL。
优选地,所述S5包括以下步骤:
S501、基于赫兹接触理论,非线性冲击碰撞力F表示为:
S502、轮齿碰撞开始接触到碰撞结束分离时的能量损耗ΔT表示为:
冲击反弹系数e定义为两物体碰撞前的相对接近速度与碰撞后的相对分离速度的比值,表示为:
公式(14)中,δ0为轮齿碰撞的最大变形量;K表示摩擦片组件轮齿单齿对接触的综合刚度;
S504、根据能量守恒定理,结合式(12)和式(14)得到迟滞阻尼系数μ,表示为:
优选地,所述S6包括以下步骤:
S601、轮齿的冲击力合力依靠内毂轮齿与摩擦片轮齿相互碰撞来传递,结合式(3)、式(9)和式(15),任意时刻内毂与摩擦片工作齿面与非工作齿面的碰撞力计算为
公式(16)中,Fi表示摩擦片与内毂第i个齿的碰撞力;p=μ/K;μ为迟滞阻尼系数;K表示摩擦片组件轮齿单齿对接触的综合刚度;为内毂角速度;为摩擦片角速度;R表示分度圆半径;δL表示内毂与摩擦片轮齿非工作齿面的变形量;δR表示内毂与摩擦片轮齿工作齿面的变形量;
S602、第i个齿正碰和反碰的法向冲击碰撞力在x方向和y方向上碰撞力分量分别表示为:
公式(17)中,Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;Fi表示摩擦片与内毂第i个齿的碰撞力;α0为压力角;i为第i个齿;Z为齿的总数;ω为浮动齿圈角速度;θ0为分度圆齿厚的半齿角,δRi表示内毂与摩擦片第i个齿的工作齿面的变形量;δLi表示内毂与摩擦片第i个齿的非工作齿面的变形量;t表示为时间。
优选地,所述S7包括以下步骤:
S701、由于内毂扭转振动造成摩擦片的转速波动,内毂输入角速度ω及输入角加速度α表示为
ω=ω0+Asin(2πft),α=2πfAcos(2πft) (18)
公式(18)中,ω表示内毂转速波动,ω0表示恒定转速;A为转速波动幅值;f为转速波动频率;t表示t时刻;α表示内毂输入角加速度;
S702、建立内毂平动-扭转动力学方程,其形式如下:
公式(19)中,m1表示内毂的质量;表示内毂在x方向的加速度;x1表示内毂在x方向的位移;表示内毂在x方向的速度;表示内毂在y方向的加速度;y1表示内毂在y方向的位移;表示内毂在y方向的速度;k1x表示内毂在x方向的支撑刚度;c1x表示内毂在x方向的支撑阻尼;Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;J1表示内毂的转动惯量;为内毂角加速度;α表示内毂输入角加速度;Fi表示摩擦片与内毂第i个齿的碰撞力;i表示第i个齿;Z为轮齿的总数;R为分度圆半径;α0为压力角;
S703、建立摩擦片平动-扭转动力学方程,其形式如下
公式(20)中,m2表示摩擦片的质量;表示摩擦片在x方向的加速度;x2表示摩擦片在x方向的位移;表示摩擦片在x方向的加速度;表示摩擦片在y方向的加速度;y2表示摩擦片在y方向的位移;表示摩擦片在y方向的速度;k2x表示摩擦片在x方向的支撑刚度;c2x表示摩擦片在x方向的支撑阻尼;Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;J2表示摩擦片的转动惯量;为摩擦片角加速度;α0为压力角;Fi表示摩擦片与内毂第i个齿的碰撞力;i表示第i个齿;Z为轮齿的总数;R为分度圆半径;
S704、对S702中构建的内毂平动-扭转动力学方程和S703中构建的摩擦片平动-扭转动力学方程时间t内迭代计算,输出每时刻内毂与摩擦片的所有齿的总冲击载荷和单个轮齿的冲击载荷。
综上所述,由于采用了上述技术方案,与现有技术相比,本发明至少具有以下有益效果:
本发明提供了一种利用能量法求解内毂与摩擦片的单齿接触的综合刚度的方法;
本发明还提供了一种摩擦片与内毂在制造过程中引起随机齿向误差的情况下摩擦片在重力作用下产生的初始偏心距计算方法,以及每时刻摩擦片的时变偏心距和时变偏心相位的迭代计算方法。
本发明还提供了一种摩擦片与内毂的轮齿间存在随机齿向误差和摩擦片产生时变偏心距和时变偏心相位的情况下,内毂与摩擦片齿面任意时刻发生正碰、反碰的确定方法以及碰撞齿数和随机碰撞齿位置的计算方法;
本发明还提供了一种摩擦片与内毂与摩擦片轮齿碰撞前后的冲击反弹系数的确定方法;
本发明还提供了一种内毂与摩擦片的任意时刻每个轮齿发生正碰、反碰的齿部碰撞力的计算方法;
本发明还提供了一种内毂与摩擦片的平动-扭转多自由度冲击动载动力学模型迭代计算方法,确定每时刻内毂与摩擦片的冲击载荷,保证较高计算效率的同时使获取的摩擦片冲击力与实际值误差较小,提高获取的摩擦片冲击力精准性。
附图说明:
图1为根据本发明示例性实施例的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法示意图。
图2为根据本发明示例性实施例的不偏心随机齿侧间隙初始状态示意图。
图3为根据本发明示例性实施例的重力作用下初始偏心距状态示意图。
图4为根据本发明示例性实施例的正碰接触状态示意图。
图5为根据本发明示例性实施例的反碰接触状态示意图。
具体实施方式
下面结合实施例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
如图1所示,本发明提供一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,包括以下步骤:
S1、获取待计算的浮动支撑摩擦片组件参数,所述摩擦片组件参数包括基本参数和随机齿形制造误差量,所述摩擦片组件包括摩擦片和内毂。
本实施例中,浮动支撑摩擦片组件的基本参数包括:摩擦片与内毂的齿数、模数、齿宽、分度圆压力角、弹性模量、泊松比、摩擦片与内毂质量、转动惯量、齿侧间隙、内毂的转速、转速波动频率和幅值等。
S2、基于摩擦片组件的基本参数,利用势能法计算摩擦片与内毂单齿接触的综合刚度。
S201、根据摩擦片组件的参数,利用势能法分别计算内毂或摩擦片的轮齿单齿对的弯曲刚度kbj、剪切刚度ksj和轴向压缩刚度kaj,对应的关系式如下:
公式(1)中,j=1时代表内毂;j=2时代表摩擦片;Uaj、Ubj和Usj分别为内毂轮齿或摩擦片轮齿的轴向压缩、弯曲和剪切变形存储的弹性势能;Fj表示内毂轮齿或摩擦片轮齿的啮合点接触载荷;Faj为内毂轮齿或摩擦片轮齿的啮合力在水平方向的分量;Ej表示内毂轮齿或摩擦片轮齿的弹性模量;Axj表示距离内毂轮齿或摩擦片轮齿的固定端x处截面的面积;dj表示内毂轮齿或摩擦片轮齿啮合位置距齿根圆距离;Mj为距离内毂轮齿或摩擦片轮齿的固定端x处的弯矩;Ixj表示距离内毂轮齿或摩擦片轮齿的固定端x处截面的惯性矩;Fbj为内毂轮齿或摩擦片轮齿的啮合力在竖直方向的分量;Gj表示内毂轮齿或摩擦片轮齿的剪切模量;
S202、根据摩擦片组件的参数,利用O’Donnell变形理论计算轮齿单齿对的基体刚度kfj,表达式如下:
公式(2)中,j=1时代表内毂;j=2时代表摩擦片;θfj为内毂或摩擦片的轮齿基体的变形;Mfj为距离内毂轮齿或摩擦片轮齿的固定端x处的单位弯矩;Vj表示距离内毂轮齿或摩擦片轮齿的固定端x处的单位剪切力;Ej表示内毂轮齿或摩擦片轮齿的弹性模量;υj表示内毂轮齿或摩擦片轮齿的泊松比;hj’表示距离内毂轮齿或摩擦片轮齿固定端x处截面;Fj表示内毂轮齿或摩擦片轮齿的啮合点接触载荷;
S203、根据S201中内毂和摩擦片轮齿的弯曲刚度kbj、剪切刚度ksj、轴向压缩刚度kaj和S202中基体刚度kfj,获取摩擦片组件轮齿单齿对接触的综合刚度K,计算如下:
S3、基于摩擦片组件的基本参数和随机齿形制造误差量,计算摩擦片在重力作用下引起的偏心距,以及每时刻摩擦片的时变偏心距和时变偏心相位。
S301、内毂与摩擦片在加工过程中存在误差,会产生随机齿侧间隙,如图2所示,不偏心状态下摩擦片与内毂的轮齿之间存在多个单侧齿侧间隙c1、c2、ci、cn,假设单侧齿侧间隙服从高斯分布,计算形式如下
ci=normrnd(μ,σ2) (4)
公式(4)中,ci表示不偏心状态第i个齿的齿侧间隙;μ表示齿侧间隙理论值;σ2表示齿侧间隙的离散程度;齿侧间隙的范围在(μ-3σ,μ+3σ)之间分布。
S302、如图3所示,因此重力作用下摩擦片会产生初始自偏心距ab;
本实施例中,重力作用下摩擦片右侧齿(工作齿面)先接触产生的偏心距计算如下:
公式(5)中,aR表示重力作用下摩擦片右侧齿(工作齿面)接触引起的偏心距;abs(…)表示绝对值;min(…)表示最小值;ci表示不偏心状态第i个齿的齿侧间隙;α0表示压力角;Z为摩擦片和内毂的齿数,i为第i个齿的齿序号;
重力作用下摩擦片左侧齿(非工作齿面)先接触产生的偏心距计算如下:
公式(6)中,aL表示重力作用下摩擦片左侧齿(工作齿面)接触引起的偏心距;abs(…)表示绝对值;min(…)表示最小值;ci表示不偏心状态第i个齿的齿侧间隙;α0表示压力角;Z为摩擦片和内毂的齿数,i为第i个齿的齿序号;
重力作用下摩擦片产生的初始自偏心距ab可表示为:
ab=(aL+aR)/2。
公式(7)中,a表示任意时刻摩擦片与内毂的偏心距;表示偏心相位;x为每时刻摩擦片的质心相对于内毂质心在水平方向(x方向)的位移;y为每时刻摩擦片的质心相对于内毂质心在竖直方向(y方向)的位移;sign(x)表示为符号函数。
S4、根据所述摩擦片与内毂的随机齿侧间隙和相对位置关系,判断摩擦片与内毂是否发生碰撞,确定内毂与摩擦片碰撞齿数和随机碰撞齿位置。
S401、利用S1所述的内毂与摩擦片的基本参数,计算任一时刻第i个齿的工作齿面与非工作齿面的齿侧间隙,计算形式如下:
公式(8)中,cRi表示自偏心状态第i个齿的工作齿面的齿侧间隙;cLi表示自偏心状态第i个齿的非工作齿面的齿侧间隙;ci表示不偏心状态第i个齿的齿侧间隙;a表示任意时刻摩擦片与内毂的偏心距;α0表示压力角;Z为摩擦片和内毂的齿数;i为第i个齿的齿序号;表示偏心相位。
S402、考虑摩擦副组件在扭矩作用下的相对位移,计算内毂与摩擦片的工作齿面与非工作齿面的变形量,计算形式如下:
δRi=R(θ1-θ2)cosα0-cRi,δLi=R(θ2-θ1)cosα0-cLi (9)
公式(9)中,δRi表示内毂与摩擦片第i个齿的工作齿面的变形量;δLi表示内毂与摩擦片第i个齿的非工作齿面的变形量;R为分度圆半径;θ1为内毂角位移;θ2为摩擦片角位移;α0为压力角;cRi表示自偏心状态第i个齿的工作齿面的齿侧间隙;cLi表示自偏心状态第i个齿的非工作齿面的齿侧间隙。
S403、根据S402的工作齿面变形量δRi和非工作齿面变形量δLi的计算结果进行碰撞判断:
若δRi>0则摩擦片与内毂的第i个工作齿面发生碰撞(正碰),如图4所示,为正碰接触状态示意图,记CRi=1;若δLi>0则摩擦片与内毂的第i个非工作齿面发生碰撞(反碰),如图5所示,为反碰接触状态示意图,记CLi=1;若δRi≤0或δLi≤0,则摩擦片与内毂的第i个齿不发生碰撞,记CRi=0或CLi=0,通过CRi和CLi可判断碰撞齿位置。
S404、根据S403的计算结果,可分别获得内毂与摩擦片正碰和反碰的碰撞齿数,可表示为:
公式(10)中,Z为摩擦片和内毂的齿数;i为第i个齿的齿序号;
则内毂与摩擦片轮齿总碰撞齿数为:
N=NR+NL。
S5、基于能量平衡理论和内毂与摩擦片的相对运动关系,计算摩擦片与内毂轮齿间碰撞的冲击反弹系数和迟滞阻尼系数。
S501、基于赫兹接触理论,非线性冲击碰撞力可表示为:
S502、轮齿碰撞开始接触到碰撞结束分离时的能量损耗可表示为:
冲击反弹系数e定义为两物体碰撞前的相对接近速度与碰撞后的相对分离速度的比值,表示为:
公式(14)中,δ0为轮齿碰撞的最大变形量。
S504、根据能量守恒定理,结合式(12)和式(14)可得到迟滞阻尼系数μ,可表示为:
S6、根据S1中确定的随机齿形制造误差量、S2中确定的摩擦片单齿综合刚度、S3中确定的摩擦片与内毂的时变偏心距和时变偏心相位、S4中确定的内毂与摩擦片碰撞齿数和S5中确定的摩擦片与内毂轮齿间的冲击反弹系数和迟滞阻尼系数,基于Hertz接触理论,计算摩擦片与内毂每个接触齿的碰撞力。
S601、轮齿的冲击力合力依靠内毂轮齿与摩擦片轮齿相互碰撞来传递,采用Hertz接触力学模型可以描述接触面之间的非线性冲击碰撞作用。基于Hertz接触理论,结合式(3)、式(9)和式(15),任意时刻内毂与摩擦片工作齿面与非工作齿面的碰撞力可计算为
公式(16)中,Fi表示摩擦片与内毂第i个齿的碰撞力;p=μ/K;μ为迟滞阻尼系数;K表示摩擦片组件轮齿单齿对接触的综合刚度;为内毂角速度;为摩擦片角速度;R表示分度圆半径;δL表示内毂与摩擦片轮齿非工作齿面的变形量;δR表示内毂与摩擦片轮齿工作齿面的变形量;
S602、第i个齿正碰和反碰的法向冲击碰撞力在x方向和y方向上碰撞力分量分别可表示为:
公式(17)中,Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;Fi表示摩擦片与内毂第i个齿的碰撞力;α0为压力角;i为第i个齿;Z为齿的总数;ω为浮动齿圈角速度;θ0为分度圆齿厚的半齿角,δRi表示内毂与摩擦片第i个齿的工作齿面的变形量;δLi表示内毂与摩擦片第i个齿的非工作齿面的变形量;t表示为时间。
S7、根据S6中确定的内毂与摩擦片每个轮齿的碰撞力,基于牛顿力学理论,构建内毂与摩擦片平动-扭转多自由度冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载。
S701、由于内毂扭转振动造成摩擦片的转速波动,内毂输入角速度ω及输入角加速度α可表示为
ω=ω0+Asin(2πft),α=2πfAcos(2πft) (18)
公式(18)中,ω表示内毂转速波动,ω0表示恒定转速;A为转速波动幅值;f为转速波动频率;t表示t时刻;α表示内毂输入角加速度。
S702、建立内毂平动-扭转动力学方程,其形式如下:
公式(19)中,m1表示内毂的质量;表示内毂在x方向(水平方向)的加速度;x1表示内毂在x方向(水平方向)的位移;表示内毂在x方向(水平方向)的速度;表示内毂在y方向(竖直方向)的加速度;y1表示内毂在y方向(竖直方向)的位移;表示内毂在y方向(竖直方向)的速度;k1x表示内毂在x方向的支撑刚度;c1x表示内毂在x方向的支撑阻尼;Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;J1表示内毂的转动惯量;为内毂角加速度;α表示内毂输入角加速度;Fi表示摩擦片与内毂第i个齿的碰撞力;i表示第i个齿;Z为齿数(轮齿的总数);R为分度圆半径;α0为压力角;
S703、建立摩擦片平动-扭转动力学方程,其形式如下
公式(20)中,m2表示摩擦片的质量;表示摩擦片在x方向(水平方向)的加速度;x2表示摩擦片在x方向(水平方向)的位移;表示摩擦片在x方向(水平方向)的加速度;表示摩擦片在y方向(竖直方向)的加速度;y2表示摩擦片在y方向(竖直方向)的位移;表示摩擦片在y方向(竖直方向)的速度;k2x表示摩擦片在x方向的支撑刚度;c2x表示摩擦片在x方向的支撑阻尼;Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;J2表示摩擦片的转动惯量;为摩擦片角加速度;α0为压力角;Fi表示摩擦片与内毂第i个齿的碰撞力;i表示第i个齿;Z为齿数(轮齿的总数);R为分度圆半径。
S704、利用龙格库塔法对S702中构建的内毂平动-扭转动力学方程和S703中构建的摩擦片平动-扭转动力学方程时间t内迭代计算,若t≤t1,更新θ1、θ2、x、y以及更新时变偏心距和时变偏心相位,并将t+1,重复S2-S6;若t>t1,输出每时刻内毂与摩擦片的所有齿的总冲击载荷和单个轮齿的冲击载荷。
本领域的普通技术人员可以理解,上述各实施方式是实现本发明的具体实施例,而在实际应用中,可以在形式上和细节上对其作各种改变,而不偏离本发明的精神和范围。
Claims (8)
1.一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,具体包括以下步骤:
S1:获取待计算的浮动支撑摩擦片组件参数,所述摩擦片组件参数包括基本参数和随机齿形制造误差量,所述摩擦片组件包括摩擦片和内毂;
S2、基于摩擦片组件的基本参数,利用势能法计算摩擦片与内毂单齿接触的综合刚度;
S3、基于摩擦片组件的基本参数和随机齿形制造误差量,计算摩擦片在重力作用下引起的偏心距,以及每时刻摩擦片的时变偏心距和时变偏心相位;
S4、根据所述摩擦片与内毂的随机齿侧间隙和相对位置关系,判断摩擦片与内毂是否发生碰撞,确定内毂与摩擦片碰撞齿数和随机碰撞齿位置;
S5、基于能量平衡理论和内毂与摩擦片的相对运动关系,计算摩擦片与内毂轮齿间碰撞的冲击反弹系数和迟滞阻尼系数;
S6、根据S1中确定的随机齿形制造误差量、S2中确定的摩擦片单齿综合刚度、S3中确定的摩擦片与内毂的时变偏心距和时变偏心相位、S4中确定的内毂与摩擦片碰撞齿数和S5中确定的摩擦片与内毂轮齿间的冲击反弹系数和迟滞阻尼系数,计算摩擦片与内毂每个接触齿的碰撞力;
S7、根据S6中确定的内毂与摩擦片每个轮齿的碰撞力,构建内毂与摩擦片平动-扭转多自由度冲击动载动力学模型,输出每时刻内毂与摩擦片的冲击动载。
2.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S1中,摩擦片组件的基本参数包括:摩擦片与内毂的齿数、模数、齿宽、分度圆压力角、弹性模量、泊松比、摩擦片与内毂质量、转动惯量、齿侧间隙、内毂的转速、转速波动频率和幅值。
3.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S2包括以下步骤:
S201、根据摩擦片组件的参数,利用势能法分别计算内毂或摩擦片的轮齿单齿对的弯曲刚度kbj、剪切刚度ksj和轴向压缩刚度kaj,对应的关系式如下:
公式(1)中,j=1时代表内毂;j=2时代表摩擦片;Uaj、Ubj和Usj分别为内毂轮齿或摩擦片轮齿的轴向压缩、弯曲和剪切变形存储的弹性势能;Fj表示内毂轮齿或摩擦片轮齿的啮合点接触载荷;Faj为内毂轮齿或摩擦片轮齿的啮合力在水平方向的分量;Ej表示内毂轮齿或摩擦片轮齿的弹性模量;Axj表示距离内毂轮齿或摩擦片轮齿的固定端x处截面的面积;dj表示内毂轮齿或摩擦片轮齿啮合位置距齿根圆距离;Mj为距离内毂轮齿或摩擦片轮齿的固定端x处的弯矩;Ixj表示距离内毂轮齿或摩擦片轮齿的固定端x处截面的惯性矩;Fbj为内毂轮齿或摩擦片轮齿的啮合力在竖直方向的分量;Gj表示内毂轮齿或摩擦片轮齿的剪切模量;
S202、根据摩擦片组件的参数,利用O’Donnell变形理论计算轮齿单齿对的基体刚度kfj,表达式如下:
公式(2)中,j=1时代表内毂;j=2时代表摩擦片;θfj为内毂或摩擦片的轮齿基体的变形;Mfj为距离内毂轮齿或摩擦片轮齿的固定端x处的单位弯矩;Vj表示距离内毂轮齿或摩擦片轮齿的固定端x处的单位剪切力;Ej表示内毂轮齿或摩擦片轮齿的弹性模量;υj表示内毂轮齿或摩擦片轮齿的泊松比;hj’表示距离内毂轮齿或摩擦片轮齿固定端x处截面;Fj表示内毂轮齿或摩擦片轮齿的啮合点接触载荷;
S203、根据S201中内毂和摩擦片轮齿的弯曲刚度kbj、剪切刚度ksj、轴向压缩刚度kaj和S202中基体刚度kfj,获取摩擦片组件轮齿单齿对接触的综合刚度K,计算如下:
公式(3)中,K表示摩擦片组件轮齿单齿对接触的综合刚度。
4.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S3包括以下步骤:
S301、内毂与摩擦片之间存在随机齿侧间隙且服从高斯分布,形式如下:
ci=normrnd(μ,σ2) (4)
公式(4)中,ci表示不偏心状态第i个齿的齿侧间隙;μ表示齿侧间隙理论值;σ2表示齿侧间隙的离散程度,齿侧间隙的范围在(μ-3σ,μ+3σ)之间分布;
S302、重力作用下摩擦片工作齿面先接触产生的偏心距计算如下:
公式(5)中,aR表示重力作用下摩擦片工作齿面接触引起的偏心距;abs(…)表示绝对值;min(…)表示最小值;ci表示不偏心状态第i个齿的齿侧间隙;α0表示压力角;Z为摩擦片和内毂的齿数,i为第i个齿的齿序号;
重力作用下摩擦片非工作齿面先接触产生的偏心距计算如下:
公式(6)中,aL表示重力作用下摩擦片非工作齿面接触引起的偏心距;abs(…)表示绝对值;min(…)表示最小值;ci表示不偏心状态第i个齿的齿侧间隙;α0表示压力角;Z为摩擦片和内毂的齿数,i为第i个齿的齿序号;
则重力作用下摩擦片产生的初始自偏心距表示为:
ab=(aL+aR)/2;
5.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S4包括以下步骤:
S401、利用S1所述的内毂与摩擦片的基本参数,计算任一时刻第i个齿的工作齿面与非工作齿面的齿侧间隙,计算形式如下:
公式(8)中,cRi表示自偏心状态第i个齿的工作齿面的齿侧间隙;cLi表示自偏心状态第i个齿的非工作齿面的齿侧间隙;ci表示不偏心状态第i个齿的齿侧间隙;a表示任意时刻摩擦片与内毂的偏心距;α0表示压力角;Z为摩擦片和内毂的齿数;i为第i个齿的齿序号;
S402、考虑摩擦副组件在扭矩作用下的相对位移,计算内毂与摩擦片的工作齿面与非工作齿面的变形量,计算形式如下:
δRi=R(θ1-θ2)cosα0-cRi,δLi=R(θ2-θ1)cosα0-cLi (9)
公式(9)中,δRi表示内毂与摩擦片第i个齿的工作齿面的变形量;δLi表示内毂与摩擦片第i个齿的非工作齿面的变形量;R为分度圆半径;θ1为内毂角位移;θ2为摩擦片角位移;α0为压力角;cRi表示自偏心状态第i个齿的工作齿面的齿侧间隙;cLi表示自偏心状态第i个齿的非工作齿面的齿侧间隙;
S403、根据S402的工作齿面变形量δRi和非工作齿面变形量δLi的计算结果进行碰撞判断:
若δRi>0则摩擦片与内毂的第i个工作齿面发生碰撞(正碰),如图4所示,为正碰接触状态示意图,记CRi=1;若δLi>0则摩擦片与内毂的第i个非工作齿面发生碰撞(反碰),如图5所示,为反碰接触状态示意图,记CLi=1;若δRi≤0或δLi≤0,则摩擦片与内毂的第i个齿不发生碰撞,记CRi=0或CLi=0;
S404、根据S403的计算结果,分别获得内毂与摩擦片正碰和反碰的碰撞齿数,表示为:
公式(10)中,Z为摩擦片和内毂的齿数;i为第i个齿的齿序号;
则内毂与摩擦片轮齿总碰撞齿数为:
N=NR+NL。
6.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S5包括以下步骤:
S501、基于赫兹接触理论,非线性冲击碰撞力F表示为:
S502、轮齿碰撞开始接触到碰撞结束分离时的能量损耗ΔT表示为:
冲击反弹系数e定义为两物体碰撞前的相对接近速度与碰撞后的相对分离速度的比值,表示为:
公式(14)中,δ0为轮齿碰撞的最大变形量;K表示摩擦片组件轮齿单齿对接触的综合刚度;
S504、根据能量守恒定理,结合式(12)和式(14)得到迟滞阻尼系数μ,表示为:
7.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S6包括以下步骤:
S601、轮齿的冲击力合力依靠内毂轮齿与摩擦片轮齿相互碰撞来传递,结合式(3)、式(9)和式(15),任意时刻内毂与摩擦片工作齿面与非工作齿面的碰撞力计算为
公式(16)中,Fi表示摩擦片与内毂第i个齿的碰撞力;p=μ/K;μ为迟滞阻尼系数;K表示摩擦片组件轮齿单齿对接触的综合刚度;为内毂角速度;为摩擦片角速度;R表示分度圆半径;δL表示内毂与摩擦片轮齿非工作齿面的变形量;δR表示内毂与摩擦片轮齿工作齿面的变形量;
S602、第i个齿正碰和反碰的法向冲击碰撞力在x方向和y方向上碰撞力分量分别表示为:
8.如权利要求1所述的一种考虑随机齿向误差的离合器摩擦片冲击力计算方法,其特征在于,所述S7包括以下步骤:
S701、由于内毂扭转振动造成摩擦片的转速波动,内毂输入角速度ω及输入角加速度α表示为
ω=ω0+A sin(2πft),α=2πfA cos(2πft) (18)
公式(18)中,ω表示内毂转速波动,ω0表示恒定转速;A为转速波动幅值;f为转速波动频率;t表示t时刻;α表示内毂输入角加速度;
S702、建立内毂平动-扭转动力学方程,其形式如下:
公式(19)中,m1表示内毂的质量;表示内毂在x方向的加速度;x1表示内毂在x方向的位移;表示内毂在x方向的速度;表示内毂在y方向的加速度;y1表示内毂在y方向的位移;表示内毂在y方向的速度;k1x表示内毂在x方向的支撑刚度;c1x表示内毂在x方向的支撑阻尼;Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;J1表示内毂的转动惯量;为内毂角加速度;α表示内毂输入角加速度;Fi表示摩擦片与内毂第i个齿的碰撞力;i表示第i个齿;Z为轮齿的总数;R为分度圆半径;α0为压力角;
S703、建立摩擦片平动-扭转动力学方程,其形式如下
公式(20)中,m2表示摩擦片的质量;表示摩擦片在x方向的加速度;x2表示摩擦片在x方向的位移;表示摩擦片在x方向的加速度;表示摩擦片在y方向的加速度;y2表示摩擦片在y方向的位移;表示摩擦片在y方向的速度;k2x表示摩擦片在x方向的支撑刚度;c2x表示摩擦片在x方向的支撑阻尼;Fxi表示第i个齿的法向冲击碰撞力在x方向上碰撞力分量;Fyi表示第i个齿的法向冲击碰撞力在y方向上碰撞力分量;J2表示摩擦片的转动惯量;为摩擦片角加速度;α0为压力角;Fi表示摩擦片与内毂第i个齿的碰撞力;i表示第i个齿;Z为轮齿的总数;R为分度圆半径;
S704、对S702中构建的内毂平动-扭转动力学方程和S703中构建的摩擦片平动-扭转动力学方程时间t内迭代计算,输出每时刻内毂与摩擦片的所有齿的总冲击载荷和单个轮齿的冲击载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210857752.4A CN115098975B (zh) | 2022-07-20 | 2022-07-20 | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210857752.4A CN115098975B (zh) | 2022-07-20 | 2022-07-20 | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115098975A true CN115098975A (zh) | 2022-09-23 |
CN115098975B CN115098975B (zh) | 2024-07-05 |
Family
ID=83298456
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210857752.4A Active CN115098975B (zh) | 2022-07-20 | 2022-07-20 | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115098975B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006045241A1 (fr) * | 2004-10-25 | 2006-05-04 | Wenming Han | Transmission mecanique a variation continue |
CN103901308A (zh) * | 2014-04-22 | 2014-07-02 | 扬州电力设备修造厂 | 一种电液制动的电动执行机构智能型测试台及其试验方法 |
CN105987094A (zh) * | 2015-02-13 | 2016-10-05 | 中国北方车辆研究所 | 一种改进的离合器 |
CN111563317A (zh) * | 2020-04-14 | 2020-08-21 | 中国北方车辆研究所 | 一种浮动支撑摩擦片喷丸强化评估方法 |
CN115098977A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种浮动支撑摩擦片组件冲击动载计算方法 |
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115130247A (zh) * | 2022-07-20 | 2022-09-30 | 重庆大学 | 一种具有轮齿裂纹的浮动支撑摩擦片冲击动载计算方法 |
CN115203957A (zh) * | 2022-07-20 | 2022-10-18 | 重庆大学 | 一种考虑轮齿碰撞激励下的离合器摩擦冲击载荷计算方法 |
-
2022
- 2022-07-20 CN CN202210857752.4A patent/CN115098975B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006045241A1 (fr) * | 2004-10-25 | 2006-05-04 | Wenming Han | Transmission mecanique a variation continue |
CN103901308A (zh) * | 2014-04-22 | 2014-07-02 | 扬州电力设备修造厂 | 一种电液制动的电动执行机构智能型测试台及其试验方法 |
CN105987094A (zh) * | 2015-02-13 | 2016-10-05 | 中国北方车辆研究所 | 一种改进的离合器 |
CN111563317A (zh) * | 2020-04-14 | 2020-08-21 | 中国北方车辆研究所 | 一种浮动支撑摩擦片喷丸强化评估方法 |
CN115098977A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种浮动支撑摩擦片组件冲击动载计算方法 |
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115130247A (zh) * | 2022-07-20 | 2022-09-30 | 重庆大学 | 一种具有轮齿裂纹的浮动支撑摩擦片冲击动载计算方法 |
CN115203957A (zh) * | 2022-07-20 | 2022-10-18 | 重庆大学 | 一种考虑轮齿碰撞激励下的离合器摩擦冲击载荷计算方法 |
Non-Patent Citations (5)
Title |
---|
MARCIN BAK: "Torque capacity of multidisc wet clutch with reference to friction occurrence on its spline connections", SCIENTIFIC REPORTS, 29 October 2021 (2021-10-29), pages 1 - 17 * |
ZHENG CAO 等: "Effects of the gear eccentricities on the dynamic performance of a planetary gears", NONLINEAR DYNAMICS, vol. 91, 6 November 2017 (2017-11-06), pages 1 - 15, XP036388752, DOI: 10.1007/s11071-017-3738-0 * |
李杰 等: "高能摩擦副齿部啮合非线性冲击损伤新算法及疲劳寿命预测", 中国公路学报, vol. 33, no. 08, 8 April 2020 (2020-04-08), pages 110 - 118 * |
王利明 等: "二十辊轧机轧辊磨床砂轮动不平衡对磨削颤振的影响", 工程科学学报, vol. 37, no. 1, 15 May 2015 (2015-05-15), pages 78 - 83 * |
赫雷 等: "导轨制动器卸载机构的设计与仿真", 机械制造与自动化, vol. 45, no. 06, 20 December 2016 (2016-12-20), pages 110 - 113 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115098976B (zh) * | 2022-07-20 | 2024-05-28 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115098975B (zh) | 2024-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115098976B (zh) | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 | |
CN115098977B (zh) | 一种浮动支撑摩擦片组件冲击动载计算方法 | |
Ma et al. | Evaluation of optimum profile modification curves of profile shifted spur gears based on vibration responses | |
CN115130247A (zh) | 一种具有轮齿裂纹的浮动支撑摩擦片冲击动载计算方法 | |
Xu et al. | Dynamic modeling and contact analysis of a cycloid-pin gear mechanism with a turning arm cylindrical roller bearing | |
CN115098975A (zh) | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 | |
CN115203957A (zh) | 一种考虑轮齿碰撞激励下的离合器摩擦冲击载荷计算方法 | |
CN110162909B (zh) | 一种渐开线直齿轮传动系统动态特性求解方法 | |
Chen et al. | Dynamics analysis of a crowned gear transmission system with impact damping: Based on experimental transmission error | |
Ooi et al. | Modal and stress analysis of gear train design in portal axle using finite element modeling and simulation | |
Yu et al. | Influence of the addendum modification on spur gear back-side mesh stiffness and dynamics | |
CN111488681B (zh) | 基于不确定性的斜齿轮副随机动力学建模方法 | |
Wu et al. | The effect of bending loads on the dynamic behaviors of a rolling guide | |
CN108416120B (zh) | 一种直齿圆柱齿轮双齿啮合区载荷分配率的确定方法 | |
CN107341313B (zh) | 基于adams的行星轮系非线性动力学建模方法 | |
Abbes et al. | Effect of transmission error on the dynamic behaviour of gearbox housing | |
Lin et al. | Nonlinear dynamic analysis of eccentric curve-face gear transmission system | |
CN110889215A (zh) | 一种基于adams的椭圆齿轮间隙的仿真分析方法 | |
CN115935694A (zh) | 一种含转动副间隙的微小齿轮动力学建模方法 | |
Chen et al. | Dynamic analysis of a geared rotor-bearing system with translational motion due to shaft deformation under residual shaft bow effect | |
Shing et al. | An improved model for the dynamics of spur gear systems with backlash consideration | |
CN103217924A (zh) | 面向实时控制应用的过约束重型并联机床动力学建模方法 | |
CN115045904A (zh) | 用于纳米定位平台的双轴柔性导向机构及其刚度建模方法 | |
Atanasiu et al. | Dynamic contact loads of spur gear pairs with addendum modifications | |
Lei | A multibody system model for meshing gears |
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 |