《1 前言》
1 前言
液体静压支承动态性能是指支承系统在所加动态干扰力作用下所产生的动力学响应。由于液体静压支承是一个具有良好粘性阻尼、质量与刚性的液体与机械系统, 因此, 分析与研究这种系统动态性能应以传递函数分析法最为简便[1]。液体静压支承动态干扰力实际上只有两种:一种是脉冲与阶跃瞬态干扰力, 另一种是稳态循环周期干扰力。因此, 它们的动力学响应就是具有粘性阻尼的自由振动与循环受迫振动的响应[2]。
分析与研究液体静压支承系统动态性能, 应以支承系统在外加动态干扰力作用下能否正常工作为唯一标准。因此, 首先应确定它是否为一个稳定系统;二是系统在外加瞬态干扰力作用下所引起的振动衰减趋于稳定状态所需的时间;三是系统在稳态周期干扰力作用下的系统动刚度, 以及自身的谐振频率, 这种频率一般称之为固有频率或共振频率。
瞬态干扰力在机床精密加工中很少出现, 但它的响应时间却能反映支承系统的抗瞬态干扰能力。衰减很慢的稳定系统并不能保证设备能够很好地正常工作, 因此, 支承系统对瞬态干扰力引起的衰减时间就成为检查系统抗瞬态干扰能力的一种手段。
稳态周期干扰力在任何机床设备上均存在。例如, 主轴驱动皮带的抖动频率, 滚动轴承滚体波动频率, 齿轮传动的轮齿啮合频率, 油泵供油的脉动频率以及主轴自身的固有频率[3]。尽管这些周期干扰力在正常工作条件下不会引起大的振动, 但支承系统固有频率与这些干扰频率接近时会引起强烈的共振, 导致机床设备不能正常工作, 甚至损坏某些构件[2]。因此, 在设计液体静压支承时, 不能不给予特别重视。
新型可变节流器具有优良的节流性能, 采用最佳节流参数可保持支承间隙不随外加载荷的增加而减小, 也就是能够保持间隙近于不变。但是, 新型可变节流器具有弹性环节, 因而不论采用传统的还是新型的可变节流液体静压支承系统, 均存在系统固有频率。只要系统固有频率与机床原已存在的稳定周期干扰频率相差20 % 以上, 就不会引起共振。固定节流液体静压支承由于节流器中不存在弹性环节, 故属于恒稳系统, 也没有系统固有频率存在。但是固定节流液体静压支承的支承间隙随外加载荷的增加而剧烈减小, 所以其静态性能不佳。
新型可变节流液体静压支承由于节流器中具有弹性环节, 与固定节流液体静压支承有明显的差异, 因此, 笔者在所建立的数学模型——传递函数中特别包含有公式 (2) 。国内外资料讨论的液体静压支承动态性能均未涉及节流器弹性环节对系统动态性能的严重影响[4,5], 因此, 其结论也会与笔者在文中讨论的结果有较大差别。采用可变节流液体静压支承系统而忽略弹性环节对系统动态性能的影响, 将在应用时产生严重后果, 这一点不能不予以特别注意。液体静压支承系统静态性能保证加工精度, 而系统动态性能保证稳定条件, 二者不能偏废, 否则很难获得良好的工作效果。
《2 液体静压支承系统数学摸型》
2 液体静压支承系统数学摸型
1) 根据静压支承系统的工作原理[6], 由油垫力平衡条件
m1y¨+c1y˙=f(t)−Aepr, (1)
式中m1为上支承质量, y为m的位移, c1为阻尼系数, f (t) 为动载荷, Ae为油腔有效面积, pr为油腔内油的压强。
2) 根据阀芯力平衡条件
m2x¨+c2x˙+k2x=Ack2, (2)
式中m2为阀芯质量, x为阀芯位移, c2为阀芯运动时的粘性阻尼系数, k2为弹簧系数, Ac为阀芯横截面积。
3) 根据支承油囊流量连续条件
q=Q−Aey¨+Acx˙+βvp˙r, (3)
式中q为节流器供油流量, Q为油腔排油流量, βv为压力变化系数。
这里 kx=∂q∂x∣∣∣pr=pr0x=0, ky=∂Q∂y∣∣∣pr=pr0y=0,
k′p=∂q∂pr∣∣∣pr=pr0x=0,k′′p=∂Q∂pr∣∣∣pr=pr0y=0,
将kx, ky, k′p与k″p及其正负号均代入式 (3) , 并对式 (1) 至式 (3) 进行拉氏变换得
m1s2Y(s)+ceYs(s)=F(s)−Aepr(s)‚ (4)m2s2X(s)+c2sX(s)+k2X(s)=Acpr(s), (5)(kx−Acs)X(s)+(ky+Acs)Y(s)=(kp+βvs)pr(s)。 (6)
将式 (4) 至式 (6) 合并, 经过适当运算可得支承系统的闭环传递函数
G(s)=D6s3+D7s2+D8s+D9D0s5+D1s4+D2s3+D3s2+D4s+D5, (7)
D0=βVm1m2,D1=kpm1m2+βV(m1c2+m2c1),D2=βV(m1k2+c1c2)+kp(m1k2+m2c1)+A2em2+A2cm1,D3=βVc1k2+kp(m1k2+c1c2)+A2ec2+A2cc1+kyAem2−kxAcm1,D4=kpc1k2+kyAec2+A2ek2−kxAcc1,D5=kyAek2,D6=βVm2,D7=kpm2+βVc2,D8=kpc2+A2c+βVk2,D9=kpk2−kxAc。
《3 液体静压支承系统动态性能新表达式》
3 液体静压支承系统动态性能新表达式
《3.1 系统稳定性判别式——霍维茨判据》
3.1 系统稳定性判别式——霍维茨判据
∣∣∣∣∣∣D4D500D2D3D4D5D0D1D2D300D0D2∣∣∣∣∣∣>0 (8)
《3.2 系统抗瞬态干扰能力》
3.2 系统抗瞬态干扰能力
G(s)=Y(s)F(s),Y(s)=G(s)⋅F(s)。
令F (s) 为单位脉冲力δ (t) , 则拉氏变换后有L[F (t) ]=L[δ (t) ]=1, 由此
y(t)=L−1[G(s)⋅F(s)]=L−1[G(s)],
已知G (s) 为有理函数, 故可用分式表达。当G (s) 分母特征方程有5个独立根时,
G(s)=∑i=15Ais+Bi, (9)
f(t)=∑i=15Aie−Bit。 (10)
G(s)=∑i=13Ais+Bi+s+a1(s+a)2+b2,
f(t)=∑i=13Aie−Bit+1b[(a1+a)2+b2]12e−atsin(bt+Ψ), (11)
这里Ψ=tan−1ba1+a。 (12)
《3.3 支承系统动刚度》
3.3 支承系统动刚度
G(s)=Y(s)F(s), 所以 J(s)=F(s)Y(s),
J(s)=F(s)Y(s)=D0s5+D1s4+D2s3+D3s2+D4s+D5D6s3+D7s2+D8s+D9。 (13)
假定液体静压支承系统中不存在保有少量空气的“死腔”, 而油液中又未溶有大量空气, 则βV=0 。一般节流器阀芯或薄板很小, 故其质量m2≈0, 由此令D0=D1=D6=D7=0, 这样, 式 (13) 可简化为三次方程
J(s)=F(s)Y(s)=D0s3+D3s2+D4s+D5D8s+D9, (14)
令s=jω并代入式 (14) , 可得支承系统动刚度计算式
J(ω)=A2+B2ω2−−−−−−−−−√‚ (15)A=D5D9+(D4D8−D3D9)ω2−D2D8ω4D28ω2+D29 (16)B=(D5D8−D4D9)−(D3D8−D2D4)ω2D28ω2+D29。 (17)
J(ω=0)=D5D9。 (18)
《3.4 液体静压支承系统固有频率》
3.4 液体静压支承系统固有频率
根据液体静压支承静态性能表达式[8], 当采用新型节流器节流时, 令WxT=0, 则J (B) →∞, 据此可得D9 = 0 。
J(ω)=(D4D8−D2D8)2+(D5D8ω2−D3D8)2ω2−−−−−−−−−−−−−−−−−−−−−−−−−√。 (19)
ω6+(D23−2D4D2D22)ω4−D252D22=0。 (20)
令Ω=ω2‚a=D23−2D4D2D22‚c=−D252D22,
令Ω=y−13a, 代入式 (21) , 得典型三次方程
y3+py+q=0, (22)
这里p=−13a3,q=227a3+c。
若Δ=(13p)3+(12q)2=127a3c+14c2>0, 则式 (22) 中有一个实根, 这个实根即系统固有频率ω*。
已知Ω=ω2=y−12a‚y=−q2+Δ−−√−−−−−−−−√3+−q2−Δ−−√−−−−−−−−√3,
ω∗=−q2+Δ−−√−−−−−−−−√3+−q2+Δ−−√−a3−−−−−−−−−−−−√3−−−−−−−−−−−−−−−−−−−−−−−−−−√。 (23)
新型节流器中有弹性环节, 故支承系统必然存在固有频率。应当指出:在设计液体静压支承时, 应当保证系统固有频率与机床设备上原已存在的各种周期干扰频率相差20 %以上, 方能保证不会引起共振。
对于传统固定节流器, 如小孔与毛细管节流, 由于节流器中不包含弹性环节, 故支承系统没有固有频率。因此, 固定节流支承系统动态性能良好, 但静态性能不佳, 而静态性能却是保证加工精度的条件, 因而应用时应当考虑这些特点。
《4 液体静压支承系统动态性能实验验证》
4 液体静压支承系统动态性能实验验证
《4.1 试验台结构与工作原理》
4.1 试验台结构与工作原理
试验台采用新型变径毛细管节流单油囊支承, 其工作原理可参阅文献[6]。应当说明的是, 为了实验更具有普遍性, 参数未按D9= 0选择, 故其固有频率不能用式 (23) 计算。
《4.2 试验台支承参数》
4.2 试验台支承参数
W0=705 N, k2=15.7 N/mm, Ae=12 780 mm2, Ac=154 mm2, kBL=30.8, h0=35×10-3 mm, hc=0.012 5 mm, γe =2.727 0, Le=60 mm, kP=2 958 mm5/N·s, kx=2 598 mm2/s, ky=12 226 mm2/s, c1=226.26 N·s /mm, c2=0.044 N·s /mm , m1=0.07 N·s2 /mm, m2=1.22 ×10-4 N·s2/mm, μ50℃=8.5×10-8 N·s /mm2 , ps=0.44 MPa, pr0=0.055 MPa 。
《4.3 支承系统动态系数》
4.3 支承系统动态系数
D2=1.79×104 N·s ·mm2 , D3=1 260.12×104 N·s ·mm3 , D4=26 345.61×104 N·mm3 , D5=24 059.10×104 mm3/s, D8=2.38×104 mm4 , D9=0.065×104 mm4/s。
《4.4 系统稳定性判别》
4.4 系统稳定性判别
《4.5 系统抗瞬态干扰能力》
4.5 系统抗瞬态干扰能力
G(s)=2.38s+0.0651.79s3+1260.12s2+26345.61s+24059.10=1.33{s+0.0115s3+703.97s2+14718.22s+13723.52},
解分母特征方程得三个负实根:0.92 , 682.44与 20.59 , 由此令
G(s)=1.33{As+0.92+Bs+6.82.44+Cδ+20.59},
A=3.11×10−3‚B=−1.52×10−3,C=−1.59×10−3‚所以G(s)=4.14×10−3{1s+0.92−0.49s+682.44−0.51s+20.59}。
《4.6 试验台液体静压支承系统动刚度》
4.6 试验台液体静压支承系统动刚度
J(ω)=[(279.94+11063.87ω2−0.85ω4ω2)2+(100029.29−529.84ω2ω2)2ω2]12‚
《表1》
ω/Hz | 0 | 1 | 4 | 4.5 | 5 | 100 | 1000 |
J (ω) /×104N/mm | ∞ | 1.499 | 1.107 | 1.106 | 1.107 | 5.295 | 99.39 |
由于支承参数未按静刚度JB→∞条件选择, 故D9≠0, 因而不能按式 (23) 求固有频率, 但从上表可看出, 动刚度最小时的频率即为系统固有频率, 由此
ω∗=4.5Hz,A∗max=11.106×104=0.09×10−3mm/N≈0.1μm/N。
根据文献[6]得知抗瞬态干扰能力实测结果与计算非常近似。由于固有频率太低, 小于5 Hz, A*max≈0.1 μm/N太小, 受条件限制, 未能在0~1 000 Hz 扫频中测出数值来。
《5 结论》
5 结论
1) 液体静压支承系统中若不存有“死腔”, 而油液中又未溶解大量空气, 则支承系统一般是稳定系统。
2) 可变节流支承系统由于节流器中含有弹性环节, 因此有固有频率存在。频率高低与弹性环节刚性有关。因此, 支承系统固有频率必须与机床上原已存在的周期干扰频率相差距20 %以上, 方能保证不会引起共振。
3) 传统固定节流液体静压支承系统由于节流器中不含弹性环节, 故不存在系统固有频率, 而动态性能良好。但由于静态性能不佳[8], 应用效果很不理想。
4) 液体静压支承系统抗瞬态干扰能力良好, 一般衰减至原有状态时不超过1.0 s。
5) 液体静压支承动态性能新表达式物理概念清晰, 公式简单, 计算结果可靠。
6) 液体静压支承系统动态性能保证设备工作条件, 而静态性能保证加工精度, 静态性能不佳很难获得令人满意的结果。
7) 实验证实, 计算与实测结果很近似, 故动态性能表达式可用于实际。