CN115098977A - 一种浮动支撑摩擦片组件冲击动载计算方法 - Google Patents
一种浮动支撑摩擦片组件冲击动载计算方法 Download PDFInfo
- Publication number
- CN115098977A CN115098977A CN202210857772.1A CN202210857772A CN115098977A CN 115098977 A CN115098977 A CN 115098977A CN 202210857772 A CN202210857772 A CN 202210857772A CN 115098977 A CN115098977 A CN 115098977A
- Authority
- CN
- China
- Prior art keywords
- friction plate
- inner hub
- gear teeth
- impact
- collision
- 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
- 238000007667 floating Methods 0.000 title claims abstract description 21
- 238000004364 calculation method Methods 0.000 title claims abstract description 11
- 238000013016 damping Methods 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims description 19
- 238000005452 bending Methods 0.000 claims description 18
- 230000001133 acceleration Effects 0.000 claims description 13
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000000926 separation method Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005381 potential energy Methods 0.000 claims description 6
- 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
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 239000011248 coating agent Substances 0.000 claims description 2
- 238000000576 coating method Methods 0.000 claims description 2
- 238000004134 energy conservation Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010998 test method Methods 0.000 description 2
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 230000001052 transient effect 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
- 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)
- 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)
- Gears, Cams (AREA)
- Mechanical Operated Clutches (AREA)
Abstract
本发明公开一种浮动支撑摩擦片组件冲击动载计算方法,具体包括以下步骤:S1、获取待计算的浮动支撑摩擦片组件的参数,所述摩擦片组件包括摩擦片和内毂;S2、基于摩擦片组件的参数,计算摩擦片组件轮齿单齿对综合刚度;S3、根据摩擦片与内毂的齿侧间隙和相对位置关系计算变形量,判断摩擦片与内毂是否发生碰撞;S4、计算摩擦片与内毂轮齿间的冲击反弹系数e和迟滞阻尼系数u;S5、根据摩擦片组件轮齿单齿对接触的综合刚度、内毂与摩擦片碰撞齿数、摩擦片与内毂轮齿间的冲击反弹系数和迟滞阻尼系数,计算摩擦片与内毂接触齿的碰撞力;S6、根据S5中确定的内毂与摩擦片的碰撞力构建冲击碰撞动力学模型,计算每时刻内毂与摩擦片的冲击动载。
Description
技术领域
本发明涉及摩擦片传动组件技术领域,特别涉及一种浮动支撑摩擦片组件冲击动载计算方法。
背景技术
摩擦片传动组件工作条件恶劣,存在大功率、高动载、高线速度、强瞬态等工作特点。摩擦片组件在浮动支撑条件下,由于变速机构的扭转振动等引起内毂转速表现为高频特征的非平稳性,造成摩擦片与内毂发生冲击碰撞,摩擦片齿部受到较大的冲击载荷,会引起摩擦片轮齿的齿面与齿根损伤等,影响摩擦片寿命和工作稳定性,直接导致离合器、制动器的工作性能下降,将严重影响车辆整体性能的提升,甚至导致恶性事故发生。
现有技术对浮动支撑摩擦片组件在工作过程的冲击载荷的计算主要集中在试验法或有限元法。但摩擦片传动组件运动形式复杂,利用有限元法难以准确设定摩擦片组件的边界条件,建模过程复杂且繁重,计算效率很低;试验法存在操作难度大、经济投入高、干扰因素多等问题。上述缺陷是本领域技术人员期望克服的。
发明内容
针对现有技术中浮动支撑摩擦片组件冲击载荷计算效率较低的问题,本发明提出一种浮动支撑摩擦片组件冲击动载计算方法,通过内毂与摩擦片的单齿对接触的综合刚度、冲击反弹系数、迟滞阻尼系数计算摩擦片与内毂接触齿的碰撞力F*,构建内毂与摩擦片扭转冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载,提高冲击载荷获取效率。
为了实现上述目的,本发明提供以下技术方案:
一种浮动支撑摩擦片组件冲击动载计算方法,具体包括以下步骤:
S1、获取待计算的浮动支撑摩擦片组件的参数,所述摩擦片组件包括摩擦片和内毂;
S2、基于摩擦片组件的参数,计算摩擦片组件轮齿单齿对综合刚度;
S3、根据摩擦片与内毂的齿侧间隙和相对位置关系计算内毂与摩擦片轮齿变形量,判断摩擦片与内毂是否发生碰撞;
S4、基于能量平衡理论和内毂与摩擦片的相对运动关系,计算摩擦片与内毂轮齿间的冲击反弹系数e和迟滞阻尼系数u;
S5、根据S2中摩擦片组件轮齿单齿对接触的综合刚度K、S3中确定的内毂与摩擦片轮齿的变形量、S4中确定的摩擦片与内毂轮齿间的冲击反弹系数e和迟滞阻尼系数u,计算摩擦片与内毂接触齿的碰撞力F*;
S6、根据S5中确定的内毂与摩擦片的碰撞力F*,构建内毂与摩擦片扭转冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载。
优选地,所述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处的单位剪切力;υj表示内毂轮齿或摩擦片轮齿的泊松比;hj’表示距离内毂轮齿或摩擦片轮齿固定端x处截面高度的1.5倍;
S203、根据S201中内毂和摩擦片轮齿的弯曲刚度kbj、剪切刚度ksj、轴向压缩刚度kaj和S202中基体刚度kfj,获取摩擦片组件轮齿单齿对接触的综合刚度K,计算如下:
优选地,所述S3包括以下步骤:
S301、利用S1所述的摩擦片组件的参数,考虑摩擦副组件在扭矩作用下的相对位移,计算内毂与摩擦片每对齿的工作齿面与非工作齿面的变形量,计算公式如下:
δR=R(θ1-θ2)·cosα0-c (4)
δL=R(θ2-θ1)·cosα0-c
公式(4)中,δR表示内毂与摩擦片在工作齿面的变形量;δL表示内毂与摩擦片在非工作齿面的变形量;θ1为内毂角位移;θ2为摩擦片角位移;c表示齿侧间隙;α0为压力角;R为分度圆半径;
S302、若δR>0,则说明内毂与摩擦片每对齿的工作齿面发生碰撞,即正碰,若δL>0,则说明内毂与摩擦片每对齿的非工作齿面发生碰撞,即反碰。
优选地,所述S4包括以下步骤:
S401、基于赫兹接触理论,非线性冲击碰撞力F表示为:
S402、轮齿碰撞开始接触到碰撞结束分离时的能量损耗ΔT表示为:
冲击反弹系数e定义为两物体碰撞前的相对接近速度与碰撞后的相对分离速度的比值,表示为:
S403、当轮齿碰撞以冲击阻尼力形式做功时,能量耗散ΔT表示为:
公式(8)中,δ0为轮齿碰撞的最大变形量,K表示摩擦片组件轮齿单齿对接触的综合刚度;
S404、根据能量守恒定理,结合式(6)和式(8)得到迟滞阻尼系数μ,表示为:
优选地,所述S5中,摩擦片与内毂单个轮齿的碰撞力计算为:
公式(10)中,F*表示摩擦片与内毂轮齿的碰撞力;p=μ/K;μ为迟滞阻尼系数;K表示摩擦片组件轮齿单齿对接触的综合刚度;为内毂角速度;为摩擦片角速度;R表示分度圆半径;δL表示内毂与摩擦片轮齿非工作齿面的变形量;δR表示内毂与摩擦片轮齿工作齿面的变形量。
优选地,所述S6包括以下步骤:
S601、由于内毂扭转振动造成摩擦片的转速波动,内毂输入转速ω及输入加速度α表示为
ω=ω0+Asin(2πft),α=2πfAcos(2πft) (11)
公式(11)中,ω表示内毂转速波动,ω0表示恒定转速;A为转速波动幅值;f为转速波动频率;t表示时间;
S602、建立内毂与摩擦片二自由度冲击碰撞动力学模型,其形式如下
公式(12)中,J1表示内毂的转动惯量;J2表示摩擦片的转动惯量;为内毂角加速度;为摩擦片角加速度;α表示内毂输入角加速度;Ftotal为内毂与摩擦片每时刻所有齿的总冲击载荷;Rb为基圆半径;z为齿数,即轮齿的总数;
S603、对S602中构建的内毂与摩擦片二自由度冲击碰撞动力学模型在时间t内迭代计算,输出每时刻内毂与摩擦片的冲击载荷。
综上所述,由于采用了上述技术方案,与现有技术相比,本发明至少具有以下有益效果:
本发明通过内毂与摩擦片的单齿对接触的综合刚度、冲击反弹系数、迟滞阻尼系数计算摩擦片与内毂接触齿的碰撞力F*,构建内毂与摩擦片扭转冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载,从而提高摩擦片组件冲击载荷获取效率。
本发明在摩擦片与内毂的轮齿间存在均匀齿侧间隙的情况下,确定内毂与摩擦片齿部发生正碰或反碰。
附图说明:
图1为根据本发明示例性实施例的一种浮动支撑摩擦片组件冲击动载计算方法流程示意图。
具体实施方式
下面结合实施例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
如图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处截面。
S203、根据S201中内毂或摩擦片轮齿的弯曲刚度kbj、剪切刚度ksj、轴向压缩刚度kaj和S202中基体刚度kfj,获取摩擦片组件轮齿单齿对接触的综合刚度K,计算如下:
S3、根据摩擦片与内毂的齿侧间隙和相对位置关系,计算内毂与摩擦片轮齿的变形量,并根据变形量判断摩擦片与内毂是否发生碰撞。
S301、利用S1所述的摩擦片组件的参数,考虑摩擦副组件在扭矩作用下的相对位移,计算内毂与摩擦片每对齿的工作齿面与非工作齿面的变形量,计算形式如下:
δR=R(θ1-θ2)·cosα0-c (4)
δL=R(θ2-θ1)·cosα0-c
公式(4)中,δR表示内毂与摩擦片在工作齿面的变形量;δL表示内毂与摩擦片在非工作齿面的变形量;θ1为内毂角位移;θ2为摩擦片角位移;c表示齿侧间隙;α0为压力角;R为分度圆半径;
S302、若δR>0,则说明内毂与摩擦片每对齿的工作齿面发生碰撞,即正碰,若δL>0,则说明内毂与摩擦片每对齿的非工作齿面发生碰撞,即反碰。
S4、基于能量平衡理论和内毂与摩擦片的相对运动关系,计算摩擦片与内毂轮齿间的冲击反弹系数和迟滞阻尼系数;
S401、基于赫兹接触理论,非线性冲击碰撞力可表示为:
S402、轮齿碰撞开始接触到碰撞结束分离时的能量损耗ΔT可表示为:
本实施例中,冲击反弹系数e定义为两物体碰撞前的相对接近速度与碰撞后的相对分离速度的比值,表示为:
S403、当轮齿碰撞以冲击阻尼力形式做功时,能量耗散ΔT可表示为:
公式(8)中,δ0为轮齿碰撞的最大变形量,K表示摩擦片组件轮齿单齿对接触的综合刚度。
S404、根据能量守恒定理,结合式(6)和式(8)可得到迟滞阻尼系数μ,可表示为:
S5、根据S2中摩擦片组件轮齿单齿对接触的综合刚度K、S3中确定的内毂与摩擦片轮齿的变形量、S4中确定的摩擦片与内毂轮齿间的冲击反弹系数e和迟滞阻尼系数u,基于Hertz接触理论,计算摩擦片与内毂轮齿的碰撞力。
本实施例中,轮齿的冲击力合力依靠内毂轮齿与摩擦片轮齿相互碰撞来传递,采用Hertz接触力学模型可以描述接触面之间的非线性冲击碰撞作用。基于Hertz接触理论,内毂与摩擦片工作齿面(正碰)与非工作齿面(反碰)的碰撞力可计算为:
公式(10)中,F*表示摩擦片与内毂轮齿的碰撞力;p=μ/K;μ为迟滞阻尼系数;K表示摩擦片组件轮齿单齿对接触的综合刚度;为内毂角速度;为摩擦片角速度;R表示分度圆半径;δL表示内毂与摩擦片轮齿非工作齿面的变形量;δR表示内毂与摩擦片轮齿工作齿面的变形量;
S6、根据S5中确定的内毂与摩擦片的碰撞力F*,基于牛顿力学理论,构建内毂与摩擦片扭转冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载。
S601、由于内毂扭转振动造成摩擦片的转速波动,内毂输入转速波动ω及输入角加速度α可表示为
ω=ω0+Asin(2πft),α=2πfAcos(2πft) (11)
公式(11)中,ω表示内毂转速波动,ω0表示恒定转速;A为转速波动幅值;f为转速波动频率;α表示内毂输入角加速度;t为时间。
S602、建立内毂与摩擦片二自由度冲击碰撞动力学模型,其形式如下
公式(12)中,J1表示内毂的转动惯量;J2表示摩擦片的转动惯量;为内毂角加速度;为摩擦片角加速度;α表示内毂输入角加速度;Ftotal为内毂与摩擦片每时刻所有齿的总冲击载荷;Rb为基圆半径;z为齿数,即轮齿的总数。
本实施例中,将S602构建的内毂与摩擦片二自由度冲击碰撞动力学模型在时间t(阈值设为t1)内迭代计算,若t≤t1,则将t+1,重复S2-S6;若t>t1,输出每时刻内毂与摩擦片的所有齿的总冲击载荷和单个轮齿的冲击载荷。
本领域的普通技术人员可以理解,上述各实施方式是实现本发明的具体实施例,而在实际应用中,可以在形式上和细节上对其作各种改变,而不偏离本发明的精神和范围。
Claims (7)
1.一种浮动支撑摩擦片组件冲击动载计算方法,其特征在于,具体包括以下步骤:
S1、获取待计算的浮动支撑摩擦片组件的参数,所述摩擦片组件包括摩擦片和内毂;
S2、基于摩擦片组件的参数,计算摩擦片组件轮齿单齿对综合刚度;
S3、根据摩擦片与内毂的齿侧间隙和相对位置关系计算内毂与摩擦片轮齿变形量,判断摩擦片与内毂是否发生碰撞;
S4、基于能量平衡理论和内毂与摩擦片的相对运动关系,计算摩擦片与内毂轮齿间的冲击反弹系数e和迟滞阻尼系数u;
S5、根据S2中摩擦片组件轮齿单齿对接触的综合刚度K、S3中确定的内毂与摩擦片轮齿的变形量、S4中确定的摩擦片与内毂轮齿间的冲击反弹系数e和迟滞阻尼系数u,计算摩擦片与内毂接触齿的碰撞力F*;
S6、根据S5中确定的内毂与摩擦片的碰撞力F*,构建内毂与摩擦片扭转冲击动载动力学模型,计算每时刻内毂与摩擦片的冲击动载。
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处的单位剪切力;υj表示内毂轮齿或摩擦片轮齿的泊松比;hj’表示距离内毂轮齿或摩擦片轮齿固定端x处截面高度的1.5倍;
S203、根据S201中内毂和摩擦片轮齿的弯曲刚度kbj、剪切刚度ksj、轴向压缩刚度kaj和S202中基体刚度kfj,获取摩擦片组件轮齿单齿对接触的综合刚度K,计算如下:
4.如权利要求1所述的一种浮动支撑摩擦片组件冲击动载计算方法,其特征在于,所述S3包括以下步骤:
S301、利用S1所述的摩擦片组件的参数,考虑摩擦副组件在扭矩作用下的相对位移,计算内毂与摩擦片每对齿的工作齿面与非工作齿面的变形量,计算公式如下:
δR=R(θ1-θ2)·cosα0-c (4)
δL=R(θ2-θ1)·cosα0-c
公式(4)中,δR表示内毂与摩擦片在工作齿面的变形量;δL表示内毂与摩擦片在非工作齿面的变形量;θ1为内毂角位移;θ2为摩擦片角位移;c表示齿侧间隙;α0为压力角;R为分度圆半径;
S302、若δR>0,则说明内毂与摩擦片每对齿的工作齿面发生碰撞,即正碰,若δL>0,则说明内毂与摩擦片每对齿的非工作齿面发生碰撞,即反碰。
5.如权利要求1所述的一种浮动支撑摩擦片组件冲击动载计算方法,其特征在于,所述S4包括以下步骤:
S401、基于赫兹接触理论,非线性冲击碰撞力F表示为:
S402、轮齿碰撞开始接触到碰撞结束分离时的能量损耗ΔT表示为:
冲击反弹系数e定义为两物体碰撞前的相对接近速度与碰撞后的相对分离速度的比值,表示为:
S403、当轮齿碰撞以冲击阻尼力形式做功时,能量耗散ΔT表示为:
公式(8)中,δ0为轮齿碰撞的最大变形量,K表示摩擦片组件轮齿单齿对接触的综合刚度;
S404、根据能量守恒定理,结合式(6)和式(8)得到迟滞阻尼系数μ,表示为:
7.如权利要求1所述的一种浮动支撑摩擦片组件冲击动载计算方法,其特征在于,所述S6包括以下步骤:
S601、由于内毂扭转振动造成摩擦片的转速波动,内毂输入转速ω及输入加速度α表示为
ω=ω0+Asin(2πft),α=2πfAcos(2πft) (11)
公式(11)中,ω表示内毂转速波动,ω0表示恒定转速;A为转速波动幅值;f为转速波动频率;t表示时间;
S602、建立内毂与摩擦片二自由度冲击碰撞动力学模型,其形式如下
公式(12)中,J1表示内毂的转动惯量;J2表示摩擦片的转动惯量;为内毂角加速度;为摩擦片角加速度;α表示内毂输入角加速度;Ftotal为内毂与摩擦片每时刻所有齿的总冲击载荷;Rb为基圆半径;z为齿数,即轮齿的总数;
S603、对S602中构建的内毂与摩擦片二自由度冲击碰撞动力学模型在时间t内迭代计算,输出每时刻内毂与摩擦片的冲击载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210857772.1A CN115098977B (zh) | 2022-07-20 | 2022-07-20 | 一种浮动支撑摩擦片组件冲击动载计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210857772.1A CN115098977B (zh) | 2022-07-20 | 2022-07-20 | 一种浮动支撑摩擦片组件冲击动载计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115098977A true CN115098977A (zh) | 2022-09-23 |
CN115098977B CN115098977B (zh) | 2024-07-05 |
Family
ID=83299671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210857772.1A Active CN115098977B (zh) | 2022-07-20 | 2022-07-20 | 一种浮动支撑摩擦片组件冲击动载计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115098977B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115098975A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4197931A (en) * | 1977-07-26 | 1980-04-15 | Norcia John A | Clutch friction plate and damper spring construction |
IT8312573A0 (it) * | 1982-07-12 | 1983-07-11 | Valeo | Dispositivo ammortizzatore di torsione in particolare frizione di innesto per veicoli automobili |
CN101034048A (zh) * | 2006-03-06 | 2007-09-12 | 中国北方车辆研究所 | 摩擦片动态强度试验装置 |
WO2013091954A2 (de) * | 2011-12-23 | 2013-06-27 | Robert Bosch Gmbh | Elektrische maschine mit einer axialen vorspannung zwischen einem wälzlager und einem gehäuseteil, sowie verfahren zum betreiben einer solchen elektrischen maschine |
CN104634676A (zh) * | 2015-02-13 | 2015-05-20 | 中国北方车辆研究所 | 一种控制摩擦片动态试验载荷的方法及实现该方法的系统 |
CN105889412A (zh) * | 2016-05-11 | 2016-08-24 | 安徽江淮汽车股份有限公司 | 一种离合器扭转减震机构及离合器 |
CN108000554A (zh) * | 2017-11-30 | 2018-05-08 | 中国矿业大学 | 一种基于片簧的可变刚度柔性关节及其控制方法 |
US20190024736A1 (en) * | 2016-02-18 | 2019-01-24 | Continental Teves Ag & Co. Ohg | Friction lining assembly having a restoring spring for limiting clearance for a motor-vehicle partially lined disk brake |
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115962240A (zh) * | 2022-10-28 | 2023-04-14 | 中国北方车辆研究所 | 一种具有齿部抗冲击结构的浮动支撑摩擦片 |
-
2022
- 2022-07-20 CN CN202210857772.1A patent/CN115098977B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4197931A (en) * | 1977-07-26 | 1980-04-15 | Norcia John A | Clutch friction plate and damper spring construction |
IT8312573A0 (it) * | 1982-07-12 | 1983-07-11 | Valeo | Dispositivo ammortizzatore di torsione in particolare frizione di innesto per veicoli automobili |
CN101034048A (zh) * | 2006-03-06 | 2007-09-12 | 中国北方车辆研究所 | 摩擦片动态强度试验装置 |
WO2013091954A2 (de) * | 2011-12-23 | 2013-06-27 | Robert Bosch Gmbh | Elektrische maschine mit einer axialen vorspannung zwischen einem wälzlager und einem gehäuseteil, sowie verfahren zum betreiben einer solchen elektrischen maschine |
CN104634676A (zh) * | 2015-02-13 | 2015-05-20 | 中国北方车辆研究所 | 一种控制摩擦片动态试验载荷的方法及实现该方法的系统 |
US20190024736A1 (en) * | 2016-02-18 | 2019-01-24 | Continental Teves Ag & Co. Ohg | Friction lining assembly having a restoring spring for limiting clearance for a motor-vehicle partially lined disk brake |
CN105889412A (zh) * | 2016-05-11 | 2016-08-24 | 安徽江淮汽车股份有限公司 | 一种离合器扭转减震机构及离合器 |
CN108000554A (zh) * | 2017-11-30 | 2018-05-08 | 中国矿业大学 | 一种基于片簧的可变刚度柔性关节及其控制方法 |
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115962240A (zh) * | 2022-10-28 | 2023-04-14 | 中国北方车辆研究所 | 一种具有齿部抗冲击结构的浮动支撑摩擦片 |
Non-Patent Citations (3)
Title |
---|
NING K: "Impacting load control of floating supported friction plate and its experimental verification", JOURNAL OF PHYSICS: CONFERENCE SERIES. IOP PUBLISHING, vol. 842, no. 1, 1 January 2017 (2017-01-01) * |
ZHANG C: "Dynamics modeling and analysis of the multistage planetary gear set-bearing-rotor-clutch coupling system considering the tooth impacts of clutches", MECHANICAL SYSTEMS AND SIGNAL PROCESSING, 1 January 2024 (2024-01-01) * |
谢陆: "基于SIMPACK仿真的盘式制动器齿部冲击力影响分析", 设备管理与维修, no. 23, 12 December 2022 (2022-12-12) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115098976A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
CN115098975A (zh) * | 2022-07-20 | 2022-09-23 | 重庆大学 | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 |
CN115098976B (zh) * | 2022-07-20 | 2024-05-28 | 重庆大学 | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115098977B (zh) | 2024-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115098977A (zh) | 一种浮动支撑摩擦片组件冲击动载计算方法 | |
CN115098976B (zh) | 一种考虑自偏心的浮动齿圈摩擦片组件冲击力计算方法 | |
CN105677980B (zh) | 汽车驱动桥系统的模态综合动力学建模与分析方法 | |
CN115130247A (zh) | 一种具有轮齿裂纹的浮动支撑摩擦片冲击动载计算方法 | |
CN115203957A (zh) | 一种考虑轮齿碰撞激励下的离合器摩擦冲击载荷计算方法 | |
CN112836319B (zh) | 一种考虑非均匀分布齿根裂纹故障的仿真方法 | |
CN115618688A (zh) | 高速薄壁齿轮传动系统切片耦合动力学建模方法 | |
CN115098975B (zh) | 一种考虑随机齿向误差的离合器摩擦片冲击力计算方法 | |
CN107341313B (zh) | 基于adams的行星轮系非线性动力学建模方法 | |
Okura et al. | Dynamic analysis of blade reversal behavior in a windshield wiper system | |
CN104590558B (zh) | 一种压电复合材料直升机桨叶结构的控制方法 | |
CN106202676A (zh) | 大功率钻井泵双侧斜齿轮啮合传动系统的受力分析方法 | |
CN104331556A (zh) | 基于车辆平顺性非线性仿真模型的仿真试验方法 | |
Zhang et al. | Contact Mechanics Analysis and Optimization of Shape Modification of Electric Vehicle Gearbox. | |
Yang et al. | Dynamics of vehicle-pavement coupled system based on a revised flexible roller contact tire model | |
CN102221402A (zh) | 低噪声微面汽车后驱动桥设计方法 | |
CN115935694A (zh) | 一种含转动副间隙的微小齿轮动力学建模方法 | |
CN115099035B (zh) | 随机位移激励下含负刚度和惯容协作的悬架减振设计方法 | |
CN106055749A (zh) | 一种提高含间隙连杆机构运动稳定性的方法 | |
Chen et al. | Residual forming effects on full vehicle frontal impact and body-in-white durability analyses | |
Abdullah et al. | Optimization of shape and design parameters of the rigid clutch disc using FEM | |
Pingale et al. | Structure design and development of engine crankshaft damper | |
CN1815229A (zh) | 一种预测汽车蒙皮表面阻尼处理结构声辐射的方法 | |
CN113010975B (zh) | 一种综合考虑加工成本与运动平稳性的齿轮间隙优化设计方法 | |
CN103793564A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |