首页 / 专利库 / 电子零件及设备 / 厚度剪切振动 / 一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法

一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法

阅读:922发布:2020-05-11

专利汇可以提供一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种抗 颤振 约束阻尼型 复合材料 镗杆的最优厚度确定方法,所述方法包括:基于Euler-Bernoulli梁理论,建立约束阻尼型复合材料镗杆的振动微分方程;确定约束阻尼型复合材料镗杆振型函数与固有 频率 ;建立考虑材料阻尼时的镗杆振动微分方程;确定结构阻尼比;确定约束阻尼型复合材料镗杆临界 主轴 转速与切削宽度;确定各材料层对应的最优厚度。本发明基于Euler-Bernoulli梁理论,建立了具有约束阻尼层的复合材料镗杆的再生颤振线性动 力 学模型,构建了系统临界稳定状态下主轴转速切削深度的表达式,并且确定了各材料层最优厚度的计算方法。,下面是一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法专利的具体信息内容。

1.一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于:所述方法以下步骤:
a.基于Euler-Bernoulli梁理论,建立约束阻尼型复合材料镗杆的振动微分方程;
b.确定约束阻尼型复合材料镗杆振型函数与固有频率
c.建立考虑材料阻尼时的镗杆振动微分方程;
d.确定结构阻尼比;
e.确定约束阻尼型复合材料镗杆临界主轴转速与切削宽度;
f.确定各材料层对应的最优厚度。
2.根据权利要求1所述的一一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于,上述步骤a中的振动微分方程通过式(1)确定:
式中,y是镗杆上距原点x处的截面在t时刻的位移(mm);(EJ)equiv是约束阻尼型结构镗杆的等效弯曲刚度;(ρA)equiv是约束阻尼型结构镗杆的等效面密度;L是镗杆的悬伸长度;δ是单位脉冲函数;
其中
(EJ)equiv=(EJ)1+(EJ)2+(EJ)3                   (2)
(ρA)equiv=(ρA)1+(ρA)2+(ρA)3                   (3)
式中,(EJ)1、(EJ)2、(EJ)3分别表示为基体、阻尼层和约束层的弯曲刚度,(ρA)1、(ρA)2、(ρA)3分别表示为基体、阻尼层和约束层的面密度;
根据材料学理论,可得约束层和阻尼层的弯曲刚度与面密度为
复合材料基体的弯曲刚度和质量分别为
其中, 是基体复合材料镗杆的偏轴刚度系数。
3.根据权利要求1所述的一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于,所述步骤b中的振型函数与固有频率确定如下:
镗杆的主振动表示为
将式(7)代入式(1)中并根据分离变量法可得振型函数为
Y(x)=C1sinβx+C2cosβx+C3shβx+C4chβx             (8)
其中,β为
悬臂梁的边界条件为
将边界条件代入式(8),得
cosβLchβL=-1                        (11)
解式(11)可得
βL=1.8751                          (13)
由式(9)与式(13)可得镗杆的固有频率为
由式(7)并根据主振函数的正交性,式(1)可化为
式中,Mi和Ki计算公式分别如下
Mi和Ki分别为模态质量和模态刚度。
4.根据权利要求1所述的一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于,所述步骤c中的考虑材料阻尼时的镗杆振动微分方程确定如下:
在考虑材料阻尼的影响(材料的阻尼并不影响材料本身的固有频率和振型)的情况下,设结构的阻尼比为ξi,则式(15)可变为
5.根据权利要求1所述的一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于,所述步骤d中的结构阻尼比确定如下:
结构阻尼比ξi为
根据振动力学理论,镗杆的阻尼系数C可表示为
式中,m为镗杆的质量,k为约束阻尼层复合镗杆的静刚度,η为约束阻尼型复合材料镗杆结构损耗因子;
k计算公式如下
式中,D1为基体直径,D2为阻尼层外圆直径,D3为约束层外圆直径;
η计算公式如下
式中,μ表示阻尼层材料的损耗因子,X是剪切参数,Y为刚度参数;
X表达式为
Y表达式为
6.根据权利要求1所述的一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于,所述步骤e中的临界主轴转速与切削宽度确定如下:
假定只考虑切削深度对切削力的影响,由于切削力的大小与动态切削深度有关,而后者的变化仅取决于再生效应,则式(16)可写为
其中:Kc是进给方向上单位切削深度上的切削力系数,b是切削深度;
在切削过程中,通常是较低阶模态被激发,颤振频率略高于系统的一阶固有频率,故在颤振分析中,只考虑镗杆的基频;
于是,具有约束型阻尼结构的复合材料镗杆的动力学方程为
对式(27)进行拉氏变换,得到
由控制工程基本原理可知,振动系统的稳定性取决于特征方程根s的性质,s可以写成s=σ+iω的形式,且切削系统稳定性与s的实部σ的取值有关:如果σ>0,系统是不稳定的;如果σ<0,系统是稳定的;当σ=0(此时ω=ωc),系统是临界稳定的;
考虑在临界稳定状态下,即σ=0,s=iωc,代入式(28)得
根据欧拉公式e-ix=cos(x)-isin(x)将上式展开,并分离实部与虚部,得
2M1ξ1ω1ωc=KcbY12(L)sin(Tωc)                 (31)
由式(30)和式(31)得
其中
得到系统临界稳定状态下主轴转速n和切削深度blim计算公式如下
7.根据权利要求1所述的一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,其特征在于,所述步骤f中的各材料层对应的最优厚度确定如下:
在刀具前和被加工材料相同的情况下,最大切削深度与结构的阻尼比和樘杆的静弯曲刚度的乘积成正比,即与动刚度的大小成正比;
动刚度的大小可由公式(36)计算
在给出镗杆整体半径尺寸及在给出R3的情况下,确定基体层R1和阻尼层R2尺寸变化范围,可将各层厚度优化问题转化为求解max{Kd}在边界(37)上的最优值问题;

