本文主要研究一个受干摩擦作用的振动系统的周期解的近似求解。此系统的周期解近似地表示Fourier级数,并通过Galerkin法把微分方程转换成非线性代数方程组,利用Broyden法来对此非线性代数方程组求解。最后数值积分的结果用来验证该方法的有效性。 In this paper, the approximate periodic solution of a vibration system subjected to dry friction is studied. The periodic solution of this system approximately represents the Fourier series, the dif-ferential equation is transformed into nonlinear algebraic equations by Galerkin method, and the nonlinear algebraic equations are solved by Broyden method. Finally, numerical integration results are used to verify the effectiveness of the method.
栗鹏
同济大学航空航天与力学学院,上海
收稿日期:2019年11月7日;录用日期:2019年11月22日;发布日期:2019年11月29日
本文主要研究一个受干摩擦作用的振动系统的周期解的近似求解。此系统的周期解近似地表示Fourier级数,并通过Galerkin法把微分方程转换成非线性代数方程组,利用Broyden法来对此非线性代数方程组求解。最后数值积分的结果用来验证该方法的有效性。
关键词 :干摩擦,Fourier级数,Galerkin法,Broyden法
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
众所周知,干摩擦广泛地存在于各类机械系统中,比如刹车系统、离合器、轮胎运动和弓弦乐器等 [
关于干摩擦动力学系统的研究方法主要有:近似解析法 [
本文中所考虑的干摩擦振动系统如图1所示,一个质量为m的物体,通过一个刚度系数k,粘性系数c的弹簧与一个固定端连接,而且该物体与地面之间存在干摩擦力,摩擦系数为 μ ,一个正弦激励u作用在该物体上,设 u = F sin ω t ,F是激励振幅, ω 是激励频率。
图1. 干摩擦振动系统的示意图
根据牛顿第二定律得到:
m x ¨ + c x ˙ + k x = F sin ω t + f (1)
其中f是干摩擦力,形式如下:
f = − μ m g sign ( x ˙ ) (2)
化简得到:
m x ¨ + c x ˙ + k x + μ m g sign ( x ˙ ) = F sin ω t (3)
为了减少系统的参数,引进如下参数:
ω 0 = k m , ξ = c ω 0 2 k , τ = ω t , Ω = ω ω 0 , α = F k , β = μ m g k
则方程(3)变为:
Ω 2 x ¨ + 2 Ω ξ x ˙ + x + β sign ( x ˙ ) = α sin τ (4)
对于强非线性的动力学系统,Galerkin法可以作为对系统的周期解求解的一种近似方法。
设该系统的周期解的Fourier级数为:
x = a 0 + ∑ n = 1 N ( a n cos ( n τ ) + b n sin ( n τ ) ) = V ( τ ) A (5)
其中 V ( τ ) = [ 1 , cos τ , ⋯ , cos N τ , sin τ , ⋯ , sin N τ ] , A = [ a 0 , a 1 , ⋯ , a N , b 1 , ⋯ , b N ] T ;N是Fourier级数的阶数, A 为解的Fourier级数的系数所组成的向量, [ ] T 表示向量的转置。
x ˙ = ( V ( τ ) A ) ′ = V ( τ ) C A (6)
x ¨ = ( V ( τ ) C A ) ′ = V ( τ ) D A (7)
其中 C = [ 0 0 1 × N 0 1 × N 0 N × 1 0 N × N − S N × N 0 N × 1 S N × N 0 N × N ] , D = [ 0 0 1 × N 0 1 × N 0 N × 1 − W N × N 0 N × N 0 N × 1 0 N × N − W N × N ] , S N × N = [ 1 0 ⋱ 0 N ] , W N × N = [ 1 0 ⋱ 0 N 2 ] 。
把(6)和(7)代入方程(4)得到:
r = Ω 2 V ( τ ) D A + 2 Ω ξ V ( τ ) C A + V ( τ ) A + β sign ( V ( τ ) C A ) − α sin τ = 0 (8)
由Galerkin过程得到系统的残差方程:
R = ∫ 0 2 π ( Ω 2 V ( τ ) D A + 2 Ω ξ V ( τ ) C A + V ( τ ) A + β sign ( V ( τ ) C A ) − α sin τ ) V ( τ ) T d τ = 0 (9)
对于一般的Galerkin法来说,求解残差方程一般利用解析的方法,进而得到Fourier级数的系数。但是对于非线性方程组(9)来说,由于符号sign函数的存在,方程(9)中的 a 0 , a 1 , ⋯ , a N , b 1 , ⋯ , b N 并不能解析求出,而且方程(9)对于系数 a 0 , a 1 , ⋯ , a N , b 1 , ⋯ , b N 的导数也不能显示表示,尤其系统中存在stick运动的情形,数值求解中的Newton法也就不再适用于此方程的求解。因此可以利用Broyden法来对方程(9)进行数值求解,迭代过程中系统的Jacobian矩阵可以数值给出。其迭代过程如下:
第一步:给出近似的初值 A 0 和Jacobian矩阵 Β 0 ,并由此初值开始,对系统进行迭代,直到求出满足精度的解;
第二步:假设方程已经进行了第i次迭代,得到方程的解为 A i ,把 A i 代入(5)中,并把(5)在一个周
期内进行离散化,即把 [ 0 , 2 π ] 等分为H个小区间,则 τ j = j Δ τ , j = 1 , 2 , ⋯ , H ,其中 Δ τ = 2 π H , H ≥ 2 N + 1 ,
则对应的向量为:
x ¯ = [ x ( τ 1 ) , x ( τ 2 ) , ⋯ , x ( τ H ) ] T = Γ A i (10)
其中 Γ = [ 1 cos τ 1 ⋯ cos N τ 1 sin τ 1 ⋯ sin N τ 1 ⋮ ⋮ ⋱ ⋮ ⋮ ⋱ ⋮ 1 cos τ H ⋯ cos N τ H sin τ H ⋯ sin N τ H ] ,把(10)代入(8)得到:
r ¯ = [ r ( τ 1 ) , r ( τ 2 ) , ⋯ , r ( τ H ) ] T = Γ r F , r F = [ r 0 i , 2 r 1 i , ⋯ , 2 r N i , 2 q 1 i , ⋯ , 2 q N i ] T 。
第三步:对于r利用离散的Fourier变换得到: r F = Γ − 1 r ¯ ,其中
Γ − 1 = 1 H ( 1 ⋯ 1 2 cos τ 1 ⋯ 2 cos ( τ H ) ⋮ ⋱ ⋮ 2 cos ( N τ 1 ) ⋯ 2 cos ( N τ H ) 2 sin τ 1 ⋯ 2 sin ( τ H ) ⋮ ⋱ ⋮ 2 sin ( N τ 1 ) ⋯ 2 sin ( N τ H ) ) ,
则:
R i = 1 2 π ∫ 0 2 π r T T d t = [ r 0 i , r 1 i , ⋯ , r N i , q 1 i , ⋯ , q N i ] T (11)
第四步: Δ A i = − B i − 1 R i ,则:
A i + 1 = A i + Δ A i (12)
若 ‖ Δ A i ‖ < δ 成立,则迭代结束,这里的 ‖ Δ A i ‖ 是当前增量的模, δ 为容许误差。
若不满足,则下一次迭代的近似Jacobian矩阵为:
B i + 1 = B i + ( R i + 1 − R i − B i Δ A i ) Δ A i T Δ A i T Δ A i (13)
一直重复上面的过程,直至收敛。
本文中考虑的系统具有强非线性和非光滑的特性,为了保证系统解的精确性,这里Fourier级数的阶数取为40。并且利用数值积分的方法来验证该方法的有效性。
当 Ω = 1 , ξ = 0.1 , β = 1 , α = 1.6 时,系统的稳态速度与时间的关系如图2所示,Galerkin表示本文的方法,*表示数值积分的方法。可以看到,系统在运动过程中不存在粘–滑状态转换,仅有滑动状态,两种方法的计算结果可以有比较好的吻合。
当 Ω = 1 , ξ = 0.1 , β = 1 , α = 1.4 时,系统的稳态速度与时间的关系如图3所示。可以看到:系统在这组参数下,其运动过程中存在粘–滑状态的转换,两种方法计算的结果也有很好的吻合。
图2. 系统的稳态速度与时间的关系
图3. 系统的稳态速度与时间的关系
本文主要研究了Galerkin法在一类干摩擦振动系统中的应用。其主要思想是:把系统的周期解设为Fourier级数,Galerkin过程把系统的微分方程转化为非线性的代数方程组,然后利用Broyden法求解非线性代数方程组,并利用数值积分的结果验证该方法的有效性。但是Boryden法有其自身的局限性,就是该方法强烈地依赖于给定的初始值和初始的Jacobain矩阵。因此寻找更有效的算法来求解此类代数方程组是将来更进一步的研究内容。
栗 鹏. 应用于干摩擦振动系统中的Galerkin法Galerkin Method Applied to the Vibration System with Dry Friction[J]. 声学与振动, 2019, 07(04): 155-160. https://doi.org/10.12677/OJAV.2019.74017