说明书全文

一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法

技术领域

[0001] 本发明涉及一种镗杆颤振稳定性分析方法,具体涉及一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法。

背景技术

[0002] 镗杆是深孔加工镗削工艺中的主要承载部件之一,多用于工件内孔的成型加工在镗削加工中,切削多为不均匀力,而镗杆的刚度有限,使得加工过程中镗杆很容易产生振动,导致工件的加工质量精度难以得到保证。
[0003] 对于大长径比约束阻尼型复合材料镗杆,传统的镗杆颤振稳定性分析方法,很少同时考虑约束层、阻尼型镗杆与复合材料镗杆共同作用对颤振稳定性的影响,并且对于约束阻尼型复合材料镗杆缺少一种计算各材料层最优厚度的计算方法。

发明内容

[0004] 本发明的目的是提供一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法。该方法为大长径比约束阻尼型复合材料镗杆规格的选取和整杆设计提供参考与依据。
[0005] 为实现上述目的,本发明所采用的技术方案为:
[0006] 本发明涉及一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,所述方法以下步骤:
[0007] a.基于Euler-Bernoulli梁理论,建立约束阻尼型复合材料镗杆的振动微分方程;
[0008] b.确定约束阻尼型复合材料镗杆振型函数与固有频率
[0009] c.建立考虑材料阻尼时的镗杆振动微分方程;
[0010] d.确定结构阻尼比;
[0011] e.确定约束阻尼型复合材料镗杆临界主轴转速与切削宽度;
[0012] f.确定各材料层对应的最优厚度;
[0013] 所述步骤a中的振动微分方程通过式(1)确定:
[0014]
[0015] 式中,y是镗杆上距原点x处的截面在t时刻的位移(mm);(EJ)equiv是约束阻尼型结构镗杆的等效弯曲刚度;(ρA)equiv是约束阻尼型结构镗杆的等效面密度;L是镗杆的悬伸长度;δ是单位脉冲函数。
[0016] 其中
[0017] (EJ)equiv=(EJ)1+(EJ)2+(EJ)3                  (2)
[0018] (ρA)equiv=(ρA)1+(ρA)2+(ρA)3                  (3)
[0019] 式中,(EJ)1、(EJ)2、(EJ)3分别表示为基体、阻尼层和约束层的弯曲刚度,(ρA)1、(ρA)2、(ρA)3分别表示为基体、阻尼层和约束层的面密度。
[0020] 复合材料基体的弯曲刚度和质量分别为
[0021]
[0022] 其中, 是基体复合材料镗杆的偏轴刚度系数。
[0023] 根据材料力学理论,可得约束层和阻尼层的弯曲刚度与面密度为
[0024]
[0025]
[0026] 所述步骤b中的振型函数与固有频率确定如下:
[0027] 镗杆的主振动表示为
[0028] 将式(7)代入式(1)中并根据分离变量法可得振型函数为
[0029] Y(x)=C1 sinβx+C2 cosβx+C3 shβx+C4 chβx            (8)
[0030] 其中,β为
[0031]
[0032] 悬臂梁的边界条件为
[0033]
[0034] 将边界条件代入式(8),得
[0035] cosβLchβL=-1                       (11)
[0036]
[0037] 解式(11)可得
[0038] βL=1.8751                         (13)
[0039] 由式(9)与式(13)可得镗杆的固有频率为
[0040]
[0041] 由式(7)并根据主振函数的正交性,式(1)可化为
[0042]
[0043]
[0044] 式中,Mi和Ki分别为模态质量和模态刚度,计算公式分别如下
[0045]
[0046]
[0047] 所述步骤c中的考虑材料阻尼时的镗杆振动微分方程确定如下:
[0048] 在考虑材料阻尼的影响(材料的阻尼并不影响材料本身的固有频率和振型)的情况下,设结构的阻尼比为ξi,则式(15)可变为
[0049]
[0050] 所述步骤d中的结构阻尼比确定如下:
[0051] 结构阻尼比ξi为
[0052]
[0053] 根据振动力学理论,镗杆的阻尼系数C可表示为
[0054]
[0055] 式中,m为镗杆的质量,k为约束阻尼层复合镗杆的静刚度,η为约束阻尼型复合材料镗杆结构损耗因子。
[0056] k计算公式如下
[0057]
[0058] 式中,D1为基体直径,D2为阻尼层外圆直径,D3为约束层外圆直径。
[0059] η计算公式如下
[0060]
[0061] 式中,μ表示阻尼层材料的损耗因子,X是剪切参数,Y为刚度参数。
[0062] X表达式为
[0063]
[0064] Y表达式为
[0065]
[0066] 所述步骤e中的临界主轴转速与切削宽度确定如下:
[0067] 假定只考虑切削深度对切削力的影响。由于切削力的大小与动态切削深度有关,而后者的变化仅取决于再生效应,则式(16)可写为:
[0068]
[0069] 其中,Kc是进给方向上单位切削深度上的切削力系数,b是切削深度。
[0070] 在切削过程中,通常是较低阶模态被激发,颤振频率略高于系统的一阶固有频率,故在颤振分析中,只考虑镗杆的基频。
[0071] 于是,具有约束型阻尼结构的复合材料镗杆的动力学方程为
[0072]
[0073] 对式(27)进行拉氏变换,得到
[0074]
[0075] 由控制工程基本原理可知,振动系统的稳定性取决于特征方程根s的性质,s可以写成s=σ+iω的形式,且切削系统稳定性与s的实部σ的取值有关:如果σ>0,系统是不稳定的;如果σ<0,系统是稳定的;当σ=0(此时ω=ωc),系统是临界稳定的。
[0076] 考虑在临界稳定状态下,即σ=0,s=iωc,代入式(28)得
[0077]
[0078] 根据欧拉公式e-ix=cos(x)-isin(x)将上式展开,并分离实部与虚部,得[0079]
[0080] 2M1ξ1ω1ωc=KcbY12(L)sin(Tωc)                 (31)
[0081] 由式(30)和式(31)得
[0082]
[0083] 其中
[0084]
[0085] 得到系统临界稳定状态下主轴转速n和切削深度blim计算公式如下
[0086]
[0087]
[0088] 所述步骤f中的各材料层对应的最优厚度确定如下:
[0089] 在刀具前和被加工材料相同的情况下,最大切削深度与结构的阻尼比和樘杆的静弯曲刚度的乘积成正比,即与动刚度的大小成正比。
[0090] 动刚度的大小可由公式(36)计算
[0091]
[0092] 在给出镗杆整体半径尺寸及在给出R3的情况下,确定基体层R1和阻尼层R2尺寸变化范围,可将各层厚度优化问题转化为求解max{Kd}在边界(37)上的最优值问题。
[0093]
[0094] 本发明基于Euler-Bernoulli梁理论,建立了具有约束阻尼层的复合材料镗杆的再生颤振线性动力学模型,构建了系统临界稳定状态下主轴转速n和切削深度blim的表达式,并且确定了各材料层最优厚度的计算方法。附图说明
[0095] 图1为本发明的约束阻尼型复合材料镗杆的结构示意图。
[0096] 图2为图1在径向方向上的截面图。
[0097] 图3为本发明中不同材料层厚度对应的颤振稳定性叶瓣曲线图。

具体实施方式

[0098] 下面结合具体实施例对本发明的实施方式进一步说明。
[0099] 约束阻尼型复合材料镗杆的结构如图1至图2所示,包括由内到外依次为基体层3、阻尼层2和约束层1。
[0100] 本发明涉及一种抗颤振约束阻尼型复合材料镗杆的最优厚度确定方法,所述方法以下步骤:
[0101] a.基于Euler-Bernoulli梁理论,建立约束阻尼型复合材料镗杆的振动微分方程;
[0102] b.确定约束阻尼型复合材料镗杆振型函数与固有频率;
[0103] c.建立考虑材料阻尼时的镗杆振动微分方程;
[0104] d.确定结构阻尼比;
[0105] e.确定约束阻尼型复合材料镗杆临界主轴转速与切削宽度;
[0106] f.确定各材料层对应的最优厚度;
[0107] 所述步骤a中的振动微分方程通过式(1)确定:
[0108]
[0109] 式中,y是镗杆上距原点x处的截面在t时刻的位移(mm);(EJ)equiv是约束阻尼结构镗杆的等效弯曲刚度;(ρA)equiv是约束阻尼结构镗杆的等效面密度;L是镗杆的悬伸长度;δ是单位脉冲函数。
[0110] 其中
[0111] (EJ)equiv=(EJ)1+(EJ)2+(EJ)3                  (2)
[0112] (ρA)equiv=(ρA)1+(ρA)2+(ρA)3                   (3)
[0113] 式中,(EJ)1、(EJ)2、(EJ)3分别表示为基体、阻尼层和约束层的弯曲刚度,(ρA)1、(ρA)2、(ρA)3分别表示为基体、阻尼层和约束层的面密度。
[0114] 复合材料基体的弯曲刚度和质量分别为:
[0115]
[0116] 其中, 是基体复合材料镗杆的偏轴刚度系数。
[0117] 根据材料力学理论,可得约束层和阻尼层的弯曲刚度与面密度为:
[0118]
[0119]
[0120] 所述步骤b中的振型函数与固有频率确定如下:
[0121] 镗杆的主振动表示为
[0122]
[0123] 将式(7)代入式(1)中并根据分离变量法可得振型函数为
[0124] Y(x)=C1 sinβx+C2 cosβx+C3 shβx+C4 chβx            (8)
[0125] 其中,β为
[0126]
[0127] 悬臂梁的边界条件为
[0128]
[0129] 将边界条件代入式(8),得
[0130] cosβLchβL=-1                       (11)
[0131]
[0132] 解式(11)可得
[0133] βL=1.8751                         (13)
[0134] 由式(9)与式(13)可得镗杆的固有频率为
[0135]
[0136] 由式(7)并根据主振函数的正交性,式(1)可化为
[0137]
[0138]
[0139] 式中,Mi和Ki分别为模态质量和模态刚度,计算公式分别如下
[0140]
[0141]
[0142] 所述步骤c中的考虑材料阻尼时的镗杆振动微分方程确定如下:
[0143] 在考虑材料阻尼的影响(材料的阻尼并不影响材料本身的固有频率和振型)的情况下,设结构的阻尼比为ξi,则式(15)可变为
[0144]
[0145] 所述步骤d中的结构阻尼比确定如下:
[0146] 结构阻尼比ξi为
[0147]
[0148] 根据振动力学理论,镗杆的阻尼系数C可表示为
[0149]
[0150] 式中,m为镗杆的质量,k为约束阻尼层复合镗杆的静刚度,η为约束阻尼型复合材料镗杆结构损耗因子。
[0151] k计算公式如下
[0152]
[0153] 式中,D1为基体直径,D2为阻尼层外圆直径,D3为约束层外圆直径。
[0154] η计算公式如下
[0155]
[0156] 式中,μ表示阻尼层材料的损耗因子,X是剪切参数,Y为刚度参数。
[0157] X表达式为
[0158]
[0159] Y表达式为
[0160]
[0161] 所述步骤e中的临界主轴转速与切削宽度确定如下:
[0162] 假定只考虑切削深度对切削力的影响。由于切削力的大小与动态切削深度有关,而后者的变化仅取决于再生效应,则式(16)可写为:
[0163]
[0164] 其中,Kc是进给方向上单位切削深度上的切削力系数,b是切削深度。
[0165] 在切削过程中,通常是较低阶模态被激发,颤振频率略高于系统的一阶固有频率,故在颤振分析中,只考虑镗杆的基频。
[0166] 于是,具有约束型阻尼结构的复合材料镗杆的动力学方程为
[0167]
[0168] 对式(27)进行拉氏变换,得到
[0169]
[0170] 由控制工程基本原理可知,振动系统的稳定性取决于特征方程根s的性质,s可以写成s=σ+iω的形式,且切削系统稳定性与s的实部σ的取值有关:如果σ>0,系统是不稳定的;如果σ<0,系统是稳定的;当σ=0(此时ω=ωc),系统是临界稳定的。
[0171] 考虑在临界稳定状态下,即σ=0,s=iωc,代入式(28)得
[0172]
[0173] 根据欧拉公式e-ix=cos(x)-isin(x)将上式展开,并分离实部与虚部,得[0174]
[0175] 2M1ξ1ω1ωc=KcbY12(L)sin(Tωc)                 (31)
[0176] 由式(30)和式(31)得
[0177]
[0178] 其中
[0179]
[0180] 得到系统临界稳定状态下主轴转速n和切削深度blim计算公式如下
[0181]
[0182]
[0183] 所述步骤f中的各材料层对应的最优厚度确定如下:
[0184] 在刀具前角和被加工材料相同的情况下,最大切削深度与结构的阻尼比和樘杆的静弯曲刚度的乘积成正比,即与动刚度的大小成正比。
[0185] 动刚度的大小可由公式(36)计算
[0186]
[0187] 在给出镗杆整体半径尺寸及在给出R3的情况下,确定基体层R1和阻尼层R2尺寸变化范围,可将各层厚度优化问题转化为求解max{Kd}在边界(37)上的最优值问题。
[0188]
[0189] 优选地,所述基体层材料选用/环树脂
[0190] 优选地,所述阻尼层材料选用聚四氟乙烯(Teflon)。
[0191] 优选地,所述约束层材料选用YG20C。
[0192] 优选地,所述L选用L=100mm。
[0193] 优选地,所述R3选用R3=1.5mm。
[0194] 求解各层材料厚度最优值,可以得到R1=0.4828mm,R2=1.0172mm。即H1=0.4828mm,H2=0.5344mm,H3=0.4828mm,转化为比例表示:H1:H2:H3=1:1.2:1,如图3所示,稳定性叶瓣曲线在坐标平面上的位置表现为先上升,后下降的趋势,说明各层材料之间存在着某个厚度值使镗杆的颤振稳定性达到最佳并且验证了优化方法的可行性。
[0195] 最后应该说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本权利要求范围当中。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